A test of invariance of halo surface density for FIRE-2 simulations with cold dark matter and self-interacting dark matter

Using FIRE-2 simulations of dwarf galaxies, this study demonstrates that the dark matter halo surface density remains nearly constant across both cold and self-interacting dark matter models (with and without baryons), yielding results consistent with observational data and established scaling relations.

The Gist

Imagine the universe is a giant, invisible ocean made of "dark matter." We can't see this water, but we know it's there because it holds galaxies together like a giant, invisible net. For decades, scientists have been arguing about what this "net" looks like. Is it a smooth, dense ball in the center (a "cusp"), or is it a fluffy, spread-out cloud (a "core")?

This paper is like a team of detectives running a massive simulation to see if a specific rule about this invisible ocean holds true, even when we change the rules of the game.

The Mystery: The "Flat Fee" of the Universe

Scientists noticed something weird about galaxies. No matter how big or small the galaxy is, the "surface density" of its dark matter halo seems to be roughly the same.

The Analogy: Imagine you are buying a pizza. Usually, a bigger pizza costs more. But imagine if, no matter if you ordered a personal slice or a giant party pie, the amount of cheese per square inch was exactly the same. That's what the "constant surface density" rule suggests about dark matter. It's like the universe charges a "flat fee" for how much dark matter is packed into the center of a galaxy, regardless of the galaxy's size.

The Experiment: The "FIRE-2" Simulator

The authors used a super-computer simulation called FIRE-2 (Feedback In Realistic Environments). Think of this as a video game engine that builds tiny galaxies from scratch, complete with stars, gas, and gravity.

They tested three different versions of the game:

1. The Standard Model (CDM): The "Cold Dark Matter" we usually believe in. It's like a ghost that doesn't bump into anything.
2. The Bouncy Model (SIDM): "Self-Interacting Dark Matter." Here, the dark matter particles can bump into each other like billiard balls, potentially smoothing out the center of the galaxy.
3. The "No-Player" Mode (DMO): They ran the Bouncy model but turned off all the stars and gas to see what happens with just the dark matter.

They focused on dwarf galaxies, which are the small, scrappy cousins of big galaxies like our Milky Way.

The Investigation: Fitting the Puzzle Pieces

To measure the "flat fee" (surface density), the scientists had to fit mathematical shapes to the data coming out of the simulation. They used three different shapes (profiles) to describe the dark matter:

• The Burkert Profile: A shape that assumes a fluffy, flat center (a "core").
• The Core-Einasto & Alpha-Beta-Gamma: More complex shapes that can handle both fluffy centers and sharp spikes.

The Challenge: In the "Standard Model" (CDM) with stars, the dark matter often forms a sharp spike (a cusp) rather than a fluffy center. Trying to force a "fluffy" shape onto a "spiky" reality is like trying to fit a square peg in a round hole. The math didn't fit well for the smaller galaxies.

The Findings: The Rule Holds Up (Mostly)

Despite the fitting challenges, the results were fascinating:

1. The "Flat Fee" is Real: For the galaxies where the math worked, the "surface density" (the amount of dark matter packed in) was indeed nearly constant. Whether they used the Standard Model or the Bouncy Model, the result was the same. The universe seems to stick to this "flat fee" rule even for tiny dwarf galaxies.
2. Agreement with Reality: When they compared their simulation results to real observations of dwarf galaxies in our neighborhood (like those orbiting the Milky Way), the numbers matched up perfectly. The simulations predicted the same "flat fee" that astronomers see in the real sky.
3. The "Bouncy" vs. "Ghost" Debate: Interestingly, the simulations showed that the "Bouncy" dark matter (SIDM) fits the data slightly better than the "Ghost" dark matter (CDM) when it comes to creating those fluffy centers. However, even the "Ghost" model, once you account for how stars push the dark matter around, produces results that look very similar to the real world.

The Big Picture: Why This Matters

Think of this paper as a stress test for our understanding of the universe.

