ABSTRACT
We study the formation of planes of dwarf galaxies around Milky Way (MW)-mass haloes in the eagle galaxy formation simulation. We focus on satellite systems similar to the one in the MW: ...spatially thin or with a large fraction of members orbiting in the same plane. To characterize the latter, we introduce a robust method to identify the subsets of satellites that have the most coplanar orbits. Out of the 11 MW classical dwarf satellites, 8 have highly clustered orbital planes whose poles are contained within a 22° opening angle centred around (l, b) = (182°, −2°). This configuration stands out when compared to both isotropic and typical ΛCDM satellite distributions. Purely flattened satellite systems are short-lived chance associations and persist for less than $1\, \rm {Gyr}$. In contrast, satellite subsets that share roughly the same orbital plane are longer lived, with half of the MW-like systems being at least $4\, \rm {Gyr}$ old. On average, satellite systems were flatter in the past, with a minimum in their minor-to-major axes ratio about $9\, \rm {Gyr}$ ago, which is the typical infall time of the classical satellites. MW-like satellite distributions have on average always been flatter than the overall population of satellites in MW-mass haloes and, in particular, they correspond to systems with a high degree of anisotropic accretion of satellites. We also show that torques induced by the aspherical mass distribution of the host halo channel some satellite orbits into the host’s equatorial plane, enhancing the fraction of satellites with coplanar orbits. In fact, the orbital poles of coplanar satellites are tightly aligned with the minor axis of the host halo.
The total satellite population of the Milky Way Newton, Oliver; Cautun, Marius; Jenkins, Adrian ...
Monthly notices of the Royal Astronomical Society,
09/2018, Letnik:
479, Številka:
3
Journal Article
The development of methods and algorithms to solve the N-body problem for classical, collisionless, non-relativistic particles has made it possible to follow the growth and evolution of cosmic dark ...matter structures over most of the universe’s history. In the best-studied case—the cold dark matter or CDM model—the dark matter is assumed to consist of elementary particles that had negligible thermal velocities at early times. Progress over the past three decades has led to a nearly complete description of the assembly, structure, and spatial distribution of dark matter haloes, and their substructure in this model, over almost the entire mass range of astronomical objects. On scales of galaxies and above, predictions from this standard CDM model have been shown to provide a remarkably good match to a wide variety of astronomical data over a large range of epochs, from the temperature structure of the cosmic background radiation to the large-scale distribution of galaxies. The frontier in this field has shifted to the relatively unexplored subgalactic scales, the domain of the central regions of massive haloes, and that of low-mass haloes and subhaloes, where potentially fundamental questions remain. Answering them may require: (i) the effect of known but uncertain baryonic processes (involving gas and stars), and/or (ii) alternative models with new dark matter physics. Here we present a review of the field, focusing on our current understanding of dark matter structure from N-body simulations and on the challenges ahead.
ABSTRACT
We examine the formation of dark matter (DM) cores in dwarf galaxies simulated with the eagle model of galaxy formation. As in earlier work, we find that the star formation (SF) gas density ...threshold (ρth) plays a critical role. At low thresholds (LT), gas is unable to reach densities high enough to dominate the gravitational potential before being dispersed by feedback from supernovae. LT runs show little effect on the inner DM profile, even in systems with extended and bursty SF, two ingredients often cited as critical for core formation. For higher thresholds, gas is able to dominate the gravitational potential before being ejected by feedback. This can lead to a substantial reduction in the inner DM content, but only if the gas is gravitationally important over an extended period of time, allowing the halo to contract before gas removal. Rapid assembly and removal of gas in short SF bursts is less effective at altering the inner DM content. Subsequent gas accretion may draw DM back in and reform a cusp, unless SF is bursty enough to prevent it, preserving the core. Thus, for the eagle SF + feedback model, there is no simple relation between core formation and SF history, contrary to recent claims. The dependence of the inner DM content of dwarfs on ρth hinders robust predictions and the interpretation of observations. A simulation of a $(12 \rm \ Mpc)^3$ volume with high ρth results in dwarfs with sizeable cores over a limited halo mass range, but with insufficient variety in mass profiles to explain the observed diversity of dwarf galaxy rotation curves.
N-body simulations suggest that the substructures that survive inside dark matter haloes follow universal distributions in mass and radial number density. We demonstrate that a simple analytical ...model can explain these subhalo distributions as resulting from tidal stripping which increasingly reduces the mass of subhaloes with decreasing halocentric distance. As a starting point, the spatial distribution of subhaloes of any given infall mass is shown to be largely indistinguishable from the overall mass distribution of the host halo. Using a physically motivated statistical description of the amount of mass stripped from individual subhaloes, the model fully describes the joint distribution of subhaloes in final mass, infall mass and radius. As a result, it can be used to predict several derived distributions involving combinations of these quantities including, but not limited to, the universal subhalo mass function, the subhalo spatial distribution, the gravitational lensing profile, the dark matter annihilation radiation profile and boost factor. This model clarifies a common confusion when comparing the spatial distributions of galaxies and subhaloes, the so-called anti-bias, as a simple selection effect. We provide a python code subgen for populating haloes with subhaloes at http://icc.dur.ac.uk/data/.
Galaxies fall into two clearly distinct types: 'blue-sequence' galaxies which are rapidly forming young stars, and 'red-sequence' galaxies in which star formation has almost completely ceased. Most ...galaxies more massive than ... follow the red sequence, while less massive central galaxies lie on the blue sequence. We show that these sequences are created by a competition between star formation-driven outflows and gas accretion on to the supermassive black hole at the galaxy's centre. We develop a simple analytic model for this interaction. In galaxies less massive than ..., young stars and supernovae drive a high-entropy outflow which is more buoyant than any tenuous corona. The outflow balances the rate of gas inflow, preventing high gas densities building up in the central regions. More massive galaxies, however, are surrounded by an increasingly hot corona. Above a halo mass of ..., the outflow ceases to be buoyant and star formation is unable to prevent the build-up of gas in the central regions. This triggers a strongly non-linear response from the black hole. Its accretion rate rises rapidly, heating the galaxy's corona, disrupting the incoming supply of cool gas and starving the galaxy of the fuel for star formation. The host galaxy makes a transition to the red sequence, and further growth predominantly occurs through galaxy mergers. We show that the analytic model provides a good description of galaxy evolution in the EAGLE hydrodynamic simulations. So long as star formation-driven outflows are present, the transition mass scale is almost independent of subgrid parameter choice. (ProQuest: ... denotes formulae/symbols omitted.)
Abstract
We estimate the systemic orbital kinematics of the Milky Way classical satellites and compare them with predictions from the Λ cold dark matter (ΛCDM) model derived from a semi-analytical ...galaxy formation model applied to high-resolution cosmological N-body simulations. We find that the Galactic satellite system is atypical of ΛCDM systems. The subset of 10 Galactic satellites with proper motion measurements has a velocity anisotropy, β = −2.2 ± 0.4, which lies in the 2.9 per cent tail of the ΛCDM distribution. Individually, the Milky Way satellites have radial velocities that are lower than expected for their proper motions, with 9 out of the 10 having at most 20 per cent of their orbital kinetic energy invested in radial motion. Such extreme values are expected in only 1.5 per cent of ΛCDM satellites systems. In the standard cosmological model, this tangential motion excess is unrelated to the existence of a Galactic ‘disc of satellites’. We present theoretical predictions for larger satellite samples that may become available as more proper motion measurements are obtained.
We use a suite of cosmological simulations to study the mass–concentration–redshift relation, c(M, z), of dark matter haloes. Our simulations include standard Λ-cold dark matter (CDM) models, and ...additional runs with truncated power spectra, consistent with a thermal warm dark matter (WDM) scenario. We find that the mass profiles of CDM and WDM haloes are self-similar and well approximated by the Einasto profile. The c(M, z) relation of CDM haloes is monotonic: concentrations decrease with increasing virial mass at fixed redshift, and decrease with increasing redshift at fixed mass. The mass accretion histories (MAHs) of CDM haloes are also scale-free, and can be used to infer concentrations directly. These results do not apply to WDM haloes: their MAHs are not scale-free because of the characteristic scale imposed by the power spectrum suppression. Further, the WDM c(M, z) relation is non-monotonic: concentrations peak at a mass scale dictated by the truncation scale, and decrease at higher and lower masses. We show that the assembly history of a halo can still be used to infer its concentration, provided that the total mass of its progenitors is considered (the ‘collapsed mass history’; CMH), rather than just that of its main ancestor. This exploits the scale-free nature of CMHs to derive a simple scaling that reproduces the mass–concentration–redshift relation of both CDM and WDM haloes over a vast range of halo masses and redshifts. Our model therefore provides a robust account of the mass, redshift, cosmology and power spectrum dependence of dark matter halo concentrations.
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.