Kyle A. Oman

The APOSTLE project: Local Group kinematic mass constraints and simulation candidate selection

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 × 1012M . 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 × 1012M for the APOSTLE project (A Project Of Simulating The Local Environment), a suite of hydrodynamical resimulations at various numerical resolution levels (reaching up to ∼104M 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.

Missing dark matter in dwarf galaxies

We use cosmological hydrodynamical simulations of the APOSTLE project to examine the fraction of baryons in CDM haloes that collect into galaxies. This ‘galaxy formation efficiency’ correlates strongly and with little scatter with halo mass, dropping steadily towards dwarf galaxies. The baryonic mass of a galaxy may thus be used to place a lower limit on total halo mass and, consequently, on its asymptotic maximum circular velocity. A number of dwarfs seem to violate this constraint, having baryonic masses up to ten times higher than expected from their rotation speeds, or, alternatively, rotating at only half the speed expected for their mass. Taking the data at face value, either these systems have formed galaxies with extraordinary efficiency ‐ highly unlikely given their shallow potential wells ‐ or they inhabit haloes with extreme deficits in their dark matter content. This ‘missing dark matter’ is reminiscent of the inner mass deficits of galaxies with slowly-rising rotation curves, but extends to regions larger than the luminous galaxies themselves, disfavouring explanations based on star formationinduced ‘cores’ in the dark matter. An alternative could be that galaxy inclination errors have been underestimated, and that these are just systems where inferred mass profiles have been compromised by systematic uncertainties in interpreting the velocity field. This should be investigated further, since it might provide a simple explanation not only for missing-dark-matter galaxies but also for other challenges to our understanding of the inner structure of cold dark matter haloes.

Shaken and stirred: the Milky Way's dark substructures

The predicted abundance and properties of the low-mass substructures embedded inside larger dark matter haloes differ sharply among alternative dark matter models. Too small to host galaxies themselves, these subhaloes may still be detected via gravitational lensing or via perturbations of the Milky Ways globular cluster streams and its stellar disc. Here, we use the Apostle cosmological simulations to predict the abundance and the spatial and velocity distributions of subhaloes in the range 106.5–108.5 M⊙ inside haloes of mass ∼1012 M⊙ in Λ cold dark matter. Although these subhaloes are themselves devoid of baryons, we find that baryonic effects are important. Compared to corresponding dark matter only simulations, the loss of baryons from subhaloes and stronger tidal disruption due to the presence of baryons near the centre of the main halo reduce the number of subhaloes by ∼1/4 to 1/2, independently of subhalo mass, but increasingly towards the host halo centre. We also find that subhaloes have non-Maxwellian orbital velocity distributions, with centrally rising velocity anisotropy and positive velocity bias that reduces the number of low-velocity subhaloes, particularly near the halo centre. We parametrize the predicted population of subhaloes in terms of mass, galactocentric distance and velocities. We discuss implications of our results for the prospects of detecting dark matter substructures and for possible inferences about the nature of dark matter.

The low-mass end of the baryonic Tully–Fisher relation

The scaling of disc galaxy rotation velocity with baryonic mass (the ‘baryonic Tully–Fisher’ relation, BTF) has long confounded galaxy formation models. It is steeper than the M ∝ V3 scaling relating halo virial masses and circular velocities and its zero-point implies that galaxies comprise a very small fraction of available baryons. Such low galaxy formation efficiencies may, in principle, be explained by winds driven by evolving stars, but the tightness of the BTF relation argues against the substantial scatter expected from such a vigorous feedback mechanism. We use the APOSTLE/EAGLE simulations to show that the BTF relation is well reproduced in _cold dark matter (CDM) simulations that match the size and number of galaxies as a function of stellarmass. In such models, galaxy rotation velocities are proportional to halo virial velocity and the steep velocity-mass dependence results from the decline in galaxy formation efficiency with decreasing halo mass needed to reconcile the CDM halo mass function with the galaxy luminosity function. The scatter in the simulated BTF is smaller than observed, even when considering all simulated galaxies and not just rotationally supported ones. The simulations predict that the BTF should become increasingly steep at the faint end, although the velocity scatter at fixed mass should remain small. Observed galaxies with rotation speeds below ∼40 km s−1 seem to deviate from this prediction. We discuss observational biases and modelling uncertainties that may help to explain this disagreement in the context of _CDM models of dwarf galaxy formation.

Knowing the unknowns: uncertainties in simple estimators of galactic dynamical masses

The observed stellar kinematics of dispersion-supported galaxies are often used to measure dynamical masses. Recently, several analytical relationships between the stellar line-of-sight velocity dispersion, the projected (2D) or deprojected (3D) half-light radius and the total mass enclosed within the half-light radius, relying on the spherical Jeans equation, have been proposed. Here, we use the APOSTLE cosmological hydrodynamical simulations of the Local Group to test the validity and accuracy of such mass estimators for both dispersion and rotation-supported galaxies, for field and satellite galaxies, and for galaxies of varying masses, shapes and velocity dispersion anisotropies. We find that the mass estimators of Walker et al. and Wolf et al. are able to recover the masses of dispersion-dominated systems with little systematic bias, but with a 1σ scatter of 25 and 23 per cent, respectively. The error on the estimated mass is dominated by the impact of the 3D shape of the stellar mass distribution, which is difficult to constrain observationally. This intrinsic scatter becomes the dominant source of uncertainty in the masses estimated for galaxies like the dwarf spheroidal (dSph) satellites of the Milky Way, where the observational errors in their sizes and velocity dispersions are small. Such scatter may also affect the inner density slopes of dSphs derived from multiple stellar populations, relaxing the significance with which Navarro–Frenk–White profiles may be excluded, depending on the degree to which the relevant properties of the different stellar populations are correlated. Finally, we derive a new optimal mass estimator that removes the residual biases and achieves a statistically significant reduction in the scatter to 20 per cent overall for dispersion-dominated galaxies, allowing more precise and accurate mass estimates.

The origin of the mass discrepancy-acceleration relation in ΛCDM.

We examine the origin of the mass discrepancy–radial acceleration relation (MDAR) of disc galaxies. This is a tight empirical correlation between the disc centripetal acceleration and that expected from the baryonic component. The MDAR holds for most radii probed by disc kinematic tracers, regardless of galaxy mass or surface brightness. The relation has two characteristic accelerations: a0, above which all galaxies are baryon dominated, and amin, an effective minimum acceleration probed by kinematic tracers in isolated galaxies. We use a simple model to show that these trends arise naturally in Λ cold dark matter (ΛCDM). This is because (i) disc galaxies in ΛCDM form at the centre of dark matter haloes spanning a relatively narrow range of virial mass; (ii) cold dark matter halo acceleration profiles are self-similar and have a broad maximum at the centre, reaching values bracketed precisely by amin and a0 in that mass range and (iii) halo mass and galaxy size scale relatively tightly with the baryonic mass of a galaxy in any successful ΛCDM galaxy formation model. Explaining the MDAR in ΛCDM does not require modifications to the cuspy inner mass profiles of dark haloes, although these may help to understand the detailed rotation curves of some dwarf galaxies and the origin of extreme outliers from the main relation. The MDAR is just a reflection of the self-similar nature of cold dark matter haloes and of the physical scales introduced by the galaxy formation process.

Mass-Discrepancy Acceleration Relation : A Natural Outcome of Galaxy Formation in Cold Dark Matter Halos

We analyze the total and baryonic acceleration profiles of a set of well-resolved galaxies identified in the eagle suite of hydrodynamic simulations. Our runs start from the same initial conditions but adopt different prescriptions for unresolved stellar and active galactic nuclei feedback, resulting in diverse populations of galaxies by the present day. Some of them reproduce observed galaxy scaling relations, while others do not. However, regardless of the feedback implementation, all of our galaxies follow closely a simple relationship between the total and baryonic acceleration profiles, consistent with recent observations of rotationally supported galaxies. The relation has small scatter: Different feedback implementations-which produce different galaxy populations-mainly shift galaxies along the relation rather than perpendicular to it. Furthermore, galaxies exhibit a characteristic acceleration g_{†}, above which baryons dominate the mass budget, as observed. These observations, consistent with simple modified Newtonian dynamics, can be accommodated within the standard cold dark matter paradigm.

Tidal stripping and the structure of dwarf galaxies in the Local Group.

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.

The core–cusp problem: a matter of perspective

The existence of two kinematically and chemically distinct stellar subpopulations in the Sculptor and Fornax dwarf galaxies offers the opportunity to constrain the density profile of their matter haloes by measuring the mass contained within the well-separated half-light radii of the two metallicity subpopulations. Walker and Penarrubia have used this approach to argue that data for these galaxies are consistent with constant-density ‘cores’ in their inner regions and rule out ‘cuspy’ Navarro–Frenk–White (NFW) profiles with high statistical significance, particularly in the case of Sculptor. We test the validity of these claims using dwarf galaxies in the APOSTLE (A Project Of Simulating The Local Environment) Λ cold dark matter cosmological hydrodynamic simulations of analogues of the Local Group. These galaxies all have NFW dark matter density profiles and a subset of them develop two distinct metallicity subpopulations reminiscent of Sculptor and Fornax. We apply a method analogous to that of Walker and Penarrubia to a sample of 50 simulated dwarfs and find that this procedure often leads to a statistically significant detection of a core in the profile when in reality there is a cusp. Although multiple factors contribute to these failures, the main cause is a violation of the assumption of spherical symmetry upon which the mass estimators are based. The stellar populations of the simulated dwarfs tend to be significantly elongated and, in several cases, the two metallicity populations have different asphericity and are misaligned. As a result, a wide range of slopes of the density profile are inferred depending on the angle from which the galaxy is viewed.

The oldest and most metal-poor stars in the APOSTLE Local Group simulations

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 (tform < 0.8 Gyr) and most metal-poor ([Fe/H]<−2.5) stars, both of which increase outwards. The most metalpoor 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 metalpoor 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.