Skip to content

Eternal Black Holes and the Thermofield Double

The eternal AdS black hole is the cleanest laboratory for separating three ideas that are often blurred together:

  1. a thermal density matrix in one quantum system,
  2. a pure entangled state in two quantum systems,
  3. a two-sided black-hole geometry with an Einstein—Rosen bridge.

The essential statement is

TFDβHLHRmaximally extended AdS black hole\boxed{ |\mathrm{TFD}_\beta\rangle \in \mathcal H_L\otimes \mathcal H_R \quad \longleftrightarrow \quad \text{maximally extended AdS black hole} }

where

TFDβ=1Z(β)neβEn/2nLnR.|\mathrm{TFD}_\beta\rangle = \frac{1}{\sqrt{Z(\beta)}} \sum_n e^{-\beta E_n/2} |n\rangle_L |n\rangle_R .

Here HL\mathcal H_L and HR\mathcal H_R are two independent copies of the same CFT Hilbert space, Hn=EnnH|n\rangle = E_n |n\rangle, and

Z(β)=TreβH.Z(\beta)=\mathrm{Tr}\, e^{-\beta H}.

More carefully, the left state is often the CPT conjugate of the right state, so one writes nLnR|n^*\rangle_L |n\rangle_R. For most computations in this page the distinction only affects how one defines left operators.

This page begins the black-hole information module because the thermofield double is the simplest exact example where a black-hole exterior looks thermal while the complete boundary state is pure. It is not yet the evaporating black-hole problem. It is the equilibrium, two-sided, AdS version of the problem: a place where horizons, entropy, entanglement, and behind-the-horizon geometry can all be discussed with sharp boundary variables.

A two-panel diagram showing the thermofield-double state of two CFTs and its dual two-sided eternal AdS black hole. The Penrose diagram has left and right AdS boundaries, four horizons meeting at a bifurcation surface, and an Einstein-Rosen bridge on the time-symmetric slice.

The thermofield double TFDβ|\mathrm{TFD}_\beta\rangle is a pure entangled state of two non-interacting CFTs. At large NN and high enough temperature for the black-hole saddle to dominate, it is dual to the two-sided eternal AdS black hole. The entropy of either CFT equals the area of the bifurcation surface at leading order, SR=AH/(4GN)S_R=A_{\mathcal H}/(4G_N).

Start with one CFT in a thermal state

ρβ=eβHZ(β).\rho_\beta = \frac{e^{-\beta H}}{Z(\beta)}.

This is a mixed state on a single Hilbert space H\mathcal H. Any mixed state can be purified by embedding it into a larger Hilbert space. The thermofield double is the canonical purification obtained by using a second identical copy:

Hdouble=HLHR.\mathcal H_{\mathrm{double}} = \mathcal H_L\otimes \mathcal H_R.

The TFD state is

TFDβ=1Z(β)neβEn/2nLnR.|\mathrm{TFD}_\beta\rangle = \frac{1}{\sqrt{Z(\beta)}} \sum_n e^{-\beta E_n/2} |n\rangle_L |n\rangle_R.

It is normalized because

TFDβTFDβ=1Z(β)neβEn=1.\langle \mathrm{TFD}_\beta|\mathrm{TFD}_\beta\rangle = \frac{1}{Z(\beta)} \sum_n e^{-\beta E_n} =1.

Tracing out the left CFT gives

ρR=TrLTFDβTFDβ=1Z(β)neβEnnRnR=eβHRZ(β).\rho_R =\mathrm{Tr}_L |\mathrm{TFD}_\beta\rangle\langle \mathrm{TFD}_\beta| = \frac{1}{Z(\beta)} \sum_n e^{-\beta E_n}|n\rangle_R\langle n|_R = \frac{e^{-\beta H_R}}{Z(\beta)}.

Thus every right-sided observable has the ordinary thermal expectation value:

TFDβORTFDβ=TrR(ρβOR).\langle \mathrm{TFD}_\beta|O_R|\mathrm{TFD}_\beta\rangle =\mathrm{Tr}_R(\rho_\beta O_R).

This is the first crucial lesson. A thermal density matrix can arise either because the system is fundamentally in a mixed state, or because it is a subsystem of a larger pure state. In AdS/CFT the latter option is explicit: the right exterior of the eternal black hole is dual to the right CFT alone, while the full two-sided geometry is dual to a pure state in HLHR\mathcal H_L\otimes\mathcal H_R.

The entanglement entropy of one side is