• The "Cusp vs. Core" Problem: For years, the standard model of the universe predicted sharp spikes of dark matter, but real galaxies look fluffy. This has been a major headache for physicists.
• The Verdict: This study suggests that even if the dark matter is the "Ghost" type (Standard Model), the presence of stars and gas (baryons) acts like a blender, smoothing out the spikes into fluffy centers. This makes the "flat fee" rule work.
• The Alternative: If dark matter is actually "Bouncy" (SIDM), it naturally creates fluffy centers without needing the stars to do the work.

Conclusion

In simple terms, the authors ran a cosmic simulation to see if the "flat fee" rule of dark matter holds up for small galaxies. It does.

Whether the dark matter is a passive ghost or an active bouncer, the universe seems to have a consistent way of packing it into the centers of dwarf galaxies. This gives scientists confidence that our current models of the universe are on the right track, even if we still need to figure out exactly how the dark matter behaves. It's a victory for the idea that the universe, in its own chaotic way, follows a surprisingly simple rule.

Technical Summary

Here is a detailed technical summary of the paper "A test of invariance of halo surface density for FIRE-2 simulations with cold dark matter and self-interacting dark matter" by Dalui and Desai.

1. Problem Statement and Motivation

The paper addresses the constancy of dark matter (DM) halo surface density, an observational anomaly where the product of the core radius ($r_c$) and central density ($\rho_c$) of cored DM halos appears nearly constant ($\log(\rho_c r_c) \approx 2.15 \pm 0.2$ in units of $M_\odot \text{pc}^{-2}$) across 18 orders of magnitude in galaxy luminosity.

While this invariance has been used to test alternatives to the standard $\Lambda$CDM model (e.g., Fuzzy Dark Matter), previous theoretical tests using mock catalogs were limited to massive halos ($10^{11} - 10^{13} M_\odot$) or lacked full baryonic physics. Furthermore, there is debate regarding whether this constancy holds for lower-mass dwarf galaxies ($\sim 10^{10} M_\odot$) and whether it persists in Self-Interacting Dark Matter (SIDM) scenarios.

Key Questions:

• Does the halo surface density remain constant for dwarf galaxy scales ($M_{halo} \approx 10^{10} M_\odot$) in high-resolution simulations?
• Does this invariance hold for both Cold Dark Matter (CDM) and Self-Interacting Dark Matter (SIDM) when baryonic feedback is included?
• How do these simulation results compare with recent empirical scaling relations proposed by Kaneda et al. (K24)?

2. Methodology

Data Source:
The authors utilized dark matter density profiles from the Feedback In Realistic Environments (FIRE-2) project, specifically the study by Straight et al. (S25). The dataset includes:

• 8 Dwarf Galaxies: Isolated halos with virial masses $\sim 10^{10} M_\odot$ (varying by $\pm 30\%$) and stellar masses $10^5 - 10^7 M_\odot$.
• Three Cosmological Scenarios:
  1. CDM + Baryons: Standard Cold Dark Matter with full baryonic physics.
  2. SIDM + Baryons: Self-Interacting Dark Matter ($\sigma/m = 1 \text{ cm}^2/\text{g}$) with full baryonic physics.
  3. SIDM DMO: SIDM without baryons (Dark Matter Only).

Fitting Procedures:
The authors fitted the simulation density profiles to three parametric models to extract structural parameters:

1. Burkert Profile: A cored profile ($\rho \propto (r^2+r_c^2)(r+r_c)^{-1}$). Used to directly compare with observational estimates of $\rho_c r_c$.
2. Core-Einasto Profile: A cored version of the Einasto profile.
3. $\alpha\beta\gamma$ Profile (gNFW): A generalized NFW profile.

Note: For CDM simulations, the Burkert profile provided poor fits (high residuals) for low-mass halos due to their "cuspy" nature. Consequently, only two CDM halos (m10k, m10m) met the quality threshold ($Q < 0.25$) for Burkert fitting.

Calculations:

