A comparative test of different pressure profile models in clusters of galaxies using recent ACT data

Imagine the universe is filled with giant, invisible bubbles of hot gas. These aren't just any bubbles; they are galaxy clusters, the largest structures in the cosmos, held together by gravity. Inside these bubbles, the gas is so hot (millions of degrees) that it glows in X-rays and leaves a specific "fingerprint" on the cosmic microwave background (the afterglow of the Big Bang). This fingerprint is called the Sunyaev-Zel'dovich (SZ) effect.

Think of the SZ effect like a shadow cast by a hot stove on a wall. The hotter the stove (the cluster), the darker and more distinct the shadow. By measuring this shadow, astronomers can map out the pressure of the gas inside these clusters.

The Big Question: Is There a "Universal" Recipe?

For a long time, scientists have tried to find a single mathematical "recipe" to describe how this gas pressure is distributed inside any galaxy cluster, regardless of its size or how far away it is. It's like asking: "Is there one single blueprint for how a cake rises, whether it's a tiny cupcake or a massive wedding cake?"

Most previous studies used a specific recipe called the gNFW (generalized Navarro-Frenk-White) model. It's the "standard" recipe in the cookbook. But is it the only recipe that works? Or are there other, perhaps more physically accurate, ways to describe the gas?

The Experiment: A Taste Test with 3,500 Clusters

In this paper, the author, Denis Tramonte, decided to put four different recipes to the test using a massive amount of new data from the Atacama Cosmology Telescope (ACT).

The Ingredients: He gathered data on 3,496 galaxy clusters. That's a huge sample size, covering a wide range of masses and distances (redshifts).
The Method (Stacking): Individual clusters are faint and hard to see clearly. So, the author used a technique called "stacking." Imagine taking 3,500 blurry photos of different people and stacking them on top of each other perfectly aligned. The result is one incredibly sharp, clear photo of an "average" person. Similarly, he stacked the signals of 3,500 clusters to create a clear picture of the "average" pressure profile.
The Contenders: He tested four different mathematical models (recipes) to see which one best matched the "average" shadow he saw:
- The gNFW (UPP): The classic, standard recipe.
- The $\beta$ -model: An older, simpler recipe based on the idea that the gas is in a calm, steady state.
- The Polytropic Model: A recipe based on how gas behaves when compressed (like air in a tire).
- The Exponential Model (EUP): A brand-new, custom recipe the author invented to fix some mathematical quirks in the others.

The Results: A Dead Heat

Here is the surprising twist: All four recipes worked almost equally well.

No Clear Winner: When the author compared the predictions of each recipe against the actual data, they all fit the measurements within the margin of error. It's like trying to guess the flavor of a mystery soup; the salt, pepper, and garlic all tasted just right. You couldn't say for sure which one was the "true" flavor.
The "Universal" Myth: The study found that while these recipes work well on average, they aren't perfectly "universal."
- The Trend: The "average" cluster isn't the same everywhere. The recipe needed to fit the biggest, oldest clusters (low redshift) was slightly different from the one needed for smaller, younger clusters (high redshift).
- The Analogy: Think of it like baking. If you try to use the exact same baking time and temperature for a cupcake and a 10-layer cake, they won't both come out perfect. The biggest, most relaxed clusters seem to have a "steeper" pressure drop-off (the gas falls off faster as you move away from the center) than the younger, messier ones.

Why Does This Matter?

Simplicity vs. Reality: The standard recipe (gNFW) is popular because it's convenient and has been used for decades. However, this study shows that simpler, more physically motivated recipes (like the $\beta$ -model) work just as well. We don't need to stick to the complicated standard just because it's tradition.
The Limits of "Averages": The study highlights a limitation in astronomy. When we look at thousands of clusters at once, we get a smooth average. But individual clusters are messy, turbulent, and unique. Trying to force every single cluster into one "Universal" box might be impossible without making the math so complex that it loses its meaning.
Future Directions: The author suggests that to truly crack the code, we might need to combine data from different sources (like X-ray and SZ data) to break the "degeneracy" (the confusion between different recipes).

The Bottom Line

