Skip to content

Unruh Temperature and Thermal Periodicity

The previous page rewrote ordinary Minkowski two-point functions in Rindler coordinates. A striking fact appeared immediately: when both points lie in the right wedge, the Minkowski Wightman function has an analytic continuation whose boundary values are related by an imaginary shift of Rindler time. This page explains why that KMS relation, rather than naive pointwise periodicity, is the defining structure of thermal equilibrium.

The route is deliberately indirect. We first review thermal correlators in ordinary quantum mechanics and field theory, because the Unruh effect is not a mysterious new kind of heat; it is the Kubo–Martin–Schwinger condition applied to the boost Hamiltonian. We then connect the KMS condition to detailed balance, Euclidean time circles, Matsubara frequencies, and fluctuation–dissipation. Finally we apply the same logic to a uniformly accelerated detector and obtain

TU=a2π.\boxed{T_U={a\over 2\pi}}.

Here aa is the detector’s proper acceleration, and kB==c=1k_B=\hbar=c=1. The result says that the Minkowski vacuum, restricted to the algebra of observables accessible to one uniformly accelerated observer, behaves as a thermal state with respect to that observer’s proper-time evolution.

Required background. Rindler coordinates and Green functions supplies boost time, wedge localization, and the Minkowski Wightman function in accelerated coordinates.

Helpful background. In–out, in–in, and Schwinger–Keldysh functionals distinguishes amplitudes from expectation values, while proper time, determinants, and thermal traces introduces Euclidean traces and periodic paths.

Thermal-time conventions. For an operator AA in the Heisenberg picture,

A(t)=eiHtAeiHt.A(t)=e^{iHt}Ae^{-iHt}.

A thermal expectation value at inverse temperature β\beta is

Oβ=1ZTr(eβHO),Z=TreβH.\langle \mathcal O\rangle_\beta={1\over Z}\operatorname{Tr}\left(e^{-\beta H}\mathcal O\right), \qquad Z=\operatorname{Tr}e^{-\beta H}.

For Rindler time we write

η=aτ,\eta=a\tau,

where τ\tau is proper time along the worldline ρ=1/a\rho=1/a. The inverse temperature conjugate to dimensionless boost time η\eta is βR=2π\beta_R=2\pi. The inverse temperature conjugate to proper time τ\tau is therefore

β=2πa.\beta={2\pi\over a}.

Let

CAB>(t)=A(t)B(0)β,CAB<(t)=B(0)A(t)β.C^>_{AB}(t)=\langle A(t)B(0)\rangle_\beta, \qquad C^<_{AB}(t)=\langle B(0)A(t)\rangle_\beta.

The two functions are not generally equal. Their difference is a commutator and controls response; their sum is noise. Thermal equilibrium relates them by a complex-time shift.

Using

A(tiβ)=eβHA(t)eβH,A(t-i\beta)=e^{\beta H}A(t)e^{-\beta H},

we find

CAB>(tiβ)=1ZTr(eβHeβHA(t)eβHB(0)).C^>_{AB}(t-i\beta) ={1\over Z}\operatorname{Tr}\left(e^{-\beta H}e^{\beta H}A(t)e^{-\beta H}B(0)\right).

Cancel the first two exponential factors and then use cyclicity of the trace:

CAB>(tiβ)=1ZTr(A(t)eβHB(0))=1ZTr(eβHB(0)A(t)).C^>_{AB}(t-i\beta) ={1\over Z}\operatorname{Tr}\left(A(t)e^{-\beta H}B(0)\right) ={1\over Z}\operatorname{Tr}\left(e^{-\beta H}B(0)A(t)\right).

Therefore

CAB>(tiβ)=CAB<(t).\boxed{ C^>_{AB}(t-i\beta)=C^<_{AB}(t). }

Equivalently,

A(t)B(0)β=B(0)A(t+iβ)β.\boxed{ \langle A(t)B(0)\rangle_\beta = \langle B(0)A(t+i\beta)\rangle_\beta. }

This is the KMS condition. It is the analytic definition of thermal equilibrium. For bosonic operators it often looks like periodicity in imaginary time, but the precise statement is stronger and more careful: shifting by iβ-i\beta changes the operator ordering.

The trace manipulation is literal for a finite system with a trace-class Gibbs operator. In continuum QFT, the thermodynamic limit may not admit a density matrix of the form Z1eβHZ^{-1}e^{-\beta H} on the vacuum Hilbert space. The intrinsic statement is then that the correlators are boundary values of functions analytic in the KMS strip and obey the same boundary relation. This is the form needed for wedge-local quantum fields.

KMS strip and imaginary-time thermal circle

The KMS condition says that a thermal Wightman function has analytic continuation in a strip of height β\beta. The lower edge is the oppositely ordered correlator. In Euclidean language, the same information becomes a thermal circle of circumference β\beta.

For a single bosonic operator JJ, define

