We present a set of self-consistent chemo-dynamical simulations of MW-type galaxies formation to study the origin of the bimodality of $\alpha$-elements in stellar populations. We explore how the bimodality is related to the geometrically and kinematically defined stellar discs, gas accretion and radial migration. We find that the two $\alpha$-sequences are formed in quite different physical environments. The high-$\alpha$ sequence is formed early from a burst of star formation (SF) in a turbulent, compact gaseous disc which forms a thick disc. The low-$\alpha$ stellar populations is the result of quiescent SF supported by the slow accretion of enriched gas onto a radially extended thin disc. Stellar feedback-driven outflows during the formation of the thick disc are responsible for the enrichment of the surrounding gaseous halo, which subsequently feeds the disc on a longer time-scale. During the thin disc phase, chemical evolution reaches an equilibrium metallicity and abundance, where the stars pile-up. This equilibrium metallicity decreases towards the outer disc, generating the ridgeline that forms the low-$\alpha$ sequence. We identify a second mechanism capable of creating a low-$\alpha$ sequence in one of our simulations. Rapid shutdown of the SF, provoked by the feedback at the end of the thick disc phase, suppresses the chemical enrichment of the halo gas, which, once accreted onto the star-forming disc, dilutes the ISM at the beginning of the thin disc formation. Both mechanisms can operate in a galaxy, but the former is expected to occur when SF efficiency ceases to be dominated by the formation of the thick disc, while the latter can occur in the inner regions. Being the result of the presence of low and high gas density environments, the bimodality is independent of any particular merger history, suggesting that it could be much more widespread than has been claimed.

Understanding radial migration is a crucial point to build relevant chemical and dynamical evolution models of the Milky Way disk. In this paper, we analyze a high-resolution N-body simulation of a Milky Way-type galaxy to study the role that the slowing down of a stellar bar has is generating migration from the inner to the outer disk. Stellar particles are trapped by the main resonances (corotation and Outer Lindblad resonance) which then propagate outwards across the disk due to the bar slowing down. Once the bar strength reaches its maximal amplitude, some of the stars, delivered to the outer disk, escape the resonances and some of them settle on nearly circular orbits. The number of the escaped stars gradually increases also due to the decrease of the bar strength when the boxy/peanut bulge forms. We show that this mechanism is not limited only to stars on nearly circular orbits: also stars initially on more eccentric orbits can be transferred outwards (out to the OLR location) and can end up on nearly circular orbits. Therefore, the propagation of the bar resonances outwards can induce the circularization of the orbits of some of the migrating stars. The mechanism investigated in this paper can explain the presence of metal-rich stars at the solar vicinity and more generally in the outer galactic disk. Our dynamical model predicts that up to 3% of stars in between of corotation and the OLR can be formed in the innermost region of the Milky Way. The epoch of the Milky Way bar formation can be potentially constrained by analyzing the age distribution of the most metal-rich stars at the solar vicinity.

We analyse a set of very metal-poor stars for which accurate chemical abundances have been obtained as part of the ESO Large Program "First stars" in the light of the Gaia DR2 data. The kinematics and orbital properties of the stars in the sample show they probably belong to the thick disc, partially heated to halo kinematics, and to the accreted Nissen & Schuster-Gaia Sausage-Enceladus satellite. The continuity of these properties with stars at both higher ($\rm [Fe/H]>-2$) and lower metallicities ($\rm [Fe/H]<-4.$) suggests that the Galaxy at $\rm [Fe/H] \lesssim -0.5$ and down to at least $\rm [Fe/H]\sim-6$ is dominated by these two populations. In particular, we show that the disc extends continuously from $\rm [Fe/H] \le -4$ (where stars with disc-like kinematics have been recently discovered) up to $\rm [Fe/H] \ge -2$, the metallicity regime of the Galactic thick disc. There exists indeed an "ultra-metal poor thick disc", which constitutes the extremely metal-poor tail of the canonical Galactic thick disc, and which extends this latter from $\rm [Fe/H] \sim -0.5$ up to the most metal-poor stars discovered in the Galaxy up to date. These results suggest that the disc may be the main, and possibly the only stellar population that has formed in the Galaxy at these metallicities. This requires that the dissipative collapse that led to the formation of the old Galactic disc must have been extremely fast. We also discuss these results in the light of recent simulation efforts made to reproduce the first stages of Milky Way-type galaxies.

In this paper we introduce a new method for analysing Milky Way phase-space which allows us to reveal the imprint left by the Milky Way bar and spiral arms on the stars with full phase-space data in Gaia Data Release 2. The unprecedented quality and extended spatial coverage of these data enable us to discover six prominent stellar density structures in the disc to a distance of 5 kpc from the Sun. Four of these structures correspond to the spiral arms detected previously in the gas and young stars (Scutum-Centaurus, Sagittarius, Local and Perseus). The remaining two are associated with the main resonances of the Milky Way bar where corotation is placed at around 6.2 kpc and the outer Lindblad resonance beyond the Solar radius, at around 9 kpc. For the first time we provide evidence of the imprint left by spiral arms and resonances in the stellar densities not relying on a specific tracer, through enhancing the signatures left by these asymmetries. Our method offers new avenues for studying how the stellar populations in our Galaxy are shaped.

As the remnants of stars with initial masses $\lesssim$ 8 M$_{\odot}$, white dwarfs contain valuable information on the formation histories of stellar populations. In this paper, we use deep, high-quality, u-band photometry from the Canada France Imaging Survey (CFIS), griz photometry from Pan-STARRS 1 (PS1), as well as proper motions from Gaia DR2, to select 25,156 white dwarf candidates over $\sim$4500 deg$^2$ using a reduced proper motion diagram. We develop a new white dwarf population synthesis code that returns mock observations of the Galactic field white dwarf population for a given star formation history, while simultaneously taking into account the geometry of the Milky Way, survey parameters, and selection effects. We use this model to derive the star formation histories of the thin disk, thick disk, and stellar halo. Our results show that the Milky Way disk began forming stars (11.3 $\pm$ 0.5) Gyr ago, with a peak rate of (8.8 $\pm$ 1.4) M$_{\odot}$yr$^{-1}$ at (9.8 $\pm$ 0.4) Gyr, before a slow decline to a constant rate until the present day --- consistent with recent results suggesting a merging event with a satellite galaxy. Studying the residuals between the data and best-fit model shows evidence for a slight increase in star formation over the past 3 Gyr. We fit the local fraction of helium-atmosphere white dwarfs to be (21 $\pm$ 3) %. Incorporating this methodology with data from future wide-field surveys such as LSST, Euclid, CASTOR, and WFIRST should provide an unprecedented view into the formation of the Milky Way at its earliest epoch through its white dwarfs.

Previous studies based on the analysis of Gaia DR2 data have revealed that accreted stars, possibly originating from a single progenitor satellite, are a significant component of the halo of our Galaxy, potentially constituting most of the halo stars at $\rm [Fe/H] < -1$ within a few kpc from the Sun and beyond. In this paper, we couple astrometric data from Gaia DR2 with elemental abundances from APOGEE DR14 to characterize the kinematics and chemistry of in-situ and accreted populations up to $\rm [Fe/H] \sim -2$. Accreted stars appear to significantly impact the Galactic chemo-kinematic relations, not only at $\rm [Fe/H] < -1$, but also at metallicities typical of the thick and metal-poor thin discs. They constitute about 60% of all stars at $\rm [Fe/H] < -1$, the remaining 40% being made of (metal-weak) thick disc stars. We find that the stellar kinematic fossil record shows the imprint left by this accretion event which heated the old Galactic disc. We are able to age-date this kinematic imprint, showing that the accretion occurred between 9 and 11 Gyr ago, and that it led to the last significant heating of the Galactic disc. An important fraction of stars with abundances typical of the (metal-rich) thick disc, and heated by this interaction, is now found in the Galactic halo. Indeed about half of the kinematically defined halo at few kpc from the Sun is composed of metal-rich thick disc stars. Moreover, we suggest that this metal-rich thick disc component dominates the stellar halo of the inner Galaxy. The new picture that emerges from this study is one where the standard non-rotating in-situ halo population, the collapsed halo, seems to be more elusive than ever.

The period 2035-50 considered in the ESA Voyage long-term plan will coincide with a series of foreseeable advances in the characterization of the stellar content of the Milky Way. The Gaia mission, combined with large-scale spectroscopic surveys, is helping to build an unprecedented census in terms of the astrometric, kinematic and chemical properties of Galactic stellar populations. Within a decade, precise measurements of such properties will be available for hundreds of millions of stars. Meanwhile, time-domain surveys initiated with CoRoT and Kepler/K2 and carried on by space missions such as TESS and PLATO or ground-based projects like the LSST, will have brought asteroseismology to a high level of maturity. The combination of precise ages from asteroseismology with astrometric and spectroscopic data, on large stellar samples, is allowing Galactic archaeologists to gain new insight into the assembly history of the Milky Way. Recent breakthroughs --- based on the detection of solar-like oscillations in tens of thousands of red-giant stars --- demonstrate the potential of such approach. Therefore, we are convinced that an all-sky, high-cadence, long-duration stellar variability survey will become a scientific priority in the 2035-50 period. The Chronos concept presented here consists in a time-domain extension to Gaia. It will allow for mass and age estimates for half a million red giants within 1.7 kpc from the Sun and hence shed a new light on our understanding of the Galactic dynamics and archaeology. In terms of the targeted pulsators, Chronos will bridge the gap between PLATO and the LSST by surveying stars all the way from the subgiant branch to the early AGB. Finally, it will surpass all previous surveys capable of conducting asteroseismology in terms of the combined sky coverage and duration of the observations (2 x 3.75 months over the whole sky and >5 years in the CVZ).

