Skip to content

Thermodynamics of the N=4 Plasma

The planar AdS5_5 black brane is dual to the thermal state of N=4\mathcal N=4 super Yang-Mills theory on flat space R3,1\mathbb R^{3,1}. This page computes its equilibrium thermodynamics.

The result is famous enough to memorize, but it is much more useful to understand where every factor comes from. In the classical supergravity regime,

N,λ=gYM2N,N\to\infty, \qquad \lambda=g_{\mathrm{YM}}^2N\to\infty,

the strongly coupled plasma has

s=π22N2T3,p=π28N2T4,ϵ=3p=3π28N2T4.s={\pi^2\over2}N^2T^3, \qquad p={\pi^2\over8}N^2T^4, \qquad \epsilon=3p={3\pi^2\over8}N^2T^4.

Here ss is the entropy density, pp is the pressure, and ϵ\epsilon is the energy density. The free-energy density is

F=FV3=p=π28N2T4.\mathcal F={F\over V_3}=-p=-{\pi^2\over8}N^2T^4.

The entire computation is a clean example of the thermal AdS/CFT dictionary.

Boundary quantityBulk quantity
temperature TTsmoothness period of Euclidean time, or surface gravity
entropy density sshorizon area density divided by 4G54G_5
free-energy density F\mathcal Frenormalized Euclidean on-shell action density
stress tensor Tμν\langle T_{\mu\nu}\ranglenormalizable metric coefficient in Fefferman-Graham form
O(N2)O(N^2) degrees of freedomclassical bulk action scale L3/G5N2L^3/G_5\sim N^2

The AdS5 black brane computes the entropy density, pressure, energy density, and free energy of strongly coupled N=4 SYM plasma.

The planar AdS5_5 black brane turns horizon data into the equation of state of the strongly coupled N=4\mathcal N=4 SYM plasma. The horizon location fixes TT, the horizon area gives ss, the Euclidean on-shell action gives F=p\mathcal F=-p, and conformality gives ϵ=3p\epsilon=3p.

The word “plasma” can mislead. At weak coupling one imagines a gas of gluons, scalars, and fermions. At strong coupling and large NN, the theory is better described by a smooth black brane. The thermodynamic variables remain perfectly sharp, but they are not describing long-lived quasiparticles.

Use Poincaré coordinates for the five-dimensional AdS black brane:

ds52=L2z2[f(z)dt2+dx2+dz2f(z)],f(z)=1(zzh)4.ds_5^2 = {L^2\over z^2} \left[ -f(z)dt^2+d\vec x^{\,2}+{dz^2\over f(z)} \right], \qquad f(z)=1-\left({z\over z_h}\right)^4.

The conformal boundary is at z=0z=0, and the horizon is at z=zhz=z_h. In Lorentzian signature, z=zhz=z_h is a regular future event horizon after changing to ingoing Eddington-Finkelstein coordinates. In Euclidean signature, regularity at z=zhz=z_h fixes the periodicity of Euclidean time and hence the temperature.

Near z=zhz=z_h, write z=zhρz=z_h-\rho with ρzh\rho\ll z_h. Since

f(z)=1(1ρzh)4=4ρzh+O(ρ2),f(z) = 1-\left(1-{\rho\over z_h}\right)^4 = {4\rho\over z_h}+O(\rho^2),

the Euclidean (tE,z)(t_E,z) part of the metric becomes a smooth plane only if tEt_E has period

β=1T=πzh.\beta={1\over T}=\pi z_h.

Thus

T=1πzh.T={1\over\pi z_h}.

Equivalently, with r=L2/zr=L^2/z, the horizon is at r=rh=L2/zhr=r_h=L^2/z_h, and

T=rhπL2.T={r_h\over\pi L^2}.

The two radial conventions are both common. The zz coordinate is convenient near the boundary; the rr coordinate makes the horizon radius look more like an ordinary black-hole radius.

At fixed Lorentzian time and at the horizon, the induced spatial metric in the three boundary directions is

dsH2=L2zh2dx2.ds_{\mathcal H}^2={L^2\over z_h^2}d\vec x^{\,2}.

Therefore the horizon area per unit boundary volume is

AHV3=L3zh3.{A_{\mathcal H}\over V_3}={L^3\over z_h^3}.

The Bekenstein-Hawking entropy density is

