Dwarf Galaxy Constraints on Interacting Fermionic Dark Matter

This study utilizes stellar kinematic data from eight Milky Way dwarf spheroidal galaxies to constrain degenerate fermionic dark matter models, finding that fermion masses between 100 and 300 eV are favored and that current observations do not significantly distinguish between interacting and non-interacting equations of state.

The Gist

The Big Picture: What is Dark Matter?

Imagine the universe is like a giant, invisible ocean. We can't see the water (Dark Matter), but we can see the boats floating on it (stars and galaxies). For a long time, scientists have assumed this "ocean" is made of invisible, non-interacting particles that just float around, bumping into each other only through gravity.

But what if the water isn't just floating? What if the particles can actually "talk" to each other, pulling or pushing slightly? This paper asks: Can we tell if Dark Matter particles interact with each other by looking at the smallest galaxies in our neighborhood?

The Laboratory: Dwarf Galaxies

The authors chose "dwarf spheroidal galaxies" as their laboratory. Think of these as tiny, lonely islands of stars.

• Why them? They are mostly made of Dark Matter (like 99% of the island is invisible water, and only 1% is the visible boat).
• The Clue: We can't see the Dark Matter, but we can watch how the stars move. If the stars are moving fast, it means there is a lot of invisible weight holding them together. If they move slowly, there is less weight.

The Two Theories: The Crowd vs. The Crowd with Handshakes

The researchers tested two different ideas about how Dark Matter behaves:

1. The "Non-Interacting" Crowd (The Standard Model):
   Imagine a crowd of people in a room who are all very shy. They don't talk to each other, they don't hold hands, and they don't push. They only care about not bumping into each other because of a rule called the "Pauli Exclusion Principle" (like how you can't sit in a seat someone else is already sitting in). This creates a "pressure" that keeps the crowd from collapsing into a tiny pile. This is the standard theory.

2. The "Interacting" Crowd (The New Idea):
   Now, imagine that same crowd, but these people have a secret ability to gently pull on each other (attractive force).

   • The Effect: If they pull on each other, the crowd can get much tighter and denser. It's like the shy crowd suddenly deciding to huddle together for warmth. This changes how the "pressure" works.

The Experiment: Measuring the "Squishiness"

The authors created a mathematical model to see what happens if Dark Matter has this "huddling" ability.

• The Analogy: Think of the Dark Matter halo as a giant, invisible balloon.
  • In the Standard Model, the balloon is stiff. It resists being squished.
  • In the Interacting Model, the balloon is softer in some spots. If you squeeze it, it collapses more easily, making the center very dense and the edges very thin.

They solved complex equations (like a very advanced weather forecast for gravity) to predict how the stars should move in these two different types of "balloons."

The Results: What the Data Says

The team took real data from eight dwarf galaxies (like Carina, Draco, and Fornax) and compared the star movements to their two models.

1. The Mass Guess: Both models agreed on one thing: The Dark Matter particles must be very light, roughly 100 to 300 electron-volts (which is incredibly light, about a million times lighter than a proton).
2. The Interaction Test: This is the big finding.
   • The data did not show a clear preference for the "huddling" (interacting) crowd over the "shy" (non-interacting) crowd.
   • The "shy" crowd model fits the data just as well as the "huddling" crowd model.
   • The Catch: If the "huddling" force were too strong, the galaxies would collapse into tiny, dense balls, and the stars would move in ways that don't match what we see.

The Conclusion

The paper concludes that current observations don't prove that Dark Matter particles interact with each other.

• The Verdict: The "shy crowd" (non-interacting) model is still the best description we have.
• The Limit: We can say that if Dark Matter does have a "huddling" force, it must be very weak. If it were strong, the galaxies would look different than they do.

In short: The authors looked at the smallest galaxies to see if Dark Matter particles hold hands. They found no evidence that they do. The particles seem to be just as shy and independent as we thought they were, and any "huddling" they might do must be very, very subtle.

Technical Summary

Technical Summary: Dwarf Galaxy Constraints on Interacting Fermionic Dark Matter

Problem Statement
The microscopic nature of dark matter (DM) remains an open question. While the standard Cold Dark Matter (CDM) model succeeds on large scales, small-scale structures like dwarf spheroidal galaxies (dSphs) offer a unique laboratory to probe DM particle properties, such as mass, spin, and self-interactions. This work focuses on fermionic DM candidates, where degeneracy pressure can stabilize gravitational collapse and generate extended cores. Previous studies have largely treated DM as a non-interacting, degenerate Fermi gas. However, this minimal scenario neglects potential attractive self-interactions within the dark sector (e.g., mediated by Yukawa-like forces), which could significantly alter the equation of state (EoS), compressibility, and stability of DM halos. The central question addressed is whether current stellar kinematic data from dwarf galaxies are sensitive to these qualitative EoS features beyond the free degenerate Fermi-gas limit.

Methodology
The authors employ a multi-step approach combining hydrostatic modeling, stellar kinematics, and statistical inference:

