Confronting fuzzy dark matter with the rotation curves of nearby dwarf irregular galaxies

While fuzzy dark matter provides an excellent fit to the rotation curves of nearby dwarf irregular galaxies with an axion mass of approximately $2\times10^{-23}$ eV, the model is ultimately ruled out by significant tensions ($\gtrsim5\sigma$) regarding the observed core scaling relations and the suppression of the linear power spectrum.

The Gist

Imagine the universe is a giant, invisible ocean. For decades, scientists have believed this ocean is made of "Cold Dark Matter" (CDM)—a substance that acts like a swarm of tiny, invisible, non-interacting marbles. These marbles clump together to form the scaffolding for galaxies.

However, there's a problem. When scientists look at small, lonely galaxies (dwarf irregulars), the "marbles" seem to pile up too sharply in the center, like a steep mountain peak. But when they measure how stars and gas move in these galaxies, it looks like the center is flat, like a gentle hill. This is the "core-cusp" puzzle.

To fix this, some scientists proposed a new idea: Fuzzy Dark Matter (FDM). Instead of tiny marbles, imagine the dark matter is made of ultra-light waves, like ripples on a pond. Because these waves are so light and spread out, they can't pile up into a sharp peak; instead, they naturally form a smooth, flat "core" in the center of a galaxy. This wave-like behavior is called a "soliton."

This paper is a reality check. The authors took a very high-quality, clean dataset of 11 nearby dwarf galaxies (from the "LITTLE THINGS" survey) and asked: "Does the Fuzzy Dark Matter wave theory actually fit the data?"

Here is what they found, broken down into simple concepts:

1. The "Goldilocks" Mass

First, they tried to find the "weight" (mass) of these invisible waves. If the waves are too heavy, they act like marbles; if they are too light, they spread out too much.

• The Result: The data fit the Fuzzy Dark Matter model very well. They found a "sweet spot" for the mass of these waves: roughly $2 \times 10^{-23}$ electron-volts.
• The Catch: If FDM were the perfect answer, every galaxy should point to this exact same mass. Instead, the authors found a strange pattern: the "heavier" the galaxy's stars were, the "lighter" the dark matter wave seemed to be. It's as if the waves were changing their weight depending on the neighborhood they lived in, which shouldn't happen if they are a fundamental particle of the universe.

2. The "Wrong Shape" of the Core

The theory predicts specific rules for how the size of the flat core relates to its density and mass. Think of it like a recipe: "If you double the size of the cake, the density must drop by a specific amount."

• The Result: The galaxies in the study broke the recipe. The relationship between the size of the core and its mass was almost the opposite of what the Fuzzy Dark Matter theory predicted.
• The Analogy: Imagine a theory that says "The bigger the balloon, the lighter the air inside." But when the scientists measured the balloons, they found "The bigger the balloon, the heavier the air inside." The data was so different from the prediction that it was a statistical mismatch of over 5 standard deviations (a very strong "no").

3. The "Missing Galaxies" Problem

This is the biggest blow to the theory. Fuzzy Dark Matter acts like a filter. Because the waves are so big, they smooth out the universe on small scales, preventing tiny clumps from forming.

• The Theory: If the mass of the wave is the "sweet spot" the authors found ($2 \times 10^{-23}$ eV), the universe should be so smooth that tiny dwarf galaxies shouldn't exist at all. The waves would have washed them away before they could form.
• The Reality: We are looking at 11 of these tiny galaxies right now. They exist.
• The Conclusion: The mass required to make the rotation curves of these galaxies look "flat" (the wave theory) is the exact same mass that would prevent these galaxies from existing in the first place. It's a "Catch-22." To explain the shape of the galaxy, you need a wave mass that erases the galaxy.

4. Did Stars and Gas Mess It Up?

The authors wondered: "Maybe the stars and gas inside the galaxies are squishing the dark matter waves, changing the results?"

• The Result: They did the math to include the gravity of stars and gas. While it did change the numbers slightly, it wasn't enough to fix the problems. The "wrong shape" of the core and the "missing galaxies" paradox remained.

The Bottom Line

The paper concludes that while Fuzzy Dark Matter looks beautiful on paper and fits the shape of the rotation curves surprisingly well, it fails the "sanity checks."

1. The properties of the galaxy cores don't match the theoretical rules.
2. The mass required to explain the curves would have prevented these galaxies from forming in the first place.

In short, the "Fuzzy" wave theory might be a nice idea, but when tested against the real, messy data of nearby dwarf galaxies, it doesn't hold up. The universe seems to be more complex than a simple wave of invisible particles.

Technical Summary

Problem Statement
The fundamental nature of dark matter (DM) remains one of the most pressing open problems in contemporary physics. While the Cold Dark Matter (CDM) paradigm within the $\Lambda$CDM model successfully describes large-scale cosmological observations, it faces significant challenges at small galactic scales. These include the "core-cusp" puzzle, where CDM simulations predict density profiles peaking as $\rho(r) \propto 1/r$ (cusps), whereas observations of dwarf galaxies suggest flat density cores; the "missing satellites" problem; and the "too-big-to-fail" problem.

