Statistical imprints of wave-like dark matter on multiply-imaged galaxies in strong cluster lenses from JWST

This paper demonstrates that the residual power spectrum of multiply-imaged galaxies in strong cluster lenses, as observed by JWST, can statistically distinguish wave-like dark matter from standard cold dark matter by quantifying sub-galactic density fluctuations across the $1 \lesssim k \lesssim 11\,\mathrm{kpc}^{-1}$ scale.

The Gist

The Big Picture: Hunting for "Ghost" Waves in the Dark

Imagine the universe is filled with invisible "dark matter" that holds galaxies together. For decades, scientists thought this stuff was made of tiny, solid particles (like invisible marbles) called Cold Dark Matter (CDM). But recent experiments haven't found these marbles, and some galaxy observations don't quite fit the "marble" theory.

This paper proposes a different idea: Wave-like Dark Matter ($\psi$DM). Instead of solid marbles, this theory suggests dark matter is a giant, fuzzy wave (like a sound wave or a ripple in a pond) that is so light it behaves like a wave over huge distances.

The authors ask: Can we tell the difference between "invisible marbles" and "fuzzy waves" by looking at how gravity bends light?

The Setup: A Cosmic Funhouse Mirror

To answer this, the team used the James Webb Space Telescope (JWST) as their tool. They focused on galaxy clusters—massive groups of galaxies that act like giant, natural magnifying glasses.

• The Analogy: Imagine looking at a distant streetlight through a funhouse mirror. The mirror (the galaxy cluster) bends the light, stretching the streetlight into long, curved arcs.
• The Goal: If the mirror is perfectly smooth, the arc looks smooth. But if the mirror has tiny bumps or ripples (caused by dark matter), the arc will have tiny wiggles or distortions.

The Method: Listening for the "Static"

The researchers simulated what JWST would see if the universe were filled with "marbles" (CDM) versus "waves" ($\psi$DM). They then created a computer program to analyze these images.

1. The Smooth Model: First, the computer tries to draw a perfect, smooth curve that matches the bent light (the arc). It assumes the dark matter is a smooth, invisible sheet.
2. The Residuals (The Leftovers): After the computer draws its perfect smooth curve, it subtracts that from the actual image. What's left over? The "residuals."
   • The Analogy: Imagine you are trying to trace a perfect circle on a piece of paper. If the paper has a tiny wrinkle, your pen will wobble. The "residual" is the wobble.
3. The Power Spectrum ($P_\delta(k)$): The team didn't just look at the wobbles with their eyes; they measured the "static" or "noise" in the wobbles. They used a mathematical tool called a Power Spectrum to see if the wobbles were random (like static on an old TV) or if they had a specific pattern (like a rhythmic hum).

The Discovery: Waves Leave a Different Fingerprint

The paper found that the "fuzzy waves" and the "solid marbles" leave very different fingerprints in the residuals:

• Cold Dark Matter (Marbles): The wobbles are small, random, and scattered. It's like static on a TV screen—chaotic and unorganized.
• Wave-like Dark Matter (Fuzzy Waves): The wobbles are coherent. Because the dark matter is a wave, it creates interference patterns (like ripples in a pond where waves crash into each other). This creates large, organized patches of wiggles that stretch across the image.

The Key Finding:
The team simulated deep observations (20 hours of looking at the sky). They found that:

• If the dark matter waves are very light (specifically, with a mass around $10^{-23}$ eV), the "organized wiggles" are so strong that JWST can easily spot them. The "static" looks different than it would if the universe were full of marbles.
• Even if the waves are a bit heavier, if the "ripples" are strong enough, JWST can still distinguish them from the marble theory.

The "Systematic Noise" Problem

The authors were very careful. They admitted that their computer models aren't perfect. Sometimes, the computer makes mistakes in drawing the smooth curve, creating "fake wiggles" that look like dark matter.

• The Analogy: Imagine you are trying to hear a whisper in a noisy room. The "noise" is the computer's imperfection.
• The Result: They found that for the "marble" theory, the real signal is hidden under the computer's noise. But for the "wave" theory (with the right mass), the signal is so loud that it pops right out above the noise, even with a single 20-hour observation.

