Skip to content

Thermal, Hydrodynamic, and Nonequilibrium Holography

Scrollable formula. Use Left and Right arrow keys to reveal hidden content.

Scrollable table. Use Left and Right arrow keys to reveal additional columns.

Scrollable figure. Use Left and Right arrow keys to inspect hidden content.

Thermal holography turns a black brane into a calculational engine for retarded response, transport, hydrodynamic constitutive relations, and nonequilibrium evolution. This chapter follows one causal chain: impose ingoing horizon data, extract a retarded correlator, take its low-frequency limit, organize long-wavelength dynamics, and then test where that organization fails during quenches, collisions, or turbulence. Frontier evidence is assessed through 2026-08-20.

Helpful background. Gauge Theory and String-Theory Preparation supplies the holographic operator map; Quantum Gravity supplies semiclassical control; Quantum Information supplies entropy diagnostics; and Black-Hole Perturbations supplies quasinormal boundary conditions and mode relaxation.

For a bulk perturbation ϕ\phi dual to an operator O{\cal O}, solve the linearized equation with ingoing horizon behavior and a fixed boundary source JJ. Let ΠJ\Pi_J be the renormalized canonical response conjugate to that source, including the field’s boundary scaling. At linear order,

GR(ω,k)=−ΠJ(ω,k)J(ω,k).G_R(\omega,\mathbf{k})=-\frac{\Pi_J(\omega,\mathbf{k})} {J(\omega,\mathbf{k})}.

Its poles are quasinormal modes; its small-frequency slope gives transport; and promoting equilibrium parameters to slowly varying fields produces fluid dynamics. Far from equilibrium, the same dictionary extracts a conserved boundary stress tensor from a nonlinear bulk evolution, but no single observable defines “thermalized.”

Read the pages in this order:

  1. Real-Time Holographic Correlators derives the causal prescription, contact-term qualifications, and pole interpretation.
  2. Kubo Formulae and Holographic Transport takes controlled low-frequency limits and checks radial flux.
  3. Fluid/Gravity Correspondence derives boundary constitutive equations from regular long-wavelength black-brane geometries.
  4. Viscosity, Universality, and Corrections separates the two-derivative result from finite-coupling and higher-derivative violations.
  5. Holographic Quenches and Thermalization compares probe clocks and provides a deterministic AdS–Vaidya benchmark.
  6. Far-From-Equilibrium Dynamics and Turbulence assesses shocks, hydrodynamization, attractors, and cascades across model and data choices.

From horizon regularity to constitutive dynamics

Section titled “From horizon regularity to constitutive dynamics”

With the global e−iωte^{-i\omega t} convention, an ingoing mode is analytic in a future-horizon coordinate. This selects the retarded, rather than advanced, solution. Son and Starinets connect the resulting boundary flux to Minkowski correlators in Son and Starinets 2002, §§ 3–4. For shear response,

η=−lim⁡ω→01ωIm⁡GRxy,xy(ω,0).\eta=-\lim_{\omega\to0} \frac{1}{\omega}\operatorname{Im}G_R^{xy,xy}(\omega,\mathbf0).

At two-derivative Einstein order, regularity relates the conserved low-frequency flux to horizon area, producing η/s=1/(4π)\eta/s=1/(4\pi) in a broad class of isotropic translationally invariant plasmas. That result is a theorem of a restricted effective action, not a universal inequality for all quantum field theories.

Fluid/gravity promotes temperature and velocity to fields and solves Einstein’s equations order by order in boundary derivatives. Constraint equations become ∇μTμν=0\nabla_\mu T^{\mu\nu}=0; dynamical radial equations determine transport coefficients. The expansion is useful when constitutive residuals are small, even if pressure anisotropy is not. It is not identical to the late-time convergence of an infinite gradient series.

The matrix records what is controlled and which axis breaks the simplest interpretation.

