Skip to content

Euclidean Gravity and Free Energy

The previous pages used black holes and black branes as geometries that encode thermal CFT states. This page explains the quantitative rule behind that statement:

ZCFT(β)=TrHΣeβHZbulkE[Sβ1×Σ]exp ⁣[IE,ren].Z_{\mathrm{CFT}}(\beta) = \mathrm{Tr}_{\mathcal H_\Sigma}\, e^{-\beta H} \quad\longleftrightarrow\quad Z_{\mathrm{bulk}}^{E}[S^1_\beta\times \Sigma] \approx \exp\!\left[-I_{E,\mathrm{ren}}\right].

In the saddle approximation,

IE,ren=βF,I_{E,{\rm ren}} = \beta F,

where FF is the Helmholtz free energy of the boundary theory in the ensemble specified by the Euclidean boundary conditions. This is the workhorse behind Hawking–Page transitions, black-brane thermodynamics, holographic pressure, and the comparison of competing saddles.

The conceptual point is simple but powerful: equilibrium thermodynamics in the boundary theory is computed by comparing smooth Euclidean bulk fillings of the same boundary thermal manifold.

Boundary setup: thermal partition functions

Section titled “Boundary setup: thermal partition functions”

For a QFT on a spatial manifold Σ\Sigma, the canonical thermal partition function is

Z(β)=TrHΣeβH,β=1T.Z(\beta)= {\rm Tr}_{\mathcal H_\Sigma} e^{-\beta H}, \qquad \beta=\frac{1}{T}.

The Euclidean path integral representation is a path integral on

Sβ1×Σ,S^1_\beta\times \Sigma,

with periodic boundary conditions for bosons and antiperiodic boundary conditions for fermions around Sβ1S^1_\beta.

If one also turns on chemical potentials, sources, or background fields, the partition function becomes a functional. For example, a source JJ for an operator O\mathcal O and a background metric g(0)ijg_{(0)ij} give

ZCFT[g(0),J;β]=Sβ1×Σ ⁣DΦexp ⁣(SE[Φ;g(0)]JO).Z_{\rm CFT}[g_{(0)},J;\beta] = \int_{S^1_\beta\times \Sigma}\!\mathcal D\Phi\, \exp\!\left(-S_E[\Phi;g_{(0)}]-\int J\mathcal O\right).

In thermal holography, this boundary Euclidean manifold is not merely decorative. It is the asymptotic boundary condition for the Euclidean bulk saddle.

The Euclidean bulk path integral is formally

ZbulkE[g(0),J]=M(g(0),J)DgDϕeIE[g,ϕ].Z_{\rm bulk}^{E}[g_{(0)},J] = \int_{\partial M\sim (g_{(0)},J)} \mathcal Dg\,\mathcal D\phi\, e^{-I_E[g,\phi]}.

In the classical gravity limit this becomes a saddle-point expansion:

ZbulkE[g(0),J]smooth saddles aexp ⁣[IE,ren(a)[g(0),J]].Z_{\rm bulk}^{E}[g_{(0)},J] \approx \sum_{\text{smooth saddles }a} \exp\!\left[-I_{E,{\rm ren}}^{(a)}[g_{(0)},J]\right].

At leading order in large NN, the dominant saddle is the one with the smallest renormalized Euclidean action. Thus

F=TlogZTIE,rendominant.F = - T\log Z \approx T\,I_{E,{\rm ren}}^{\rm dominant}.

Equivalently,

IE,rendominant=βF.I_{E,{\rm ren}}^{\rm dominant}=\beta F.

Euclidean gravity and free energy workflow

The Euclidean thermal CFT partition function is represented by a bulk filling whose conformal boundary contains the same thermal circle. In the saddle approximation, the renormalized Euclidean action gives IE,ren=βFI_{E,{\rm ren}}=\beta F; thermodynamic derivatives of IE,renI_{E,{\rm ren}} give EE, SS, pressure, and charge densities.

For pure Einstein gravity with negative cosmological constant, a common Euclidean convention is

IE=116πGd+1Mdd+1xg(R2Λ)18πGd+1MddxγK+Ict.I_E = -\frac{1}{16\pi G_{d+1}} \int_M d^{d+1}x\sqrt g\,(R-2\Lambda) -\frac{1}{8\pi G_{d+1}} \int_{\partial M} d^d x\sqrt\gamma\,K +I_{\rm ct}.

Here:

  • γij\gamma_{ij} is the induced metric on the cutoff boundary;
  • KK is the trace of the extrinsic curvature;
  • IctI_{\rm ct} is a sum of local counterterms on the cutoff surface;
  • Λ=d(d1)/(2L2)\Lambda=-d(d-1)/(2L^2) for AdSd+1_{d+1}.