• Surface Density ($S$): Calculated as $S = \rho_c \times r_c$ using the best-fit Burkert parameters.
• Column Density ($S^*$): Calculated by integrating the density profile along the line of sight up to a specific radius $R$ (defined differently for each profile, e.g., $R=1.66r_c$ for Burkert).
• Scaling Relations: The authors reconstructed the relation between mean surface density within the radius of maximum circular velocity ($\bar{\Sigma}(<r_{max})$) and $V_{max}$, comparing it to observational data from Milky Way and M31 dwarfs.

3. Key Contributions

1. Extension to Dwarf Scales: This is one of the first systematic tests of halo surface density invariance specifically at the dwarf galaxy scale ($10^{10} M_\odot$) using simulations that include full baryonic feedback (stellar winds, supernovae, etc.), which is crucial for creating cored profiles in CDM.
2. Comparison of DM Models: It provides a direct comparison between CDM and SIDM (with and without baryons) regarding the invariance of surface density, testing if the "constancy" is a generic feature of DM physics or specific to certain interaction models.
3. Validation of K24 Relations: The study validates the empirical scaling relations proposed by Kaneda et al. (K24) using high-fidelity FIRE-2 data, specifically testing the "cusp-to-core" transformation model against standard NFW relations.

4. Results

A. Constancy of Surface and Column Density

• Invariance: For all three simulation sets (CDM+baryons, SIDM+baryons, SIDM DMO), both the halo surface density and column density were found to be nearly constant across the mass range of the simulated dwarfs.
• Agreement with Observations: The median surface density derived from the Burkert profile fits is consistent with the observational estimate of $\rho_c r_c = (141 \pm 65) M_\odot \text{pc}^{-2}$ within $1\sigma$ for all models.
  • SIDM+baryons: $82^{+8}{-10} M\odot \text{pc}^{-2}$
  • SIDM DMO: $95^{+13}{-13} M\odot \text{pc}^{-2}$
  • CDM+baryons: $85 M_\odot \text{pc}^{-2}$ (mean of two valid halos).
• Profile Dependence: While the Burkert profile fits poorly for low-mass CDM halos (due to cusps), the column density calculated using Core-Einasto and $\alpha\beta\gamma$ profiles also showed near-constancy, suggesting the invariance is robust to the choice of density profile.

B. Scaling Relations ($\bar{\Sigma}$ vs. $V_{max}$)

• The authors derived the relation between mean surface density and maximum circular velocity ($V_{max}$).
• Observational Match: The simulated data points for CDM and SIDM dwarfs align well with observational data for Milky Way and M31 dwarfs.
• Model Preference: The simulated scaling relations agree more closely with the cusp-to-core transformation model (proposed by K24) than with the standard cuspy NFW profile based on the concentration-mass ($c-M$) relation. This suggests that baryonic feedback or self-interactions effectively transform cusps to cores, driving the observed scaling.

5. Significance and Conclusion

• Robustness of Invariance: The study confirms that the near-constancy of dark matter surface density is not limited to massive galaxies or specific DM models. It holds for dwarf galaxies in both CDM (with strong baryonic feedback) and SIDM scenarios.
• Baryonic Feedback Role: The results imply that baryonic feedback processes in CDM simulations are sufficient to create cores and maintain the surface density invariance, even without invoking exotic DM physics like SIDM.
• Theoretical Implications: The agreement with the K24 cusp-to-core transformation model supports the idea that the observed scaling relations are a natural outcome of galaxy formation physics (feedback) rather than a fundamental property of the dark matter particle itself.
• Limitations: The authors note that the dynamic range of halo masses in the FIRE-2 sample is limited, preventing a definitive test of mass dependence scaling. Future work will involve 2D projections of the full simulation data to refine surface density calculations.

In summary, the paper demonstrates that the "universal" dark matter surface density is a robust feature emerging from realistic galaxy formation simulations, consistent across different dark matter paradigms when baryonic physics is properly accounted for.