We explore the origin of phase-space substructures revealed by the second Gaia data release in the disc of the Milky Way, such as the ridges in the $V_{\phi}$-$r$ plane, the undulations in the $V_{\phi}$-$r$-$V_r$ space and the streams in the $V_{\phi}$-$V_r$ plane. We use a collisionless N-body simulation with co-spatial thin and thick discs, along with orbit integration, to study the orbital structure close to the Outer Lindblad Resonance (OLR) of the bar. We find that a prominent, long-lived ridge is formed in the $V_{\phi}$-$r$ plane due to the OLR which translates to streams in the $V_{\phi}$-$V_r$ plane and examine which closed periodic and trapped librating orbits are responsible for these features. We find that orbits which carry out small librations around the $x_1(1)$ family are preferentially found at negative $V_r$, giving rise to a `horn'-like feature, while orbits with larger libration amplitudes, trapped around the $x_1(2)$ and $x_1(1)$ families, constitute the positive $V_r$ substructure, i.e. the Hercules-like feature. This changing libration amplitude of orbits will translate to a changing ratio of thin/thick disc stars, which could have implications on the metallicity distribution in this plane. We find that a scenario in which the Sun is placed close to the OLR gives rise to a strong asymmetry in $V_r$ in the $V_{\phi}$-$V_r$ plane (i.e. Hercules vs. `the horn') and subsequently to undulations in the $V_{\phi}$-$r$-$V_r$ space. We also explore a scenario in which the Sun is placed closer to the bar corotation and find that the bar perturbation $alone$ cannot give rise to the these features.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

We present a scenario of the chemical enrichment of the solar neighborhood that solves the G-dwarf problem by taking into account constraints on a larger scale. We argue that the Milky Way disk within 10 kpc has been enriched to solar metallicity by a massive stellar population: the thick disk, which itself formed from a massive turbulent gaseous disk. The key new consideration is that the pre-enrichment provided by the thick disk is not related to the mass fraction of this stellar population at the solar radius, as is classically assumed in inside-out scenarios, but is actually related to the formation of the entire massive thick disk, due to the vigorous gas phase mixing that occurred during its formation. Hence, the fact that this population represents only 15-25% of the local stellar surface density today is irrelevant for `solving' the G-dwarf problem. The only condition for this scenario to work is that the thick disk was formed from a turbulent gaseous disk that permitted a homogeneous -- not radially dependent -- distribution of metals, allowing the solar ring to be enriched to solar metallicity. At the solar radius, the gas flowing from the outer disk combined with the solar metallicity gas left over from thick disk formation, providing the fuel necessary to form the thin disk at the correct metallicity to solve the G-dwarf problem. Chemical evolution at R$>$6 kpc, and in particular beyond the solar radius, can be reproduced with the same scheme. These results imply that the local metallicity distribution is not connected to the gas accretion history of the Milky Way. Finally, we argue that the Sun is the result of the evolution typical of stars in the disk beyond $\sim$6 kpc (i.e., also undergoing dilution), and has none of the characteristics of inner disk stars. [Abridged]

Using a single N-body simulation ($N=0.14\times 10^9$) we explore the formation, evolution and spatial variation of the phase-space spirals similar to those recently discovered by Antoja et al. in the Milky Way disk, with Gaia DR2. For the first time in the literature, we use a self-consistent N-body simulation of an isolated Milky Way-type galaxy to show that the phase-space spirals develop naturally from vertical oscillations driven by the buckling of the stellar bar. We claim that the physical mechanism standing behind the observed incomplete phase-space mixing process can be internal and not necessarily due to the perturbation induced by a massive satellite. In our model, the bending oscillations propagate outwards and produce axisymmetric variations of the mean vertical coordinate and of the vertical velocity component. As a consequence, the phase-space wrapping results in the formation of patterns with various morphology across the disk, depending on the bar orientation, distance to the galactic center and time elapsed since the bar buckling. Once bending waves appear, they are supported for a long time via disk self-gravity. The underlying physical mechanism implies the link between in-plane and vertical motion that leads directly to phase-space structures whose amplitude and shape are in remarkable agreement with those of the phase-space spirals observed in the Milky Way disk. In our isolated galaxy simulation, phase-space spirals are still distinguishable, at the solar neighbourhood, 3 Gyr after the buckling phase. The long-lived character of the phase-space spirals generated by the bar buckling instability cast doubts on the timing argument used so far to get back at the time of the onset of the perturbation: phase-space spirals may have been caused by perturbations originated several Gyrs ago, and not as recent as suggested so far.

In the MW bulge, metal-rich stars form a strong bar and are more peanut-shaped than metal-poor stars. It has been recently claimed that this behavior is driven by the initial in-plane radial velocity dispersion of these populations, rather than by their initial vertical random motions. This has led to the suggestion that a thick disc is not necessary to explain the characteristics of the MW bulge. We rediscuss this issue by analyzing two dissipationless N-body simulations of boxy/peanut (b/p)-shaped bulges formed from composite stellar discs, made of kinematically cold and hot stellar populations, and we conclude that initial vertical random motions are as important as in-plane random motions in determining the relative contribution of cold and hot disc populations with height above the plane, the metallicity and age trends. Previous statements emphasizing the dominant role of in-plane motions in determining these trends are not confirmed. However, differences exist in the morphology and strength of the resulting b/p-shaped bulges: a model where disc populations have initially only different in-plane random motions, but similar thickness, results into a b/p bulge where all populations have a similar peanut shape, independently on their initial kinematics, or metallicity. We discuss the reasons behind these differences, and also predict the signatures that these two extreme initial conditions would leave on the vertical age and metallicity gradients of disc stars, outside the bulge region. We conclude that a metal-poor, kinematically (radial and vertical) hot component, that is a thick disc, is necessary in the MW before bar formation, supporting the scenario traced in previous works. [abridged]

Globular clusters (GCs), the oldest stellar systems observed in the Milky Way, have for long been considered single stellar populations. As such, they provided an ideal laboratory to understand stellar dynamics and primordial star formation processes. However, during the last two decades, observations unveiled their real, complex nature. Beside their pristine stars, GCs host one or more helium enriched and possibly younger stellar populations whose formation mechanism is still unknown. Even more puzzling is the existence of GCs showing star by star iron spreads. Using detailed N-body simulations we explore the hypothesis that these anomalies in metallicity could be the result of mutual stripping and mergers between a primordial population of disc GCs. In the first paper of this series we proved, both with analytical arguments and short-term N-body simulations, that disc GCs have larger fly-by and close encounter rates with respect to halo clusters. These interactions lead to mass exchange and even mergers that form new GCs, possibly showing metallicity spreads. Here, by means of long-term direct N-body simulations, we provide predictions on the dynamical properties of GCs that underwent these processes. The comparison of our predictions with available and future observational data could provide insights on the origin of GCs and on the Milky Way build-up history as a whole.

We present the results of a self-consistent $N$-body simulation following the evolution of a primordial population of thick disc globular clusters (GCs). We study how the internal properties of such clusters evolve under the action of mutual interactions, while they orbit a Milky Way-like galaxy. For the first time, through analytical and numerical considerations, we find that physical encounters between disc GCs are a crucial factor that contributed to the shape of the current properties of the Galactic GC system. Close passages or motion on similar orbits may indeed have a significant impact on the internal structure of clusters, producing multiple gravitationally bound sub-populations through the exchange of mass and even mergers. Our model produces two major mergers and a few small mass exchanges between pairs of GCs. Two of our GCs accrete stars from two companions, ending up with three internal sub-populations. We propose these early interactions and mergers between thick disc GCs with slightly different initial chemical compositions as a possible explanation for the presence of the spreads in metallicity observed in some of the massive Milky Ways GCs.

We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.

We investigate the nature of the double color-magnitude sequence observed in the Gaia DR2 HR diagram of stars with high transverse velocities. The stars in the reddest-color sequence are likely dominated by the dynamically-hot tail of the thick disk population. Information from Nissen & Schuster (2010) and from the APOGEE survey suggests that stars in the blue-color sequence have elemental abundance patterns that can be explained by this population having a relatively low star-formation efficiency during its formation. In dynamical and orbital spaces, such as the `Toomre diagram', the two sequences show a significant overlap, but with a tendency for stars on the blue-color sequence to dominate regions with no or retrograde rotation and high total orbital energy. In the plane defined by the maximal vertical excursion of the orbits versus their apocenters, stars of both sequences redistribute into discrete wedges. We conclude that stars which are typically assigned to the halo in the solar vicinity are actually both accreted stars lying along the blue sequence in the HR diagram, and the low velocity tail of the old Galactic disk, possibly dynamically heated by past accretion events. Our results imply that a halo population formed in situ and responsible for the early chemical enrichment prior to the formation of the thick disk is yet to be robustly identified, and that what has been defined as the stars of the in situ stellar halo of the Galaxy may be in fact fossil records of its last significant merger.

