Cold Dark Matter Simulations Research Articles

Fast simulation mapping: From standard to modified gravity cosmologies using the bias assignment method

Context. We assess the effectiveness of a non-parametric bias model in generating mock halo catalogues for modified gravity (MG) cosmologies, relying on the distribution of dark matter from either MG or Λ cold dark matter (ΛCDM) simulations. Aims. We aim to generate halo catalogues that effectively capture the distinct impact of MG, ensuring high accuracy in both two- and three-point statistics for a comprehensive analysis of large-scale structures. We investigated the inclusion of MG in non-local bias to directly map the tracers onto ΛCDM fields, which would significantly reduce computational costs. Methods. We employed the bias assignment method (BAM) to model halo distribution statistics by leveraging seven high-resolution COLA simulations of MG cosmologies. Taking cosmic-web dependences into account when learning the bias relations, we designed two experiments to map the MG effects: one utilising the consistent MG density fields and the other employing the benchmark ΛCDM density field. Results. BAM generates MG halo catalogues from both calibration experiments with excellent summary statistics, achieving a ~1% accuracy in the power spectrum across a wide range of k modes, with minimal differences well below 10% for modes subject to cosmic variance, particularly below k < 0.07 h Mpc−1. The reduced bispectrum remains consistent with the reference catalogues within 10% for the studied configuration. Our results demonstrate that a non-linear and non-local bias description can model the effects of MG starting from a ΛCDM field.

Nearest neighbour analysis as a new probe for fuzzy dark matter

ABSTRACT Fuzzy dark matter (FDM) is a promising candidate for dark matter (DM), characterized by its ultra-light mass, which gives rise to wave effects at astrophysical scales. These effects offer potential solutions to the small-scale issues encountered within the standard cold dark matter (CDM) paradigm. In this paper, we investigate the large-scale structure of the cosmic web using FDM simulations, comparing them to CDM-only simulations and a simulation incorporating baryonic effects. Our study employs the nearest neighbour (NN) analysis as a new statistical tool for examining the structure and statistics of the cosmic web in an FDM universe. This analysis could capture the information absent in the two-point correlation functions. In particular, we analyse data related to the spherical contact, nearest neighbour distances (NND), and the angle between the first and second nearest neighbours of haloes (NNA). Specifically, we utilize probability distribution functions, statistical moments, and fitting parameters, as well as G(x), F(x), and J(x) functions to analyse the above data. Remarkably, the results from the FDM simulations differ significantly from the others across these analyses, while no noticeable distinction is observed between the baryonic and CDM-only simulations. Moreover, the lower FDM mass leads to more significant deviations from the CDM simulations. These compelling results highlight the efficiency of the NN analysis – mainly through the use of the J(x) function, $s_3$, $l_{3}$, and $a_4$ parameters – as a prominent new tool for investigating FDM on large scales and making observational predictions.

Satellite group infall into the Milky Way: Exploring the Crater-Leo case with new HST proper motions

Context. Within Λ cold dark matter (ΛCDM) simulations, Milky Way-like galaxies accrete some of their satellite galaxies in groups of 3–5 members rather than individually, and this has been suggested as a possible mechanism driving the formation of satellite planes. Objects accreted in groups are expected to share similar specific total energy and angular momentum, and to also have identical orbital planes and directions. Aims. Looking at observatio ns of Milky Way satellites, the dwarf galaxies Leo II, IV, V, and Crater II, and the star cluster Crater 1 were proposed to be a vestige of group infall. The suggested ‘Crater-Leo group’ shows a coherent distance gradient and all these objects align along a great circle on the sky. We used proper motion data to investigate whether the phase-space distribution of the members of the proposed group are indeed consistent with group infall. Methods. To further investigate this possibility, we used Gaia Data Release 3 (DR3) and new Hubble Space Telescope (HST) proper motions – namely, (μα*, μδ) = (−0.1921 ± 0.0514, −0.0686 ± 0.0523) mas yr−1 for Leo IV and (μα*, μδ) = (0.1186 ± 0.1943, −0.1183 ± 0.1704) mas yr−1 for Leo V – to derive accurate orbital properties for the proposed group objects. In addition, we explored other possible members of this putative association. Results. Leo II, Leo IV, and Crater 1 show orbital properties consistent with those we predict from assuming group infall. However, our results suggest that Crater II was not accreted with the rest of the objects. If confirmed with increasingly accurate proper motions in the future, the Crater-Leo objects would appear to constitute the first identified case of a cosmologically expected, typical group infall event, as opposed to the highly hierarchical Magellanic Cloud system.

