$X+y$: insights on gas thermodynamics from the combination of X-ray and thermal Sunyaev-Zel'dovich data cross-correlated with cosmic shear

By cross-correlating cosmic shear, thermal Sunyaev-Zel'dovich, and X-ray data, this study develops a physical model that jointly constrains the spatial distribution and thermodynamic properties of hot gas while highlighting the necessity of accounting for X-ray AGN contamination and non-thermal pressure support to resolve parameter tensions.

The Gist

Imagine the universe as a giant, invisible ocean. We can see the islands (galaxies) and the ships sailing on them, but the water itself—the hot gas that makes up most of the "normal" matter in the universe—is invisible to our eyes. Scientists have been trying to figure out how this invisible ocean behaves: how deep it is, how hot it gets, and how much of it gets pushed out of the way by the islands.

This paper is like a detective story where three different types of detectives team up to solve the mystery of this cosmic ocean.

The Three Detectives

1. The Lens Detective (Cosmic Shear): Imagine looking at a distant galaxy through a funhouse mirror. The gravity of the invisible gas and dark matter in between bends the light, distorting the galaxy's shape. This detective (using data from the Dark Energy Survey) maps out where the total mass is, but it can't tell the difference between the heavy dark matter and the lighter hot gas.
2. The Heat Detective (tSZ): Hot gas doesn't just sit there; it's boiling with energy. When light from the Big Bang (the Cosmic Microwave Background) passes through this hot gas, it gets a little "kick," changing its color slightly. This detective (using data from the Planck satellite) sees the heat of the gas.
3. The Glow Detective (X-ray): When gas is super hot and dense, it glows in X-rays, like a hot stove element. This detective (using data from the ROSAT satellite) sees the brightness of the gas.

The Problem: The "Blind Spot"

Individually, these detectives have blind spots:

• The Heat Detective knows the gas is hot, but doesn't know exactly how much of it is there. It's like feeling a warm breeze but not knowing if it's a small fan or a giant heater.
• The Glow Detective sees the brightness, which depends on how dense the gas is. But if the gas is hot, it glows differently than if it's cool.
• The Lens Detective sees the total weight but can't separate the gas from the invisible dark matter.

If you only use one detective, you get stuck in a loop of guesses. You can't tell if a signal is strong because there's a lot of gas, or because the gas is very hot.

The Solution: The Team-Up

The authors of this paper decided to make these three detectives work together. They looked at how the "Lensing" distortions matched up with the "Heat" and the "Glow" across the sky.

By combining the data, they could finally break the code. It's like having a recipe that tells you both the temperature of the oven and the weight of the cake batter. Now you can figure out exactly how much flour (gas) and how much heat you have.

What They Found

Using this team approach, they built a simple model to describe the cosmic gas and found some key numbers:

• The "Halfway" Mass: They found the size of the galaxy cluster where half of the gas has been kicked out. It's like finding the exact size of a pot where the boiling water starts splashing out the sides.
• The Gas "Stiffness": They measured how the gas pressure changes as you move away from the center of a galaxy cluster.
• The Temperature: They figured out how hot the gas is in the center compared to what physics predicts it should be.

The Plot Twist: The "Imposter"

Here is where it gets tricky. The detectives realized there was a third, sneaky character in the room: Active Galactic Nuclei (AGN). These are super-bright black holes at the centers of galaxies that shoot out huge amounts of energy.

• The Confusion: The "Glow Detective" (X-rays) was seeing light from these black holes and thinking, "Wow, there must be a huge amount of hot gas here!" But actually, it was just the black hole screaming.
• The Fix: The team had to teach their model to ignore the black holes. Once they filtered out the "imposter" signals, the measurements of the gas became much clearer.

They also realized that the gas isn't just hot; it has some "non-thermal" energy (like turbulence or magnetic fields) pushing on it, which they had to account for to get the math right.

Why This Matters

Why do we care about invisible gas?