We use planetary nebulae (PNe) as probes to determine the Galactic radial oxygen gradients, and other abundance patterns. We select data homogeneously from recent data sets, including PNe at large Galactocentric distances. The radial oxygen gradient calculated for the general PN population, which probes the region between the Galactic center and out to $\sim$28 kpc, is shallow, with slope $\sim$-0.02 dex kpc$^{-1}$, in agreement with previous findings. We looked for time evolution of the metallicity gradient using PNe with different age progenitors as metallicity probes. We identify PNe whose progenitor stars are younger than 1 Gyr (YPPNe), and those whose progenitor stars are older than 7.5 Gyr (OPPNe), based on the comparison between evolutionary yields and elemental abundances of the PNe. By studying OPPNe and YPPNe separately we found that: (i) The OPPNe oxygen gradient is shallower ($\sim-0.015$ dex kpc$^{-1}$) than that derived from YPPNe ($\sim-0.027$ dex kpc$^{-1}$); (ii) the OPPNe inner radial distribution of oxygen is compatible with no gradient to the radial extent of the thick disk population ($\sim$10 kpc), similarly to what has been observed in thick disk stars; (iii) planetary nebulae (especially OPPNe) indicate that significant gradient slope is limited to Galactocentric distances between $\sim$10 to $\sim$13.5 kpc, as observed for open clusters and field stars. Outside this range, the distribution is almost flat. We found that the radial oxygen gradient is steeper for a PN population closer to the Galactic disk, similarly to what is observed in the general stellar population by the SEGUE survey. We use our novel population dating to compare our results with current chemical evolutionary models, and with gradients from other Galactic populations, for insight on galaxy chemical evolution.

There is a long-standing debate on the origin of the metal-poor stellar populations of the Milky Way (MW) bulge, with the two leading scenarios being that these populations are either i) part of a classical metal-poor spheroid or ii) the same population as the chemically defined thick disc seen at the Solar neighbourhood. Here we test whether the latter scenario can reproduce the observed chemical properties of the MW bulge. To do so we compare an N-body simulation of a composite (thin+thick) stellar disc -- which evolves secularly to form a bar and a boxy/peanut (b/p) bulge -- to data from APOGEE DR13. This model, in which the thick disc is massive and centrally concentrated, can reproduce the morphology of the metal-rich and metal-poor stellar populations in the bulge, as well as the mean metallicity and [$\alpha$/Fe] maps as obtained from the APOGEE data. It also reproduces the trends, in both longitude and latitude, of the bulge metallicity distribution function (MDF). Additionally, we show that the model predicts small but measurable azimuthal metallicity variations in the inner disc due to the differential mapping of the thin and thick disc in the bar. We therefore see that the chemo-morphological relations of stellar populations in the MW bulge are naturally reproduced by mapping the thin and thick discs of the inner MW into a b/p.

Despite the huge amount of photometric and spectroscopic efforts targetting the Galactic bulge over the past few years, its age distribution remains controversial owing to both the complexity of determining the age of individual stars and the difficult observing conditions. Taking advantage of the recent release of very deep, proper-motion-cleaned colour--magnitude diagrams (CMDs) of four low reddening windows obtained with the Hubble Space Telescope (HST), we used the CMD-fitting technique to calculate the star formation history (SFH) of the bulge at -2deg > b > -4deg along the minor axis. We find that over 80 percent of the stars formed before 8 Gyr ago, but that a significant fraction of the super-solar metallicity stars are younger than this age. Considering only the stars that are within reach of the current generation of spectrographs (i.e. V < 21), we find that 10 percent of the bulge stars are younger than 5 Gyr, while this fraction rises to 20-25 percent in the metal-rich peak. The age-metallicity relation is well parametrized by a linear fit implying an enrichment rate of dZ/dt ~ 0.005 Gyr$^{-1}$. Our metallicity distribution function accurately reproduces that observed by several spectroscopic surveys of Baade's window, with the bulk of stars having metal-content in the range [Fe/H] ~ -0.7 to ~0.6, along with a sparse tail to much lower metallicities.

Although there is consensus that metal-rich stars in the Milky Way bulge are formed via secular evolution of the thin disc, the origin of their metal-poor counterparts is still under debate. Two different origins have been invoked for metal-poor stars: they might be classical bulge stars or stars formed via internal evolution of a massive thick disc. We use N-body simulations to calculate the kinematic signature given by the difference in the mean Galactocentric radial velocity ($\Delta\rm V_{\rm GC}$) between metal-rich stars ([Fe/H] $\ge$ 0) and moderately metal-poor stars (-1.0 $\le$ [Fe/H] $<$ 0) in two models, one containing a thin disc and a small classical bulge (B/D=0.1), and the other containing a thin disc and a massive centrally concentrated thick disc. We reasonably assume that thin-disk stars in each model may be considered as a proxy of metal-rich stars. Similarly, bulge stars and thick-disc stars may be considered as a proxy of metal-poor stars. We calculate $\Delta\rm V_{\rm GC}$ at different latitudes ($b=0^\circ$, $-2^\circ$, $-4^\circ$,$-6^\circ$, $-8^\circ$ and $-10^\circ$) and longitudes ($l=0^\circ$, $\pm5^\circ$, $\pm10^\circ$ and $\pm15^\circ$) and show that the $\Delta\rm V_{\rm GC}$ trends predicted by the two models are different. We compare the predicted results with ARGOS data and APOGEE DR13 data and show that moderately metal-poor stars are well reproduced with the co-spatial stellar discs model, which has a massive thick disc. Our results give more evidence against the scenario that most of the metal-poor stars are classical bulge stars. If classical bulge stars exists, most of them probably have metallicities [Fe/H] $<$ -1 dex, and their contribution to the mass of the bulge should be a small percentage of the total bulge mass.

We study how migration affects stars of a galaxy with a thin stellar disc and thicker stellar components. The simulated galaxy has a strong bar and lasting spiral arms. We find that the amplitude of the churning (change in angular momentum) is similar for thin and thick components, and of limited amplitude, and that stars of all components can be trapped at the corotation of the bar. At the exception of those stars trapped at the corotation, we find that stars that are far from their initial guiding radius are more likely so due to blurring rather than churning effects. We compare the simulation to orbits integration with a fixed gravitational potential rotating at a constant speed. In the latter case, stars trapped at corotation are churned periodically outside and inside the corotation radius, with a zero net average. However, as the bar speed of the simulated galaxy decreases and its corotation radius increases, stars trapped at corotation for several Gyrs can be churned outwards on average. We study the location of extreme migrators (stars experimenting the largest churning) and find that extreme migrators come from regions on the leading side of the effective potential local maxima.

We show that the bulge and the disk of the Milky Way (MW) at R$\lesssim$7~kpc are well described by a unique chemical evolution and a two-phase star-formation history (SFH). We argue that the populations within this inner disk, not the entire disk, are the same, and that the outer Lindblad resonance (OLR) of the bar plays a key role in explaining this uniformity. In our model of a two-phase star formation history, the metallicity, [$\alpha$/Fe] and [$\alpha$/H] distributions, and age-metallicity relation are all compatible with the observations of both the inner disk and bulge. The dip at [Fe/H]$\sim$0 dex seen in the metallicity distributions of the bulge and inner disk reflects the quenching episode in the SFH of the inner MW at age $\sim$8 Gyr, and the common evolution of the bulge and inner disk stars. We show that at z$\le$1.5, when the MW was starting to quench, transitioning between the end of the $\alpha$-enhanced thick disk formation to the start of the thin disk, and yet was still gas rich, the gas accretion rate could not have been significant. The [$\alpha$/Fe] abundance ratio before and after this quenching phase would be different, which is not observed. The present analysis suggests that Milky Way history, and in particular at the transition from the thick to the thin disk -- the epoch of the quenching -- must have been driven by a decrease of the star formation efficiency. We argue that the decline in the intensity of gas accretion, the formation of the bar, and the quenching of the SFR at the same epoch may be causally connected thus explaining their temporal coincidence (abridged).

This letter studies the formation of azimuthal metallicity variations in the disks of spiral galaxies in the absence of initial radial metallicity gradients. Using high-resolution $N$-body simulations, we model composite stellar discs, made of kinematically cold and hot stellar populations, and study their response to spiral arm perturbations. We find that, as expected, disk populations with different kinematics respond differently to a spiral perturbation, with the tendency for dynamically cooler populations to show a larger fractional contribution to spiral arms than dynamically hotter populations. By assuming a relation between kinematics and metallicity, namely the hotter the population, the more metal-poor it is, this differential response to the spiral arm perturbations naturally leads to azimuthal variations in the mean metallicity of stars in the simulated disk. Thus, azimuthal variations in the mean metallicity of stars across a spiral galaxy are not necessarily a consequence of the reshaping, by radial migration, of an initial radial metallicity gradient. They indeed arise naturally also in stellar disks which have initially only a negative vertical metallicity gradient.

We examine the metallicity trends in the Milky Way (MW) bulge - using APOGEE DR13 data - and explore their origin by comparing two N-body models of isolated galaxies which develop a bar and a boxy/peanut (b/p) bulge. Both models have been proposed as scenarios for reconciling a disc origin of the MW bulge with a negative vertical metallicity gradient. The first model is a superposition of co-spatial, i.e. overlapping, disc populations with different scaleheights, kinematics and metallicities. In this model the thick, metal-poor, and centrally concentrated disc populations contribute significantly to the stellar mass budget in the inner galaxy. The second model is a single disc with an initial steep radial metallicity gradient, which is mapped by the bar into the b/p bulge in such a way that the vertical metallicity gradient of the MW bulge is reproduced -- as shown already in previous works in the literature. However, as we show here, the latter model does not reproduce the positive longitudinal metallicity gradient of the inner disc, nor the metal-poor innermost regions of the Bulge seen in the data. On the other hand, the model with co-spatial thin and thick disc populations reproduces all the aforementioned trends. We therefore see that it is possible to reconcile a (primarily) disc origin for the MW bulge with the observed trends in metallicity by mapping the inner thin and thick discs of the MW into a b/p. For this scenario to reproduce the observations, the $\alpha$-enhanced, metal-poor, thick disc populations must have a significant mass contribution in the inner regions -- as has been suggested for the Milky Way.

