Sunyaev–Zeldovich y Parameter Calculator
The Thermal Sunyaev–Zeldovich Effect in Galaxy Clusters
The Sunyaev–Zeldovich (SZ) effect is one of the clearest examples of how tiny interactions can leave measurable signatures across cosmological distances. It appears when cosmic microwave background photons pass through hot ionized gas, especially the intracluster medium inside galaxy clusters, and are inverse-Compton scattered by energetic electrons. Each scattering transfers only a small amount of energy, but over a long path through a dense, hot plasma the cumulative effect becomes observable. The standard quantity used to describe this distortion is the dimensionless Compton y parameter. In practical terms, y measures the integrated electron pressure along the line of sight.
This calculator is designed as a quick-look tool for the thermal SZ effect in a simple uniform slab model. You enter an electron density, an electron temperature, and a path length. The page then computes three related outputs: the Compton y parameter, the Thomson optical depth τ, and the approximate Rayleigh–Jeans brightness temperature shift of the cosmic microwave background. That makes the tool useful for students learning the physics, for astronomy readers checking order-of-magnitude estimates, and for anyone who wants to connect cluster gas properties to an observable CMB signal.
Introduction
In a hot cluster atmosphere, electrons typically have temperatures of several keV, corresponding to tens of millions of kelvin. CMB photons are much lower in energy, so when they scatter from these electrons they gain energy on average. The result is not simply a dimming or brightening of the background. Instead, the CMB spectrum is distorted: at low frequencies there is a decrement, while at sufficiently high frequencies there is an increment. In the low-frequency Rayleigh–Jeans limit, the temperature change is approximately proportional to −2y. That simple relation is what this calculator uses to report an estimated temperature decrement.
The importance of the SZ effect in astronomy comes from its unusual observational behavior. Ordinary surface brightness fades strongly with distance, but the SZ surface brightness does not suffer the same redshift dimming in the usual way. Because of that, distant clusters can still produce detectable SZ signals. This is one reason why SZ surveys have become powerful tools for finding galaxy clusters across a wide range of cosmic time. The same physics also makes the effect valuable for studying pressure profiles, cluster energetics, feedback, and cosmological structure growth.
For a slab of plasma with uniform electron density ne, temperature Te, and physical depth L, the classical expression for y is
This compact expression contains the main physical ingredients. The Thomson cross section σT sets the scattering probability for low-energy photons. The factor kBTe represents the thermal energy scale of the electrons. The denominator mec² normalizes that thermal energy to the electron rest energy, giving the characteristic fractional energy transfer per scattering in the non-relativistic limit. Multiplying by density and path length turns the local plasma properties into a line-of-sight integral. In other words, the thermal SZ effect traces pressure rather than just density.
How to Use
Using the calculator is straightforward, but it helps to know exactly what each field means. The first input is the electron density ne in cm⁻³. This is the number density of free electrons in the plasma. In cluster outskirts it may be around 10⁻⁴ cm⁻³, while in denser central regions it can be closer to 10⁻³ or 10⁻² cm⁻³. The second input is the electron temperature Te in keV. X-ray astronomers often quote temperatures in keV because that unit directly reflects the thermal energy scale. Typical cluster gas temperatures range from a few keV to above 10 keV in very massive or merging systems. The third input is the path length L in kiloparsecs, representing the effective thickness of the plasma along the line of sight.
After entering values, select the compute button. The script converts the density from cm⁻³ to m⁻³, converts the temperature from keV into joules for the quantity kBT, and converts the path length from kiloparsecs to meters. It then evaluates the thermal SZ expression and displays the results in scientific notation where appropriate. Because the model is linear in density, temperature, and path length, the output is easy to scale mentally. If you double the density, both y and τ double. If you double the temperature, y doubles but τ does not change. If you double the path length, both quantities double again.
When interpreting the result, remember that y is dimensionless and usually small. Values around 10⁻⁶ to 10⁻⁴ are common for diffuse or moderate cluster gas, while larger values can occur in dense, hot regions. The optical depth τ is also dimensionless and is usually much less than one, which is why the single-scattering approximation is generally appropriate. The reported ΔT value is the approximate Rayleigh–Jeans temperature shift in kelvin. A negative value means a decrement in the low-frequency CMB brightness temperature.
Formula
The numerical implementation follows the standard non-relativistic thermal SZ approximation. The code computes
and also computes the Thomson optical depth
followed by the Rayleigh–Jeans temperature relation
with TCMB taken to be 2.725 K. This means the displayed temperature shift is
The formula is intentionally simple. It assumes a uniform plasma slab, non-relativistic thermal electrons, and the Rayleigh–Jeans limit for the observed frequency. In real observations, the full SZ spectrum depends on frequency, and the gas properties vary with radius and substructure. Even so, this approximation is excellent for building intuition and for checking whether a chosen set of cluster parameters gives a plausible signal level.
The small spectral distortion caused by a non-relativistic population of electrons produces a decrement in brightness temperature at radio and millimeter wavelengths and an increment in the Wien tail. In the Rayleigh–Jeans regime, the fractional temperature shift is approximately ΔT/TCMB = −2y. Our calculator adopts this approximation and reports both y and the implied ΔT for the canonical TCMB = 2.725 K. It also returns the optical depth τ = neσTL to help users gauge how many scatterings typically occur along the path.
The numerical implementation multiplies the user-supplied density (converted from cm⁻³ to m⁻³), temperature (keV to joules), and path length (kpc to meters). Together with fundamental constants, these inputs produce y, τ, and the Rayleigh–Jeans temperature decrement. The results scale linearly with each parameter, so doubling the density doubles both y and τ, whereas doubling the temperature doubles y but leaves τ unchanged because the optical depth does not depend on thermal energy.
Example
A useful worked example is a cluster atmosphere with electron density 1×10⁻³ cm⁻³, electron temperature 5 keV, and path length 500 kpc. These are reasonable order-of-magnitude values for hot intracluster gas. Entering those numbers into the calculator gives a Compton y of about 3.2×10⁻⁵. The corresponding optical depth is small, confirming that the plasma is optically thin to Thomson scattering. The Rayleigh–Jeans temperature shift is about −0.00017 K, or roughly −0.17 mK. That may sound tiny, but modern CMB instruments are designed to detect signals at this level.
This example also shows how to think about scaling. If the same gas were twice as hot, the optical depth would stay the same because the number of electrons and the path length would be unchanged, but the Compton y value would double because the pressure doubled. If instead the density were five times larger while the temperature stayed fixed, both y and τ would increase by a factor of five. These simple proportionalities are one reason the thermal SZ effect is so useful as a pressure probe.
To illustrate typical numbers encountered in cluster astrophysics, the table below provides examples for idealized plasma slabs. These demonstrate how modest variations in electron pressure modify the y parameter by orders of magnitude:
| ne (cm⁻³) | Te (keV) | L (kpc) | y | ΔT (mK) |
|---|---|---|---|---|
| 1×10⁻³ | 5 | 500 | 3.2×10⁻⁵ | −0.17 |
| 5×10⁻³ | 8 | 1000 | 2.1×10⁻⁴ | −1.14 |
| 1×10⁻⁴ | 10 | 1500 | 8.0×10⁻⁶ | −0.04 |
Though the table lists idealized uniform slabs, real clusters exhibit radial gradients and substructure. Observers often combine resolved SZ maps with X-ray measurements to infer three-dimensional pressure distributions. Those reconstructions support studies of cluster scaling relations, baryon content, feedback from active galactic nuclei, and the thermodynamic history of the intracluster medium. Because the SZ signal remains observable at high redshift, surveys such as ACT, SPT, and Planck have discovered large cluster samples that are valuable for cosmology.
Limitations and Assumptions
This calculator intentionally uses a simplified physical model. The first limitation is the assumption of a uniform slab. Real galaxy clusters are not uniform: density and temperature vary with radius, shocks can heat localized regions, and mergers can create strong asymmetries. A single density, temperature, and path length therefore represent an effective average rather than a full physical description. The result is best interpreted as an order-of-magnitude estimate or a pedagogical approximation.
The second limitation is the use of the non-relativistic thermal SZ formula. Relativistic corrections become increasingly important when the electron temperature reaches the high-keV regime, especially in very hot or merging clusters. Those corrections alter the detailed frequency dependence of the SZ spectrum and can shift the inferred signal if high precision is required. This page does not include those refinements, so users doing precision analysis should treat the output as a baseline estimate rather than a final observational model.
A third limitation is that the displayed temperature shift uses the Rayleigh–Jeans approximation. That relation is appropriate at low observing frequencies, but the full thermal SZ effect is frequency dependent and crosses through a null near 217 GHz in the non-relativistic limit. Above that frequency the sign of the thermal distortion changes. Therefore, the reported ΔT should be read specifically as a low-frequency brightness-temperature approximation, not as a universal temperature shift at every observing band.
There are also physical effects beyond the thermal SZ signal. The kinetic SZ effect arises from the bulk motion of the plasma relative to the CMB rest frame and depends on line-of-sight velocity rather than temperature. Non-thermal electron populations, magnetic fields, clumping, and projection effects can further complicate interpretation. None of those are included here. Even with those caveats, the calculator remains useful because it captures the leading dependence of the thermal SZ signal on electron pressure and optical depth.
Relativistic corrections become non-negligible when kBTe approaches tens of keV. They modify the spectral shape of the SZ effect and must be accounted for in high-precision work, particularly for merging clusters with shock-heated gas. Our calculator does not include these corrections, but the output y still conveys the integrated pressure regardless of electron rest frame velocity, making it a useful quick-look diagnostic.
Beyond clusters, the SZ effect reveals the broader cosmic web. The warm-hot intergalactic medium predicted by cosmological simulations should produce faint y signatures at the level of 10⁻⁷–10⁻⁶. Detecting this diffuse component requires stacking observations of many large-scale structures or performing cross-correlations with galaxy surveys. Future experiments like CMB-S4 and space missions dedicated to spectral distortions may map this low-level signal, helping to resolve the long-standing missing baryon problem.
From a theoretical perspective, the SZ effect can be derived by solving the Kompaneets equation, a Fokker-Planck approximation to the Boltzmann equation for Compton scattering. The differential form of the Kompaneets equation is
where x = hν / (kBTe) is the dimensionless photon frequency and n is the occupation number. Integrating this equation along a line of sight under the assumption of small energy exchange per scattering yields the y parameter. Although solving the full equation is beyond the scope of this web tool, mentioning it underscores the deep kinetic theory underpinning the deceptively simple y formula.
Importantly, the SZ effect allows direct measurements of the Hubble constant when combined with X-ray data. The combination yields the angular diameter distance to a cluster without relying on intermediate standard candles. Such geometric determinations help resolve tensions among various cosmological measurements. Additionally, the kinetic SZ effect—arising from the bulk motion of clusters relative to the CMB rest frame—offers a way to map large-scale velocity fields, potentially constraining dark energy and modifications to gravity.
Technological advances are rapidly improving SZ observations. Multi-band imagers can separate the thermal, kinetic, and relativistic components of the effect. Interferometers deliver arcsecond resolution, revealing subcluster merger shocks and turbulence. Polarization measurements, though challenging, could someday map transverse cluster motions. As data quality rises, accurate modeling of relativistic corrections and feedback processes becomes essential to avoid biases in cosmological inference.
The SZ effect also has applications in fundamental physics. Because it probes the hot electron content of clusters independent of redshift, comparisons between SZ and gravitational lensing masses test the validity of hydrostatic equilibrium and thereby inform theories of modified gravity. Constraints on the integrated y signal across the sky, often called the Compton y background, set limits on energy injection from early universe processes such as primordial black hole evaporation or dark matter annihilation.
In summary, the Sunyaev–Zeldovich effect transforms the cosmic microwave background into a backlight for the hot universe. By entering densities, temperatures, and path lengths into this calculator, you can quickly estimate the Compton y parameter, optical depth, and resulting low-frequency brightness temperature decrement. These quantities are small, but they connect directly to cluster pressure, scattering probability, and observable CMB distortions. That makes the calculator a practical bridge between basic plasma properties and one of the most important probes of hot gas in modern cosmology.