🔭 astrophysics

Spatially Resolved Kinematics of SLACS Lens Galaxies. II: Breaking Degeneracies with Lensing and Dynamical Models

By combining spatially resolved kinematics from Keck KCWI with Hubble Space Telescope lensing data for 14 SLACS galaxies, this study confirms that nearly isothermal power-law mass profiles are a robust first-order description of these systems that do not introduce measurable bias in time-delay cosmography, while demonstrating that more flexible models are required to fully capture the detailed kinematic data.

Original authors: Shawn Knabel, Tommaso Treu, Michele Cappellari, Simon Birrer, Xiang-Yu Huang, Anowar J. Shajib, William Sheu

Published 2026-04-15

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Shawn Knabel, Tommaso Treu, Michele Cappellari, Simon Birrer, Xiang-Yu Huang, Anowar J. Shajib, William Sheu

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

The Big Picture: Solving the "Hubble Tension"

Imagine the universe is a giant clock, and astronomers are trying to figure out exactly how fast it is ticking. This speed is called the Hubble Constant ( $H_0$ ).

Right now, there is a massive argument in the scientific community. One group of astronomers (looking at the early universe) says the clock is ticking at one speed. Another group (looking at nearby stars and galaxies) says it's ticking faster. They disagree so much that it's like two people arguing over whether a car is going 60 mph or 70 mph, and the difference is so huge it suggests the laws of physics might be broken.

To settle this, scientists are using gravitational lenses. These are massive galaxies that act like a cosmic magnifying glass, bending light from a distant quasar behind them. By measuring how long it takes for that light to travel different paths around the galaxy, they can calculate the speed of the universe's expansion.

The Problem: To get an accurate time, you need to know exactly how heavy the "magnifying glass" (the lens galaxy) is. But there's a trick in the math called the Mass-Sheet Degeneracy (MSD).

The Analogy: The Foggy Window

Imagine you are looking at a friend through a window.

