⚛️ phenomenology

Using Strong Lensing to Detect Subhalos with Steep Inner Density Profiles

This study demonstrates that strong gravitational lensing observations can detect dark matter subhalos with steep inner density profiles at masses over an order of magnitude lower than standard NFW subhalos, a capability that remains robust even when accounting for complex lens galaxy mass models, thereby offering a powerful tool for distinguishing between Cold Dark Matter and alternative dark matter theories.

Original authors: Kassidy E. Kollmann, James W. Nightingale, Mariangela Lisanti, Andrew Robertson, Oren Slone

Published 2026-02-25

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Kassidy E. Kollmann, James W. Nightingale, Mariangela Lisanti, Andrew Robertson, Oren Slone

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.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

Imagine the universe is a giant, invisible ocean made of Dark Matter. We can't see this ocean directly, but we know it's there because its gravity acts like a giant lens, bending the light from distant galaxies behind it. This phenomenon is called Gravitational Lensing.

Usually, when astronomers look for clumps of this dark matter (called subhalos) inside these giant lenses, they assume the clumps have a standard shape: a "hazy" cloud that gets denser toward the center, but not too dense. This is the standard "NFW" profile, named after the scientists who described it.

The Big Question:
What if the dark matter isn't just a hazy cloud? What if, due to exotic physics, some of these clumps are actually super-dense spikes at their very core? This paper asks: Does the shape of the dark matter clump change how easy it is to find?

The Core Discovery: The "Sharp Needle" vs. The "Soft Pillow"

The authors ran thousands of computer simulations to answer this. They created fake telescope images (like those from Hubble, Euclid, or the James Webb Space Telescope) containing a dark matter clump. They tested three types of clumps:

The Cored Clump (The Soft Pillow): The center is flat and empty.
The NFW Clump (The Standard Haze): The center gets denser, but gradually.
The Steep Clump (The Sharp Needle): The center is incredibly dense and spikes up sharply.

The Result:
They found that the "Sharp Needle" (Steep) clumps are incredibly easy to spot, even if they are very small and light. In fact, they can be detected at masses 10 times smaller than the standard "Haze" clumps.

The Analogy:
Imagine you are trying to find a hidden object in a dark room by looking at how it distorts a sheet of fabric stretched over a frame.

If you put a soft pillow (Cored) or a standard rock (NFW) under the fabric, it makes a gentle, broad dip. If the object is small, the dip is so subtle you might miss it entirely.
If you put a sharp needle (Steep) under the fabric, it creates a tiny, incredibly deep, and sharp spike. Even if the needle is very light, that sharp spike is impossible to ignore. The fabric reacts violently to the sharpness of the tip.

Why This Matters for "Steep" Clumps

The paper highlights a unique superpower of these "Steep" clumps: They are immune to confusion.

Usually, when astronomers try to find these dark matter clumps, they have to subtract the complex shape of the main galaxy (the lens) from the image. If the galaxy isn't a perfect circle or has weird bumps (multipole perturbations), it's hard to tell if a bump in the image is a dark matter clump or just a weird shape of the galaxy itself. It's like trying to hear a whisper in a noisy room; the noise drowns out the whisper.

For the "Soft Pillow" and "Standard Rock": When the galaxy has complex shapes, the "noise" drowns them out. They become invisible.
For the "Sharp Needle": The spike is so sharp and distinct that even a noisy, complex galaxy background can't hide it. The needle's signature is so unique that the galaxy's "noise" can't mimic it.

The Telescope Factor

The authors also tested how different telescopes would see these clumps:

Hubble (Current): Can see the "Sharp Needles" if they are on the ring of light (Einstein ring), but misses the softer ones.
Euclid (Upcoming): Will find even more, but still struggles with the soft ones.
James Webb (JWST): With its incredible sharpness and sensitivity, it can find the "Sharp Needles" even when they are very tiny, and it might even start to see the "Soft Pillows" if they are massive enough.

The Bigger Picture: What is Dark Matter Made Of?

This isn't just about finding clumps; it's about solving the mystery of what dark matter is.