s=1V3AH4G5=L34G5zh3.s={1\over V_3}{A_{\mathcal H}\over4G_5} ={L^3\over4G_5z_h^3}.

For type IIB string theory on AdS5×S5\mathrm{AdS}_5\times S^5,

L3G5=2N2π.{L^3\over G_5}={2N^2\over\pi}.

Using zh=1/(πT)z_h=1/(\pi T) gives

s=142N2π(πT)3=π22N2T3.s ={1\over4}{2N^2\over\pi}(\pi T)^3 ={\pi^2\over2}N^2T^3.

This is the sharpest version of the slogan “black-hole area counts field-theory entropy.” The black-brane horizon is extended along the three boundary spatial directions, so its area is proportional to V3V_3. The entropy is extensive in the boundary volume but controlled by an area in the bulk.

The N2N^2 scaling comes from 1/G51/G_5. In the CFT, the same scaling is the number of adjoint color degrees of freedom. This is not a coincidence; in holographic CFTs, the stress-tensor central charge is proportional to L3/G5L^3/G_5.

Since the boundary theory on flat space is conformal, the thermal stress tensor must be traceless:

Tμμ=ϵ+3p=0.\langle T^\mu{}_{\mu}\rangle=-\epsilon+3p=0.

Thus

ϵ=3p.\epsilon=3p.

The thermodynamic identity

s=pTs={\partial p\over\partial T}

then integrates to

p(T)=0TdTs(T)=π28N2T4.p(T)=\int_0^T dT'\,s(T')={\pi^2\over8}N^2T^4.

The integration constant is zero for the homogeneous plasma on R3\mathbb R^3. Hence

FV3=p=π28N2T4,{F\over V_3}=-p=-{\pi^2\over8}N^2T^4,

and

ϵ=3p=3π28N2T4.\epsilon=3p={3\pi^2\over8}N^2T^4.

Equivalently, the holographic stress tensor gives the mixed-index expectation value

Tμν=π2N2T48diag(3,1,1,1).\langle T^\mu{}_{\nu}\rangle ={\pi^2N^2T^4\over8} \operatorname{diag}(-3,1,1,1).

With lower indices in the rest frame, this is Tμν=diag(ϵ,p,p,p)T_{\mu\nu}=\operatorname{diag}(\epsilon,p,p,p).

Several thermodynamic quantities follow immediately:

w=ϵ+p=Ts=π22N2T4,w=\epsilon+p=Ts={\pi^2\over2}N^2T^4, CV=ϵT=3π22N2T3=3s,C_V={\partial\epsilon\over\partial T} ={3\pi^2\over2}N^2T^3=3s, cs2=pϵ=13.c_s^2={\partial p\over\partial\epsilon}={1\over3}.

The value cs2=1/3c_s^2=1/3 is fixed by conformal invariance in four spacetime dimensions. What holography computes is the overall normalization of the equation of state at strong coupling.

The same free energy can be found directly from the renormalized Euclidean gravitational action. Schematically,

ZCFT(β)=ZIIB[Sβ1×R3]exp(IEren[gbb]),Z_{\mathrm{CFT}}(\beta) = Z_{\mathrm{IIB}}[S^1_\beta\times\mathbb R^3] \simeq \exp\left(-I_E^{\mathrm{ren}}[g_{\mathrm{bb}}]\right),

so

F=TIEren.F=T I_E^{\mathrm{ren}}.

The five-dimensional Euclidean action is

IE=116πG5d5xg(R+12L2)18πG5Md4xγK+Ict.I_E = -{1\over16\pi G_5}\int d^5x\sqrt g\left(R+{12\over L^2}\right) -{1\over8\pi G_5}\int_{\partial M}d^4x\sqrt\gamma\,K +I_{\mathrm{ct}}.

The counterterm action IctI_{\mathrm{ct}} removes UV divergences as the cutoff surface z=ϵz=\epsilon is taken to the boundary. For the planar black brane, the finite renormalized answer is

FV3=L316πG5zh4=π28N2T4.{F\over V_3} =-{L^3\over16\pi G_5z_h^4} =-{\pi^2\over8}N^2T^4.

The entropy follows from

S=FT.S=-{\partial F\over\partial T}.

This agreement between the horizon-area calculation and the Euclidean-action calculation is a valuable normalization check. It tells us that the temperature, ensemble, counterterms, and Newton constant have been matched consistently.

