Pair Correlation Functions Overview

• Pair correlation functions are two-point statistics that quantify the probability of finding a second particle at a given separation, normalized to one in uncorrelated regimes.
• They are expressed in various forms—such as radial distribution functions in liquids, correlation functions in heavy-ion femtoscopy, and spectral statistics in number theory—across diverse scientific fields.
• PCFs are pivotal for interpreting experimental data and validating models through methods like integral equations, machine learning, and obstacle corrections in empirical setups.

Pair correlation functions (PCFs) are two-point statistics that quantify how pairs are distributed relative to an uncorrelated reference. In liquids, plasmas, glasses, and electron systems they are commonly written as radial distribution functions $g(r)$ or $g_{\alpha\beta}(r)$; in spatial point-process theory they appear as $g_0(h)$; in heavy-ion femtoscopy they are written as two-particle correlation functions $C(q)$; in one-dimensional model sets they are averaged pair correlations $\nu(y)$; and in the analytic theory of the Riemann zeta function they appear as pair-correlation sums such as $F(\alpha,T)$ (Ayush et al., 2023, Shaw et al., 2020, Xiao et al., 27 Oct 2025, Baake et al., 27 Feb 2025, Banks, 27 Feb 2025). Across these settings, the normalization is typically chosen so that the statistic tends to $1$ in an uncorrelated or asymptotic regime, while deviations from $1$ encode ordering, exclusion, clustering, screening, or oscillatory structure.

1. Definitions, normalizations, and basic objects

In atomistic statistical mechanics, the PCF measures the probability of finding a second particle at a separation $r$ relative to an ideal gas at the same density. For glasses, the pair correlation function or radial distribution function $g_{\alpha\beta}(r)$ is defined for atom types $g_{\alpha\beta}(r)$0 and $g_{\alpha\beta}(r)$1 through interatomic separations $g_{\alpha\beta}(r)$2 and an ensemble average over molecular-dynamics trajectories (Ayush et al., 2023). For the unpolarized uniform electron gas, the spin-averaged pair correlation $g_{\alpha\beta}(r)$3 is defined so that $g_{\alpha\beta}(r)$4 is the probability to find a second electron between $g_{\alpha\beta}(r)$5 and $g_{\alpha\beta}(r)$6, with spin-resolved functions $g_{\alpha\beta}(r)$7, $g_{\alpha\beta}(r)$8, and $g_{\alpha\beta}(r)$9 entering the unpolarized average (Larkin et al., 2022). In crystals and in classical density-functional theory, one starts from the one-particle density $g_0(h)$0, the two-particle density $g_0(h)$1, the pair distribution function $g_0(h)$2, and the total pair correlation $g_0(h)$3 (Jaiswal et al., 2014, Singh, 18 May 2026).

For spatial point processes under second-order intensity-reweighted stationarity, the central relation is

$g_0(h)$4

where translation invariance is imposed only after reweighting by the intensity field $g_0(h)$5 (Shaw et al., 2020). The associated inhomogeneous $g_0(h)$6-function is

$g_0(h)$7

and in isotropic settings one writes $g_0(h)$8 (Shaw et al., 2020).

In heavy-ion femtoscopy, the object of interest is the correlation function of emitted fragments as a function of relative momentum,

$g_0(h)$9

with

$C(q)$0

where $C(q)$1 is the same-event histogram, $C(q)$2 is the mixed-event reference, and $C(q)$3 fixes $C(q)$4 (Xiao et al., 27 Oct 2025). This formulation makes explicit that femtoscopic PCFs are experimentally constructed from coincidence yields rather than directly from spatial coordinates.

In one-dimensional regular model sets, the averaged pair correlation is

$C(q)$5

and it is equivalently the mass of the autocorrelation measure at $C(q)$6 (Baake et al., 27 Feb 2025). For $C(q)$7-type decompositions $C(q)$8, one refines this to cross-correlations $C(q)$9 (Baake et al., 27 Feb 2025).

In Montgomery-type pair correlation for zeta zeros, assuming RH so that nontrivial zeros are $\nu(y)$0, one defines

$\nu(y)$1

with $\nu(y)$2 in the classical case (Banks, 27 Feb 2025). The terminology is the same, but the “pairs” are pairs of zero ordinates rather than particles or points in Euclidean space.

2. Physical interpretation and representative regimes

