Low-mass Subhaloes Research Articles

Survival of gas in subhalos and its impact on the 21 cm forest signals: insights from hydrodynamic simulations

Understanding the survival of gas within subhalos under various astrophysical processes is crucial for elucidating cosmic structure formation and evolution. We study the resilience of gas in subhalos, focusing on the impact of tidal and ram pressure stripping through hydrodynamic simulations. Our results uncover significant gas stripping primarily driven by ram pressure effects, which also profoundly influence the gas distribution within these subhalos. Notably, despite their vulnerability to ram pressure effects, the low-mass subhalos can play a pivotal role in influencing the observable characteristics of cosmic structures due to their large abundance. Specifically, we explore the application of our findings to the 21 cm forest, showing how the survival dynamics of gas in subhalos can modulate the 21 cm optical depth, a key probe for detecting minihalos in the pre-reionization era. Our previous study demonstrated that the 21 cm optical depth can be enhanced by the subhalos, but the effects of tidal and ram pressure stripping on the subhalo abundance have not been fully considered. In this work, we further investigate the contribution of subhalos to the 21 cm optical depth with hydrodynamic simulations, particularly highlighting the trajectories and fates of subhalos within mass ranges of 104-6 M ⊙ h -1 in a host halo of 107 M ⊙h-1, and subhalos within mass range of 104-5 M ⊙h-1 in a host halo of 106 M ⊙h-1. Despite their susceptibility to ram pressure stripping, the contribution of abundant low-mass subhalos to the 21 cm optical depth is more significant than that of their massive counterparts primarily due to their greater abundance. We find that the 21 cm optical depth can be increased by a factor of approximately two due to the abundant low-mass subhalos. However, this enhancement is about twice as low as previously estimated in our earlier study, a discrepancy attributed to the effects of ram pressure stripping. Our work provides critical insights into the gas dynamics within subhalos in the early universe, highlighting their resilience against environmental stripping effects, and their impact on observable 21 cm signals.

THE THREE HUNDRED: Contrasting cluster galaxy density in hydrodynamical and dark matter simulations

Cluster number counts at visible and IR wavelengths will be a key cosmological probe in the next decade thanks to the Euclid satellite mission. For this purpose, the performance of cluster detection algorithms, which at these wavelengths are sensitive to the spatial distributions of the cluster galaxy members and their luminosity functions, need to be accurately characterized. Using THE THREE HUNDRED hydrodynamical and dark-matter-only simulations, we studied a complete sample of massive clusters beyond 7 (5) × 1014 M⊙ at redshift 0 (1) on a (1.48 Gpc)3 volume. We find that the mass resolution of the current hydrodynamical simulations (1.5 × 109 M⊙) is not enough to characterize the luminosity function of the sample in the perspective of Euclid data. Nevertheless, these simulations are still useful to characterize the spatial distribution of the cluster substructures assuming a common relative mass threshold for the different flavours and resolutions. By comparing with the dark-matter-only version of these simulations, we demonstrate that baryonic physics preserves significantly low-mass subhalos (galaxies), as has also been observed in previous studies with less statistics. Furthermore, by comparing the hydro simulations with higher resolution dark-matter-only simulations of the same objects and taking the same limit in subhalo mass, we find galaxy density profiles that are significantly more cuspy towards the centre of the clusters, where the low-mass substructures tend to concentrate. We conclude that using a dark-matter-only simulation may lead to some biases on the spatial distribution and density of galaxy cluster members. Based on the preliminary analysis of few high-resolution hydro simulations we conclude that a mass resolution of 1.8 × 108h−1 M⊙ will be needed for THE THREE HUNDRED simulations to approach the expected magnitude limits for the Euclid survey. These simulations are currently under way.

The viability of low-mass subhaloes as targets for gamma-ray dark matter searches