Conclusion: A New Way to Listen to the Universe

The paper concludes that by looking at the "wiggles" in the light of distant galaxies, we can statistically tell if dark matter is made of particles or waves.

• If the wiggles are organized and large: It supports the Wave-like Dark Matter theory.
• If the wiggles are random and small: It supports the standard Cold Dark Matter theory.

This method doesn't require finding a single, specific dark matter particle. Instead, it listens to the "hum" of the entire population of dark matter, offering a new, independent way to solve one of the biggest mysteries in physics.

Technical Summary

Technical Summary: Statistical Imprints of Wave-like Dark Matter on Multiply-Imaged Galaxies in Strong Cluster Lenses from JWST

Problem Statement
The standard Cold Dark Matter (CDM) paradigm faces persistent challenges at sub-galactic scales, including the cusp–core, missing-satellite, and too-big-to-fail problems. While baryonic feedback may alleviate some tensions, the non-detection of Weakly Interacting Massive Particles (WIMPs) has motivated alternative theories. Wave-like dark matter ($\psi$DM), also known as fuzzy or ultralight dark matter, proposes that dark matter consists of extremely light bosons ($m \sim 10^{-22}$ eV). This model predicts a kpc-scale de Broglie wavelength, leading to wave-interference density fluctuations and solitonic cores that differ fundamentally from the discrete subhalo population of CDM.

Existing constraints on $\psi$DM particle mass ($m_\psi$) rely on individual systems (e.g., flux-ratio anomalies in quads, VLBI lensing) and are often vulnerable to systematic uncertainties in baryonic physics or incomplete wave-interference simulations. There is a need for independent, statistical probes that can distinguish $\psi$DM from CDM across a population of systems without requiring the resolution of individual subhalos.

Methodology
The authors propose using the residual power spectrum, $P_\delta(k)$, as a statistical probe. This metric quantifies deviations between observed lensed surface brightness and the prediction of a smooth lens model. The core hypothesis is that the amplitude and shape of $P_\delta(k)$ encode the population statistics of dark matter substructure along the line of sight, allowing differentiation between the granular, wave-interference fluctuations of $\psi$DM and the discrete subhalos of CDM.

The study employs a comprehensive simulation and analysis pipeline:

1. Mock Observations: The authors generate JWST-quality mock observations of strong gravitational lensing in galaxy clusters. They use real COSMOS galaxy images as sources and model the cluster lens as an elliptical Navarro-Frenk-White (NFW) profile ($M_{200} = 10^{15} M_\odot$, $z_{lens}=0.4$).
2. Substructure Modeling: Using the pyHalo package, they simulate line-of-sight substructure under three paradigms:
   • CDM: A population of subhalos and line-of-sight halos following the Sheth-Tormen mass function, subject to tidal stripping.
   • $\psi$DM (Case 1): $m_\psi = 10^{-22}$ eV with a fluctuation amplitude $100\times$ the physical prediction (to test detectability limits).
   • $\psi$DM (Case 2): $m_\psi = 10^{-23}$ eV with the physical fluctuation amplitude ($\log_{10} A_\psi = -0.8$) derived from numerical simulations.
3. Observational Simulation: Mocks simulate 20-hour deep exposures with JWST NIRCam (F150W filter), incorporating realistic noise, Point Spread Function (PSF), and detector specifications.
4. Lens and Source Modeling: The authors perform a joint Bayesian fit to the mock data using:
   • Source: A shapelet basis ($n_{max}=20$) to model the complex morphology of background galaxies.
   • Lens: The Curved Arc Basis (CAB) formalism, a local lensing approximation that describes distortions near each arc without requiring a global parametric cluster mass model.
5. Power Spectrum Analysis: After fitting, normalized residual maps are generated. The 2D power spectrum of these residuals, $P_\delta(k)$, is computed and averaged over 20 source realizations (60 residual maps total) to suppress variance from individual source morphologies and stochastic halo realizations.