The counterterms are not optional. The raw on-shell action diverges because the AdS boundary has infinite volume. Holographic renormalization defines

IE,ren=limϵ0(IEzϵ),I_{E,{\rm ren}} = \lim_{\epsilon\to 0} \left( I_E^{z\ge \epsilon}\right),

where IEzϵI_E^{z\ge \epsilon} includes the bulk action, the Gibbons–Hawking–York boundary term, and the counterterms on the cutoff surface z=ϵz=\epsilon.

A useful warning: different sign conventions for the Euclidean action and extrinsic curvature exist in the literature. The invariant statement is the saddle relation

ZEeIE,ren.Z_E \approx e^{-I_{E,{\rm ren}}}.

In the canonical ensemble,

F=IE,renβ,E=IE,renβ,S=βEIE,ren.F = \frac{I_{E,{\rm ren}}}{\beta}, \qquad E = \frac{\partial I_{E,{\rm ren}}}{\partial \beta}, \qquad S = \beta E - I_{E,{\rm ren}}.

Equivalently,

S=(ββ1)IE,ren.S = \left(\beta\frac{\partial}{\partial \beta}-1\right) I_{E,{\rm ren}}.

If the ensemble includes a chemical potential μ\mu for a charge QQ, then the Euclidean source is usually the boundary value of a bulk gauge field component, and the action computes the grand potential:

IE,ren=βΩ,Ω=ETSμQ.I_{E,{\rm ren}}=\beta \Omega, \qquad \Omega = E-TS-\mu Q.

Boundary conditions matter. Fixing the boundary value of AτA_\tau computes a grand-canonical ensemble, while fixing electric flux computes a canonical ensemble. These two choices differ by a boundary Legendre transform.

Consider the Euclidean planar AdSd+1_{d+1} black brane

dsE2=L2z2[f(z)dτ2+dx2+dz2f(z)],f(z)=1(zzh)d.ds_E^2 = \frac{L^2}{z^2} \left[ f(z)d\tau^2+d\vec x^{\,2}+\frac{dz^2}{f(z)} \right], \qquad f(z)=1-\left(\frac{z}{z_h}\right)^d.

Smoothness at the Euclidean horizon requires

β=4πzhd,T=d4πzh.\beta=\frac{4\pi z_h}{d}, \qquad T=\frac{d}{4\pi z_h}.

For this solution, holographic renormalization gives the renormalized Euclidean action density

IE,renβVd1=Ld116πGd+11zhd.\frac{I_{E,{\rm ren}}}{\beta V_{d-1}} = -\frac{L^{d-1}}{16\pi G_{d+1}}\frac{1}{z_h^d}.

Therefore the free energy density is

fthermFVd1=Ld116πGd+11zhd.f_{\rm therm} \equiv \frac{F}{V_{d-1}} = -\frac{L^{d-1}}{16\pi G_{d+1}}\frac{1}{z_h^d}.

The pressure is

p=ftherm=Ld116πGd+11zhd,p=-f_{\rm therm} = \frac{L^{d-1}}{16\pi G_{d+1}}\frac{1}{z_h^d},

and conformal invariance gives

ε=(d1)p=(d1)Ld116πGd+11zhd.\varepsilon=(d-1)p = \frac{(d-1)L^{d-1}}{16\pi G_{d+1}}\frac{1}{z_h^d}.

The entropy density is either the thermodynamic derivative

s=fthermT,s=-\frac{\partial f_{\rm therm}}{\partial T},

or the Bekenstein–Hawking area density

s=14Gd+1(Lzh)d1.s = \frac{1}{4G_{d+1}}\left(\frac{L}{z_h}\right)^{d-1}.

Using T=d/(4πzh)T=d/(4\pi z_h), these agree. The first law also works:

dε=Tds,ε+p=Ts.d\varepsilon = T ds, \qquad \varepsilon + p = Ts.

This is one of the cleanest demonstrations that the horizon is not merely a geometric feature. It carries the entropy of the thermal CFT state.

It is useful to see where the free energy comes from. On shell,

R=d(d+1)L2,R2Λ=2dL2.R=-\frac{d(d+1)}{L^2}, \qquad R-2\Lambda=-\frac{2d}{L^2}.

For the planar black brane,

g=Ld+1zd+1.\sqrt g = \frac{L^{d+1}}{z^{d+1}}.

Thus the regulated bulk action per βVd1\beta V_{d-1} is

IbulkβVd1=dLd18πGd+1ϵzhdzzd+1=Ld18πGd+1(1ϵd1zhd).\frac{I_{\rm bulk}}{\beta V_{d-1}} = \frac{dL^{d-1}}{8\pi G_{d+1}} \int_\epsilon^{z_h}\frac{dz}{z^{d+1}} = \frac{L^{d-1}}{8\pi G_{d+1}} \left(\frac{1}{\epsilon^d}-\frac{1}{z_h^d}\right).

