Search

WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, will see first light in late 2022. WEAVE comprises a new 2-degree field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable 'mini' integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrograph covering the wavelength range 366$-$959\,nm at $R\sim5000$, or two shorter ranges at $R\sim20\,000$. After summarising the design and implementation of WEAVE and its data systems, we present the organisation, science drivers and design of a five- to seven-year programme of eight individual surveys to: (i) study our Galaxy's origins by completing Gaia's phase-space information, providing metallicities to its limiting magnitude for $\sim$3 million stars and detailed abundances for $\sim1.5$ million brighter field and open-cluster stars; (ii) survey $\sim0.4$ million Galactic-plane OBA stars, young stellar objects and nearby gas to understand the evolution of young stars and their environments; (iii) perform an extensive spectral survey of white dwarfs; (iv) survey $\sim400$ neutral-hydrogen-selected galaxies with the IFUs; (v) study properties and kinematics of stellar populations and ionised gas in $z1$ million spectra of LOFAR-selected radio sources; (viii) trace structures using intergalactic/circumgalactic gas at $z>2$. Finally, we describe the WEAVE Operational Rehearsals using the WEAVE Simulator., 41 pages, 27 figures, accepted for publication by MNRAS

We present the e-TidalGCs Project which aims at modeling and predicting the extra-tidal features surrounding all Galactic globular clusters for which 6D phase space information, masses and sizes are available (currently 159 globular clusters). We focus the analysis and presentation of the results on the distribution of extra-tidal material on the sky, and on the different structures found at different heliocentric distances. We emphasize the wide variety of morphologies found: beyond the canonical tidal tails, our models reveal that the extra-tidal features generated by globular clusters take a wide variety of shapes, from thin and elongated shapes, to thick, and complex halo-like structures. We also compare some of the most well studied stellar streams found around Galactic globular clusters to our model predictions, namely those associated to the clusters NGC 3201, NGC 4590, NGC 5466 and Pal 5. Additionally, we investigate how the distribution and extension in the sky of the simulated streams vary with the Galactic potential by making use of three different models, containing or not a central spheroid, or a stellar bar. Overall, our models predict that the mass lost by the current globular cluster population in the field from the last 5 Gyrs is between $0.3-2.1\times10^{7}M_{\odot}$, an amount comparable between 7-55 % of current mass. Most of this lost mass is found in the inner Galaxy, with the half-mass radius of this population being between 4-6 kpc. The outputs of the simulations will be publicly available, at a time when the ESA Gaia mission and complementary spectroscopic surveys are delivering exquisite data to which these models can be compared., Comment: 51 pages, 34 figures, Accepted for publication on A&A

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., Comment: 22 pages, 19 figures, accepted for publication in MNRAS

Understanding our Galaxy’s structure, formation, and evolution will, over the next decades, continue to benefit from the wonderful large survey by Gaia, for astrometric, kinematic, and spectroscopic characterization, and by large spectroscopic surveys for chemical characterization. The weak link for full exploitation of these data is age characterization, and stellar age estimation relies predominantly on mass estimates. The ideas presented in this White Paper shows that a seismology survey is the way out of this situation and a natural complement to existing and planned surveys. These ideas are strongly rooted in the past decade’s experience of the so-called Seismology revolution, initiated with CoRoT and Kepler. The case of red giant stars is used here as the best current illustration of what we can expect from seismology for large samples, but premises for similar developments exist in various other classes of stars covering other ranges of age or mass. Whatever the star considered, the first information provided by stellar pulsations is always related to the mean density and thus to the mass (and age). In order to satisfy the need for long-duration and allsky coverage, we rely on a new instrumental concept which decouples integration time and sampling time. We thus propose a long (~1 year) all-sky survey which would perfectly fit between TESS, PLATO, and the Rubin Observatory (previously known as LSST) surveys to offer a time domain complement to the current and planned astrometric and spectroscopic surveys. The fine characterization of host stars is also a key aspect for the interpretation and exploitation of the various projects – anticipated in the framework of the Voyage 2050 programme – searching for atmospheric characterization of terrestrial planets or, more specifically, looking for a signature of life, in distant planets.

The Galactic globular cluster system went and is still going through dynamical processes that require to be explored in detail. Here we illustrate how primordial massive globular clusters born in the Milky Way’s disc evolved by stripping material from each other or even merging very early during their lives. These processes might explain the puzzling presence of star-by-star spreads in iron content observed in massive globular clusters and should be taken into account when studying globular cluster stellar populations. In this context, we show how the direct comparison between the predictions provided by our direct N-body simulations and observations can shed light on the origin and chemo-dynamical evolution of globular clusters.