Key Contributions

• Application of $P_\delta(k)$ to Cluster Lenses: While the substructure power spectrum has been proposed for galaxy-galaxy lenses, this work extends the methodology to strong cluster lenses, leveraging the high magnification and multiple images available in cluster fields to reduce statistical variance.
• Integration of Realistic Systematics: The study explicitly accounts for modeling systematics by comparing results against "smooth" mocks (no substructure). It demonstrates that the CAB + shapelet pipeline preserves the $\psi$DM signal while establishing a "modeling-systematics floor" that serves as the baseline for detection.
• Direct Measurement from Data: The paper demonstrates that $P_\delta(k)$ can be measured directly from data using local lens modeling (CAB), avoiding the complexities and degeneracies of global parametric cluster modeling.

Results

1. Sensitivity to Parameters: $P_\delta(k)$ is sensitive to both the boson mass ($m_\psi$) and the fluctuation amplitude ($A_\psi$).
   • Amplitude Dependence: Increasing $A_\psi$ monotonically raises the power spectrum amplitude without changing its spectral shape. At high amplitudes ($\log_{10} A_\psi = 1.2$), $\psi$DM produces a broadband power excess $\sim 3-10\times$ above the CDM baseline.
   • Mass Dependence: The detectability depends on the interplay between subhalo suppression and wave-fluctuation scale.
     • At $m_\psi = 10^{-21}$ eV, the de Broglie wavelength is sub-pixel, and the signal is indistinguishable from CDM.
     • At $m_\psi = 10^{-22}$ eV, subhalos are suppressed, but fluctuations remain small-scale and low-amplitude; the signal falls below or is indistinguishable from the CDM + systematics baseline in the full fitting analysis.
     • At $m_\psi = 10^{-23}$ eV, the de Broglie wavelength ($\sim 0.78$ kpc) is resolved, and fluctuations are strong enough to produce a clear excess over CDM.
2. Statistical Separation:
   • For the physical amplitude scenario ($m_\psi = 10^{-23}$ eV), the $\psi$DM and CDM $P_\delta(k)$ bands are non-overlapping across the range $1 \lesssim k \lesssim 11$ kpc$^{-1}$.
   • For the high-amplitude scenario ($m_\psi = 10^{-22}$ eV, $100\times$ amplitude), a statistically significant separation is also achieved across $k \sim 2-11$ kpc$^{-1}$.
3. Nature of the Signal: The $\psi$DM signature manifests as a broadband power-law modification rather than a sharp spectral peak at $k \sim 2\pi/\lambda_{dB}$. This is due to the projection of 3D density fluctuations onto the 2D lens plane, which smears the characteristic scale into a broadband excess.

Significance and Claims
The paper claims that the residual power spectrum is a viable, independent statistical probe for the wave-like nature of dark matter.

• Complementarity: This approach complements existing constraints from flux-ratio anomalies and VLBI lensing by probing the population-level statistics of substructure in cluster environments rather than individual systems.
• Feasibility: The authors demonstrate that a single 20-hour JWST observation of a strong cluster lens is sufficient to statistically distinguish $\psi$DM with $m_\psi \lesssim 10^{-23}$ eV (at physical amplitudes) or $m_\psi \sim 10^{-22}$ eV (at enhanced amplitudes) from CDM.
• Robustness: The method is robust against source morphology variations and local lens modeling uncertainties, provided the modeling systematics floor is well-characterized.
• Limitations: The authors note that the fluctuation amplitude in cluster environments is uncertain (extrapolated from galaxy-scale simulations), and the current analysis is limited to a single cluster configuration. They emphasize that for $m_\psi \gtrsim 10^{-22}$ eV, the signal may be suppressed below the current modeling-systematics floor, requiring improved source/lens modeling or compound lens configurations (e.g., galaxy-in-cluster lenses) to probe higher masses.

In conclusion, the work establishes $P_\delta(k)$ as a powerful tool for testing the wave-like nature of dark matter, offering a path to constrain $\psi$DM parameters using the rich statistical data expected from JWST cluster lensing programs.