Effect of the Large Magellanic Cloud on the kinematics of Milky Way satellites and virial mass estimate

We present a study illustrating the effects of the passage of a Large Magellanic Cloud (LMC) mass satellite on the distance and velocity distributions of satellites in Λ+Cold Dark Matter simulations of Milky Way (MW) sized halos. In agreement with previous studies, we find that during such a passage the velocity distribution develops a high-velocity tail, which can bias velocity-based virial halo mass estimates. When the velocity distribution of MW satellites is corrected for effects of the LMC passage, it is consistent with the distributions in halos of masses as low as and as high as 1.5×1012M⊙. We present a new halo mass estimator , where c is the coefficient calibrated using satellite systems in the simulated MW-sized halos, is the variance of 3D velocities taken with the sign of the radial velocity of each satellite, and is the median halocentric distance of the satellites. We show that the estimator has only s=8% scatter around the median relation of the estimated and true halo masses and deviates by <2s from the median during the pericentric passage of an LMC-like subhalo. This is because and deviate in the opposite directions during such passages. We apply the estimator to the MW satellite system and estimate the virial mass of the Milky Way of , in good agreement with several recent estimates using other methods.

Varying primordial state fractions in exo- and endothermic SIDM simulations of Milky Way-mass haloes

Abstract Self-interacting dark matter (SIDM) is increasingly studied as a potential solution to small-scale discrepancies between simulations of cold dark matter (CDM) and observations. We examine a physically motivated two-state SIDM model with both elastic and inelastic scatterings. In particular, endothermic, exothermic, and elastic scattering have equal transfer cross sections at high relative velocities (vrel ≳ 400 km s−1). In a suite of cosmological zoom-in simulation of Milky Way-size haloes, we vary the primordial state fractions to understand the impact of inelastic dark matter self-interactions on halo structure and evolution. In particular, we test how the initial conditions impact the present-day properties of dark matter haloes. Depending on the primordial state fraction, scattering reactions will be dominated by either exothermic or endothermic effects for high and low initial excited state fractions respectively. We find that increasing the initial excited fraction reduces the mass of the main halo, as well as the number of subhaloes on all mass scales. The main haloes are cored, with lower inner densities and higher outer densities compared with CDM. Additionally, we find that the shape of the main halo becomes more spherical the higher the initial excited state fraction is. Finally, we show that the number of satellites steadily decreases with initial excited state fraction across all satellite masses.

Dissipative Dark Substructure: The Consequences of Atomic Dark Matter on Milky Way Analog Subhalos

Abstract Using cosmological hydrodynamical zoom-in simulations, we explore the properties of subhalos in Milky Way analogs that contain a subcomponent of atomic dark matter (ADM). ADM differs from cold dark matter (CDM) due to the presence of self-interactions that lead to energy dissipation, analogous to standard model baryons. This model can arise in dark sectors that are natural and theoretically motivated extensions to the standard model. The simulations used in this work were carried out using GIZMO and utilize the FIRE-2 galaxy formation physics in the standard model baryonic sector. For the parameter points we consider, the ADM gas cools efficiently, allowing it to collapse to the center of subhalos. This increases a subhalo’s central density and affects its orbit, with more subhalos surviving small pericentric passages. The subset of subhalos that host satellite galaxies have cuspier density profiles and smaller stellar half-mass radii relative to CDM. The entire population of dwarf galaxies produced in the ADM simulations is more compact than those seen in CDM simulations, unable to reproduce the entire diversity of observed dwarf galaxy structures. Additionally, we also identify a population of highly compact subhalos that consist nearly entirely of ADM and form in the central region of the host, where they can leave distinctive imprints in the baryonic disk. This work presents the first detailed exploration of subhalo properties in a strongly dissipative dark matter scenario, providing intuition for how other regions of ADM parameter space, as well as other dark sector models, would impact galactic-scale observables.

