🔭 astrophysics

A Salpeter IMF and an NFW halo: Disentangling the dark and stellar mass of an elliptical galaxy through precise lens modelling of a double-source-plane system

By applying precise strong lensing modeling to the double-source-plane system J0946+1006, the study disentangles stellar and dark matter components to confirm a Salpeter initial mass function and an NFW dark matter halo, yielding a stellar mass of $4.4 \times 10^{11} M_{\odot}$ and demonstrating a robust methodology for future dark matter studies with Euclid.

Original authors: Tian Li, Thomas E. Collett, Coleman M. Krawczyk, Giovanni Granata, Wolfgang J. R. Enzi, Daniel J. Ballard, Natalie E. P. Lines, Ana Sainz de Murieta, Luke Weisenbach, Dan Ryczanowski

Published 2026-07-10

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Tian Li, Thomas E. Collett, Coleman M. Krawczyk, Giovanni Granata, Wolfgang J. R. Enzi, Daniel J. Ballard, Natalie E. P. Lines, Ana Sainz de Murieta, Luke Weisenbach, Dan Ryczanowski

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 as a giant, cosmic funhouse mirror. Usually, when a massive galaxy sits in front of a distant one, its gravity acts like a lens, bending the light of the background object into a perfect ring or a set of distorted arcs. This is called "strong gravitational lensing." Scientists use these rings to weigh the galaxy, but there's a tricky problem: it's like trying to figure out how much of a cake's weight comes from the fluffy sponge and how much comes from the heavy chocolate frosting just by looking at the whole thing. You can't easily tell the "dark matter" (the invisible frosting) from the "stars" (the visible sponge).

Enter the "Jackpot" lens, a rare cosmic treasure found in the galaxy SDSSJ0946+1006. This isn't just a one-ring show; it's a double-decker. There are two background galaxies at different distances being warped by the same foreground galaxy. This setup is the universe's way of handing scientists a second clue, helping them untangle the mystery of what the galaxy is actually made of.

The Great Cosmic Weigh-In

The team of astronomers, led by Tian Li, decided to build a super-detailed digital model of this galaxy to see if they could finally separate the stars from the dark matter. They treated the galaxy like a complex recipe:

The Stars: They modeled the visible light using a "multi-Gaussian expansion," which is basically stacking 15 different fuzzy, elliptical blobs of light on top of each other to match the galaxy's shape perfectly.
The Dark Matter: They assumed the invisible dark matter halo followed a specific shape called a "generalised NFW" profile. Think of this as a cloud of invisible particles that is denser in the middle and gets thinner as you go out, but they allowed the "steepness" of that density to change.

They ran this model through a massive computer simulation using a tool called Herculens, which is built to run super-fast on graphics cards (the same kind gamers use). They didn't just guess; they let the data speak, adjusting the model millions of times to see what fit best.

What They Found (The Real Deal)

After crunching the numbers, the data told a very specific story:

