No Arabic abstract
We analyse from an observational perspective the formation history and kinematics of a Milky Way-like galaxy from a high-resolution zoom-in cosmological simulation that we compare to those of our Galaxy as seen by Gaia DR2 to better understand the origin and evolution of the Galactic thin and thick discs. The cosmological simulation was carried out with the GADGET-3 TreePM+SPH code using the MUlti Phase Particle Integrator (MUPPI) model. We disentangle the complex overlapping of stellar generations that rises from the top-down and inside-out formation of the galactic disc. We investigate cosmological signatures in the phase-space of mono-age populations and highlight features stemming from past and recent dynamical perturbations. In the simulation, we identify a satellite with a stellar mass of $1.2 times 10^9$ M$_odot$, i.e. stellar mass ratio $Delta sim 5.5$ per cent at the time, accreted at $z sim 1.6$, which resembles the major merger Gaia-Sausage-Enceladus that produced the Galactic thick disc, i.e. $Delta sim 6$ per cent. We found at $z sim 0.5-0.4$ two merging satellites with a stellar mass of $8.8 times 10^8$ M$_odot$ and $5.1 times 10^8$ M$_odot$ that are associated to a strong starburst in the Star Formation History, which appears fairly similar to that recently found in the Solar Neighbourhood. Our findings highlight that detailed studies of coeval stellar populations kinematics, which are made available by current and future Gaia data releases and in synergy with simulations, are fundamental to unravel the formation and evolution of the Milky Way discs.
We present the analysis of a suite of simulations run with different particle-and grid-based cosmological hydrodynamical codes and compare them with observational data of the Milky Way. This is the first study to make comparisons of properties of galaxies simulated with particle and grid-based codes. Our analysis indicates that there is broad agreement between these different modelling techniques. We study the velocity dispersion - age relation for disc stars at z=0 and find that four of the simulations are more consistent with observations by Holmberg et al. (2008) in which the stellar disc appears to undergo continual/secular heating. Two other simulations are in better agreement with the Quillen & Garnett (2001) observations that suggest a saturation in the heating profile for young stars in the disc. None of the simulations have thin discs as old as that of the Milky Way. We also analyse the kinematics of disc stars at the time of their birth for different epochs in the galaxies evolution and find that in some simulations old stars are born cold within the disc and are subsequently heated, while other simulations possess old stellar populations which are born relatively hot. The models which are in better agreement with observations of the Milky Ways stellar disc undergo significantly lower minor-merger/assembly activity after the last major merger - i.e. once the disc has formed. All of the simulations are significantly hotter than the Milky Way disc; on top of the effects of mergers, we find a floor in the dispersion that is related to the underlying treatment of the heating and cooling of the interstellar medium, and the low density threshold which such codes use for star formation. This finding has important implications for all studies of disc heating that use hydrodynamical codes.
We develop a chemical evolution model in order to study the star formation history of the Milky Way. Our model assumes that the Milky Way is formed from a closed box-like system in the inner regions, while the outer parts of the disc experience some accretion. Unlike the usual procedure, we do not fix the star formation prescription (e.g. Kennicutt law) in order to reproduce the chemical abundance trends. Instead, we fit the abundance trends with age in order to recover the star formation history of the Galaxy. Our method enables one to recover with unprecedented accuracy the star formation history of the Milky Way in the first Gyrs, in both the inner (R<7-8kpc) and outer (R>9-10kpc) discs as sampled in the solar vicinity. We show that, in the inner disc, half of the stellar mass formed during the thick disc phase, in the first 4-5 Gyr. This phase was followed by a significant dip in the star formation activity (at 8-9 Gyr) and a period of roughly constant lower level star formation for the remaining 8 Gyr. The thick disc phase has produced as many metals in 4 Gyr as the thin disc in the remaining 8 Gyr. Our results suggest that a closed box model is able to fit all the available constraints in the inner disc. A closed box system is qualitatively equivalent to a regime where the accretion rate, at high redshift, maintains a high gas fraction in the inner disc. In such conditions, the SFR is mainly governed by the high turbulence of the ISM. By z~1 it is possible that most of the accretion takes place in the outer disc, while the star formation activity in the inner disc is mostly sustained by the gas not consumed during the thick disc phase, and the continuous ejecta from earlier generations of stars. The outer disc follows a star formation history very similar to that of the inner disc, although initiated at z~2, about 2 Gyr before the onset of the thin disc formation in the inner disc.
The formation of the Galactic disc is an enthusiastically debated issue. Numerous studies and models seek to identify the dominant physical process(es) that shaped its observed properties. Taking advantage of the improved coverage of the inner Milky Way provided by the SDSS DR16 APOGEE catalogue and of the ages published in the APOGEE-AstroNN Value Added Catalogue (VAC), we examine the radial evolution of the chemical and age properties of the Galactic stellar disc, with the aim to better constrain its formation. Using a sample of 199,307 giant stars with precise APOGEE abundances and APOGEE-astroNN ages, selected in a +/-2 kpc layer around the galactic plane, we assess the dependency with guiding radius of: (i) the median metallicity, (ii) the ridge lines of the [Fe/H]-[Mg/Fe] and age-[Mg/Fe] distributions and (iii) the Age Distribution Function (ADF). The giant star sample allows us to probe the radial behaviour of the Galactic disc from Rg = 0 to 14-16 kpc. The thick disc [Fe/H]-[Mg/Fe] ridge lines follow closely grouped parallel paths, supporting the idea that the thick disc did form from a well-mixed medium. However, the ridge lines present a small drift in [Mg/Fe], which decreases with increasing guiding radius. At sub-solar metallicity, the intermediate and outer thin disc [Fe/H]-[Mg/Fe] ridge lines follow parallel sequences shifted to lower metallicity as the guiding radius increases. We interpret this pattern, as the signature of a dilution of the inter-stellar medium from Rg~6 kpc to the outskirt of the disc, which occured before the onset of the thin disc formation. The APOGEE-AstroNN VAC provides stellar ages for statistically significant samples of thin disc stars from the Galactic centre up to Rg~14 kpc. An important result provided by this dataset, is that the thin disc presents evidence of an inside-out formation up to R_g~10-12 kpc.(Abridged)
We apply a semi-analytic galaxy formation model to two high resolution cosmological N-body simulations to investigate analogues of the Milky Way system. We select these according to observed properties of the Milky Way rather than by halo mass as in most previous work. For disk-dominated central galaxies with stellar mass (5--7) x 10d10Msun, the median host halo mass is 1.4 x 10d12Msun, with 1 sigma dispersion in the range [0.86, 3.1] x 10d12Msun, consistent with dynamical measurements of the Milky Way halo mass. For any given halo mass, the probability of hosting a Milky Way system is low, with a maximum of ~20% in haloes of mass ~10d12Msun. The model reproduces the V-band luminosity function and radial profile of the bright (MV < -9) Milky Way satellites. Galaxy formation in low mass haloes is found to be highly stochastic, resulting in an extremely large scatter in the relation between MV (or stellar mass) for satellites and the depth of the subhalo potential well in which they live, as measured by the maximum of the rotation curve, Vmax. We conclude that the too big to fail problem is an artifact of selecting satellites in N-body simulations according to subhalo properties: in 10% of cases we find that three or fewer of the brightest (or most massive) satellites have Vmax > 30 km/s. Our model predicts that around half of the dark matter subhaloes with Vmax > 20 km/s host satellites fainter than MV = -9 and so may be missing from existing surveys.
Binary black holes are the primary endpoint of massive stellar evolution. Their properties provide a unique opportunity to constrain binary evolution, which is still poorly understood. In this paper, we predict the inventory of binary black holes and their merger products in/around the Milky Way, and detail their main properties. We present the first combination of a high-resolution cosmological simulation of a Milky Way-mass galaxy with a binary population synthesis model. The hydrodynamic simulation, taken from the FIRE project, provides a cosmologically realistic star formation history for the galaxy and its stellar halo and satellites. We apply a metallicity-dependent evolutionary model to the star particles to produce individual binary black holes. We find that a million binary black holes have merged in the model Milky Way, and 3 million binaries are still present, with an average mass of 28 Msun per binary. Because the black hole progenitors are biased towards low metallicity stars, half reside in the stellar halo and satellites and 40 per cent of the binaries were formed outside the main galaxy. This trend increases with the masses of the black holes. The numbers and mass distribution of the merged systems is compatible with the LIGO/Virgo detections. Observations of these black holes will be challenging, both with electromagnetic methods and LISA. We find that a cosmologically realistic star formation history, with self-consistent metal enrichment and Galactic accretion history, are key ingredients for determining binary black hole rates that can be compared with observations to constrain massive binary evolution.