Abstract
Merging galaxy clusters have been touted as one of the best probes for constraining self-interacting dark matter, but few simulations exist to back up this claim. We simulate equal-mass ...mergers of 1015 M⊙ haloes, like the El Gordo and Sausage clusters, with cosmologically motivated halo and merger parameters, and with velocity-independent dark-matter self-interactions. Although the standard lore for merging clusters is that self-interactions lead to large separations between the galaxy and dark-matter distributions, we find that maximal galaxy–dark matter offsets of ≲20 kpc form for a self-interaction cross-section of σSI/mχ = 1 cm2 g−1. This is an order of magnitude smaller than those measured in observed equal-mass and near-equal-mass mergers, and is likely to be even smaller for lower mass systems. While competitive cross-section constraints are thus unlikely to emerge from offsets, we find other signatures of self-interactions that are more promising. Intriguingly, we find that after dark-matter haloes coalesce, the collisionless galaxies and especially the brightest cluster galaxy (BCG) oscillate around the centre of the merger remnant on stable orbits of 100 kpc for σSI/mχ = 1 cm2 g−1 for at least several Gyr, well after the clusters have relaxed. If BCG miscentring in relaxed clusters remains a robust prediction of self-interacting dark matter under the addition of gas physics, substructure, merger mass ratios (e.g. 10:1 like the Bullet Cluster) and complex cosmological merger histories, the observed BCG offsets may constrain σSI/mχ to ≲0.1 cm2 g−1 – the tightest constraint yet.
A critical challenge to the cold dark matter (CDM) paradigm is that there are fewer satellites observed around the Milky Way than found in simulations of dark matter substructure. We show that there ...is a match between the observed satellite counts corrected by the detection efficiency of the Sloan Digital Sky Survey (for luminosities L≳340 L_{⊙}) and the number of luminous satellites predicted by CDM, assuming an empirical relation between stellar mass and halo mass. The "missing satellites problem," cast in terms of number counts, is thus solved. We also show that warm dark matter models with a thermal relic mass smaller than 4 keV are in tension with satellite counts, putting pressure on the sterile neutrino interpretation of recent x-ray observations. Importantly, the total number of Milky Way satellites depends sensitively on the spatial distribution of satellites, possibly leading to a "too many satellites" problem. Measurements of completely dark halos below 10^{8} M_{⊙}, achievable with substructure lensing and stellar stream perturbations, are the next frontier for tests of CDM.
ABSTRACT
In the standard Lambda cold dark matter paradigm, pure dark matter simulations predict dwarf galaxies should inhabit dark matter haloes with a centrally diverging density ‘cusp’. This is in ...conflict with observations that typically favour a constant density ‘core’. We investigate this ‘cusp-core problem’ in ‘ultra-faint’ dwarf galaxies simulated as part of the ‘Engineering Dwarfs at Galaxy formation’s Edge’ project. We find, similarly to previous work, that gravitational potential fluctuations within the central region of the simulated dwarfs kinematically heat the dark matter particles, lowering the dwarfs’ central dark matter density. However, these fluctuations are not exclusively caused by gas inflow/outflow, but also by impulsive heating from minor mergers. We use the genetic modification approach on one of our dwarf’s initial conditions to show how a delayed assembly history leads to more late minor mergers and, correspondingly, more dark matter heating. This provides a mechanism by which even ultra-faint dwarfs ($M_* \lt 10^5\, \text{M}_{\odot }$), in which star formation was fully quenched at high redshift, can have their central dark matter density lowered over time. In contrast, we find that late major mergers can regenerate a central dark matter cusp, if the merging galaxy had sufficiently little star formation. The combination of these effects leads us to predict significant stochasticity in the central dark matter density slopes of the smallest dwarfs, driven by their unique star formation and mass assembly histories.
Abstract
We present the first detailed comparison of populations of dwarf galaxy stellar streams in cosmological simulations and the Milky Way. In particular, we compare streams identified around 13 ...Milky Way analogs in the FIRE-2 simulations to streams observed by the Southern Stellar Stream Spectroscopic Survey (
S
5
). For an accurate comparison, we produce mock Dark Energy Survey (DES) observations of the FIRE streams and estimate the detectability of their tidal tails and progenitors. The number and stellar mass distributions of detectable stellar streams is consistent between observations and simulations. However, there are discrepancies in the distributions of pericenters and apocenters, with the detectable FIRE streams, on average, forming at larger pericenters (out to >110 kpc) and surviving only at larger apocenters (≳40 kpc) than those observed in the Milky Way. We find that the population of high-stellar-mass dwarf galaxy streams in the Milky Way is incomplete. Interestingly, a large fraction of the FIRE streams would only be detected as intact satellites in DES-like observations, since their tidal tails have too low surface brightness to be detectable. We thus predict a population of yet-undetected tidal tails around Milky Way satellites, as well as a population of fully undetected low-surface-brightness stellar streams, and estimate their detectability with the Rubin Observatory. Finally, we discuss the causes and implications of the discrepancies between the stream populations in FIRE and the Milky Way, and explore future avenues for tests of satellite disruption in cosmological simulations.
Abstract
We build a statistical framework to infer the global properties of the satellite system of the Andromeda galaxy (M31) from the properties of individual dwarf galaxies located in the ...Pan-Andromeda Archaelogical Survey (PAndAS) and the previously determined completeness of the survey. Using forward modeling, we infer the slope of the luminosity function of the satellite system, the slope of its spatial density distribution, and the size–luminosity relation followed by the dwarf galaxies. We find that the slope of the luminosity function is
β
= −1.5 ± 0.1. Combined with the spatial density profile, it implies that, when accounting for survey incompleteness, M31 hosts
92
−
26
+
19
dwarf galaxies with
M
V
< −5.5 and a sky-projected distance from M31 between 30 and 300 kpc. We conclude that many faint or distant dwarf galaxies remain to be discovered around Andromeda, especially outside the PAndAS footprint. Finally, we use our model to test if the higher number of satellites situated in the hemisphere facing the Milky Way could be explained simply by the detection limits of dwarf galaxy searches. We rule this out at >99.9% confidence and conclude that this anisotropy is an intrinsic feature of the M31 satellite system. The statistical framework we present here is a powerful tool to robustly constrain the properties of a satellite system and compare those across hosts, especially considering the upcoming start of the Euclid or Rubin large photometric surveys that are expected to uncover a large number of dwarf galaxies in the Local Volume.
ABSTRACT
We show how the interplay between feedback and mass-growth histories introduces scatter in the relationship between stellar and neutral gas properties of field faint dwarf galaxies ...($M_{\star }\lessapprox 10^{6} \, \mbox{M}_\mathrm{\odot }$). Across a suite of cosmological, high-resolution zoomed simulations, we find that dwarf galaxies of stellar masses $10^5 \le M_{\star }\le 10^{6} \, \mbox{M}_\mathrm{\odot }$ are bimodal in their cold gas content, being either H i-rich or H i-deficient. This bimodality is generated through the coupling between (i) the modulation of H i contents by the background of ultraviolet radiation (UVB) at late times and (ii) the significant scatter in the stellar-to-halo mass relationship induced by reionization. Furthermore, our H i-rich dwarfs exhibit disturbed and time-variable neutral gas distributions primarily due to stellar feedback. Over the last four billion years, we observe order-of-magnitude changes around the median $M_{\mathrm{H\,\small {I} }}$, factor-of-a-few variations in H i spatial extents, and spatial offsets between H i and stellar components regularly exceeding the galaxies’ optical sizes. Time variability introduces further scatter in the $M_{\star }\!-\! M_{\mathrm{H\,\small {I} }}$ relation and affects a galaxy’s detectability in H i at any given time. These effects will need to be accounted for when interpreting observations of the population of faint, H i-bearing dwarfs by the combination of optical and radio wide, deep surveys.
ABSTRACT
The Eridanus II (EriII) ‘ultra-faint’ dwarf has a large (15 pc) and low-mass (4.3 × 103 M⊙) star cluster (SC) offset from its centre by 23 ± 3 pc in projection. Its size and offset are ...naturally explained if EriII has a central dark matter core, but such a core may be challenging to explain in a ΛCDM cosmology. In this paper, we revisit the survival and evolution of EriII’s SC, focusing for the first time on its puzzlingly large ellipticity ($0.31^{+0.05}_{-0.06}$). We perform a suite of 960 direct N-body simulations of SCs, orbiting within a range of spherical background potentials fit to ultra-faint dwarf (UFD) galaxy simulations. We find only two scenarios that come close to explaining EriII’s SC. In the first scenario, EriII has a low-density dark matter core (of size ${\sim}70\, \text{pc}$ and density $\lesssim 2\times 10^8\, \text{M}_{\odot }\, \text{kpc}^{-3}$). In this model, the high ellipticity of EriII’s SC is set at birth, with the lack of tidal forces in the core allowing its ellipticity to remain frozen for long times. In the second scenario, EriII’s SC orbits in a partial core, with its high ellipticity owing to its imminent tidal destruction. However, this latter model struggles to reproduce the large size of EriII’s SC, and it predicts substantial tidal tails around EriII’s SC that should have already been seen in the data. This leads us to favour the cored model. We discuss potential caveats to these findings, and the implications of the cored model for galaxy formation and the nature of dark matter.
ABSTRACT
Collisionless dark matter only (DMO) structure formation simulations predict that dark matter (DM) haloes are prolate in their centres and triaxial towards their outskirts. The addition of ...gas condensation transforms the central DM shape to be rounder and more oblate. It is not clear, however, whether such shape transformations occur in ‘ultra-faint’ dwarfs, which have extremely low baryon fractions. We present the first study of the shape and velocity anisotropy of ultra-faint dwarf galaxies that have gas mass fractions of fgas(r < Rhalf) < 0.06. These dwarfs are drawn from the Engineering Dwarfs at Galaxy formation’s Edge (EDGE) project, using high-resolution simulations that allow us to resolve DM halo shapes within the half-light radius (∼100 pc). We show that gas-poor ultra-faints (M200c ≤ 1.5 × 109 M⊙; fgas < 10−5) retain their pristine prolate DM halo shape even when gas, star formation, and feedback are included. This could provide a new and robust test of DM models. By contrast, gas-rich ultra-faints (M200c > 3 × 109 M⊙; fgas > 10−4) become rounder and more oblate within ∼10 half-light radii. Finally, we find that most of our simulated dwarfs have significant radial velocity anisotropy that rises to $\tilde{\beta } \gt 0.5$ at R ≳ 3Rhalf. The one exception is a dwarf that forms a rotating gas/stellar disc because of a planar, major merger. Such strong anisotropy should be taken into account when building mass models of gas-poor ultra-faints.
ABSTRACT
Ultra-faint dwarf galaxies (UFDs) are commonly found in close proximity to the Milky Way and other massive spiral galaxies. As such, their projected stellar ellipticity and extended light ...distributions are often thought to owe to tidal forces. In this paper, we study the projected stellar ellipticities and faint stellar outskirts of tidally isolated ultra-faints drawn from the ‘Engineering Dwarfs at Galaxy Formation’s Edge’ (EDGE) cosmological simulation suite. Despite their tidal isolation, our simulated dwarfs exhibit a wide range of projected ellipticities (0.03 < ε < 0.85), with many possessing anisotropic extended stellar haloes that mimic tidal tails, but owe instead to late-time accretion of lower mass companions. Furthermore, we find a strong causal relationship between ellipticity and formation time of a UFD, which is robust to a wide variation in the feedback model. We show that the distribution of projected ellipticities in our suite of simulated EDGE dwarfs matches well with a sample of 19 Local Group dwarf galaxies and a sample of 11 isolated dwarf galaxies. Given ellipticity in EDGE arises from an ex-situ accretion origin, the agreement in shape indicates the ellipticities of some observed dwarfs may also originate from a non-tidal scenario. The orbital parameters of these observed dwarfs further support that they are not currently tidally disrupting. If the baryonic content in these galaxies is still tidally intact, then the same may be true for their dark matter content, making these galaxies in our Local Group pristine laboratories for testing dark matter and galaxy formation models.
ABSTRACT
Low-mass dwarf galaxies are expected to reside within dark matter haloes that have a pristine, ‘cuspy’ density profile within their stellar half-light radii. This is because they form too ...few stars to significantly drive dark matter heating through supernova-driven outflows. Here, we study such simulated faint systems ($10^4 \le M_{\star }\le 2 \times 10^6 \, \mbox{M}_\mathrm{\odot }$) drawn from high-resolution (3 pc) cosmological simulations from the ‘Engineering Dwarf Galaxies at the Edge of galaxy formation’ (EDGE) project. We confirm that these objects have steep and rising inner dark matter density profiles at z = 0, little affected by galaxy formation effects. But five dwarf galaxies from the suite also showcase a detectable H i reservoir ($M_{\mathrm{H\, {\small I} }}\approx 10^{5}-10^{6} \, \mbox{M}_\mathrm{\odot }$), analogous to the observed population of faint, H i-bearing dwarf galaxies. These reservoirs exhibit episodes of ordered rotation, opening windows for rotation curve analysis. Within actively star-forming dwarfs, stellar feedback easily disrupts the tenuous H i discs ($v_{\phi , g}\approx 10\, \mathrm{km} \, \mathrm{s}^{-1}$), making rotation short-lived ($\ll 150 \, \mathrm{Myr}$) and more challenging to interpret for dark matter inferences. In contrast, we highlight a long-lived ($\ge 500 \, \mathrm{Myr}$) and easy-to-interpret H i rotation curve extending to $\approx 2\, r_{1/2, \text{3D}}$ in a quiescent dwarf, that has not formed new stars since z = 4. This stable gas disc is supported by an oblate dark matter halo shape that drives high-angular momentum gas flows. Our results strongly motivate further searches for H i in rotation curves in the observed population of H i-bearing low-mass dwarfs, that provide a key regime to disentangle the respective roles of dark matter microphysics and galaxy formation effects in driving dark matter heating.