SR=TrρRlogρR=βH+logZ(β).S_R = -\mathrm{Tr}\,\rho_R\log\rho_R = \beta \langle H\rangle + \log Z(\beta).

Equivalently,

SR=(1ββ)logZ(β).S_R = \left(1-\beta\frac{\partial}{\partial \beta}\right)\log Z(\beta).

For a holographic large-NN CFT in the deconfined phase, this entropy is order N2N^2 and is reproduced by the black-hole area.

The TFD state has a simple path-integral construction. The thermal partition function of one CFT is the Euclidean path integral on a circle of circumference β\beta:

Z(β)=TreβH.Z(\beta) =\mathrm{Tr}\,e^{-\beta H}.

Cut this circle into two equal intervals of Euclidean length β/2\beta/2. The path integral on the half-circle prepares a state in two copies of the Hilbert space. In the energy basis it gives precisely

nLmRTFDβneβH/2m=eβEm/2δnm,\langle n|_L \langle m|_R |\mathrm{TFD}_\beta\rangle \propto \langle n|e^{-\beta H/2}|m\rangle =e^{-\beta E_m/2}\delta_{nm},

up to the left-right conjugation convention.

This Euclidean construction is also the cleanest way to see the bulk geometry. In the high-temperature saddle of global AdS, or for planar black branes at any nonzero temperature, the Euclidean bulk is a smooth black-hole cigar: the Euclidean time circle smoothly caps off in the interior. Cutting the Euclidean cigar across a time-reflection-symmetric slice and analytically continuing gives the two-sided Lorentzian black hole.

So the TFD is not an arbitrary entangled state. It is prepared by Euclidean evolution for a time β/2\beta/2. This strongly constrains its entanglement structure, and that special structure is what allows a simple semiclassical wormhole interpretation.

A convenient representative of the eternal AdS-Schwarzschild family is

ds2=f(r)dt2+dr2f(r)+r2dΩd12,ds^2 = -f(r)dt^2 +\frac{dr^2}{f(r)} +r^2 d\Omega_{d-1}^2,

with

f(r)=1+r2L2μrd2.f(r)=1+\frac{r^2}{L^2}-\frac{\mu}{r^{d-2}}.

For the planar black brane one replaces dΩd12d\Omega_{d-1}^2 by dxd12d\vec x_{d-1}^{\,2} and uses the planar form of f(r)f(r). The horizon is at

f(rh)=0.f(r_h)=0.

The exterior metric above only covers one asymptotic region. Its maximal analytic extension has two AdS boundaries, two exterior regions, a future interior, a past interior, and a bifurcation surface where the future and past horizons meet. Each asymptotic boundary has its own CFT:

right boundaryCFTR,left boundaryCFTL.\text{right boundary}\leftrightarrow \mathrm{CFT}_R, \qquad \text{left boundary}\leftrightarrow \mathrm{CFT}_L.

The two CFTs are not coupled. The Hamiltonian is simply

Hdouble=HL+HR.H_{\mathrm{double}}=H_L+H_R.

The wormhole is therefore not a channel through which the two field theories exchange signals. It is a geometric representation of entanglement in a particular state of two decoupled systems.

At tL=tR=0t_L=t_R=0, the time-symmetric spatial slice connects the two exteriors through an Einstein—Rosen bridge. The bifurcation surface has area AHA_{\mathcal H}, and the Bekenstein-Hawking entropy is

SBH=AH4Gd+1.S_{\mathrm{BH}} =\frac{A_{\mathcal H}}{4G_{d+1}}.

From the boundary viewpoint this is the thermal entropy of one CFT, which is also the entanglement entropy between the two CFTs:

S(ρR)=S(ρL)=SBH+O(N0).S(\rho_R)=S(\rho_L)=S_{\mathrm{BH}}+O(N^0).

The O(N0)O(N^0) term is the bulk quantum correction discussed in the previous module: the leading area term is corrected by bulk entanglement and counterterm contributions.

The RT/HRT formula gives a particularly sharp interpretation of the area law. Take the boundary region AA to be the entire right boundary at the time-symmetric slice. The complementary boundary region is the entire left boundary. The extremal surface homologous to the right boundary is the bifurcation surface H\mathcal H.

Thus, at leading order,

SR=Area(H)4GN.S_R =\frac{\mathrm{Area}(\mathcal H)}{4G_N}.

This is one of the most beautiful checks of the TFD dictionary:

thermal entropy of one CFT=entanglement entropy between the two CFTs=horizon area in Planck units\boxed{ \text{thermal entropy of one CFT} = \text{entanglement entropy between the two CFTs} = \text{horizon area in Planck units} }

There is a useful conceptual distinction here. The total state TFD|\mathrm{TFD}\rangle is pure, so

SLR=0.S_{LR}=0.

But either side alone is mixed, so

SL=SR=Sthermal.S_L=S_R=S_{\mathrm{thermal}}.

The mutual information between the two complete CFTs is therefore

I(L:R)=SL+SRSLR=2Sthermal.I(L:R)=S_L+S_R-S_{LR}=2S_{\mathrm{thermal}}.

This large mutual information is the boundary signature of the connected two-sided geometry.

The TFD state has a simple but important symmetry. Since the same energy eigenvalue appears on both sides,

(HLHR)TFDβ=0.(H_L-H_R)|\mathrm{TFD}_\beta\rangle=0.

Therefore the state is invariant under

eia(HLHR)TFDβ=TFDβ.e^{-ia(H_L-H_R)}|\mathrm{TFD}_\beta\rangle =|\mathrm{TFD}_\beta\rangle.

This is the boundary version of the boost Killing symmetry of the eternal black hole. In the Penrose diagram this Killing vector is future-directed on the right exterior and past-directed on the left exterior. That is the source of a common sign confusion: the geometric time coordinate in the left exterior is naturally oriented oppositely to the physical future time of the left CFT.

Now evolve both CFTs forward with their ordinary Hamiltonians:

TFDβ(t)=eit(HL+HR)TFDβ.|\mathrm{TFD}_\beta(t)\rangle =e^{-it(H_L+H_R)}|\mathrm{TFD}_\beta\rangle.

In the energy basis,

TFDβ(t)=1Z(β)neβEn/2e2iEntnLnR.|\mathrm{TFD}_\beta(t)\rangle = \frac{1}{\sqrt{Z(\beta)}} \sum_n e^{-\beta E_n/2}e^{-2iE_n t} |n\rangle_L |n\rangle_R.

This state is not equal to the original TFD. However, the reduced density matrix on either side is unchanged:

ρR(t)=TrLTFDβ(t)TFDβ(t)=eβHRZ(β).\rho_R(t) =\mathrm{Tr}_L |\mathrm{TFD}_\beta(t)\rangle\langle \mathrm{TFD}_\beta(t)| =\frac{e^{-\beta H_R}}{Z(\beta)}.

So one-sided observables remain thermal and time-translation invariant, while two-sided observables change. In the bulk this corresponds to the fact that the exterior regions remain stationary, but the Einstein—Rosen bridge grows with the two-sided time. This is a first glimpse of a major modern theme: black-hole interiors can encode time-dependent information that is invisible in strictly one-sided thermal observables.

For a right-sided operator, the TFD gives ordinary thermal correlators:

TFDβOR(t)OR(0)TFDβ=Tr(ρβO(t)O(0)).\langle \mathrm{TFD}_\beta|O_R(t)O_R(0)|\mathrm{TFD}_\beta\rangle = \mathrm{Tr}\left(\rho_\beta O(t)O(0)\right).

The right exterior of the black hole is enough to compute such correlators. This is the Lorentzian continuation of the finite-temperature dictionary studied in the previous modules.

Two-sided correlators are more interesting. Schematically,

TFDβOL(tL)OR(tR)TFDβ\langle \mathrm{TFD}_\beta|O_L(t_L)O_R(t_R)|\mathrm{TFD}_\beta\rangle

measures correlations between the two boundaries. With the standard left-right conjugation convention, this correlator is equivalent to a thermal correlator with an imaginary time shift by β/2\beta/2:

OL(tL)OR(tR)TFDTr(ρβO(tR)O(tL+iβ/2)),\langle O_L(t_L)O_R(t_R)\rangle_{\mathrm{TFD}} \sim \mathrm{Tr}\left( \rho_\beta\, O(t_R)O(t_L+i\beta/2) \right),

up to operator-ordering and conjugation conventions. The shift by iβ/2i\beta/2 is the algebraic trace of the Euclidean half-circle that prepares the TFD.

In the bulk, heavy-operator two-sided correlators can be approximated by spacelike geodesics passing through the Einstein—Rosen bridge. More generally, two-sided correlators are sensitive to the wormhole geometry in a way that one-sided thermal correlators are not.