Galaxy surveys have suggested that rapid and sustained decrease in the star-formation rate, "quenching", in massive disk galaxies is frequently related to the presence of a bar. Optical and near-IR observations reveal that nearly 60% of disk galaxies in the local universe are barred, thus it is important to understand the relationship between bars and star formation in disk galaxies. Recent observational results imply that the Milky Way quenched about 9-10 Gyr ago, at the transition between the cessation of the growth of the kinematically hot, old, metal-poor thick disk and the kinematically colder, younger, and more metal-rich thin disk. Although perhaps coincidental, the quenching episode could also be related to the formation of the bar. To explore the relation between bar formation and star formation in gas rich galaxies quantitatively, we simulated gas-rich disk isolated galaxies. Our simulations include prescriptions for star formation, stellar feedback, and for regulating the multi-phase interstellar medium. We find that the action of stellar bar efficiently quenches star formation, reducing the star-formation rate by a factor of 10 in less than 1 Gyr. Analytical and self-consistent galaxy simulations with bars suggest that the action of the stellar bar increases the gas random motions within the co-rotation radius of the bar. Indeed, we detect an increase in the gas velocity dispersion up to 20-35 km/s at the end of the bar formation phase. The star-formation efficiency decreases rapidly, and in all of our models, the bar quenches the star formation in the galaxy. The star-formation efficiency is much lower in simulated barred compared to unbarred galaxies and more rapid bar formation implies more rapid quenching.

The Canada-France Imaging Survey (CFIS) will map the northern high Galactic latitude sky in the $u$-band ("CFIS-u", 10,000$\, {\rm deg^2}$) and in the $r$-band ("CFIS-r", 5,000$\, {\rm deg^2}$), enabling a host of stand-alone science investigations, and providing some of the ground-based data necessary for photometric redshift determination for the Euclid mission. In this first contribution we present the $u$-band component of the survey, describe the observational strategy, and discuss some first highlight results, based on approximately one third of the final area. We show that the Galactic anticenter structure is distributed continuously along the line of sight, out to beyond 20 kpc, and possesses a metallicity distribution that is essentially identical to that of the outer disk sampled by APOGEE. This suggests that it is probably a buckled disk of old metal-rich stars, rather than a stream or a flare. We also discuss the future potential for CFIS-u in discovering star-forming dwarf galaxies around the Local Group, the characterization of the white dwarf and blue straggler population of the Milky Way, as well as its sensitivity to low-surface brightness structures in external galaxies.

We present the chemical distribution of the Milky Way, based on 2,900$\, {\rm deg^2}$ of $u$-band photometry taken as part of the Canada-France Imaging Survey. When complete, this survey will cover 10,000$\, {\rm deg^2}$ of the Northern sky. By combing the CFHT $u$-band photometry together with SDSS and Pan-STARRS $g,r,$ and $i$, we demonstrate that we are able to measure reliably the metallicities of individual stars to $\sim 0.2$ dex, and hence additionally obtain good photometric distance estimates. This survey thus permits the measurement of metallicities and distances of the dominant main-sequence population out to approximately 30 kpc, and provides much higher number of stars at large extraplanar distances than have been available from previous surveys. We develop a non-parametric distance-metallicity decomposition algorithm and apply it to the sky at $30\deg < |b| < 70\deg$ and to the North Galactic Cap. We find that the metallicity-distance distribution is well-represented by three populations whose metallicity distributions do not vary significantly with vertical height above the disk. As traced in main-sequence stars, the stellar halo component shows a vertical density profile that is close to exponential, with a scale height of around 3 kpc. This may indicate that the inner halo was formed partly from disk stars ejected in an ancient minor merger.

We investigate here the effect of the selection function on the metallicity distribution function (MDF) and on the vertical metallicity gradient by studying similar lines of sight using four different spectroscopic surveys (APOGEE, LAMOST, RAVE, and Gaia-ESO), which have different targeting strategies and therefore different selection functions. We use common fields between the spectroscopic surveys of APOGEE, LAMOST, RAVE (ALR) and APOGEE, RAVE, Gaia-ESO (AGR) and use two stellar population synthesis models, GALAXIA and TRILEGAL, to create mock fields for each survey. We apply the selection function in the form of colour and magnitude cuts of the respective survey to the mock fields to replicate the observed source sample. We make a basic comparison between the models to check which best reproduces the observed sample distribution. We carry out a quantitative comparison between the synthetic MDF from the mock catalogues using both models to understand the effect of the selection function on the MDF and on the vertical metallicity gradient. Using both models, we find a negligible effect of the selection function on the MDF for APOGEE, LAMOST, and RAVE. We find a negligible selection function effect on the vertical metallicity gradients as well, though GALAXIA and TRILEGAL have steeper and shallower slopes, respectively, than the observed gradient. After applying correction terms on the metallicities of RAVE and LAMOST with respect to our reference APOGEE sample, our observed vertical metallicity gradients between the four surveys are consistent within 1-sigma. We also find consistent gradient for the combined sample of all surveys in ALR and AGR. We estimated a mean vertical metallicity gradient of -0.241+/-0.028 dex kpc-1. There is a significant scatter in the estimated gradients in the literature, but our estimates are within their ranges.

Deciphering the assembly history of the Milky Way is a formidable task, which becomes possible only if one can produce high-resolution chrono-chemo-kinematical maps of the Galaxy. Data from large-scale astrometric and spectroscopic surveys will soon provide us with a well-defined view of the current chemo-kinematical structure of the Milky Way, but will only enable a blurred view on the temporal sequence that led to the present-day Galaxy. As demonstrated by the (ongoing) exploitation of data from the pioneering photometric missions CoRoT, Kepler, and K2, asteroseismology provides the way forward: solar-like oscillating giants are excellent evolutionary clocks thanks to the availability of seismic constraints on their mass and to the tight age-initial-mass relation they adhere to. In this paper we identify five key outstanding questions relating to the formation and evolution of the Milky Way that will need precise and accurate ages for large samples of stars to be addressed, and we identify the requirements in terms of number of targets and the precision on the stellar properties that are needed to tackle such questions. By quantifying the asteroseismic yields expected from PLATO for red-giant stars, we demonstrate that these requirements are within the capabilities of the current instrument design, provided that observations are sufficiently long to identify the evolutionary state and allow robust and precise determination of acoustic-mode frequencies. This will allow us to harvest data of sufficient quality to reach a 10% precision in age. This is a fundamental pre-requisite to then reach the more ambitious goal of a similar level of accuracy, which will only be possible if we have to hand a careful appraisal of systematic uncertainties on age deriving from our limited understanding of stellar physics, a goal which conveniently falls within the main aims of PLATO's core science.

We explore trends in the morphology and line-of-sight (los) velocity of stellar populations in the inner regions of disc galaxies, using N-body simulations with both a thin (kinematically cold) and a thick (kinematically hot) disc which form a bar and boxy/peanut (b/p) bulge. The bar in the thin disc component is $\sim$50\% stronger than the thick disc bar and is more elongated, with an axis ratio almost half that of the thick disc bar. The thin disc b/p bulge has a pronounced X-shape, while the thick disc b/p is weaker with a rather boxy shape. This leads to the signature of the b/p bulge in the thick disc to be weaker and further away from the plane than in the thin disc. Regarding the kinematics, we find that the los velocity of thick disc stars in the outer parts of the b/p bulge can be \emph{larger} than that of thin disc stars, by up to 40\% and 20\% for side-on and Milky Way-like orientations of the bar respectively. This is due to the different orbits followed by thin and thick disc stars in the bar-b/p region, which are affected by the fact that: i) thin disc stars are trapped more efficiently in the bar - b/p instability and thus lose more angular momentum than their thick disc counterparts and ii) thick disc stars have large radial excursions and therefore stars from large radii with high angular momenta can be found in the bar region. We also find that the difference between the los velocities of the thin and thick disc in the b/p bulge ($\Delta v_{los}$) correlates with the initial difference between the radial velocity dispersions of the two discs ($\Delta \sigma$) . We therefore conclude that stars in the bar - b/p bulge will have considerably different morphologies and kinematics depending on the kinematic properties of the disc population they originate from.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

In this work, two new axisymmetric models for the Galactic mass distribution are presented. Motivated by recent results, these two models include the contribution of a stellar thin disc and of a thick disc, as massive as the thin counterpart but with a shorter scale-length. Both models satisfy a number of observational constraints: stellar densities at the solar vicinity, thin and thick disc scale lengths and heights, rotation curve(s), and the absolute value of the perpendicular force Kz as a function of distance to the Galactic centre. We numerically integrate into these new models the motion of all Galactic globular clusters for which distances, proper motions, and radial velocities are available, and the orbits of about one thousand stars in the solar vicinity. The retrieved orbital characteristics are compared to those obtained by integrating the clusters and stellar orbits in pure thin disc models. We find that, due to the possible presence of a thick disc, the computed orbital parameters of disc stars can vary by as much as 30-40%. We also show that the systematic uncertainties that affect the rotation curve still plague computed orbital parameters of globular clusters by similar amounts.