The missing quasar image in the gravitationally lensed quasar HE0230−2130: Implications for the cored lens mass distribution and dark satellites

Strongly lensed systems with peculiar configurations allow us to probe the local properties of the deflecting lens mass while simultaneously testing general profile assumptions. The quasar HE0230−2130 is lensed by two galaxies at similar redshifts (Δz ∼ 0.003) into four observed images. Using modeled quasar positions from fitting the brightness of the quasar images in ground-based imaging data from the Magellan telescope, we find that lens-mass models where each of these two galaxies is parametrized with a singular power-law (PL) profile predict five quasar images. One of the predicted images is unobserved despite it being distinctively offset from the lensing galaxies and likely bright enough to be observable. This missing image gives rise to new opportunities to study the mass distribution of these galaxies. To interpret the quad configuration of the system, we tested 12 different profile assumptions with the aim of obtaining lens-mass models that correctly predict only four observed images. We tested the effects of adopting: cored profiles for the lensing galaxies; external shear; and additional profiles to represent a dark matter clump. We find that half of our model classes can produce the correct image multiplicity. By comparing the Bayesian evidence of different model parametrizations, we favor two model classes: (i) one that incorporates two singular PL profiles for the lensing galaxies and a cored isothermal sphere in the region of the previously predicted fifth image (rNIS profile), and (ii) one with a bigger lensing galaxy parametrized by a singular PL profile and the smaller galaxy by a cored PL profile with external shear. We estimated the mass of the rNIS clump for each candidate model of our final Markov chain Monte Carlo sample, and find that only 2% are in the range of 106 M⊙ ≤ MrNIS ≤ 109 M⊙, which is the predicted mass range of dark matter subhalos in cold dark matter simulations, or the mass of dark-matter-dominated and low-surface-brightness galaxies. We therefore favor the models with a cored mass distribution for the lens galaxy close to the predicted fifth image. Our study further demonstrates that lensed quasar images are sensitive to the dark matter structure in the gravitational lens. We are able to describe this exotic lensing configuration with relatively simple models, which demonstrates the power of strong lensing for studying galaxies and lens substructure.

Lopsided Galaxy Distributions in Massive TNG300 Clusters

Studies of isolated central galaxies and their satellites have shown that the spatial distributions of the satellites are distinctly “lopsided” with respect to the locations of the central galaxies. Here we extend this type of analysis to larger systems by analyzing the lopsidedness of 280 massive galaxy clusters in the IllustrisTNG300 Λ Cold Dark Matter simulation. Using a pairwise clustering statistic, we compute the probability that pairs of cluster galaxies are separated by a given polar angle difference, Δϕ, in the plane of the sky. Relative to the location of the central cluster galaxy, we find a statistically significant excess of galaxy pairs that are located on the same side of the cluster. The lopsidedness of the galaxy distribution is most pronounced for large clusters and for pairs of intrinsically red galaxies. The results summarized here were presented at the 243rd meeting of the American Astronomical Society.

The Effective Field Theory of Large-scale Structures of a Fuzzy Dark Matter Universe

Ultralight scalar fields and their noninteracting class, i.e., the so-called fuzzy dark matter (FDM), are dark matter candidates introduced to solve the small-scale problems of the standard cold dark matter. In this paper, we investigate whether the physics of FDM, particularly the quantum pressure that leads to the suppression of structure formation on small scales, could leave significant imprints on the large-scale statistics of matter fluctuations. For this purpose, we utilize the Effective Field Theory of Large Scale Structures (EFT of LSS), wherein small-scale physics is integrated and represented on large scales by only a set of free parameters. These parameters can be determined by fitting them into the cosmological simulations. By fitting the EFT predictions to the simulation data, we determine the value of the speed of sound as a quantitative measure of how UV physics affects large-scale perturbation. We use the Gadget-2 code to study the evolution of 5123 particles in a box with a side length 250 h −1 Mpc. We exploit the suppressed FDM initial power for the FDM universe and perform N-body simulation sufficient to produce accurate—enough for our purpose—results on large scales. In particular, we perform three FDM simulations with different masses and compare their sound speed with the standard cold dark matter (CDM) simulation. We found no difference between the FDM and CDM sound speeds beyond the confidence intervals. However, a consistently increasing trend can be seen in the sound speed for lower masses. This result suggests further investigations using higher-resolution simulations.