This is why the eternal black hole is a useful information-theoretic laboratory. The one-sided observer sees thermal physics; the doubled system contains the purification; and two-sided observables know about the entanglement that purifies the thermal state.

The two-sided eternal AdS black hole contains a wormhole, but it is not traversable. A causal signal sent from the right boundary into the black hole cannot emerge at the left boundary. It hits the future singularity.

The boundary explanation is simple. The two CFTs are decoupled:

Hdouble=HL+HR.H_{\mathrm{double}}=H_L+H_R.

Therefore a right-sided perturbation cannot causally affect a left-sided observable. For tL,tRt_L,t_R in a causal ordering where no explicit left-right interaction has been turned on,

[OL(tL),OR(tR)]=0.[O_L(t_L),O_R(t_R)]=0.

Entanglement creates correlations, not communication. This distinction is absolutely central. A Bell pair has correlations, but it does not transmit a signal. The eternal black hole is the gravitational version of the same principle, except that in the large-NN and strong-coupling regime the correlations are geometrized as an Einstein—Rosen bridge.

There are controlled deformations that can make related wormholes traversable. The most famous construction couples the two boundaries by a brief double-trace interaction of the schematic form

δH(t)=h(t)OL(t)OR(t),\delta H(t)=h(t)\,O_L(t)O_R(t),

which can generate negative averaged null energy in the bulk and open a causal window through the wormhole. That is a deformation of the system. It is not a property of the undeformed TFD state.

Low temperature, high temperature, and saddle dominance

Section titled “Low temperature, high temperature, and saddle dominance”

The TFD state exists for every β\beta in any quantum system. The statement

TFDβtwo-sided AdS black hole|\mathrm{TFD}_\beta\rangle \longleftrightarrow \text{two-sided AdS black hole}

is a statement about the dominant semiclassical bulk saddle in a suitable regime.

For a CFT on Sd1S^{d-1} at large NN, low-temperature thermal physics may be dominated by thermal AdS rather than by an AdS-Schwarzschild black hole. At high temperature the large black hole dominates, and the two-sided black-hole saddle is the natural bulk dual of the TFD. This is the Lorentzian counterpart of the Hawking—Page discussion in the finite-temperature module.

For a planar black brane, there is no finite-temperature Hawking—Page transition in the same sense; the black brane describes the deconfined thermal plasma for any nonzero temperature in the infinite-volume limit.

The moral is simple: the TFD is exact, but the semiclassical geometry it is approximated by depends on the phase and on the large-NN saddle.

Eternal black hole versus black-hole microstate

Section titled “Eternal black hole versus black-hole microstate”

The TFD is not a typical pure state of one CFT. It is a pure state of two CFTs:

TFDHLHR.|\mathrm{TFD}\rangle\in \mathcal H_L\otimes\mathcal H_R.

A one-sided black-hole microstate, by contrast, is a state

ΨHR|\Psi\rangle\in \mathcal H_R

with energy of order that of a black hole. Such a state may reproduce thermal one-point functions and simple two-point functions for a long time, but it is not literally the TFD and it does not have a second asymptotic boundary.

This distinction matters for the information problem. The eternal black hole shows how a black-hole exterior can look thermal while the complete quantum state remains pure. But it does not by itself solve the problem of an evaporating black hole formed from collapse. For evaporation, one must understand time-dependent geometry, Hawking radiation, Page curves, islands, and the fine-grained entropy of radiation. Those topics appear later in this module.

Boundary statementBulk statement
Two decoupled CFTs, HLHR\mathcal H_L\otimes\mathcal H_RTwo asymptotic AdS boundaries
$\mathrm{TFD}_\beta\rangle$
$\rho_R=\mathrm{Tr}_L\mathrm{TFD}\rangle\langle\mathrm{TFD}
S(ρR)S(\rho_R)Horizon area AH/(4GN)A_{\mathcal H}/(4G_N) at leading order
HLHRH_L-H_R symmetryBoost Killing symmetry of the eternal geometry
Evolution by HL+HRH_L+H_RGrowth of the Einstein—Rosen bridge
One-sided correlatorsExterior black-hole perturbation theory
Two-sided correlatorsProbes of the bridge and behind-horizon geometry
No left-right couplingNon-traversable wormhole
Double-trace left-right couplingPossible traversable deformation in special setups

Mistake 1: thinking the TFD is just a thermal state

Section titled “Mistake 1: thinking the TFD is just a thermal state”

The right CFT is thermal after tracing out the left CFT, but the doubled state is pure:

SLR=0,SR=SL=Sthermal.S_{LR}=0, \qquad S_R=S_L=S_{\mathrm{thermal}}.

Confusing the purified state with the reduced state is the most common conceptual error.

Mistake 2: treating the two CFTs as interacting

Section titled “Mistake 2: treating the two CFTs as interacting”

The two CFTs in the undeformed TFD setup are decoupled. The connected bulk geometry is not produced by an interaction term in the boundary Hamiltonian. It is produced by entanglement in the state.

Mistake 3: calling the wormhole traversable

Section titled “Mistake 3: calling the wormhole traversable”

The eternal black-hole wormhole is non-traversable. Entanglement gives correlations, not signals. Traversability requires additional physics, such as a carefully chosen coupling between the two boundaries.

Mistake 4: ignoring the left-time sign convention

Section titled “Mistake 4: ignoring the left-time sign convention”

The TFD is invariant under HLHRH_L-H_R, not under HL+HRH_L+H_R. Geometrically, the eternal black-hole Killing time is future-directed on one boundary and past-directed on the other.

Mistake 5: identifying every entangled state with a smooth wormhole

Section titled “Mistake 5: identifying every entangled state with a smooth wormhole”

The TFD has a very special entanglement pattern: energy eigenstates are paired with Boltzmann weights. Generic entangled states of two CFTs do not automatically have a smooth semiclassical Einstein—Rosen bridge.

Mistake 6: confusing the eternal black hole with an evaporating black hole

Section titled “Mistake 6: confusing the eternal black hole with an evaporating black hole”

The eternal AdS black hole is an equilibrium object. It is not the same as a black hole formed from collapse and then evaporating into an external bath. It is an essential warm-up, not the whole information paradox.

Starting from

TFDβ=1ZneβEn/2nLnR,|\mathrm{TFD}_\beta\rangle = \frac{1}{\sqrt{Z}} \sum_n e^{-\beta E_n/2}|n\rangle_L|n\rangle_R,

show explicitly that

ρR=TrLTFDβTFDβ=eβHRZ.\rho_R=\mathrm{Tr}_L |\mathrm{TFD}_\beta\rangle\langle \mathrm{TFD}_\beta| =\frac{e^{-\beta H_R}}{Z}.

Then show that

SR=βH+logZ.S_R=\beta\langle H\rangle +\log Z.
Solution

The projector is

TFDTFD=1Zm,neβ(Em+En)/2mLmRLnRn.|\mathrm{TFD}\rangle\langle \mathrm{TFD}| = \frac{1}{Z} \sum_{m,n} e^{-\beta(E_m+E_n)/2} |m\rangle_L|m\rangle_R {}_L\langle n|{}_R\langle n|.

Tracing over the left Hilbert space gives

ρR=1Zm,neβ(Em+En)/2LnmLmRRn.\rho_R = \frac{1}{Z} \sum_{m,n} e^{-\beta(E_m+E_n)/2} {}_L\langle n|m\rangle_L |m\rangle_R{}_R\langle n|.

Using LnmL=δnm{}_L\langle n|m\rangle_L=\delta_{nm},

ρR=1ZneβEnnRRn=eβHRZ.\rho_R = \frac{1}{Z}\sum_n e^{-\beta E_n}|n\rangle_R{}_R\langle n| =\frac{e^{-\beta H_R}}{Z}.

The entropy is

SR=TrρRlogρR.S_R =-\mathrm{Tr}\,\rho_R\log\rho_R.

Since

logρR=βHRlogZ,\log \rho_R=-\beta H_R-\log Z,

we get

SR=βTr(ρRHR)+logZ=βH+logZ.S_R =\beta\,\mathrm{Tr}(\rho_R H_R)+\log Z =\beta\langle H\rangle +\log Z.

Exercise 2: Which Hamiltonian leaves the TFD invariant?

Section titled “Exercise 2: Which Hamiltonian leaves the TFD invariant?”

Show that

(HLHR)TFDβ=0,(H_L-H_R)|\mathrm{TFD}_\beta\rangle=0,

but that HL+HRH_L+H_R does not annihilate the TFD unless the state has support only at zero energy.

Solution

Acting on a basis term gives

HLnLnR=EnnLnR,H_L |n\rangle_L|n\rangle_R =E_n |n\rangle_L|n\rangle_R,

and

HRnLnR=EnnLnR.H_R |n\rangle_L|n\rangle_R =E_n |n\rangle_L|n\rangle_R.

Therefore