The $\Lambda$CDM cosmological scenario predicts that our Galaxy should contain hundreds of stellar streams at the solar vicinity, fossil relics of the merging history of the Milky Way and more generally of the hierarchical growth of galaxies. Because of the mixing time scales in the inner Galaxy, it has been claimed that these streams should be difficult to detect in configuration space but can still be identifiable in kinematic-related spaces like the energy/angular momenta spaces, E-Lz and Lperp-Lz, or spaces of orbital/velocity parameters. By means of high-resolution, dissipationless N-body simulations, containing between 25$\times10^6$ and 35$\times10^6$ particles, we model the accretion of a series of up to four 1:10 mass ratio satellites then up to eight 1:100 satellites and we search systematically for the signature of these accretions in these spaces. In all spaces considered (1) each satellite gives origin to several independent overdensities; (2) overdensities of multiple satellites overlap; (3) satellites of different masses can produce similar substructures; (4) the overlap between the in-situ and the accreted population is considerable everywhere; (5) in-situ stars also form substructures in response to the satellite(s) accretion. These points are valid even if the search is restricted to kinematically-selected halo stars only. As we are now entering the 'Gaia era', our results warn that an extreme caution must be employed before interpreting overdensities in any of those spaces as evidence of relics of accreted satellites. Reconstructing the accretion history of our Galaxy will require a substantial amount of accurate spectroscopic data, that, complemented by the kinematic information, will possibly allow us to (chemically) identify accreted streams and measure their orbital properties. (abridged)

The determination of the age of the bulge has led to two contradictory results. On the one side, the color-magnitude diagrams in different bulge fields seem to indicate a uniformly old ($>$10 Gyr) population. On the other side, individual ages derived from dwarfs observed through microlensing events seem to indicate a large spread, from $\sim$ 2 to $\sim$ 13 Gyr. Because the bulge is now recognised as being mainly a boxy peanut-shaped bar, it is suggested that disk stars are one of its main constituents, and therefore also stars with ages significantly younger than 10 Gyr. Other arguments as well point to the fact that the bulge cannot be exclusively old, and in particular cannot be a burst population, as it is usually expected if the bulge was the fossil remnant of a merger phase in the early Galaxy. In the present study, we show that given the range of metallicities observed in the bulge, a uniformly old population would be reflected into a significant spread in color at the turn-off which is not observed. Inversely, we demonstrate that the correlation between age and metallicity expected to hold for the inner disk would conspire to form a color-magnitude diagram with a remarkably small spread in color, thus mimicking the color-magnitude diagram of a uniformly old population. If stars younger than 10 Gyr are part of the bulge, as must be the case if the bulge has been mainly formed through dynamical instabilities in the disk, then a very small spread at the turn-off is expected, as seen in the observations.

Quenching, the cessation of star formation, is one of the most significant events in the life cycle of galaxies. We show here the first evidence that the Milky Way experienced a generalised quenching of its star formation at the end of its thick disk formation $\sim$9 Gyr ago. Elemental abundances of stars studied as part of the APOGEE survey reveal indeed that in less than $\sim$2 Gyr the star formation rate in our Galaxy dropped by an order-of-magnitude. Because of the tight correlation between age and alpha abundance, this event reflects in the dearth of stars along the inner disk sequence in the [Fe/H]-[$\alpha$/Fe] plane. Before this phase, which lasted about 1.5 Gyr, the Milky Way was actively forming stars. Afterwards, the star formation resumed at a much lower level to form the thin disk. These events are very well matched by the latest observation of MW-type progenitors at high redshifts. In late type galaxies, quenching is believed to be related to a long and secular exhaustion of gas. In our Galaxy, it occurred on a much shorter time scale, while the chemical continuity before and after the quenching indicates that it was not due to the exhaustion of the gas. While quenching is generally associated with spheroids, our results show that it also occurs in galaxies like the Milky Way, possibly when they are undergoing a morphological transition from thick to thin disks. Given the demographics of late type galaxies in the local universe, in which classical bulges are rare, we suggest further that this may hold true generally in galaxies with mass lower than or approximately $M^*$, where quenching could be directly a consequence of thick disk formation. We emphasize that the quenching phase in the Milky Way could be contemporaneous with, and related to, the formation of the bar. We sketch a scenario on how a strong bar may inhibit star formation.

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

(Abridged) By using an N-body simulation of a bulge that was formed via a bar instability mechanism, we analyse the imprints of the initial (i.e. before bar formation) location of stars on the bulge kinematics, in particular on the heliocentric radial velocity distribution of bulge stars. Four different latitudes were considered: $b=-4^\circ$, $-6^\circ$, $-8^\circ$, and $-10^\circ$, along the bulge minor axis as well as outside it, at $l=\pm5^\circ$ and $l=\pm10^\circ$. The bulge X-shaped structure comprises stars that formed in the disk at different locations. Stars formed in the outer disk, beyond the end of the bar, which are part of the boxy peanut-bulge structure may show peaks in the velocity distributions at positive and negative heliocentric radial velocities with high absolute values that can be larger than 100 $\rm km$ $\rm s^{-1}$, depending on the observed direction. In some cases the structure of the velocity field is more complex and several peaks are observed. Stars formed in the inner disk, the most numerous, contribute predominantly to the X-shaped structure and present different kinematic characteristics. Our results may enable us to interpret the cold high-velocity peak observed in the APOGEE commissioning data, as well as the excess of high-velocity stars in the near and far arms of the X-shaped structure at $l$=$0^\circ$ and $b$=$-6^\circ$. When compared with real data, the kinematic picture becomes more complex due to the possible presence in the observed samples of classical bulge and/or thick disk stars. Overall, our results point to the existence of complex patterns and structures in the bulge velocity fields, which are generated by the bar. This suggests that caution should be used when interpreting the bulge kinematics: the presence of substructures, peaks and clumps in the velocity fields is not necessarily a sign of past accretion events.

Radial stellar migration in galactic discs has received much attention in studies of galactic dynamics and chemical evolution, but remains a dynamical phenomenon that needs to be fully quantified. In this work, using a Tree-SPH simulation of an Sb-type disc galaxy, we quantify the effects of blurring (epicyclic excursions) and churning (change of guiding radius). We quantify migration (either blurring or churning) both in terms of flux (the number of migrators passing at a given radius), and by estimating the population of migrators at a given radius at the end of the simulation compared to non-migrators, but also by giving the distance over which the migration is effective at all radii. We confirm that the corotation of the bar is the main source of migrators by churning in a bar-dominated galaxy, its intensity being directly linked to the episode of a strong bar, in the first 1-3 Gyr of the simulation. We show that within the outer Lindblad resonance (OLR), migration is strongly dominated by churning, while blurring gains progressively more importance towards the outer disc and at later times. Most importantly, we show that the OLR limits the exchange of angular momentum, separating the disc in two distinct parts with minimal or null exchange, except in the transition zone, which is delimited by the position of the OLR at the epoch of the formation of the bar, and at the final epoch. We discuss the consequences of these findings for our understanding of the structure of the Milky Way disc. Because the Sun is situated slightly outside the OLR, we suggest that the solar vicinity may have experienced very limited churning from the inner disc.

We analyse the chemical properties of a set of solar vicinity stars, and show that the small dispersion in abundances of \alpha-elements at all ages provides evidence that the SFH has been uniform throughout the thick disk. In the context of long time scale infall models, we suggest that this result points either to a limited dependence of the gas accretion on the Galactic radius in the inner disk (R<10 kpc), or to a decoupling of the accretion history and star formation history due to other processes governing the ISM in the early disk, suggesting that infall cannot be a determining parameter of the chemical evolution at these epochs. We argue however that these results and other recent observational constraints -- namely the lack of radial metallicity gradient and the non-evolving scale length of the thick disk -- are better explained if the early disk is viewed as a pre-assembled gaseous system, with most of the gas settled before significant star formation took place -- formally the equivalent of a closed-box model. In any case, these results point to a weak, or non-existent inside-out formation history in the thick disk, or in the first 3-5 Gyr of the formation of the Galaxy. We argue however that the growing importance of an external disk whose chemical properties are distinct from those of the inner disk would give the impression of an inside-out growth process when seen through snapshots at different epochs. However, the progressive, continuous process usually invoked may not have actually existed in the Milky Way.

The apparent correlation between the specific star formation rate (sSFR) and total stellar mass (M_star) of galaxies is a fundamental relationship indicating how they formed their stellar populations. To attempt to understand this relation, we hypothesize that the relation and its evolution is regulated by the increase in the stellar and gas mass surface density in galaxies with redshift, which is itself governed by the angular momentum of the accreted gas, the amount of available gas, and by self-regulation of star formation. With our model, we can reproduce the specific SFR-M_star relations at z~1-2 by assuming gas fractions and gas mass surface densities similar to those observed for z=1-2 galaxies. We further argue that it is the increasing angular momentum with cosmic time that causes a decrease in the surface density of accreted gas. The gas mass surface densities in galaxies are controlled by the centrifugal support (i.e., angular momentum), and the sSFR is predicted to increase as, sSFR(z)=(1+z)^3/t_H0, as observed (where t_H0 is the Hubble time and no free parameters are necessary). At z>~2, we argue that star formation is self-regulated by high pressures generated by the intense star formation itself. The star formation intensity must be high enough to either balance the hydrostatic pressure (a rather extreme assumption) or to generate high turbulent pressure in the molecular medium which maintains galaxies near the line of instability (i.e. Toomre Q~1). The most important factor is the increase in stellar and gas mass surface density with redshift, which allows distant galaxies to maintain high levels of sSFR. Without a strong feedback from massive stars, such galaxies would likely reach very high sSFR levels, have high star formation efficiencies, and because strong feedback drives outflows, ultimately have an excess of stellar baryons (abridged).

