Measuring the mass profile of galaxy clusters beyond their virial radius

Traditional estimators of the mass of galaxy clusters assume that the cluster components (galaxies, intracluster medium, and dark matter) are in dynamical equilibrium. Two additional estimators, that do not require this assumption, were proposed in the 1990s: gravitational lensing and the caustic technique. With these methods, we can measure the cluster mass within radii much larger than the virial radius. In the caustic technique, the mass measurement is only based on the celestial coordinates and redshifts of the galaxies in the cluster field of view; therefore, unlike lensing, it can be, in principle, applied to clusters at any redshift. Here, we review the origin, the basics and the performance of the caustic method.

💡 Research Summary

**
The paper reviews two mass‑estimation techniques for galaxy clusters that do not rely on the assumption of dynamical equilibrium: gravitational lensing and the caustic method. While lensing measures the deflection of background light and works well for low‑to‑moderate redshift clusters, it becomes observationally expensive at high redshift because sufficient background sources are required. The caustic technique, introduced by Diaferio and Geller in the early 1990s, overcomes this limitation by using only the projected positions and redshifts of galaxies in the cluster field.

The core idea is that galaxies in the infall region populate a characteristic “trumpet‑shaped” envelope in the phase‑space diagram (projected radius versus line‑of‑sight velocity). The envelope’s outer edge corresponds to the escape velocity at each radius. By measuring this edge, denoted (\mathcal{A}(r)), and integrating its square, one obtains the cumulative mass profile under the assumption of spherical symmetry:
(GM(<r)=\mathcal{F}{\beta}\int{0}^{r}\mathcal{A}^{2}(R),dR).
(\mathcal{F}_{\beta}) is a correction factor that accounts for the unknown orbital anisotropy (\beta); simulations suggest a typical value around 0.5, but the factor can be refined for individual clusters.

Implementation proceeds in four steps: (1) collect galaxy coordinates and redshifts, (2) construct a smoothed two‑dimensional density map in the ((R,\Delta v)) plane, (3) identify the caustic edge by tracing a chosen density contour, and (4) compute (\mathcal{A}(R)) and integrate to derive (M(<r)). Careful treatment of interlopers and adequate sampling in the outer regions are essential for reliable edge detection.

Performance tests using N‑body simulations show that the method recovers the true mass to within 20–30 % out to (3–5,r_{200}). Comparisons with X‑ray hydrostatic masses and weak‑lensing measurements for real clusters (e.g., CIRS, HeCS samples) reveal consistent results, confirming that the caustic technique can provide accurate mass estimates well beyond the virial radius. Its reliance solely on optical spectroscopy makes it especially valuable for high‑redshift clusters where lensing data are scarce.

Limitations are also discussed. Strong departures from spherical symmetry, ongoing major mergers, or the presence of massive substructures can distort the caustic envelope, leading to systematic biases. Sparse galaxy sampling reduces the robustness of the edge determination, and the assumption of a universal (\beta) may not hold for all clusters, necessitating cluster‑specific calibration of (\mathcal{F}_{\beta}). Recent advances aim to mitigate these issues through machine‑learning edge‑finding algorithms, joint analyses with X‑ray or Sunyaev‑Zel’dovich data, and Bayesian frameworks that incorporate simulation‑based priors on anisotropy.

The authors conclude that the caustic method is a powerful, observationally economical tool for probing cluster mass distributions out to several times the virial radius. When combined with forthcoming large‑scale spectroscopic surveys such as DESI and LSST, it promises to improve measurements of the cluster mass function and to tighten constraints on cosmological parameters, while also offering insights into the radial distribution of dark matter and the dynamics of cluster assembly.