This paper is a reality check for astronomers. It says: "Stop obsessing over finding the one perfect universal recipe. It might not exist."

Instead, we should accept that different mathematical models can describe the universe just as well as each other. The "Universal Pressure Profile" is a useful tool for general cosmology, but if you want to understand the specific physics of a single cluster, you need to remember that nature is messy, and no single equation can perfectly capture every nuance of a galaxy cluster's hot gas.

In short: We tested four different ways to describe the hot gas in galaxy clusters. They all passed the test. The "standard" way isn't the only way, and the universe is a bit more flexible (and less uniform) than we hoped.

Here is a detailed technical summary of the paper "A comparative test of different pressure profile models in clusters of galaxies using recent ACT data" by Denis Tramonte.

1. Problem Statement

Galaxy clusters are the largest virialized structures in the universe, and their intracluster medium (ICM) contains the majority of their baryonic mass. Characterizing the thermodynamical state of the ICM is crucial for cosmology and astrophysics. The electron pressure profile is the primary tool for this characterization.

While the "Universal Pressure Profile" (UPP), based on a generalized Navarro-Frenk-White (gNFW) functional form, is the standard model adopted in the literature, it suffers from several issues:

Physical Limitations: It exhibits unphysical divergence at the cluster center and strong degeneracies between its slope parameters ( $\alpha, \beta, \gamma$ ), making individual parameter values difficult to interpret physically.
Sample Bias: Previous constraints often relied on high-mass, well-resolved clusters or heterogeneous optical catalogs, potentially introducing selection biases.
Universality Assumption: It is unclear if a single functional form can accurately describe clusters across a wide range of masses and redshifts without residual dependencies.

This study aims to test the efficacy of the UPP against three alternative functional forms (Beta-model, Polytropic, and a new Exponential model) using a large, homogeneous sample of galaxy clusters from the Atacama Cosmology Telescope (ACT).

2. Methodology

Data Sets

Compton Parameter Map: The study utilizes the ACT DR6 Compton parameter ( $y$ ) map, covering $\sim13,000$ deg $^2$ with a resolution of 1.6 arcmin. This map was generated via a blind component separation pipeline combining Planck and ACT data.
Cluster Catalog: The sample consists of 3,496 galaxy clusters from the ACT DR5 catalog (Hilton et al. 2021). These were identified via a blind search on ACT DR5 maps using a multifrequency matched filter.
- Mass Range: $10^{14} $to$ 10^{15.1} M_\odot$.
- Redshift Range: $z \in [0, 2]$ .
- Selection: The catalog is homogeneous (derived from the same SZ data as the map) and uses bias-uncorrected masses derived directly from SZ scaling relations.

Data Processing

Stacking: The 3,496 clusters were stacked on the $y$ -map to extract average angular Compton profiles.
Correction: A critical methodological upgrade over previous work (e.g., Tramonte et al. 2023) involved correcting for pixel area distortion in the equirectangular projection. Pixels at high declinations subtend smaller solid angles; the authors corrected the Right Ascension (RA) separation by $\cos(\delta)$ before stacking to prevent artificial ellipticity in the profiles.
Sub-sampling: The full sample was divided into 16 subsamples based on three mass bins and three redshift bins to test for evolutionary trends.

Theoretical Modeling

The study models the line-of-sight (LoS) integrated electron pressure. The pressure profile $P_e(r)$ is factorized into a global scaling function (dependent on mass $M_{500}$ and redshift $z$ ) and a scaled universal shape function $P(x)$ , where $x = r/r_s$ and $r_s = R_{500}/c_{500}$ .

Four distinct functional forms for $P(x)$ were tested:

Universal Pressure Profile (UPP): A gNFW form: $P(x) = P_0 x^{-\gamma} [1 + x^\alpha]^{(\beta-\gamma)/\alpha}$ . (Fixed $\gamma=0.31$ ).
$\beta$ -Model Profile (BMP): Based on hydrostatic equilibrium and isothermality: $P(x) = P_0 [1 + x^2]^{-3\beta/2}$ .
Polytropic Profile (PTP): Based on a polytropic equation of state with NFW density: $P(x) = P_0 [\frac{1}{x}\ln(1+x)]^{n+1}$ .
Exponential Universal Profile (EUP): A new ansatz: $P(x) = P_0 e^{-\gamma x} (1+x)^{-\zeta}$ . This aims to solve the central divergence of the UPP and reduce parameter degeneracy.