ABSTRACT In this work, we investigate the discovery potential of low-mass Galactic dark matter (DM) subhaloes for indirect searches of DM. We use data from the Via Lactea II (VL-II) N-body cosmological simulation, which resolves subhaloes down to $\mathcal {O}(10^4)$ solar masses and it is thus ideal for this purpose. First, we characterize the abundance, distribution and structural properties of the VL-II subhalo population in terms of both subhalo masses and maximum circular velocities. Then, we repopulate the original simulation with millions of subhaloes of masses down to about five orders of magnitude below the minimum VL-II subhalo mass (more than one order of magnitude in velocities). We compute subhalo DM annihilation astrophysical ‘J-factors’ and angular sizes for the entire subhalo population, by placing the Earth at a random position but at the right Galactocentric distance in the simulation. Thousands of these realizations are generated in order to obtain statistically meaningful results. We find that some nearby low-mass Galactic subhaloes, not massive enough to retain stars or gas, may indeed yield DM annihilation fluxes comparable to those expected from other, more massive and acknowledgeable DM targets like dwarf satellite galaxies. Typical angular sizes are of the order of the degree, thus subhaloes potentially appearing as extended sources in gamma-ray telescopes, depending on instrument angular resolution and sensitivity. Our work shows that low-mass Galactic subhaloes with no visible counterparts are expected to play a relevant role in current and future indirect DM search searches and should indeed be considered as excellent DM targets.

Prospects for Detecting Gaps in Globular Cluster Stellar Streams in External Galaxies with the Nancy Grace Roman Space Telescope

Stellar streams form through the tidal disruption of satellite galaxies or globular clusters orbiting a host galaxy. Globular cluster streams are exciting since they are thin (dynamically cold) and therefore sensitive to perturbations from low-mass subhalos. Since the subhalo mass function differs depending on the dark matter composition, these gaps can provide unique constraints on dark matter models. However, current samples are limited to the Milky Way. With its large field of view, deep imaging sensitivity, and high angular resolution, the upcoming Nancy Grace Roman Space Telescope (Roman) presents a unique opportunity to increase the number of observed streams and gaps significantly. This paper presents a first exploration of the prospects for detecting gaps in streams in M31 and other nearby galaxies with resolved stars. We simulate the formation of gaps in a Palomar 5–like stream and generate mock observations of these gaps with background stars in M31 and foreground Milky Way stellar fields. We assess Roman's ability to detect gaps out to 10 Mpc through visual inspection and with the gap-finding tool FindTheGap. We conclude that gaps of ≈1.5 kpc in streams that are created from subhalos of masses ≥5 × 106 M ⊙ are detectable within a 2–3 Mpc volume in exposure times of 1000 s to 1 hr. This volume contains ≈150 galaxies, including ≈eight galaxies with luminosities >109 L ⊙. Large samples of stream gaps in external galaxies will open up a new era of statistical analyses of gap characteristics in stellar streams and help constrain dark matter models.

Extending the unified subhalo model to warm dark matter

ABSTRACT Using a set of high-resolution N-body simulations, we extend the unified distribution model of cold dark matter (CDM) subhaloes to the warm dark matter (WDM) case. The same model framework combining the unevolved mass function, unevolved radial distribution, and tidal stripping can predict the mass function and spatial distribution of subhaloes in both CDM and WDM simulations. The dependence of the model on the DM particle property is universally parametrized through the half-mode mass of the initial power spectrum. Compared with the CDM model, the WDM model differs most notably in two aspects. (1) In contrast to the power-law form in CDM, the unevolved subhalo mass function for WDM is scale-dependent at the low mass end due to the cut-off in the initial power spectrum. (2) WDM subhaloes are more vulnerable to tidal stripping and disruption due to their lower concentrations at accretion time. Their survival rate is also found to depend on the infall mass. Accounting for these differences, the model predicts a final WDM subhalo mass function that is also proportional to the unevolved subhalo mass function. The radial distribution of WDM subhaloes is predicted to be mass-dependent. For low mass subhaloes, the radial distribution is flatter in the inner halo and steeper in the outer halo compared to the CDM counterpart, due to the scale-dependent unevolved mass function and the enhanced tidal stripping. The code for sampling subhaloes according to our generalized model is available at https://github.com/fhtouma/subgen2.

Soft scattering evaporation of dark matter subhalos by inner galactic gases

