Universal Pressure Profile (UPP) Overview

• Universal Pressure Profile (UPP) is a self-similar model describing intracluster electron pressure using scaled radius and pressure for cross-cluster comparison.
• It employs a generalized NFW (gNFW) parameterization, facilitating thermal SZ analysis, X-ray/SZ profile comparisons, and halo-model forecasts.
• Research reveals that UPP’s universality is approximate, with variations arising from entropy differences, hydrostatic bias, and mass calibration corrections.

Searching arXiv for the cited UPP papers and related work. {"query":"(Aslanbeigi et al., 2013) Universal Pressure Profile intracluster pressure profile methodology and first applications", "max_results": 5} {"query":"(Lapi et al., 2011) Intracluster Plasma Universal Pressure Profile", "max_results": 5} The Universal Pressure Profile (UPP) is a self-similar description of the intracluster electron pressure profile of the hot intracluster medium (ICM), usually expressed in scaled radius and scaled pressure so that clusters of different mass and redshift can be compared within a common dimensionless framework. In its standard usage, the UPP combines a characteristic pressure normalization with a generalized Navarro–Frenk–White (gNFW) shape, and it has become a central model for thermal Sunyaev–Zeldovich (tSZ) analyses, X-ray/SZ profile comparisons, matched-filter cluster detection, mass calibration, and halo-model forecasts. At the same time, the literature treats “universality” as approximate rather than exact: some works regard the UPP as a statistically useful mean profile, whereas others emphasize entropy-class bimodality, hydrostatic-mass bias, or residual mass and redshift dependence (Aslanbeigi et al., 2013, Lapi et al., 2011, He et al., 2020).

1. Canonical definition and self-similar scaling

In the self-similar formulation, different clusters are assumed to share the same dimensionless radial shape once radius is scaled by a characteristic size and pressure by a characteristic amplitude. One explicit statement of this assumption is

$P_e^{(a)}(r)=P\!\left(r/r_c^{(a)}\right)P_c^{(a)},$

with $r_c^{(a)}=R_{500}^{(a)}$. The overdensity radius is defined through

$$M_\Delta=\frac{4\pi}{3}\Delta R_\Delta^3 \rho_{\rm crit}(z), \qquad \rho_{\rm crit}(z)=\frac{3H(z)^2}{8\pi G}.$$

Using $R_{500}$ as the scaling radius, the characteristic pressure normalization is parameterized as