This is divergent. The Gibbons–Hawking–York term and counterterms cancel the ϵd\epsilon^{-d} divergence and adjust the finite term. The final renormalized answer is

IE,renβVd1=Ld116πGd+1zhd.\frac{I_{E,{\rm ren}}}{\beta V_{d-1}} = -\frac{L^{d-1}}{16\pi G_{d+1}z_h^d}.

The lesson is not the numerical factor alone. The lesson is that finite thermodynamics comes from the renormalized on-shell action, not from the bulk volume integral by itself.

For a CFT on Sd1S^{d-1}, the relevant Euclidean boundary is

Sβ1×Sd1.S^1_\beta\times S^{d-1}.

Two important smooth bulk fillings are:

thermal global AdS,Euclidean AdS-Schwarzschild black hole.\text{thermal global AdS}, \qquad \text{Euclidean AdS-Schwarzschild black hole}.

The global AdS-Schwarzschild metric has

dsE2=f(r)dτ2+dr2f(r)+r2dΩd12,ds_E^2 = f(r)d\tau^2+\frac{dr^2}{f(r)}+r^2 d\Omega_{d-1}^2,

with

f(r)=1+r2L2μrd2,μ=rhd2(1+rh2L2).f(r)=1+\frac{r^2}{L^2}-\frac{\mu}{r^{d-2}}, \qquad \mu=r_h^{d-2}\left(1+\frac{r_h^2}{L^2}\right).

The temperature is

T=drh2+(d2)L24πL2rh.T = \frac{d r_h^2+(d-2)L^2}{4\pi L^2 r_h}.

The free energy is

F=Ωd1rhd216πGd+1L2(L2rh2).F = \frac{\Omega_{d-1} r_h^{d-2}}{16\pi G_{d+1}L^2} \left(L^2-r_h^2\right).

So large black holes with rh>Lr_h>L have F<0F<0 and dominate over thermal AdS, while small black holes with rh<Lr_h<L have F>0F>0 and do not dominate the canonical ensemble. The transition occurs at

rh=L,THP=d12πL.r_h=L, \qquad T_{\rm HP}=\frac{d-1}{2\pi L}.

This is the Hawking–Page transition. On the boundary, it is interpreted as a large-NN thermal phase transition on compact space.

The formula

S=(ββ1)IES=\left(\beta\frac{\partial}{\partial\beta}-1\right)I_E

reproduces the Bekenstein–Hawking entropy. One intuitive derivation uses a conical defect.

Near a Euclidean horizon, the metric locally looks like

ds2dR2+R2dθ2+dsH2.ds^2 \approx dR^2 + R^2 d\theta^2 + ds_{\mathcal H}^2.

Smoothness requires θθ+2π\theta\sim \theta+2\pi. If instead the angular period is 2πα2\pi\alpha, the origin has a conical defect with localized curvature. Evaluating the derivative of the action with respect to α\alpha at α=1\alpha=1 gives

S=Area(H)4Gd+1.S=\frac{\mathrm{Area}(\mathcal H)}{4G_{d+1}}.

In higher-derivative gravity, this result is replaced by Wald entropy or its appropriate generalization. In two-derivative Einstein gravity, it is the area law.

At leading large NN, the saddle with smallest IEI_E dominates:

ZeIEmin.Z\approx e^{-I_{E}^{\rm min}}.

If two saddles exchange dominance as β\beta changes, the boundary theory has a large-NN phase transition. The Hawking–Page transition is the canonical example:

ΔIE(β)=IEblack holeIEthermal AdS.\Delta I_E(\beta)=I_E^{\rm black\ hole}-I_E^{\rm thermal\ AdS}.

Then

ΔIE<0black hole dominates,\Delta I_E<0 \quad\Rightarrow\quad \text{black hole dominates},

while

ΔIE>0thermal AdS dominates.\Delta I_E>0 \quad\Rightarrow\quad \text{thermal AdS dominates}.

At finite NN, the transition is rounded or corrected by subleading saddles and quantum fluctuations. In the strict classical gravity limit, the saddle competition becomes sharp.

Boundary quantityEuclidean bulk quantity
thermal circle Sβ1S^1_\betaasymptotic Euclidean time circle
partition function Z(β)Z(\beta)bulk Euclidean path integral
free energy FFIE,ren/βI_{E,{\rm ren}}/\beta
entropy SSββIEIE\beta\partial_\beta I_E-I_E or horizon area
pressure ppF/V-F/V for homogeneous states
chemical potential μ\muboundary value of Euclidean gauge field
phase transitionexchange of dominant Euclidean saddle

