Galaxy Clusters and Hydrostatic Dynamics

Scrollable formula. Use Left and Right arrow keys to reveal hidden content.

Scrollable table. Use Left and Right arrow keys to reveal additional columns.

Scrollable figure. Use Left and Right arrow keys to inspect hidden content.

Galaxy clusters offer several gravitational tracers in the same object: galaxy velocities respond to the potential, hot gas responds through pressure support, and lensing responds to light deflection. Agreement can test a common weak-field metric map; disagreement can signal nonthermal pressure, triaxial geometry, line-of-sight structure, selection, or modified gravity. A cluster “mass” is therefore not one directly observed number. It is a model-dependent reconstruction from each tracer, and their covariance and validity domains determine whether the comparison tests gravity. The evidence assessment is current through 2026-08-19.

Required background. Classical Field Theory supplies potentials and stress-supported matter, while Probability and Statistical Inference supplies joint likelihoods, population selection, covariance, and model checks.

Helpful background. Weak-Field Gravity supplies dynamical and lensing potentials, and Black-Hole Astrophysics supplies gas, plasma, and galaxy-population context.

Hydrostatic, dynamical, and lensing reconstructions

For a spherical ideal gas in hydrostatic equilibrium, thermal pressure obeys $dP_{\rm th}/dr=-\rho_g\,d\Phi/dr$ if it supplies all support. Using $P_{\rm th}=n_gk_{\rm B}T$ and $\rho_g=\mu m_p n_g$ gives

M_{\rm HSE}(<r) =-\frac{rk_{\rm B}T(r)}{G\mu m_p} \left[\frac{d\ln n_g}{d\ln r}+\frac{d\ln T}{d\ln r}\right],

where $n_g$ , $T$ , and $\mu m_p$ are particle number density, temperature, and mean particle mass; composition is assumed constant over the derivative. The global convention sets $k_{\rm B}=1$ , but it is restored here because cluster temperatures are reported instrumentally as $k_{\rm B}T$ . The factor $1/G$ converts a measured force into a Newtonian mass. In GR’s weak-field spherical limit this is enclosed gravitating mass; in another force law it is an effective dynamical mass, not automatically the actual enclosed matter.

An isothermal gas with $n_g\propto r^{-2}$ illustrates the result: $M_{\rm HSE}=2rk_{\rm B}T/(G\mu m_p)$ , positive because density and pressure decrease outward. If a scalar nonthermal pressure also supports a stationary spherical atmosphere, write $P_{\rm tot}=P_{\rm th}+P_{\rm nt}$ and $f_{\rm nt}=P_{\rm nt}/P_{\rm tot}$ . Differentiation gives

M_{\rm HSE} =(1-f_{\rm nt})M_{\rm dyn} +\frac{r^2P_{\rm tot}}{G\rho_g}\frac{df_{\rm nt}}{dr}, \qquad M_{\rm dyn}=\frac{r^2}{G}\frac{d\Phi}{dr}.

For constant $f_{\rm nt}=0.2$ , thermal pressure alone returns $0.8M_{\rm dyn}$ . A radially varying fraction adds the derivative term, so the fraction itself is not a universal multiplicative mass correction. Turbulence, bulk flow, cosmic rays, magnetic stress, and asphericity can require still more general support or geometric terms; a merging cluster need not obey this equilibrium reduction.

Galaxy velocities obey a collisionless Jeans equation with orbital anisotropy and membership uncertainties. Weak and strong lensing instead infer projected combinations of metric potentials and matter along the line of sight. In general relativity with negligible anisotropic stress, the dynamical and lensing maps are tied; in modified gravity, the potential governing nonrelativistic motion need not equal the light-deflection combination. Comparing them is consequently valuable only when gas physics, tracer anisotropy, projection, and source-redshift calibration are independently controlled.

The X-ray route from surface brightness and temperature to mass, together with its equilibrium assumptions, is developed in Ettori et al. 2013, §§2–4 and Eqs. (1)–(8). Cluster lensing reconstruction and its mass-sheet, redshift, shear, and projection systematics are reviewed in Umetsu 2020, §§3–6. These sources make clear that agreement among catalog masses does not provide an independent test if one proxy was calibrated against another under the same gravitational assumptions.

Data products and joint likelihood

A reproducible comparison pins a named cluster sample and release: X-ray event files and response matrices, background and point-source masks, temperature pipeline; Sunyaev–Zeldovich maps and bandpass or beam calibration; shear and strong-lensing catalogs with photometric-redshift and multiplicative-shear calibration; galaxy spectroscopy and membership; and cluster centers and redshifts. Instrument cross-calibration, PSF, effective area, beam, absolute astrometry, source-redshift distribution, and velocity zero points enter explicitly.