1. Precision Cosmology: To understand the universe's expansion and the "Dark Energy" pushing it apart, we need to know exactly how much normal matter is hiding in the gas. If we get the gas wrong, our calculations for the whole universe are wrong.
2. The "S8 Tension": There is a famous disagreement in physics. Some measurements say the universe is "clumpier" than others. This paper suggests that gas feedback (gas being pushed out by black holes) might be the reason for this disagreement. If gas is pushed out of galaxy clusters, it smooths things out, changing how we see the universe's structure.

The Bottom Line

This paper is a success story of cross-referencing. By combining three different ways of looking at the sky, and by being careful to filter out the "noise" from black holes, the team created a much clearer picture of the invisible hot gas that fills our universe. They proved that a relatively simple model can explain complex data, provided you account for the messy reality of black holes and turbulent gas.

It's like finally solving a jigsaw puzzle where three people were holding different pieces, and once they put them together, the picture of the cosmic ocean finally made sense.

Technical Summary

Here is a detailed technical summary of the paper "X + y: insights on gas thermodynamics from the combination of X-ray and thermal Sunyaev-Zel'dovich data cross-correlated with cosmic shear."

1. Problem Statement

The distribution and thermodynamic properties of cosmic baryons (specifically hot intergalactic gas) remain a major uncertainty in precision cosmology. While baryons constitute ~5% of the universe's energy density, their complex physics—governed by radiative cooling, gravitational virialization, and non-gravitational heating from stars and Active Galactic Nuclei (AGN)—creates significant uncertainties in:

• Weak Lensing: Baryonic feedback suppresses the matter power spectrum on small scales, potentially mimicking or masking cosmological parameters (e.g., contributing to the $S_8$ tension).
• CMB Secondary Anisotropies: Uncertainties in cluster mass-observable relations and optical depth degeneracies.
• X-ray Science: Hydrostatic mass biases in cluster measurements.

Current probes often suffer from degeneracies. For instance, the thermal Sunyaev-Zel'dovich (tSZ) effect depends on the line-of-sight integral of thermal pressure ($P_{th} \propto \rho T$), creating a degeneracy between gas density ($\rho$) and temperature ($T$). Similarly, X-ray emissivity depends on $\rho^2 \Lambda(T)$, offering complementary but distinct sensitivities. The challenge is to develop a unified physical model that can simultaneously describe X-ray and tSZ data while constraining the gas thermodynamics and accounting for observational contaminants like unresolved AGN and non-thermal pressure.

2. Methodology

Data Sources:
The authors perform a joint analysis of cross-correlations between:

1. Cosmic Shear: From the Dark Energy Survey (DES) Year 3 release (4 redshift bins).
2. X-ray Maps: From the ROSAT All-Sky Survey (RASS) in the 0.5–2.0 keV band.
3. tSZ Maps: Compton-$y$ parameter maps from the Planck 2015 release (MILCA component separation).

Theoretical Framework:
The study employs a hydrodynamical halo model to describe the gas distribution within dark matter halos. The model decomposes the total matter density into Cold Dark Matter (CDM), stars, bound gas, and ejected gas.

• Bound Gas: Modeled with a Komatsu-Seljak (KS) profile derived from hydrostatic equilibrium with a polytropic equation of state ($P \propto \rho^\Gamma$).
• Ejected Gas: Modeled with a Gaussian profile representing gas expelled by AGN feedback.
• Thermodynamics: The gas temperature profile is parametrized by a central temperature ($T_c$) scaled by a factor $\alpha_T$ relative to the virial temperature, and a polytropic index $\Gamma$. Non-thermal pressure support is introduced via a parameter $\gamma_T$.

Key Parameters Constrained:

• $\log_{10} M_c$: The mass scale at which half of the gas content has been expelled from the halo (halfway mass).
• $\Gamma$: The polytropic index describing the scale dependence of the gas distribution.
• $\alpha_T$: The ratio of central gas temperature to virial temperature.

Analysis Strategy:
The authors compute angular power spectra ($C_\ell$) for the cross-correlations ($\gamma \times X$, $\gamma \times y$) and the auto-correlation ($X \times y$). They use a Gaussian likelihood to constrain the model parameters. Crucially, they test the robustness of the model against:

• AGN Contamination: Masking resolved point sources and marginalizing over the contribution of unresolved AGN using a luminosity function model.
• Non-thermal Pressure: Allowing for non-thermal pressure support in the gas.
• Systematics: Testing different tSZ maps (Planck PR4, NILC, CIB-deprojected) and scale cuts to assess foreground contamination.

3. Key Contributions

1. First Joint Constraints: This work presents the first simultaneous constraints on gas thermodynamics derived from the cross-correlation of cosmic shear with both X-ray and tSZ data within a single physical model.
2. Resolution of Parameter Tension: The authors identify a significant tension ($\sim 3.5\sigma$) between the polytropic index $\Gamma$ derived from X-ray-shear correlations versus tSZ-shear correlations when using a minimal model. They demonstrate that this tension is resolved by marginalizing over unresolved AGN contamination and accounting for non-thermal pressure.
3. Unified Model Validation: They show that a relatively simple hydrodynamic model, once extended to include these astrophysical effects, can simultaneously describe the cross-correlations ($\gamma X$, $\gamma y$) and successfully predict external observables not used in the fit (tSZ auto-correlation, bias-weighted mean pressure, and the $X \times y$ cross-correlation).
4. Quantification of Baryonic Suppression: Using the derived constraints on $M_c$, the authors calculate the suppression factor of the matter power spectrum due to baryonic feedback, comparing it favorably with hydrodynamical simulations (FLAMINGO).

4. Key Results

• Parameter Constraints: After marginalizing over AGN contamination and non-thermal pressure, the joint analysis yields:
  • $\log_{10} M_c = 14.83^{+0.16}_{-0.23}$ (Mass scale where 50% of gas is expelled).
  • $\Gamma = 1.144^{+0.016}_{-0.013}$ (Polytropic index).
  • $\alpha_T = 1.30^{+0.15}_{-0.28}$ (Central temperature ratio).
• Tension Resolution: The tension between the $\Gamma$ values preferred by X-ray data ($\sim 1.22$) and tSZ data ($\sim 1.15$) disappears when the model accounts for the fact that unresolved AGN contribute $\sim 20\%$ to the X-ray signal and when non-thermal pressure is included.
• Predictive Power: The best-fit model derived from the cross-correlations accurately predicts:
  • The tSZ auto-correlation ($C_{yy}$).
  • The bias-weighted mean gas pressure $\langle b P_e \rangle$ measured in other studies.
  • The $X \times y$ cross-correlation (with a PTE of 77.9% when AGN marginalization is applied).
• Impact of Systematics: The study highlights that the choice of tSZ map and scale cuts significantly affects the inferred $\Gamma$. However, the inclusion of AGN marginalization is the most critical step in reconciling the datasets.

5. Significance

• Cosmological Implications: By providing a self-calibrated model for baryonic feedback, this work offers a pathway to reduce systematic uncertainties in Stage-IV weak lensing surveys (like Euclid and LSST). Accurate modeling of baryonic suppression is essential for resolving the $S_8$ tension and measuring cosmological parameters with high precision.
• Astrophysical Insights: The results provide direct observational constraints on the efficiency of AGN feedback (via $M_c$) and the thermodynamic state of the intracluster medium. The finding that $\Gamma \approx 1.14$ suggests gas profiles are slightly steeper than isothermal but less steep than adiabatic, consistent with simulations including feedback.
• Methodological Advance: The paper demonstrates the power of multi-wavelength cross-correlations. By combining tracers with different dependencies on density and temperature (X-ray $\propto \rho^2 T$, tSZ $\propto \rho T$, Shear $\propto \rho$), the authors break degeneracies that plague single-probe analyses.
• Future Directions: The authors outline necessary improvements for future work, including better modeling of metallicity gradients, more flexible AGN clustering models, and corrections for the transition between 1-halo and 2-halo regimes in the halo model.

In conclusion, this paper establishes that while simple hydrodynamic models face tensions when applied to multi-probe data, incorporating realistic astrophysical contaminants (AGN) and physics (non-thermal pressure) allows for a consistent, predictive, and physically motivated description of the hot gas universe.