The primary difficulty in measuring dynamical masses of galaxy clusters from galaxy data lies in the separation between true cluster members from interloping galaxies along the line of sight. We study the impact of membership contamination and incompleteness on cluster mass estimates obtained with 25 commonly used techniques applied to nearly 1000 mock clusters. We show that all methods overestimate or underestimate cluster masses when applied to contaminated or incomplete galaxy samples respectively. This appears to be the main source of the intrinsic scatter in the mass scaling relation. Applying corrections based on a prior knowledge of contamination and incompleteness can reduce the scatter to the level of shot noise expected for poorly sampled clusters. We establish an empirical model quantifying the effect of imperfect membership on cluster mass estimation and discuss its universal and method-dependent features. We find that both imperfect membership and the response of the mass estimators depend on cluster mass, effectively causing a flattening of the estimated - true mass relation. Imperfect membership thus alters cluster counts determined from spectroscopic surveys, hence the cosmological parameters that depend on such counts.

We use N-body simulations to quantify how the escape velocity in cluster-sized halos maps to the gravitational potential in a LambdaCDM universe. Using spherical density-potential pairs and the Poisson equation, we find that the matter density inferred gravitational potential profile predicts the escape velocity profile to within a few percent accuracy for group and cluster-sized halos (10^13 < M_200 < 10^15 M_sun, with respect to the critical density). The accuracy holds from just outside the core to beyond the virial radius. We show the importance of explicitly incorporating a cosmological constant when inferring the potential from the Poisson equation. We consider three density models and find that the Einasto and Gamma profiles provide a better joint estimate of the density and potential profiles than the Navarro, Frenk and White profile, which fails to accurately represent the escape velocity. For individual halos, the 1 sigma scatter between the measured escape velocity and the density-inferred potential profile is small (<5%). Finally, while the sub-halos show 15% biases in their representation of the particle velocity dispersion profile, the sub-halo escape velocity profile matches the dark matter escape velocity profile to high accuracy with no evidence for velocity bias outside 0.4r_200.

Ongoing and future spectroscopic surveys will measure numerous galaxy redshifts within tens of thousands of galaxy clusters. However, the sampling within these clusters will be low, 15 < N < 50 per cluster. With such data, it will be difficult to achieve accurate and precise mass estimates for individual clusters using phase-space mass estimation techniques. We develop and test a new stacking algorithm based on the caustic technique, which reduces the mass scatter in <log M_caustic | M_200 > for ensemble clusters from 70% for individual clusters to less than 10% for ensemble clusters with only 15 galaxies per cluster and 100 clusters per ensemble. With > 1000 galaxies per ensemble phase-space, the escape-velocity edge becomes readily identifiable and the presence of interloping galaxies is minimized. We develop and test an algorithm to trace the projected phase-space surface directly, which results in minimally biased dynamical mass estimates. We then quantify how binning and sampling affect the phase-space-based mass estimates when using an observational proxy that incorporates realistic mass scatter, like richness, and find the added uncertainty in the binning procedure has minimal influence on the resulting bias and scatter of the stacked mass estimates.

We derive the escape velocity profile for an Einasto density field in an accelerating universe and demonstrate its physical viability by comparing theoretical expectations to both light-cone data generated from N-body simulations and archival data on 20 galaxy clusters. We demonstrate that the projection function ($g(\beta )$) is deemed physically viable only for the theoretical expectation that includes a cosmology-dependent term. Using simulations, we show that the inferred velocity anisotropy is more than 6{\sigma} away from the expected value for the theoretical profile that ignores the acceleration of the universe. In the archival data, we constrain the average velocity anisotropy parameter of a sample of 20 clusters to be $\beta ={0.248}_{-0.360}^{+0.164}$ at the 68% confidence level. Lastly, we briefly discuss how our analytic model may be used as a novel cosmological probe based on galaxy clusters.

We measure the evolution of the velocity dispersion--temperature ($\sigma_{\rm v}$--$T_{\rm X}$) relation up to $z = 1$ using a sample of 38 galaxy clusters drawn from the \textit{XMM} Cluster Survey. This work improves upon previous studies by the use of a homogeneous cluster sample and in terms of the number of high redshift clusters included. We present here new redshift and velocity dispersion measurements for 12 $z > 0.5$ clusters observed with the GMOS instruments on the Gemini telescopes. Using an orthogonal regression method, we find that the slope of the relation is steeper than that expected if clusters were self-similar, and that the evolution of the normalisation is slightly negative, but not significantly different from zero ($\sigma_{\rm v} \propto T^{0.86 \pm 0.14} E(z)^{-0.37 \pm 0.33}$). We verify our results by applying our methods to cosmological hydrodynamical simulations. The lack of evolution seen in our data is consistent with simulations that include both feedback and radiative cooling.

Modified theories of gravity provide us with a unique opportunity to generate innovative tests of gravity. In Chameleon f(R) gravity, the gravitational potential differs from the weak-field limit of general relativity (GR) in a mass dependent way. We develop a probe of gravity which compares high mass clusters, where Chameleon effects are weak, to low mass clusters, where the effects can be strong. We utilize the escape velocity edges in the radius/velocity phase space to infer the gravitational potential profiles on scales of 0.3-1 virial radii. We show that the escape edges of low mass clusters are enhanced compared to GR, where the magnitude of the difference depends on the background field value |fR0|. We validate our probe using N-body simulations and simulated light cone galaxy data. For a DESI (Dark Energy Spectroscopic Instrument) Bright Galaxy Sample, including observational systematics, projection effects, and cosmic variance, our test can differentiate between GR and Chameleon f(R) gravity models, |fR0| = 4e-6 (2e-6) at > 5{\sigma} (> 2{\sigma}), more than an order of magnitude better than current cluster-scale constraints.

