Tracing the Evolution of $\Omega_m(z)$ over the Last 10 Billion Years with Non-parametric Methods

This study uses Gaussian Process Regression to non-parametrically reconstruct the redshift evolution of the matter density parameter $\Omega_m(z)$ over the last 10 billion years, finding that while its evolution aligns with the standard $\Lambda$CDM model, the inferred present-day value $\Omega_{m0}$ is highly sensitive to uncertainties in galaxy cluster mass calibration.

The Gist

Imagine the universe as a giant, expanding loaf of raisin bread. As the dough rises (the universe expands), the raisins (galaxies and galaxy clusters) move further apart. But how much "dough" (matter) is actually in that loaf? And does the amount of dough change as the loaf rises, or does it just get spread thinner?

This paper is a team of cosmologists trying to answer that question by looking at the universe's history over the last 10 billion years. They are specifically trying to measure $\Omega_m$, which is essentially the "density of matter" in our cosmic loaf.

Here is the breakdown of their journey, using simple analogies:

1. The Problem: We Can't Weigh the Universe Directly

To know how much matter is in the universe, scientists usually look at Galaxy Clusters. Think of these clusters as the "heavyweights" of the universe—giant groups of galaxies held together by gravity. They are mostly made of invisible "dark matter" (about 80%) and hot gas (about 15%), with just a tiny bit of stars.

The tricky part? We can't put a galaxy cluster on a scale.

• The Scale Issue: To weigh them, scientists have to guess their mass based on how they move or how hot their gas is. But these guesses often have a "systematic error," like a scale that is slightly off because it was calibrated in a different room. This is called Mass Bias.

2. The Solution: A "Model-Free" Detective Approach

Instead of forcing the data to fit a specific theory (like the standard "Big Bang" recipe), the authors used a technique called Gaussian Process Regression (GPR).

• The Analogy: Imagine you are trying to draw a smooth line through a bunch of scattered dots on a piece of paper.
  • Old Way: You assume the dots must form a perfect parabola (a specific curve) and force them to fit that shape, even if they look a bit weird.
  • This Paper's Way: You use a flexible, stretchy rubber band (GPR) to connect the dots. You let the data tell you what the shape is, without forcing it into a pre-made mold. This is called a non-parametric method. It's "model-independent," meaning they aren't assuming the universe works a specific way; they are just letting the observations speak.

3. The Ingredients: Three Clues

To reconstruct the history of matter density, they combined three different types of cosmic clues:

1. Galaxy Cluster Gas: They measured the ratio of hot gas to total mass in two different groups of clusters (one with 44 clusters, one with 103). This is their main "weight" indicator.
2. Cosmic Chronometers: These are like "cosmic stopwatches." By looking at the ages of old galaxies, they can tell how fast the universe was expanding at different times in the past.
3. Type Ia Supernovae: These are "standard candles" (like lightbulbs of known brightness). By seeing how dim they look, scientists can measure how far away they are, which helps map the expansion of the universe.

4. The Results: The Rubber Band Stretches Perfectly

When they ran their flexible rubber band (GPR) through the data, they found something reassuring:

• The Shape: The data followed the standard "recipe" perfectly. As the universe expanded, the matter density dropped exactly as predicted by the standard model ($\Lambda$CDM). It's like watching the raisin bread rise, and the raisins spreading out exactly as physics says they should.
• The "Now" Value: They calculated what the matter density is today ($\Omega_{m0}$).
  • Using the smaller group (44 clusters): ~0.30
  • Using the larger group (103 clusters): ~0.27 (depending on how they corrected the scale).

5. The Catch: The Scale is Still Wobbly

Here is the most important takeaway: The shape of the curve is perfect, but the exact number depends on how you calibrate the scale.

• The Calibration Problem: The authors found that the final number for matter density changed significantly depending on which "calibration" they used for the galaxy clusters.
  • If they used one calibration method (based on weak gravitational lensing), they got a value around 0.27.
  • If they used another (based on the Cosmic Microwave Background), they got a value around 0.21.
• The Metaphor: Imagine you are measuring the height of a growing tree. Your ruler is flexible and accurate (the GPR method), but the zero point on the ruler is slightly different depending on which brand of ruler you use. The growth rate (the evolution over time) is clear, but the exact height today depends on which ruler you trust.