The Euclidean action computes equilibrium thermodynamics. It does not, by itself, compute real-time response. For that we need the Lorentzian prescription.

“The free energy is just the bulk action without boundary terms.”

Section titled ““The free energy is just the bulk action without boundary terms.””

No. The Gibbons–Hawking–York term and counterterms are essential. Without them, the variational problem is not correct and the action is divergent.

“A Euclidean black hole has a Lorentzian horizon.”

Section titled ““A Euclidean black hole has a Lorentzian horizon.””

The Euclidean geometry has a smooth origin where the thermal circle caps off. The Lorentzian horizon is recovered after analytic continuation. The Euclidean cap is why the temperature is fixed.

“All saddles with the same boundary are equally important.”

Section titled ““All saddles with the same boundary are equally important.””

They all contribute in principle, but in the classical large-NN limit the smallest renormalized action dominates. Subdominant saddles are exponentially suppressed by factors of order eN2e^{-N^2} in standard large-NN examples.

“The Euclidean action always computes Helmholtz free energy.”

Section titled ““The Euclidean action always computes Helmholtz free energy.””

Only if the boundary conditions define a canonical ensemble. If sources such as chemical potentials are fixed, the action computes an appropriate thermodynamic potential, usually a grand potential.

“Euclidean correlators are enough for transport.”

Section titled ““Euclidean correlators are enough for transport.””

Euclidean correlators determine Matsubara data. Transport coefficients require real-time retarded correlators and a prescription for horizon boundary conditions. This is the subject of the next page.

Exercise 1: Planar black-brane equation of state

Section titled “Exercise 1: Planar black-brane equation of state”

Starting from

ftherm=Ld116πGd+1zhd,T=d4πzh,f_{\rm therm} = -\frac{L^{d-1}}{16\pi G_{d+1}z_h^d}, \qquad T=\frac{d}{4\pi z_h},

show that ε=(d1)p\varepsilon=(d-1)p and ε+p=Ts\varepsilon+p=Ts.

Solution

The pressure is

p=ftherm=Ld116πGd+1zhd.p=-f_{\rm therm} = \frac{L^{d-1}}{16\pi G_{d+1}z_h^d}.

Conformal invariance gives Tii=0T^i{}_i=0, so for a homogeneous thermal state

ε+(d1)p=0,-\varepsilon+(d-1)p=0,

hence

ε=(d1)p.\varepsilon=(d-1)p.

The entropy density from the horizon area is

s=14Gd+1(Lzh)d1.s=\frac{1}{4G_{d+1}}\left(\frac{L}{z_h}\right)^{d-1}.

Then

Ts=d4πzhLd14Gd+1zhd1=dLd116πGd+1zhd=dp.Ts = \frac{d}{4\pi z_h}\frac{L^{d-1}}{4G_{d+1}z_h^{d-1}} = \frac{dL^{d-1}}{16\pi G_{d+1}z_h^d} =dp.

Since ε+p=(d1)p+p=dp\varepsilon+p=(d-1)p+p=dp, we get

ε+p=Ts.\varepsilon+p=Ts.

Exercise 2: Hawking–Page free energy sign

Section titled “Exercise 2: Hawking–Page free energy sign”

For a global AdSd+1_{d+1} black hole,

F=Ωd1rhd216πGd+1L2(L2rh2).F = \frac{\Omega_{d-1} r_h^{d-2}}{16\pi G_{d+1}L^2} \left(L^2-r_h^2\right).

For which values of rhr_h does the black hole dominate over thermal AdS?

Solution

With the usual normalization in which thermal AdS has F=0F=0, the black hole dominates when F<0F<0. Since the prefactor

Ωd1rhd216πGd+1L2\frac{\Omega_{d-1} r_h^{d-2}}{16\pi G_{d+1}L^2}

is positive, the sign is determined by L2rh2L^2-r_h^2. Therefore

F<0rh>L.F<0 \quad\Longleftrightarrow\quad r_h>L.

The transition occurs at rh=Lr_h=L, corresponding to

THP=d12πL.T_{\rm HP}=\frac{d-1}{2\pi L}.

Suppose a saddle has

IE(β)=βMS0,I_E(\beta)=\beta M-S_0,

where MM and S0S_0 are independent of β\beta along the family considered. Use

S=(ββ1)IES=\left(\beta\partial_\beta-1\right)I_E

to find the entropy.

Solution

First compute

βIE=M.\partial_\beta I_E = M.

Then

(ββ1)IE=βM(βMS0)=S0.\left(\beta\partial_\beta-1\right)I_E = \beta M-(\beta M-S_0) =S_0.

For an Einstein black hole, S0S_0 is the Bekenstein–Hawking entropy

S0=Area(H)4G.S_0=\frac{\mathrm{Area}(\mathcal H)}{4G}.