Skip to content

The Hydrodynamic Limit and Slow Variables

Hydrodynamics applies when a disturbance evolves on scales much longer than microscopic relaxation scales because conservation laws or another declared mechanism prevent its rapid decay. The hydrodynamic limit is therefore a joint limit in frequency, wave number, and any small symmetry-breaking rate—not a synonym for “locally thermal” and not a fixed list of fluid variables.

Required background. Spurions, Local Counterterms, and Symmetry Response supplies the source-response Ward identities from which conserved-density poles follow.

Helpful background. Linearized Kinetics and Relaxation Modes shows the microscopic separation between collision zero modes and gapped modes. Memory Functions and Slow-Mode Projection develops the complementary projection formalism.

Consider a conserved scalar density nn and current j\mathbf j,

tn+ ⁣ ⁣j=0.\partial_t n+\boldsymbol\nabla\!\cdot\!\mathbf j=0.

If all nonconserved local variables relax in a time τmicro\tau_{\mathrm{micro}}, locality and rotational invariance permit the constitutive expansion

j=Dn+O(2,n2).\mathbf j=-D\boldsymbol\nabla n+O(\partial^2,n^2).

For eiωt+ikxe^{-i\omega t+i\mathbf k\cdot\mathbf x}, conservation closes to

iω+Dk2+O(k4)=0,ω(k)=iDk2+O(k4).-i\omega+Dk^2+O(k^4)=0, \qquad \omega(k)=-iDk^2+O(k^4).

The rate vanishes as k0k\to0 because a spatially uniform conserved excess cannot decay. A retarded density correlator consequently has the low-energy form

GnnR(ω,k)=χDk2iω+Dk2+terms analytic near (0,0),G^R_{nn}(\omega,k) = -\chi\,\frac{Dk^2}{-i\omega+Dk^2} +\text{terms analytic near }(0,0),

where χ\chi is the static susceptibility. This uses the site convention GR=iθ(t)[n(t),n(0)]G^R=-i\theta(t)\langle[n(t),n(0)]\rangle: for a Hamiltonian source HHμnH\to H-\int\mu n, the physical response is GR-G^R and GR(0,k0)=χG^R(0,k\to0)=-\chi. Kadanoff and Martin derive this connection between local conservation, constitutive closure, and hydrodynamic singularities Kadanoff and Martin 1963, §§2–3, pp. 421–430.

Energy and longitudinal momentum instead form a reversible pair. Conservation gives a matrix linear in kk, producing ω=±csk+O(k2)\omega=\pm c_s k+O(k^2) rather than diffusion. The mechanism is the same: exact conservation forces eigenvalues to approach zero, while reversible couplings determine whether the approach is linear or quadratic in kk.

Let τmicro\tau_{\mathrm{micro}} and micro\ell_{\mathrm{micro}} characterize the least-damped omitted mode and the shortest length resolved by a local constitutive law. Ordinary hydrodynamics requires

ωτmicro1,kmicro1.|\omega|\tau_{\mathrm{micro}}\ll1, \qquad k\ell_{\mathrm{micro}}\ll1.

These inequalities concern derivatives, not amplitudes. Nonlinear hydrodynamics can describe order-one changes spread over long distances; linear response imposes a separate small-amplitude approximation. Fluctuating hydrodynamics adds a loop or noise expansion, and near criticality that expansion can fail even while kmicrok\ell_{\mathrm{micro}} is small.

The constitutive relation integrates out modes with nonzero relaxation rates at k=0k=0. If the smallest omitted rate is not well separated from the target frequency, the resulting kernels become strongly frequency dependent and a local derivative expansion ceases to be economical. Kovtun gives this pole-based formulation and its charged and stress-tensor examples Kovtun 2012, §2, pp. 19–29, Open PDF.

Suppose a density obeys

tn+ ⁣ ⁣j=Γn,Γτmicro1.\partial_t n+\boldsymbol\nabla\!\cdot\!\mathbf j=-\Gamma n, \qquad \Gamma\tau_{\mathrm{micro}}\ll1.

Then

ω=iΓiDk2+.\omega=-i\Gamma-iDk^2+\cdots .

This mode is not hydrodynamic in the strict k0k\to0 sense because its gap is nonzero. It is nevertheless quasihydrodynamic on times tΓ1t\ll\Gamma^{-1} or whenever Γ\Gamma is part of the controlled small-parameter counting. Integrating it out at ωΓ\omega\sim\Gamma produces a memory kernel; retaining it restores a local enlarged system.