The large gap between a galactic dark matter subhalo’s velocity and its own gravitational binding velocity creates the situation that small subhalos can be evaporated before dark matter thermalize with baryons due to the low binding velocity. In case dark matter acquires an electromagnetic dipole moment, the survival of low-mass subhalos requires stringent limits on the photon-mediated soft scattering. The current stringent direct detection limits indicate for a small dipole moment, which lets DM decouple early and allows small subhalos to form. We calculate the DM kinetic decoupling temperature in the Early Universe and evaluate the smallest protohalo mass. In the late Universe, low-mass subhalos can be evaporated via soft collision by ionized gas and accelerated cosmic rays. We calculate the subhalos evaporation rate and show that subhalos lighter than 10^{-5}M_{odot } in the gaseous inner galactic region are subject to evaporation via dark matter’s effective electric and magnetic dipole moments below current direct detection limits, which potentially affects the low-mass subhalos distribution in the galactic center.

Matching the mass function of Milky Way satellites in competing dark matter models

ABSTRACT Any successful model of dark matter must explain the diversity of observed Milky Way (MW) satellite density profiles, from very dense ultrafaints to low-density satellites so large that they could be larger than their inferred dark matter haloes. Predictions for these density profiles are complicated by the limitations of simulation resolution in the stripping of subhaloes by the MW system. We consider cold dark matter (CDM), warm dark matter (WDM, 3.3 keV thermal relic power spectrum), and a self-interacting dark matter model (SIDM) that induces gravothermal collapse in low-mass subhaloes. Using N-body simulations combined with a halo stripping algorithm, we find that most CDM and WDM subhaloes of mass >108 ${\, \rm M_\odot }$ are large enough after stripping to fit most satellites; however, the required amount of stripping often requires a stronger tidal field than is available on the subhalo’s orbit. The lower concentrations of WDM subhaloes enable more stripping to take place, even on orbits with large pericentres. SIDM cores offer the best fits to massive, low-density satellites at the expense of predicting >109 ${\, \rm M_\odot }$ subhaloes to host low-density satellites with no observed analogue. The agreement of the total number of satellites with observations in CDM and WDM depends strongly on the assumptions made to draw the observational estimates. We conclude that an SIDM model must have a very high velocity-dependent cross-section in order to match all satellites, and that WDM offers a marginally better fit than CDM to the MW satellite mass function.

Shedding light on low-mass subhalo survival and annihilation luminosity with numerical simulations

ABSTRACT In this work, we carry out a suite of specially designed numerical simulations to shed light on dark matter (DM) subhalo survival at mass scales relevant for gamma-ray DM searches, a topic subject to intense debate nowadays. We have employed an improved version of DASH, a GPU N-body code, to study the evolution of low-mass subhaloes inside a Milky-Way-like halo with unprecedented accuracy, reaching solar-mass and sub-parsec resolution. We simulate subhaloes with varying mass, concentration, and orbital properties, and consider the effect of baryons in the host. We analyse the evolution of the bound mass fraction and annihilation luminosity, finding that most subhaloes survive until present, yet losing in some cases more than 99 per cent of their initial mass. Baryons induce a much greater mass-loss, especially when the subhalo orbit is more parallel to the Galactic disc. Many of these subhaloes cross the solar Galactocentric radius, making it easier to detect their annihilation fluxes from Earth. We find subhaloes orbiting a DM-only halo with a pericentre in the solar vicinity to lose 70–90 per cent of their initial annihilation luminosity at present, which increases up to 99 per cent when including baryons. We find a strong relation between subhalo’s mass-loss and the effective tidal field at pericentre. Indeed, much of the dependence on all considered parameters can be explained through this single parameter. In addition to shedding light on the survival of low-mass Galactic subhaloes, our results can provide detailed predictions that will aid current and future quests for the nature of DM.

Milky Way Satellite Census. IV. Constraints on Decaying Dark Matter from Observations of Milky Way Satellite Galaxies

