Search

The presence of dark gaps, a preferential light deficit along the bar minor axis, is observationally well known. The properties of dark gaps are thought to be associated with the properties of bars, and their spatial locations are often associated with bar resonances. However, a systematic study, testing the robustness and universality of these assumptions, is still largely missing. Here, we investigate the formation and evolution of bar-induced dark gaps using a suite of N-body models of (kinematically cold) thin and (kinematically hot) thick discs with varying thick disc mass fraction, and different thin-to-thick disc geometry. We find that dark gaps are a natural consequence of the trapping of disc stars by the bar. The properties of dark gaps (such as strength and extent) are well correlated with the properties of bars. For stronger dark gaps, the fractional mass loss along the bar minor axis can reach up to ~60-80 percent of the initial mass contained, which is redistributed within the bar. These trends hold true irrespective of the mass fraction in the thick disc and the assumed disc geometry. In all our models harbouring slow bars, none of the resonances (corotation, Inner Lindblad resonance, and 4:1 ultra-harmonic resonance) associated with the bar correspond to the location of dark gaps, thereby suggesting that the location of dark gaps is not a universal proxy for these bar resonances, in contrast with earlier studies., Comment: 13 pages, 12 figures, 1 table (including appendix), accepted for publication in MNRAS

Previous studies of the chemo-kinematic properties of stars in the Galactic bulge have revealed a puzzling trend. Along the bulge minor axis, and close to the Galactic plane, metal-rich stars display a higher line-of-sight velocity dispersion compared to metal-poor stars, while at higher latitudes metal-rich stars have lower velocity dispersions than metal-poor stars, similar to what is found in the Galactic disc. In this work, we re-examine this issue, by studying the dependence of line-of-sight velocity dispersions on metallicity and latitude in the latest APOGEE Data Release 17, confirming the results of previous works. We then analyse an N-body simulation of a Milky Way-like galaxy, also taking into account observational biases introduced by the APOGEE selection function. We show that the inversion in the line-of-sight velocity dispersion-latitude relation observed in the Galactic bulge can be reproduced by our model. We show that this inversion is a natural consequence of a scenario in which the bulge is a boxy/peanut-shaped structure, whose metal-rich and metal-poor stars mainly originate from the thin and thick disc of the Milky Way, respectively. Due to their cold kinematics, metal-rich, thin disc stars, are efficiently trapped in the boxy/peanut bulge, and, at low latitudes, show a strong barred morphology, which results in high velocity dispersions which are larger than those attained by the metal-poor populations. Extremely metal-rich stars in the Galactic bulge, which have received renewed attention in the literature, do follow the same trends as those of the metal-rich populations. The line-of-sight velocity-latitude relation observed in the Galactic bulge for metal-poor and metal-rich stars are thus both an effect of the intrinsic nature of the Galactic bulge and of the angle at which we observe it from the Sun., Comment: 12 pages, 9 figures

We derive the metallicity profile of the Milky Way low-$\alpha$ disc population from 2 to 20 kpc from the Galactic centre in 1 Gyr age bins using the astroNN catalogue, and show that it is highly structured, with a plateau between 4 and 7 kpc and a break at 10-12 kpc. We argue that these features result from the two main bar resonances, the corotation and the Outer Lindblad Resonance (OLR), respectively. We show that the break in the metallicity profile is most visible in stars having 7-8 Gyr, reaching an amplitude of about 0.4 dex, and is the signpost of the position of the bar OLR. The bar formation was accompanied by an episode of radial migration triggered by its slowing down and is responsible for spreading old metal-rich stars up to the OLR. The data show that the slowdown of the bar ended 6-7 Gyr ago. Based on numerical simulations that reproduce well the break observed in the metallicity profile, we argue that this implies that the bar formed in our Galaxy 8-10 Gyr ago. Analysis of the metallicity distribution as a function of radius shows no evidence of significant systematic outward radial migration after this first episode. We argue that the variation of the metallicity dispersion as a function of the guiding radius is dominated by the migration triggered by the bar, but also that the libration of orbits around the bar resonances induces a mixing that may have a significant impact on the observed metallicity dispersion. In contrast, the absence of a break in the metallicity profile of populations younger than about $\sim$6 Gyr and the flattening of the gradient at younger ages is interpreted as evidence that the strength of the bar has decreased, loosening its barrier effect and allowing the gas and metals on both sides of the OLR to mix, erasing the break. Beyond the OLR, stars younger than 7 Gyr show very small metallicity dispersion, suggesting no or limited migration., Comment: 20 pages, 23 figures. Accepted in A&A