G>(t)=J(t)J(0)β,G<(t)=J(0)J(t)β.G^>(t)=\langle J(t)J(0)\rangle_\beta, \qquad G^<(t)=\langle J(0)J(t)\rangle_\beta.

Then KMS gives

G>(tiβ)=G<(t).G^>(t-i\beta)=G^<(t).

For a time-translation-invariant thermal state, G<(t)=G>(t)G^<(t)=G^>(-t) if JJ is Hermitian. Thus a common one-operator form is

G>(tiβ)=G>(t).G^>(t-i\beta)=G^>(-t).

This is exactly the form that appears for the Minkowski vacuum restricted to a Rindler wedge.

The KMS condition becomes especially physical after a Fourier transform. Take a Hermitian operator JJ and insert a complete set of energy eigenstates,

Hn=Enn.H|n\rangle=E_n|n\rangle.

With the convention

G>(ω)=dteiωtG>(t),G^>(\omega)=\int_{-\infty}^{\infty}dt\,e^{i\omega t}G^>(t),

we obtain

G>(ω)=2πZm,neβEmJmn2δ(ω(EnEm)).G^>(\omega) ={2\pi\over Z}\sum_{m,n}e^{-\beta E_m}|J_{mn}|^2 \delta\left(\omega-(E_n-E_m)\right).

The same calculation gives

G<(ω)=2πZm,neβEmJmn2δ(ω(EmEn)).G^<(\omega) ={2\pi\over Z}\sum_{m,n}e^{-\beta E_m}|J_{mn}|^2 \delta\left(\omega-(E_m-E_n)\right).

Relabeling mnm\leftrightarrow n in the second expression yields

G<(ω)=eβωG>(ω).G^<(\omega)=e^{-\beta\omega}G^>(\omega).

For a Hermitian operator, G<(ω)=G>(ω)G^<(\omega)=G^>(-\omega), so

G>(ω)=eβωG>(ω).\boxed{ G^>(-\omega)=e^{-\beta\omega}G^>(\omega). }

This is detailed balance. Processes related by reversing the direction of energy flow differ by the Boltzmann factor.

For a two-level detector with gap Ω>0\Omega>0, weak coupling and stationary long-time evolution give the Golden Rule rates

Γ=λ2μeg2F(Ω),Γ=λ2μeg2F(Ω),\Gamma_\uparrow=\lambda^2|\mu_{eg}|^2\mathcal F(\Omega), \qquad \Gamma_\downarrow=\lambda^2|\mu_{eg}|^2\mathcal F(-\Omega),

where

F(Ω)=dseiΩsG+(s).\mathcal F(\Omega)=\int_{-\infty}^{\infty}ds\,e^{-i\Omega s}G^+(s).

The KMS condition implies

ΓΓ=eβΩ.\boxed{ {\Gamma_\uparrow\over\Gamma_\downarrow}=e^{-\beta\Omega}. }

Equivalently, for a bosonic oscillator mode of frequency Ω\Omega,

Γn(Ω),Γ1+n(Ω),\Gamma_\uparrow\propto n(\Omega), \qquad \Gamma_\downarrow\propto 1+n(\Omega),

with

n(Ω)=1eβΩ1.\boxed{ n(\Omega)={1\over e^{\beta\Omega}-1}. }

The 11 in 1+n1+n is spontaneous emission; the nn term is stimulated emission. Equilibrium is the condition that upward and downward transitions balance after including the Boltzmann weights of the detector states.

Detailed balance for a two-level detector in a thermal bath

A two-level detector in a thermal bath has upward and downward transition rates related by detailed balance. The ratio Γ/Γ=eβΩ\Gamma_\uparrow/\Gamma_\downarrow=e^{-\beta\Omega} is the operational content of temperature.

This is the right language for acceleration. An accelerated detector does not infer temperature by reading a coordinate label. It infers temperature because its transition rates obey detailed balance with a definite β\beta.

Thermal traces can be written as Euclidean path integrals. For a bosonic coordinate qq, the matrix element of eβHe^{-\beta H} is an imaginary-time transition amplitude. Taking the trace identifies the endpoints, so

Z(β)=TreβH=q(β)=q(0)DqeSE[q].Z(\beta)=\operatorname{Tr}e^{-\beta H} =\int_{q(\beta)=q(0)}\mathcal Dq\,e^{-S_E[q]}.

For a scalar field,

Z(β)=ϕ(τE+β,x)=ϕ(τE,x)DϕeSE[ϕ].Z(\beta)=\int_{\phi(\tau_E+\beta,\mathbf x)=\phi(\tau_E, \mathbf x)}\mathcal D\phi\,e^{-S_E[\phi]}.

Fermions are antiperiodic rather than periodic,

ψ(τE+β,x)=ψ(τE,x),\psi(\tau_E+\beta, \mathbf x)=-\psi(\tau_E, \mathbf x),