Context. Recent observations of the Milky Way and galaxies at high redshifts suggest that galaxy discs were already in place soon after the Big Bang. While the gas infall history of the Milky Way in the inner disc has long been assumed to be characterised by a short accretion timescale, this has not been directly constrained using observations. Aims. Using data for the inner regions of the Milky Way recently produced by APOGEE and Gaia and of unprecedented quantity and quality, we aim to derive strong constraints on the infall history of the inner (< 6 kpc) Galaxy (with a focus on stars between 4 and 6 kpc, which we show is an appropriate proxy for the entire inner disc). Methods. We implemented gas infall into a chemical evolution model of the Galaxy disc, and used a Schmidt–Kennicutt law to connect the infall to the star formation. We explore a number of models, and two different formulations of the infall law. In one formulation, the infall is non-parametric, and in the other the infall has an explicitly exponential form. We fit the model parameters to the time–[Si/Fe] distribution of solar vicinity stars, and the metallicity and [Si/Fe] distribution function of stars with a galactocentric radius of between 4 and 6 kpc from APOGEE. Results. Our results point to a fast, early gas accretion, and an upper limit on the accretion timescale of around 2 Gyr in the inner disc of the Milky Way. This suggests that at least half the baryons were in place within 2−3 Gyr of the Big Bang, and that half the stars of the inner disc formed within the first 5 Gyr, during the thick disc formation phase. This implies that the stellar mass of the inner disc is dominated by the thick disc, supporting our previous work, and that the gas accretion onto the inner disc was rapid and early.

As the remnants of stars with initial masses ≲8 M⊙, 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, griz photometry from Pan-STARRS1, as well as proper motions from Gaia DR2, to select 25,156 white dwarf candidates over ∼4500 deg2 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 (MW), 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 MW disk began forming stars (11.3 ± 0.5) Gyr ago, with a peak rate of (8.8 ± 1.4) M⊙ yr −1 at (9.8 ± 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 ± 3)%. Incorporating this methodology with data from future wide-field surveys such as the Large Synoptic Survey Telescope, Euclid, The Cosmological Advanced Survey Telescope for Optical and ultraviolet Research, and the Wide Field Infrared Survey Telescope should provide an unprecedented view into the formation of the MW 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., 25 pages, 24 figures, accepted for publication on A&A

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., 9 pages, 8 figures, A&A in press

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], Comment: 30 pages, 18 figures, A&A, in press

International audience; In the Milky Way bulge, metal-rich stars form a strong bar and are more peanut-shaped than metal-poor stars. It has recently been claimed that this behavior is driven by the initial (i.e., before bar formation) 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 disk is not necessary to explain the characteristics of the Milky Way bulge. We discuss this issue again by analyzing two dissipationless N-body simulations of boxy or peanut-shaped bulges formed from composite stellar disks that consist of kinematically cold and hot stellar populations. These two models represent two extreme cases: one where all three components of the disk have a fixed vertical velocity dispersion and different in-plane radial dispersion, and another where they all have a fixed radial dispersion and different vertical random motions (thickness). This is intended to quantify the drivers of the main features that are observed in composite boxy or peanut-shaped bulges and their origin. We quantify the mapping into a boxy or peanut-shaped bulge of disk populations in these two cases, and we conclude that initial vertical random motions are as important as in-plane random motions in determining the relative contribution of cold-and hotdisk 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, significant differences exist in the morphology and strength of the resulting boxy or peanut-shaped bulges. In particular, the model where disk populations initially have only different in-plane random motions, but similar thickness, results in a boxy or peanut-shaped bulge where all populations have a similar peanut shape, independent of their initial kinematics or metallicity. This is at odds with the trends observed in the Milky Way bulge. 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 disk stars outside the bulge region. As a consequence of this analysis, we conclude that given our current knowledge of the Milky Way bulge and of the properties of its main stellar components, a metal-poor, kinematically (radial and vertical) hot component, that is, a thick disk, is necessary in the Milky Way before bar formation. This supports the scenario that has been traced in previous works. Boxy or peanut-shaped bulges and their surrounding regions are fossil records of the conditions present at early times in disk galaxies, and by dissecting their stellar components by chemical compositions and/or age, it may be possible to reconstruct their early state.

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., Comment: 17 pages, accepted for publication in MNRAS