Boxy/peanut (b/p) bulges, the vertically extended inner parts of bars, are ubiquitous in barred galaxies in the local Universe, including our own Milky Way. At the same time, a majority of external galaxies and the Milky Way also possess a thick-disc. However, the dynamical effect of thick-discs in the b/p formation and evolution is not fully understood. Here, we investigate the effect of thick-discs in the formation and evolution of b/ps by using a suite of N-body models of (kinematically cold) thin and (kinematically hot) thick discs. Within the suite of models, we systematically vary the mass fraction of the thick disc, and the thin-to-thick disc scale length ratio. The b/ps form in almost all our models via a vertical buckling instability, even in the presence of a massive thick disc. The thin disc b/p is much stronger than the thick disc b/p. With increasing thick disc mass fraction, the final b/p structure gets progressively weaker in strength and larger in extent. Furthermore, the time-interval between the bar formation and the onset of buckling instability gets progressively shorter with increasing thick-disc mass fraction. The breaking and restoration of the vertical symmetry (during and after the b/p formation) show a spatial variation -- the inner bar region restores vertical symmetry rather quickly (after the buckling) while in the outer bar region, the vertical asymmetry persists long after the buckling happens. Our findings also predict that at higher redshifts, when discs are thought to be thicker, b/ps would have more 'boxy-shaped' appearance than more 'X-shaped' appearance. This remains to be tested from future observations at higher redshifts., Comment: 19 pages, 1 table, 20 figures (including appendices), accepted for publication in A&A

One of the challenges related to stellar bars is to accurately determine the length of the bar in a disc galaxy. In the literature, a wide variety of methods have been employed to measure the extent of a bar. However, a systematic study on determining the robustness and accuracy of different bar length estimators is still beyond our grasp. Here, we investigate the accuracy and the correlation (if any) between different bar length measurement methods while using an N-body model of a barred galaxy, where the bar evolves self-consistently in the presence of a live dark matter halo. We investigate the temporal evolution of the bar length, using different estimators (involving isophotal analysis of de-projected surface brightness distribution and Fourier decomposition of surface density), and we study their robustness and accuracy. We made further attempts to determine correlations among any two of these bar length estimators used here. In the presence of spirals, the bar length estimators that only consider the amplitudes of different Fourier moments (and do not take into account the phase-angle of m=2 Fourier moment) systematically overestimate the length of the bar. The strength of dark-gaps (produced by bars) is strongly correlated with the bar length in early rapid growth phase and is only weakly anti-correlated during subsequent quiescent phase of bar evolution. However, the location of dark-gaps is only weakly correlated to the bar length, hence, this information cannot be used as a robust proxy for determining the bar length. In addition, the bar length estimators, obtained using isophotal analysis of de-projected surface brightness distribution, systematically overestimate the bar length. The implications of bar length over(under)estimation in the context of determining fast and slow bars are further discussed in this work., Comment: 14 pages, 11 figures (including appendices), accepted for publication in A&A

Young alpha-rich (YAR) stars have been detected in the past as outliers to the local age $\rm-$ [$\alpha$/Fe] relation. These objects are enhanced in $\alpha$-elements but apparently younger than typical thick disc stars. We study the global kinematics and chemical properties of YAR giant stars in APOGEE DR17 survey and show that they have properties similar to those of the standard thick disc stellar population. This leads us to conclude that YAR are rejuvenated thick disc objects, most probably evolved blue stragglers. This is confirmed by their position in the Hertzsprung-Russel diagram (HRD). Extending our selection to dwarfs allows us to obtain the first general straggler distribution in an HRD of field stars. We also compare the elemental abundances of our sample with those of standard thick disc stars, and find that our YAR stars are shifted in oxygen, magnesium, sodium, and the slow neutron-capture element cerium. Although we detect no sign of binarity for most objects, the enhancement in cerium may be the signature of a mass transfer from an asymptotic giant branch companion. The most massive YAR stars suggest that mass transfer from an evolved star may not be the only formation pathway, and that other scenarios, such as collision or coalescence should be considered., Comment: 18 Pages, 20 Figures, 1 Table; accepted for publication in Astronomy & Astrophysics

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

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 $z<0.5$ cluster galaxies; (vi) survey stellar populations and kinematics in $\sim25\,000$ field galaxies at $0.3\lesssim z \lesssim 0.7$; (vii) study the cosmic evolution of accretion and star formation using $>1$ 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., Comment: 41 pages, 27 figures, accepted for publication by MNRAS; updated version including information on individual grants in a revised Acknowledgements section, corrections to the affiliation list, and an updated references list