We use a recent census of the Milky Way (MW) satellite galaxy population to constrain the lifetime of particle dark matter (DM). We consider two-body decaying dark matter (DDM) in which a heavy DM particle decays with lifetime τ comparable to the age of the universe to a lighter DM particle (with mass splitting ϵ) and to a dark radiation species. These decays impart a characteristic “kick velocity,” V kick = ϵ c, on the DM daughter particles, significantly depleting the DM content of low-mass subhalos and making them more susceptible to tidal disruption. We fit the suppression of the present-day DDM subhalo mass function (SHMF) as a function of τ and V kick using a suite of high-resolution zoom-in simulations of MW-mass halos, and we validate this model on new DDM simulations of systems specifically chosen to resemble the MW. We implement our DDM SHMF predictions in a forward model that incorporates inhomogeneities in the spatial distribution and detectability of MW satellites and uncertainties in the mapping between galaxies and DM halos, the properties of the MW system, and the disruption of subhalos by the MW disk using an empirical model for the galaxy–halo connection. By comparing to the observed MW satellite population, we conservatively exclude DDM models with τ < 18 Gyr (29 Gyr) for V kick = 20 kms−1 (40 kms−1) at 95% confidence. These constraints are among the most stringent and robust small-scale structure limits on the DM particle lifetime and strongly disfavor DDM models that have been proposed to alleviate the Hubble and S 8 tensions.

An accurate comprehensive approach to substructure – II. Stripped subhaloes

ABSTRACT In Salvador-Solé, Manrique &amp; Botella (Paper I), we used the ConflUent System of Peak trajectories (CUSP) formalism to derive from first principles and no single free parameter the accurate abundance and radial distribution of both diffuse dark matter (dDM) and subhaloes accreted on to haloes and their progenitors at all previous times. Here we use those results as initial conditions for the monitoring of the evolution of subhaloes and dDM within the host haloes. Specifically, neglecting dynamical friction, we accurately calculate the effects of repetitive tidal stripping and heating on subhaloes as they orbit inside the host halo and infer the amount of dDM and subsubhaloes they release into the intrahalo medium. We then calculate the expected abundance and radial distribution of stripped subhaloes and dDM. This derivation clarifies the role of halo concentration in substructure and unravels the origin of some key features found in simulations including the dependence of substructure on halo mass. In addition, it unveils the specific effects of dynamical friction on substructure. The results derived here are for purely accreting haloes. In Salvador-Solé et al. (Paper III), we complete the study by addressing the case of low-mass subhaloes, unaffected by dynamical friction, in ordinary haloes having suffered major mergers.

Sensitivity of the Cherenkov Telescope Array to dark subhalos

In this work, we study the potential of the Cherenkov Telescope Array (CTA) for the detection of Galactic dark matter (DM) subhalos. We focus on low-mass subhalos that do not host any baryonic content and therefore lack any multiwavelength counterpart. If the DM is made of weakly interacting massive particles (WIMPs), these dark subhalos may thus appear in the gamma-ray sky as unidentified sources. A detailed characterization of the instrumental response of CTA to dark subhalos is performed, for which we use the {\it ctools} analysis software and simulate CTA observations under different array configurations and pointing strategies, such as the scheduled extragalactic survey. This, together with information on the subhalo population as inferred from N-body cosmological simulations, allows us to predict the CTA detectability of dark subhalos, i.e., the expected number of subhalos in each of the considered observational scenarios. In the absence of detection, for each observation strategy we set competitive limits to the annihilation cross section as a function of the DM particle mass, that are at the level of $\langle\sigma v\rangle\sim4\times10^{-24}$ ($7\times10^{-25}$) $\mathrm{cm^3s^{-1}}$ for the $b\bar{b}$ ($\tau^+\tau^-$) annihilation channel in the best case scenario. Interestingly, we find the latter to be reached with no dedicated observations, as we obtain the best limits by just accumulating exposure time from all scheduled CTA programs and pointings over the first 10 years of operation. This way CTA will offer the most constraining limits from subhalo searches in the intermediate range between $\sim 1-3$ TeV, complementing previous results with \textit{Fermi}-LAT and HAWC at lower and higher energies, respectively.

