---
title: Void–Galaxy Cross-Correlation Function
url: https://www.emergentmind.com/topics/void-galaxy-cross-correlation-function
type: topic
---

# Void–Galaxy Cross-Correlation Function

The void–galaxy cross-correlation function is the two-point statistic between cosmic-void centres and galaxies, defined as the excess or deficit of galaxies at a given void-centric separation relative to a random distribution. In configuration space it is equivalent to the stacked radial density profile of tracers around void centres; in redshift space it becomes anisotropic and encodes coherent outflows, the Alcock–Paczynski effect, and higher-order relativistic contributions. It is used as a diagnostic of void environments in SDSS, as a growth-rate test in 6dF, as a high-precision AP/RSD observable in BOSS, and as a forecasted Euclid and DESI probe in recent survey analyses [1606.03092] [1904.01030] [2302.05302] [2504.08221].

## 1. Definitions and mathematical representations

In its most direct form, the void–galaxy cross-correlation function is the tracer density contrast around void centres,
\[
\xi_{vg}(r) \equiv \frac{\rho_g(r)}{\bar{\rho}_g}-1 \equiv \delta_g(r),
\]
or, more generally for a tracer population \(T\),
\[
1+\xi_{vT}(r)=\frac{\langle n_T(r)\rangle}{\bar n_T}.
\]
This identifies the cross-correlation with the stacked density profile around voids. The generic profile recalled in void analyses is a strong underdensity near the centre, a rise toward the void edge, an overdense ridge or wall at the boundary, and an asymptotic approach to the mean density at larger radii [2009.14751].

| Representation | Expression | Typical use |
|---|---|---|
| 3D isotropic CCF | \(\xi_{vg}(r)\) | Stacked tracer profile around void centres |
| Anisotropic CCF | \(\xi_{vg}(r_p,\pi)\) or \(\xi^s(s,\mu)\) | RSD and AP analyses |
| Projected CCF | \(w_{p,vg}(r_p)=2\int_0^{\Pi_{\max}}\xi_{vg}(r_p,\pi)\,d\pi\) | Real-space or weak-RSD compressed observable |
| Angular tomographic form | \(C^{vg}_{ij}(\ell)\) | Photometric and harmonic-space analyses |

The projected statistic used in forward-model DESI forecasts is
\[
w_{p,vg}(r_p)=2\int_0^{\Pi_{\max}}\xi_{vg}(r_p,\pi)\,d\pi,
\]
with \(\Pi_{\max}=200\,h^{-1}\,\mathrm{Mpc}\) in that analysis, and with forecast fits restricted to \(25<r_p<150\,h^{-1}\,\mathrm{Mpc}\) to reduce sensitivity to noisy derivatives and modelling challenges [2504.08221]. On large scales, the Fourier-space cross-spectrum is commonly written as
\[
P_{vg}(k,z)\simeq b_v(z)\,b_g(z)\,P_{mm}(k,z),
\]
with \(b_v\) and \(b_g\) the void and galaxy biases [2206.14211].

A tomographic angular counterpart follows the same projection structure used for void–lensing. Replacing the lensing kernel by a galaxy tracer kernel yields
\[
C^{vg}_{ij}(\ell)=\frac{c}{H_0}\int dz\,
\frac{W_i^v(z)\,W_j^g(z)}{E(z)\,r^2(z)}
\,\hat P_{mm}\!\left(\frac{\ell+1/2}{r(z)},z\right),
\]
with \(W_j^g(z)=\frac{H(z)}{c}n_j^g(z)b_g(z)\). This is mathematically parallel to the void–lensing cross-spectrum and makes explicit that the void–galaxy CCF is one instance of a broader void–tracer cross-correlation framework [2206.14211].

## 2. Redshift-space formulation and multipole structure

The central observable in spectroscopic analyses is the anisotropic redshift-space cross-correlation. A standard streaming description writes
\[
1+\xi^{\rm rs}(\mathbf{s})=\int [1+\xi^{\rm rr}(\mathbf{r})]\,P(v_\parallel,\mathbf{r})\,dv_\parallel,
\]
where \(\xi^{\rm rr}\) is the real-space CCF, \(P(v_\parallel,\mathbf{r})\) is the line-of-sight velocity PDF, and the mapping between \(\mathbf{r}\) and \(\mathbf{s}\) is supplied by the galaxy peculiar velocity field. The coherent radial outflow is linked to the enclosed density contrast through the linear continuity equation,
\[
v_r(r)=-\frac{1}{3}\,f\,a\,H\,r\,\Delta(r),
\qquad
\Delta(r)=\frac{3}{r^3}\int_0^r \delta(y)\,y^2\,dy,
\]
so the anisotropy of \(\xi_{vg}\) directly probes \(f\sigma_8\) once the profile amplitude is fixed [2302.05302].