because the trace over fermionic states and coherent-state variables produces an extra sign.

For bosons, the allowed Euclidean frequencies are

ωn=2πnβ,nZ.\boxed{ \omega_n={2\pi n\over\beta}, \qquad n\in\mathbb Z. }

For fermions,

ωn=2π(n+12)β,nZ.\boxed{ \omega_n={2\pi\left(n+{1\over2}\right)\over\beta}, \qquad n\in\mathbb Z. }

As the simplest example, consider a harmonic oscillator of frequency ω0\omega_0. Its Euclidean thermal Green function satisfies

(d2dτE2+ω02)GE(τE)=δβ(τE),\left(-{d^2\over d\tau_E^2}+\omega_0^2\right)G_E(\tau_E)=\delta_\beta(\tau_E),

where δβ\delta_\beta is the periodic delta function on the circle. Therefore

GE(τE)=1βnZeiωnτEωn2+ω02.\boxed{ G_E(\tau_E) ={1\over\beta}\sum_{n\in\mathbb Z} {e^{i\omega_n\tau_E}\over \omega_n^2+\omega_0^2}. }

For 0τEβ0\le \tau_E\le \beta, the same answer can be written as

GE(τE)=12ω0cosh[ω0(β2τEβ)]sinh(βω02),\boxed{ G_E(\tau_E) ={1\over2\omega_0} {\cosh\left[\omega_0\left({\beta\over2}-|\tau_E|_\beta\right)\right] \over \sinh\left({\beta\omega_0\over2}\right)}, }

where τEβ|\tau_E|_\beta means the shortest distance to 00 on the thermal circle. In the zero-temperature limit β\beta\to\infty, the sum becomes an integral and

GE(τE)12ω0eω0τE.G_E(\tau_E)\to {1\over2\omega_0}e^{-\omega_0|\tau_E|}.

Euclidean thermal circle and Matsubara modes

A thermal trace becomes a Euclidean path integral on a circle of circumference β\beta. Bosonic fields are periodic and therefore have Matsubara frequencies ωn=2πn/β\omega_n=2\pi n/\beta.

The connection with Rindler space is now almost visible by inspection. Euclidean Rindler time is an angle. If the geometry is smooth at the origin, that angle must have period 2π2\pi. Thermal periodicity is therefore built into the regularity of the Euclidean plane.

Thermal field theory distinguishes two real-time questions. One asks how much the system fluctuates in equilibrium. The other asks how the system responds to a small external perturbation.

For a Hermitian operator JJ, define the spectral function

ρ(ω)=G>(ω)G<(ω)=dteiωt[J(t),J(0)]β,\rho(\omega)=G^>(\omega)-G^<(\omega) =\int_{-\infty}^{\infty}dt\,e^{i\omega t} \langle[J(t),J(0)]\rangle_\beta,

and the symmetrized noise

S(ω)=12(G>(ω)+G<(ω))=12dteiωt{J(t),J(0)}β.S(\omega)={1\over2}\left(G^>(\omega)+G^<(\omega)\right) ={1\over2}\int_{-\infty}^{\infty}dt\,e^{i\omega t} \langle\{J(t),J(0)\}\rangle_\beta.

Using

G<(ω)=eβωG>(ω),G^<(\omega)=e^{-\beta\omega}G^>(\omega),

we get

ρ(ω)=(1eβω)G>(ω),\rho(\omega)=\left(1-e^{-\beta\omega}\right)G^>(\omega),

and

S(ω)=12(1+eβω)G>(ω).S(\omega)={1\over2}\left(1+e^{-\beta\omega}\right)G^>(\omega).

Therefore

S(ω)=12coth(βω2)ρ(ω).\boxed{ S(\omega)={1\over2}\coth\left({\beta\omega\over2}\right)\rho(\omega). }

This is the fluctuation–dissipation relation. In a convention where the retarded function is

GR(t)=iθ(t)[J(t),J(0)]β,G_R(t)=-i\theta(t)\langle[J(t),J(0)]\rangle_\beta,

the dissipative part is proportional to the spectral function,

ρ(ω)=2ImGR(ω)\rho(\omega)=-2\operatorname{Im}G_R(\omega)

for real ω\omega, with the overall sign tied to the definition of GRG_R and the Fourier transform. The convention-independent point is that the commutator measures response, while the anticommutator measures equilibrium noise, and KMS fixes their ratio.

This relation is why the same mathematics appears in apparently different problems: blackbody radiation, current fluctuations, hot matter in a box, detector clicks, and accelerated motion in the vacuum.

Uniform acceleration and the Unruh response

Section titled “Uniform acceleration and the Unruh response”

A uniformly accelerated detector with proper acceleration aa follows

T(τ)=1asinh(aτ),X(τ)=1acosh(aτ),T(\tau)={1\over a}\sinh(a\tau), \qquad X(\tau)={1\over a}\cosh(a\tau),