The most standard interpretation of a PCF is that $\nu(y)$3 signals the ideal-gas limit or decay of correlations, while peaks signal preferred length scales. This is explicit in the classical radial PCF, where $\nu(y)$4 gives the expected number of particles in a shell at distance $\nu(y)$5, and where peaks in $\nu(y)$6 encode preferred separations (Worlitzer et al., 2022). In cosmology, the galaxy pair-correlation function $\nu(y)$7 is the excess probability, relative to a Poisson distribution, of finding two galaxies separated by $\nu(y)$8; oscillations in the range $\nu(y)$9–$F(\alpha,T)$0 are identified with baryon acoustic oscillations (Zaninetti, 2014).

In warm dense matter, the onset of short-range order is visible as a first peak in $F(\alpha,T)$1. In the unpolarized uniform electron gas, SMPIMC calculations show that when the classical coupling $F(\alpha,T)$2 exceeds roughly $F(\alpha,T)$3–$F(\alpha,T)$4, depending on degeneracy, a pronounced first peak develops, and this appearance of short-range order is accompanied by a high-momentum “quantum tail” in the momentum distribution (Larkin et al., 2022). The paper interprets these as two aspects of the same many-body effect: transient effective wells formed by neighboring electrons and tunneling through those wells (Larkin et al., 2022).

In two-temperature strongly coupled plasmas, PCFs resolve interspecies screening and the validity of closure theories. Molecular-dynamics comparisons show that the Seuferling–Vogel–Teopffer model reproduces electron-electron, ion-ion, and electron-ion radial distribution functions well over a wide range of mass ratios, temperature ratios, and couplings, while the Yukawa OCP description of ion-ion correlations is accurate only up to $F(\alpha,T)$5, beyond which it rapidly breaks down (Shaffer et al., 2017).

In fractional quantum Hall states, the pair-correlation function $F(\alpha,T)$6 is a reduced radial object constructed from the trial wavefunction after fixing one particle at the origin and another at separation $F(\alpha,T)$7. Its thermodynamic-limit form feeds directly into the projected static structure factor and the single-mode-approximation estimate of the magneto-roton gap (Fulsebakke et al., 2022). In the two-dimensional Hubbard model, the phrase “pair correlation” refers instead to superconducting pairing correlators such as $F(\alpha,T)$8, where enhancement away from half-filling and growth with system size are used as a criterion for superconductivity (Yanagisawa, 2013). The shared terminology therefore covers both geometric pair distributions and operator-valued pairing correlators.

3. Integral equations, closures, and broken symmetry

A large part of PCF theory is organized by the Ornstein–Zernike equation. In a general inhomogeneous system,

$F(\alpha,T)$9

with $1$0 the direct pair correlation function (Jaiswal et al., 2014, Singh, 18 May 2026). In homogeneous fluids this reduces to the usual convolution form, and a closure relation such as Percus–Yevick, hypernetted-chain, Rogers–Young, or Zerah–Hansen is then required (Singh, 18 May 2026).

For crystals and other ordered phases, the essential step is to decompose one- and two-particle quantities into symmetry-conserving and symmetry-broken parts. The density is written as $1$1, and the pair correlations as

$1$2

where $1$3 and $1$4 are homogeneous-fluid correlations at the average density, and $1$5, $1$6 are periodic in the center-of-mass coordinate and encode the broken translational symmetry (Jaiswal et al., 2014). The broken direct correlation function is then expanded in powers of the order parameters by using higher-body direct correlation functions of the uniform fluid (Jaiswal et al., 2014, Singh, 18 May 2026). In the two-dimensional hexagonal soft-disk example, truncation at the three-body term already captures more than $1$7 of the full $1$8 near freezing (Jaiswal et al., 2014).

Within classical density-functional theory, these PCFs are not merely diagnostic. The excess free-energy functional generates the direct correlations through functional differentiation, and the grand-potential difference between competing phases contains both a symmetry-conserving contribution and a symmetry-breaking contribution involving $1$9 (Singh, 18 May 2026). This is the basis for the “exact” DFT described in the review, where directly computed PCFs of broken-symmetry phases are used to predict freezing, nematic–isotropic transitions, and crystal–crystal transitions (Singh, 18 May 2026).