We present an initial study of the mass and evolutionary state of a massive and distant cluster, RCS2 J232727.6-020437. This cluster, at z=0.6986, is the richest cluster discovered in the RCS2 project. The mass measurements presented in this paper are derived from all possible mass proxies: X-ray measurements, weak-lensing shear, strong lensing, Sunyaev Zel'dovich effect decrement, the velocity distribution of cluster member galaxies, and galaxy richness. While each of these observables probe the mass of the cluster at a different radius, they all indicate that RCS2 J232727.6-020437 is among the most massive clusters at this redshift, with an estimated mass of M_200 ~3 x10^15 h^-1 Msun. In this paper, we demonstrate that the various observables are all reasonably consistent with each other to within their uncertainties. RCS2 J232727.6-020437 appears to be well relaxed -- with circular and concentric X-ray isophotes, with a cool core, and no indication of significant substructure in extensive galaxy velocity data.

This article is the second in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our aim is to quantify the scatter, systematic bias and completeness of cluster masses derived from a diverse set of 25 galaxy-based methods using two contrasting mock galaxy catalogues based on a sophisticated halo occupation model and a semi-analytic model. Analysing 968 clusters, we find a wide range in the RMS errors in log M200c delivered by the different methods (0.18 to 1.08 dex, i.e., a factor of ~1.5 to 12), with abundance matching and richness methods providing the best results, irrespective of the input model assumptions. In addition, certain methods produce a significant number of catastrophic cases where the mass is under- or over-estimated by a factor greater than 10. Given the steeply falling high-mass end of the cluster mass function, we recommend that richness or abundance matching-based methods are used in conjunction with these methods as a sanity check for studies selecting high mass clusters. We see a stronger correlation of the recovered to input number of galaxies for both catalogues in comparison with the group/cluster mass, however, this does not guarantee that the correct member galaxies are being selected. We do not observe significantly higher scatter for either mock galaxy catalogues. Our results have implications for cosmological analyses that utilise the masses, richnesses, or abundances of clusters, which have different uncertainties when different methods are used.

This paper is the first in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our primary aim is to test the performance of these cluster mass estimation techniques on a diverse set of models that will increase in complexity. We begin by providing participating methods with data from a simple model that delivers idealised clusters, enabling us to quantify the underlying scatter intrinsic to these mass estimation techniques. The mock catalogue is based on a Halo Occupation Distribution (HOD) model that assumes spherical Navarro, Frenk and White (NFW) haloes truncated at R_200, with no substructure nor colour segregation, and with isotropic, isothermal Maxwellian velocities. We find that, above 10^14 M_solar, recovered cluster masses are correlated with the true underlying cluster mass with an intrinsic scatter of typically a factor of two. Below 10^14 M_solar, the scatter rises as the number of member galaxies drops and rapidly approaches an order of magnitude. We find that richness-based methods deliver the lowest scatter, but it is not clear whether such accuracy may simply be the result of using an over-simplistic model to populate the galaxies in their haloes. Even when given the true cluster membership, large scatter is observed for the majority non-richness-based approaches, suggesting that mass reconstruction with a low number of dynamical tracers is inherently problematic.

We quantify the expected observed statistical and systematic uncertainties of the escape velocity as a measure of the gravitational potential and total mass of galaxy clusters. We focus our attention on low redshift (z < 0.15) clusters, where large and deep spectroscopic datasets currently exist. Utilizing a suite of Millennium Simulation semi-analytic galaxy catalogs, we find that the dynamical mass, as traced by either the virial relation or the escape velocity, is robust to variations in how dynamical friction is applied to "orphan" galaxies in the mock catalogs (i.e., those galaxies whose dark matter halos have fallen below the resolution limit). We find that the caustic technique recovers the known halo masses (M_200) with a third less scatter compared to the virial masses. The bias we measure increases quickly as the number of galaxies used decreases. For N_gal > 25, the scatter in the escape velocity mass is dominated by projections along the line-of-sight. Algorithmic uncertainties from the determination of the projected escape velocity profile are negligible. We quantify how target selection based on magnitude, color, and projected radial separation can induce small additional biases into the escape velocity masses. Using N_gal = 150 (25), the caustic technique has a per cluster scatter in ln(M|M_200) of 0.3 (0.5) and bias 1+/-3% (16+/-5%) for clusters with masses > 10^14M_solar at z<0.15.

We use the Millennium Simulation to quantify the statistical accuracy and precision of the escape velocity technique for measuring cluster-sized halo masses at z~0.1. We show that in 3D, one can measure nearly unbiased (<4%) halo masses (>1.5x10^14 M_solar h^-1) with 10-15% scatter. Line-of-sight projection effects increase the scatter to ~25%, where we include the known velocity anisotropies. The classical "caustic" technique incorporates a calibration factor which is determined from N-body simulations. We derive and test a new implementation which eliminates the need for calibration and utilizes only the observables: the galaxy velocities with respect to the cluster mean v, the projected positions r_p, an estimate of the Navarro-Frenk-White (NFW) density concentration and an estimate of the velocity anisotropies, beta. We find that differences between the potential and density NFW concentrations induce a 10% bias in the caustic masses. We also find that large (100%) systematic errors in the observed ensemble average velocity anisotropies and concentrations translate to small (5%-10%) biases in the inferred masses.