1. Equation of State (EoS) Formulation:

   • Baseline: The standard non-relativistic degenerate Fermi gas EoS, $P(\rho) = K\rho^{5/3}$.
   • Interacting Model: A phenomenological EoS incorporating attractive self-interactions:
     $$P(\rho) = K\rho^{5/3} - \alpha \frac{\rho^2}{\rho_s^2 + \rho^2}$$
     Here, $\alpha$ parametrizes the interaction strength, and $\rho_s$ defines the density scale where finite-density effects saturate. This form allows for a "softening" of the EoS and, in certain parameter regions, a negative compressibility ($dP/d\rho < 0$), signaling mechanical instability.
   • Maxwell Construction: For regions where the EoS becomes mechanically unstable, the authors apply a Maxwell construction to enforce phase coexistence, replacing the unstable branch with a constant-pressure plateau. This ensures physically admissible hydrostatic solutions.

2. Hydrostatic Equilibrium:
   The authors solve the non-relativistic hydrostatic equations (coupled with the mass continuity equation) to determine the DM density profiles $\rho(r)$ and enclosed mass profiles $M(r)$. For the interacting case, solutions are computed numerically for each set of parameters $(m_\chi, \alpha, \rho_s, \rho_c)$, accounting for the Maxwell construction where necessary.

3. Kinematic Modeling:
   The predicted line-of-sight (LOS) velocity dispersion profiles, $\sigma_{\text{los}}(R)$, are derived using the spherical Jeans equation. The models assume a Plummer profile for the stellar tracer distribution and treat the stellar velocity anisotropy ($\beta$) as a free parameter.

4. Statistical Analysis:
   Markov Chain Monte Carlo (MCMC) fits are performed on LOS velocity dispersion data from eight classical Milky Way dwarf spheroidals: Carina, Draco, Fornax, Leo I, Leo II, Sculptor, Sextans, and Ursa Minor. The analysis compares the non-interacting model (parameters: $m_\chi, \rho_c, \beta$) against the interacting model (parameters: $m_\chi, \alpha, \rho_s, \rho_c, \beta$).

   • Constraints: Physical cuts are applied post-sampling to ensure dynamical stability (excluding regions where $dM/d\rho_c < 0$) and to satisfy the Tremaine-Gunn bound (phase-space density constraints).
   • Priors: Broad uniform priors are used for all parameters, with logarithmic sampling for densities and interaction scales.

Key Results

• Fermion Mass: The data robustly constrain the fermion mass to the range $m_\chi \sim 100\text{--}300$ eV. The preferred mass values are nearly identical for both the interacting and non-interacting models, indicating that the data primarily constrain the halo extent via degeneracy pressure.
• Interaction Strength: The analysis yields upper limits on the attractive interaction strength. For the phenomenological EoS, the 95% confidence level upper bound is typically $\log_{10}(\alpha/\text{eV}^4) \lesssim -13$. The strongest constraint comes from Leo II ($\le -13.80$).
• Equation of State Preference: The posterior distributions for fermion mass ($m_\chi$), central density ($\rho_c$), and stellar anisotropy ($\beta$) are broadly similar between the interacting and non-interacting models. The data do not strongly prefer an interacting EoS over the free degenerate Fermi gas. Instead, the interacting model is constrained by the absence of large deviations from the non-interacting limit; strong attractive interactions would render halos too compact or mechanically unstable, leading to poor fits with observed velocity dispersions.
• Degeneracy and Stability: The study highlights a degeneracy where a lighter fermion mass (which naturally creates larger cores) can be compensated by attractive interactions (which steepen profiles and reduce core sizes). However, the stability analysis reveals that strong interactions can lead to dynamically unstable configurations or require Maxwell constructions that exclude certain central densities.

Significance and Claims
The paper claims that current stellar kinematic data from classical dwarf spheroidals are sufficient to test and constrain phenomenological equations of state for interacting fermionic dark matter. The primary significance lies in demonstrating that:

1. Exclusion of Large Deviations: Current data exclude large deviations from the non-interacting limit. While attractive interactions are theoretically possible, they must be weak enough that the resulting halo structure remains consistent with the free Fermi-gas predictions.
2. Methodological Rigor: The work provides a consistent framework for solving hydrostatic equilibrium in the presence of phase transitions (Maxwell construction) and integrating these solutions into Jeans analysis for parameter estimation.
3. Complementarity: The approach complements previous phase-space arguments by starting directly from the DM EoS rather than assuming halo profiles, offering a direct probe of the dark sector's microphysics (specifically the interplay between degeneracy pressure and self-interactions).

The authors conclude modestly that while the data do not rule out interactions entirely, they place stringent upper bounds on their strength, effectively favoring the non-interacting degenerate Fermi gas scenario within the tested parameter space. Future improvements in observational data (e.g., larger spectroscopic samples, Gaia proper motions) and theoretical refinements (e.g., connecting the phenomenological EoS more explicitly to particle physics models) are suggested to sharpen these constraints.