with fixed transverse coordinates. This is the Rindler trajectory ρ=1/a\rho=1/a, η=aτ\eta=a\tau.

For a massless scalar field in four-dimensional Minkowski space,

GM+(X,X)=14π21(TTi0)2+XX2.G_M^+(X,X')={1\over4\pi^2} {1\over-(T-T'-i0)^2+|\mathbf X-\mathbf X'|^2}.

Along the accelerated worldline, set s=ττs=\tau-\tau'. The invariant denominator becomes

(T(τ)T(τ)i0)2+(X(τ)X(τ))2=4a2sinh2(a(si0)2).-(T(\tau)-T(\tau')-i0)^2+(X(\tau)-X(\tau'))^2 =-{4\over a^2}\sinh^2\left({a(s-i0)\over2}\right).

Thus

G+(s)=a216π21sinh2(a(si0)2).\boxed{ G^+(s) =-{a^2\over16\pi^2} {1\over\sinh^2\left({a(s-i0)\over2}\right)}. }

This correlator satisfies the KMS relation

G+(siβ)=G+(s),β=2πa.G^+(s-i\beta)=G^+(-s), \qquad \beta={2\pi\over a}.

Therefore a detector following this trajectory must satisfy detailed balance at temperature

TU=1β=a2π.\boxed{ T_U={1\over\beta}={a\over2\pi}. }

For an eternally coupled detector, stationarity reduces the long-time transition rate to the response kernel

F(Ω)=dseiΩsG+(s).\mathcal F(\Omega)=\int_{-\infty}^{\infty}ds\,e^{-i\Omega s}G^+(s).

Evaluating the contour integral gives

F(Ω)=Ω2π1e2πΩ/a1.\boxed{ \mathcal F(\Omega)={\Omega\over2\pi} {1\over e^{2\pi\Omega/a}-1}. }

This compact formula is valid for either sign of Ω\Omega with the usual i0i0 prescription. Coupling constants and detector matrix elements have been stripped off. For Ω>0\Omega>0, it gives the stationary excitation rate of a detector initially in its ground state. For Ω<0\Omega<0, it gives de-excitation and includes spontaneous emission.

Finite-time detectors require a switching function χ(τ)\chi(\tau). Their probability involves a double integral,