Modelling post-reionization H i distributions in fuzzy dark matter cosmologies using conditional normalizing flows

ABSTRACT Upcoming 21-cm intensity mapping experiments like the Square Kilometer Array (SKA) hold significant potential to constrain the properties of dark matter. In this work, we model neutral hydrogen (H i) distributions using high-resolution hydrodynamical N-body simulations of both cold dark matter (CDM) and fuzzy dark matter (FDM) cosmologies in the post-reionization redshift range of z = 3.42−4.94. We show that the H i abundance decreases in FDM-like cosmologies. Extreme FDM models with m ∼ 10−22 eV are at odds with a range of measurements. Due to the increased halo bias, the H i bias increases, paralleled by the damped Lyman-α (DLA) bias which we infer from the cross-section of DLAs. The distribution of the latter in extreme FDM models has a high median at the low-mass end, which can be traced to the high column density of cosmic filaments. FDM models exhibit a very similar abundance of DLAs compared to CDM while sub-DLAs are already less abundant. We study the prospects of detecting the brightest H i peaks with SKA1-Low at z = 4.94, indicating moderate signal-to-noise ratios (SNRs) at angular resolution θA = 2 arcmin with a rapidly declining SNR for lower values of θA. After training the conditional normalizing flow network HIGlow on 2D H i maps, we interpolate its latent space of axion masses to predict the peak flux for a new, synthetic FDM cosmology, finding good agreement with expectations. This work thus underscores the potential of normalizing flows in capturing complex, non-linear structures within H i maps, offering a versatile tool for conditional sample generation and prediction tasks.

Halo assembly in cold and warm dark matter during the JWST frontier epoch

ABSTRACT The JWST mission is in the process of probing the galaxy mass function at z > 10, when conceivably any delay in halo assembly due to the presence of a dwarf galaxy-scale power spectrum cutoff may drastically suppress the number of galaxies relative to the cold dark matter (CDM) expectation. We employ N-body simulations of CDM and warm dark matter (WDM) to explore how the difference in halo collapse time between these models scales with z = 0 descendant halo mass. We demonstrate that collapse begins first for the most massive haloes, and the delay in collapse time between CDM and WDM haloes correlates inversely with descendant mass. We thus infer that only present-day dwarf galaxies exhibit any difference in their assembly history between CDM and WDM at z = 10, and therefore support previous studies that have found JWST is unlikely to determine whether our Universe is better described by the CDM cosmology or the WDM cosmology.

Exploring the correlation between dark matter, intracluster light, and globular cluster distribution in SMACS0723

We present a free-form model of SMACS0723, the first cluster observed with JWST. This model does not make any strong assumptions on the distribution of mass (mostly made up of dark matter) in the cluster and we use it to study the possible correlation between dark matter with the intracluster light and distribution of globular clusters (GCs). To explore the uncertainty in mass modeling, we derived three lens models based on spectroscopically confirmed systems and new candidate systems with redshifts predicted by the lens model derived from the spectroscopic systems. We find that beyond the radius of influence for the brightest cluster galaxy (BCG), the total mass does not trace the intracluster light (ICL), implying the need for a dark component (dark matter). Two loop-like structures observed in the intracluster light do not have any obvious correspondence with the total mass (of mostly dark matter) distribution. The radial profiles of the ICL and the distribution of GCs are similar to each other, but they are steeper than the profile of the lens model. More specifically, we find that the total mass is shallower by 1 dex in log scale than both ICL and GC profiles. This is in excellent agreement with current N-body simulations of cold dark matter.

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.

Densities and mass assembly histories of the Milky Way satellites are not a challenge to ΛCDM

