Skip to content

Correlations, Susceptibilities, and Correlation Lengths

Static response and spatial correlation are two views of the same equilibrium fluctuations. A susceptibility is a source derivative and, after contact terms and normalizations are matched, the integral of a connected correlation function. A correlation length can be read from asymptotic decay, a screening pole, or a second moment, but these estimators agree only when one isolated long-distance scale dominates. None is automatically a real-time quasiparticle mass.

The scaling relations among correlations, susceptibilities, and correlation lengths are reviewed with their critical assumptions in Fisher 1967, pp. 615–730.

Required background. Partition Functions and Thermodynamic Response supplies the source derivatives. Connected Correlators and Cumulants fixes connected normalization. Helpful background. Critical Exponents, Scaling Relations, and Hyperscaling Caveats develops the scaling limit used near criticality.

For a translation-invariant scalar observable O(x)\mathcal O(\mathbf x), define

C(x)=O(x)O(0)O2.C(\mathbf x)= \langle\mathcal O(\mathbf x)\mathcal O(\mathbf0)\rangle -\langle\mathcal O\rangle^2.

Couple a uniform source through Hh=HhVdd1xOH_h=H-h\int_V\mathrm d^{d-1}x\,\mathcal O. Then

χ=Oh=βVdd1xC(x),\chi =\frac{\partial\langle\mathcal O\rangle}{\partial h} =\beta\int_V\mathrm d^{d-1}x\,C(\mathbf x),

where O\langle\mathcal O\rangle denotes the spatially averaged density. If the source also changes the operator or action explicitly, a contact term must be added. On a finite periodic lattice the same identity holds as a sum, with the zero mode and volume normalization stated explicitly.

The structure factor is

S(k)=dd1xeikxC(x).S(\mathbf k)=\int\mathrm d^{d-1}x\, e^{-i\mathbf k\cdot\mathbf x}C(\mathbf x).

Although the global Fourier convention uses e+ipxe^{+ip\cdot x}, this spatial definition is written explicitly and is even in k\mathbf k for a parity-invariant real channel. With the convention above, χ=βS(0)\chi=\beta S(\mathbf0) after any required contact subtraction.

Exponential, screening, and second-moment lengths

Section titled “Exponential, screening, and second-moment lengths”

If one isolated screening state controls large separation,

C(r)Ar(d2)/2er/ξexp.C(r)\sim A\,r^{-(d-2)/2}e^{-r/\xi_{\mathrm{exp}}}.

Equivalently S(k)S(\mathbf k) has a nearest singularity at k2=ξexp2\mathbf k^2=-\xi_{\mathrm{exp}}^{-2}. At small real momentum an Ornstein–Zernike form reads

S(k)S(0)1+ξ22k2+O(k4),S(\mathbf k)\simeq \frac{S(0)}{1+\xi_2^2\mathbf k^2+O(k^4)},

with second-moment length

ξ22=12(d1)dd1xx2C(x)dd1xC(x).\xi_2^2 =\frac{1}{2(d-1)} \frac{\int\mathrm d^{d-1}x\,\mathbf x^2 C(\mathbf x)} {\int\mathrm d^{d-1}x\,C(\mathbf x)}.

On a periodic lattice of extent LL a common estimator avoids coordinate wrapping:

ξ2,L=12sin(π/L)[S(0)S(kmin)1]1/2,kmin=2πL(1,0,).\xi_{2,L} =\frac{1}{2\sin(\pi/L)} \left[ \frac{S(\mathbf0)}{S(\mathbf k_{\min})}-1 \right]^{1/2}, \qquad \mathbf k_{\min}=\frac{2\pi}{L}(1,0,\ldots).

This formula assumes an approximately isotropic single-pole small-momentum form. Direction-dependent estimators are required on anisotropic lattices; a multi-exponential channel need not yield ξ2=ξexp\xi_2=\xi_{\mathrm{exp}}.

For the static Gaussian functional

F[ϕ]=12dd1x[(ϕ)2+ms2ϕ2],\mathcal F[\phi] =\frac12\int\mathrm d^{d-1}x \left[(\nabla\phi)^2+m_s^2\phi^2\right],

the covariance is

S(k)=Tk2+ms2.S(\mathbf k)=\frac{T}{\mathbf k^2+m_s^2}.

Therefore

ξ=1ms,χ=βS(0)=1ms2\xi=\frac1{m_s}, \qquad \chi=\beta S(0)=\frac1{m_s^2}

for the displayed source normalization. This is a screening mass: it controls static spatial decay. A pole mass comes from a real-time frequency singularity, and a curvature mass from a local effective-potential Hessian. Interactions can separate all three.

Near a continuous transition, a scaling form is

C(r;t)1rd3+ηC(r/ξ),ξtν,χξ2η.C(r;t)\sim \frac{1}{r^{d-3+\eta}}\, \mathcal C(r/\xi), \qquad \xi\sim\lvert t\rvert^{-\nu}, \qquad \chi\sim\xi^{2-\eta}.

Finite volume cuts off the divergence when ξ\xi approaches LL. Extracting a critical exponent therefore requires a joint finite-size and scaling-window analysis; fitting a large but finite ξ\xi at one volume does not establish a critical point.

  • Confirm that connected rather than full correlations are integrated.
  • State the source sign and any contact or zero-mode subtraction.
  • Compare exponential and second-moment estimators over increasing volumes.
  • Vary the fit window and include more than one screening state when required.
  • Test anisotropy by rotating kmin\mathbf k_{\min}.
  • Keep screening, curvature, pole, and damping scales distinct.
  • Near criticality, demonstrate L/ξL/\xi control rather than extrapolating a single volume.

The finite-volume and limit distinctions are summarized in the shared equilibrium table. Gaussian determinants and the failure of quadratic control continue on Gaussian Fluctuations and Effective Free Energy.

The schematic below organizes the relationships used on this page. Inspect it with this question in mind: How do Euclidean configurations, correlations, Gaussian fluctuations, phases, and thermodynamic limits relate?

A Euclidean weight defines correlations; a positive Hessian licenses Gaussian fluctuations locally, while phase coexistence and spontaneous breaking require a separately ordered thermodynamic and source limit.

A Euclidean weight defines correlations; a positive Hessian licenses Gaussian fluctuations locally, while phase coexistence and spontaneous breaking require a separately ordered thermodynamic and source limit. Solid connections show the primary relation; dashed outlines or arrows mark qualifications and failure boundaries. The diagram is schematic and not to scale.

The surrounding discussion supplies the relevant equations and checks in text form; the figure is a navigational summary.

For S(k)=A/(k2+m12)+B/(k2+m22)S(k)=A/(k^2+m_1^2)+B/(k^2+m_2^2) with 0<m1<m20<m_1<m_2 and A,B>0A,B>0, compare the exponential and second-moment lengths.

Solution

The large-distance decay is controlled by the nearest pole, so ξexp=1/m1\xi_{\mathrm{exp}}=1/m_1. Expanding at small kk gives ξ22=(A/m14+B/m24)/(A/m12+B/m22)\xi_2^2=(A/m_1^4+B/m_2^4)/(A/m_1^2+B/m_2^2). Unless the light pole dominates the zero-momentum weight, ξ2<1/m1\xi_2<1/m_1. The difference is a diagnostic of multiple scales, not a contradiction.