Fuzzy Dark Matter (FDM) has emerged as a prominent alternative candidate, proposing that DM consists of ultralight axions (mass $m_a \sim 10^{-23} - 10^{-21}$ eV) behaving as coherent waves. A key prediction of FDM is the formation of a solitonic core at the center of halos due to quantum pressure effects, which could naturally resolve the core-cusp problem. However, FDM also predicts a suppression of the linear matter power spectrum on small scales, potentially inhibiting the formation of low-mass halos. This paper aims to confront the predictions of FDM with high-resolution rotation curves of nearby dwarf irregular galaxies to determine if FDM is a viable solution to these small-scale issues or to set exclusion ranges for the axion mass.

Methodology
The authors utilize a homogeneous and robust sample of 11 isolated, dark matter-dominated dwarf-irregular galaxies from the "LITTLE THINGS in 3D" (LTs) catalog. These galaxies were selected from the LITTLE THINGS survey and re-analyzed using state-of-the-art 3D reconstruction techniques to obtain high-resolution HI rotation curves, excluding galaxies with biased dynamical analyses ("rogues").

The theoretical framework models the DM density profile as a piecewise function: a solitonic core (derived from the ground-state solution of the Schrödinger-Poisson equations) in the inner region, matched to a Navarro-Frenk-White (NFW) halo in the outer region. The model depends on four primary parameters: the axion mass ($m_a$), the solitonic core mass ($M_c$), the halo concentration parameter ($c$), and the total halo mass ($M_h$).

The authors employ a statistical framework based on $\chi^2$ analysis using Markov Chain Monte Carlo (MCMC) techniques (emcee). They fit the theoretical rotation curves (including contributions from DM, stars, and gas) to the observed data. The analysis incorporates loose priors for most parameters but enforces the "core-halo relation" predicted by FDM simulations (linking $M_c$ and $M_h$) and ensures the transition radius between the soliton and the halo is physically consistent. The study also investigates the effects of baryonic feedback by solving for modified soliton solutions in the presence of a static baryonic potential.

Key Contributions and Results

1. Excellent Fit to Rotation Curves: The FDM model provides an excellent fit to the rotation curves of the LTs sample. The inferred axion masses for individual galaxies cluster around a narrow range, with a weighted average of $m_a \approx 1.90^{+0.24}_{-0.21} \times 10^{-23}$ eV.
2. Inconsistency with Soliton Scaling Relations: Despite the good fit to the velocity curves, the authors identify a major tension regarding the internal properties of the fitted cores. The data follow scaling relations between core radius ($r_c$), core mass ($M_c$), and central density ($\rho_c$) that are inconsistent with FDM predictions.
   • Specifically, the $r_c - M_c$ relation in the data shows a correlation almost orthogonal to the FDM prediction (heavier cores are bigger, whereas FDM predicts heavier cores are smaller). This discrepancy has a statistical significance of $\gtrsim 5\sigma$.
   • Similarly, the $r_c - \rho_c$ relation disagrees with FDM predictions at a significance $> 3\sigma$.
   • Furthermore, the derived axion masses show a mild but significant ($\sim 3.3\sigma$) negative correlation with stellar mass ($M_\star$), suggesting that astrophysical processes beyond a pure FDM soliton are influencing the observed cores.
3. Tension with Halo Abundance (HMF): The second major problem arises from cosmological constraints. The axion mass preferred by the rotation curve fits ($m_a \approx 2 \times 10^{-23}$ eV) predicts a strong suppression of the linear power spectrum, leading to a severe cutoff in the Halo Mass Function (HMF) at low masses.
   • Comparing the predicted number of halos in a Local Group volume with the observed census of dwarf galaxies, the authors find a tension exceeding $5\sigma$.
   • Using a conservative lower bound based on the sample galaxies, they derive $m_a \gtrsim 4.3 \times 10^{-23}$ eV. Using a broader census of 116 galaxies, the lower bound tightens to $m_a \gtrsim 5.5 \times 10^{-22}$ eV. This creates a "catch-22": the mass required to fit the rotation curves is too low to allow the existence of the observed number of dwarf galaxies.
4. Baryonic Effects: The authors estimate the impact of baryons (stars and gas) on the soliton structure. They find that including baryonic potentials leads to more compact solitons, requiring a systematic negative correction to the axion mass (typically $\lesssim 15\%$) to reproduce the original density profiles. However, these corrections are insufficient to resolve the tensions regarding the scaling relations or the halo abundance suppression.

Significance and Conclusion
The paper concludes that while the FDM soliton+NFW model can mathematically fit the rotation curves of dwarf irregular galaxies, the physical parameters required to do so are inconsistent with the broader predictions of the FDM model. The study highlights a "catch-22" scenario: the axion mass that produces the observed core structures in dwarf galaxies simultaneously predicts a suppression of small-scale structure formation that is too strong to be consistent with the observed abundance of such galaxies.

The authors state that their analysis poses a "serious challenge" to FDM as a hypothesis for solving the core-cusp problem or other small-scale issues of $\Lambda$CDM. They suggest that while baryonic physics might alleviate some tensions, it is unlikely to resolve the fundamental discrepancies found. The paper implies that alternative approaches, such as self-interacting dark matter or models where FDM constitutes only a fraction of the total dark matter, may be necessary. The work underscores the need for further observational reanalysis of HI data from other isolated dwarf galaxies to increase the statistical power of such cosmological tests.