🔭 astrophysics

Wavelet-Scattering Signatures of Fuzzy Dark Matter in Simulated 21 cm Brightness-Temperature Maps

This paper demonstrates that applying the wavelet scattering transform to simulated 21 cm brightness-temperature maps reveals complementary, non-Gaussian morphological signatures of fuzzy dark matter that, when combined with the power spectrum, yield tighter cosmological constraints than the power spectrum alone, even under realistic foreground contamination scenarios.

Original authors: Hayato Shimabukuro, Shihang Liu, Bohua Li

Published 2026-06-09

📖 4 min read☕ Coffee break read

CC0 1.0

Original authors: Hayato Shimabukuro, Shihang Liu, Bohua Li

Original paper dedicated to the public domain under CC0 1.0 (http://creativecommons.org/publicdomain/zero/1.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Big Picture: A Cosmic Mystery

Imagine the universe as a giant, dark ocean. For a long time, scientists believed this ocean was filled with invisible, heavy rocks (called "Cold Dark Matter") that clumped together to form islands (galaxies). But there's a problem: when we look at the smallest islands, there aren't as many as our rock theory predicts.

Enter Fuzzy Dark Matter (FDM). Instead of heavy rocks, this theory suggests the ocean is filled with a ghostly, ultra-light mist made of tiny waves. Because these waves are so light, they spread out and smooth things over, preventing the formation of tiny islands. This paper asks: Can we see the difference between the "rocky" universe and the "misty" universe?

The Tool: Listening to the Cosmic Dawn

To answer this, the scientists looked at the "Cosmic Dawn"—the time when the very first stars turned on. They used a special kind of radio signal called the 21 cm line. Think of this signal like a thermal camera for the early universe. It shows us where the gas is cold (absorbing light) and where it is hot (emitting light).

Usually, scientists look at this signal by measuring its "loudness" at different sizes (like checking the volume of a song at different frequencies). This is called the Power Spectrum. It's like listening to a song and just counting how many high notes vs. low notes there are.

The New Idea: The Wavelet Scattering Transform (WST)

The authors of this paper wanted to look deeper. They used a new tool called the Wavelet Scattering Transform (WST).

The Analogy:
Imagine you are trying to recognize a friend in a crowded room.

The Power Spectrum is like counting how many people are wearing red shirts vs. blue shirts. It gives you a general idea of the crowd's color, but it doesn't tell you who is standing next to whom.
The WST is like looking at the patterns of the crowd. It notices that "the person in the red shirt is always standing next to the person in the blue shirt," or "the group of red shirts is clustered in a specific shape." It captures the texture and the relationships between different parts of the image, not just the total amount of color.

What They Found

The scientists ran computer simulations of the early universe using both the "Rocky" (Cold Dark Matter) and "Misty" (Fuzzy Dark Matter) theories. They then applied their new "Pattern Scanner" (WST) to the radio maps.

The "Mist" Delays the Party: In the Fuzzy Dark Matter universe, the tiny waves smooth out the gas. This means the first stars form later than in the Rocky universe. The scientists saw this delay clearly in their data. The "party" (the formation of stars and heat) started later in the misty version.
Changing the Shape: The "Mist" doesn't just delay things; it changes the shape of the cosmic structures. The WST tool was able to spot these shape changes, especially in the fine details (the small-scale structures) that the standard "loudness" check (Power Spectrum) missed.
Better Together: When they combined the standard "loudness" check with their new "pattern" check, they got the best possible result. It's like having both a microphone (to hear the volume) and a camera (to see the shape). Together, they give a much clearer picture of what the universe is made of.

The "Noise" Problem

In real life, radio telescopes (like the SKA) have to deal with a lot of static and interference from Earth (like cell phones and TV signals). This is called "foregrounds."

The paper tested a "cleaning" method to remove this interference. They found that while the cleaning process changes the "loudness" numbers a bit, the pattern relationships (the ratios between different scales) remain surprisingly stable. This suggests that even with a messy telescope, the "pattern scanner" (WST) can still find the clues about Fuzzy Dark Matter.

The Conclusion

This paper doesn't claim to have found Fuzzy Dark Matter yet. Instead, it proves that our new "pattern scanner" (WST) is a powerful tool.

If we build the next generation of radio telescopes, using this new method alongside the old one will help us tell the difference between a universe made of heavy rocks and one made of ghostly waves. It's a new way to look at the cosmic map that reveals details we couldn't see before.

Technical Summary: Wavelet-Scattering Signatures of Fuzzy Dark Matter in Simulated 21 cm Brightness-Temperature Maps

Problem Statement
The nature of dark matter remains a central open question in cosmology. While the Cold Dark Matter (CDM) paradigm explains large-scale structure, it faces challenges on sub-galactic scales, such as the cusp-core problem and the missing-satellite problem. Fuzzy Dark Matter (FDM), composed of ultralight bosons ( $m \sim 10^{-22}$ eV), offers a theoretically motivated alternative where quantum pressure suppresses small-scale structure formation. Current constraints on FDM mass from Lyman- $\alpha$ forests and satellite counts are tightening, yet complementary probes at high redshift ( $10 \lesssim z \lesssim 30$ ) are needed. The redshifted 21 cm line of neutral hydrogen traces the thermal and ionization history of the intergalactic medium (IGM) during Cosmic Dawn and the Epoch of Reionization (EoR). Because the formation of the first stars depends on low-mass halos, the 21 cm signal is sensitive to the suppression of small-scale structure in FDM cosmologies. However, conventional analyses relying solely on the 21 cm power spectrum may not fully capture the non-Gaussian, multiscale morphological changes induced by FDM.

Methodology
This study employs a morphological extension of existing FDM 21 cm calculations rather than a new simulation framework. The authors utilize the semi-numerical code 21cmFAST, modified to incorporate FDM effects as detailed in Liu et al. (2025). These modifications include:

Linear Regime: A sharp cutoff in the linear matter power spectrum derived from the Schrödinger-Poisson system, characterized by a Jeans scale dependent on the particle mass $m_{22} = m_{\text{FDM}} / 10^{-22} \text{ eV}$ .
Nonlinear Regime: A modified halo mass function (HMF) that suppresses halos below the Jeans mass and a phenomenological fitting formula for the HMF shape, including an environmental modulation prescription for local halo abundance.

The authors generate simulated 21 cm brightness-temperature cubes ( $\delta T_b$ ) for CDM and various FDM models (varying $m_{22}$ and astrophysical parameters). They apply the two-dimensional Wavelet Scattering Transform (WST) to 2D slices of these maps. The WST provides translation-invariant, deformation-stable statistics:

First-order coefficients ( $S_1$ ): Measure localized, scale-dependent variance (analogous to a localized power spectrum).
Second-order coefficients ( $S_2$ ): Measure cross-scale couplings (non-Gaussian phase information).
Normalized ratios ( $R = S_2/S_1$ ): Quantify the relative strength of cross-scale interactions, reducing dependence on overall amplitude.

The analysis includes:

Fisher Forecasts: Comparing constraints from power-spectrum-only, WST-only ( $S_1$ , $R$ , or combined), and combined power-spectrum+WST data vectors. These use identical mock-observation assumptions (SKA1-Low-like), redshift bins, nuisance parameters, and thermal noise models.
Foreground Treatment: An idealized test applying wedge masks to 3D light-cone chunks to remove contaminated low- $k_\parallel$ modes before extracting WST summaries.

Key Contributions

Morphological Characterization: The paper introduces the WST as a tool to characterize how FDM-induced delays in source formation and small-scale suppression manifest in the morphology of 21 cm maps, moving beyond two-point statistics.
Complementarity Analysis: It explicitly tests whether WST summaries carry information complementary to the power spectrum. The study constructs Fisher forecasts to quantify the incremental information gain when combining WST with the power spectrum.
Foreground Robustness Check: It provides a diagnostic showing how the normalized second-order ratio ( $R$ ) responds to idealized foreground-wedge filtering, demonstrating that $R$ is less strongly reshaped by horizon-wedge filtering than first-order amplitudes ( $S_1$ ).

Results

WST Coefficients:
- $S_1$ (First-order): The redshift evolution of $S_1(j)$ mirrors power-spectrum trends. FDM models exhibit a systematic shift of $S_1$ peaks toward lower redshifts and a reshaping of peak profiles compared to CDM. This reflects the delayed build-up of Ly $\alpha$ coupling, X-ray heating, and reionization due to the suppression of low-mass halos. The effect is most pronounced at smaller wavelet scales (finer structures) and for lighter FDM masses ( $m_{22}=10$ ).
- $R$ (Second-order ratio): The normalized ratio $R(j_1, j_2)$ highlights cross-scale phase coupling. In light FDM models, the suppression of small-scale power ( $S_1$ ) is stronger than the reduction in cross-scale modulation ( $S_2$ ), leading to an increase in the normalized ratio $R$ at specific redshifts (intermediate $z \sim 10-15$ ). This indicates that while absolute power is damped, the relative coupling between scales remains distinct.
- Astrophysical Parameters: Variations in virial temperature ( $T_{\text{vir}}$ ), X-ray luminosity ( $L_X/\text{SFR}$ ), and ionizing efficiency ( $\zeta$ ) produce coherent shifts in the redshift evolution of both $S_1$ and $R$ , allowing them to serve as compact summaries of the timing of Cosmic Dawn and EoR.
Fisher Forecasts:
- The combined power-spectrum + WST data vector yields the tightest marginalized constraints on FDM parameters and astrophysical nuisance parameters.
- This confirms that WST summaries carry information complementary to the power spectrum within the adopted local Fisher setup.
Foreground Impact:
- The idealized wedge-avoidance test shows that while low- $k_\parallel$ modes are removed, the normalized second-order ratio $R$ is less severely reshaped by this filtering than the first-order amplitudes. However, the authors note this is not a full foreground-subtraction pipeline and does not constitute a foreground-aware detectability forecast.

Significance and Claims
The paper claims to establish the Wavelet Scattering Transform as an interpretable morphological summary that remains useful when compared directly with the power spectrum. It demonstrates that WST captures non-Gaussian, multiscale features of the 21 cm signal that are altered by FDM physics, specifically the delayed emergence of luminous sources and the suppression of small-scale structure.

The authors are modest about the scope of their work:

They do not claim to provide a final SKA mass reach or an end-to-end foreground-subtraction pipeline.
They utilize a minimal astrophysical parameter set and fixed cosmological parameters.
The observational model is a controlled mock setup rather than a full visibility-domain pipeline.

The primary significance lies in establishing that WST offers a unified, higher-order statistical framework capable of probing the microphysics of dark matter across the spatial and redshift scales traced by the 21 cm signal, providing constraints that are complementary to conventional power-spectrum analyses.