Redshift-Space Distortion Measurement

• Redshift-space distortion (RSD) measurement is a technique that quantifies anisotropic imprints from galaxy peculiar velocities, enabling inference of the cosmic growth rate.
• Modern analyses employ sophisticated models, statistical estimators, and realistic mock catalogues to correct nonlinear effects and systematics like the Finger-of-God.
• By combining RSD with external probes and advanced modeling, researchers rigorously test gravitational theories and constrain key cosmological parameters.

Redshift-space distortion (RSD) measurement quantifies the anisotropic imprint of galaxy (or matter) peculiar velocities on observed large-scale structure statistics, enabling direct inference of the cosmic growth rate and rigorous tests of the underlying theory of gravity. RSDs manifest as scale- and angle-dependent deformations of clustering signals, arising from the mapping of galaxy positions from real to redshift space due to line-of-sight velocity components. Modern RSD analyses draw from a sophisticated suite of theoretical models, statistical estimators, and validation pipelines to extract unbiased constraints on structure growth and test gravitational physics from spectroscopic and, increasingly, photometric galaxy surveys.

1. Theoretical Foundations and Mathematical Framework

The origin of RSD lies in the fact that redshift-space positions are measured as

$$\mathbf{s} = \mathbf{r} + \frac{v_\parallel(\mathbf{r})}{aH} \hat{\mathbf{r}}$$

where $v_\parallel$ is the peculiar velocity component along the line-of-sight. The observed clustering is thus a convolution of the real-space density field and the large scale coherent flows and random virialized motions in biased tracer populations.

In linear theory (“Kaiser limit”), the Fourier-space galaxy power spectrum is modified to

$$P_g^s(k, \mu; z) = \left[b(z) + f(z) \mu^2\right]^2 P_m^r(k; z)$$

where $b(z)$ is the linear galaxy bias, $f(z) = d\ln D/d\ln a$ is the logarithmic growth rate, $\mu$ is the cosine of the angle between $\mathbf{k}$ and the line-of-sight, and $P_m^r(k; z)$ is the real-space matter power spectrum (Raccanelli et al., 2015, Howlett et al., 2014).

Nonlinear effects, notably the “Finger-of-God” (FoG) suppression due to small-scale velocity dispersions, are treated via phenomenological damping factors such as

$$D_{\text{FoG}}(k, \mu) = \exp\left[-(k\mu\sigma_v)^2\right] \quad \text{(Gaussian)} \ D_{\text{FoG}}(k, \mu) = \left[1 + (k\mu\sigma_v)^2\right]^{-1} \quad \text{(Lorentzian)}$$

with $\sigma_v$ the one-dimensional velocity dispersion (Raccanelli et al., 2015, Zheng et al., 2016).

Beyond linear theory, redshift-space statistics encode higher-order velocity-moment correlators, nonlinear density-velocity couplings, and bias stochasticity (Okumura et al., 2011, Howlett et al., 2014, Zhang et al., 2012). Model extensions include the TNS (Taruya–Nishimichi–Saito) approach and the Gaussian streaming model (GSM) combined with convolved Lagrangian perturbation theory (CLPT).

2. RSD Measurement Pipelines: Data, Estimators, and Fitting Procedures

RSD measurement requires carefully constructed data vectors and likelihoods:

• Sample selection: High-completeness spectroscopic samples (e.g., SDSS MGS: $z_{\mathrm{eff}}=0.15$, 63,163 galaxies; (Howlett et al., 2014)) are typically used, with well-defined radial and angular footprints and controlled selection functions.
• Mock catalogues: Realistic mock realizations are essential for error estimation and pipeline validation. For the MGS RSD measurement, 1000 mocks were synthesized via the PICOLA fast COLA N-body scheme (box size $1280\,h^{-1}\text{Mpc}$, $1536^3$ particles), populated with galaxies using a five-parameter HOD prescription tuned to match observed clustering (Howlett et al., 2014).
• Measurables: The two-point statistics are compressed into the configuration-space multipoles,

$$\xi_\ell(s) = \frac{2\ell+1}{2} \int_{-1}^{1} d\mu\, \xi(s, \mu) P_\ell(\mu)$$

with $\ell=0$ (monopole), $\ell=2$ (quadrupole), and sometimes $\ell=4$ (hexadecapole) (Howlett et al., 2014, Samushia et al., 2011).

• Theoretical modeling: The configuration-space correlation function is modeled with a GSM+CLPT framework, predicting $\xi(r)$, the mean infall velocity $v_{12}(r)$, and pairwise velocity dispersion $\sigma_{12}^2(r, \mu)$, properly incorporating scale-dependent bias and small-scale damping.
• Parameter inference: A Markov Chain Monte Carlo (MCMC) procedure (e.g., "emcee") is employed, varying $f\sigma_8$, $b\sigma_8$, higher-order bias, nuisance parameters for non-linear corrections, the BAO shape, the Alcock–Paczynski (AP) dilation and anisotropy parameters, and integral constraint offsets. Covariances are robustly estimated from mock samples using the Hartlap correction (Howlett et al., 2014).
• Range and binning: The RSD fits are stable with respect to choice of fitting range ($25<s<160\,h^{-1}\mathrm{Mpc}$), bin width ($\Delta s=8\,h^{-1}\mathrm{Mpc}$), and binning scheme. Results are tested for robustness to fitting/analysis choices (Howlett et al., 2014).

3. Statistical and Systematic Error Budget

