Constraints on Self-Interacting Fuzzy Dark Matter from the Stellar Kinematics of the Dwarf Galaxy Leo II

This study uses stellar kinematics from the dwarf galaxy Leo II to constrain the two-dimensional parameter space of self-interacting fuzzy dark matter, revealing that attractive (repulsive) self-interactions lead to more concentrated (diffuse) density profiles and establishing 95% confidence lower limits on the particle mass within the range of $(1-10)\times10^{-22}\,\mathrm{eV}$ for interaction strengths up to $f_a^{-1}\lesssim 10^{-14}\,\mathrm{GeV}^{-1}$.

The Gist

Imagine the universe is filled with invisible "dark matter" that holds galaxies together. For a long time, scientists thought this matter was made of heavy, slow-moving particles (like cold, invisible rocks). But a newer theory, called Fuzzy Dark Matter (FDM), suggests it's actually made of incredibly light, wave-like particles. Think of these particles not as rocks, but as a giant, invisible fog or a vibrating string that stretches across the galaxy.

However, this "Fuzzy" theory has been getting into trouble. When scientists looked at small, lonely galaxies (dwarf galaxies) like Leo II, the math didn't add up. The "fog" was supposed to be spread out, but the stars inside Leo II were moving in a way that suggested the dark matter was clumped up tighter than the simple Fuzzy theory allowed. It was like trying to fit a fluffy cloud into a tiny jar; the cloud kept getting squished in ways the theory said it shouldn't.

The New Idea: The "Social" Fog
The authors of this paper asked: "What if these fuzzy particles aren't just floating around alone? What if they can talk to each other?"

In physics, this is called Self-Interaction (SI).

• Repulsive Interaction: Imagine the particles are like magnets with the same pole facing each other. They push away from one another. This makes the "fog" spread out even more, becoming very diffuse and fluffy.
• Attractive Interaction: Imagine the particles are like magnets with opposite poles. They pull toward each other. This makes the "fog" clump together, becoming denser and more compact in the center.

The Experiment: Testing Leo II
The team used the dwarf galaxy Leo II as their laboratory. They looked at how the stars in Leo II were moving (their "kinematics"). By measuring these speeds, they could map out exactly how much dark matter was where.

They then ran a simulation with three scenarios:

1. No Interaction: The standard Fuzzy theory (just the wave).
2. Repulsive Interaction: The particles push each other apart.
3. Attractive Interaction: The particles pull each other together.

The Results: Finding the Sweet Spot
Here is what they found, using simple analogies:

• The "No Interaction" Problem: Without any social interaction, the Fuzzy dark matter creates a core that is too big and too fluffy for Leo II. The stars are moving too fast for such a spread-out cloud. To fix this, the particles would have to be heavier than the theory allows, which creates a conflict with other observations.
• The "Repulsive" Worsening: If the particles push each other away, the cloud gets even fluffier. This makes the mismatch with the star movements even worse. It's like trying to stuff a giant beach ball into a shoebox; it just won't fit.
• The "Attractive" Solution: If the particles pull toward each other, the cloud shrinks and gets denser in the middle. This "clumping" matches the star movements in Leo II much better. It's like compressing that beach ball until it fits perfectly into the shoebox.

The Conclusion
The paper concludes that the simple "Fuzzy Dark Matter" theory is too rigid. However, if we add a "social" element where the particles attract each other, the theory can survive the test.

They found a "Goldilocks zone" for the strength of this attraction. If the attraction is too weak, the theory still fails. If it's too strong, it might break other rules. But within a specific range of attraction strength, the theory works, and the mass of the particles can be in a range that fits the data from Leo II.

In a Nutshell:
The universe's invisible "fog" might be too fluffy on its own. But if those invisible particles have a tendency to huddle together (attract), they can form a dense enough core to explain why the stars in Leo II move the way they do. This study doesn't prove the theory is right, but it shows that adding a little bit of "attraction" between the particles saves the theory from being ruled out by this specific galaxy.

Technical Summary

Technical Summary: Constraints on Self-Interacting Fuzzy Dark Matter from the Stellar Kinematics of the Dwarf Galaxy Leo II

Problem Statement
The one-parameter Fuzzy Dark Matter (FDM) model, characterized by an ultralight boson mass $m_a$, has faced increasingly stringent constraints from Lyman-$\alpha$ forest observations and local dwarf galaxy measurements. These constraints often challenge the viability of the standard non-interacting FDM paradigm. A natural theoretical extension to mitigate these limits is the inclusion of self-interactions (SI) within the FDM field, introducing a coupling parameter $f_a$. While self-interacting FDM (SIFDM) offers a broader parameter space to evade existing bounds, systematic and robust constraints on the two-dimensional parameter space $(m_a, f_a)$ remain limited. This study addresses the gap in constraining SIFDM using stellar kinematic data from the dwarf spheroidal galaxy Leo II, a system previously employed to constrain non-interacting FDM.

