The AIDA-TNG project. Abundance, radial distribution, and clustering properties of halos in alternative dark matter models

The AIDA-TNG project utilizes cosmological simulations to demonstrate that alternative dark matter models, specifically warm and self-interacting variants, significantly alter halo abundance, radial distribution, and clustering properties, with small-scale clustering emerging as a key discriminator for future observational comparisons.

The Gist

Imagine the universe as a giant, invisible ocean. In this ocean, there are islands of matter that we can see (stars, galaxies) and a vast, invisible "fog" that holds everything together. This invisible fog is Dark Matter.

For a long time, scientists thought this fog was made of "Cold" particles—like heavy, slow-moving snowflakes that just clump together under gravity. This is the standard model, called $\Lambda$CDM. It works great for big things, but when we look at the small islands (galaxies), the standard model predicts there should be way more tiny satellite galaxies than we actually see. It's like a recipe that says you should have 1,000 chocolate chips in your cookie, but you only count 50.

This paper is like a team of cosmic chefs (the AIDA-TNG project) running a massive simulation kitchen to test two new recipes for that invisible fog: Warm Dark Matter (WDM) and Self-Interacting Dark Matter (SIDM). They want to see if these new ingredients fix the "missing chocolate chip" problem.

Here is what they found, explained simply:

1. The Kitchen Setup (The Simulations)

The scientists didn't just guess; they built a virtual universe inside a computer. They created different "boxes" (simulations) of space, ranging from small neighborhoods to large cities. Inside these boxes, they ran the universe forward in time, but with different rules for how the dark matter fog behaves:

• The Standard (CDM): Heavy, non-interacting snowflakes.
• The Warm Version (WDM): Lighter, faster-moving particles (like warm air). Because they move fast, they can't clump together easily in small groups.
• The Social Version (SIDM): Particles that bump into each other (like people at a crowded party). When they collide, they swap energy and spread out, smoothing out the dense centers of galaxies.

2. Counting the Guests (Abundance)

The first thing they checked was: How many "satellite" galaxies (sub-halos) are there?

• The Warm Recipe (WDM): Because these particles are fast and "warm," they wash away the tiny clumps. It's like trying to build a sandcastle in a strong wind; the small towers get blown away. The result? Fewer small galaxies. This actually helps fix the "missing satellites" problem because the simulation now predicts fewer small islands, closer to what we see in real life.
• The Social Recipe (SIDM): The number of galaxies stays about the same as the standard model. The particles don't disappear; they just move around differently.

3. Where Do They Live? (Radial Distribution)

Next, they looked at where these galaxies sit inside their larger host galaxies. Imagine a host galaxy as a big city, and the satellites as suburbs.

• Standard & Warm Models: The satellites tend to crowd right into the city center. It's like a dense, crowded downtown. The paper found that the standard "crowded downtown" model (NFW profile) wasn't quite right; the satellites were even more concentrated than expected in the Warm models.
• The Social Model (SIDM): Because the dark matter particles bump into each other, they act like a fluid that smooths out the center. Instead of a crowded downtown, you get a "suburban sprawl." The satellites are pushed out to the edges, creating a shallower density profile. The center is less crowded, and the suburbs are more spread out.

4. How They Clump Together (Clustering)

Finally, they looked at how these galaxies group together on a cosmic scale. Do they hang out in tight packs, or are they spread out?

• Warm Dark Matter: Because the small clumps are missing, the remaining galaxies are forced to live in bigger, more massive groups. This makes the "big groups" clump together more tightly. It's like if you remove all the small families from a town, leaving only the large estates; the town looks very different and more clustered.
• Self-Interacting Dark Matter: These galaxies are less clustered at small scales. Because the centers are "smoothed out," they don't stick together as tightly as the standard model predicts.

The Big Takeaway

The paper concludes that how galaxies cluster together is a powerful detective tool.

