Skip to content

Modular KMS Relations and Modular Correlators

Every faithful normal state is an equilibrium state for its own modular dynamics: two ordered correlators are the boundary values of one analytic function across a unit strip. The statement is intrinsic to the algebra–state pair. Its sign, strip, and operator order become unambiguous only after the modular-flow convention has been fixed.

Required background. Thermal density operators and the KMS condition supply the Gibbs calculation and physical-time KMS condition used for comparison; Tomita–Takesaki flow supplies the standard pair and modular automorphism group used in the theorem.

Helpful background. Relative modular operators and cocycles explain how the same analytic methods compare two states rather than one state with itself.

The chapter’s route from standard pairs to geometric flow locates the intrinsic KMS theorem. The representative-flow table keeps algebraic, wedge, ball, and regulated formulas separate, while the claim-licensing checklist tests any proposed physical interpretation.

Use the site convention

σs(A)=ΔisAΔ−is,K=−log⁡Δ,\sigma_s(A)=\Delta^{is}A\Delta^{-is}, \qquad K=-\log\Delta,

for the modular group of a standard pair (A,Ω)(\mathcal A,\Omega), with ω(A)=⟨Ω,AΩ⟩\omega(A)=\langle\Omega,A\Omega\rangle. For every bounded A,B∈AA,B\in\mathcal A, there is a bounded function FA,BF_{A,B}, continuous on the closed strip and analytic inside it, whose boundary values are

FA,B(s)=ω ⁣(Aσs(B)),FA,B(s−i)=ω ⁣(σs(B)A),s∈R,\begin{aligned} F_{A,B}(s)&=\omega\!\left(A\sigma_s(B)\right),\\ F_{A,B}(s-i)&=\omega\!\left(\sigma_s(B)A\right), \qquad s\in\mathbb R, \end{aligned}

with −1<Im⁡z<0-1<\operatorname{Im}z<0. Entire analytic elements are needed only if one wants to write the interior value literally as ω(Aσz(B))\omega(A\sigma_z(B)); the theorem for arbitrary bounded elements asserts the existence of the strip function. This is the modular KMS boundary theorem in Takesaki 1970, “The Modular Automorphism Group and the Kubo–Martin–Schwinger Boundary Condition,” pp. 63–72.

The same statement appears in three common forms:

evolutioncorrelator on the real edgeanalytic stripopposite edge
site flow σs\sigma_sF(s)=ω(Aσs(B))F(s)=\omega(A\sigma_s(B))−1<Im⁡z<0-1<\operatorname{Im}z<0F(s−i)=ω(σs(B)A)F(s-i)=\omega(\sigma_s(B)A)
reversed flow αs=σ−s\alpha_s=\sigma_{-s}G(s)=ω(Aαs(B))G(s)=\omega(A\alpha_s(B))0<Im⁡z<10<\operatorname{Im}z<1G(s+i)=ω(αs(B)A)G(s+i)=\omega(\alpha_s(B)A)
Gibbs Heisenberg flow τt\tau_tH(t)=ω(Aτt(B))H(t)=\omega(A\tau_t(B))0<Im⁡z<β0<\operatorname{Im}z<\betaH(t+iβ)=ω(τt(B)A)H(t+i\beta)=\omega(\tau_t(B)A)

For ρ=Z−1e−βH\rho=Z^{-1}e^{-\beta H} and τt(A)=eitHAe−itH\tau_t(A)=e^{itH}Ae^{-itH},

σs(A)=ρisAρ−is=τ−βs(A).\sigma_s(A)=\rho^{is}A\rho^{-is}=\tau_{-\beta s}(A).

Thus the lower modular strip maps to the upper physical-time strip because physical time is reversed. Modular ss is dimensionless; a temperature in physical units exists only after a particular physical evolution has been identified. The original infinite-system equilibrium formulation of the KMS boundary condition is developed in Haag, Hugenholtz, and Winnink 1967, §§1–3, pp. 215–229.

Let ρ=∑npn∣n⟩⟨n∣\rho=\sum_n p_n\lvert n\rangle\langle n\rvert with every pn>0p_n>0. Direct matrix multiplication gives

FA,B(z)=∑m,npm 1−izpn izAmnBnm.F_{A,B}(z)= \sum_{m,n}p_m^{\,1-iz}p_n^{\,iz}A_{mn}B_{nm}.

At z=s−iz=s-i, the coefficient becomes pn(pn/pm)isp_n(p_n/p_m)^{is}; relabeling mm and nn yields the required reversed order.

For a numerical convention check, choose