Methodology
The analysis proceeds through a three-step framework combining observational data analysis with theoretical reconstruction:

1. Jeans Analysis of Leo II: The authors utilize high-quality kinematic data from 175 confirmed member stars of Leo II. They perform a Jeans analysis using the general coreNFWtides parameterization for the dark matter (DM) density profile. This approach allows for the inference of a distribution of allowed DM density profiles without assuming a specific cosmological or galaxy evolution history. The analysis accounts for tidal stripping via a modified coreNFW profile and treats parameters such as virial mass ($M_{200}$), concentration ($c_{200}$), core radius ($r_c$), and velocity anisotropy ($\beta$) as free parameters in a Markov Chain Monte Carlo (MCMC) sampling. The result is a statistical ensemble of $\sim 1.6 \times 10^5$ density profiles, with a median profile serving as the fiducial input.

2. Reconstruction of SIFDM Wavefunctions: For a given set of parameters $(m_a, f_a)$, the authors reconstruct the corresponding SIFDM wavefunction $\psi$ that satisfies the Gross-Pitaevskii-Poisson equations. This is achieved via an eigenstate decomposition method. The input density profile (the median from the Jeans analysis) is used to define the effective gravitational and interaction potentials. The wavefunction is expressed as a linear combination of eigenstates, with expansion coefficients determined by fitting the time-averaged density profile of the reconstructed wavefunction to the input profile using a non-negative least-squares algorithm.

3. Statistical Constraints: Three distinct statistical methods are employed to quantify the consistency between the reconstructed SIFDM density profile ($\rho_{out}$) and the input profile ($\rho_{in}$) derived from the Jeans analysis:

   • Method 1 & 2 (Distance Measures): These methods calculate a symmetric logarithmic distance $d$ between $\rho_{out}$ and $\rho_{in}$ over specific radial intervals ($[r_{1/2}/4, 4r_{1/2}]$ and $[r_{1/2}/4, r_{1/2}]$). A parameter set is excluded at the 95% confidence level if the distance exceeds the 95th percentile of the distribution of distances calculated among the Jeans-inferred profiles.
   • Method 3 (Lower Bound): This method defines a lower bound curve corresponding to the lowest 5% of the density distribution from the Jeans analysis. A model is considered consistent only if $\rho_{out}$ remains entirely above this curve within the radial interval.

Key Contributions and Results
The study maps the constraints on the SIFDM parameter space $(m_a, f_a)$, revealing the interplay between particle mass and interaction strength:

• Effect of Self-Interactions: The analysis demonstrates that attractive SI leads to a more concentrated (compact) soliton core and higher central density, thereby improving agreement with the relatively cuspy density profiles inferred from the Leo II kinematics. Conversely, repulsive SI results in more diffuse cores and lower central densities, worsening the agreement with the data.
• Constraints on Non-Interacting FDM: In the absence of SI, the 95% confidence-level lower limits on the dimensionless mass parameter $m_{22} \equiv m_a / 10^{-22}$ eV are found to be 4.1, 6.3, and 7.5, depending on the statistical method used. These limits are consistent in order of magnitude with previous studies but slightly weaker, attributed to variations in statistical methodology.
• Constraints on SIFDM:
  • For a fixed particle mass, stronger attractive SI improves the fit, while stronger repulsive SI degrades it.
  • For a fixed SI strength, increasing the particle mass improves the fit.
  • The results indicate that for SI strengths satisfying $f_a^{-1} \lesssim 10^{-14} \text{ GeV}^{-1}$, the allowed range of particle masses lies broadly within $(1 - 10) \times 10^{-22}$ eV.
  • Specifically, for smaller particle masses, sufficiently strong attractive SI is required to maintain consistency with the data, whereas for larger masses, both weakly repulsive and attractive interactions remain allowed.

Significance
The paper claims that this analysis offers a complementary probe to existing constraints on FDM. By simultaneously constraining the two parameters $(m_a, f_a)$ without relying on assumptions regarding cosmological or galaxy evolution histories, the study provides a model-independent test of the SIFDM framework. The results highlight that self-interactions play a qualitatively distinct and physically transparent role in shaping the inner halo structure, effectively shifting the viable parameter space of FDM models. The work underscores the sensitivity of dynamical constraints to the weighting of inner versus outer halo regions and the specific statistical methods employed.