The Stars are "Normal": The team found that the stars in this galaxy have a very consistent weight-to-light ratio. In plain English, the stars aren't hiding any secret heavyweights. The data prefers a "Salpeter-like" Initial Mass Function (IMF). Imagine the IMF as a rulebook for how many big, heavy stars versus small, light stars a galaxy makes. A Salpeter rulebook means the galaxy is full of standard, heavy stars, not a weird mix of tiny, light ones.
- The Numbers: They calculated the total stellar mass to be 4.4 +0.25 −0.39 × 10¹¹ M⊙ (that's 440 billion times the mass of our Sun, give or take a bit).
The Dark Matter is "Standard": The invisible dark matter halo turned out to be very close to the standard "NFW" shape predicted by the Cold Dark Matter theory. It's not weirdly steep or flat; it's just right.
- The Numbers: The inner slope of this dark matter cloud (how fast it gets dense near the center) is 1.04 +0.10 −0.14. This is very close to 1, which is the classic NFW prediction.
- The total mass of the dark matter halo is 1.11 +0.37 −0.32 × 10¹³ M⊙.
The Ratio: When they compared the star mass to the halo mass, they found a ratio of log10(M200/M★) = 1.41 +0.13 −0.14. This means the dark matter halo is about 26 times heavier than the stars inside it. This is slightly higher than what is typically seen in other galaxies at this distance, but it's still within the normal range of variation (about 0.1 dex higher than average).

What They Ruled Out (The "Not This" List)

The paper is very clear about what this galaxy is NOT:

It's not a "Steep Cusp": Previous studies on this same galaxy suggested the dark matter was incredibly steep in the center (a slope of about 1.7). This new, more flexible model rules that out. The data does not support a super-steep dark matter spike. The apparent steepness in older studies was likely an illusion caused by using a simpler, rigid model (a single power law) that couldn't bend enough to fit the real, slightly different shape of the dark matter.
It's not a "Bottom-Heavy" Star Factory: The data argues against the idea that the center of this galaxy is packed with tiny, low-mass stars that make the galaxy heavier than it looks (a "Chabrier" IMF). If that were true, the stars would be much lighter, and the math wouldn't add up with the lensing data. The stars are heavy, just like the Salpeter rulebook says.
It's not a "Flat" Mass-to-Light Gradient: The team allowed the possibility that the stars get heavier or lighter as you move away from the center (a gradient). However, the data prefers a flat profile. The stars weigh the same everywhere relative to how bright they are.

How Sure Are They?

The authors are quite confident, but they are careful not to overstate it.

The Method: They didn't just simulate this; they measured it using real telescope data from the Hubble Space Telescope (HST) and the Very Large Telescope (VLT).
The "Double-Source" Trick: The biggest strength of their result is the "double-source-plane" geometry. Because there are two background galaxies, the "Mass-Sheet Transformation" (a sneaky mathematical trick that usually lets scientists hide the true mass profile) is broken. The model has to fit both rings simultaneously, which locks the answer down much tighter than usual.
The Confidence: They state that even when they added extra data (like the speed of the stars moving inside the galaxy) to their model, the results barely changed. This suggests the lensing data alone was already doing a fantastic job. The inner slope of the dark matter is measured as 1.04 +0.10 −0.14, which is a solid measurement, though the uncertainty range (the plus/minus numbers) shows there is still a little wiggle room.

The Big Picture

This paper is like finding a perfect puzzle piece that finally fits the picture of how massive galaxies are built. It confirms that for this "Jackpot" galaxy, the stars are normal, the dark matter is standard, and the whole thing follows the rules of the Cold Dark Matter theory.

The authors suggest that this method—using double-source lenses—is the future. With upcoming telescopes like Euclid and the Nancy Grace Roman Space Telescope, we might find thousands of these "Jackpot" systems. That will allow scientists to stop guessing and start mapping the dark matter of the universe with incredible precision, finally solving the mystery of how much "invisible frosting" is on our cosmic cakes.

Technical Summary: Disentangling Dark and Stellar Mass in the Jackpot Lens

Problem Statement
Determining the distribution of dark matter (DM) in massive galaxies is challenging because DM does not emit light, requiring inference from its gravitational influence. While strong gravitational lensing provides a direct probe of total mass, it suffers from the mass-sheet transformation (MST) degeneracy, where different mass profiles can produce identical lensing observables. Furthermore, disentangling the stellar and dark matter components is complicated by uncertainties in the stellar mass-to-light ratio ( $M/L$ ), which depends on the Initial Mass Function (IMF) and potential radial gradients. Previous studies of the "Jackpot" lens (J0946+1006), a double-source-plane lens (DSPL), have yielded conflicting results regarding the inner DM slope ( $\gamma_{in}$ ), with some suggesting very steep cusps ( $\gamma \approx 1.7$ ) and others finding profiles consistent with the Navarro-Frenk-White (NFW) prediction. These discrepancies often arise from simplified modeling assumptions, such as constant stellar $M/L$ ratios, spherical symmetry, or the use of single power-law profiles that may not capture the true complexity of composite mass distributions.

Methodology
The authors present a strong lensing analysis of J0946+1006 using a composite mass model that explicitly separates stellar and dark matter components, leveraging the unique geometry of the DSPL system to break the MST degeneracy.

Data: The analysis utilizes Hubble Space Telescope (HST) ACS F814W imaging and Multi Unit Spectroscopic Explorer (MUSE) integral-field spectroscopy for stellar velocity dispersion measurements. The system comprises a foreground lens at $z=0.222$ and two background sources at $z=0.609$ (S1) and $z=2.045$ (S2).
Light Modeling: The lens light is modeled using a Multi-Gaussian Expansion (MGE) with 15 components. Source light is modeled using both parametric (elliptical Gaussians) and non-parametric (pixelated Gaussian Processes with Matérn covariance) approaches.
Mass Modeling:
- Stellar Component: The stellar mass is tied to the light model via a global mass-to-light ratio ( $M/L$ ) and allows for a free radial gradient ( $\nabla M/L$ ).
- Dark Matter Component: The primary model uses an elliptical generalised NFW (gNFW) profile, approximated by a 3D MGE of 30 Gaussians to accelerate computation. The inner slope ( $\gamma_{in}$ ) and scale radius ( $r_s$ ) are free parameters.
- Alternative Model: For comparison, a model using an Elliptical Power Law (EPL) for the dark matter component (Star+EPL) is also fitted.
Inference Strategy: The analysis employs the Herculens code, built on JAX for GPU acceleration and automatic differentiation. The workflow involves Stochastic Variational Inference (SVI) to find a posterior mode, followed by Hamiltonian Monte Carlo (NUTS) within a Gibbs sampling scheme to map the full posterior distribution.
Breaking Degeneracies: The DSPL geometry is used to constrain the MST, as a mass sheet that fits the first source plane generally fails to reproduce the imaging of the second source. Additionally, the authors incorporate kinematic data (line-of-sight velocity dispersion $\sigma_{los} = 249.1 \pm 2.1$ km s $^{-1}$ ) and a concentration-mass ( $c-M$ ) prior from $\Lambda$ CDM simulations to reweight the posteriors.

Key Contributions

Composite Modeling with DSPL Constraints: The paper applies a flexible, two-component (star + gNFW) mass model to a DSPL system, utilizing the second source plane to partially break the MST without relying solely on kinematic assumptions.
Flexible Stellar $M/L$ : Unlike previous analyses of this system that assumed a constant stellar $M/L$ , this work allows for a radial gradient and a free global normalization, testing the robustness of IMF inferences against these variations.
Re-evaluation of Inner Slope: The study challenges previous claims of a very steep inner DM cusp ( $\gamma \approx 1.7$ ) by demonstrating that such steepness may be an artifact of power-law modeling assumptions rather than a physical reality.

Results

Stellar Mass and IMF: The data prefer an approximately constant stellar mass-to-light ratio ( $\nabla M/L \approx 0$ ). The inferred stellar mass is $M_\star = 4.4^{+0.25}_{-0.39} \times 10^{11} M_\odot$ . The IMF mismatch parameter is $\alpha_{SPS} = 0.96^{+0.16}_{-0.14}$ (with kinematic and $c-M$ priors), indicating consistency with a Salpeter IMF and tension with a Chabrier IMF at the $\sim 3\sigma$ level.
Dark Matter Profile: The inner slope of the dark matter halo is constrained to $\gamma_{in} = 1.04^{+0.10}_{-0.14}$ , consistent with the standard NFW prediction ( $\gamma=1$ ). This contrasts with previous findings of $\gamma \approx 1.7$ ; the authors attribute the earlier steep results to the limitations of EPL models, which cannot access the low-slope parameter space allowed by gNFW profiles.
Halo Mass and SHMR: The total halo mass is $M_{200} = 1.11^{+0.37}_{-0.32} \times 10^{13} M_\odot$ , yielding a stellar-to-halo mass ratio of $\log_{10}(M_{200}/M_\star) = 1.41^{+0.13}_{-0.14}$ . This value lies $\sim 0.1$ dex above typical literature relations at similar redshift, a difference comparable to the intrinsic scatter of those relations.
Model Comparison: When an EPL profile is forced to model the dark matter (Star+EPL), the fit recovers a steep slope ( $\gamma \approx 1.61$ ), reproducing previous results. However, the gNFW model shows that this steepness is an approximation artifact over the narrow radial range probed by the arcs.

Significance and Claims
The paper claims that the DSPL geometry provides constraining power on the mass profile comparable to that of stellar kinematics, offering a practical alternative for high-redshift lenses where kinematic data is difficult to obtain. The results suggest that the inner dark matter profiles of massive early-type galaxies are consistent with NFW predictions when modeled with sufficient flexibility and when the MST is constrained by multiple source planes. The authors conclude that the Jackpot lens does not require an intrinsically steep dark matter cusp, and that the apparent steepness in previous studies was likely a consequence of modeling choices. They posit that this approach serves as a template for future large-scale DSPL analyses with upcoming surveys like Euclid, the Chinese Space Station Telescope (CSST), and the Nancy Grace Roman Space Telescope, which will enable population-level studies to statistically recover intrinsic dark matter halo properties and resolve selection biases.