Redshift-Space Distortions in Cosmology

• Redshift-space distortions (RSD) are anisotropies in observed galaxy clustering caused by line-of-sight peculiar velocities, linking real-space to redshift-space measurements.
• Advanced modeling techniques, from linear Kaiser effects to non-perturbative streaming models, capture both large-scale coherent infall and small-scale damping (Fingers-of-God).
• RSD analyses provide actionable constraints on cosmic structure growth, galaxy-halo connection, and tests of modified gravity, crucial for next-generation cosmological surveys.

Redshift-space distortions (RSD) are observational anisotropies in large-scale structure and astrophysical signals caused by the mapping between physical positions and observed redshifts in the presence of line-of-sight peculiar velocities. These distortions fundamentally alter clustering statistics by introducing correlations along the line of sight, providing a unique probe of the growth of cosmic structure, the galaxy–halo connection, and tests of gravity and fundamental physics. Due to their sensitivity to coherent infall and small-scale virialized motions, RSD have become a central cosmological observable in galaxy and intensity mapping surveys, as well as in 21-cm experiments.

1. Theoretical Foundations of RSD

In the standard plane-parallel approximation, the observed redshift-space position $\mathbf{s}$ of a tracer is related to its real-space position $\mathbf{r}$ and its line-of-sight peculiar velocity $v_\parallel$ via

$$\mathbf{s} = \mathbf{r} + \frac{v_\parallel}{aH} \, \hat{\mathbf{n}}$$

where $a$ is the cosmological scale factor and $H$ is the Hubble parameter (Raccanelli et al., 2015). The associated effect on the overdensity field leads, even in the linear theory with a uniform bias $b$, to a characteristic anisotropy in clustering. The Fourier-space galaxy power spectrum becomes

$P^s_g(k,\mu) = [b + f \mu^2]^2 P_m(k)$

with $f\equiv d\ln D/d\ln a$ the logarithmic growth rate, $D(a)$ the linear growth factor, $P_m(k)$ the real-space matter power spectrum, and $\mu = \cos\theta$ the cosine of the angle between the wavevector and the line of sight (Raccanelli et al., 2015, Seljak et al., 2011). This is known as the Kaiser effect and dominates large, quasi-linear scales.

On small (nonlinear) scales, random galaxy motions within halos introduce additional line-of-sight smearing ("Fingers-of-God" or FoG), typically modeled phenomenologically as a Gaussian or Lorentzian velocity dispersion: $$P^s_{\text{nl}}(k,\mu) = P^s_{\text{lin}}(k,\mu) \exp[-(k\mu\sigma_v)^2]$$ with $\sigma_v$ the 1D velocity dispersion (Raccanelli et al., 2015).

RSD also appear in other observables, e.g., the 21-cm intensity mapping power spectrum or velocity fields, where analogous anisotropic enhancements and damping are present (Ross et al., 2020, Dam et al., 2021).

2. Advanced Modelling and Moment Expansions

Beyond the linear regime, RSD require modeling the full mapping of the real-space density and velocity fields into redshift space. The exact transformation can be expanded as

$$\delta_s(\mathbf{k}) = \sum_{L=0}^{\infty}\frac{1}{L!} \left(\frac{ik\mu}{H}\right)^L T^L_{\parallel}(\mathbf{k})$$

with $T^L_{\parallel}(\mathbf{k})$ the Fourier transform of the density-weighted $L$-th moment of the line-of-sight velocity (Seljak et al., 2011, Okumura et al., 2011, Vlah et al., 2012). The redshift-space power spectrum is then written as a double sum over velocity moment correlators

