The redshift-space cluster-galaxy cross-correlation function - I. Modelling galaxy infall on to Millennium simulation clusters and SDSS groups

The large-scale infall of galaxies around massive clusters provides a potentially powerful diagnostic of structure growth, dark energy and cosmological deviations from General Relativity. We develop and test a method to recover galaxy infall kinematics (GIK) from measurements of the redshift-space cluster-galaxy cross-correlation function ξ ^s_cg(r_p,r_π ). Using galaxy and halo samples from the Millennium simulation, we calibrate an analytic model of the galaxy kinematic profiles comprising a virialized component with an isotropic Gaussian velocity distribution and an infall component described by a skewed 2D t-distribution with a characteristic infall velocity v_{r, c} and separate radial and tangential dispersions. We show that convolving the real-space cross-correlation function with this velocity distribution accurately predicts the redshift-space ξ ^s_cg, and we show that measurements of ξ ^s_cg can be inverted to recover the four distinct elements of the GIK profiles. These in turn provide diagnostics of cluster mass profiles, and we expect the characteristic infall velocity v_{r, c}(r) in particular to be insensitive to galaxy formation physics that can affect velocity dispersions within haloes. As a proof of concept we measure ξ ^s_cg for rich galaxy groups in the Sloan Digital Sky Survey and recover GIK profiles for groups in two bins of central galaxy stellar mass. The higher mass bin has a v_{r, c}(r) curve very similar to that of 10¹⁴ h^-1 M_⊙ haloes in the Millennium simulation, and the recovered kinematics follow the expected trends with mass. GIK modelling of cluster-galaxy cross-correlations can be a valuable complement to stacked weak lensing analyses, allowing novel tests of modified gravity theories that seek to explain cosmic acceleration.