
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 phasespace 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 phasespace, the escapevelocity edge becomes readily identifiable and
the presence of interloping galaxies is minimized. We develop and test an
algorithm to trace the projected phasespace surface directly, which results in
minimally biased dynamical mass estimates. We then quantify how binning and
sampling affect the phasespacebased 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 use Nbody simulations to quantify how the escape velocity in
clustersized halos maps to the gravitational potential in a LambdaCDM
universe. Using spherical densitypotential 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
clustersized 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 densityinferred potential profile is small (<5%). Finally,
while the subhalos show 15% biases in their representation of the particle
velocity dispersion profile, the subhalo escape velocity profile matches the
dark matter escape velocity profile to high accuracy with no evidence for
velocity bias outside 0.4r_200.

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 lightcone data generated from Nbody
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 cosmologydependent 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 dispersiontemperature
($\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 selfsimilar, 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 weakfield 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.31 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 Nbody
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 = 4e6 (2e6)
at > 5{\sigma} (> 2{\sigma}), more than an order of magnitude better than
current clusterscale constraints.

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 semianalytic 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 lineofsight. 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(MM_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 clustersized 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 1015% scatter. Lineofsight 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 Nbody 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 NavarroFrenkWhite (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.