Parameter Estimation

MCMC Approach: A multi-stage Monte Carlo Markov Chain (MCMC) approach was used to constrain parameters ( $c_{500}$ , $P_0$ , and shape parameters).
Degeneracy Breaking: To address strong parameter degeneracies, a two-step fitting process was employed:
1. Fit $c_{500}$ alone while fixing other parameters.
2. Fix the new $c_{500}$ and fit the remaining parameters.
Covariance: To ensure stability against spurious correlations in outer bins, a diagonal-only covariance matrix was used in the final fitting stage, with Fisher-matrix corrections applied to $c_{500}$ uncertainties.

3. Key Contributions

Large-Scale Homogeneous Sample: The use of a single, blind SZ-derived catalog (ACT DR5) avoids the systematic biases associated with merging heterogeneous optical catalogs used in previous population-level studies.
Methodological Refinement: The implementation of declination-dependent pixel area corrections in the stacking procedure significantly improves the accuracy of the reconstructed angular profiles compared to previous studies using ACT DR4.
Comparative Model Testing: Unlike previous works that focused primarily on the UPP, this study rigorously compares the UPP against physically motivated alternatives (BMP, PTP) and a new physically motivated ansatz (EUP).
Resolution of Degeneracies: The multi-stage MCMC strategy effectively breaks the degeneracy between the concentration parameter and the profile shape, providing tighter constraints than previous single-stage fits.

4. Results

Model Performance: All four functional forms (UPP, BMP, PTP, EUP) successfully reproduce the measured angular Compton profiles within their error bars. There is no single model that is statistically favored over the others based on reduced chi-squared ( $\chi^2_r$ ) values.
Parameter Constraints:
- Concentration ( $c_{500}$ ): Best-fit values are generally $>2.0$ , consistent with some recent literature but higher than early UPP estimates.
- Slopes: The $\beta$ -model and Polytropic models yield slopes consistent with previous X-ray and SZ studies.
- UPP Parameters: The best-fit UPP parameters are in broad agreement with literature, though the study confirms the strong degeneracy between parameters.
Residual Trends (Non-Universality):
- While the models fit the data well, there are hints of residual dependencies on mass and redshift.
- Trend: High-mass, low-redshift clusters (typically relaxed, cool-core systems) favor higher amplitudes and steeper profiles (larger slope parameters) compared to low-mass, high-redshift clusters.
- This suggests that a strictly "universal" profile may not capture the full physical diversity of the ICM, particularly the differences between relaxed and disturbed systems.
EUP Viability: The new Exponential Universal Profile (EUP) performs as well as the standard UPP but offers a more physically motivated form (no central divergence) and potentially reduced parameter degeneracy.

5. Significance and Conclusions

Limitations of SZ-Only Studies: The study concludes that population-level studies using SZ data alone lack the constraining power to definitively discriminate between different pressure profile models. The line-of-sight integration of the Compton parameter allows different functional forms to yield nearly identical predictions.
Re-evaluating the UPP: The gNFW (UPP) parametrization, while convenient and historically standard, is not physically superior to models like the $\beta$ -model or Polytropic profile. Its ad-hoc nature and intrinsic degeneracies limit its utility for extracting deep physical insights into the ICM.
Future Directions: The authors suggest that combining SZ data (proportional to integrated density) with X-ray data (proportional to integrated density squared) is necessary to break degeneracies and constrain the radial temperature profile.
Universal Model Validity: The assumption of a single universal pressure profile is challenged by the observed trends in amplitude and slope with mass and redshift. While universal models remain useful for cosmological scaling relations, they may fail to capture the specific thermodynamical states of individual cluster populations when high precision is required.

In summary, this paper demonstrates that while the UPP remains a valid descriptive tool, alternative, more physically motivated models are equally effective. The observed trends suggest that the "universal" nature of the ICM pressure profile is an approximation that breaks down when analyzing specific mass and redshift regimes in detail.