The statistical precision of RSD-derived growth parameters is controlled by survey volume, tracer bias, and number density. For the Kaiser+exponential model on scales $3<r_p,\pi<35\,h^{-1}\mathrm{Mpc}$, statistical errors on $\beta=f/b$ are encapsulated by

$$\frac{\Delta\beta}{\beta} \simeq C\,b^{0.7} V^{-1/2} \exp\left[\frac{n_0}{b^2 n}\right]$$

with $C=4.9 \times 10^{2} h^{-1.5}\textrm{Mpc}^{1.5}$, $n_0=1.7\times10^{-4}h^3\textrm{Mpc}^{-3}$, $V$ the survey volume, $n$ the galaxy number density, and $b$ the bias (Bianchi et al., 2012). However, systematic errors—especially the systematic underestimation of $\beta$ due to inadequacies of the standard Kaiser+exponential-damping model even on quasi-linear scales—outstrip statistical errors, reaching $\sim 10\%$ or more unless nonlinear corrections (e.g., TNS, GSM) are included.

Systematic uncertainties are diagnosed and controlled using comparison to mocks, fit-variant ensembles, and detailed theory-model systematics validation. For the SDSS MGS at $z=0.15$, the primary result $f\sigma_8=0.49^{+0.15}_{-0.14}$ is unchanged within errors across algorithmic and modeling choices (Howlett et al., 2014).

4. Extensions and Applications: Gravity Testing, Alcock–Paczynski Effect, and Joint Constraints

RSD measurements, in conjunction with external geometric probes (CMB, BAO), enable precise tests of General Relativity (GR) and parameterizations of modified gravity via the growth index $\gamma$:

$f(z)=\Omega_m(z)^{\gamma}$

The SDSS MGS “free AP” RSD fit ($f\sigma_8=0.53\pm0.19$) can be combined with CMB and other BAO+RSD datasets to yield constraints $\gamma=0.64\pm0.09$, consistent at $1\sigma$ with the GR value $\gamma\approx0.55$, but with mild indication of a preference for weaker gravity (higher $\gamma$) (Howlett et al., 2014).

The inclusion of the AP effect—allowing for deviations in the fiducial cosmological distance scale—introduces additional degeneracies but can be robustly fit by marginalizing over isotropic dilation ($\alpha$) and anisotropy ($\epsilon$) parameters, with results validated on mocks and against BAO/cosmological priors (Howlett et al., 2014).

5. Model Development: Nonlinearities, Velocity Decomposition, and Advanced RSD Modeling

Contemporary research has established that percent-level RSD cosmology is unattainable without accounting for nonlinear evolution, bias, and multi-component velocity-field physics:

• Velocity decomposition: The total peculiar velocity is decomposed into three eigen-components: $v_\delta$ (correlated with density), $v_S$ (stochastic irrotational), and $v_B$ (rotational, post shell-crossing). Each contributes distinct terms to the redshift-space power spectrum and damping profiles (Zhang et al., 2012, Zheng et al., 2013).
• Streaming models and higher-order kernels: Accurate models, such as the GSM+CLPT adopted for SDSS MGS, incorporate scale-dependent Lagrangian bias, pairwise velocity statistics, and small-scale stochasticity, computed consistently to $O(P_\mathrm{lin}^2)$. These developments are essential for modeling the velocity-propagation and higher-point correlators that shape the RSD signal on relevant scales (Howlett et al., 2014).
• Multi-streaming effect: Random virial motions (“multi-streaming”) inside halos, especially in high-mass systems, introduce additional “Finger-of-God” damping that must be separated from coherent bulk flows for unbiased $f\sigma_8$ inference and breaking the $f\sigma_8$-$\sigma_v$ degeneracy in cosmological fits (Zheng et al., 2016).

6. Impact, Survey Prospects, and Future Directions

The anchoring of the $f\sigma_8(z)$ curve at low redshift ($z_\mathrm{eff}=0.15$ from SDSS MGS), together with measurements from surveys such as 6dFGS and BOSS LOWZ, provides a robust low-$z$ baseline for growth-of-structure tests and enables high-precision consistency checks against predictions from Planck CMB $\Lambda$CDM cosmology and GR (Howlett et al., 2014).

Prospects for future work include:

• Percent-level precision: Upcoming surveys (eBOSS, DESI) are expected to achieve percent-level precision on $f\sigma_8$ and $\gamma$, necessitating further theoretical advances in modeling nonlinear RSD and bias.
• Comprehensive error control: Continued use of large mock catalog suites, flexible bias parameterizations, robust AP marginalization, and improved velocity modeling are required to control both statistical and systematic uncertainties.
• Synergy with lensing and BAO: Joint lensing+RSD analyses will maximize information and break parameter degeneracies in $f$ and $b$, enabling unbiased extraction of structure growth in the low-$z$ universe.
• Constraints on fundamental physics: As statistical and systematics limits on $f\sigma_8$ and $\gamma$ tighten, RSD will continue to provide leading constraints on the behavior of gravity on cosmological scales, complementing geometric and weak lensing probes.

---

For a comprehensive technical reference on practical RSD measurement pipelines—including the MGS DR7 analysis, error budgets, and the impact of systematic effects—see (Howlett et al., 2014). For a systematic investigation of survey parameter dependence and scaling of RSD errors, see (Bianchi et al., 2012). For discussions of nonlinear corrections, advanced modeling, and the cosmological implications of RSD, see (Zhang et al., 2012, Taruya et al., 2013), and the references within.