Standard Theory (Cold Dark Matter): Predicts mostly "Standard Rocks" (NFW).
Alternative Theory (Self-Interacting Dark Matter): Predicts that dark matter particles can bump into each other. This can cause the center of a clump to either flatten out (a "Soft Pillow") or collapse into a super-dense spike (a "Sharp Needle").

The Conclusion:
If we look at the sky with our best telescopes and find a lot of these tiny, "Sharp Needle" dark matter clumps, it would be a smoking gun that Self-Interacting Dark Matter is real. If we only find the "Standard Rocks," then the old theory holds up.

Because these "Sharp Needles" are so much easier to find than the others, this paper suggests that we might have been missing a huge population of dark matter subhalos all along, simply because we were looking for the wrong shape. By looking for the "spikes," we might finally crack the code of the universe's invisible ocean.

1. Problem Statement

The Cold Dark Matter (CDM) model predicts a universal Navarro-Frenk-White (NFW) density profile for dark matter subhalos, characterized by an inner density slope of $\beta \approx 1$ . However, alternative dark matter models, such as Self-Interacting Dark Matter (SIDM), predict a diverse range of inner density profiles. SIDM can lead to:

Cored profiles ( $\beta \sim 0$ ): During the core-expansion phase.
Steeply rising (cuspy) profiles ( $\beta > 2$ ): During the gravothermal core-collapse phase.

Current strong gravitational lensing analyses typically assume NFW profiles or treat subhalo concentration as a free parameter but often fix the inner slope to $\beta=1$ . This assumption may bias mass inferences and detection limits. The central problem addressed is: How does the inner density slope ( $\beta$ ) of a subhalo affect its detectability and parameter recovery in strong lensing observations, and how robust are these detections against modeling complexities (e.g., multipole perturbations and source morphology)?

2. Methodology

Mock Data Generation

The authors generated mock strong lensing observations using the open-source software PyAutoLens.

System Configuration: A standard "cross quad" lens system with a source galaxy at $z=1.0$ and a lens galaxy at $z=0.2$ . The lens galaxy is modeled with an Elliptical Power-Law (EPL) mass profile and cored-Sérsic light profiles (bulge + disk).
Subhalo Perturbers: A single subhalo was inserted into the lens galaxy using a generalized NFW (gNFW) profile:
$\rho(r) = \frac{\rho_0}{(r/r_s)^\beta (1 + r/r_s)^{3-\beta}}$
Three distinct inner slopes were tested:
- Cored: $\beta = 0.2$
- NFW: $\beta = 1.0$
- Steep: $\beta = 2.2$
  Subhalo masses ( $M_{200}$ ) ranged from $10^{7.5}$ to $10^{10} M_\odot$ .
Observational Suites: Five suites of mock data were created to test different conditions:
1. HST On-Ring: Subhalo on the Einstein ring; SNR $\sim 30$ .
2. HST Offset: Subhalo $\sim 1''$ off the Einstein ring; SNR $\sim 30$ .
3. HST On-Ring (Higher SNR): SNR $\sim 60$ .
4. Euclid On-Ring: Simulating lower resolution ( $0.1''$ pixel scale); SNR $\sim 30$ .
5. JWST On-Ring: Simulating high resolution ( $0.031''$ ) and high SNR ( $\sim 70$ ).

Modeling and Detection Procedure

Lens Modeling: The authors used a Bayesian non-linear search (Nautilus/Dynesty) to fit the data.
- Smooth Model: Fits the macro-lens (mass and light) without a subhalo.
- Subhalo Model: Adds a gNFW subhalo with free parameters for mass, concentration, inner slope ( $\beta$ ), and position.
Detection Metric: Detection significance was quantified using the difference in log-Bayesian evidence ( $\Delta \ln \mathcal{E}$ ) between the model with a subhalo and the smooth model. A threshold of $\Delta \ln \mathcal{E} \geq 50$ was used to define a "strong detection."
Complexity Tests:
- Multipole Perturbations: To test degeneracy with macro-model complexity, 1st, 3rd, and 4th-order multipoles were added to the lens mass model.
- Pixelized Source Reconstruction: In Section 5, a subset of data was re-analyzed using a flexible, non-parametric (Voronoi mesh) source reconstruction and including the lens galaxy's light in the fit, rather than subtracting it.

3. Key Contributions