The Milky Way as well as a majority of external galaxies possess a thick disc. However, the dynamical role of the (geometrically) thick disc on the bar formation and evolution is not fully understood. Here, we investigate the effect of thick discs in bar formation and evolution by means of a suite of N-body models of (kinematically cold) thin-(kinematically hot) thick discs. We systematically vary the mass fraction of the thick disc, the thin-to-thick disc scale length ratio as well as thick disc's scale height to examine the bar formation under diverse dynamical scenarios. Bars form almost always in our models, even in presence of a massive thick disc. The part of the bar constituted by the thick disc closely follows the overall growth and temporal evolution of the part of the bar constituted by the thin disc, only the part of the bar in the thick disc is weaker than the part of the bar in the thin disc. The formation of stronger bars is associated with a simultaneous larger loss of angular momentum and a larger radial heating. In addition, we demonstrate a preferential loss of angular momentum and a preferential radial heating of disc stars, along the azimuthal direction within the extent of the bar, in both thin and thick disc stars. For purely thick disc models (without any thin disc), the bar formation critically depends on the disc scale length and scale height. A larger scale length and/or a larger vertical scale height delays the bar formation time and/or suppresses the bar formation almost completely in thick-disc-only models. We find that the Ostriker-Peeble criterion predicts the bar instability scenarios in our models better than the Efstathiou-Lake-Negroponte criterion., Comment: 20 pages, 14 figures, 1 table (including appendix), accepted for publication in A&A

Reconstructing how all the stellar components of the Galaxy formed and assembled over time, by studying the properties of the stars which make it, is the aim of Galactic archeology. In these last years, thanks to the launch of the ESA Gaia astrometric mission, and the development of many spectroscopic surveys, we are for the first time in the position to delve into the layers of the past of our galaxy. Globular clusters (GCs) play a fundamental role in this research field since they are among the oldest stellar systems in the Milky Way (MW) and so bear witness of its entire past. In the recent years, there have been several attempts to constrain the nature of clusters (accreted or formed in the MW itself) through the analysis of kinematic spaces and to reconstruct from this the properties of the accretions events experienced by the MW through time. This work aims to test a widely-used assumption about the clustering of the accreted populations of GCs in the integrals of motions space. We analyze a set of dissipation-less N-body simulations that reproduce the accretion of one or two satellites with their GC population on a MW-type galaxy. Our results demonstrate that a significant overlap between accreted and "kinematically-heated" in-situ GCs is expected in kinematic spaces, for mergers with mass ratios of 1:10. In contrast with standard assumptions made in the literature so far, we find that accreted GCs do not show dynamical coherence, that is they do not cluster in kinematic spaces. In addition, GCs can also be found in regions dominated by stars which have a different origin (i.e. different progenitor). This casts doubt on the association between GCs and field stars that is generally made in the literature to assign them to a common origin. Our findings severely question the recovered accretion history of the MW based on the phase-space clustering of the GC population., Comment: Accepted for publication on A&A. Major additions implemented to the first version

In the Milky Way, recent progress in the exploration of its assembly history is driven by the tremendous amount of high-quality data delivered by Gaia, which has revealed a number of substructures potentially linked to several ancient accretion events. In this work, aiming to explore the phase-space structure of accreted stars, we analyze six M31/MW analogues from the HESTIA suite of cosmological hydrodynamics zoom-in simulations of the Local Group. We found that all the HESTIA galaxies experience a few dozen mergers but only 1-4 mergers have the stellar mass ratio >0.2 where, depending on the halo definition, the most massive merger contributes from 20% to 70% of the total stellar halo. Individual merger remnants show diverse density distributions at z=0, significantly overlapping with each other and with the in-situ stars in the ELz, UV and RVphi coordinates. The mergers debris often change their position in the ELz with time due to the galactic mass growth and the non-axisymmetry of the potential. In agreement with previous works, we show that even individual merger debris exhibit a number of distinct ELz features. In the UV plane, all HESTIA galaxies reveal radially hot, non-rotating or weakly counter-rotating, Gaia-Sausage-like features. We found an age gradient in Elz space for the individual debris, where the youngest stars, formed in the inner regions of accreting systems, deposit to the innermost regions of the host. The bulk of these stars is being formed during the last stages of accretion, making it possible to date the merger. In actions space (Jr, Jz, J\phi), the mergers debris do not appear as isolated substructures but are instead scattered over a large parameters area and overlapping with the in-situ stars. We also introduce a purely kinematic space (Jz/Jr-eccentricity), where different merger debris can be disentangled better from each other and from the in-situ stars., Comment: Paper II in a series of III, 22 pages, 16 figures, A&A in press

