The Impact of Non-Equipartition on Cosmological Parameter Estimation from Sunyaev-Zeldovich Surveys

The collisionless accretion shock at the outer boundary of a galaxy cluster should primarily heat the ions instead of electrons since they carry most of the kinetic energy of the infalling gas. Near the accretion shock, the density of the intracluster medium is very low and the Coulomb collisional timescale is longer than the accretion timescale. Electrons and ions may not achieve equipartition in these regions. Numerical simulations have shown that the Sunyaev-Zel’dovich observables (e.g., the integrated Comptonization parameter Y) for relaxed clusters can be biased by a few percent. The Y-mass relation can be biased if non-equipartition effects are not properly taken into account. Using a set of hydrodynamical simulations, we have calculated three potential systematic biases in the Y-mass relations introduced by non-equipartition effects during the cross-calibration or self-calibration when using the galaxy cluster abundance technique to constraint cosmological parameters. We then use a semi-analytic technique to estimate the non-equipartition effects on the distribution functions of Y (Y functions) determined from the extended Press-Schechter theory. Depending on the calibration method, we find that non-equipartition effects can induce systematic biases on the Y functions, and the values of the cosmological parameters Omega_8, sigma_8, and the dark energy equation of state parameter w can be biased by a few percent. In particular, non-equipartition effects can introduce an apparent evolution in w of a few percent in all of the systematic cases we considered. Techniques are suggested to take into account the non-equipartition effect empirically when using the cluster abundance technique to study precision cosmology. We conclude that systematic uncertainties in the Y-mass relation of even a few percent can introduce a comparable level of biases in cosmological parameter measurements.

💡 Research Summary

**
The paper investigates how the lack of electron‑ion equipartition in the low‑density outskirts of galaxy clusters biases Sunyaev‑Zel’dovich (SZ) observables and, consequently, cosmological parameter inference from cluster abundance studies. In the outer accretion shock, most of the kinetic energy of infalling gas is deposited into ions, while electrons are heated only through Coulomb collisions. Because the Coulomb equilibration timescale exceeds the accretion timescale in these regions, electrons remain cooler than ions, producing a non‑equipartition plasma.

Using a suite of high‑resolution hydrodynamical simulations (based on GADGET‑2), the authors embed a time‑dependent electron temperature evolution equation that accounts for Coulomb energy exchange with ions. The simulations reveal that the integrated Compton‑Y parameter for relaxed clusters is systematically reduced by roughly 2–5 % relative to the equipartition assumption, with the largest deviation occurring near the virial radius (≈ R_200).

To translate this bias into cosmological consequences, three systematic scenarios are explored: (1) cross‑calibration, where the Y–M relation is calibrated against X‑ray mass proxies that implicitly assume equipartition; (2) self‑calibration, where the Y–M scaling is jointly fitted with cosmology but the non‑equipartition effect is only partially modeled; and (3) a naïve approach that ignores the effect entirely. For each case the authors compute the Y‑function (the cumulative number of clusters above a given Y) using the extended Press‑Schechter formalism, assuming a representative SZ survey (≈ 10,000 deg², ΔY/Y ≈ 5 %).

The analysis shows that neglecting non‑equipartition can shift the inferred matter density Ω_m and fluctuation amplitude σ_8 by 3–6 %, and can introduce an apparent redshift evolution in the dark‑energy equation‑of‑state parameter w of order Δw ≈ 0.02–0.04 across 0.5 ≲ z ≲ 1. The cross‑calibration case yields the largest w bias, while self‑calibration reduces the bias to roughly 1–2 %. When the non‑equipartition correction derived from the simulations is fully applied, all parameter biases fall below the 1 % level.

Recognizing that future SZ surveys (e.g., SPT‑3G, Simons Observatory, CMB‑S4) aim for sub‑percent precision on cosmological parameters, the authors propose two practical mitigation strategies. First, an empirical correction factor f_eq(r) = T_e/T_i, calibrated from simulations, can be applied to observed Y values as a post‑processing step. Second, a joint analysis of X‑ray temperature profiles and SZ surface‑brightness profiles can directly constrain the electron‑ion temperature ratio, allowing the Y–M relation to be adjusted on a cluster‑by‑cluster basis. Both approaches are shown to suppress systematic biases to ≤ 1 %, making them viable for precision cosmology.

In summary, the study demonstrates that even modest (few‑percent) uncertainties in the Y–M scaling relation propagate into comparable biases in key cosmological parameters. Accurate modeling of electron‑ion non‑equipartition is therefore essential for extracting unbiased constraints on Ω_m, σ_8, and the dark‑energy equation of state from SZ‑selected cluster samples. The paper provides both a quantitative assessment of the bias and concrete, implementable methods to correct for it, underscoring the importance of plasma physics in the era of high‑precision cosmology.