The observed anisotropy is usually compressed into Legendre multipoles,
\[
\xi_\ell(s)=\frac{2\ell+1}{2}\int_{-1}^{1} d\mu\,\xi(s,\mu)\,P_\ell(\mu),
\]
with the monopole and quadrupole carrying most of the cosmological information in BOSS analyses, and the hexadecapole often added in Euclid forecasts [1904.01030] [2302.05302]. In the BOSS configuration-space treatment where void centres are approximately in real space and galaxy velocities are radial and spherically symmetric about the void centre, the base redshift-space mapping can be written as
\[
1+\xi^s(\mathbf{s})=
\left(1+\xi^r(r)\right)
\left[1-\frac{f\Delta(r)}{3}-f\mu^2\left(\delta(r)-\Delta(r)\right)\right]^{-1},
\]
which isolates the anisotropic RSD contribution in the term proportional to \(\mu^2\) [1904.01030].

Beyond the even multipoles, the dipole \(\xi_1^{(s)}\) is a distinct higher-order observable. When the real-to-redshift mapping is extended to second order in peculiar velocity and gravitational potential, the dipole is dominated by gravitational redshift inside voids, whereas the monopole is almost unaffected and the quadrupole receives smaller higher-order corrections. In that formulation the dipole becomes a probe of the gravitational potential profile of voids rather than merely of the density or coherent outflow field [1805.05708].

## 3. Void definition and measurement practice

The measured CCF depends on how voids are defined. In VIDE/ZOBOV-like methods, the tracer field is tessellated with Voronoi cells, local densities are taken as inversely proportional to cell volumes, and a watershed transform merges low-density basins into voids; the resulting centre may be the macrocenter or volume-weighted barycentre, and the effective radius is usually
\[
R_v=\left(\frac{3V}{4\pi}\right)^{1/3}.
\]
Alternative choices include the centre of the largest empty sphere in a watershed basin, the minimum-density voxel in a gridded watershed, the area-weighted centre of a 2D photometric void, and the centre of a spherical underdensity in excursion-set-like finders. These centre definitions are not equivalent, and later systematics analyses show that the choice matters for AP and dipole measurements [2009.14751] [1904.01030] [2407.03221] [2511.04099].

The CCF is estimated from pair counts. Survey analyses often use a Landy–Szalay form,
\[
\xi(r,\mu)=
\frac{{\rm DD}_{12}-{\rm DR}_{12}-{\rm DR}_{21}+{\rm RR}_{12}}
{{\rm RR}_{12}},
\]
with voids and galaxies as the two species, while periodic-box analyses can use natural estimators such as \({\rm DD}/{\rm RR}-1\) because no survey mask is present [2302.05302] [2107.01314]. Earlier SDSS void work also used a Davis–Peebles estimator,
\[
\xi(\sigma,\pi)=\frac{DD(\sigma,\pi)}{DR(\sigma,\pi)}-1,
\]
and verified that a symmetric Landy–Szalay-type estimator gives essentially the same result for the data considered [1306.5799].

Binning conventions are analysis-specific but conceptually stable. BOSS measured \(\xi^s(s,\mu)\) in 30 bins over \(0<s<120\,h^{-1}\mathrm{Mpc}\) and 80 \(\mu\)-bins, then compressed to monopole and quadrupole [1904.01030]. The DESI forward-model study measured \(\xi_{vg}(r_p,\pi)\) in 20 logarithmic transverse bins over \(15<r_p<150\,h^{-1}\mathrm{Mpc}\) and 200 linear bins over \(0<\pi<200\,h^{-1}\mathrm{Mpc}\), then projected to \(w_{p,vg}\) [2504.08221]. Euclid forecasts use 30 radial bins up to \(120\,h^{-1}\mathrm{Mpc}\), 200 angular bins in \(\mu\), and the multipoles \(\ell=0,2,4\) as the redshift-space data vector [2302.05302].

## 4. Alcock–Paczynski distortions, reconstruction, and systematic effects

A central assumption in many AP analyses has been that void centres transform under a change of fiducial cosmology in the same way as galaxies. Stretched-box tests show that this is not generally correct: applying AP distortions to the tracer field and then running the void finder is not equivalent to finding voids first and then rescaling their centres. The response of the void finder reduces the amplitude of the AP signal in the CCF, the size of the effect depends on the void-finding package, and incorrect treatment biases recovered parameters for \(\texttt{revolver}\), \(\texttt{vide}\), \(\texttt{voxel}\), and the spherical finder in \(\texttt{Pylians3}\) [2407.02699].

A more complete redshift-space picture identifies several distinct effects on voids and on the void–galaxy CCF: a systematic expansion induced by galaxy dynamics, the AP volume effect, a systematic off-centring along the line of sight caused by void dynamics, and distortions associated with void ellipticity. Off-centring and ellipticity are detectable in projected versions of the CCF and generate additional anisotropies beyond the standard Gaussian streaming description unless they are explicitly included. In the simplified tests where all these effects are accounted for, the Gaussian streaming model remains robust [2107.01314].