We develop a chemical evolution model in order to study the star formation history of the Milky Way. Our model assumes that the Milky Way is formed from a closed box-like system in the inner regions, while the outer parts of the disc experience some accretion. Unlike the usual procedure, we do not fix the star formation prescription (e.g. Kennicutt law) in order to reproduce the chemical abundance trends. Instead, we fit the abundance trends with age in order to recover the star formation history of the Galaxy. Our method enables one to recover with unprecedented accuracy the star formation history of the Milky Way in the first Gyrs, in both the inner (R<7-8kpc) and outer (R>9-10kpc) discs as sampled in the solar vicinity. We show that, in the inner disc, half of the stellar mass formed during the thick disc phase, in the first 4-5 Gyr. This phase was followed by a significant dip in the star formation activity (at 8-9 Gyr) and a period of roughly constant lower level star formation for the remaining 8 Gyr. The thick disc phase has produced as many metals in 4 Gyr as the thin disc in the remaining 8 Gyr. Our results suggest that a closed box model is able to fit all the available constraints in the inner disc. A closed box system is qualitatively equivalent to a regime where the accretion rate, at high redshift, maintains a high gas fraction in the inner disc. In such conditions, the SFR is mainly governed by the high turbulence of the ISM. By z~1 it is possible that most of the accretion takes place in the outer disc, while the star formation activity in the inner disc is mostly sustained by the gas not consumed during the thick disc phase, and the continuous ejecta from earlier generations of stars. The outer disc follows a star formation history very similar to that of the inner disc, although initiated at z~2, about 2 Gyr before the onset of the thin disc formation in the inner disc.

By analyzing a N-body simulation of a bulge formed simply via a bar instability mechanism operating on a kinematically cold stellar disk, and by comparing the results of this analysis with the structural and kinematic properties of the main stellar populations of the Milky Way bulge, we conclude that the bulge of our Galaxy is not a pure stellar bar formed from a pre-existing thin stellar disk, as some studies have recently suggested. On the basis of several arguments emphasized in this paper, we propose that the bulge population which, in the Milky Way, is observed not to be part of the peanut structure corresponds to the old galactic thick disk, thus implying that the Milky Way is a pure thin+thick disk galaxy, with only a possible limited contribution of a classical bulge.

We compare the star-formation history and dynamics of the Milky Way (MW) with the properties of distant disk galaxies. During the first ~4 Gyr of its evolution, the MW formed stars with a high star-formation intensity (SFI), Sigma_SFR~0.6 Msun/yr/kpc2 and as a result, generated outflows and high turbulence in its interstellar medium. This intense phase of star formation corresponds to the formation of the thick disk. The formation of the thick disk is a crucial phase which enables the MW to have formed approximately half of its total stellar mass by z~1 which is similar to "MW progenitor galaxies" selected by abundance matching. This agreement suggests that the formation of the thick disk may be a generic evolutionary phase in disk galaxies. Using a simple energy injection-kinetic energy relationship between the 1-D velocity dispersion and SFI, we can reproduce the average perpendicular dispersion in stellar velocities of the MW with age. This relationship, its inferred evolution, and required efficiency are consistent with observations of galaxies from z~0-3. The high turbulence generated by intense star formation naturally resulted in a thick disk, a chemically well-mixed ISM, and is the mechanism that links the evolution of MW to the observed characteristics of distant disk galaxies.

By means of idealized, dissipationless N-body simulations which follow the formation and subsequent buckling of a stellar bar, we study the characteristics of boxy/peanut-shaped bulges and compare them with the properties of the stellar populations in the Milky Way bulge. The main results of our modeling, valid for the general family of boxy/peanut shaped bulges, are the following: (i) because of the spatial redistribution in the disk initiated at the epoch of bar formation, stars from the innermost regions to the outer Lindblad resonance of the stellar bar are mapped into a boxy bulge; (ii) the contribution of stars to the local bulge density depends on their birth radius: stars born in the innermost disk tend to dominate the innermost regions of the boxy bulge, while stars originating closer to the OLR are preferentially found in the outer regions of the boxy/peanut structure; (iii) stellar birth radii are imprinted in the bulge kinematics, the larger the birth radii of stars ending up in the bulge, the greater their rotational support and the higher their line-of- sight velocity dispersions (but note that this last trend depends on the bar viewing angle); (iv) the higher the classical bulge-over-disk ratio, the larger its fractional contribution of stars at large vertical distance from the galaxy mid-plane. (ABRIDGED) On the basis of their chemical and kinematic characteristics, the results of our modeling suggests that the populations A, B and C, as defined by the ARGOS survey, can be associated, respectively, with the inner thin disk, to the young thick and to the old thick disk, following the nomenclature recently suggested for stars in the solar neighborhood by Haywood et al. (2013).

The understanding of Galaxy evolution can be facilitated by the use of population synthesis models, which allows us to test hypotheses on the star formation history, star evolution, and chemical and dynamical evolution of the Galaxy. The new version of the Besancon Galaxy model (hereafter BGM) aims to provide a more flexible and powerful tool to investigate the initial mass function (IMF) and star formation rate (SFR) of the Galactic disc. We present a new strategy for the generation of thin disc stars, which assumes the IMF, SFR and evolutionary tracks as free parameters. We have updated most of the ingredients for the star count production and, for the first time, binary stars are generated in a consistent way. The local dynamical self-consistency is maintained in this new scheme. We then compare simulations from the new model with Tycho-2 data and the local luminosity function, as a first test to verify and constrain the new ingredients. The effects of changing thirteen different ingredients of the model are systematically studied. For the first time, a full sky comparison is performed between BGM and data. This strategy allows us to constrain the IMF slope at high masses, which is found to be close to 3.0 and excludes a shallower slope such as Salpeter's one. The SFR is found decreasing whatever IMF is assumed. The model is compatible with a local dark matter density of 0.011 Mo/pc^3 implying that there is no compelling evidence for the significant amount of dark matter in the disc. While the model is fitted to Tycho-2 data, which is a magnitude limited sample with V<11, we check that it is still consistent with fainter stars. The new model constitutes a new basis for further comparisons with large scale surveys and is being prepared to become a powerful tool for the analysis of the Gaia mission data.

Constraints on the Galactic bulge/bar structure and formation history from stellar kinematics and metallicities mainly come from relatively high-latitude fields (|b|>4) where a complex mix of stellar population is seen. We aim here to constrain the formation history of the Galactic bar by studying the radial velocity and metallicity distributions of stars in-situ (|b|<1). We observed red clump stars in four fields along the bar's major axis (l=10,6,-6 and b=0 plus a field at l=0,b=1) with low-resolution spectroscopy from VLT/FLAMES, observing around the CaII triplet. We developed robust methods for extracting radial velocity and metallicity estimates from these low signal-to-noise spectra. We derived distance probability distributions using Bayesian methods rigorously handling the extinction law. We present radial velocities and metallicity distributions, as well as radial velocity trends with distance. We observe an increase in the radial velocity dispersion near the Galactic plane. We detect the streaming motion of the stars induced by the bar in fields at l=+/-6, the highest velocity components of this bar stream being metal-rich ([Fe/H]~0.2 dex). Our data is consistent with a bar inclined at 26+/-3 from the Sun-Galactic centre line. We observe a significant fraction of metal-poor stars, in particular in the field at l=0,b=1. We confirm the flattening of the metallicity gradient along the minor axis when getting closer to the plane, with a hint that it could actually be inverted. Our stellar kinematics corresponds to the expected behaviour of a bar issued from the secular evolution of the Galactic disc. The mix of several populations, seen further away from the plane, is also seen in the bar in-situ since our metallicity distributions highlight a different spatial distribution between metal-poor and metal-rich stars, the more metal-poor stars being more centrally concentrated.

We report the first robust measurement of the Milky Way star formation history using the imprint left on chemical abundances of long-lived stars. The formation of the Galactic thick disc occurs during an intense star formation phase between 9.0 (z~1.5) and 12.5 Gyr (z~4.5) ago and is followed by a dip (at z~1.1) lasting about 1 Gyr. Our results imply that the thick disc is as massive as the Milky Way's thin disc, suggesting a fundamental role of this component in the genesis of our Galaxy, something that had been largely unrecognized. This new picture implies that huge quantities of gas necessary to feed the building of the thick disc must have been present at these epochs, in contradiction with the long-term infall assumed by chemical evolution models in the last two decades. These results allow us to fit the Milky Way within the emerging features of the evolution of disc galaxies in the early Universe.

Standard chemical evolution models based on long-term infall are affected by a number of problems, evidenced by the analysis of the most recent data. Among these: (1) models rely on the local metallicity distribution, assuming its shape is valid for the entire Galaxy, which it is not; (2) they assume that the solar vicinity abundance patterns resulted from a unique chemical evolution, which it does not; (3) they assume the disk is a single structure with chemical properties that are a smooth function of the distance to the galactic center, which it is not. Moreover, new results point to a thick disk being as massive as the thin disk, leading to a change of paradigm in the way we see the formation of these structures. I discuss these various issues, and, commenting on Snaith et al. (2014), how a closed box model offers an interesting approximation to the galactic chemical evolution, by providing the conditions in which large amounts of gas are available in the disk at high redshift. The novel way presented in Snaith et al. (2014) to derive SFH from stellar abundances is also discussed, providing a measurement of the SFH of old populations that is valid for the entire Galaxy. The derived SFH shows that the formation of the thick disk has been the dominant epoch of star formation in our Galaxy.