assumptions and settingcontrolled evidence or claimadverse limitation or discriminatorprimary source
Linear perturbation of a stationary AdS black brane, fixed counterterms, ingoing future-horizon conditionBoundary canonical momentum gives the retarded correlator; quasinormal poles encode relaxationGauge-invariant variables, contact terms, alternate quantization, and coupled channels must be handled before interpreting a raw ratioSon and Starinets 2002, §§ 3–4
Isotropic translationally invariant plasma with a two-derivative Einstein dualThe shear channel gives η/s=1/(4π)\eta/s=1/(4\pi) and a regular membrane fluxHigher-derivative interactions can change the graviton coupling and violate the proposed universal lower boundPolicastro, Son, and Starinets 2001, complete Letter, pp. 1–4
Slowly varying boosted AdS black brane with regular horizon and derivative expansionBulk constraint equations reproduce conformal hydrodynamics and determine nonlinear transportLarge gradients, extra slow modes, phase transitions, or singular horizon data lie outside the controlled expansionBhattacharyya et al. 2008, §§ 3–5
Homogeneous thin-shell AdS–Vaidya probes in the large-dimension saddle approximationShorter-scale correlators and surfaces approach thermal values before longer-scale entanglement probes; interval entropy saturates a causal time in AdS3Probe ordering depends on geometry, observable, subtraction, and tolerance; event-horizon formation is not a boundary thermalization timeBalasubramanian et al. 2011, §§ III.1–III.3, pp. 23–36
Long-wavelength AdS4 black-hole flows with effectively two-dimensional boundary turbulenceGauge-invariant boundary spectra support an inverse cascade in the hydrodynamic regimeThe cascade direction and inertial range depend on dimension and scale separation; this is not evidence for one universal turbulent black-hole cascadeAdams, Chesler, and Liu 2014, complete Letter, pp. 1–5

Competing positions and discriminating tests

Section titled “Competing positions and discriminating tests”

One position emphasizes universality: horizon regularity fixes low-frequency dissipation, and many strongly coupled plasmas hydrodynamize rapidly. The evidence for this is robust within two-derivative large-N models. The competing position emphasizes effective-theory dependence: higher-derivative terms, charge, broken translation symmetry, anisotropy, or new slow modes alter transport and sometimes the very set of hydrodynamic variables. The decisive test is not whether a number is close to 1/(4π)1/(4\pi), but whether all omitted couplings and symmetry assumptions are parametrically controlled.

A similar distinction governs nonequilibrium claims. A small constitutive residual supports hydrodynamization for a chosen stress tensor and truncation; it does not prove isotropization, entanglement saturation, or a thermal density matrix. An event horizon is teleological, an apparent horizon is slicing dependent, and a nonlocal probe samples a scale-dependent bulk region. These clocks answer different questions.

Turbulent and attractor claims require the strongest controls: constraint convergence, energy conservation, gauge-invariant spectra, ultraviolet refinement, a genuine flux plateau, and variation across initial-data families. It remains open how much of the observed rapid hydrodynamization survives finite coupling and finite N, and which cascade or attractor properties persist in fully three-dimensional inhomogeneous flows. Current evidence does not establish a model-independent universal ordering of all equilibration diagnostics.

The conclusions could change when an independent nonlinear simulation supplies a new counterexample initial state, when a theorem or analytic calculation controls a wider transport class, or when finite-coupling data change the hierarchy of nonhydrodynamic modes. A reproducible higher-resolution cascade with a longer gauge-invariant inertial range, or a change in field content or boundary condition that introduces a new slow sector, could also change the interpretation. A new plot without constraint, cutoff, and constitutive-residual convergence does not.

Compute a Kubo ratio and state its scope. Suppose a renormalized shear correlator has Im⁡GRxy,xy=−0.8ω+O(ω2)\operatorname{Im}G_R^{xy,xy}=-0.8\omega+O(\omega^2) and the entropy density is s=3.2πs=3.2\pi. Find η/s\eta/s. Does agreement with the Einstein value establish universal transport?

Solution

The Kubo slope gives η=0.8\eta=0.8. Dividing by s=3.2πs=3.2\pi yields η/s=1/(4π)\eta/s=1/(4\pi). This checks the normalization and is consistent with an isotropic two-derivative Einstein dual. It does not establish a universal bound or show that higher-derivative, finite-density, translation-breaking, or anisotropic corrections vanish. One must inspect the effective action, symmetries, order of limits, and any additional slow modes before exporting the value to another theory or regime.

  • Adams, Allan, Paul M. Chesler, and Hong Liu. “Holographic Turbulence.” Physical Review Letters 112 (2014): 151602. DOI. Open PDF.
  • Balasubramanian, Vijay, et al. “Holographic Thermalization.” Physical Review D 84 (2011): 026010. DOI. Open PDF.
  • Bhattacharyya, Sayantani, et al. “Nonlinear Fluid Dynamics from Gravity.” Journal of High Energy Physics 02 (2008): 045. DOI. Open PDF.
  • Policastro, Giuseppe, Dam T. Son, and Andrei O. Starinets. “The Shear Viscosity of Strongly Coupled N=4 Supersymmetric Yang–Mills Plasma.” Physical Review Letters 87 (2001): 081601. DOI. Open PDF.
  • Son, Dam T., and Andrei O. Starinets. “Minkowski-Space Correlators in AdS/CFT Correspondence: Recipe and Applications.” Journal of High Energy Physics 09 (2002): 042. DOI. Open PDF.