skip to main content


Title: Dark and luminous satellites of LMC-mass galaxies in the FIRE simulations
ABSTRACT Within lambda cold dark matter ($\Lambda$CDM), dwarf galaxies like the Large Magellanic Cloud (LMC) are expected to host numerous dark matter subhaloes, several of which should host faint dwarf companions. Recent Gaia proper motions confirm new members of the LMC system in addition to the previously known SMC, including two classical dwarf galaxies ($M_\ast$$\gt 10^5$ M$_{\odot }$; Carina and Fornax) as well as several ultrafaint dwarfs (Car2, Car3, Hor1, and Hyd1). We use the Feedback In Realistic Environments (FIRE) simulations to study the dark and luminous (down to ultrafaint masses, $M_\ast$$\sim$6$\times 10^ {3}$ M$_{\odot }$) substructure population of isolated LMC-mass hosts ($M_{\text{200m}}$ = 1–3$\times 10^ {11}$ M$_{\odot }$) and place the Gaia  + DES results in a cosmological context. By comparing number counts of subhaloes in simulations with and without baryons, we find that, within 0.2 $r_{\text{200m}}$, LMC-mass hosts deplete $\sim$30 per cent of their substructure, significantly lower than the $\sim$70 per cent of substructure depleted by Milky Way (MW) mass hosts. For our highest resolution runs ($m_\text{bary}$  = 880 M$_{\odot }$), $\sim 5\!-\!10$ subhaloes form galaxies with $M_\ast$$\ge 10^{4}$ M$_{\odot }$ , in agreement with the seven observationally inferred pre-infall LMC companions. However, we find steeper simulated luminosity functions than observed, hinting at observation incompleteness at the faint end. The predicted DM content for classical satellites in FIRE agrees with observed estimates for Carina and Fornax, supporting the case for an LMC association. We predict that tidal stripping within the LMC potential lowers the inner dark matter density of ultrafaint companions of the LMC. Thus, in addition to their orbital consistency, the low densities of dwarfs Car2, Hyd1, and Hyd2 reinforce their likelihood of Magellanic association.  more » « less
Award ID(s):
1752913
PAR ID:
10174834
Author(s) / Creator(s):
; ; ; ; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
489
Issue:
4
ISSN:
0035-8711
Page Range / eLocation ID:
5348 to 5364
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    We study a suite of extremely high-resolution cosmological Feedback in Realistic Environments simulations of dwarf galaxies ($M_{\rm halo} \lesssim 10^{10}\rm \, M_{\odot }$), run to z = 0 with $30\, \mathrm{M}_{\odot }$ resolution, sufficient (for the first time) to resolve the internal structure of individual supernovae remnants within the cooling radius. Every halo with $M_{\rm halo} \gtrsim 10^{8.6}\, \mathrm{M}_{\odot }$ is populated by a resolved stellar galaxy, suggesting very low-mass dwarfs may be ubiquitous in the field. Our ultra-faint dwarfs (UFDs; $M_{\ast }\lt 10^{5}\, \mathrm{M}_{\odot }$) have their star formation (SF) truncated early (z ≳ 2), likely by reionization, while classical dwarfs ($M_{\ast }\gt 10^{5}\, \mathrm{M}_{\odot }$) continue forming stars to z < 0.5. The systems have bursty star formation histories, forming most of their stars in periods of elevated SF strongly clustered in both space and time. This allows our dwarf with M*/Mhalo > 10−4 to form a dark matter core ${\gt}200\rm \, pc$, while lower mass UFDs exhibit cusps down to ${\lesssim}100\rm \, pc$, as expected from energetic arguments. Our dwarfs with $M_{\ast }\gt 10^{4}\, \mathrm{M}_{\odot }$ have half-mass radii (R1/2) in agreement with Local Group (LG) dwarfs (dynamical mass versus R1/2 and stellar rotation also resemble observations). The lowest mass UFDs are below surface brightness limits of current surveys but are potentially visible in next-generation surveys (e.g. LSST). The stellar metallicities are lower than in LG dwarfs; this may reflect pre-enrichment of the LG by the massive hosts or Pop-III stars. Consistency with lower resolution studies implies that our simulations are numerically robust (for a given physical model).

     
    more » « less
  2. null (Ed.)
    ABSTRACT We present the first set of cosmological baryonic zoom-in simulations of galaxies including dissipative self-interacting dark matter (dSIDM). These simulations utilize the Feedback In Realistic Environments galaxy formation physics, but allow the dark matter to have dissipative self-interactions analogous to standard model forces, parametrized by the self-interaction cross-section per unit mass, (σ/m), and the dimensionless degree of dissipation, 0 < fdiss < 1. We survey this parameter space, including constant and velocity-dependent cross-sections, and focus on structural and kinematic properties of dwarf galaxies with $M_{\rm halo} \sim 10^{10-11}{\, \rm M_\odot }$ and $M_{\ast } \sim 10^{5-8}{\, \rm M_\odot }$. Central density profiles (parametrized as ρ ∝ rα) of simulated dwarfs become cuspy when $(\sigma /m)_{\rm eff} \gtrsim 0.1\, {\rm cm^{2}\, g^{-1}}$ (and fdiss = 0.5 as fiducial). The power-law slopes asymptote to α ≈ −1.5 in low-mass dwarfs independent of cross-section, which arises from a dark matter ‘cooling flow’. Through comparisons with dark matter only simulations, we find the profile in this regime is insensitive to the inclusion of baryons. However, when $(\sigma /m)_{\rm eff} \ll 0.1\, {\rm cm^{2}\, g^{-1}}$, baryonic effects can produce cored density profiles comparable to non-dissipative cold dark matter (CDM) runs but at smaller radii. Simulated galaxies with $(\sigma /m) \gtrsim 10\, {\rm cm^{2}\, g^{-1}}$ and the fiducial fdiss develop significant coherent rotation of dark matter, accompanied by halo deformation, but this is unlike the well-defined thin ‘dark discs’ often attributed to baryon-like dSIDM. The density profiles in this high cross-section model exhibit lower normalizations given the onset of halo deformation. For our surveyed dSIDM parameters, halo masses and galaxy stellar masses do not show appreciable difference from CDM, but dark matter kinematics and halo concentrations/shapes can differ. 
    more » « less
  3. ABSTRACT

    A variety of observational campaigns seek to test dark matter models by measuring dark matter subhaloes at low masses. Despite their predicted lack of stars, these subhaloes may be detectable through gravitational lensing or via their gravitational perturbations on stellar streams. To set measurable expectations for subhalo populations within Lambda cold dark matter, we examine 11 Milky Way (MW)-mass haloes from the FIRE-2 baryonic simulations, quantifying the counts and orbital fluxes for subhaloes with properties relevant to stellar stream interactions: masses down to $10^{6}\, \text{M}_\odot$, distances ≲50 kpc of the galactic centre, across z = 0 − 1 (tlookback = 0–8 Gyr). We provide fits to our results and their dependence on subhalo mass, distance, and lookback time, for use in (semi)analytical models. A typical MW-mass halo contains ≈16 subhaloes $\gt 10^{7}\, \text{M}_\odot$ (≈1 subhalo $\gt 10^{8}\, \text{M}_\odot$) within 50 kpc at z ≈ 0. We compare our results with dark matter-only versions of the same simulations: because they lack a central galaxy potential, they overpredict subhalo counts by 2–10×, more so at smaller distances. Subhalo counts around a given MW-mass galaxy declined over time, being ≈10× higher at z = 1 than at z ≈ 0. Subhaloes have nearly isotropic orbital velocity distributions at z ≈ 0. Across our simulations, we also identified 4 analogues of Large Magellanic Cloud satellite passages; these analogues enhance subhalo counts by 1.4–2.1 times, significantly increasing the expected subhalo population around the MW today. Our results imply an interaction rate of ∼5 per Gyr for a stream like GD-1, sufficient to make subhalo–stream interactions a promising method of measuring dark subhaloes.

     
    more » « less
  4. ABSTRACT

    We use analytical and N-body methods to study the capture of field stars by gravitating substructures moving across a galactic environment. The majority of stars captured by a substructure move on temporarily bound orbits that are lost to galactic tides after a few orbital revolutions. In numerical experiments where a substructure model is immersed into a sea of field particles on a circular orbit, we find a population of particles that remain bound to the substructure potential for indefinitely long times. This population is absent from substructure models, initially placed outside the galaxy on an eccentric orbit. We show that gravitational capture is most efficient in dwarf spheroidal galaxies (dSphs) on account of their low velocity dispersions and high stellar phase-space densities. In these galaxies, ‘dark’ sub-subhaloes, which do not experience in situ star formation, may capture field stars and become visible as stellar overdensities with unusual properties: (i) they would have a large size for their luminosity, (ii) contain stellar populations indistinguishable from the host galaxy, and (iii) exhibit dark matter (DM)-dominated mass-to-light ratios. We discuss the nature of several ‘anomalous’ stellar systems reported as star clusters in the Fornax and Eridanus II dSphs that exhibit some of these characteristics. DM sub-subhaloes with a mass function ${\rm d}N/{\rm d}M_\bullet \sim M_\bullet ^{-\alpha }$ are expected to generate stellar systems with a luminosity function, ${\rm d}N/{\rm d}M_\star \sim M_\star ^{-\beta }$, where $\beta =(2\alpha +1)/3=1.6$ for $\alpha =1.9$. Detecting and characterizing these objects in dSphs would provide unprecedented constraints on the particle mass and cross-section of a large range of DM particle candidates.

     
    more » « less
  5. ABSTRACT

    Small-scale dark matter structures lighter than a billion solar masses are an important probe of primordial density fluctuations and dark matter microphysics. Due to their lack of starlight emission, their only guaranteed signatures are gravitational in nature. We report on results of a search for astrometric weak lensing by compact dark matter subhaloes in the Milky Way with Gaia DR3 data. Using a matched-filter analysis to look for correlated imprints of time-domain lensing on the proper motions of background stars in the Magellanic Clouds, we exclude order-unity substructure fractions in haloes with masses Ml between 107 and $10^9 \, {\rm M}_\odot$ and sizes of one parsec or smaller. We forecast that a similar approach based on proper accelerations across the entire sky with data from Gaia DR4 may be sensitive to substructure fractions of fl ≳ 10−3 in the much lower mass range of $10 \, {\rm M}_\odot \lesssim M_l \lesssim 3 \times 10^3 \, {\rm M}_\odot$. We further propose an analogous technique for stacked star–star lensing events in the regime of large impact parameters. Our first implementation is not yet sufficiently sensitive but serves as a useful diagnostic and calibration tool; future data releases should enable average stellar mass measurements using this stacking method.1

     
    more » « less