The biggest splash Belokurov, Vasily; Sanders, Jason L; Fattahi, Azadeh ...
Monthly notices of the Royal Astronomical Society,
05/2020, Letnik:
494, Številka:
3
Journal Article
Recenzirano
Odprti dostop
ABSTRACT
Using a large sample of bright nearby stars with accurate Gaia Data Release 2 astrometry and auxiliary spectroscopy we map out the properties of the principle Galactic components such as the ...‘thin’ and ‘thick’ discs and the halo. We confirm previous claims that in the Solar neighbourhood, there exists a large population of metal-rich (Fe/H > −0.7) stars on highly eccentric orbits. By studying the evolution of elemental abundances, kinematics, and stellar ages in the plane of azimuthal velocity vϕ and metallicity Fe/H, we demonstrate that this metal-rich halo-like component, which we dub the Splash, is linked to the α-rich (or ‘thick’) disc. Splash stars have little to no angular momentum and many are on retrograde orbits. They are predominantly old, but not as old as the stars deposited into the Milky Way (MW) in the last major merger. We argue, in agreement with several recent studies, that the Splash stars may have been born in the MW’s protodisc prior to the massive ancient accretion event which drastically altered their orbits. We cannot, however, rule out other (alternative) formation channels. Taking advantage of the causal connection between the merger and the Splash, we put constraints of the epoch of the last massive accretion event to have finished 9.5 Gyr ago. The link between the local metal-rich and metal-poor retrograde stars is confirmed using a large suite of cutting-edge numerical simulations of the MW’s formation.
Abstract
We use 12 cosmological N-body simulations of Local Group systems (the apostle models) to inspect the relation between the virial mass of the main haloes (Mvir,1 and Mvir,2), the mass derived ...from the relative motion of the halo pair (Mtim), and that inferred from the local Hubble flow (Mlhf). We show that within the spherical collapse model (SCM), the correspondence between the three mass estimates is exact, i.e. Mlhf = Mtim = Mvir,1 + Mvir,2. However, comparison with apostle simulations reveals that, contrary to what the SCM states, a relatively large fraction of the mass that perturbs the local Hubble flow and drives the relative trajectory of the main galaxies is not contained within Rvir, and that the amount of ‘extravirial’ mass tends to increase in galaxies with a slow accretion rate. In contrast, modelling the peculiar velocities around the Local Group returns an unbiased constraint on the virial mass ratio of the main galaxy pair. Adopting the outer halo profile found in N-body simulations, which scales as ρ ∼ R−4 at R ≳ Rvir, indicates that the galaxy masses perturbing the local Hubble flow roughly correspond to the asymptotically convergent (total) masses of the individual haloes. We show that estimates of Mvir based on the dynamics of tracers at R ≫ Rvir require a priori information on the internal matter distribution and the growth rate of the main galaxies, both of which are typically difficult to quantify.
ABSTRACT
We determine the Milky Way (MW) mass profile inferred from fitting physically motivated models to the Gaia DR2 Galactic rotation curve and other data. Using various hydrodynamical ...simulations of MW-mass haloes, we show that the presence of baryons induces a contraction of the dark matter (DM) distribution in the inner regions, r ≲ 20 kpc. We provide an analytic expression that relates the baryonic distribution to the change in the DM halo profile. For our galaxy, the contraction increases the enclosed DM halo mass by factors of roughly 1.3, 2, and 4 at radial distances of 20, 8, and 1 kpc, respectively compared to an uncontracted halo. Ignoring this contraction results in systematic biases in the inferred halo mass and concentration. We provide a best-fitting contracted NFW halo model to the MW rotation curve that matches the data very well.1 The best-fit has a DM halo mass, $M_{200}^{\rm DM}=0.97_{-0.19}^{+0.24}\times 10^{12}\,\mathrm{M}_\odot$, and concentration before baryon contraction of $9.4_{-2.6}^{+1.9}$, which lie close to the median halo mass–concentration relation predicted in ΛCDM. The inferred total mass, $M_{200}^{\rm total}=1.08_{-0.14}^{+0.20} \times 10^{12}\,\mathrm{M}_\odot$, is in good agreement with recent measurements. The model gives an MW stellar mass of $5.04_{-0.52}^{+0.43}\times 10^{10}\,\mathrm{M}_\odot$ and infers that the DM density at the Solar position is $\rho _{\odot }^{\rm DM}=8.8_{-0.5}^{+0.5}\times 10^{-3}\,\mathrm{M}_\odot \,\mathrm{pc}^{-3}\equiv 0.33_{-0.02}^{+0.02}\,\rm {GeV}\,\rm {cm}^{-3}$. The rotation curve data can also be fitted with an uncontracted NFW halo model, but with very different DM and stellar parameters. The observations prefer the physically motivated contracted NFW halo, but the measurement uncertainties are too large to rule out the uncontracted NFW halo.
Abstract
We model the fastest moving ($v_{\rm tot} \gt 300 \, {\rm km}\, \, {\rm s}^{-1}$) local (D ≲ 3 kpc) halo stars using cosmological simulations and six-dimensional Gaia data. Our approach is ...to use our knowledge of the assembly history and phase-space distribution of halo stars to constrain the form of the high-velocity tail of the stellar halo. Using simple analytical models and cosmological simulations, we find that the shape of the high-velocity tail is strongly dependent on the velocity anisotropy and number density profile of the halo stars – highly eccentric orbits and/or shallow density profiles have more extended high-velocity tails. The halo stars in the solar vicinity are known to have a strongly radial velocity anisotropy, and it has recently been shown the origin of these highly eccentric orbits is the early accretion of a massive ($M_{\rm star}\sim 10^9 \, \mathrm{M}_\odot$) dwarf satellite. We use this knowledge to construct a prior on the shape of the high-velocity tail. Moreover, we use the simulations to define an appropriate outer boundary of 2r200, beyond which stars can escape. After applying our methodology to the Gaia data, we find a local (r0 = 8.3 kpc) escape speed of $v_{\rm esc}(r_0) = 528^{+24}_{-25} \, {\rm km}\, \, {\rm s}^{-1}$. We use our measurement of the escape velocity to estimate the total Milky Way mass, and dark halo concentration: $M_{200, \rm tot} = 1.00^{+0.31}_{-0.24} \times 10^{12}\, \mathrm{M}_\odot$, $c_{200}=10.9^{+4.4}_{-3.3}$. Our estimated mass agrees with recent results in the literature that seem to be converging on a Milky Way mass of $M_{200, \rm tot} \sim 10^{12}\, \mathrm{M}_\odot$.
We use a large sample of isolated dark matter halo pairs drawn from cosmological N-body simulations to identify candidate systems whose kinematics match that of the Local Group (LG) of galaxies. We ...find, in agreement with the ‘timing argument’ and earlier work, that the separation and approach velocity of the Milky Way (MW) and Andromeda (M31) galaxies favour a total mass for the pair of ∼5 × 1012 M⊙. A mass this large, however, is difficult to reconcile with the small relative tangential velocity of the pair, as well as with the small deceleration from the Hubble flow observed for the most distant LG members. Halo pairs that match these three criteria have average masses a factor of ∼2 times smaller than suggested by the timing argument, but with large dispersion. Guided by these results, we have selected 12 halo pairs with total mass in the range 1.6–3.6 × 1012 M⊙ for the apostle project (A Project Of Simulating The Local Environment), a suite of hydrodynamical resimulations at various numerical resolution levels (reaching up to ∼104 M⊙ per gas particle) that use the subgrid physics developed for the eagle project. These simulations reproduce, by construction, the main kinematics of the MW–M31 pair, and produce satellite populations whose overall number, luminosities, and kinematics are in good agreement with observations of the MW and M31 companions. The apostle candidate systems thus provide an excellent testbed to confront directly many of the predictions of the Λ cold dark matter cosmology with observations of our local Universe.
ABSTRACT
We use a distribution function analysis to estimate the mass of the Milky Way (MW) out to 100 kpc using a large sample of halo stars. These stars are compiled from the literature, and the ...vast majority (${\sim } 98{{\ \rm per\ cent}}$) have 6D phase-space information. We pay particular attention to systematic effects, such as the dynamical influence of the Large Magellanic Cloud (LMC), and the effect of unrelaxed substructure. The LMC biases the (pre-LMC infall) halo mass estimates towards higher values, while realistic stellar haloes from cosmological simulations tend to underestimate the true halo mass. After applying our method to the MW data, we find a mass within 100 kpc of M (<100 kpc) = 6.07 ± 0.29 (stat.) ± 1.21 (sys.) × 1011 M⊙. For this estimate, we have approximately corrected for the reflex motion induced by the LMC using the Erkal et al. model, which assumes a rigid potential for the LMC and MW. Furthermore, stars that likely belong to the Sagittarius stream are removed, and we include a 5 per cent systematic bias, and a 20 per cent systematic uncertainty based on our tests with cosmological simulations. Assuming the mass–concentration relation for Navarro–Frenk–White haloes, our mass estimate favours a total (pre-LMC infall) MW mass of M200c = 1.01 ± 0.24 × 1012 M⊙, or (post-LMC infall) mass of M200c = 1.16 ± 0.24 × 1012 M⊙ when a 1.5 × 1011 M⊙ mass of a rigid LMC is included.
ABSTRACT
N-body simulations make unambiguous predictions for the abundance of substructures within dark matter haloes. However, the inclusion of baryons in the simulations changes the picture because ...processes associated with the presence of a large galaxy in the halo can destroy subhaloes and substantially alter the mass function and velocity distribution of subhaloes. We compare the effect of galaxy formation on subhalo populations in two state-of-the-art sets of hydrodynamical Λcold dark matter (ΛCDM) simulations of Milky Way mass haloes, Apostle and Auriga. We introduce a new method for tracking the orbits of subhaloes between simulation snapshots that gives accurate results down to a few kiloparsecs from the centre of the halo. Relative to a dark matter-only simulation, the abundance of subhaloes in Apostle is reduced by 50 per cent near the centre and by 10 per cent within r200. In Auriga, the corresponding numbers are 80 per cent and 40 per cent. The velocity distributions of subhaloes are also affected by the presence of the galaxy, much more so in Auriga than in Apostle. The differences on subhalo properties in the two simulations can be traced back to the mass of the central galaxies, which in Auriga are typically twice as massive as those in Apostle. We show that some of the results from previous studies are inaccurate due to systematic errors in the modelling of subhalo orbits near the centre of haloes.
The Local Group galaxies offer some of the most discriminating tests of models of cosmic structure formation. For example, observations of the Milky Way (MW) and Andromeda satellite populations ...appear to be in disagreement with N-body simulations of the ‘lambda cold dark matter’ (ΛCDM) model: there are far fewer satellite galaxies than substructures in CDM haloes (the ‘missing satellites’ problem); dwarf galaxies seem to avoid the most massive substructures (the ‘too-big-to-fail’ problem); and the brightest satellites appear to orbit their host galaxies on a thin plane (the ‘planes of satellites’ problem). Here we present results from apostle (A Project Of Simulating The Local Environment), a suite of cosmological hydrodynamic simulations of 12 volumes selected to match the kinematics of the Local Group (LG) members. Applying the eagle code to the LG environment, we find that our simulations match the observed abundance of LG galaxies, including the satellite galaxies of the MW and Andromeda. Due to changes to the structure of haloes and the evolution in the LG environment, the simulations reproduce the observed relation between stellar mass and velocity dispersion of individual dwarf spheroidal galaxies without necessitating the formation of cores in their dark matter profiles. Satellite systems form with a range of spatial anisotropies, including one similar to the MWs, confirming that such a configuration is not unexpected in ΛCDM. Finally, based on the observed velocity dispersion, size, and stellar mass, we provide estimates of the maximum circular velocity for the haloes of nine MW dwarf spheroidals.
How unusual is the Milky Way’s assembly history? Evans, Tilly A; Fattahi, Azadeh; Deason, Alis J ...
Monthly notices of the Royal Astronomical Society,
10/2020, Letnik:
497, Številka:
4
Journal Article
Recenzirano
Odprti dostop
ABSTRACT
In the lambda cold dark matter (ΛCDM) model of structure formation galactic haloes build-up by accretion of mass and mergers of smaller haloes. The most recent massive merger event ...experienced by the Milky Way (MW) halo was the accretion of the Large Magellanic Cloud (LMC; which has a stellar mass of ∼109M⊙). Recent analyses of galactic stellar data from the Gaia satellite have uncovered an earlier massive accretion event, the Gaia-Enceladus Sausage (GES), which merged with the MW around 10 Gyr ago. Here, we use the EAGLE cosmological hydrodynamics simulation to study properties of simulated MW-mass haloes constrained to have accretion histories similar to that of the MW, specifically the recent accretion of an ‘LMC’ galaxy and a ‘GES’ merger, with a quiescent period between the GES merger and the infall of the LMC (the ‘LMC and GES’ category). We find that ∼16 per cent of MW-mass haloes have an LMC; ∼5 per cent have a GES event and no further merger with an equally massive object since z = 1; and only 0.65 per cent belong to the LMC and GES category. The progenitors of the MWs in this last category are much less massive than average at early times but eventually catch up with the mean. The LMC and GES category of galaxies naturally end up in the ‘blue cloud’ in the colour–magnitude diagram at z = 0, tend to have a disc morphology and have a larger than average number of satellite galaxies.