• If we look at the universe and see that small galaxies are missing and the big ones are very tightly packed, it might point to Warm Dark Matter.
• If we see that the centers of galaxies are "puffy" and less dense, it might point to Self-Interacting Dark Matter.

The authors are saying: "We have these new recipes, and we've cooked up some virtual universes to test them. The results show that the way galaxies are distributed and clustered can tell us which recipe is the real one."

The Future: Right now, they are only looking at the invisible dark matter "fog." In the future, they plan to add the "visible" ingredients (stars, gas, black holes) to the simulation to see how the real, visible galaxies behave. This will help them compare their virtual cookies with the real ones we see in the sky, finally solving the mystery of what dark matter really is.

Technical Summary

Here is a detailed technical summary of the paper "The AIDA-TNG project. Abundance, radial distribution, and clustering properties of halos in alternative dark matter models."

1. Problem Statement

The standard $\Lambda$ Cold Dark Matter ($\Lambda$CDM) model successfully explains cosmological observations at large scales but faces persistent tensions at kiloparsec (kpc) scales, such as the "missing satellites," "too-big-to-fail," and "cusp-core" problems. These discrepancies suggest that the nature of dark matter (DM) might differ from the standard cold, collisionless assumption.
Alternative models, specifically Warm Dark Matter (WDM) and Self-Interacting Dark Matter (SIDM), have been proposed to resolve these issues. However, distinguishing between these scenarios requires precise characterization of how they affect the abundance, radial distribution, and clustering of dark matter halos and subhalos. While previous studies have touched on these aspects, a comprehensive analysis using the Halo Occupation Distribution (HOD) formalism across multiple resolutions and redshifts in a unified simulation suite was lacking.

2. Methodology

The authors utilized the AIDA-TNG (Alternative Dark Matter in the TNG universe) project, a suite of cosmological N-body simulations designed to compare different DM scenarios.

• Simulations: The study focused on dark matter-only runs of boxes with side lengths of 35 $h^{-1}$ Mpc (labeled 50/A and 50/B) and 75 $h^{-1}$ Mpc (100/A).
  • Models Tested:
    • CDM: Standard Cold Dark Matter.
    • WDM: Thermal relic models with particle masses $m_X = 1, 3, 5$ keV.
    • SIDM: Constant cross-section ($\sigma/m = 1$ cm$^2$ g$^{-1}$) and velocity-dependent cross-section (vSIDM) models.
  • Redshifts: Analyses were conducted at $z = 0, 0.5, 1, 2$.
• Halo Identification: Dark matter halos and subhalos were identified using the Friends-of-Friends (FoF) and Subfind algorithms.
• Analytical Framework:
  • Halo Occupation Distribution (HOD): The authors modeled the average number of halos ($\langle N_h \rangle$) as a sum of central ($\langle N_c \rangle$) and satellite ($\langle N_s \rangle$) components. They fitted the parameters $M_1$ (normalization) and $\alpha$ (slope) to describe the satellite population.
  • Radial Density Profiles: Recognizing that the standard Navarro-Frenk-White (NFW) profile is insufficient for subhalos, the authors employed a generalized NFW (gNFW) profile:
    $$\rho(r) = \frac{\rho_s}{(r/r_s)^\gamma (1 + r/r_s)^{3-\gamma}}$$
    where $\gamma$ corrects the inner slope and $f_c$ scales the concentration relative to the host halo.
  • Clustering Analysis: The two-point correlation function $\xi(r)$ was decomposed into one-halo ($\xi_{1h}$) and two-halo ($\xi_{2h}$) terms. The analysis included corrections for integral constraints (finite box size effects) and nonlinear bias.
• Statistical Approach: A Bayesian inference framework with Gaussian likelihoods was used to constrain HOD parameters ($M_1, \alpha$) and gNFW parameters ($\gamma, f_c$).

3. Key Contributions