We analyze high quality abundances data of solar neighborhood stars and show that there are two distinct regimes of [alpha/Fe] versus age which we identify as the epochs of the thick and thin disk formation. A tight correlation between metallicity and [alpha/Fe] versus age is clearly identifiable on thick disk stars, implying that this population formed from a well mixed ISM, over a time scale of 4-5 Gyr. Thick disk stars vertical velocity dispersion correlate with age, with the youngest objects having as small scale heights as those of thin disk stars. A natural consequence of these two results is that a vertical metallicity gradient is expected in this population. We suggest that the thick disk set the initial conditions for the formation of the inner thin disk. This provides also an explanation of the apparent coincidence between the step in metallicity at 7-10 kpc in the thin disk and the confinment of the thick disk at about R<10 kpc. We suggest that the outer thin disk developped outside the influence of the thick disk, but also that the high alpha-enrichment of the outer regions may originate from a primordial pollution by the gas expelled from the thick disk. Local metal-poor thin disk stars, whose properties are best explained by an origin in the outer disk, are shown to be as old as the youngest thick disk (9-10 Gyr), implying that the outer thin disk started to form while the thick disk formation was still on-going in the inner Galaxy. We point out that, given the tight age-abundance relations in the thick disk, an inside-out process would give rise to a radial gradient in abundances in this population which is not observed. Finally, we argue that the data discussed here leave little room for radial migration, either to have contaminated the solar vicinity, or, to have redistributed stars in significant proportion across the solar annulus.

In this paper we report a study (see Mastrobuono-Battisti et al., 2012) about the formation and characteristics of the tidal tails around Palomar 5 along its orbit in the Milky Way potential, by means of direct N-body simulations and simplified numerical models. Unlike previous findings, we are able to reproduce the substructures observed in the stellar streams of this cluster, without including any lumpiness in the dark matter halo. We show that overdensities similar to those observed in Palomar 5 can be reproduced by the epicyclic motion of stars along its tails, i.e. a simple local accumulation of orbits of stars that escaped from the cluster with very similar positions and velocities. This process is able to form stellar clumps at distances of several kiloparsecs from the cluster, so it is not a phenomenon confined to the inner part of Palomar 5's tails, as previously suggested.

By means of N-body simulations, we show that radial migration in galaxy disks, induced by bar and spiral arms, leads to significant azimuthal variations in the metallicity distribution of old stars at a given distance from the galaxy center. Metals do not show an axisymmetric distribution during phases of strong migration. Azimuthal variations are visible during the whole phase of strong bar phase, and tend to disappear as the effect of radial migration diminishes, together with a reduction in the bar strength. These results suggest that the presence of inhomogeneities in the metallicity distribution of old stars in a galaxy disk can be a probe of ongoing strong migration. Such signatures may be detected in the Milky Way by Gaia (and complementary spectroscopic data), as well as in external galaxies, by IFU surveys like CALIFA and ATLAS3D. Mixing - defined as the tendency toward a homogeneous, azimuthally symmetric, stellar distribution in the disk - and migration turns out to be two distinct processes, the effects of mixing starting to be visible when strong migration is over.

By means of direct N-body simulations and simplified numerical models, we study the formation and characteristics of the tidal tails around Palomar 5, along its orbit in the Milky Way potential. Unlike previous findings, we are able to reproduce the substructures observed in the stellar streams of this cluster, without including any lumpiness in the dark matter halo. We show that overdensities similar to those observed in Palomar 5 can be reproduced by the epicyclic motion of stars along its tails, i.e. a simple local accumulation of orbits of stars that escaped from the cluster with very similar positions and velocities. This process is able to form stellar clumps at distances of several kiloparsecs from the cluster, so it is not a phenomenon confined to the inner part of Palomar 5's tails, as previously suggested. Our models can reproduce the density contrast between the clumps and the surrounding tails found in the observed streams, without including any lumpiness in the dark halo, suggesting new upper limits on its granularity.

We investigate the characteristics of the thick disk in the Canada - France - Hawaii - Telescope Legacy Survey (CFHTLS) fields, complemented at bright magnitudes with Sloan Digital Sky Survey (SDSS) data. The ([Fe/H], Z) distributions are derived in the W1 and W3 fields, and compared with simulated maps produced using the Besan\c{c}on model. It is shown that the thick disk, represented in star-count models by a distinct component, is not an adequate description of the observed ([Fe/H], Z) distributions in these fields.

We have chosen the name of GYES, one of the mythological giants with one hundred arms, offspring of Gaia and Uranus, for our instrument study of a multifibre spectrograph for the prime focus of the Canada-France-Hawaii Telescope. Such an instrument could provide an excellent ground-based complement for the Gaia mission and a northern complement to the HERMES project on the AAT. The CFHT is well known for providing a stable prime focus environment, with a large field of view, which has hosted several imaging instruments, but has never hosted a multifibre spectrograph. Building upon the experience gained at GEPI with FLAMES-Giraffe and X-Shooter, we are investigating the feasibility of a high multiplex spectrograph (about 500 fibres) over a field of view 1 degree in diameter. We are investigating an instrument with resolution in the range 15000 to 30000, which should provide accurate chemical abundances for stars down to 16th magnitude and radial velocities, accurate to 1 km/s for fainter stars. The study is led by GEPI-Observatoire de Paris with a contribution from Oxford for the study of the positioner. The financing for the study comes from INSU CSAA and Observatoire de Paris. The conceptual study will be delivered to CFHT for review by October 1st 2010.

We use an extended and homogeneous data set of Galactic planetary nebulae (PNe) to study the metallicity gradients and the Galactic structure and evolution. The most up-to-date abundances, distances (calibrated with Magellanic Cloud PNe) have been employed, together with a novel homogeneous morphological classification, to characterize the different PN populations. We confirm that morphological classes have a strong correlation with PN Peimbert's Type, and also with their distribution on the Galactic landscape. We studied the alpha-element distribution within the Galactic disk, and found that the best selected disk population, together with the most reliable PN distance scale yields to a radial oxygen gradient of d[log(O/H)]/dR=-0.023 +- 0.006 dex/ kpc for the whole disk sample, and of d[log(O/H)]/dR= -0.035+-0.024, -0.023+-0.005, and -0.011+-0.013 dex/kpc respectively for Type I, II, and III PNe. Neon gradients for the same PN types confirm the trend. Accurate statistical analysis show moderately high uncertainties in the slopes, but also confirm the trend of steeper gradient for PNe with more massive progenitors, indicating a possible steepening with time of the Galactic disk metallicity gradient. The PN metallicity gradients presented here are consistent with the local metallicity distribution; furthermore, oxygen gradients determined with young and intermediate age PNe show good consistency with oxygen gradients derived respectively from other young (OB stars, HII regions) and intermediate (open cluster) Galactic populations. We also extend the Galactic metallicity gradient comparison by revisiting the open cluster [Fe/H] data from high resolution spectroscopy. The analysis suggests that they could be compliant with the same general picture of a steepening of gradient with time.

The correlation between stellar metallicity and the presence of giant planets is well established. It has been tentatively explained by the possible increase of planet formation probability in stellar disks with enhanced amount of metals. However, there are two caveats to this explanation. First, giant stars with planets do not show a metallicity distribution skewed towards metal-rich objects, as found for dwarfs. Second, the correlation with metallicity is not valid at intermediate metallicities, for which it can be shown that giant planets are preferentially found orbiting thick disk stars. None of these two peculiarities is explained by the proposed scenarios of giant planet formation. We contend that they are galactic in nature, and probably not linked to the formation process of giant planets. It is suggested that the same dynamical effect, namely the migration of stars in the galactic disk, is at the origin of both features, with the important consequence that most metal-rich stars hosting giant planets originate from the inner disk, a property that has been largely neglected until now. We illustrate that a planet-metallicity correlation similar to the observed one is easily obtained if stars from the inner disk have a higher percentage of giant planets than stars born at the solar radius, with no specific dependence on metallicity. We propose that the density of molecular hydrogen in the inner galactic disk (the molecular ring) could play a role in setting the high percentage of giant planets that originate from this region.

Using test particle simulations we examine the structure of the outer Galactic disk as it is perturbed by a satellite in a tight eccentric orbit about the Galaxy. A satellite of mass a few times 10^9 Msol can heat the outer Galactic disk, excite spiral structure and a warp and induce streams in the velocity distribution. We examine particle eccentricity versus the change in mean radius between initial and current orbits. Correlations between these quantities are reduced after a few satellite pericenter passages. Stars born in the outer galaxy can be moved in radius from their birth positions and be placed in low eccentricity orbits inside their birth radii. We propose that mergers and perturbations from satellite galaxies and subhalos can induce radial mixing in the stellar metallicity distribution.

The analysis of the kinematics of solar neighbourhood stars shows that the low and high metallicity tails of the thin disc are populated by objects which orbital properties suggest an origin in the outer and inner galactic disc, respectively. Signatures of radial migration are identified in various recent samples, and are shown to be responsible for the high metallicity dispersion in the age-metallicity distribution. Most importantly, it is shown that the population of low metallicity wanderers of the thin disc (-0.7<[Fe/H]<-0.3 dex) is also responsible for the apparent hiatus in metallicity with the thick disc (which terminal metallicity is about -0.2 dex). It implies that the thin disc at the solar circle has started to form stars at about this same metallicity. This is also consistent with the fact that 'transition' objects, which have alpha-element abundance intermediate between that of the thick and thin discs, are found in the range [-0.4,-0.2] dex. Once the metal-poor thin disc stars are recognised for what they are - wanderers from the outer thin disc - the parenthood between the two discs can be identified on stars genuinely formed at the solar circle through an evolutionary sequence in [alpha/Fe] and [Fe/H] . Another consequence is that stars that can be considered as truly resulting of the chemical evolution at the solar circle have a metallicity restricted to about [-0.2,+0.2] dex, confirming an old idea that most chemical evolution in the Milky Way have preceded the thin disc formation.