Heavy-ion femtoscopy uses a different dynamical framework but addresses the same two-point object. The classical trajectory approximation CTA-I combines a thermal equilibrium emission source, sequential emission with an exponential time delay, a self-consistent mean field that is harmonic at small $1$0 and Coulombic at large $1$1, and the Coulomb interaction between the emitted fragments (Xiao et al., 27 Oct 2025). The resulting three-body Monte Carlo propagates the source nucleus and two emitted particles to asymptotic momenta before constructing $1$2 from accepted events and mixed events (Xiao et al., 27 Oct 2025).

4. Estimation, normalization, and corrections in empirical settings

The numerical value of a PCF depends critically on the reference measure used in its denominator. In inhomogeneous spatial point processes, the “global” estimator replaces pointwise intensity products by the overlap integral

$1$3

leading to the unbiased estimators

$1$4

Under SOIRS, these remain unbiased, and the simulations in the paper show reduced bias and variance relative to local normalization, especially when the intensity must itself be estimated (Shaw et al., 2020).

On discrete lattices, exact normalization requires counting all site pairs at a given lattice distance. For occupied sites $1$5, metric $1$6, and separation $1$7, the discrete PCF is

$1$8

where $1$9 counts occupied-occupied pairs and $r$0 counts site-site pairs at distance $r$1 (Gavagnin et al., 2018). On square lattices, taxicab and uniform metrics yield isotropic, correctly normalized PCFs, unlike the annular PCF, which is poorly calibrated on a lattice, or the rectilinear PCF, which is anisotropic and misses diagonal patterns such as chessboards and diagonal stripes (Gavagnin et al., 2018).

When obstacles are present, naive normalization is no longer valid because straight-line lattice distances do not respect connectivity. The corrected pair-correlation function replaces the standard pair count $r$2 by

$r$3

thereby removing accessible-obstacle and obstacle-obstacle artifacts and reassigning pairs whose shortest paths are altered by obstacle avoidance (Johnston et al., 2018). The paper shows that standard PCFs can produce spurious oscillations and false long-range aggregation in obstacle-strewn domains, while the corrected PCF recovers the genuine short-range peak at $r$4 associated with birth clustering (Johnston et al., 2018).

A related critique applies to the classical radial PCF more generally. Because it averages all pairs in metric shells, it can wash out subtle local rearrangements, rare defects, or local symmetries. The Voronoi PCF addresses this by replacing Euclidean shells with topological shells in the Delaunay graph: $r$5 is the minimum number of Delaunay edges in a path from $r$6 to $r$7, $r$8 is the mean number of $r$9-neighbors, and

$g_{\alpha\beta}(r)$0

This produces a discrete, topology-based PCF that is scale-invariant, robust under infinitesimal perturbations, and sensitive to local connectivity (Worlitzer et al., 2022). The same paper also notes a limitation: averaging over all particles still loses higher-order distributional information, so one may need the variance or full distribution of shell populations for finer discrimination (Worlitzer et al., 2022).

5. Aperiodic order and analytic number theory

In regular model sets, PCFs are linked directly to covariograms of internal-space windows. For a cut-and-project set with windows $g_{\alpha\beta}(r)$1, one has

$g_{\alpha\beta}(r)$2

with $g_{\alpha\beta}(r)$3 outside the corresponding difference set (Baake et al., 27 Feb 2025). This identification converts pair-correlation questions into convolution questions in internal space.

When the model set comes from a primitive inflation, the cross-correlations satisfy exact renormalisation equations. In the paper’s formulation,

$g_{\alpha\beta}(r)$4

and the window covariograms satisfy the corresponding internal-space contraction equation (Baake et al., 27 Feb 2025). In the sister silver mean example, recursion on a self-consistent finite set determines all pair-correlation values, yielding a continuous but extremely spiky graph. In the “splitting” Cantorval example, the numerical sample reveals two disjoint branches associated with an even or odd number of $g_{\alpha\beta}(r)$5-tiles in the digital decomposition, and convergence to the continuous curve is very slow (Baake et al., 27 Feb 2025). The paper characterizes the resulting covariograms as unexpectedly complex and wild.

In analytic number theory, the PCF language refers to distributions of differences between zero ordinates. Banks introduces a two-parameter family of weights

$g_{\alpha\beta}(r)$6