Harnessing the Population Statistics of Subhalos to Search for Annihilating Dark Matter

Abstract The Milky Way’s dark matter halo is expected to host numerous low-mass subhalos with no detectable associated stellar component. Such subhalos are invisible unless their dark matter annihilates to visible states such as photons. One of the established methods for identifying candidate subhalos is to search for individual unassociated gamma-ray sources with properties consistent with the dark matter expectation. However, robustly ruling out an astrophysical origin for any such candidate is challenging. In this work, we present a complementary approach that harnesses information about the entire population of subhalos—such as their spatial and mass distribution in the Galaxy—to search for a signal of annihilating dark matter. Using simulated data, we show that the collective emission from subhalos can imprint itself in a unique way on the statistics of observed photons, even when individual subhalos may be too dim to be resolved on their own. Additionally, we demonstrate that, for the models we consider, the signal can be identified even in the face of unresolved astrophysical point-source emission of extragalactic and Galactic origin. This establishes a new search technique for subhalos that is complementary to established methods, and that could have important ramifications for gamma-ray dark matter searches using observatories such as the Fermi Large Area Telescope and the Cerenkov Telescope Array.

The effect of hydrodynamics alone on the subhalo population in a ΛCDM rich cluster sized dark matter halo

We perform a set of non-radiative hydro-dynamical (NHD) simulations of a rich cluster sized dark matter halo from the Phoenix Project with three different numerical resolutions, to investigate the effect of hydrodynamics alone on the subhalo population in the halo. Compared to dark matter only (DMO) simulations of the same halo, subhaloes are less abundant for relatively massive subhaloes (M sub > 2.5 × 109 h −1 M ⊙, or V max > 70 km s−1) but more abundant for less massive subhaloes in the NHD simulations. This results in different shapes in the subhalo mass/V max function in two different sets of simulations. At given subhalo mass, the subhaloes less massive than 1010 h −1 M ⊙ have larger V max in the NHD than DMO simulations, while V max is similar for the subhaloes more massive than the mass value. This is mainly because the progenitors of present day low mass subhaloes have larger concentration parameters in the NHD than DMO simulations. The survival number fraction of the accreted low mass progenitors of the main halo at redshift 2 is about 50 percent higher in the NHD than DMO simulations.

Dynamical Simulations of the First Globular Clusters

Abstract A Milky Way–like halo is simulated with tidally limited star clusters inserted in the dark matter halos present at high redshift. An n-body code augmented with velocity relaxation in the star clusters evolves the system. The stripped stars and remnant clusters that began in the lower mass subhalos have a distribution somewhat more extended than the dark matter halo, with a mean galactic radius of about 60 kpc inside 150 kpc, whereas the distribution of all stripped stars and clusters is more centrally concentrated than the dark matter. The clusters from low mass subhalos, those with a peak circular velocity of 12–18 , also produce most of the population’s thin stellar streams. Together these suggest a search strategy for extremely metal-poor cluster stars and their remnants. The dependence of the stellar population distribution on subhalo mass is not seen in simulations that start clusters at lower redshift. The half mass radii of the clusters are set by the tidal fields of the initial subhalo, “dwarf galaxy” location, which causes the average half mass–radius to decrease with increasing redshift of formation. Starting clusters at greater than redshift 4 leads to cluster half mass radii approximately as seen in the Milky Way.

The Density Structure of Simulated Stellar Streams

Abstract Star particles in a set of dense clusters are self-consistently evolved within an LCDM dark matter distribution with an n-body code. The clusters are started on nearly circular orbits in the more massive sub-halos. Each cluster develops a stellar tidal stream, initially within its original sub-halo. When a sub-halo merges into the main halo the early time stream is dispersed as a somewhat chaotic thick stream, roughly the width of the orbit of the cluster in the sub-halo. Once the cluster orbits freely in the main halo the star stream forms a thin stream again, usually resulting in a thin stream surrounded by a wider distribution of star particles lost at earlier times. To examine the role of the lower-mass dark matter sub-halos in the creation of density variations along the thin tidal star streams two realizations of the simulation are run with and without a normal cold dark matter sub-halo population below 4 × 108 . About 70(40)% of thin streams show density variations that are 2(5) times the star count noise level, irrespective of the presence or absence of low-mass sub-halos. A counts-in-cells analysis (related to the two-point correlation function and power spectrum) of the density along nearly 8000° of streams in the two well matched models finds that the full sub-halo population leads to slightly larger, but statistically significant, density fluctuations on scales of 2°–6°.