In the boundary theory, all elementary fields are in the adjoint representation, so a thermal gas of weakly coupled fields has O(N2)O(N^2) degrees of freedom. In the bulk, the same scaling comes from

L3G5N2.{L^3\over G_5}\sim N^2.

At large NN, the planar black brane is the dominant homogeneous thermal saddle on R3\mathbb R^3. There is no finite-temperature Hawking-Page transition in the planar case, because the black-brane free-energy density scales as T4-T^4 for every nonzero TT.

This differs from the theory on S3S^3, where the radius of the sphere supplies an additional scale and thermal AdS can dominate at low temperature. In the planar case, conformal invariance leaves no dimensionless ratio except those involving sources or chemical potentials.

The free-field comparison and the factor 3/43/4

Section titled “The free-field comparison and the factor 3/43/43/4”

At zero coupling, N=4\mathcal N=4 SYM is a gas of adjoint massless fields. Per adjoint color generator, the physical degrees of freedom are

2gauge-boson polarizations,6real scalars,4Weyl fermions.2 \quad \text{gauge-boson polarizations}, \qquad 6 \quad \text{real scalars}, \qquad 4 \quad \text{Weyl fermions}.

The bosonic degeneracy is nb=8n_b=8, and the fermionic degeneracy is nf=8n_f=8. Therefore the free pressure is

pfree=π2T490(nb+78nf)(N21)=π26(N21)T4.p_{\mathrm{free}} ={\pi^2T^4\over90}\left(n_b+{7\over8}n_f\right)(N^2-1) ={\pi^2\over6}(N^2-1)T^4.

At large NN,

sfree=pfreeT=2π23N2T3.s_{\mathrm{free}} ={\partial p_{\mathrm{free}}\over\partial T} ={2\pi^2\over3}N^2T^3.

The strong-coupling entropy is

sλ==π22N2T3.s_{\lambda=\infty}={\pi^2\over2}N^2T^3.

Thus

sλ=sλ=0=34.{s_{\lambda=\infty}\over s_{\lambda=0}} ={3\over4}.

This 3/43/4 is famous, but it is not an exact theorem about all holographic plasmas. It compares two different limits of the same large-NN theory:

λ=0andλ=.\lambda=0 \qquad \text{and} \qquad \lambda=\infty.

The striking fact is not merely that 3/43/4 is close to 11. The striking fact is that a classical gravitational area predicts a definite normalization for a strongly coupled, non-BPS thermal quantity.

Finite-coupling and finite-NN corrections

Section titled “Finite-coupling and finite-NNN corrections”

The supergravity answer is the first term in a controlled expansion. For the canonical duality,