Stars with planets at intermediate metallicities ([-0.7,-0.2] dex) exhibit properties that differ from the general field stars. Thirteen stars with planets reported in this metallicity range belong to the thick disc, while only one planet have been detected among stars of the thin disc. Although this statistics is weak, it contradicts the known correlation between the presence of planet and metallicity. We relate this finding to the specific property of the thin disc in this metallicity range, where stars are shown to rotate around the Galaxy faster than the Sun. Their orbital parameters are conveniently explained if they are contaminants coming from the outer Galactic disc, as a result of radial mixing. This must be considered together with the fact that metal-rich stars ([Fe/H]>+0.1 dex) found in the solar neighbourhood, which are the hosts of most of the detected planets, are suspected of being wanderers from the inner Galactic disc. It is then questionned why stars that originate in the inner and outer thin disc show respectively the highest and lowest rate of detected planets. It is suggested that the presence of giant planets might be primarily a function of a parameter linked to galactocentric radius, but not metallicity. Combined with the existing radial metallicity gradient, then radial mixing explains the correlation at high metallicity observed locally, but also the peculiarity found at low metallicity, which cannot be accounted for by a simple correlation between metallicity and planet probability.

I review the uncertainties in two observational constraints of the Galactic disc chemical evolution: the metallicity distribution (MD) of long-lived dwarfs and the age-metallicity relation (AMR). It is shown first that the observed MD, designed with standard methods, is more fit to a closed-box model than to the infall MD. We argue that this is due to the specific contribution of the thick disc population, which has been overlooked both in the derivation of the observed MD and in the standard chemical evolution models. Although this agreement disqualifies the MD as the best supportive evidence for infall, we argue that the evolution must be more complex than described by either the closed-box or the standard infall models. A new age-metallicity distribution is obtained, where particularities of the previous recent determinations are phased out. The new AMR shows a mean increase limited to about a factor of 2 in Z over the disc age. It is shown that below 3 Gyrs, the dispersion in metallicity is about 0.1 dex, which, (given the observational uncertainties), is compatible with the small cosmic dispersion measured on the interstellar medium and meteoritic pre-solar dust grains. A population that is progressively older and more metal rich arises at a metallicity greater than that of the Hyades, to reach [Fe/H] ~ +0.5 dex at ages > 5 Gyr. We suggest that this is best explained by radial migration. A symmetrical widening towards lower metallicities is seen at about the same age. Finally, the new derived ages are sufficiently consistent that an AMR within the thick disc is confirmed. These new features altogether draw a picture of the chemical evolution where dynamical effects and complexity in the AMR dominate, rather than a generalized high dispersion at all ages. (abridged)

The Gaia mission is designed as a Galaxy explorer, and will measure simultaneously, in a survey mode, the five or six phase space parameters of all stars brighter than 20th magnitude, as well as providing a description of their astrophysical characteristics. These measurements are obtained by combining an astrometric instrument with micro-arcsecond capabilities, a photometric system giving the magnitudes and colours in 15 bands and a medium resolution spectrograph named the Radial Velocity Spectrometer (RVS). The latter instrument will produce spectra in the 848 to 874 nm wavelength range, with a resolving power R = 11 500, from which radial velocities, rotational velocities, atmospheric parameters and abundances can be derived. A companion paper (Katz et al. 2004) presents the characteristics of the RVS and its performance. This paper details the outstanding scientific impact of this important part of the Gaia satellite on some key open questions in present day astrophysics. The unbiased and simultaneous acquisition of multi-epoch radial velocities and individual abundances of key elements in parallel with the astrometric parameters is essential for the determination of the dynamical state and formation history of our Galaxy. Moreover, for stars brighter than V=15, the resolving power of the RVS will give information about most of the effects which influence the position of a star in the Hertzsprung-Russell diagram, placing unprecedented constraints on the age, internal structure and evolution of stars of all types. Finally, the RVS multi-epoch observations are ideally suited to the identification, classification and characterisation of the many types of double, multiple and variable stars.

The definition and optimisation studies for the Gaia satellite spectrograph, the Radial Velocity Spectrometer (RVS), converged in late 2002 with the adoption of the instrument baseline. This paper reviews the characteristics of the selected configuration and presents its expected performance. The RVS is a 2.0 by 1.6 degree integral field spectrograph, dispersing the light of all sources entering its field of view with a resolving power R=11 500 over the wavelength range [848, 874] nm. The RVS will continuously and repeatedly scan the sky during the 5 years of the Gaia mission. On average, each source will be observed 102 times over this period. The RVS will collect the spectra of about 100-150 million stars up to magnitude V~17-18. At the end of the mission, the RVS will provide radial velocities with precisions of ~2 km/s at V=15 and \~15-20 km/s at V=17, for a solar metallicity G5 dwarf. The RVS will also provide rotational velocities, with precisions (at the end of the mission) for late type stars of sigma_vsini ~5 km/s at V~15 as well as atmospheric parameters up to V~14-15. The individual abundances of elements such as Silicon and Magnesium, vital for the understanding of Galactic evolution, will be obtained up to V~12-13. Finally, the presence of the 862.0 nm Diffuse Interstellar Band (DIB) in the RVS wavelength range will make it possible to derive the three dimensional structure of the interstellar reddening.

The mean metallicity of the Milky Way thin disc in the solar neighbourhood is still a matter of debate, and has recently been subject to upward revision (Haywood, 2001). Our star sample was drawn from a set of solar neighbourhood dwarfs with photometric metallicities. In a recent study, Reid (2002) suggests that our metallicity calibration, based on Geneva photometry, is biased. We show here that the effect detected by Reid is not a consequence of our adopted metallicity scale, and we confirm that our findings are robust. On the contrary, the application to Stromgren photometry of the Schuster & Nissen metallicity scale is problematic. Systematic discrepancies of about 0.1 to 0.3 dex affect the photometric metallicity determination of metal rich stars, on the colour interval 0.22< b-y <0.59, i.e including F and G stars. For F stars, it is shown that this is a consequence of a mismatch between the standard sequence m_1(b-y) of the Hyades used by Schuster & Nissen to calibrate their metallicity scale, and the system of Olsen (1993, 1994ab). It means that although Schuster & Nissen calibration and Olsen photometry are intrinsically correct, there are mutually incompatible for metal rich, F-type stars. For G stars, the discrepancy is most probably the continuation of the same problem, albeit worthen by the lack of spectroscopic calibrating stars. A corrected calibration is proposed which renders the calibration of Schuster & Nissen applicable to the catalogues of Olsen. We also give a simpler calibration referenced to the Hyades sequence, valid over the same color and metallicity ranges.

We present a revised metallicity distribution of dwarfs in the solar neighbourhood. This distribution is centred on solar metallicity. We show that previous metallicity distributions, selected on the basis of spectral type, are biased against stars with solar metallicity or higher. A selection of G-dwarf stars is inherently biased against metal rich stars and is not representative of the solar neighbourhood metallicity distribution. Using a sample selected on colour, we obtain a distribution where approximately half the stars in the solar neighbourhood has a metallicity higher than [Fe/H]=0. The percentage of mid-metal-poor stars ([Fe/H]$<$-0.5) is approximately 4 per cent, in agreement with present estimates of the thick disc. In order to have a metallicity distribution comparable to chemical evolution model predictions, we convert the star fraction to mass fraction, and show that another bias against metal-rich stars affects dwarf metallicity distributions, due to the colour (or spectral type) limits of the samples. Reconsidering the corrections due to the increasing thickness of the stellar disc with age, we show that the Simple Closed-Box model with no instantaneous recycling approximation gives a reasonable fit to the observed distribution. Comparisons with the age-metallicity relation and abundance ratios suggest that the Simple Closed-Box model may be a viable model of the chemical evolution of the Galaxy at solar radius.

Accurate characterization of thick disc properties from recent kinematic and photometric surveys provides converging evidences that this intermediate population is a sequel of the violent heating of early disc populations by a merging satellite galaxy. The thick disc population is revisited under the light of new data in a number of galactic sample fields. Various thick disc hypotheses are fitted to observational data through a maximum likelihood technique. The resulting characteristics of the thick disc are the following : a scale height of 760 $\pm$ 50 pc, with a local density of 5.6 $\pm$ 1 \% of the thin disc. The scale length is constrained to be 2.8 $\pm$ 0.8 kpc, well in agreement with the disc scale length (2.5 $\pm$ 0.3 kpc). The mean metallicity of the thick disc is found to be -0.7 $\pm$ 0.2 dex, with no {significant} metallicity gradients. These photometric constraints in combination with kinematic data give new constraints on the thick disc formation. We show that thick disc characteristics are {hardly compatible with a top-down formation scenario} but fully compatible with a violent merging event arising at the early thin disc life time as described by Quinn, Hernquist \& Fullagar (1993).

We have been doing a photometric and astrometric sample survey programme on Schmidt plates in a few selected directions in the Galaxy with an aim to study galactic stellar populations. In order to calibrate Schmidt plates photometrically, it is necessary to have a number of photometric standards which cover the entire range of magnitudes to be studied on the plates. As Schmidt plates cover a wide area of the sky, the standards should also be spread over the surface of the plate so as to minimise geometrical effects present on the plate. This paper is first in a series of papers aimed at to provide photometric standards in the selected galactic directions to be studied. In this paper we present photometric magnitudes obtained through CCD observation in the galactic direction (l=170deg, b=45deg) centered at alpha(2000) = 9h 41m 20s and delta(2000) = +49deg 54' 20'' and covering an area of 14 sq degs. B and V magnitudes have been obtained for a total of 214 stars out of which for 73 stars U magnitudes also have been obtained. The magnitude range covered is 11<V<21 and 0.18<B-V<1.85.