ABSTRACT We use the GRUMPY galaxy formation model based on a suite of zoom-in, high-resolution, dissipationless Λ Cold Dark Matter (ΛCDM) simulations of the Milky Way (MW) sized haloes to examine total matter density within the half-mass radius of stellar distribution, ρtot(< r1/2), of satellite dwarf galaxies around the MW hosts and their mass assembly histories. We compare model results to ρtot(< r1/2) estimates for observed dwarf satellites of the Milky Way spanning their entire luminosity range. We show that observed MW dwarf satellites exhibit a trend of decreasing total matter density within a half-mass radius, ρtot(< r1/2), with increasing stellar mass. This trend is in general agreement with the trend predicted by the model. None of the observed satellites are overly dense compared to the results of our ΛCDM-based model. We also show that although the halo mass of many satellite galaxies is comparable to the halo mass of the MW progenitor at z ≳ 10, at these early epochs halos that survive as satellites to z = 0 are located many virial radii away from the MW progenitors and thus do not have a chance to merge with it. Our results show that neither the densities estimated in observed Milky Way satellites nor their mass assembly histories pose a challenge to the ΛCDM model. In fact, the broad agreement between density trends with the stellar mass of the observed and model galaxies can be considered as yet another success of the model.

Endothermic self-interacting dark matter in Milky Way-like dark matter haloes

ABSTRACT Self-interacting dark matter (SIDM) offers the potential to mitigate some of the discrepancies between simulated cold dark matter (CDM) and observed galactic properties. We introduce a physically motivated SIDM model to understand the effects of self interactions on the properties of Milky Way and dwarf galaxy sized haloes. This model consists of dark matter with a nearly degenerate excited state, which allows for both elastic and inelastic scattering. In particular, the model includes a significant probability for particles to up-scatter from the ground state to the excited state. We simulate a suite of zoom-in Milky Way-sized N-body haloes with six models with different scattering cross sections to study the effects of up-scattering in SIDM models. We find that the up-scattering reaction greatly increases the central densities of the main halo through the loss of kinetic energy. However, the physical model still results in significant coring due to the presence of elastic scattering and down-scattering. These effects are not as apparent in the subhalo population compared to the main halo, but the number of subhaloes is reduced compared to CDM.

Time evolution of the local gravitational parameters and gravitational wave polarizations in a relativistic MOND theory

The recently proposed Skordis-Z\l{}o\'{s}nik theory is the first relativistic MOND theory that can recover the success of the standard $\Lambda$CDM model at matching observations of the cosmic microwave background. This paper aims to revisit the Newtonian and MOND approximations and the gravitational wave analysis of the theory. For the local gravitational parameters, we show that one could obtain both time-varying effective Newtonian gravitational \textit{constant} $G_\textrm{N}$ and time-varying characteristic MOND acceleration scale $a_\textrm{MOND}$, by relaxing the static assumption extensively adopted in the literature. Specially, we successfully demonstrate how to reproduce the redshift dependence of $a_\textrm{MOND}$ observed in the \textit{Magneticum} cold dark matter simulations. For the gravitational waves, we show that there are only two tensor polarizations, and reconfirm that its speed is equal to the speed of light.

The local dark matter distribution in self-interacting dark matter halos

We study the effects of dark matter self-interactions on the local dark matter distribution in selected Milky Way-like galaxies in the EAGLE hydrodynamical simulations. The simulations were run with two different self-interacting dark matter models, a constant and velocity-dependent self-interaction cross-section. We find that the local dark matter velocity distribution of the Milky Way-like halos in the simulations with dark matter self-interactions and baryons are generally similar to those extracted from cold collisionless dark matter simulations with baryons. In both cases, the local dark matter speed distributions agree well with their best fit Maxwellian distributions. Including baryons in the simulations with or without dark matter self-interactions increases the local dark matter density and shifts the dark matter speed distributions to higher speeds. To study the implications for direct detection, we compute the dark matter halo integrals obtained directly from the simulations and compare them to those obtained from the best fit Maxwellian velocity distribution. We find that a Maxwellian distribution provides a good fit to the halo integrals of most halos, without any significant difference between the results of different dark matter self-interaction models.

Real and counterfeit cores: how feedback expands haloes and disrupts tracers of inner gravitational potential in dwarf galaxies