Searching for dark matter sub-structure with HAWC

Numerical simulations show that the dark matter halos surrounding galaxies are expected to contain many over-densities or sub-halos. The most massive of these sub-halos can be optically observed in the form of dwarf galaxies. However, most lower mass sub-halos are predicted to exist as dark dwarf galaxies: sub-halos like dwarf galaxies with no luminous counterpart. It may be possible to detect these unseen sub-halos from gamma-ray signals originating from dark matter annihilation. The High Altitude Water Cherenkov Observatory (HAWC) is a very high energy (500 GeV to >100 TeV) gamma ray detector with a wide field-of-view and near continuous duty cycle, making HAWC ideal for unbiased sky surveys. We perform a search for gamma ray signals from dark dwarfs in the Milky Way halo with HAWC. We perform a targeted search of HAWC gamma-ray sources which have no known association with lower-energy counterparts, based on an unbiased survey of the entire sky. With no sources found to strongly prefer dark matter models, we calculate the ability of HAWC to observe dark dwarfs. We also compute the HAWC sensitivity to potential future detections for a given model of dark matter substructure. Assuming thermal dark matter, we find the corresponding J-factor of a dark dwarf required to reach the HAWC detection criterion is 5.79× 1020 GeV2 cm−5 sr for one particular set of dark matter assumptions. HAWC is found to be able to competitively constrain dark matter annihilation from discovered halos with J-factors on the scale of 1019 GeV2 cm−5 sr or greater, with better constraints obtained on dark matter models with >10 TeV masses and sources that transit overhead.

Unidentified gamma-ray sources as targets for indirect dark matter detection with the Fermi-Large Area Telescope

One of the predictions of the ΛCDM cosmological framework is the hierarchical formation of structure, giving rise to dark matter (DM) halos and subhalos. When the latter are massive enough they retain gas (i.e., baryons) and become visible. This is the case of the dwarf satellite galaxies in the Milky Way (MW) . Below a certain mass, halos may not accumulate significant amounts of baryons and remain completely dark. However, if DM particles are Weakly Interacting Massive Particles (WIMPs), we expect them to annihilate in subhalos, producing gamma rays which can be detected with the Fermi satellite. Using the three most recent point-source Fermi Large Area Telescope (LAT) catalogs (3FGL, 2FHL and 3FHL), we search for DM subhalo candidates among the unidentified sources, i.e., sources with no firm association to a known astrophysical object. We apply several selection criteria based on the expected properties of the DM-induced emission from subhalos, which allow us to significantly reduce the list of potential candidates. Then, by characterizing the minimum detection flux of the instrument and comparing our sample to predictions from the Via Lactea II (VL-II) N-body cosmological simulation, we place conservative and robust constraints on the ⟨σ v⟩−mDM parameter space. For annihilation via the τ+τ− channel, we put an upper limit of 4 × 10−26 (5 × 10−25) cm3s−1 for a mass of 10 (100) GeV . A critical improvement over previous treatments is the repopulation we made to include low-mass subhalos below the VL\ obreakdash-II mass resolution. With more advanced subhalo candidate filtering the sensitivity reach of our method can potentially improve these constraints by a factor 3 (2) for τ+τ− (bb̄) channel.

The suppression of star formation on the smallest scales: what role does environment play?