For cluster $j$ with latent mass profile $M_j$ , shape $q_j$ , orientation $\Omega_j$ , and nonthermal fraction $f_{{\rm nt},j}$ , a schematic joint likelihood is

\begin{aligned} p(d_j\mid\Theta_j,\mathcal M) &=p(d_j^{\rm X}\mid\Theta_j,\mathcal M) \,p(d_j^{\rm SZ}\mid\Theta_j,\mathcal M)\\ &\quad\times p(d_j^{\rm lens}\mid\Theta_j,\mathcal M) \,p(d_j^{\rm gal}\mid\Theta_j,\mathcal M), \end{aligned}

only when the measurement errors are conditionally independent after including all shared calibration and source variables in $\Theta_j$ . The same geometry, orientation, center, gas profile, and potential enter every relevant factor; $\beta_j$ denotes galaxy orbital anisotropy. Integrating shared nuisance variables usually correlates the probes. Residual cross-probe noise covariance requires a joint likelihood instead of this product. A modified-gravity model must map one matter distribution into each potential and gas or galaxy response; a GR analysis must specify halo, baryon, and pressure models. Priors on triaxiality, nonthermal support, anisotropy, concentration, and line-of-sight structure can drive the comparison and must be reported.

Selection includes X-ray flux or morphology thresholds, SZ signal-to-noise, optical richness, lensing detectability, redshift, masks, source quality, and the definition of “relaxed.” A relaxed subsample improves hydrostatic validity but may not represent the parent population. Completeness and the observable–mass scatter belong in a population-normalized likelihood. Posterior-predictive profiles, not only integrated masses, should be checked.

Results and limitations through 2026-08-19

Position or evidence	Result or claim	Limitation or counterevidence
Relaxed-cluster X-ray, SZ, lensing, and dynamical reconstructions can be mutually consistent under GR plus dark matter	Supports one gravitational potential after calibrated baryonic support	Shared mass priors, cross-calibration, triaxial projection, and hydrostatic bias reduce independence
Disturbed clusters exhibit significant proxy disagreement	Nonequilibrium gas and geometry can explain apparent dynamical–lensing mismatch	A flexible disturbance correction can also hide a true force-law difference unless validated on simulations
Modified gravity can change dynamical versus lensing potentials or screening with radius	Multi-tracer profiles can in principle discriminate gravitational slip or screening	Cluster baryons, neutrinos or additional unseen mass, and nonlinear screening remain theory- and environment-dependent

The conservative position is that calibrated multi-probe cluster observations are broadly modeled by general relativity with dark matter and imperfect gas support. The competing position seeks a coherent radial mismatch between matter motion and light deflection predicted by a complete modified theory. Individual discrepant clusters are not persuasive because mergers and orientation supply known alternatives. It remains unresolved whether future percent-level multi-probe profiles can separate those effects without simulation-dependent calibration dominating the result.

Adversarial multi-probe checks

An independent check applies every estimator to mock relaxed and disturbed clusters with known three-dimensional potential, instrument response, and line-of-sight structure. It must recover the input profile and achieve interval coverage separately for morphology classes. Leave-one-probe-out prediction then asks whether X-ray plus SZ predicts lensing, whether lensing predicts galaxy kinematics, and whether discrepancies correlate with disturbance indicators.

An adversarial analysis changes the center, triaxiality prior, nonthermal-pressure model, source-redshift distribution, membership algorithm, and halo or screening prescription. A gravity signal must appear with the correlated radial form predicted by the theory across multiple probes; a constant mass offset is more naturally a calibration or support degeneracy. The rerun package needs released maps or catalogs, masks, response files, selection function, covariance, code version, priors, and mock controls.

What could change the conclusion. The conclusion could change after a cross-calibrated multi-instrument release, a blind cluster-simulation challenge validates pressure and geometry recovery, a larger selection-characterized sample yields a coherent radial dynamical–lensing mismatch, or a complete screened-gravity model predicts the same pattern without object-by-object tuning. A single high-significance mass-ratio outlier is not enough.

Exercise

If unmodeled nonthermal pressure supplies part of the radial gas support, does the thermal hydrostatic estimator overestimate or underestimate the true enclosed mass?

Solution

The full equilibrium equation contains the gradient of thermal plus nonthermal pressure. Using only the thermal gradient attributes too little support to the gas. Substituting that smaller gradient into the hydrostatic formula therefore gives an enclosed mass below the true value, provided the omitted nonthermal pressure decreases outward in the usual supporting sense. The sign can fail for an arbitrary transient flow, which is why disturbed systems require a dynamical model rather than a universal correction factor.

References

Ettori, S., et al. (2013), “Mass Profiles of Galaxy Clusters from X-ray Analysis,” Space Science Reviews 177, 119–154. DOI.
Umetsu, K. (2020), “Cluster–galaxy weak lensing,” Astronomy and Astrophysics Review 28, 7. DOI.