We present a high-resolution numerical study of the phase-space diversity in an isolated Milky Way-type galaxy. Using a single N-body simulation (N ≈ 0.14 × 109) 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 Data Release 2 (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. Thus, 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 outward and produce axisymmetric variations of the mean vertical coordinate and vertical velocity component of about 100 − 200 pc and 1 − 2 km s−1, respectively. As a consequence, the phase-space wrapping results in the formation of patterns with various morphologies 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. Such vertical oscillations trigger the formation of various time-dependent phase-space spirals in the entire disk. 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 neighborhood 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 to the time of the onset of the perturbation: phase-space spirals may have been caused by perturbations originated several gigayearrs ago, and not as recent as suggested so far.

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., Comment: 13 pages, 12 figures, 3 tables. Accepted for publication in A&A

0, info:eu-repo/semantics/published

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., 10 pages, 6 figures, accepted for publication in A&A Letters

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., Comment: 15 pages, 14 figures. Accepted for publication in A&A

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., MNRAS, in press (accepted March 27th, 2018)

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., minor corrections, 10 pages, 12 figures, accepted in ApJ

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] $, Comment: 7 pages, 4 figures, accepted for publication in A&A

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)., Comment: 26 pages, 14 figures, accepted in A&A

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., Comment: 17 pages, 17 figure, accepted for publication in A&A

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., Comment: 13 pages, 9 figures, 6 embedded animations (for Acrobat only); high resolution movies available at https://goo.gl/YC546A; accepted by A&A

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., Comment: Accepted for publication in Astronomy and Astrophysics (acceptance date: 06/03/18),15 pages, 22 figures

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., 12 pages, 11 figures, accepted for publication in ApJ

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., Comment: Accepted for publication in A&A. 15 pages (2 page appendix). 16 figures

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��< |b| < 70��$ 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., 21 pages, 17 figures, accepted for publication in ApJ

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., 20 pages, 14 figures, 4 tables

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), Comment: Accepted on A&A. A high-resolution version of the paper is available at http://aramis.obspm.fr/~paola/ELZ/Elz.pdf

One year before the first release of the first data from Gaia, how robust are our views of the Milky Way stellar populations? Recent results have shown that limits, differences and/or continuities between populations are not where we thought they were just a few years ago. The outer disk (> 10kpc) has properties essentially different from the inner (thin+thick) disk, while the bulge is best explained in terms of disk populations, with a negligible or inexistent classical bulge, suggesting that the Milky Way is a pure disk galaxy. Much less contingent than previously envisaged, the thick disk is probably the main phase of stellar mass creation in the MW, and the parent population of the thin disk. These results lead to fundamental changes in our views on the stellar mass growth of the Galaxy, secular mass redistribution in the disk, and imply a change of paradigm of the chemical evolution. I review these different advances, and discuss some of the key questions.

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., 11 pages, 11 figures. Accepted for publication in A&A

To explore the relation between bar formation and star formation in Milky Way-type galaxies quantitatively, we simulated gas-rich disk isolated galaxies. We find that the action of the 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 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.

We explore morphological, kinematic and chemical trends of boxy/peanut (b/p) bulges of Milky Way (MW)-type galaxies, to better understand the formation history of the MW’s bulge. We show, using N-body simulations with both a kinematically cold and a kinematically hot disc, that colder populations develop a more prominent bar and X-shaped peanut as compared to their hotter counterpart. Colder discs also exhibit lower line-of-sight velocities (when viewed edge-on) at the edges of the b/p compared to hot discs, in agreement with what is seen for the MW bulge. Furthermore, we explore an N-body model which has three co-spatial discs with metallicities which correspond to the stellar populations of the inner Milky Way, where the α-enhanced thick disc populations are massive and centrally concentrated. The metallicity trends seen in observations of the Bulge can be reproduced in the model without the need of adding any additional components, which hints to the disc origin of the MW’s bulge.

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., The Astrophysical Journal, in press

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., 5 pages, 6 figures, Accepted for publication in A&A Letters

From the Gaia mission, kinematic and star count data, together with the physical parameters of the stars - ages and metallicities-, will allow to characterise our Galaxy populations and, from that, the overall Galactic gravitational potential. One of the promising procedures to reach such goal will be to optimise the present Population Synthesis models by fitting, through robust statistical techniques, the large and small scale structure and kinematics parameters that best will reproduce Gaia data. As a first step we present here the comparison between the data from the Tycho-2 catalogue and the Besancon Galaxy Model simulations, as well as our present and future work on the model optimization.

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., minor revision

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 (R9-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., 30 pages, 18 figures, 3 tables, accepted by A&A - minor stylistic changes

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., Comment: Accepted for publication in Astronomy and Astrophysics (acceptance date: 27/04/15), 24 pages, 24 figures

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)., A&A accepted, 13 pages

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., 9 pages, 4 figures, accepted for publication on A&A

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., 5 pages, 4 figures; accepted to ApJ Letters

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)., 21 pages, 14 figures, A&A accepted

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, accepted in Astronomy and Astrophysics, 20 pages, 16 figures

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., Accepted for publication in ApJL

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, 18 pages, 17 figures, accepted for publication in A&A

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]$, Comment: 21 pages, 15 figures, accepted for publication in MNRAS

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, Accepted in A&A, Revised version with new figures and extended discussion

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., 8 pages, 10 figures, accepted for publication on Astronomy and Astrophysics