Theory suggests that mergers play an important role in shaping galactic discs and stellar haloes, which was observationally confirmed in the MW thanks to Gaia data. In this work, aiming to probe the contribution of mergers to the in situ stellar halo formation, we analyse six M31/MW analogues from the HESTIA suite of cosmological hydrodynamical zoom-in simulations of the LG. We found that all the HESTIA galaxies experience between one to four mergers with stellar mass ratios between 0.2 and 1 relative to the host at the time of the merger. These significant mergers, with a single exception, happened 7-11Gyr ago. The overall impact of the most massive mergers in HESTIA is clearly seen as a sharp increase in the orbital eccentricity (and a corresponding decrease in the rotational velocity Vphi of pre-existing disc stars of the main progenitor, thus nicely reproducing the Splash-, Plume-like feature that was discovered in the MW. We do find a correlation between mergers and close pericentric passages of massive satellites and bursts of star formation in the in situ component. Massive mergers sharply increase the disc velocity dispersion of the in situ stars; however, the latest significant merger often heats up the disc up to the numbers when the contribution of the previous ones is less prominent in the age-velocity dispersion relation. In HESTIA galaxies, the in situ halo is an important component of the inner stellar halo where its fraction is about 30-40%, while in the outer parts it typically does not exceed ~5% beyond 15 kpc. The simulations suggest that this component of the stellar haloes continues to grow well after mergers conclude; however, the most significant contribution comes from stars that formed recently before the merger. The orbital analysis of the HESTIA galaxies suggests that wedges in Rmax-Zmax space are mainly populated by the stars born between significant mergers., Comment: Paper I in a series of III, 17 pages, 15 figures, A&A in press

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 time scale, this has not been directly constrained using observations. Using the unprecedented amount and quality of data of the inner regions of the Milky Way that has recently been produced by APOGEE and Gaia, we aim to derive strong constraints on the infall history of the inner (less than 6 kpc) Galaxy (with a focus on stars between 4-6 kpc, which we show is an appropriate proxy for the entire inner disc). We have 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 between 4-6 kpc from APOGEE. Our results point to a fast early gas accretion, with an upper limit of 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., Comment: 23 pages, 17 figures, 4 tables, submitted to A&A

An $m=1$ lopsided asymmetry is common in disc galaxies. Here, we investigate the excitation of an $m=1$ lopsidedness in host galaxies during minor mergers while choosing a set of 1:10 merger models (with varying orbital configurations, morphology of the host galaxy) from the GalMer galaxy merger library. We show that a minor merger triggers a prominent $m=1$ lopsidedness in stars of the host galaxy. The strength of the $m=1$ lopsidedness undergoes a transient amplification phase after each pericenter passage of the satellite, in concordance with past findings of exciting an $m=1$ lopsidedness by tidal encounters. However, once the merger happens, and the post-merger remnant readjusts itself, the lopsidedness disappears in short time-scale ($\sim$ 500-850 Myr). Furthermore, a delayed merger can drive a prolonged ($\sim$2 Gyr) lopsidedness in the host galaxy. We demonstrate that the $m=1$ lopsidedness rotates with a well-defined pattern speed which is much slower than the $m=2$ bar pattern speed, and is retrograde with respect to the bar. This gives rise to a dynamical scenario where the Inner Lindblad resonance of the $m=1$ lopsidedness falls in between the corotation and the Outer Lindblad resonance of the $m=2$ bar mode. A kinematic lopsidedness also arises in the host galaxy; the resulting temporal variation closely follows that of the density lopsidedness. The minor merger also triggers a transient off-centred stellar disc-dark matter halo configuration due to the tidal encounter with the satellite., Comment: 19 pages, 1 table, 16 figures (including appendices), accepted for publication in MNRAS

Minor merger of galaxies are common during the evolutionary phase of galaxies. Here, we investigate the dynamical impact of a minor merger (mass ratio 1:10) event on the final fate of a stellar bar in the merger remnant. To achieve that, we choose a set of minor merger models from the publicly available GalMer library of galaxy merger simulations. The models differ in terms of their orbital energy, orientation of the orbital spin vector, and morphology of the satellite galaxy (discy/spheroidal). We demonstrate that the central stellar bar, initially present in the host galaxy, undergoes a transient bar amplification phase after each pericentre passage of the satellite; in concordance with past studies of bar excitation due to tidal encounter. However, once the merger happens, the central stellar bar weakens substantially in the post-merger remnants. The accumulation of satellite's stars in the central region of merger remnant plays a key role in the bar weakening process; causing a net increase in the central mass concentration as well as in the specific angular momentum content. We find that the efficiency of mass accumulation from the satellite in the central parts of merger remnants depends on the orbital parameters as well as on the satellite's morphology. Consequently, different minor merger models display different degrees of bar weakening event. This demonstrates that minor merger of galaxies is a plausible avenue for bar weakening in disc galaxies., Comment: 16 pages, 16 figures, 1 table, accepted for publication in MNRAS

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

Boxy, peanut- or X-shaped "bulges" are observed in a large fraction of barred galaxies viewed in, or close to, edge-on projection, as well as in the Milky Way. They are the product of dynamical instabilities occurring in stellar bars, which cause the latter to buckle and thicken vertically. Recent studies have found nearby galaxies that harbour two such features arising at different radial scales, in a nested configuration. In this paper we explore the formation of such double peanuts, using a collisionless N-body simulation of a pure disc evolving in isolation within a live dark matter halo, which we analyse in a completely analogous way to observations of real galaxies. In the simulation we find a stable double configuration consisting of two X/peanut structures associated to the same galactic bar - rotating with the same pattern speed - but with different morphology, formation time, and evolution. The inner, conventional peanut-shaped structure forms early via the buckling of the bar, and experiences little evolution once it stabilises. This feature is consistent in terms of size, strength and morphology, with peanut structures observed in nearby galaxies. The outer structure, however, displays a strong X, or "bow-tie", morphology. It forms just after the inner peanut, and gradually extends in time (within 1 to 1.5 Gyr) to almost the end of the bar, a radial scale where ansae occur. We conclude that, although both structures form, and are dynamically coupled to, the same bar, they are supported by inherently different mechanisms., Comment: 12 pages, 11 figures

(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., Comment: 9 chapters, 301 pages, 100 figures. This version of the DSC is a comprehensive update of the original version, released in 2016, which can be downloaded at arXiv:1606.00043. A detailed summary of the design of MSE is available in the MSE Book 2018, available at arXiv:1810.08695

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., Comment: 7 pages, 6 figures; accepted for publication by A&A Letters on Dec 14, 2018

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

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 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., Comment: 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] $<$ 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., 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

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

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

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

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

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), Comment: 210 pages, 91 figures. Exposure draft. Appendices to the Detailed Science Case can be found at http://mse.cfht.hawaii.edu/docs/

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

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