The same decision arises for a weakly broken symmetry, a critical order parameter with critical slowing down, spin density with a long relaxation time, axial charge with a small violation rate, or a nonhydrodynamic stress mode promoted in a transient theory. The page Hydro+ and Parametrically Slow Critical Modes treats the critical example; the present criterion is the universal selection rule.

A charged normal fluid away from criticality usually retains energy density, momentum density, and each independent conserved charge density. This set is incomplete when:

  • a continuous symmetry is spontaneously broken, requiring a Goldstone phase;
  • a critical mode relaxes on the hydrodynamic time scale;
  • integrability supplies infinitely many relevant conservation laws;
  • a background field changes which momentum components are conserved;
  • an external drive removes energy conservation or creates a new slow phase;
  • an initial state supports a parametrically narrow quasiparticle or coherent mode.

No pole calculation can repair an omitted variable. A fit may hide the omission in a frequency-dependent “transport coefficient,” but that coefficient is then not the local constitutive datum it was assumed to be.

Three checks catch most mistakes. First, at k=0k=0 every exactly conserved density must have zero relaxation rate. Second, residues must obey the relevant Ward identity and static susceptibility. Third, adding a small breaking rate must move the pole continuously from iDk2-iDk^2 to iΓiDk2-i\Gamma-iDk^2, not remove it discontinuously.

The hydrodynamic pole is an infrared statement. Extending iDk2-iDk^2 to arbitrarily large kk gives a parabolic PDE with instantaneous tails, but microscopic causality is not decided by extrapolating a truncated low-kk series. Whether a chosen macroscopic formulation itself has a causal well-posed initial-value problem is the separate subject of Strong Hyperbolicity, Stability, and Causal Propagation.

The first two boxes in the schematic isolate this page’s central step: conservation selects slow densities, while coarse graining is the additional physical assumption that permits a local-equilibrium description. Follow the remaining boxes only as a preview of the closures that later pages will test.

A left-to-right chain runs from conserved currents and slow densities through local equilibrium, hydrostatic constraints, and constitutive tensors to ideal linear modes; a dashed branch from the constitutive box says that a frame change redefines fields rather than physical transport data.

Conservation provides the parametrically slow densities but does not by itself prove local equilibration. The later source, constitutive, and mode stages require the scale separation and completeness hypotheses stated on this page. The diagram is a qualitative dependency map, not a plot of relaxation times or a claim that every conserved density is sufficient for closure.

In text: identify all exact and approximate conservation laws; verify a gap to omitted modes; assume local equilibration only inside the resulting frequency and wavelength window; then construct and test a constitutive closure. If another mode becomes parametrically slow, it must be added before the later stages in the diagram are valid.

A nearly conserved charge has constitutive current j=Dn\mathbf j=-D\nabla n and relaxation rate Γ\Gamma. Couple a source μ\mu so that the equilibrium response is δn=χμ\delta n=\chi\mu. Find a retarded correlator with the correct static limit.

Solution

The sourced relaxation equation can be written

(t+ΓD2)δn=(ΓD2)χμ.(\partial_t+\Gamma-D\nabla^2)\delta n = (\Gamma-D\nabla^2)\chi\mu.

Therefore

GnnR(ω,k)=χΓ+Dk2iω+Γ+Dk2.G^R_{nn}(\omega,k) = -\chi\, \frac{\Gamma+Dk^2}{-i\omega+\Gamma+Dk^2}.

At ω=0\omega=0 it equals χ-\chi, so the response GR-G^R equals χ\chi; at vanishing source its pole is ω=i(Γ+Dk2)\omega=-i(\Gamma+Dk^2). Contact-term conventions can redistribute analytic pieces, but not this pole and its static normalization.

Conservation Laws and Hydrodynamic Fields now turns the slow densities into covariant stress and current variables. Generalized Hydrodynamics of Integrable Systems treats the distinct case in which the conserved set is extensive or infinite.

  • Kadanoff, Leo P., and Paul C. Martin. 1963. “Hydrodynamic Equations and Correlation Functions.” Annals of Physics 24: 419–469. DOI.

  • Kovtun, Pavel. 2012. “Lectures on Hydrodynamic Fluctuations in Relativistic Theories.” Journal of Physics A: Mathematical and Theoretical 45: 473001. DOI. Open PDF.