$P_c^{(a)}=1.65\times 10^{-3}E(z_a)^{8/3}h_{70}^2 \left[\frac{M_{500}^{(a)}{3\times10^{14}h_{70}^{-1}M_\odot}\right]^{2/3+\delta}\ {\rm keV\,cm^{-3},$

with $\delta=0$ for standard self-similar scaling and $\delta=0.12$ for the modified Arnaud et al. scaling used to compare directly to Planck (Aslanbeigi et al., 2013).

The dimensionless UPP shape is usually written in gNFW form,

$$\mathbb{P}(x)=\frac{P_0}{(c_{500}x)^\gamma\left[1+(c_{500}x)^\alpha\right]^{(\beta-\gamma)/\alpha}},$$

where $x=r/R_{500}$ or, in some formulations, $x=r/r_s$ with $r_c^{(a)}=R_{500}^{(a)}$0. In this parameterization, $r_c^{(a)}=R_{500}^{(a)}$1 is the normalization, $r_c^{(a)}=R_{500}^{(a)}$2 is the concentration-like scale parameter, $r_c^{(a)}=R_{500}^{(a)}$3 controls the inner slope, $r_c^{(a)}=R_{500}^{(a)}$4 controls the transition region, and $r_c^{(a)}=R_{500}^{(a)}$5 controls the outer slope. This gNFW family is the standard analytic backbone of UPP work; flexible fitting codes such as PreProFit use the full gNFW model and treat the UPP as a special case rather than as a hard-coded default (Castagna et al., 2019).

The observational connection is through the tSZ effect, which probes the line-of-sight integral of electron pressure. One standard expression is

$r_c^{(a)}=R_{500}^{(a)}$6

with

$r_c^{(a)}=R_{500}^{(a)}$7

Because the SZ observable is linear in pressure and redshift-independent in surface brightness, it is particularly well suited to testing the scaled-pressure hypothesis underlying the UPP (Lapi et al., 2011).

2. Measurement strategies and empirical profile reconstruction

A key methodological advance in UPP measurement was the replacement of simple stacking by an optimal multi-frequency template fit. In one explicit implementation, the tSZ signal is decomposed into radial bins,

$r_c^{(a)}=R_{500}^{(a)}$8

and the data are fitted with a Gaussian likelihood

$r_c^{(a)}=R_{500}^{(a)}$9

This framework retains the full frequency-channel resolution and sensitivity by fitting all WMAP Q/V/W channels simultaneously rather than degrading the data to a single low-resolution map, and it is designed to account explicitly for primary CMB correlations and instrumental noise correlations (Aslanbeigi et al., 2013).

Applied to WMAP 9-year data and the MCXC catalog, this method fitted 8 radial bins from $$M_\Delta=\frac{4\pi}{3}\Delta R_\Delta^3 \rho_{\rm crit}(z), \qquad \rho_{\rm crit}(z)=\frac{3H(z)^2}{8\pi G}.$$0 to $$M_\Delta=\frac{4\pi}{3}\Delta R_\Delta^3 \rho_{\rm crit}(z), \qquad \rho_{\rm crit}(z)=\frac{3H(z)^2}{8\pi G}.$$1. For all 1743 MCXC clusters, the pressure profile was detected at $$M_\Delta=\frac{4\pi}{3}\Delta R_\Delta^3 \rho_{\rm crit}(z), \qquad \rho_{\rm crit}(z)=\frac{3H(z)^2}{8\pi G}.$$2 for $$M_\Delta=\frac{4\pi}{3}\Delta R_\Delta^3 \rho_{\rm crit}(z), \qquad \rho_{\rm crit}(z)=\frac{3H(z)^2}{8\pi G}.$$3 and $$M_\Delta=\frac{4\pi}{3}\Delta R_\Delta^3 \rho_{\rm crit}(z), \qquad \rho_{\rm crit}(z)=\frac{3H(z)^2}{8\pi G}.$$4 for $$M_\Delta=\frac{4\pi}{3}\Delta R_\Delta^3 \rho_{\rm crit}(z), \qquad \rho_{\rm crit}(z)=\frac{3H(z)^2}{8\pi G}.$$5; for the resolved subsample of 162 clusters, the significance was $$M_\Delta=\frac{4\pi}{3}\Delta R_\Delta^3 \rho_{\rm crit}(z), \qquad \rho_{\rm crit}(z)=\frac{3H(z)^2}{8\pi G}.$$6 or $$M_\Delta=\frac{4\pi}{3}\Delta R_\Delta^3 \rho_{\rm crit}(z), \qquad \rho_{\rm crit}(z)=\frac{3H(z)^2}{8\pi G}.$$7. The profile was constrained out to $$M_\Delta=\frac{4\pi}{3}\Delta R_\Delta^3 \rho_{\rm crit}(z), \qquad \rho_{\rm crit}(z)=\frac{3H(z)^2}{8\pi G}.$$8 in the stacked/fit sense, but the measurements were only individually significant out to about $$M_\Delta=\frac{4\pi}{3}\Delta R_\Delta^3 \rho_{\rm crit}(z), \qquad \rho_{\rm crit}(z)=\frac{3H(z)^2}{8\pi G}.$$9, and beyond $R_{500}$0 the best-fit bins became consistent with zero and could even be negative, leading to the conclusion that there was essentially no signal beyond that scale. The same analysis compared its binned measurements to the Planck best-fit GNFW/UPP profile and described the agreement as good (Aslanbeigi et al., 2013).

Later SZ-only work extended the empirical test of the UPP to much larger samples and higher-quality multi-instrument combinations. A joint Fourier-space analysis of 461 SPT and Planck clusters performed a deconvolved average-profile reconstruction and fitted a gNFW model with $R_{500}$1 fixed. For the full sample it reported

$R_{500}$2

with $R_{500}$3 for 68 degrees of freedom and reduced $R_{500}$4. The inner profile agreed well with the classic UPP, while the outer profile agreed better with the Planck intermediate paper V profile. The same study found no strong indication of redshift evolution but did report a significant difference between low- and high-mass subsamples, indicating that universality is accurate only to a certain level of approximation (Melin et al., 2023).

3. Physical interpretation and the limits of universality

The strongest conceptual challenge to a strictly universal low-redshift pressure law comes from entropy-based interpretations of the ICM. One such framework argues that nearby clusters do not occupy a single pressure family but instead separate into low-entropy (LE) and high-entropy (HE) systems. In that description, LE clusters have $R_{500}$5 and $R_{500}$6, with steeper pressure profiles in both core and outskirts, whereas HE clusters have $R_{500}$7 and $R_{500}$8 to $R_{500}$9, with shallower pressure profiles throughout. The underlying entropy run is modeled as

$P_c^{(a)}=1.65\times 10^{-3}E(z_a)^{8/3}h_{70}^2 \left[\frac{M_{500}^{(a)}{3\times10^{14}h_{70}^{-1}M_\odot}\right]^{2/3+\delta}\ {\rm keV\,cm^{-3},$0

and the pressure follows from hydrostatic equilibrium,

$P_c^{(a)}=1.65\times 10^{-3}E(z_a)^{8/3}h_{70}^2 \left[\frac{M_{500}^{(a)}{3\times10^{14}h_{70}^{-1}M_\odot}\right]^{2/3+\delta}\ {\rm keV\,cm^{-3},$1

Within this picture, the standard UPP is interpreted as an ensemble average over two physically distinct cluster populations rather than a single exact law valid for all rich clusters at all epochs (Lapi et al., 2011).

This entropy-based critique is coupled to a specific evolutionary hypothesis. The fraction of HE clusters is argued to increase with redshift, so the cluster population is expected to converge toward a truly universal HE-like template at $P_c^{(a)}=1.65\times 10^{-3}E(z_a)^{8/3}h_{70}^2 \left[\frac{M_{500}^{(a)}{3\times10^{14}h_{70}^{-1}M_\odot}\right]^{2/3+\delta}\ {\rm keV\,cm^{-3},$2. That is a claim of redshift-dependent universality: non-universal at low $P_c^{(a)}=1.65\times 10^{-3}E(z_a)^{8/3}h_{70}^2 \left[\frac{M_{500}^{(a)}{3\times10^{14}h_{70}^{-1}M_\odot}\right]^{2/3+\delta}\ {\rm keV\,cm^{-3},$3, increasingly universal at high $P_c^{(a)}=1.65\times 10^{-3}E(z_a)^{8/3}h_{70}^2 \left[\frac{M_{500}^{(a)}{3\times10^{14}h_{70}^{-1}M_\odot}\right]^{2/3+\delta}\ {\rm keV\,cm^{-3},$4. A plausible implication is that some apparent tensions in UPP fits reflect sample composition rather than failure of the scaled-profile paradigm itself (Lapi et al., 2011).

A complementary physical interpretation treats the observed UPP as the end state of entropy injection applied to an initially gravitational entropy profile. In that formulation,

$P_c^{(a)}=1.65\times 10^{-3}E(z_a)^{8/3}h_{70}^2 \left[\frac{M_{500}^{(a)}{3\times10^{14}h_{70}^{-1}M_\odot}\right]^{2/3+\delta}\ {\rm keV\,cm^{-3},$5

and the final entropy is modeled in Lagrangian gas-fraction space as

$P_c^{(a)}=1.65\times 10^{-3}E(z_a)^{8/3}h_{70}^2 \left[\frac{M_{500}^{(a)}{3\times10^{14}h_{70}^{-1}M_\odot}\right]^{2/3+\delta}\ {\rm keV\,cm^{-3},$6

The best-fit $P_c^{(a)}=1.65\times 10^{-3}E(z_a)^{8/3}h_{70}^2 \left[\frac{M_{500}^{(a)}{3\times10^{14}h_{70}^{-1}M_\odot}\right]^{2/3+\delta}\ {\rm keV\,cm^{-3},$7 entropy-injection prescription is

$P_c^{(a)}=1.65\times 10^{-3}E(z_a)^{8/3}h_{70}^2 \left[\frac{M_{500}^{(a)}{3\times10^{14}h_{70}^{-1}M_\odot}\right]^{2/3+\delta}\ {\rm keV\,cm^{-3},$8

The explicit conclusion is that neither a constant entropy floor nor a purely multiplicative global entropy increase is sufficient; both an overall entropy increase and extra central entropy injection are required, suggesting a combination of preheating and in situ AGN/SN feedback. In the same model, self-similar evolution of the UPP corresponds to

$P_c^{(a)}=1.65\times 10^{-3}E(z_a)^{8/3}h_{70}^2 \left[\frac{M_{500}^{(a)}{3\times10^{14}h_{70}^{-1}M_\odot}\right]^{2/3+\delta}\ {\rm keV\,cm^{-3},$9

so the effective entropy injection decreases with redshift (Nath et al., 2011).

Observational tests of redshift evolution remain mixed. A Planck tSZ stacking analysis around SDSS DR7 luminous red galaxies found that the original UPP redshift scaling $\delta=0$0 did not fit all three redshift bins equally well at small angular scales, and that introducing

$\delta=0$1

improved the joint fit with $\delta=0$2. This suggests that the original UPP model with less redshift-dependence may provide a better fit to the stacked thermal SZ profile in that data set (Gong et al., 2019).

4. Mass calibration, hydrostatic bias, and debiased profiles

A major revision of the UPP concerns the mass scale used to define $\delta=0$3 and $\delta=0$4. Because the original X-ray-based calibration used hydrostatic masses, any systematic underestimation of mass propagates directly into the normalized pressure profile. A simulation-based recalibration using IllustrisTNG, BAHAMAS, and MACSIS fitted the relation between “true” and X-ray masses and found a mass-dependent hydrostatic bias: for REXCESS clusters, hydrostatic masses are underestimated by about 7% on average, rising from roughly 0% at $\delta=0$5 to about 15% at $\delta=0$6. After correction, the inferred masses are $\delta=0$7 higher and the normalized pressures are $\delta=0$8 lower than in the original UPP calibration (He et al., 2020).

The resulting Debiased Pressure Profile (DPP) is written as

$\delta=0$9

with best-fit GNFW parameters

$\delta=0.12$0

A central claim of this recalibration is that the extra mass-dependent correction term used in the original UPP is no longer required: after debiasing, $\delta=0.12$1 is consistent with zero at all radii. The recalibrated profile is about 5% lower in the center, about 20% lower by $\delta=0.12$2, and nearly 35% lower in the outskirts relative to the original UPP. When propagated to SZ statistics, the DPP yields an integrated $\delta=0.12$3–$\delta=0.12$4 relation with only minor deviations from self-similar scaling and lowers the thermal SZ angular power spectrum by about 30–40%, bringing it broadly into agreement with Planck under nominal cosmological parameters (He et al., 2020).

Independent analyses likewise identify mass calibration as the dominant systematic. In the WMAP/MCXC study, MCXC and ESZ mass estimates differed systematically by about 12% on average. That discrepancy affected the pressure profile at roughly the $\delta=0.12$5 level but induced a systematic shift of about 20% in the inferred gas mass fraction, large enough to dominate future high-precision analyses unless mass calibration improves (Aslanbeigi et al., 2013).

The mass-profile covariance has also been treated directly as a joint inference problem. A CHEX-MATE analysis fitted the UPP shape, pressure normalization, intrinsic scatter, and individual cluster $\delta=0.12$6 values simultaneously, emphasizing that changes in $\delta=0.12$7 alter both the vertical normalization through $\delta=0.12$8 and the horizontal scaling through $\delta=0.12$9. For the case with free $$\mathbb{P}(x)=\frac{P_0}{(c_{500}x)^\gamma\left[1+(c_{500}x)^\alpha\right]^{(\beta-\gamma)/\alpha}},$$0, the best-fit parameters were approximately

$$\mathbb{P}(x)=\frac{P_0}{(c_{500}x)^\gamma\left[1+(c_{500}x)^\alpha\right]^{(\beta-\gamma)/\alpha}},$$1

The fitted $$\mathbb{P}(x)=\frac{P_0}{(c_{500}x)^\gamma\left[1+(c_{500}x)^\alpha\right]^{(\beta-\gamma)/\alpha}},$$2 was reported as consistent with self-similar evolution, and fixing the UPP to different literature profiles was shown to shift inferred masses by roughly 10–50%. This suggests that a self-consistent UPP requires joint modeling of the absolute mass scale, the $$\mathbb{P}(x)=\frac{P_0}{(c_{500}x)^\gamma\left[1+(c_{500}x)^\alpha\right]^{(\beta-\gamma)/\alpha}},$$3 scaling, and the profile shape (Muñoz-Echeverría et al., 21 Oct 2025).

5. Role in SZ surveys, forward modeling, and cross-correlation analyses

In survey pipelines, the UPP is often not merely a descriptive profile but an operational template. The ACTPol two-season SZ catalog used the UPP both in the matched filter employed for cluster detection and in the mass–SZ scaling relation used to infer $$\mathbb{P}(x)=\frac{P_0}{(c_{500}x)^\gamma\left[1+(c_{500}x)^\alpha\right]^{(\beta-\gamma)/\alpha}},$$4. The chosen reference filter corresponded to a UPP cluster with

$$\mathbb{P}(x)=\frac{P_0}{(c_{500}x)^\gamma\left[1+(c_{500}x)^\alpha\right]^{(\beta-\gamma)/\alpha}},$$5

or equivalently $$\mathbb{P}(x)=\frac{P_0}{(c_{500}x)^\gamma\left[1+(c_{500}x)^\alpha\right]^{(\beta-\gamma)/\alpha}},$$6. The observable $$\mathbb{P}(x)=\frac{P_0}{(c_{500}x)^\gamma\left[1+(c_{500}x)^\alpha\right]^{(\beta-\gamma)/\alpha}},$$7 was converted into $$\mathbb{P}(x)=\frac{P_0}{(c_{500}x)^\gamma\left[1+(c_{500}x)^\alpha\right]^{(\beta-\gamma)/\alpha}},$$8 using the Arnaud et al. scaling relation with intrinsic scatter $$\mathbb{P}(x)=\frac{P_0}{(c_{500}x)^\gamma\left[1+(c_{500}x)^\alpha\right]^{(\beta-\gamma)/\alpha}},$$9 and a Tinker et al. halo-mass-function prior, and the mass-function correction lowered catalog masses by about 16%. A comparison with richness-based weak-lensing masses yielded

$x=r/R_{500}$0

showing that UPP-based SZ masses were on average lower than the weak-lensing-calibrated masses (Hilton et al., 2017).

For profile fitting at the map level, PreProFit provides a Bayesian forward-modeling framework based on the gNFW pressure law,

$x=r/R_{500}$1

followed by Abel projection,

$x=r/R_{500}$2

beam convolution, transfer-function filtering, and comparison to radial surface-brightness data with a Gaussian $x=r/R_{500}$3 likelihood. The code does not hard-code a fixed UPP parameter set, but it is explicitly suited to testing whether data are consistent with a UPP-like gNFW shape (Castagna et al., 2019).

Cross-correlation studies use the UPP as the gas-sector input to halo-model calculations. A Planck–RCSLenS analysis fitted UPP parameters to $x=r/R_{500}$4 and $x=r/R_{500}$5 and found that the concentration parameter $x=r/R_{500}$6 was the best-constrained quantity; for the $x=r/R_{500}$7 estimator, it obtained $x=r/R_{500}$8 for WMAP-9 and $x=r/R_{500}$9 for Planck-2018. Propagating the posterior to pressure space produced a factor of 3 uncertainty in the profile shape and magnitude. The same study found that most of the cross-correlation signal came from low redshift, halo masses $x=r/r_s$0–$x=r/r_s$1, and inner halo radii $x=r/r_s$2, indicating that tSZ–lensing cross-correlations test the regime where UPP assumptions are most astrophysically consequential (Ma et al., 2020).

6. Recent comparative tests and broader modeling uses

Large stacked-sample analyses continue to test the breadth of the UPP’s validity. Using Planck 2015 MILCA and ACT-DR4 $x=r/r_s$3 maps together with merged optical cluster catalogs, one study analyzed 23,820 clusters over

$x=r/r_s$4

and fitted a halo-model UPP including hydrostatic bias and miscentering. For the full cluster sample it obtained

$x=r/r_s$5

with Planck and

$x=r/r_s$6

with ACT, together with loose hydrostatic-bias constraints $x=r/r_s$7. Its main conclusion was that there is no compelling evidence for a residual mass or redshift dependence within current uncertainties, extending the empirical validity of the UPP over a much broader domain than earlier tests (Tramonte et al., 2023).

That conclusion is not uncontested. A later ACT-based comparison of four profile families—gNFW/UPP, $x=r/r_s$8-model, polytropic, and exponential—reported that all four reproduce the stacked $x=r/r_s$9-profiles within their error bars, with no compelling indication of a preferred model. In that analysis, the UPP fit the data well but was not clearly favored, and the authors emphasized strong parameter degeneracies, the fact that SZ-only population studies cannot sharply distinguish 3D pressure laws, and hints of residual subsample dependence favoring higher amplitudes and steeper profiles in high-mass, low-redshift clusters. This suggests that the universal pressure assumption is useful as an average approximation but can break down whenever high precision is required (Tramonte, 13 Jan 2026).

The UPP has also moved beyond cluster-profile fitting into numerical ICM environment modeling. Relativistic MHD simulations of radio-galaxy jets have used UPP-based hydrostatic, spherically symmetric atmospheres in which the cluster is described primarily by $r_c^{(a)}=R_{500}^{(a)}$00, with

$r_c^{(a)}=R_{500}^{(a)}$01

and a realistic temperature profile to derive density, magnetic field, Faraday rotation, and synthetic polarization observables. In this context, the UPP is treated as a representative general cluster atmosphere, and comparisons with $r_c^{(a)}=R_{500}^{(a)}$02-models indicate that the chosen atmospheric prescription affects lobe morphology, core clearing, and the partition of energy between radio lobes and shocked ambient gas (Stimpson et al., 2023). A subsequent polarimetric extension derived synthetic Stokes $r_c^{(a)}=R_{500}^{(a)}$03 and $r_c^{(a)}=R_{500}^{(a)}$04 maps and rotation-measure products from the same UPP atmospheres, arguing that such models provide the most realistic spherically symmetric baseline yet for interpreting cluster Faraday screens in AGN environments (Stimpson et al., 12 Aug 2025).

Taken together, these developments establish the UPP less as an immutable law than as a standard organizing model: a self-similar gNFW pressure prescription that remains central to SZ and X-ray cluster science, but whose interpretation depends sensitively on entropy structure, sample selection, mass calibration, and the precision threshold of the application.