🔭 astrophysics

Cluster Infall for Mass Calibration in the Stage-IV Era

This paper presents a velocity distribution model for galaxy cluster outskirts that, when applied to DESI spectroscopic data, enables sub-percent precision mass calibration comparable to Stage IV weak lensing surveys, offering a robust method insensitive to baryonic physics.

Original authors: Connor Sweeney, Eduardo Rozo

Published 2026-03-25

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Connor Sweeney, Eduardo Rozo

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 ocean. In this ocean, galaxy clusters are like massive whirlpools or whirlpools of stars, gas, and dark matter. These whirlpools are the heaviest things in the universe, and understanding how heavy they are helps scientists figure out the rules of the universe itself (like how fast it's expanding).

However, measuring the weight of these cosmic whirlpools is tricky. The "water" right in the center is turbulent, churning with explosions and hot gas (baryonic physics), which makes it hard to get a clean reading.

The Big Idea:
Instead of looking at the messy center, Connor Sweeney and Eduardo Rozo decided to look at the outskirts of these whirlpools. Think of it like trying to weigh a spinning top. If you look at the center, it's a blur. But if you look at the very edge, where the top is just starting to pull in nearby dust, you can see exactly how hard it's pulling.

They call this the "Infall Region." It's the zone where galaxies are just starting to fall into the cluster, like leaves being sucked toward a drain.

How They Did It (The Analogy)

1. The "Traffic Camera" Model
Imagine you are standing on a highway overpass (the "Line of Sight"). You see cars (galaxies) driving toward a massive traffic jam (the cluster).

Some cars are already stuck in the jam (orbiting).
Some are far away, just cruising (expansion).
The Infallers: These are the cars in the middle, speeding up as they get closer to the jam.

The authors built a sophisticated mathematical traffic model. They didn't just count how many cars there were; they looked at the speed and direction of every single car.

Radial Velocity: How fast are they speeding toward the center?
Tangential Velocity: How much are they swerving sideways?

They realized that if you know the speed distribution of these "infalling" cars, you can calculate exactly how heavy the traffic jam (the cluster) must be to pull them in that fast.

2. The "Smoothie" Problem
In the past, scientists tried to describe this traffic with a very complicated recipe that had too many ingredients (parameters). It was like trying to make a smoothie with 50 different fruits; it was hard to know which fruit made the taste what it was.

Sweeney and Rozo found a shortcut. They realized that the "shape" of the speed distribution follows a very specific, predictable pattern (like a bell curve, but with a twist). They simplified the recipe down to just two main ingredients:

The Peak Speed: How fast are the cars moving at their fastest point?
The Spread: How much do the speeds vary?

By focusing on these two things, they could describe the whole traffic pattern with much less math, making their model much more accurate.

3. The "DESI" Telescope
To test this, they used data from a massive digital map of the universe called DESI (Dark Energy Spectroscopic Instrument). Imagine DESI as a giant, high-speed camera that takes a picture of the speed of millions of galaxies at once.

They ran their "traffic model" against the DESI data. The result? Their model could predict the mass of the galaxy clusters with sub-percent precision.

Why This Matters

Think of it this way:

Old Way (Weak Lensing): Trying to weigh a person by looking at how much their shadow stretches on the ground. It works, but it's blurry and affected by the wind (baryonic physics).
New Way (Infall Calibration): Weighing the person by measuring how fast a feather falls toward them. It's cleaner, happens further away from the messy center, and gives a much sharper reading.

The Bottom Line:
This paper proves that by watching how galaxies fall into clusters from a safe distance, we can weigh those clusters with incredible accuracy. This is a game-changer for cosmology because it gives us a new, super-precise ruler to measure the expansion of the universe and the nature of dark energy, without getting confused by the messy physics happening inside the clusters themselves.

In short: They found a way to weigh the universe's heaviest objects by watching the dust dance around the edges, and they did it with a level of precision that rivals the best methods we have today.

1. Problem Statement

Galaxy clusters are powerful cosmological probes, particularly through their abundance as a function of mass. However, accurately calibrating cluster masses is challenging due to nonlinear baryonic physics (e.g., active galactic nucleus feedback) that affects the inner regions of clusters.

The Challenge: Traditional mass proxies (like weak lensing or X-ray emission) often rely on regions sensitive to baryonic effects or require complex systematics modeling.
The Opportunity: The outskirts of galaxy clusters ( $r \gtrsim 5 h^{-1}$ Mpc) are dominated by infalling matter and are relatively insensitive to baryonic physics.
The Gap: While previous studies (e.g., Aung et al. 2023) have used the dispersion of line-of-sight (LOS) velocities in these regions as a standard ruler, there is a lack of models utilizing the full shape of the LOS velocity distribution $P(v_{LOS}|R, M)$ to constrain masses with high precision. Existing models often struggle to parameterize the complex dependence of velocity distributions on both projected radius and mass simultaneously.

2. Methodology

The authors developed a comprehensive statistical framework to model the velocity distribution of infalling galaxies and forecast mass constraints using future spectroscopic data (DESI).

A. Simulation Data

Simulation: Used the Multi-Dark Planck 2 (MDPL2) simulation ( $3840^3$ particles, Planck cosmology).
Galaxy Population: Populated halos using the UniverseMachine algorithm, applying a stellar mass cut ( $M_* > 10^{10} h^{-1} M_\odot$ ) and selecting massive halos ( $M_{vir} \ge 8 \times 10^{13} h^{-1} M_\odot$ ).
Selection: Analyzed galaxy-halo pairs at redshift $z \approx 0.194$ , focusing on projected radii $R \in [5, 30] h^{-1}$ Mpc to exclude orbiting galaxies.

B. Modeling the 3D Velocity Distribution

The authors modeled the joint distribution of radial ( $v_r$ ) and tangential ( $v_t$ ) peculiar velocities, $P(v_r, v_t | r, M)$ , by decomposing it into:
$P(v_r, v_t | r, M) = P(v_r | r, M) \times P(v_t | v_r, r, M)$

Radial Velocity Distribution $P(v_r | r, M)$ :
- Fit: Modeled using the Johnson's SU (JSU) distribution, which captures skewness and kurtosis.
- Parameterization: Reduced the 4 JSU parameters to 2 interpretable parameters: the peak velocity ( $v_{r,peak}$ ) and variance ( $\sigma^2_{vr}$ ).
- Scaling: These parameters were fitted as power laws of radius ( $r$ ) and mass ( $M$ ). An additional mass-dependent correction term was introduced for the shape parameter $\delta$ to improve fit accuracy at high masses.
Conditional Tangential Velocity Distribution $P(v_t | v_r, r, M)$ :
- Fit: Modeled as a Student's t-distribution with fixed degrees of freedom ($dof=5$).
- Variance: The scale parameter $\sigma_{vt}$ was modeled as a cubic function of $v_r$ , with coefficients dependent on radius and mass.

C. Projection to Line-of-Sight (LOS)

To compare with spectroscopic surveys, the 3D model was projected along the LOS:

Weighting: The projection was weighted by the infalling halo-galaxy correlation function $\xi_{inf}(r|M)$ , based on the model from Salazar et al. (2025). This ensures the model accounts for the density profile of infalling galaxies while excluding orbiting populations.
Hubble Flow: The model incorporates the Hubble flow contribution ( $aHr_{LOS}$ ) to convert peculiar velocities to physical LOS velocities ( $v_{LOS}$ ).
Result: An analytical prediction for the observable distribution $P(v_{LOS}|R, M)$ .

D. Forecasting Mass Constraints

Fisher Matrix Analysis: The authors calculated the Fisher information matrix to estimate the uncertainty in recovered cluster masses ( $\sigma_M$ ) based on the number of galaxy pairs expected in DESI surveys (BGS and LRG samples).
Comparison: Results were compared against forecasts for Stage-IV stacked weak lensing (WL) surveys.

3. Key Contributions

Unified Velocity Model: Developed a robust, parameterized model for the joint 3D velocity distribution of infalling galaxies that accurately captures non-Gaussian features (skewness, tails) and their dependence on both radius and halo mass.
Mass-Dependent Calibration: Explicitly parameterized the mass dependence of velocity statistics, moving beyond simple scaling relations to allow for precise mass calibration.
Infall Region Focus: Demonstrated that utilizing the full velocity distribution (not just dispersion) in the cluster outskirts ( $R > 5 h^{-1}$ Mpc) provides a mass calibration method robust to baryonic physics.
DESI Forecasts: Provided the first detailed forecast showing that DESI spectroscopic data can achieve sub-percent level precision in cluster mass calibration.

4. Results

Model Accuracy: The smooth model fits the simulation data with high fidelity, reproducing the peak positions and skewed tails of the velocity distributions. The projected LOS velocity distributions are predicted with $\sim 5\%$ accuracy.
Mass Constraints:
- For the highest redshift bins ( $z \sim 0.8$ ) and high-mass clusters, the method achieves mass uncertainties comparable to stacked weak lensing.
- For lower mass bins and specific redshift ranges, the method outperforms weak lensing, achieving sub-percent precision (better than 1%).
- The precision improves at higher redshifts, contrasting with weak lensing where signal-to-noise typically degrades.
Robustness: By operating at $R \ge 5 h^{-1}$ Mpc, the method avoids the "virialized" and "splashback" regions where baryonic feedback is most significant.

5. Significance

This work establishes cluster infall kinematics as a premier tool for the "Stage-IV" era of cosmology (DESI, Euclid, Rubin).

Systematics Mitigation: It offers a pathway to calibrate cluster masses that is theoretically cleaner than methods relying on the cluster core, potentially reducing systematic errors in cosmological parameter estimation (e.g., $\sigma_8$ , $\Omega_m$ ).
Synergy with Spectroscopy: It leverages the unique capability of spectroscopic surveys (like DESI) to measure individual galaxy velocities, transforming them from mere tracers into dynamical mass probes.
Future Outlook: The authors conclude that with perfect model calibration, this approach could provide mass constraints superior to weak lensing, making it a critical component for future precision cosmology analyses.

Limitations noted: The current forecasts assume perfect model calibration and neglect cosmology dependence. Future work must address potential biases between infalling galaxies and the underlying dark matter infall.