Pgedτdτχ(τ)χ(τ)eiΩ(ττ)G+(τ,τ),P_{g\to e}\propto \int d\tau\,d\tau'\, \chi(\tau)\chi(\tau')e^{-i\Omega(\tau-\tau')} G^+(\tau,\tau'),

and contains switching transients. Exact Planckian stationarity is recovered in the appropriate long-interaction limit; a short measurement need not look exactly thermal.

For clarity, separate the two signs of the detector gap. When Ω>0\Omega>0, the detector is excited and the response is Planck suppressed:

F(Ω)=Ω2π1e2πΩ/a1.\mathcal F(\Omega)={\Omega\over2\pi}{1\over e^{2\pi\Omega/a}-1}.

For de-excitation, write Ω=Ω\Omega=-|\Omega|. The same compact formula gives

F(Ω)=Ω2π(1+1e2πΩ/a1),\mathcal F(-|\Omega|)={|\Omega|\over2\pi} \left(1+{1\over e^{2\pi |\Omega|/a}-1}\right),

where the first term is spontaneous emission and the second is stimulated emission by the thermal Rindler bath.

Uniformly accelerated detector and thermal response

An eternally and uniformly accelerated detector samples a stationary KMS Wightman function with imaginary proper-time separation 2π/a2\pi/a. Its long-time excitation rate is Planckian at TU=a/(2π)T_U=a/(2\pi); finite switching adds transients.

Nothing in this calculation says that the Minkowski vacuum is a thermal state for inertial observers. It says that the pair consisting of a state and a time evolution is thermal: the Minkowski vacuum, restricted to right-wedge observables, is thermal with respect to the boost generator.

Euclidean smoothness and the geometric temperature

Section titled “Euclidean smoothness and the geometric temperature”

There is a geometric way to see the same inverse temperature without calculating detector rates. Start from the right-wedge metric

ds2=ρ2dη2dρ2dy2.ds^2=\rho^2d\eta^2-d\rho^2-d\mathbf y^2.

Continue to Euclidean boost time by setting

η=iθ.\eta=-i\theta.

and define the positive Euclidean line element by dsE2=ds2ds_E^2=-ds^2. Then

dsE2=dρ2+ρ2dθ2+dy2.ds_E^2=d\rho^2+\rho^2d\theta^2+d\mathbf y^2.

The (ρ,θ)(\rho,\theta) part is just the flat Euclidean plane in polar coordinates. The point ρ=0\rho=0 is smooth only if

θθ+2π.\theta\sim\theta+2\pi.

If an observer uses proper time τ=η/a\tau=\eta/a, then Euclidean proper time is τE=θ/a\tau_E=\theta/a, and the period is

τEτE+2πa.\tau_E\sim \tau_E+{2\pi\over a}.

Thermal QFT identifies Euclidean time period with inverse temperature, so again

β=2πa,TU=a2π.\beta={2\pi\over a}, \qquad T_U={a\over2\pi}.

This argument is the flat-spacetime prototype of Hawking’s black-hole temperature calculation. Near any nonextremal horizon, the Euclidean metric looks like a plane in polar coordinates. Avoiding a conical singularity fixes the Euclidean period, and therefore fixes the temperature.

Wedge thermality and the regulated density-matrix picture

Section titled “Wedge thermality and the regulated density-matrix picture”

The detector calculation gives the operational meaning of the temperature. The operator-algebra statement is even sharper. Let KRK_R be the dimensionless generator of Rindler time translations in the right wedge. With a UV regulator that factorizes left- and right-wedge degrees of freedom, one may write

ρR(reg)=TrL0M0M=1ZRe2πKR.\boxed{ \rho_R^{\rm(reg)} =\operatorname{Tr}_L |0_M\rangle\langle0_M| ={1\over Z_R}e^{-2\pi K_R}. }

Here TrL\operatorname{Tr}_L traces over degrees of freedom in the left wedge. If the physical Hamiltonian for proper time is

Hτ=aKR,H_\tau=aK_R,

then

ρR(reg)=1ZReβHτ,β=2πa.\rho_R^{\rm(reg)}={1\over Z_R}e^{-\beta H_\tau}, \qquad \beta={2\pi\over a}.

Mode by mode, the same regulated statement has the schematic thermofield-double form

0Mω,k(1e2πω)1/2n=0eπωnnωkRnωkL.|0_M\rangle \sim \prod_{\omega,\mathbf k} \left(1-e^{-2\pi\omega}\right)^{1/2} \sum_{n=0}^{\infty}e^{-\pi\omega n} |n_{\omega\mathbf k}\rangle_R|n_{\omega\mathbf k}\rangle_L.

Tracing over the left wedge gives

ρR,reg(ω,k)=(1e2πω)n=0e2πωnnωkRRnωk,\rho_{R,\rm reg}^{(\omega,\mathbf k)} =\left(1-e^{-2\pi\omega}\right) \sum_{n=0}^{\infty}e^{-2\pi\omega n} |n_{\omega\mathbf k}\rangle_R\,{}_R\langle n_{\omega\mathbf k}|,

and hence

Nωk=1e2πω1.\langle N_{\omega\mathbf k}\rangle ={1\over e^{2\pi\omega}-1}.

Minkowski vacuum as an entangled state of left and right Rindler modes

The global Minkowski vacuum is pure but entangled across the Rindler horizons. A modewise or UV-regulated factorization gives ρR(reg)e2πKR\rho_R^{\rm(reg)}\propto e^{-2\pi K_R}. Intrinsically, the continuum wedge state is KMS for boost flow rather than a trace-class density matrix.

The exact construction of the modes requires a careful choice of Unruh modes, but the thermal weights are fixed by analyticity and wedge localization. In continuum algebraic QFT, a wedge algebra is not generally represented by a tensor factor with a trace-class reduced density matrix; the literal partial trace above is a regulated mnemonic. The regulator-independent content is the KMS property of the vacuum under modular boost flow, as established by the Bisognano–Wichmann theorem under its standard axioms.

Thermal equilibrium can be characterized without mentioning particles: it is the KMS analyticity condition

A(t)B(0)β=B(0)A(t+iβ)β.\langle A(t)B(0)\rangle_\beta = \langle B(0)A(t+i\beta)\rangle_\beta.

In Fourier space, KMS becomes detailed balance. For a detector with level spacing Ω\Omega,

ΓΓ=eβΩ.{\Gamma_\uparrow\over\Gamma_\downarrow}=e^{-\beta\Omega}.

In Euclidean field theory, the same condition becomes a compact imaginary-time circle of circumference β\beta, with bosonic Matsubara frequencies ωn=2πn/β\omega_n=2\pi n/\beta.

For a uniformly accelerated observer, Rindler time is boost time. The Minkowski vacuum restricted to a single wedge is KMS with inverse temperature 2π2\pi in dimensionless boost time. Along a worldline of proper acceleration aa, this becomes

β=2πa,TU=a2π.\beta={2\pi\over a}, \qquad T_U={a\over2\pi}.

The Unruh effect is therefore not a statement that the inertial vacuum is globally mixed. It is a statement that the right-wedge restriction of the inertial vacuum is KMS with respect to boost evolution. The Planck formula is a stationary long-time detector rate; finite switching adds protocol-dependent transients. A reduced density matrix is a useful regulated or modewise representation, while the continuum theorem is intrinsically a statement about the wedge algebra and modular flow.

Replacing KMS by real-time periodicity. The shift is in imaginary time and, for real-time Wightman functions, changes operator ordering. The analytic strip and its boundary values are part of the statement.

Using the Feynman propagator to count detector clicks. Transition rates are controlled by the Wightman function evaluated along the detector trajectory, with the detector’s switching and i0i0 prescription specified.

Calling the Minkowski vacuum a global thermal bath. The vacuum remains pure and Poincaré invariant. Thermality appears after restricting to one wedge and using boost evolution.

Forgetting which time generates the Hamiltonian. The inverse temperature is 2π2\pi for dimensionless boost time η\eta. Along the trajectory η=aτ\eta=a\tau, the proper inverse temperature is 2π/a2\pi/a.

Treating a finite-time response as exactly Planckian. The compact response formula is a stationary long-time rate for uniform acceleration. A finite switching function produces transients that depend on the measurement protocol.

Confusing the Rindler and Minkowski vacua. The Rindler vacuum is empty for KRK_R modes but singular at the horizons. The Minkowski vacuum is regular there and KMS when restricted to a wedge.

Taking the continuum partial trace literally. A UV regulator can supply a left–right tensor factor and a reduced density matrix. The intrinsic continuum statement is the KMS/modular-flow relation for the wedge algebra.

Exercise 1: Derive KMS from trace cyclicity

Section titled “Exercise 1: Derive KMS from trace cyclicity”

Let

CAB>(t)=Z1Tr(eβHA(t)B(0)).C^>_{AB}(t)=Z^{-1}\operatorname{Tr}\left(e^{-\beta H}A(t)B(0)\right).

Using A(tiβ)=eβHA(t)eβHA(t-i\beta)=e^{\beta H}A(t)e^{-\beta H}, prove

CAB>(tiβ)=CAB<(t),CAB<(t)=Z1Tr(eβHB(0)A(t)).C^>_{AB}(t-i\beta)=C^<_{AB}(t), \qquad C^<_{AB}(t)=Z^{-1}\operatorname{Tr}\left(e^{-\beta H}B(0)A(t)\right).
Solution

Start from

CAB>(tiβ)=1ZTr(eβHA(tiβ)B(0)).C^>_{AB}(t-i\beta) ={1\over Z}\operatorname{Tr}\left(e^{-\beta H}A(t-i\beta)B(0)\right).

Substitute

A(tiβ)=eβHA(t)eβH.A(t-i\beta)=e^{\beta H}A(t)e^{-\beta H}.

Then

CAB>(tiβ)=1ZTr(eβHeβHA(t)eβHB(0))=1ZTr(A(t)eβHB(0)).C^>_{AB}(t-i\beta) ={1\over Z}\operatorname{Tr}\left(e^{-\beta H}e^{\beta H}A(t)e^{-\beta H}B(0)\right) ={1\over Z}\operatorname{Tr}\left(A(t)e^{-\beta H}B(0)\right).

By cyclicity of the trace,

Tr(A(t)eβHB(0))=Tr(eβHB(0)A(t)).\operatorname{Tr}\left(A(t)e^{-\beta H}B(0)\right) =\operatorname{Tr}\left(e^{-\beta H}B(0)A(t)\right).

Therefore

CAB>(tiβ)=CAB<(t).C^>_{AB}(t-i\beta)=C^<_{AB}(t).

Exercise 2: Detailed balance in frequency space

Section titled “Exercise 2: Detailed balance in frequency space”

For a Hermitian operator JJ, define

G>(t)=J(t)J(0)β,G>(ω)=dteiωtG>(t).G^>(t)=\langle J(t)J(0)\rangle_\beta, \qquad G^>(\omega)=\int dt\,e^{i\omega t}G^>(t).

Show that

G>(ω)=eβωG>(ω).G^>(-\omega)=e^{-\beta\omega}G^>(\omega).
Solution

Insert energy eigenstates:

G>(t)=1Zm,neβEmei(EmEn)tJmn2.G^>(t)={1\over Z}\sum_{m,n}e^{-\beta E_m}e^{i(E_m-E_n)t}|J_{mn}|^2.

Fourier transforming gives

G>(ω)=2πZm,neβEmJmn2δ(ω(EnEm)).G^>(\omega)={2\pi\over Z}\sum_{m,n}e^{-\beta E_m}|J_{mn}|^2 \delta\left(\omega-(E_n-E_m)\right).

Now compute G>(ω)G^>(-\omega):

G>(ω)=2πZm,neβEmJmn2δ(ω(EnEm)).G^>(-\omega)={2\pi\over Z}\sum_{m,n}e^{-\beta E_m}|J_{mn}|^2 \delta\left(-\omega-(E_n-E_m)\right).

Relabel mnm\leftrightarrow n:

G>(ω)=2πZm,neβEnJmn2δ(ω(EmEn)).G^>(-\omega)={2\pi\over Z}\sum_{m,n}e^{-\beta E_n}|J_{mn}|^2 \delta\left(-\omega-(E_m-E_n)\right).

The delta function is now

δ(ω(EnEm)).\delta\left(\omega-(E_n-E_m)\right).

On its support,

En=Em+ω,E_n=E_m+\omega,

so

eβEn=eβωeβEm.e^{-\beta E_n}=e^{-\beta\omega}e^{-\beta E_m}.

Therefore

G>(ω)=eβωG>(ω).G^>(-\omega)=e^{-\beta\omega}G^>(\omega).

Exercise 3: Matsubara frequencies from periodicity

Section titled “Exercise 3: Matsubara frequencies from periodicity”

Let GE(τ)G_E(\tau) be a bosonic Euclidean thermal correlator satisfying

GE(τ+β)=GE(τ).G_E(\tau+\beta)=G_E(\tau).

Show that its Fourier series contains only frequencies ωn=2πn/β\omega_n=2\pi n/\beta. Then solve

(d2dτ2+ω02)GE(τ)=δβ(τ)\left(-{d^2\over d\tau^2}+\omega_0^2\right)G_E(\tau)=\delta_\beta(\tau)

in frequency space.

Solution

A periodic function on a circle of circumference β\beta has Fourier expansion

GE(τ)=nZcneiωnτ.G_E(\tau)=\sum_{n\in\mathbb Z}c_n e^{i\omega_n\tau}.

The condition GE(τ+β)=GE(τ)G_E(\tau+\beta)=G_E(\tau) requires

eiωnβ=1,e^{i\omega_n\beta}=1,

so

ωn=2πnβ.\omega_n={2\pi n\over\beta}.

Normalize the expansion as

GE(τ)=1βnZGneiωnτ,δβ(τ)=1βnZeiωnτ.G_E(\tau)={1\over\beta}\sum_{n\in\mathbb Z}G_n e^{i\omega_n\tau}, \qquad \delta_\beta(\tau)={1\over\beta}\sum_{n\in\mathbb Z}e^{i\omega_n\tau}.

Substitution into

(d2dτ2+ω02)GE(τ)=δβ(τ)\left(-{d^2\over d\tau^2}+\omega_0^2\right)G_E(\tau)=\delta_\beta(\tau)

gives

(ωn2+ω02)Gn=1.(\omega_n^2+\omega_0^2)G_n=1.

Hence

Gn=1ωn2+ω02,G_n={1\over\omega_n^2+\omega_0^2},

and

GE(τ)=1βnZeiωnτωn2+ω02.G_E(\tau)={1\over\beta}\sum_{n\in\mathbb Z} {e^{i\omega_n\tau}\over\omega_n^2+\omega_0^2}.

Starting from

ds2=ρ2dη2dρ2,ds^2=\rho^2d\eta^2-d\rho^2,

set η=iθ\eta=-i\theta and show that smoothness at ρ=0\rho=0 requires θθ+2π\theta\sim\theta+2\pi. Then derive the Unruh temperature for an observer whose proper time is τ=η/a\tau=\eta/a.

Solution

The continuation η=iθ\eta=-i\theta, together with dsE2=ds2ds_E^2=-ds^2 for the mostly-minus convention, gives

dsE2=dρ2+ρ2dθ2.ds_E^2=d\rho^2+\rho^2d\theta^2.

This is the Euclidean plane in polar coordinates. The point ρ=0\rho=0 is smooth only if the angular coordinate has its standard period,

θθ+2π.\theta\sim\theta+2\pi.

If η=aτ\eta=a\tau, then after continuation θ=aτE\theta=a\tau_E. The period of Euclidean proper time is therefore

aτEaτE+2π,a\tau_E\sim a\tau_E+2\pi,

or

τEτE+2πa.\tau_E\sim\tau_E+{2\pi\over a}.

Thermal field theory identifies the Euclidean time period with inverse temperature:

β=2πa.\beta={2\pi\over a}.

Thus

TU=1β=a2π.T_U={1\over\beta}={a\over2\pi}.

Exercise 5: Detailed balance for an accelerated detector

Section titled “Exercise 5: Detailed balance for an accelerated detector”

Assume the detector response function satisfies the KMS consequence

F(Ω)=eβΩF(Ω)\mathcal F(\Omega)=e^{-\beta\Omega}\mathcal F(-\Omega)

for Ω>0\Omega>0. Show that a two-level detector with energies Eg=0E_g=0 and Ee=ΩE_e=\Omega relaxes to thermal occupation probabilities.

Solution

Let pgp_g and pep_e be the probabilities for the ground and excited states. In equilibrium, the upward and downward probability currents balance:

pgΓ=peΓ.p_g\Gamma_\uparrow=p_e\Gamma_\downarrow.

The rates are proportional to the response functions,

ΓF(Ω),ΓF(Ω).\Gamma_\uparrow\propto\mathcal F(\Omega), \qquad \Gamma_\downarrow\propto\mathcal F(-\Omega).

Using detailed balance,

ΓΓ=F(Ω)F(Ω)=eβΩ.{\Gamma_\uparrow\over\Gamma_\downarrow} ={\mathcal F(\Omega)\over\mathcal F(-\Omega)} =e^{-\beta\Omega}.

Thus

pepg=ΓΓ=eβΩ.{p_e\over p_g}={\Gamma_\uparrow\over\Gamma_\downarrow}=e^{-\beta\Omega}.

This is exactly the thermal Boltzmann ratio for a two-level system:

pg=11+eβΩ,pe=eβΩ1+eβΩ.p_g={1\over1+e^{-\beta\Omega}}, \qquad p_e={e^{-\beta\Omega}\over1+e^{-\beta\Omega}}.

For a uniformly accelerated detector in the Minkowski vacuum, β=2π/a\beta=2\pi/a.

Exercise 6: De-excitation from the stationary response

Section titled “Exercise 6: De-excitation from the stationary response”

The stationary accelerated-detector response rate for a massless scalar in four dimensions can be written as

F(Ω)=Ω2π1e2πΩ/a1.\mathcal F(\Omega)={\Omega\over2\pi}{1\over e^{2\pi\Omega/a}-1}.

Show that for Ω=E<0\Omega=-E<0, with E>0E>0, this becomes

F(E)=E2π(1+1e2πE/a1).\mathcal F(-E)={E\over2\pi}\left(1+{1\over e^{2\pi E/a}-1}\right).
Solution

Substitute Ω=E\Omega=-E:

F(E)=E2π1e2πE/a1.\mathcal F(-E)={-E\over2\pi}{1\over e^{-2\pi E/a}-1}.

Multiply numerator and denominator by e2πE/ae^{2\pi E/a}:

F(E)=E2πe2πE/ae2πE/a1.\mathcal F(-E)={E\over2\pi}{e^{2\pi E/a}\over e^{2\pi E/a}-1}.

Since

exex1=1+1ex1,{e^x\over e^x-1}=1+{1\over e^x-1},

we get

F(E)=E2π(1+1e2πE/a1).\mathcal F(-E)={E\over2\pi}\left(1+{1\over e^{2\pi E/a}-1}\right).

The 11 is spontaneous emission; the Bose–Einstein term is stimulated emission.

  • J. J. Bisognano and E. H. Wichmann, “On the Duality Condition for a Hermitian Scalar Field,” Journal of Mathematical Physics 16 (1975), 985–1007, doi:10.1063/1.522605.
  • J. J. Bisognano and E. H. Wichmann, “On the Duality Condition for Quantum Fields,” Journal of Mathematical Physics 17 (1976), 303–321, doi:10.1063/1.522898.
  • P. C. W. Davies, “Scalar Particle Production in Schwarzschild and Rindler Metrics,” Journal of Physics A: Mathematical and General 8 (1975), 609–616, doi:10.1088/0305-4470/8/4/022.
  • S. A. Fulling, “Nonuniqueness of Canonical Field Quantization in Riemannian Space-Time,” Physical Review D 7 (1973), 2850–2862, doi:10.1103/PhysRevD.7.2850.
  • R. Kubo, “Statistical-Mechanical Theory of Irreversible Processes. I. General Theory and Simple Applications to Magnetic and Conduction Problems,” Journal of the Physical Society of Japan 12 (1957), 570–586, doi:10.1143/JPSJ.12.570.
  • P. C. Martin and J. Schwinger, “Theory of Many-Particle Systems. I,” Physical Review 115 (1959), 1342–1373, doi:10.1103/PhysRev.115.1342.
  • W. G. Unruh, “Notes on Black-Hole Evaporation,” Physical Review D 14 (1976), 870–892, doi:10.1103/PhysRevD.14.870.
  • N. D. Birrell and P. C. W. Davies, Quantum Fields in Curved Space, Cambridge Monographs on Mathematical Physics, Cambridge University Press, Cambridge, 1982.
  • R. Haag, Local Quantum Physics: Fields, Particles, Algebras, 2nd ed., Springer, Berlin, 1996.
  • M. Le Bellac, Thermal Field Theory, Cambridge Monographs on Mathematical Physics, Cambridge University Press, Cambridge, 1996.
  • S. Takagi, “Vacuum Noise and Stress Induced by Uniform Acceleration: Hawking–Unruh Effect in Rindler Manifold of Arbitrary Dimension,” Progress of Theoretical Physics Supplement 88 (1986), 1–142, doi:10.1143/PTPS.88.1.
  • R. M. Wald, Quantum Field Theory in Curved Spacetime and Black Hole Thermodynamics, Chicago Lectures in Physics, University of Chicago Press, Chicago, 1994.