Why Does This Matter?

This study helps solve the "Cosmic Tension" puzzle. Currently, different ways of measuring the universe give slightly different answers (e.g., the "Hubble Tension" or "S8 Tension").

• This paper says: "The universe is expanding and evolving exactly as we thought it would over the last 10 billion years."
• However, it also warns: "We still don't know the exact amount of matter in the universe today because our 'scales' (mass calibrations) aren't perfect yet."

In Summary:
The authors used a flexible, "no-assumptions" math tool to trace the history of matter in the universe. They confirmed that the universe is behaving exactly as the standard model predicts. However, they highlighted that our biggest uncertainty isn't the theory—it's our ability to accurately weigh the galaxy clusters. Until we fix the "scale," the exact number for how much matter is in the universe today will remain a bit of a moving target.

Technical Summary

1. Problem Statement

The paper addresses the challenge of determining the redshift evolution of the matter density parameter, $\Omega_m(z)$, in a way that is largely independent of specific cosmological models.

• Context: Galaxy clusters are powerful probes for cosmology, particularly for constraining $\Omega_m$ and the amplitude of matter fluctuations ($\sigma_8$). However, their use is limited by the difficulty of accurately measuring total cluster masses, leading to systematic uncertainties (e.g., hydrostatic mass bias).
• The Tension: There is an ongoing tension between early-Universe probes (CMB, e.g., Planck) and late-Universe probes (large-scale structure, e.g., weak lensing) regarding the values of $H_0$ and $S_8$.
• Goal: To reconstruct the evolution of $\Omega_m(z)$ over the last ~10 billion years ($z \approx 0$ to $1.2$) using a non-parametric approach. This allows for a direct test of the standard $\Lambda$CDM prediction ($\rho_m \propto (1+z)^3$) without assuming a specific background cosmology or a parametric form for the evolution of matter density.

2. Methodology

The authors employ a Gaussian Process Regression (GPR) technique to reconstruct cosmological quantities directly from observational data, minimizing model dependence.

Theoretical Framework

The core of the analysis relies on the cluster gas mass fraction ($f_{gas}$), defined as the ratio of baryonic gas mass to total cluster mass. The relationship is derived as:
$$\Omega_m(z) = \Omega_{b0} \frac{(1+z)^3 H_0^2}{H^2(z)} \frac{\Upsilon(z) K(z)}{f_{gas}(z)} \left( \frac{d_L^*(z)}{d_L(z)} \right)^{3/2}$$
Where:

• $\Omega_{b0}$: Present-day baryon density (constrained by primordial helium abundance).
• $H(z)$: Hubble parameter derived from Cosmic Chronometers (CC).
• $d_L(z)$: Luminosity distance derived from Type Ia Supernovae (SNe Ia).
• $d_L^*(z)$: Luminosity distance in a fiducial cosmology (used to correct for bias).
• $\Upsilon(z)$: Gas depletion factor (accounts for gas not in the intracluster medium).
• $K(z)$: Calibration term accounting for mass estimation biases (hydrostatic mass bias).

Data Sets

The study combines three distinct data sources:

1. Cosmic Chronometers: 32 measurements of $H(z)$ used to reconstruct the expansion history via GPR.
2. Type Ia Supernovae: The Union3 compilation (2087 SNe Ia, binned into 22 measurements) used to reconstruct the dimensionless luminosity distance $D_L(z)$ via GPR.
3. Galaxy Cluster Gas Mass Fractions ($f_{gas}$): Two independent samples were analyzed:
   • Sample A (44 clusters): From Mantz et al. (2022), covering $0.078 < z < 1.063$. These are massive, relaxed clusters observed by Chandra, measured in a shell ($0.8-1.2 r_{2500c}$) to minimize core physics effects.
   • Sample B (103 clusters): A compilation by Corasaniti et al. covering $0.047 \le z \le 1.235$. These are integrated measurements out to $r_{500c}$, including a larger volume but higher sensitivity to baryonic physics and calibration systematics.

Reconstruction Strategy