ρ=(3/4001/4),A=∣0⟩⟨1∣,B=∣1⟩⟨0∣.\rho=\begin{pmatrix}3/4&0\\0&1/4\end{pmatrix}, \qquad A=\lvert0\rangle\langle1\rvert, \qquad B=\lvert1\rangle\langle0\rvert.

Then

F(z)=34 3−iz,F(s−i)=14 3−is=Tr⁡ ⁣(ρ σs(B)A).F(z)=\frac34\,3^{-iz}, \qquad F(s-i)=\frac14\,3^{-is} =\operatorname{Tr}\!\left(\rho\,\sigma_s(B)A\right).

The equality is exact, so a reproducible implementation can use

RKMS(s)=∣F(s−i)−Tr⁡(ρ σs(B)A)∣max⁡(1,∣F(s−i)∣)R_{\mathrm{KMS}}(s)= \frac{\left\lvert F(s-i)- \operatorname{Tr}(\rho\,\sigma_s(B)A)\right\rvert} {\max(1,\lvert F(s-i)\rvert)}

with analytic target RKMS=0R_{\mathrm{KMS}}=0. Report the eigenvalue floor min⁡pn\min p_n, the sampled ss interval, arithmetic precision, and the maximum residual. Faithfulness is not cosmetic: if a pnp_n vanishes, the complex powers and the modular group must first be defined on the support corner.

The wedge gives a QFT realization in which modular time has a geometric meaning. For the massless scalar field in four dimensions, use point fields only as shorthand for a limit of spacetime-smeared fields supported in the right wedge. Along a uniformly accelerated orbit x(τ)x(\tau) of acceleration aa, write ϕ(τ)=ϕ(x(τ))\phi(\tau)=\phi(x(\tau)) and set

βU=2πa,G(z)=−a216π2sinh⁡2(az/2).\beta_{\mathrm U}=\frac{2\pi}{a}, \qquad \mathcal G(z)=-\frac{a^2}{16\pi^2\sinh^2(az/2)}.

The function is analytic for −βU<Im⁡z<0-\beta_{\mathrm U}<\operatorname{Im}z<0. With the site convention, the Bisognano–Wichmann trajectory moves the second field by proper time −βUs-\beta_{\mathrm U}s; boost invariance then gives

F(s)=⟨Ω,ϕ(0)σs(ϕ(0))Ω⟩=⟨Ω,ϕ(βUs)ϕ(0)Ω⟩=lim⁡ϵ↓0G(βUs−iϵ).\begin{aligned} F(s) &=\langle\Omega,\phi(0)\sigma_s(\phi(0))\Omega\rangle\\ &=\langle\Omega,\phi(\beta_{\mathrm U}s)\phi(0)\Omega\rangle =\lim_{\epsilon\downarrow0}\mathcal G(\beta_{\mathrm U}s-i\epsilon). \end{aligned}

Approaching the other edge from inside the strip,

lim⁡ϵ↓0G ⁣(βUs−iβU+iϵ)=⟨Ω,σs(ϕ(0))ϕ(0)Ω⟩.\lim_{\epsilon\downarrow0} \mathcal G\!\left( \beta_{\mathrm U}s-i\beta_{\mathrm U}+i\epsilon \right) =\langle\Omega,\sigma_s(\phi(0))\phi(0)\Omega\rangle.

The equality follows from sinh⁡(w−iπ)=−sinh⁡w\sinh(w-i\pi)=-\sinh w; the change from −i0-i0 to +i0+i0 carries the operator-order reversal. It verifies the unit modular strip and the factor 2π2\pi simultaneously. The same periodicity underlies the accelerated-detector relation in Unruh 1976, §III, especially Eq. (3.10), pp. 883–885.

For a reproducible boundary test, choose a>0a>0, a dimensionless regulator aϵa\epsilon, and a coincidence exclusion a∣u∣≥aumin⁡>0a\lvert u\rvert\ge a u_{\min}>0. Normalize the dimension-two correlator as G^=G/a2\widehat{\mathcal G}=\mathcal G/a^2 and compare

Rϵ(u)=∣G^(u−iβU+iϵ)−G^(u+iϵ)∣max⁡(1,∣G^(u+iϵ)∣).R_{\epsilon}(u)= \frac{\left\lvert \widehat{\mathcal G}(u-i\beta_{\mathrm U}+i\epsilon)- \widehat{\mathcal G}(u+i\epsilon) \right\rvert} {\max(1,\lvert\widehat{\mathcal G}(u+i\epsilon)\rvert)}.