1. Unified HOD Analysis: The paper provides the first comprehensive application of the HOD formalism to distinguish between WDM, SIDM, and CDM using the AIDA-TNG simulation suite, explicitly separating central and satellite contributions.
2. Subhalo Profile Characterization: It demonstrates that the standard NFW profile fails to accurately describe subhalo distributions and establishes the necessity of the gNFW model, quantifying the parameters $\gamma$ and $f_c$ across different DM models.
3. Clustering as a Discriminator: The study identifies small-scale clustering (specifically the one-halo term) as a powerful probe for distinguishing alternative DM models, often outperforming simple abundance counts in statistical significance.
4. Redshift Evolution: The work details how the differences between DM models evolve with redshift, highlighting that WDM effects are more pronounced at high $z$ due to delayed structure formation, while SIDM effects are driven by velocity-dependent thermalization.

4. Key Results

A. Abundance and Halo Occupation

• Satellite Suppression: Both WDM and SIDM models show a reduction in the average number of satellites compared to CDM, particularly in massive halos.
• HOD Parameters:
  • WDM: The suppression of low-mass halos due to free-streaming leads to a higher normalization parameter $M_1$ (requiring more massive hosts to host the same number of satellites) and a slight decrease in the slope $\alpha$.
  • SIDM: The reduction in satellites is interpreted as a consequence of core thermalization, which redistributes satellites to outer regions.
  • Distinction: While abundance counts alone struggle to distinguish models (especially at lower resolutions), the HOD parameters derived from clustering show a $>2\sigma$ difference between CDM and WDM3 for the parameter $M_1$ at $z=0$ and $z=2$.

B. Radial Distribution (gNFW Parameters)

• WDM (Cuspy): WDM models exhibit a steeper inner slope ($\gamma$ is higher) and a more concentrated distribution of satellites. This is attributed to the delayed formation of structures in WDM, allowing subhalos to form in denser environments and sink further via dynamical friction.
• SIDM (Cored): SIDM models exhibit a shallower inner slope ($\gamma$ is lower) and a more extended core. This is a direct result of particle self-interactions thermalizing the core and making subhalos more susceptible to tidal disruption.
• Mass Dependence: The inner slope $\gamma$ increases with host halo mass ($M_{200c}$), indicating that massive groups have steeper satellite profiles. However, in SIDM, this trend is reversed for low-mass groups which appear more peaked than massive ones due to earlier thermalization in massive systems.
• Concentration ($f_c$): The concentration factor $f_c$ remains roughly constant across redshifts and models, generally $>1$, reflecting the orbital decay of subhalos due to dynamical friction.

C. Clustering Properties

• One-Halo Term: The small-scale clustering ($r \lesssim 1$ Mpc $h^{-1}$) is dominated by satellite-satellite pairs.
  • WDM: Shows an excess of clustering at small scales compared to CDM due to the peaked subhalo distribution.
  • SIDM: Shows a deficit of clustering (less clustered) at small scales due to the cored, shallower distribution.
• Two-Halo Term: WDM models show enhanced clustering on larger scales due to the suppression of low-mass halos, which increases the effective bias of the remaining sample.
• Redshift Dependence:
  • WDM differences are more visible at high redshift ($z=2$) where structure formation delays are most acute.
  • SIDM differences are relatively constant with redshift, though velocity-dependent models (vSIDM) show stronger deviations at high $z$ where velocity dispersions are lower.

5. Significance and Future Outlook

This study establishes that small-scale clustering statistics, when combined with HOD modeling and generalized density profiles, are a robust tool for constraining the nature of dark matter. The ability to distinguish between WDM and SIDM scenarios using the two-point correlation function suggests that future observational surveys (e.g., Euclid, LSST) focusing on galaxy clustering can test these alternative models.

The authors note that current results are based on dark-matter-only simulations. The next step involves incorporating baryonic physics (using the full-physics AIDA-TNG runs) to disentangle the effects of dark matter properties from baryonic feedback (e.g., supernova winds, AGN), which also alter halo density profiles. This work lays the necessary theoretical groundwork for interpreting upcoming high-precision observational data.