(HLHR)nLnR=0(H_L-H_R)|n\rangle_L|n\rangle_R=0

term by term in the TFD sum, so

(HLHR)TFD=0.(H_L-H_R)|\mathrm{TFD}\rangle=0.

On the other hand,

(HL+HR)nLnR=2EnnLnR.(H_L+H_R)|n\rangle_L|n\rangle_R =2E_n |n\rangle_L|n\rangle_R.

Thus

(HL+HR)TFD=1Zn2EneβEn/2nLnR,(H_L+H_R)|\mathrm{TFD}\rangle =\frac{1}{\sqrt Z}\sum_n 2E_n e^{-\beta E_n/2}|n\rangle_L|n\rangle_R,

which is generally nonzero. Evolution by HL+HRH_L+H_R changes the phases of the TFD components while leaving the one-sided density matrices thermal.

Exercise 3: Euclidean half-circle preparation

Section titled “Exercise 3: Euclidean half-circle preparation”

Use the matrix element of Euclidean time evolution to show that a Euclidean path integral over an interval of length β/2\beta/2 prepares the TFD state.

Solution

The Euclidean evolution operator is

eβH/2.e^{-\beta H/2}.

Insert complete sets of energy eigenstates at the two ends of the interval:

neβH/2m=eβEm/2δnm.\langle n|e^{-\beta H/2}|m\rangle =e^{-\beta E_m/2}\delta_{nm}.

Interpreting the two ends of the interval as the two Hilbert-space factors gives a state proportional to

meβEm/2mLmR.\sum_m e^{-\beta E_m/2}|m\rangle_L|m\rangle_R.

Normalizing by Z(β)1/2Z(\beta)^{-1/2} gives

TFDβ=1Z(β)meβEm/2mLmR.|\mathrm{TFD}_\beta\rangle =\frac{1}{\sqrt{Z(\beta)}} \sum_m e^{-\beta E_m/2}|m\rangle_L|m\rangle_R.

Depending on the orientation of the Euclidean path integral, the left factor may be represented by the conjugate state mL|m^*\rangle_L.

Consider the time-symmetric slice of the two-sided eternal black hole. Let AA be the entire right boundary. Explain why the RT surface computing S(A)S(A) is the bifurcation surface of the horizon.

Solution

The RT surface must be codimension two, extremal on the time-symmetric slice, and homologous to the boundary region AA. Since AA is the entire right boundary, the bulk region homologous to AA is the entire right exterior plus the part of the bridge ending at the bifurcation surface. The surface that separates the right exterior from the left exterior is precisely the bifurcation surface H\mathcal H.

On the time-symmetric slice it is minimal by symmetry and by the standard black-hole area theorem intuition for the stationary horizon cross-section. Therefore

S(A)=Area(H)4GNS(A)=\frac{\mathrm{Area}(\mathcal H)}{4G_N}

at leading order. This equals the thermal entropy of the right CFT and the entanglement entropy between the two CFTs.

Exercise 5: One-sided observables remain thermal under two-sided time evolution

Section titled “Exercise 5: One-sided observables remain thermal under two-sided time evolution”

Let

TFD(t)=eit(HL+HR)TFD.|\mathrm{TFD}(t)\rangle=e^{-it(H_L+H_R)}|\mathrm{TFD}\rangle.

Show that the reduced density matrix ρR(t)\rho_R(t) is independent of tt.

Solution

The evolved state is

TFD(t)=1ZneβEn/2e2iEntnLnR.|\mathrm{TFD}(t)\rangle =\frac{1}{\sqrt Z}\sum_n e^{-\beta E_n/2}e^{-2iE_n t}|n\rangle_L|n\rangle_R.

Then

TFD(t)TFD(t)=1Zm,neβ(Em+En)/2e2it(EmEn)mLmRLnRn.|\mathrm{TFD}(t)\rangle\langle\mathrm{TFD}(t)| = \frac{1}{Z}\sum_{m,n} e^{-\beta(E_m+E_n)/2}e^{-2it(E_m-E_n)} |m\rangle_L|m\rangle_R{}_L\langle n|{}_R\langle n|.

Tracing over the left side sets m=nm=n, so the phase cancels:

ρR(t)=1ZneβEnnRRn=ρR(0).\rho_R(t) =\frac{1}{Z}\sum_n e^{-\beta E_n}|n\rangle_R{}_R\langle n| =\rho_R(0).

Thus all right-sided one-point functions and density-matrix observables remain thermal, even though two-sided correlators can depend on tt.