We explore a second order Hamiltonian vertical resonance model for X-shaped or peanut-shaped galactic bulges. The X-shape is caused by the 2:1 vertical Lindblad resonance with the bar, with two vertical oscillation periods per orbital period in the bar frame. We examine N-body simulations and find that due to the bar slowing down and disk thickening during bar buckling, the resonance and associated peanut-shape moves outward. The peanut-shape is consistent with the location of the vertical resonance, independent of whether the bar buckled or not. We estimate the resonance width from the potential m=4 Fourier component and find that the resonance is narrow, affecting orbits over a narrow range in the angular momentum distribution, dL/L ~ 0.05. As the resonance moves outward, stars originally in the mid plane are forced out of the mid plane into orbits just within the resonance separatrix. The height of the separatrix orbits, estimated from the Hamiltonian model, is approximately consistent with the peanut-shape height. The X-shape is comprised of stars in the vicinity of the resonance separatrix. The velocity distributions from the simulations illustrate that low inclination orbits are depleted within resonance. Within resonance, the vertical velocity distribution is broad, consistent with resonant heating caused by the passage of the resonance through the disk. In the Milky Way bulge we relate the azimuthally averaged mid-plane mass density near the vertical resonance to the rotation curve and bar pattern speed. At an estimated vertical resonance galactocentric radius of ~1.3 kpc, we confirm a mid-plane density of ~5x10^8 Msol/kpc^3, consistent with recently estimated mass distributions. We find that the rotation curve, bar pattern speed, 2:1 vertical resonance location, X-shape tips, and mid-plane mass density, are all self-consistent in the Milky Way galaxy bulge., Comment: accepted for publication in MNRAS

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., Comment: 4 pages, 3 figures. To appear in the proceedings of the conference 'Reading the book of globular clusters with the lens of stellar evolution', Rome, Nov. 26-28, 2012, to be published in Memorie della Societa Astronomica Italiana, Eds. Paolo Ventura and Corinne Charbonnel

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., Comment: 6 pages, 7 figures. A&A Letters, accepted. Top panel of Fig. A1 replaced, minor typos corrected. High resolution version available at http://mygepi.obspm.fr/~paola/Pal5/