ABSTRACT The tension between the diverging density profiles in Lambda cold dark matter simulations and the constant-density inner regions of observed galaxies is a long-standing challenge known as the ‘core–cusp’ problem. We demonstrate that the SMUGGLE galaxy formation model implemented in the arepo moving mesh code forms constant-density cores in idealized dwarf galaxies of M⋆ ≈ 8 × 107 Msun with initially cuspy dark matter (DM) haloes of M200 ≈ 1010 Msun. Identical initial conditions run with an effective equation of state interstellar medium model preserve cuspiness. Literature on the subject has pointed to the low density threshold for star formation, ρth, in such effective models as an obstacle to baryon-induced core formation. Using a SMUGGLE run with equal ρth, we demonstrate that core formation can proceed at low density thresholds, indicating that ρth is insufficient on its own to determine whether a galaxy develops a core. We reaffirm that the ability to resolve a multiphase interstellar medium at sufficiently high densities is a more reliable indicator of core formation than any individual model parameter. In SMUGGLE, core formation is accompanied by large degrees of non-circular motion, with gas rotational velocity profiles that consistently fall below the circular velocity $v_\text{circ} = \sqrt{GM/R}$ out to ∼2 kpc. Asymmetric drift corrections help recover the average underlying DM potential for some of our less efficient feedback runs, but time-variations in the instantaneous azimuthal gas velocity component are substantial, highlighting the need for careful modelling in the inner regions of dwarfs to infer the true distribution of DM.

A universal profile for stacked filaments from cold dark matter simulations

ABSTRACT We study the stacked filaments connecting group-mass halo pairs, using dark-matter-only N-body simulations. We calculate the dark matter overdensity profile of these stacked filaments at different redshifts as a function of the distance perpendicular to the filament axis. A four-parameter universal functional form, including three comoving scale radii and one amplitude parameter (core density), provides a good fit out to a radius of 20h−1cMpc for stacked filaments over a range of redshifts, lengths, and masses. The scale radii are approximately independent of redshift but increase as power laws with the comoving filament length. Lastly, we compare the scaling of the filament mass measured directly from the simulations to the predicted scaling from the halo–halo–matter three-point correlation function as a function of redshift and of the mass of the halo pairs. We find that both measured scalings are similar to, but somewhat shallower than the predictions, by 10 per cent and 30 per cent, respectively. These results provide a template to interpret present and upcoming observational results based on stacking, for example, weak lensing, thermal and kinetic Sunyaev-Zel’dovich, or X-ray observations.

X-ray morphology of cluster-mass haloes in self-interacting dark matter

ABSTRACT We perform cosmological zoom-in simulations of 19 relaxed cluster-mass haloes with the inclusion of adiabatic gas in the cold dark matter (CDM) and self-interacting dark matter (SIDM) models. These clusters are selected as dynamically relaxed clusters from a parent simulation with $M_{\rm 200} \simeq (1\!-\!3)\times 10^{15}{\, \rm M_\odot }$. Both the dark matter and the intracluster gas distributions in SIDM appear more spherical than their CDM counterparts. Mock X-ray images are generated based on the simulations and are compared to the real X-ray images of 84 relaxed clusters selected from the Chandra and ROSAT archives. We perform ellipse fitting for the isophotes of mock and real X-ray images and obtain the ellipticities at cluster-centric radii of $r\simeq 0.1\!-\!0.2R_{\rm 200}$. The X-ray isophotes in SIDM models with increasing cross-sections are rounder than their CDM counterparts, which manifests as a systematic shift in the distribution function of ellipticities. Unexpectedly, the X-ray morphology of the observed non-cool-core clusters agrees better with SIDM models with cross-section $(\sigma /m)= 0.5\!-\!1\, {\rm cm}^2\, {\rm g}^{-1}$ than CDM and SIDM with $(\sigma /m)=0.1\, {\rm cm}^2\, {\rm g}^{-1}$. Our statistical analysis indicates that the latter two models are disfavoured at the $68{{\ \rm per\ cent}}$ confidence level (as conservative estimates). This conclusion is not altered by shifting the radial range of measurements or applying a temperature selection criterion. However, the primary uncertainty originates from the lack of baryonic physics in the adiabatic model, such as cooling, star formation and feedback effects, which still have the potential to reconcile CDM simulations with observations.