The Lensing Effect: You can see your friend clearly, but the window distorts their shape slightly.
The Mass-Sheet Degeneracy: Now, imagine someone puts a very thin, slightly foggy sheet of plastic over the window.
- Your friend still looks like your friend (the shape doesn't change).
- The distance to your friend seems the same.
- But, the fog makes everything look slightly dimmer and changes the perceived scale.

In cosmology, this "fog" is an invisible sheet of mass along our line of sight. Because the math allows for this "fog" to exist without changing the shape of the images, scientists can't tell if the galaxy is heavy and close, or light and far away. This uncertainty throws off the calculation of the universe's expansion speed.

The Paper's Solution: The "Two-Pronged" Detective Work

The authors of this paper (Knabel et al.) decided to solve this mystery by using two different detective tools at the same time: Lensing and Dynamics.

Tool 1: The Lensing Map (The Shape)

They used the Hubble Space Telescope to take high-resolution pictures of the lensed galaxies. This tells them the shape of the gravity well (how the light bends).

Analogy: This is like looking at the shadow a tree casts. You know the tree is there, and you know the general shape of the shadow.

Tool 2: The Kinematic Map (The Dance)

They used the Keck KCWI telescope to watch the stars inside the galaxy move. They didn't just look at the average speed; they made a 2D map of how every single star was dancing (moving left, right, fast, slow).

Analogy: This is like watching the leaves on the tree blowing in the wind. If the tree is heavy, the leaves move slowly. If the tree is light, the leaves whip around fast. By watching the "dance" of the stars, they can weigh the galaxy directly.

The Secret Weapon: Breaking the "Fog"

When they combined the Shadow (Lensing) and the Dance (Kinematics), something magical happened.

The "fog" (Mass-Sheet Degeneracy) that confused the Lensing tool doesn't affect the Dance. The stars move based on the actual gravity, not the mathematical trickery of the lensing equations. By forcing the model to satisfy both the Shadow and the Dance simultaneously, the "fog" disappears. They can finally weigh the galaxy accurately.

The "Anisotropy" Puzzle: The Direction of the Dance

There was one catch. To weigh the galaxy using the stars' dance, you have to guess how the stars are moving. Are they moving in perfect circles? Or are they zooming in and out like a swarm of bees? This is called Anisotropy.

The Problem: If you guess wrong about the direction of the dance, you get the weight wrong. This is called the Mass-Anisotropy Degeneracy (MAD).
The Fix: The authors used data from nearby galaxies (where we have super-clear data) to create a "rulebook." They found that most of these heavy, round galaxies are "isotropic"—meaning the stars dance in all directions equally, like a swarm of bees in a jar. They used this rulebook as a "prior" (a helpful hint) to guide their calculations for the distant galaxies.

The Results: What Did They Find?

The "Fog" is Gone: They successfully measured the "Mass-Sheet" parameter ( $\lambda_{int}$ ) for 14 galaxies with incredible precision (about 5.8% error).
The Average is Perfect: On average, these galaxies act exactly like a simple, smooth "power-law" (a specific mathematical curve). The average result was 1.01, which is basically 1.0.
- Translation: The "fog" isn't actually there on average. The simple math models used in the past were actually doing a pretty good job!
The Scatter: However, when they looked at individual galaxies, some were a bit "foggier" than others. The results varied slightly from galaxy to galaxy.
- Translation: While the simple model works for the average galaxy, if you want to be super precise for one specific galaxy, you need a more complex model.

The Conclusion: Why This Matters

This paper is a huge step forward for cosmology.

Validation: It confirms that the "Time-Delay Cosmography" method (using lensing to measure the universe's speed) is robust. The "fog" doesn't introduce a systematic error that would ruin the measurement of the Hubble Constant.
Refinement: It shows that while simple models work for the big picture, we need to be careful with individual galaxies.
Future: It paves the way for using these galaxies to solve the "Hubble Tension." If we can weigh these galaxies perfectly, we can finally tell if the universe is expanding faster than we thought, or if our understanding of physics needs a rewrite.

In a nutshell: The authors built a super-accurate scale by combining a shadow and a dance. They proved that the "fog" in our cosmic measurements is mostly an illusion, and that the universe's expansion rate can be measured with high confidence, bringing us closer to solving one of the biggest mysteries in physics.

1. Problem Statement

The measurement of the Hubble constant ( $H_0$ ) via time-delay cosmography of lensed quasars is currently limited by systematic uncertainties, most notably the Mass-Sheet Degeneracy (MSD). The MSD allows for a transformation of the lens mass profile (adding a uniform mass sheet) that preserves image positions and shapes but scales the inferred time delays and $H_0$ .

To break the MSD, joint modeling of gravitational lensing and stellar dynamics is required. However, stellar dynamics introduces its own fundamental limitation: the Mass-Anisotropy Degeneracy (MAD). The second moment of stellar velocities (velocity dispersion) cannot uniquely determine the mass profile without assumptions about the orbital anisotropy of the stars. While previous studies (e.g., TDCOSMO 2025, "TDC-25") utilized spherical dynamical models and informative priors from local galaxies to address these issues, they relied on azimuthally averaged kinematics. This paper aims to break these degeneracies using spatially resolved 2D kinematic maps and axisymmetric Jeans modeling to test the validity of spherical approximations and quantify intrinsic scatter in mass profiles.

2. Methodology

The authors analyzed a sample of 14 strong gravitational lens galaxies from the Sloan Lens ACS (SLACS) survey, utilizing high-quality integral-field spectroscopy (IFS) data from the Keck Cosmic Webb Imager (KCWI).

Data Processing: Kinematic maps ( $V_{rms}$ ) were extracted using the Penalized Pixel-Fitting (pPXF) method. A full covariance matrix was constructed to account for spatial correlations and systematic uncertainties (estimated at ~1.4% uncorrelated systematics).
Modeling Framework: The study employed Jeans Anisotropic Modeling (JAM) via the jampy software.
- Mass Profile: Modeled as a Power-Law Elliptical Mass Density (PEMD) profile ( $\rho \propto r^{-\gamma}$ ) subject to the Mass-Sheet Transform (MST).
- Degeneracy Breaking:
  - MSD: Broken by combining lensing constraints (HST imaging) with dynamical constraints. The internal MST parameter ( $\lambda_{int}$ ) was fitted freely.
  - MAD: Broken by applying informative priors on orbital anisotropy ( $\beta_{ani}$ ) and intrinsic shape ( $q_{intr}$ ), derived from high-resolution local galaxy studies (e.g., SAURON, ATLAS $^{3D}$ ).
Model Variations: The authors compared six distinct model configurations for each galaxy:
1. Spherical vs. Axisymmetric geometry.
2. Spherically-aligned vs. Cylindrically-aligned velocity ellipsoids.
3. Dynamics-only vs. Joint Lensing/Dynamical constraints.
Statistical Analysis: Bayesian inference was performed using MCMC (emcee). Model selection was conducted using the Bayesian Information Criterion (BIC) to weight the contributions of different model configurations to the final results.

3. Key Contributions

First Application of 2D Axisymmetric JAM to SLACS: This work moves beyond the radially binned, spherical approximations used in TDC-25, utilizing full 2D kinematic maps to test the robustness of the spherical assumption for massive early-type galaxies (ETGs).
Quantification of Intrinsic Scatter: The study explicitly measures the intrinsic scatter in the mass-sheet parameter ( $\lambda_{int}$ ) and the power-law slope ( $\gamma$ ) across the sample, moving beyond mean values to assess galaxy-to-galaxy variations.
Validation of Priors: It rigorously tests the assumption that local galaxy anisotropy priors are valid for lens galaxies at $z \sim 0.2-0.3$ , confirming that massive ETGs are effectively isotropic within one effective radius.
Systematic Error Budget: The paper provides a detailed analysis of how systematics (e.g., background source contamination, low signal-to-noise outer bins) affect the inference of $\lambda_{int}$ , demonstrating that the data quality is sufficient to constrain deviations from pure power laws.

4. Key Results

Mass Density Slope: The sample exhibits nearly isothermal total mass density profiles with a mean slope of $\langle \gamma \rangle = 2.04 \pm 0.02$ and an intrinsic scatter of $0.08^{+0.03}_{-0.02}$ . This is consistent with previous lensing-only studies.
Mass-Sheet Parameter: The mean internal mass-sheet parameter is $\langle \lambda_{int} \rangle = 1.01 \pm 0.03$ , with an intrinsic scatter of $0.11 \pm 0.03$ .
- The mean value is consistent with unity (no MSD) within $1\sigma$ , suggesting that on average, power-law profiles do not bias time-delay cosmography.
- However, the significant intrinsic scatter ( $\sim 11\%$ ) indicates that individual galaxies deviate from the pure power-law form, necessitating flexible models for precision cosmology.
Comparison with TDC-25: The results for $\lambda_{int}$ are consistent with the TDC-25 analysis (which used spherical models and azimuthally averaged data) to within $1\sigma$ . This confirms that the scatter reported in TDC-25 is real and not an artifact of the spherical approximation.
Spherical vs. Axisymmetric: Axisymmetric models (both spherical and cylindrical alignment) yield results consistent with spherical models within uncertainties. No significant bias was detected in $\lambda_{int}$ or $\gamma$ due to the geometric assumption, though axisymmetric models showed a slight preference for cylindrical alignment in fast-rotator candidates.
Dynamics-Only Constraints: Even without lensing priors, the spatially resolved kinematics alone constrain $\lambda_{int}$ to an average precision of 8.2%, improving to 5.8% when joint lensing constraints are applied. This demonstrates the power of 2D kinematics to detect deviations from power-law profiles.

5. Significance

Cosmological Impact: The finding that the mean $\lambda_{int} \approx 1$ with no measurable bias validates the use of the SLACS sample as an external dataset for time-delay cosmography. However, the intrinsic scatter highlights a limit to the precision achievable with simple power-law models; future cosmography will require more flexible mass models to reduce this scatter.
Methodological Validation: The study confirms that the strategy of using local galaxy priors to break the MAD is robust for lens galaxies up to $z \sim 1$ . It also validates the use of spherical Jeans models as a first-order approximation, provided that correction factors (as proposed by Huang et al. 2025) are applied when high-resolution 2D data is unavailable.
Future Directions: The paper argues that to fully exploit strong lenses for precision cosmology, future work must decouple stellar mass from dark matter using flexible mass-to-light ratios and dark matter profiles. Furthermore, it emphasizes the need for higher-resolution data (e.g., from the ELT or JWST) to constrain Schwarzschild orbit models directly, rather than relying on priors.

In conclusion, this paper successfully breaks the MSD and MAD for a sample of 14 lenses using spatially resolved kinematics, confirming that while power-law profiles are a good first-order description, the intrinsic scatter in individual galaxy mass profiles must be accounted for in high-precision cosmological measurements.