Analytically Rϵ(u)=0R_\epsilon(u)=0; numerical uncertainty is therefore only roundoff and cancellation near a pole. The smearing widths, support margin from the wedge boundary, aϵa\epsilon, and aumin⁡a u_{\min} are part of the result, not optional implementation details.

Let EK(dk)E_K(dk) be the projection-valued measure of KK. Since KΩ=0K\Omega=0,

FA,B(s)=∫Re−isk dνA,B(k),νA,B(E)=⟨A∗Ω,EK(E)BΩ⟩.F_{A,B}(s) =\int_{\mathbb R}e^{-isk}\,d\nu_{A,B}(k), \qquad \nu_{A,B}(E)= \langle A^*\Omega,E_K(E)B\Omega\rangle.

This is generally a complex cross measure, not a probability measure. For analytic elements, the lower boundary inserts e−ke^{-k}, while modular invariance gives

FA,B(s−i)=FB,A(−s).F_{A,B}(s-i)=F_{B,A}(-s).

Fourier uniqueness therefore yields the weak measure identity

e−k dνA,B(k)=dνB,A(−k).\boxed{ e^{-k}\,d\nu_{A,B}(k)=d\nu_{B,A}(-k) }.

The minus sign in the reflected argument and the factor e−ke^{-k} both follow from the chosen transform e−iske^{-isk} and the lower strip. Changing either convention changes the displayed formula. For nonsmooth boundary data, interpret it after smearing in kk rather than as a pointwise identity between spectral densities.

The modular theorem concerns bounded algebra elements. A useful entire-analytic approximation is

Aϵ=1πϵ∫−∞∞dt  e−t2/ϵσt(A),ϵ>0.A_\epsilon= \frac{1}{\sqrt{\pi\epsilon}} \int_{-\infty}^{\infty}dt\; e^{-t^2/\epsilon}\sigma_t(A), \qquad \epsilon>0.

The integral is strong-∗* convergent, ∥Aϵ∥≤∥A∥\lVert A_\epsilon\rVert\le\lVert A\rVert, AϵA_\epsilon is entire analytic for the modular group, and Aϵ→AA_\epsilon\to A strongly-∗* as ϵ↓0\epsilon\downarrow0. This mollification controls modular time; it does not replace spacetime smearing of a local field. Field correlators additionally require test functions, a common invariant domain, and distributional boundary values.

The two-level example detects a missing KMS shift immediately. At s=0s=0,

F(0)=ω(AB)=34,ω(BA)=14.F(0)=\omega(AB)=\frac34, \qquad \omega(BA)=\frac14.

If one reverses the order without moving to z=−iz=-i, the absolute residual is 1/21/2, even though the correctly shifted residual is zero. What survives is ordinary stationarity; the cyclic boundary identity has been discarded.

The wedge exposes a different failure. The analytic function has poles at z=2πin/az=2\pi i n/a and behaves near coincidence as −1/(4π2z2)-1/(4\pi^2z^2). A path that crosses z=0z=0 is not an analytic continuation inside the KMS strip. Approaching the pole from the wrong side replaces (u−i0)−2(u-i0)^{-2} by (u+i0)−2(u+i0)^{-2} and hence changes the commutator/contact distribution. A pointwise check at u≠0u\ne0 can miss this failure, which is why the test must retain spacetime smearing and vary umin⁡u_{\min} and ϵ\epsilon together.

Giving a strip width without a flow sign. A unit width does not determine whether the strip is above or below the real axis. Display the adjoint action, the correlator order, and the shifted boundary together.

Calling modular temperature a thermometer reading. Modular time is dimensionless. The Unruh and Gibbs temperatures arise only after a boost or Hamiltonian evolution has been identified with the modular flow.

Treating point fields as bounded analytic elements. A Wightman function is a distributional boundary value. State the spacetime smearing, operator domain, and pole prescription before using the bounded-element theorem.

1. A two-level strip and a failed shortcut

Section titled “1. A two-level strip and a failed shortcut”

For the two-level data above, derive F(z)F(z) from matrix multiplication. Verify the lower boundary for every real ss, and evaluate the residual produced by replacing F(s−i)F(s-i) with F(s)F(s) at s=0s=0.

Solution

Only A01B10A_{01}B_{10} is nonzero. The general finite-dimensional formula therefore gives

F(z)=p01−izp1iz=34(13)iz=34 3−iz.F(z)=p_0^{1-iz}p_1^{iz} =\frac34\left(\frac13\right)^{iz} =\frac34\,3^{-iz}.

At z=s−iz=s-i,

F(s−i)=34 3−is−1=14 3−is.F(s-i)=\frac34\,3^{-is-1} =\frac14\,3^{-is}.