leading to weighted pair-correlation functions $g_{\alpha\beta}(r)$7 (Banks, 27 Feb 2025). Although these weights alter the asymptotic behavior of the weighted PCF, the paper states that they do not yield any new information about the simplicity of the zeros (Banks, 27 Feb 2025). The same work extends Montgomery’s approach from differences of ordinates to sums $g_{\alpha\beta}(r)$8, defining $g_{\alpha\beta}(r)$9 and formulating a generalized strong pair-correlation conjecture for any integer $g_{\alpha\beta}(r)$00 (Banks, 27 Feb 2025). This suggests that “pair correlation” in the number-theoretic literature is best understood as a spectral two-point statistic rather than a geometric radial function.

6. Computational representations, surrogate models, and domain-specific sensitivities

PCFs are increasingly treated as high-dimensional objects that can be learned, compressed, or stably parametrized. In glass modeling, the machine-learning pipeline of a CNN autoencoder plus random-forest regression compresses each $g_{\alpha\beta}(r)$01 grayscale PCF image into an $g_{\alpha\beta}(r)$02-dimensional latent vector and predicts that vector from a $g_{\alpha\beta}(r)$03-dimensional composition-and-pair fingerprint (Ayush et al., 2023). Direct regression on the full $g_{\alpha\beta}(r)$04-point $g_{\alpha\beta}(r)$05 curves yields $g_{\alpha\beta}(r)$06, whereas encoding plus regression plus decoding attains $g_{\alpha\beta}(r)$07 in latent-space prediction and reconstructed-curve coefficients of determination exceeding $g_{\alpha\beta}(r)$08 on all nine unseen glass compositions (Ayush et al., 2023). A plausible implication is that low-dimensional latent structure captures the dominant peak-position and peak-height variability more efficiently than naive coordinatewise regression.

For fractional quantum Hall wavefunctions, stable representation rather than generic regression is the central issue. Girvin’s original basis is physically appealing but numerically unstable because its overlap matrix is ill-conditioned; the orthogonalized spherical basis built from modified Jacobi polynomials yields coefficients that can be obtained stably by projection, and these coefficients extrapolate linearly in $g_{\alpha\beta}(r)$09 to the thermodynamic limit (Fulsebakke et al., 2022). The paper further reports that, for all states considered, the large-$g_{\alpha\beta}(r)$10 expansion coefficients are fit remarkably well by a cosine oscillation with exponentially decaying amplitude (Fulsebakke et al., 2022).

In femtoscopic heavy-ion applications, the most important practical sensitivity in CTA-I is not the thermal spectrum but the source geometry. The model shows that $g_{\alpha\beta}(r)$11 is highly sensitive to the Gaussian source size $g_{\alpha\beta}(r)$12 and only weakly dependent on the temperature parameter $g_{\alpha\beta}(r)$13. For Ar+Au$g_{\alpha\beta}(r)$14B–B at $g_{\alpha\beta}(r)$15, changing $g_{\alpha\beta}(r)$16 from $g_{\alpha\beta}(r)$17 to $g_{\alpha\beta}(r)$18 visibly modifies the Coulomb anti-correlation dip around $g_{\alpha\beta}(r)$19–$g_{\alpha\beta}(r)$20, whereas varying $g_{\alpha\beta}(r)$21 from $g_{\alpha\beta}(r)$22 to $g_{\alpha\beta}(r)$23 leaves the correlation nearly unchanged; similar behavior appears in $g_{\alpha\beta}(r)$24–$g_{\alpha\beta}(r)$25 and $g_{\alpha\beta}(r)$26–$g_{\alpha\beta}(r)$27 at $g_{\alpha\beta}(r)$28 (Xiao et al., 27 Oct 2025). The proposed workflow is therefore to determine $g_{\alpha\beta}(r)$29 from single-particle spectra and use the PCF to scan $g_{\alpha\beta}(r)$30 and extract the spatio-temporal size of the emission region (Xiao et al., 27 Oct 2025).

A recurrent methodological point is that a PCF is powerful precisely because it is low-order and broadly measurable, but that same averaging can be a limitation. This is explicit in the contrast between shell-averaged radial PCFs and Voronoi-topological PCFs (Worlitzer et al., 2022), in the need for obstacle corrections on lattices (Johnston et al., 2018), and in the need for orthogonal bases or renormalisation equations when direct numerical representations become unstable or fractal-like (Fulsebakke et al., 2022, Baake et al., 27 Feb 2025). The cumulative picture is that PCFs remain foundational descriptors of two-point structure, yet their most informative use depends on choosing the normalization, geometry, closure, and representation that match the underlying domain.