L4α2=λ,gsλN.{L^4\over\alpha'^2}=\lambda, \qquad g_s\sim {\lambda\over N}.

Finite λ\lambda corrections come from higher-derivative stringy terms, especially the type IIB α3R4\alpha'^3R^4 interaction. The leading correction to the entropy density has the form

s=π22N2T3[1+158ζ(3)λ3/2+].s ={\pi^2\over2}N^2T^3 \left[ 1+{15\over8}\zeta(3)\lambda^{-3/2}+\cdots \right].

Equivalently,

ssfree=34+4532ζ(3)λ3/2+.{s\over s_{\mathrm{free}}} ={3\over4} +{45\over32}\zeta(3)\lambda^{-3/2} +\cdots.

The positive correction means that as the coupling is decreased from infinity, the entropy begins to move upward toward the free-field value.

Finite NN corrections are quantum-gravity corrections in the bulk. They arise from string loops and are suppressed by powers of 1/N21/N^2 at fixed large λ\lambda. Schematic expansions of thermodynamic quantities therefore look like

FN2T4=classical supergravity+λ3/2(α3 correction)+1N2(bulk loops)+.{F\over N^2T^4} = \text{classical supergravity} + \lambda^{-3/2}(\alpha'^3\text{ correction}) + {1\over N^2}(\text{bulk loops}) + \cdots.

The full type IIB higher-derivative expansion is subtler than this schematic line, but the conceptual lesson is stable: large λ\lambda suppresses string-scale corrections, while large NN suppresses quantum-gravity loops.

Universal versus model-dependent information

Section titled “Universal versus model-dependent information”

Some features of the equation of state are universal for any four-dimensional CFT in a homogeneous thermal state:

ϵ=3p,pT4,sT3,cs2=13.\epsilon=3p, \qquad p\propto T^4, \qquad s\propto T^3, \qquad c_s^2={1\over3}.

These statements follow from conformal invariance and thermodynamics, not from the detailed bulk solution.

Other statements are special to strongly coupled large-NN N=4\mathcal N=4 SYM in the classical gravity regime:

p=π28N2T4,s=π22N2T3,ssfree=34.p={\pi^2\over8}N^2T^4, \qquad s={\pi^2\over2}N^2T^3, \qquad {s\over s_{\mathrm{free}}}={3\over4}.

Those coefficients know about the precise normalization of G5G_5, hence about the S5S^5 compactification and the number of D3-branes.

Still other quantities become universal only after imposing stronger assumptions on the gravitational dual. The famous shear-viscosity result

ηs=14π{\eta\over s}={1\over4\pi}

is not merely a consequence of conformal invariance. It depends on two-derivative Einstein gravity and will be derived later from real-time response.

The phrase “N=4\mathcal N=4 plasma” is precise but easy to misuse. The theory is supersymmetric at zero temperature, conformal, has adjoint scalars and fermions, and has no running coupling. The thermal state breaks supersymmetry because the fermions are antiperiodic around the thermal circle, but the microscopic theory is still not QCD.

The comparison to the quark-gluon plasma is therefore structural rather than literal. Holography gives a controlled example of a strongly coupled non-Abelian plasma with entropy of order N2N^2, rapid thermal relaxation governed by black-brane quasinormal modes, no weakly coupled quasiparticle expansion at large λ\lambda, and hydrodynamics emerging from long-wavelength gravitational perturbations.

The equilibrium equation of state is only the beginning. The real power of the black-brane geometry appears when one perturbs it: sound modes, shear diffusion, conductivity, spectral functions, and nonequilibrium relaxation all become classical wave problems in the same background.

A common mistake is to say that 3/43/4 is a universal prediction of holography. It is not. It is the ratio between the infinite-coupling and free limits of the canonical large-NN N=4\mathcal N=4 theory.

Another mistake is to identify the black-brane entropy with a microscopic derivation of every thermal microstate. The horizon-area computation gives the thermal entropy in the gravitational saddle. AdS/CFT supplies a nonperturbative unitary definition through the CFT, but the classical area calculation itself is a coarse thermodynamic computation.

A third mistake is to forget the ensemble. The planar black brane computes the canonical ensemble at fixed TT and zero chemical potentials. Charged or rotating black holes compute different ensembles with R-charge chemical potentials or angular velocities.

Finally, do not confuse the planar black brane with the global AdS black hole. On R3\mathbb R^3, the planar black brane dominates for all T>0T>0. On S3S^3, thermal AdS and global AdS-Schwarzschild compete, producing the Hawking-Page transition discussed in the previous page.

Starting from

ds52=L2z2[f(z)dt2+dx2+dz2f(z)],f(z)=1(zzh)4,ds_5^2 = {L^2\over z^2} \left[-f(z)dt^2+d\vec x^{\,2}+{dz^2\over f(z)}\right], \qquad f(z)=1-\left({z\over z_h}\right)^4,

show that

s=π22N2T3.s={\pi^2\over2}N^2T^3.

Use T=1/(πzh)T=1/(\pi z_h) and L3/G5=2N2/πL^3/G_5=2N^2/\pi.

Solution

At the horizon, the induced spatial metric along the three boundary directions is

dsH2=L2zh2dx2.ds_{\mathcal H}^2={L^2\over z_h^2}d\vec x^{\,2}.

Therefore

AHV3=L3zh3.{A_{\mathcal H}\over V_3}={L^3\over z_h^3}.

The Bekenstein-Hawking entropy density is

s=14G5L3zh3.s={1\over4G_5}{L^3\over z_h^3}.

Using zh=1/(πT)z_h=1/(\pi T) and L3/G5=2N2/πL^3/G_5=2N^2/\pi gives

s=142N2π(πT)3=π22N2T3.s={1\over4}{2N^2\over\pi}(\pi T)^3 ={\pi^2\over2}N^2T^3.

Assume conformal invariance in four dimensions and the entropy density

s=π22N2T3.s={\pi^2\over2}N^2T^3.

Derive pp, ϵ\epsilon, F/V3F/V_3, and cs2c_s^2.

Solution

Use

s=pT.s={\partial p\over\partial T}.

Then

p=0TdTπ22N2T3=π28N2T4.p=\int_0^T dT'\,{\pi^2\over2}N^2T'^3 ={\pi^2\over8}N^2T^4.

The free-energy density is

FV3=p=π28N2T4.{F\over V_3}=-p=-{\pi^2\over8}N^2T^4.

Conformal invariance gives

ϵ+3p=0,-\epsilon+3p=0,

so

ϵ=3p=3π28N2T4.\epsilon=3p={3\pi^2\over8}N^2T^4.

Finally,

cs2=pϵ=13.c_s^2={\partial p\over\partial\epsilon}={1\over3}.

Exercise 3: The free-field entropy and the factor 3/43/4

Section titled “Exercise 3: The free-field entropy and the factor 3/43/43/4”

Count the free-field degrees of freedom of N=4\mathcal N=4 SYM and show that

sfree=2π23N2T3s_{\mathrm{free}}={2\pi^2\over3}N^2T^3

at leading order in large NN. Then show that sλ=/sλ=0=3/4s_{\lambda=\infty}/s_{\lambda=0}=3/4.

Solution

Per adjoint color generator, the bosons are two gauge-boson polarizations plus six real scalars:

nb=2+6=8.n_b=2+6=8.

The fermions are four Weyl fermions, each with two helicity states:

nf=8.n_f=8.

The free pressure is

pfree=π2T490(nb+78nf)(N21)=π2T490(8+7)(N21).p_{\mathrm{free}} ={\pi^2T^4\over90}\left(n_b+{7\over8}n_f\right)(N^2-1) ={\pi^2T^4\over90}(8+7)(N^2-1).

Thus

pfree=π26(N21)T4π26N2T4.p_{\mathrm{free}}={\pi^2\over6}(N^2-1)T^4 \simeq {\pi^2\over6}N^2T^4.

Differentiating with respect to TT gives

sfree=2π23N2T3.s_{\mathrm{free}} ={2\pi^2\over3}N^2T^3.

The strong-coupling result is

sλ==π22N2T3.s_{\lambda=\infty}={\pi^2\over2}N^2T^3.

Therefore

sλ=sλ=0=π2/22π2/3=34.{s_{\lambda=\infty}\over s_{\lambda=0}} ={\pi^2/2\over2\pi^2/3} ={3\over4}.

The holographic stress tensor of the planar black brane can be written as

Tμν=Adiag(3,1,1,1).\langle T^\mu{}_{\nu}\rangle =A\operatorname{diag}(-3,1,1,1).

Use the pressure found above to determine AA and check the trace Ward identity.

Solution

For mixed indices in the rest frame,

Tμν=diag(ϵ,p,p,p).\langle T^\mu{}_{\nu}\rangle =\operatorname{diag}(-\epsilon,p,p,p).

Comparing with

Adiag(3,1,1,1)A\operatorname{diag}(-3,1,1,1)

gives

A=p=π28N2T4.A=p={\pi^2\over8}N^2T^4.

The trace is

Tμμ=3A+A+A+A=0,\langle T^\mu{}_{\mu}\rangle=-3A+A+A+A=0,

as required for a CFT on flat space.

Exercise 5: Leading finite-coupling correction

Section titled “Exercise 5: Leading finite-coupling correction”

Assume

s=π22N2T3[1+158ζ(3)λ3/2+].s ={\pi^2\over2}N^2T^3 \left[ 1+{15\over8}\zeta(3)\lambda^{-3/2}+\cdots \right].

Compute the leading strong-coupling expansion of s/sfrees/s_{\mathrm{free}}.

Solution

Since

sfree=2π23N2T3,s_{\mathrm{free}}={2\pi^2\over3}N^2T^3,

we have

ssfree=34[1+158ζ(3)λ3/2+].{s\over s_{\mathrm{free}}} ={3\over4} \left[ 1+{15\over8}\zeta(3)\lambda^{-3/2}+\cdots \right].

Thus

ssfree=34+4532ζ(3)λ3/2+.{s\over s_{\mathrm{free}}} ={3\over4} +{45\over32}\zeta(3)\lambda^{-3/2} +\cdots.