Moreover σs(B)=3−isB\sigma_s(B)=3^{-is}B and BA=∣1⟩⟨1∣BA=\lvert1\rangle\langle1\rvert, so

Tr⁡(ρ σs(B)A)=3−isTr⁡(ρBA)=14 3−is.\operatorname{Tr}(\rho\,\sigma_s(B)A) =3^{-is}\operatorname{Tr}(\rho BA) =\frac14\,3^{-is}.

At s=0s=0, the unshifted same-order value is F(0)=3/4F(0)=3/4, whereas the reversed-order value is 1/41/4. The failed shortcut has absolute residual 1/21/2.

2. Derive the modular detailed-balance sign

Section titled “2. Derive the modular detailed-balance sign”

Starting from the lower-strip boundary relation and the transform FA,B(s)=∫e−iskdνA,B(k)F_{A,B}(s)=\int e^{-isk}d\nu_{A,B}(k), derive the relation between νA,B\nu_{A,B} and νB,A\nu_{B,A}. Then state how the result changes if the Fourier kernel is e+iske^{+isk}.

Solution

At the lower edge,

FA,B(s−i)=∫e−iske−k dνA,B(k).F_{A,B}(s-i) =\int e^{-isk}e^{-k}\,d\nu_{A,B}(k).

Modular invariance gives

ω(σs(B)A)=ω(Bσ−s(A))=FB,A(−s)=∫eisk dνB,A(k).\omega(\sigma_s(B)A) =\omega(B\sigma_{-s}(A)) =F_{B,A}(-s) =\int e^{isk}\,d\nu_{B,A}(k).

Changing k↦−kk\mapsto-k in the last integral and using uniqueness of Fourier transforms of finite complex measures gives

e−kdνA,B(k)=dνB,A(−k).e^{-k}d\nu_{A,B}(k)=d\nu_{B,A}(-k).

With Fourier kernel e+iske^{+isk}, the variable called kk is replaced by −k-k; the exponential is correspondingly written e+ke^{+k}. This is a convention change, not a different physical law.

3. The accelerated strip and its forbidden crossing

Section titled “3. The accelerated strip and its forbidden crossing”

List the poles of G(z)=−a2/[16π2sinh⁡2(az/2)]\mathcal G(z)=-a^2/[16\pi^2\sinh^2(az/2)]. Show that the strip −2π/a<Im⁡z<0-2\pi/a<\operatorname{Im}z<0 contains no pole, derive its two boundary prescriptions, and explain why a contour crossing z=0z=0 cannot prove the KMS identity.

Solution

The zeros of sinh⁡(az/2)\sinh(az/2) occur at az/2=iπnaz/2=i\pi n, so the double poles are

z=2πina,n∈Z.z=\frac{2\pi i n}{a}, \qquad n\in\mathbb Z.

Only the boundary lines at imaginary parts 00 and −2π/a-2\pi/a contain poles; the open strip is analytic. Its upper boundary is

lim⁡ϵ↓0G(u−iϵ)=G+(u).\lim_{\epsilon\downarrow0}\mathcal G(u-i\epsilon)=G^+(u).

Approaching the lower edge from inside and using sinh⁡(w−iπ)=−sinh⁡w\sinh(w-i\pi)=-\sinh w gives

lim⁡ϵ↓0G(u−2πi/a+iϵ)=lim⁡ϵ↓0G(u+iϵ)=G−(u).\lim_{\epsilon\downarrow0} \mathcal G(u-2\pi i/a+i\epsilon) =\lim_{\epsilon\downarrow0}\mathcal G(u+i\epsilon) =G^-(u).

The two i0i0 prescriptions encode opposite operator orders. A contour through z=0z=0 crosses a double pole, so Cauchy continuation is unavailable and an additional commutator/contact distribution is picked up. Excluding ∣u∣<umin⁡\lvert u\rvert<u_{\min} in a pointwise plot does not remove this distributional issue; a valid test keeps smeared observables and approaches each edge from inside.

  • Haag, Rudolf, Nico M. Hugenholtz, and Marinus Winnink. “On the Equilibrium States in Quantum Statistical Mechanics.” Communications in Mathematical Physics 5 (1967): 215–236. DOI.
  • Takesaki, Masamichi. Tomita’s Theory of Modular Hilbert Algebras and Its Applications. Lecture Notes in Mathematics 128. Berlin: Springer, 1970. DOI.
  • Unruh, William G. “Notes on Black-Hole Evaporation.” Physical Review D 14 (1976): 870–892. DOI.

Original QFT.org content:CC BY 4.0, unless an item supplies different terms. Third-party material retains its own terms.