ABSTRACT
We make the most precise determination to date of the number density of extragalactic 21-cm radio sources as a function of their spectral line widths – the H i velocity width function (H i ...WF) – based on 21 827 sources from the final $7000\, \mathrm{deg}^2$ data release of the Arecibo Legacy Fast ALFA (ALFALFA) survey. The number density of sources as a function of their neutral hydrogen masses – the H i mass function (H i MF) – has previously been reported to have a significantly different low-mass slope and ‘knee mass’ in the two sky regions surveyed during ALFALFA. In contrast with this, we find that the shape of the H i WF in the same two sky regions is remarkably similar, consistent with being identical within the confidence intervals implied by the data (but the overall normalization differs). The spatial uniformity of the H i WF implies that it is likely a stable tracer of the mass function of dark matter haloes, in spite of the environmental processes to which the measured variation in the H i MF are attributed, at least for galaxies containing enough neutral hydrogen to be detected. This insensitivity of the H i WF to galaxy formation and evolution can be exploited to turn it into a powerful constraint on cosmological models as future surveys yield increasingly precise measurements. We also report on the possible influence of a previously overlooked systematic error affecting the H i WF, which may plausibly see its low-velocity slope steepen by ∼40 per cent in analyses of future, deeper surveys. Finally, we provide an updated estimate of the ALFALFA completeness limit.
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.
We measure the star formation quenching efficiency and time-scale in cluster environments. Our method uses N-body simulations to estimate the probability distribution of possible orbits for a sample ...of observed Sloan Digital Sky Survey galaxies in and around clusters based on their position and velocity offsets from their host cluster. We study the relationship between their star formation rates and their likely orbital histories via a simple model in which star formation is quenched once a delay time after infall has elapsed. Our orbit library method is designed to isolate the environmental effect on the star formation rate due to a galaxy's present-day host cluster from 'pre-processing' in previous group hosts. We find that quenching of satellite galaxies of all stellar masses in our sample ... by massive ... clusters is essentially 100 per cent efficient. Our fits show that all galaxies quench on their first infall, approximately at or within a Gyr of their first pericentric passage. There is little variation in the onset of quenching from galaxy-to-galaxy: the spread in this time is at most ~2 Gyr at fixed M*. Higher mass satellites quench earlier, with very little dependence on host cluster mass in the range probed by our sample. (ProQuest: ... denotes formulae/symbols omitted.)
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
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.
Abstract
The shallow faint-end slope of the galaxy mass function is usually reproduced in Λ cold dark matter (ΛCDM) galaxy formation models by assuming that the fraction of baryons that turn into ...stars drops steeply with decreasing halo mass and essentially vanishes in haloes with maximum circular velocities Vmax < 20–30 km s−1. Dark-matter-dominated dwarfs should therefore have characteristic velocities of about that value, unless they are small enough to probe only the rising part of the halo circular velocity curve (i.e. half-mass radii, r1/2 ≪ 1 kpc). Many dwarfs have properties in disagreement with this prediction: they are large enough to probe their halo Vmax but their characteristic velocities are well below 20 km s−1. These ‘cold faint giants’ (an extreme example is the recently discovered Crater 2 Milky Way satellite) can only be reconciled with our ΛCDM models if they are the remnants of once massive objects heavily affected by tidal stripping. We examine this possibility using the APOSTLE cosmological hydrodynamical simulations of the Local Group. Assuming that low-velocity-dispersion satellites have been affected by stripping, we infer their progenitor masses, radii, and velocity dispersions, and find them in remarkable agreement with those of isolated dwarfs. Tidal stripping also explains the large scatter in the mass discrepancy–acceleration relation in the dwarf galaxy regime: tides remove preferentially dark matter from satellite galaxies, lowering their accelerations below the amin ∼ 10−11 m s−2 minimum expected for isolated dwarfs. In many cases, the resulting velocity dispersions are inconsistent with the predictions from Modified Newtonian Dynamics, a result that poses a possibly insurmountable challenge to that scenario.
ABSTRACT
We combine orbital information from N-body simulations with an analytic model for star formation quenching and SDSS observations to infer the differential effect of the group/cluster ...environment on star formation in satellite galaxies. We also consider a model for gas stripping, using the same input supplemented with H i fluxes from the ALFALFA survey. The models are motivated by and tested on the Hydrangea cosmological hydrodynamical simulation suite. We recover the characteristic times when satellite galaxies are stripped and quenched. Stripping in massive ($M_{\rm vir}\sim 10^{14.5}\, {\rm M}_\odot$) clusters typically occurs at or just before the first pericentric passage. Lower mass ($\sim 10^{13.5}\, {\rm M}_\odot$) groups strip their satellites on a significantly longer (by $\sim 3\, {\rm Gyr}$) time-scale. Quenching occurs later: Balmer emission lines typically fade $\sim 3.5\, {\rm Gyr}$ ($5.5\, {\rm Gyr}$) after first pericentre in clusters (groups), followed a few hundred Myr later by reddenning in (g − r) colour. These ‘delay time-scales’ are remarkably constant across the entire satellite stellar mass range probed (∼109.5–$10^{11}\, {\rm M}_\odot$), a feature closely tied to our treatment of ‘group pre-processing’. The lowest mass groups in our sample ($\sim 10^{12.5}\, {\rm M}_\odot$) strip and quench their satellites extremely inefficiently: typical time-scales may approach the age of the Universe. Our measurements are qualitatively consistent with the ‘delayed-then-rapid’ quenching scenario advocated for by several other studies, but we find significantly longer delay times. Our combination of a homogeneous analysis and input catalogues yields new insight into the sequence of events leading to quenching across wide intervals in host and satellite mass.
We examine the spatial distribution of the oldest and most metal-poor stellar populations of Milky Way-sized galaxies using the A Project Of Simulating The Local Environment (APOSTLE) cosmological ...hydrodynamical simulations of the Local Group. In agreement with earlier work, we find strong radial gradients in the fraction of the oldest (t sub( form) < 0.8 Gyr) and most metal-poor (Fe/H < -2.5) stars, both of which increase outwards. The most metal-poor stars form over an extended period of time; half of them form after z = 5.3, and the last 10 per cent after z = 2.8. The age of the metal-poor stellar population also shows significant variation with environment; a high fraction of them are old in the galaxy's central regions and an even higher fraction in some individual dwarf galaxies, with substantial scatter from dwarf to dwarf. We investigate the dependence of these results on the assumptions made for metal mixing. Overall, over half of the stars that belong to both the oldest and most metal-poor population are found outside the solar circle. Somewhat counter-intuitively, we find that dwarf galaxies with a large fraction of metal-poor stars that are very old are systems where metal-poor stars are relatively rare, but where a substantial old population is present. Our results provide guidance for interpreting the results of surveys designed to hunt for the earliest and most pristine stellar component of our Milky Way.