Discovery of Slope-Dependent Detectability: The paper demonstrates that the minimum detectable mass of a subhalo is strongly dependent on its inner density slope $\beta$ .
Robustness Against Macro-Model Degeneracy: A unique finding is that subhalos with steep inner profiles ( $\beta \sim 2.2$ ) remain detectable even when complex multipole perturbations are added to the lens model. In contrast, NFW and Cored subhalos become essentially undetectable under the same conditions.
Instrumental Forecasts: The study provides specific detection limits for current (HST) and upcoming (Euclid, JWST) facilities regarding subhalos with diverse density profiles.

4. Key Results

Detectability Limits (HST On-Ring Benchmark)

Steep Subhalos ( $\beta=2.2$ ): Detectable down to $M_{200} \approx 5.4 \times 10^8 M_\odot$ .
NFW Subhalos ( $\beta=1.0$ ): Detectable only down to $M_{200} \approx 6.0 \times 10^9 M_\odot$ (over an order of magnitude higher mass limit).
Cored Subhalos ( $\beta=0.2$ ): Undetectable up to the maximum simulated mass of $10^{10} M_\odot$ .

Impact of Position

When the subhalo is moved off the Einstein ring (Offset suite), the sensitivity to the inner slope diminishes. The deflection strength difference between Steep and NFW profiles decreases as the impact parameter increases.
In the Offset suite, Steep subhalos become significantly harder to detect (limit rises to $\sim 3.2 \times 10^9 M_\odot$ ), and NFW/Cored subhalos remain undetectable.

Impact of Data Quality

Higher SNR (HST): Improves limits slightly (Steep limit $\to 2.3 \times 10^8 M_\odot$ ).
Euclid: Lower resolution reduces sensitivity. Steep subhalos detectable down to $\sim 2.2 \times 10^9 M_\odot$ .
JWST: High resolution and SNR provide the best constraints. Steep subhalos detectable down to $\sim 8.3 \times 10^7 M_\odot$ , and Cored subhalos become detectable down to $\sim 2.7 \times 10^9 M_\odot$ .

Impact of Modeling Complexity

Multipoles: Adding multipoles to the lens model drastically reduces the detectability of NFW and Cored subhalos (pushing them below the detection threshold). However, Steep subhalos on the Einstein ring remain highly detectable (limit only shifts slightly to $8.3 \times 10^8 M_\odot$ ).
Pixelized Sources: Re-analyzing with pixelized source reconstruction and lens light fitting reduces the overall $\Delta \ln \mathcal{E}$ values (making detection harder). However, Steep subhalos remain detectable and their properties ( $\beta$ ) are robustly recovered, whereas NFW subhalos often fall below the detection threshold.

Parameter Recovery

For detectable Steep subhalos, the inner slope $\beta$ is robustly inferred (e.g., $\beta \approx 2.2 \pm 0.1$ for massive systems).
There is a known degeneracy between concentration ( $c$ ) and inner slope ( $\beta$ ), but $\beta$ is generally constrained more tightly than $c$ .
NFW subhalos are often misidentified as having shallower (cored) profiles due to the broadness of the posterior, but they are rarely misidentified as Steep.

5. Significance

Testing Dark Matter Physics: This work highlights that strong lensing is a powerful probe for distinguishing between CDM and SIDM, particularly in the gravothermal core-collapse regime. The detection of low-mass subhalos with steep inner profiles would be strong evidence for SIDM or other non-CDM physics.
Overcoming Systematics: The finding that Steep subhalos break the degeneracy with macro-model complexity (multipoles) is crucial. It suggests that future surveys (Euclid, JWST) can robustly detect these specific subhalo types even when accounting for the complex mass distributions of real elliptical galaxies.
Future Observations: The results suggest that JWST will be capable of detecting subhalos an order of magnitude less massive than current HST limits, specifically for steep profiles. This opens a new window for constraining the microphysics of dark matter through the abundance and properties of subhalos.

In summary, the paper establishes that the inner density slope is a critical factor in subhalo detectability. Subhalos with steep profiles ( $\beta \sim 2.2$ ) are uniquely resilient to modeling degeneracies and can be detected at significantly lower masses than standard NFW subhalos, offering a promising avenue for testing alternative dark matter theories with upcoming astronomical surveys.