The predominantly ancient stellar populations observed in the lowest-mass galaxies (i.e. ultra-faint dwarfs) suggest that their star formation was suppressed by reionization. Most of the well-studied ultra-faint dwarfs, however, are within the central half of the Milky Way dark matter halo, such that they are consistent with a population that was accreted at early times and thus potentially quenched via environmental processes. To study the potential role of environment in suppressing star formation on the smallest scales, we utilize the Exploring the Local Volume in Simulations (ELVIS) suite of $N$-body simulations to constrain the distribution of infall times for low-mass subhalos likely to host the ultra-faint population. For the ultra-faint satellites of the Milky Way with star-formation histories inferred from $Hubble~Space~Telescope$ imaging, we find that environment is highly unlikely to play a dominant role in quenching their star formation. Even when including the potential effects of pre-processing, there is a $\lesssim 0.1\%$ probability that environmental processes quenched all of the known ultra-faint dwarfs early enough to explain their observed star-formation histories. Instead, we argue for a mass floor in the effectiveness of satellite quenching at roughly $M_{\star} \sim 10^{5}~M_{\odot}$, below which star formation in surviving galaxies is globally suppressed by reionization. We predict a large population of quenched ultra-faint dwarfs in the Local Field ($1 < R/R_{\rm{vir}} < 2$), with as many as $\sim250$ to be discovered by future wide-field imaging surveys.

Power spectrum of dark matter substructure in strong gravitational lenses

Studying the smallest self-bound dark matter structure in our Universe can yield important clues about the fundamental particle nature of dark matter. Galaxy-scale strong gravitational lensing provides a unique way to detect and characterize dark matter substructures at cosmological distances from the Milky Way. Within the cold dark matter (CDM) paradigm, the number of low-mass subhalos within lens galaxies is expected to be large, implying that their contribution to the lensing convergence field is approximately Gaussian and could thus be described by their power spectrum. We develop here a general formalism to compute from first principles the substructure convergence power spectrum for different populations of dark matter subhalos. As an example, we apply our framework to two distinct subhalo populations: a truncated Navarro-Frenk-White subhalo population motivated by standard CDM, and a truncated cored subhalo population motivated by self-interacting dark matter (SIDM). We study in detail how the subhalo abundance, mass function, internal density profile, and concentration affect the amplitude and shape of the substructure power spectrum. We determine that the power spectrum is mostly sensitive to a specific combination of the subhalo abundance and moments of the mass function, as well as to the average tidal truncation scale of the largest subhalos included in the analysis. Interestingly, we show that the asymptotic slope of the substructure power spectrum at large wave number reflects the internal density profile of the subhalos. In particular, the SIDM power spectrum exhibits a characteristic steepening at large wave number absent in the CDM power spectrum, opening the possibility of using this observable, if at all measurable, to discern between these two scenarios.

Subhalo demographics in the Illustris simulation: effects of baryons and halo-to-halo variation

We study the abundance of subhaloes in the hydrodynamical cosmological simulation Illustris, which includes both baryons and dark matter in a $\Lambda$CDM volume 106.5 Mpc a side. We compare Illustris to its dark matter-only (DMO) analogue, Illustris-Dark, and quantify the effects of baryonic processes on the demographics of subhaloes in the host mass range $10^{11}$ to $3 \times 10^{14} \msun$. We focus on both the evolved ($z=0$) subhalo cumulative mass functions (SHMF) and the statistics of subhaloes ever accreted, i.e. infall subhalo mass function. We quantify the variance in subhalo abundance at fixed host mass and investigate the physical reasons responsible for such scatter. We find that in Illustris, baryonic physics impacts both the infall and $z=0$ subhalo abundance by tilting the DMO function and suppressing the abundance of low-mass subhaloes. The breaking of self-similarity in the subhalo abundance at $z=0$ is enhanced by the inclusion of baryonic physics. The non-monotonic alteration of the evolved subhalo abundances can be explained by the modification of the concentration--mass relation of Illustris hosts compared to Illustris-Dark.Interestingly, the baryonic implementation in Illustris does not lead to an increase in the halo-to-halo variation compared to Illustris-Dark. In both cases, the normalized intrinsic scatter today is larger for Milky Way-like haloes than for cluster-sized objects. For Milky Way-like haloes, it increases from about eight per cent at infall to about 25 per cent at the current epoch. In both runs, haloes of fixed mass formed later host more subhaloes than early formers.