Velocity-field reconstruction is therefore used to mitigate anisotropic selection bias before void finding. In BOSS, reconstruction was introduced specifically to remove the complicating effects of RSD in the void centre positions themselves, enabling a joint AP and RSD fit to the anisotropic CCF [1904.01030]. Euclid forecasts adopt the same principle: remove large-scale redshift-space distortions from the galaxy field, find voids in the reconstructed catalogue, and interpolate the void catalogue and the CCF over a grid of reconstruction efficiencies parametrized by \(\beta=f/b\). This yields nearly isotropic post-reconstruction galaxy clustering on large scales and stabilizes the modelling of both \(f\sigma_8\) and \(D_M/D_H\) [2302.05302].

## 5. Cosmological information content

The CCF already provided a low-redshift consistency test in 6dF. Fitting a self-consistent RSD model to the 2D galaxy–galaxy and void–galaxy correlation functions recovered
\[
f\sigma_8 = 0.42 \pm 0.06
\]
from galaxy clustering and
\[
f\sigma_8 = 0.39 \pm 0.11
\]
from the void–galaxy CCF, indicating consistency of the growth rate measured in overdense and underdense environments within the same dataset [1606.03092].

In BOSS CMASS, the anisotropic void–galaxy CCF became a precision AP/RSD probe. A joint fit to RSD and AP, combined with velocity-field reconstruction to remove void-centre RSD, yielded
\[
D_A(z)H(z)/c = 0.4367\pm 0.0045
\]
at \(z=0.57\), a \(1\%\) AP measurement, and
\[
f\sigma_8(z=0.57)=0.501\pm0.051,
\]
a \(10\%\) growth-rate measurement from voids alone. Combining void information with BAO and galaxy RSD in the same sample improved \(D_A/r_s\), \(H r_s\), and \(f\sigma_8\) constraints and reduced uncertainties on extended-model parameters such as \(H_0\), \(\Omega_m\), and \(w\) when combined with Planck [1904.01030].

Forecasts for Euclid spectroscopic voids push the same observable into the Stage-IV regime. With voids identified after RSD reconstruction, the expected precision is about \(0.3\%\) on the ratio \(D_M/D_H\) and between \(5\%\) and \(8\%\) on \(f\sigma_8\) in each of four redshift bins covering \(0.9\le z<1.8\). In flat \(\Lambda\)CDM this translates to
\[
\Delta\Omega_\mathrm{m}=\pm0.0028,
\]
and in flat \(w\)CDM to a dark-energy equation-of-state constraint of about \(6\%\) from voids alone [2302.05302].

Forward-model approaches reinforce the same conclusion but change the modelling strategy. An HOD-based DESI Year 5 forecast using the joint data vector \((n_v,\,w_{p,vg},\,w_{p,gg})\) predicts \(1.5\%\) and \(0.8\%\) constraints on \(\Omega_m\) and \(\sigma_8\), and shows that the combination of void and galaxy summary statistics breaks degeneracies in the galaxy–halo connection and cosmology relative to galaxy clustering alone [2504.08221]. A neural-network emulator trained on AbacusSummit and fit to the BOSS void–galaxy CCF recovered
\[
\Omega_{\rm m} = 0.330\pm 0.020,\qquad
\sigma_8 = 0.777^{+0.047}_{-0.062},
\]
together with a \(28\%\) reduction in errors for \(\Omega_{\rm m}\) relative to a template-based method, while also showing that AP errors had previously been underestimated if void centres were assumed to respond to AP in the same way as galaxies [2407.03221].

## 6. Variants, selection effects, and related extensions

The CCF is not determined by cosmology alone; it is also sensitive to how the tracer population is selected. A direct study of galaxy-property dependence found that void catalogues built from luminosity-selected galaxies and halos are consistent within errors in both the void size function and density profiles, whereas star-formation-rate-selected catalogues can differ, especially at tracer densities \(\bar n \simeq 1.9\times10^{-3}\,h^3\,\mathrm{Mpc}^{-3}\). In particular, SFR-selected voids can show a lower ridge amplitude in the density profile, consistent with the lower large-scale bias of SFR-selected galaxies relative to luminosity-selected galaxies at fixed number density [2009.14751].

A broader unifying perspective comes from void–tracer cross-correlations. In void–lensing forecasts, the same large-scale machinery is written in terms of the void–matter cross-spectrum and its projection with a tracer kernel. Replacing the lensing kernel by a galaxy kernel reproduces the tomographic void–galaxy angular spectrum, and recent CSST work models the underlying void–matter cross-power with the Halo Void Dust Model and an HSW density profile before projecting to

Source: https://www.emergentmind.com/topics/void-galaxy-cross-correlation-function