• Non-parametric GPR: Used to reconstruct smooth functions for $H(z)$ and $d_L(z)$ from the discrete data points.
• Propagation: These reconstructed functions are substituted into the equation for $\Omega_m(z)$ along with the $f_{gas}$ measurements.
• Systematics Handling: The nuisance parameters $\Upsilon(z)$ and $K(z)$ are modeled with specific priors based on hydrodynamical simulations and weak-lensing calibrations. The study explicitly tests different mass calibration scenarios ($K$) to quantify systematic uncertainties.

3. Key Contributions

• Model-Independent Reconstruction: The paper provides a direct, non-parametric reconstruction of $\Omega_m(z)$ evolution, avoiding the assumption of a specific dark energy equation of state or a fixed $\Lambda$CDM background during the reconstruction phase.
• Quantification of Mass Bias: It rigorously demonstrates that the inferred value of $\Omega_{m0}$ is highly sensitive to the assumed mass calibration ($K$), identifying mass bias as the dominant source of uncertainty, exceeding statistical errors.
• Cross-Validation: By using two distinct cluster samples (shell-based vs. integrated) and three different calibration schemes, the study isolates the impact of astrophysical systematics on cosmological parameters.

4. Results

Evolution of $\Omega_m(z)$

• The reconstructed evolution of $\Omega_m(z)$ for both cluster samples is consistent with the standard $\Lambda$CDM prediction ($\rho_m \propto (1+z)^3$) across the redshift range probed ($z \approx 0$ to $1.2$).
• The data-driven reconstruction shows no statistically significant deviation from the standard scaling law over the last ~10 billion years of cosmic history.

Constraints on $\Omega_{m0}$ (Present-day Matter Density)

The inferred value of $\Omega_{m0}$ varies significantly depending on the cluster sample and the mass calibration assumption:

Sample	Calibration Scheme	$\Omega_{m0}$
44 Clusters	Mantz et al. (2022) Priors	$0.296 \pm 0.044$
103 Clusters	CCCP (Weak Lensing)	$0.271 \pm 0.016$
103 Clusters	CLASH (Weak Lensing)	$0.253 \pm 0.017$
103 Clusters	CMB-based (Planck/SPT)	$0.210 \pm 0.013$

• 44-Cluster Sample: Yields a value ($0.296$) consistent with Planck CMB results ($\sim 0.31$) and large-scale structure probes, though with larger uncertainties.
• 103-Cluster Sample: Shows a strong dependence on calibration.
  • CCCP/CLASH: Values ($0.271, 0.253$) align with "low-$S_8$" preferences of weak lensing surveys (DES, KiDS), potentially alleviating the $S_8$ tension.
  • CMB-based: Yields a very low value ($0.210$). If combined with standard assumptions, this implies a strong suppression of structure growth, potentially overcorrecting the $S_8$ tension and creating tension with CMB data unless other parameters (like $H_0$ or $\sigma_8$) are adjusted.

5. Significance and Conclusion

• Validation of $\Lambda$CDM: The study reinforces the robustness of the standard cosmological model regarding the evolution of matter density, as the non-parametric reconstruction matches the $\rho_m \propto (1+z)^3$ scaling.
• Systematics are Key: The primary takeaway is that astrophysical systematics (mass calibration) are the limiting factor in cluster cosmology, not statistical precision. The choice of calibration ($K$) shifts the central value of $\Omega_{m0}$ by more than $2\sigma$.
• Cosmological Tensions: The results highlight the multidimensional nature of the $H_0$ and $S_8$ tensions.
  • Calibrations favoring lower $\Omega_{m0}$ (like CMB-based) tend to push $H_0$ higher (alleviating the Hubble tension) but may exacerbate the $S_8$ tension.
  • Calibrations favoring intermediate $\Omega_{m0}$ (CCCP/CLASH) align well with the "low-$S_8$" universe preferred by late-time structure probes.
• Future Outlook: The authors conclude that upcoming X-ray and Sunyaev–Zel'dovich surveys, combined with improved weak-lensing mass calibrations, are essential to reduce these systematics. This will allow cluster cosmology to play a more decisive role in resolving current cosmological tensions.

In summary, the paper successfully demonstrates that while the evolution of matter density is robustly consistent with $\Lambda$CDM, the normalization ($\Omega_{m0}$) remains heavily dependent on how cluster masses are calibrated, underscoring the need for independent, high-precision mass measurements to fully leverage galaxy clusters as precision cosmological tools.