$$P^s(k,\mu) = \sum_{L,L'} \frac{(-1)^{L'}}{L!L'!} \left( \frac{ik\mu}{H} \right)^{L+L'} P_{LL'}(k)$$

where $P_{LL'}(k)$ are the auto- and cross-power spectra of the moments.

The expansion can be grouped into powers of $\mu^2$, with known terms up to $\mathcal{O}(\mu^4)$ carrying the dominant contributions on large and intermediate scales. On large scales, only a few low-order terms contribute (density–momentum, momentum–momentum), but at $k \gtrsim 0.1\,h\,\text{Mpc}^{-1}$, higher moments and non-Gaussian FoG-like corrections (from quadratic and higher velocity moments) become significant (Okumura et al., 2011, Vlah et al., 2012, Vlah et al., 2013).

These techniques also allow for a self-consistent treatment of scale-dependent velocity bias in galaxies and halos, the treatment of nonlinear bias models, and the calculation of all relevant multipoles (Vlah et al., 2013).

3. Statistical Measures and Observables

RSD induce characteristic angular and radial dependencies in statistics such as:

• Power spectrum multipoles ($P_\ell(k)$, $\ell=0,2,4$): decomposing $P^s(k,\mu)$ into Legendre series (Raccanelli et al., 2015, Collaboration et al., 28 Jan 2026).
• Configuration-space correlation function multipoles ($\xi_\ell(s)$): used for fitting the monopole and quadrupole components (Hernández-Aguayo et al., 2018).
• Projected correlation function $w(r_p)$: integrating out the redshift distortion along the line of sight gives access to the real-space clustering amplitude (Tonegawa et al., 2020).
• "Density-split" correlations: cross-correlating densities in different local environments (e.g., by quantiles) recovers non-Gaussian clustering information and tightens cosmological constraints (Paillas et al., 2021).
• 21-cm anisotropy multipoles: RSD in the 21-cm power spectrum are captured by enhancement of the quadrupole and hexadecapole during the Cosmic Dawn (Ross et al., 2020).

In angular (photometric) surveys, the impact of RSD is encoded in the angular correlation function $w(\theta)$, which for top-hat redshift bins retains RSD information, while specialized pair-center binning can suppress the distortion for validation and bias-breaking (Ross et al., 2011).

A notable advanced analysis is the use of optimally redshift-weighted statistics, where galaxies are assigned redshift-dependent weights tailored to minimize the variance for cosmological parameters, compressing the full redshift evolution of growth and geometry in a small number of weighted spectra and typically yielding better precision than redshift binning (Ruggeri et al., 2016).

4. Nonlinearities, Fingers-of-God, and Modeling Advances

Precise RSD modeling at the percent level for next-generation surveys requires capturing:

• Nonlinear evolution: Standard perturbation theory (SPT) and Effective Field Theory (EFT) approaches extend predictions beyond the linear regime, with correction kernels and bias expansions (Song et al., 2018, Eggemeier et al., 30 Jan 2025, Collaboration et al., 28 Jan 2026, Vlah et al., 2012).
• Non-perturbative streaming/FoG effects: Resummed moment expansion or streaming models, such as the Velocity Difference Generator (VDG) analytic resummation, capture virialized random motions non-perturbatively, facilitating unbiased parameter estimation to significantly smaller scales than EFT (Eggemeier et al., 30 Jan 2025, Collaboration et al., 28 Jan 2026).
• Hybrid/Emulator approaches: N-body simulation-calibrated emulators such as BACCO combine Lagrangian bias theory with N-body measured corrections, reaching high precision at intermediate and small scales (Collaboration et al., 28 Jan 2026).
• Wide-angle corrections: At large pair separations, the plane-parallel approximation fails, demanding modeling using tripolar spherical harmonics or pair-based weighting to reconstruct the fully non-parallel RSD signal (Raccanelli et al., 2010).

Table: Comparison of Modern RSD Modeling Approaches (Editor’s term)

Model Type	Valid $k$ Range [$h$/Mpc]	Key Nonlinear Feature
EFT (1-loop)	$\lesssim 0.20-0.25$	Perturbative corrections and counterterms (Eggemeier et al., 30 Jan 2025, Collaboration et al., 28 Jan 2026)
VDG/Streaming	$\lesssim 0.35-0.40$	Non-perturbative FoG kernel resummation (Eggemeier et al., 30 Jan 2025, Collaboration et al., 28 Jan 2026)
BACCO Emulator	$\lesssim 0.30$	Simulation-calibrated Lagrangian bias and nonlinear mapping (Collaboration et al., 28 Jan 2026)

Joint power spectrum and bispectrum fits incorporating the non-perturbative streaming model can deliver up to a factor two improvement in cosmological parameter uncertainties relative to EFT over the same $k$-range (Eggemeier et al., 30 Jan 2025).

5. Cosmological Information, Degeneracies, and Parameter Forecasts

RSD are uniquely sensitive to the rate of structure growth, parameterized through $f$ or $f\sigma_8$. In surveys, care is required to break degeneracies between galaxy bias $b$ and $f$. Techniques include:

• Top-hat vs. pair-center binning in $w(\theta)$: measuring $b\sigma_8$ from pair-center bins, fitting $f\sigma_8$ from top-hat bins (Ross et al., 2011).
• Nonlinear weakly-perturbative modeling: including higher-order terms in the redshift-space power spectrum and bias expansions breaks the $f-\sigma_8$ degeneracy, allowing constraints on each at few percent precision in Stage-IV spectroscopic surveys (Song et al., 2021).
• Multi-tracer analyses: Joint fits to density and velocity spectra as in peculiar-velocity surveys further break degeneracies (Dam et al., 2021).
• Optimal redshift weighting: Redshift-dependent galaxy weights compress redshift evolution while retaining optimal sensitivity to deviations in $f(z)$ or $\Omega_m(z)$ (Ruggeri et al., 2016).

Justifications, forecasts, and parameter errors:

• For next-generation photometric surveys (e.g., DES), RSD measurements enable $f(z)\sigma_8(z)$ determinations to 17% per narrow redshift slice, and growth index $\gamma$ to 25% precision (Ross et al., 2011).
• In SKA and Euclid-like spectroscopic surveys, sub-percent level constraints on $f\sigma_8$, $H(z)$, $D_A(z)$ are achievable over $z \sim 0.7$–1.8, contingent on precise RSD modeling (Raccanelli et al., 2015, Collaboration et al., 28 Jan 2026).

6. RSD as Probes of Gravity, Fundamental Physics, and Systematics

RSD are sensitive to modifications of general relativity, as the growth rate $f$ depends on the theory of gravity:

• Modified gravity: Different screening mechanisms in $f(R)$ and nDGP models lead to small but distinguishable deviations in RSD observables, especially in the quadrupole and monopole of the redshift-space correlation function on large scales (Hernández-Aguayo et al., 2018, Taruya et al., 2013).
• Velocity bias and galaxy–halo connection: Small-scale RSD, including FoG anisotropy, provide sensitivity to satellite galaxy velocity bias, merger time-scale parameters, and baryonic physics inside halos. Empirically, a satellite velocity bias $b_v^s \lesssim 1$ relative to MBPs is found, with gravitational dynamics (rather than baryonic effects) dominating the suppression of velocity dispersion (Tonegawa et al., 2020).
• Wide-angle and selection effects: For surveys with large solid-angle coverage, wide-angle and mode-coupling corrections must be explicitly modeled, as neglect leads to systematic biases on BAO scales at the percent level (Raccanelli et al., 2010).
• Nonlinear iterative reconstruction: Nonperturbative, iterative real-to-redshift-space reconstructions allow direct removal of RSD from galaxy catalogs, accurately recovering real-space clustering and BAO features and offering an internal calibration of $f/b$ (Wang et al., 2019).

7. RSD in 21-cm Cosmology and Alternative Tracers

In intensity mapping of the 21-cm line, RSD enhance the signal power (Kaiser effect), induce strong anisotropy, and modulate the detectability of higher-order statistics across reionization and Cosmic Dawn astrophysical regimes (Ross et al., 2020). During inhomogeneous X-ray heating, RSD-induced anisotropy is suppressed, while at earlier and later cosmic epochs they provide a nearly scale-independent boost to the power spectrum by a factor ≃1.87 (linear expectation).

In summary, redshift-space distortions encode rich information about cosmic structure growth, gravity, and astrophysical processes. Their accurate modeling—spanning advanced perturbation theory, velocity-moment expansions, nonperturbative resummation, hybrid emulation, and iterative reconstruction—is essential for extracting precision cosmological constraints from current and future large-scale structure surveys (Eggemeier et al., 30 Jan 2025, Collaboration et al., 28 Jan 2026, Ross et al., 2011).