Skip to content

From Euclidean Time to Euclidean QFT

Imaginary time performs two closely related jobs in quantum mechanics. Evolution by e−Hτ/ℏe^{-H\tau/\hbar} suppresses high-energy components, while a trace of e−βHe^{-\beta H} turns the imaginary-time interval into a closed path of circumference βℏ\beta\hbar. In quantum field theory, the same operator statements remain true, but each time slice now carries an entire field configuration.

That replacement produces the basic Euclidean-QFT dictionary:

q(τ)⟶ϕ(τ,x),∫Dq e−SE[q]/ℏ⟶∫Dϕ e−SE[ϕ]/ℏ,τ∈[0,βℏ]⟶Sβℏ1×Σ.\begin{aligned} q(\tau) &\longrightarrow \phi(\tau,\mathbf x), \\ \int\mathcal Dq\, e^{-S_E[q]/\hbar} &\longrightarrow \int\mathcal D\phi\, e^{-S_E[\phi]/\hbar}, \\ \tau\in[0,\beta\hbar] &\longrightarrow S^1_{\beta\hbar}\times\Sigma . \end{aligned}

The imaginary-time kernel, Wick rotation, and ground-state projection belong to Euclidean and Imaginary-Time Path Integrals. The closed thermal trace, cyclic coordinate measure, permutation sectors, ring-polymer map, and periodic oscillator determinant belong to Path Integrals for Statistical Mechanics. This page takes those results as input and explains what changes when the dynamical variable is a field, when particle statistics determine temporal boundary conditions, and when Euclidean correlation functions must be related back to Lorentzian observables.

Throughout, β=1/(kBT)\beta=1/(k_B T) and c=1c=1. Factors of ℏ\hbar are retained where they clarify the length βℏ\beta\hbar of the thermal circle.

Consider a real scalar field with Lorentzian action

SM[ϕ]=∫dt ddx [12(∂tϕ)2−12(∇ϕ)2−V(ϕ)].S_M[\phi] = \int dt\,d^d x\, \left[ \frac12(\partial_t\phi)^2 - \frac12(\boldsymbol\nabla\phi)^2 - V(\phi) \right].

Under a justified contour rotation t=−iτt=-i\tau, the exponent transforms schematically as

iℏSM⟶−1ℏSE,\frac{i}{\hbar}S_M \quad\longrightarrow\quad -\frac{1}{\hbar}S_E,

with Euclidean action

SE[ϕ]=∫dτ ddx [12(∂τϕ)2+12(∇ϕ)2+V(ϕ)].S_E[\phi] = \int d\tau\,d^d x\, \left[ \frac12(\partial_\tau\phi)^2 + \frac12(\boldsymbol\nabla\phi)^2 + V(\phi) \right].

For

V(ϕ)=12m2ϕ2+λ4!ϕ4,V(\phi) = \frac12m^2\phi^2 + \frac{\lambda}{4!}\phi^4,

the Euclidean kinetic and potential terms are nonnegative when m2≥0m^2\ge0 and λ≥0\lambda\ge0. This makes the formal weight damping rather than oscillatory:

eiSM/ℏ⟶e−SE/ℏ.e^{iS_M/\hbar} \quad\longrightarrow\quad e^{-S_E/\hbar}.

The word “formal” matters. A Wick rotation is a contour deformation in complex time or energy, not a universal symbol substitution. Singularities, boundary conditions, the state being prepared, and the large-field behavior of the action determine whether the deformation is legitimate. The construction of the field measure and its regulator belongs to From Path Integrals in QM to Field Path Integrals.

Vacuum Projection Becomes a Statement About Wavefunctionals

Section titled “Vacuum Projection Becomes a Statement About Wavefunctionals”

Let ∣φ⟩\lvert\varphi\rangle denote a generalized field eigenstate on a spatial slice:

ϕ^(x)∣φ⟩=φ(x)∣φ⟩.\widehat\phi(\mathbf x) \lvert\varphi\rangle = \varphi(\mathbf x) \lvert\varphi\rangle.

The Euclidean field kernel is

KE[φf,T;φi,0]:=⟨φf∣e−HT/ℏ∣φi⟩=∫ϕ(0,x)=φi(x)ϕ(T,x)=φf(x)Dϕ e−SE[ϕ]/ℏ.\begin{aligned} K_E[\varphi_f,\mathcal T;\varphi_i,0] &:= \langle\varphi_f\rvert e^{-H\mathcal T/\hbar} \lvert\varphi_i\rangle \\ &= \int_{\substack{ \phi(0,\mathbf x)=\varphi_i(\mathbf x)\\ \phi(\mathcal T,\mathbf x)=\varphi_f(\mathbf x) }} \mathcal D\phi\, e^{-S_E[\phi]/\hbar}. \end{aligned}

In a regulated finite-volume theory with energy eigenstates ∣n⟩\lvert n\rangle, its spectral expansion is

KE=∑nΨn[φf] Ψn∗[φi] e−EnT/ℏ,K_E = \sum_n \Psi_n[\varphi_f]\, \Psi_n^*[\varphi_i]\, e^{-E_n\mathcal T/\hbar},

where

Ψn[φ]:=⟨φ∣n⟩\Psi_n[\varphi] := \langle\varphi\vert n\rangle

is a Schrödinger-picture wavefunctional. If the vacuum is nondegenerate and the boundary data overlap it, then

eE0T/ℏKE⟶Ψ0[φf]Ψ0∗[φi](T→∞).e^{E_0\mathcal T/\hbar} K_E \longrightarrow \Psi_0[\varphi_f]\Psi_0^*[\varphi_i] \qquad (\mathcal T\to\infty).

Thus a long Euclidean cylinder prepares the vacuum just as long imaginary-time evolution projects onto the ground state in ordinary quantum mechanics. A half-space path integral can prepare a vacuum wavefunctional on its boundary. Degenerate vacua, gapless spectra, infinite volume, and topological sectors require more care: the limit may project onto a subspace, depend on boundary conditions, or converge without a simple isolated exponential.

The Thermal Trace Becomes a Thermal Circle

Section titled “The Thermal Trace Becomes a Thermal Circle”

For a scalar field, the canonical partition function is

Z(β)=Tr⁡e−βH.Z(\beta) = \operatorname{Tr}e^{-\beta H}.

Taking the trace identifies the initial and final field configurations:

Z(β)=∫Dφ ⟨φ∣e−βH∣φ⟩=∫ϕ(βℏ,x)=ϕ(0,x)Dϕ e−SE[ϕ]/ℏ.\begin{aligned} Z(\beta) &= \int\mathcal D\varphi\, \langle\varphi\rvert e^{-\beta H} \lvert\varphi\rangle \\ &= \int_{\phi(\beta\hbar,\mathbf x) = \phi(0,\mathbf x)} \mathcal D\phi\, e^{-S_E[\phi]/\hbar}. \end{aligned}

Euclidean time is therefore compact:

τ∼τ+βℏ.\tau \sim \tau+\beta\hbar.

If space is Σ\Sigma, finite-temperature QFT lives on the Euclidean geometry

Sβℏ1×Σ.S^1_{\beta\hbar}\times\Sigma.

This statement cleanly separates three common constructions:

Euclidean constructionGeometry or boundary dataQuantity prepared
Fixed-boundary kernelfinite interval with φi,φf\varphi_i,\varphi_ftransition kernel in imaginary time
Vacuum projectionlong cylinder or half-spaceground-state wavefunctional or vacuum correlator
Thermal tracecompact circle of length βℏ\beta\hbarpartition function and thermal correlators

The zero-temperature limit sends β→∞\beta\to\infty, so the thermal circle decompactifies. Vacuum Euclidean QFT and thermal Euclidean QFT are therefore related, but they are not identical formulations at finite β\beta.

Statistics Fixes the Temporal Boundary Condition

Section titled “Statistics Fixes the Temporal Boundary Condition”

For bosonic fields, the trace produces periodic boundary conditions:

ϕ(τ+βℏ,x)=ϕ(τ,x).\phi(\tau+\beta\hbar,\mathbf x) = \phi(\tau,\mathbf x).

For fermionic fields represented by Grassmann coherent states, the trace produces antiperiodic boundary conditions:

ψ(τ+βℏ,x)=−ψ(τ,x),ψ‾(τ+βℏ,x)=−ψ‾(τ,x).\begin{aligned} \psi(\tau+\beta\hbar,\mathbf x) &= -\psi(\tau,\mathbf x), \\ \overline\psi(\tau+\beta\hbar,\mathbf x) &= -\overline\psi(\tau,\mathbf x). \end{aligned}

The minus sign is not an arbitrary convention and does not mean that a physical fermion state literally changes sign after elapsed imaginary time. It arises from the graded coherent-state trace and is what gives thermal fermion propagators their Fermi–Dirac structure. Coherent-State Path Integrals Preview derives the finite-mode trace formula and its determinant check. Bosonic and fermionic composite operators inherit boundary behavior from their total fermion parity.

Gauge fields are bosonic and are ordinarily periodic around the thermal circle, but gauge transformations, holonomy around the circle, and gauge fixing add structure not present for a scalar field. Those issues belong to finite-temperature gauge theory rather than to this bridge.

Compact Euclidean time makes temporal frequency discrete. A periodic bosonic field has Fourier expansion

ϕ(τ,x)=∑n∈Ze−iωnτϕn(x),\phi(\tau,\mathbf x) = \sum_{n\in\mathbb Z} e^{-i\omega_n\tau} \phi_n(\mathbf x),

with bosonic Matsubara frequencies

ωn=2πnβℏ.\omega_n = \frac{2\pi n}{\beta\hbar}.

An antiperiodic fermionic field has

ψ(τ,x)=∑n∈Ze−iνnτψn(x),\psi(\tau,\mathbf x) = \sum_{n\in\mathbb Z} e^{-i\nu_n\tau} \psi_n(\mathbf x),

where

νn=(2n+1)πβℏ.\nu_n = \frac{(2n+1)\pi}{\beta\hbar}.

Two consequences are immediate:

  • bosons possess a zero-frequency mode n=0n=0;
  • fermions have no zero-frequency mode.

As β→∞\beta\to\infty, the frequency spacing vanishes and a bosonic sum approaches a continuous integral:

1βℏ∑nf(ωn)⟶∫−∞∞dω2π f(ω).\frac{1}{\beta\hbar} \sum_n f(\omega_n) \longrightarrow \int_{-\infty}^{\infty} \frac{d\omega}{2\pi}\, f(\omega).

At high temperature or long spatial distance, the bosonic zero mode can dominate static observables. This is the seed of dimensional reduction: under suitable scale-separation assumptions, a thermal field theory in d+1d+1 dimensions admits a dd-dimensional effective description for its longest-distance modes. Interactions, gauge constraints, infrared divergences, and matching coefficients decide whether that description is quantitatively reliable.

For an operator OO, define imaginary-time evolution by

O(τ)=eτH/ℏOe−τH/ℏ.O(\tau) = e^{\tau H/\hbar} O e^{-\tau H/\hbar}.

A thermal Euclidean two-point function is

GE(τ)=1Z(β)Tr⁡[e−βHTτO(τ)O(0)].G_E(\tau) = \frac{1}{Z(\beta)} \operatorname{Tr} \left[ e^{-\beta H} \mathcal T_\tau O(\tau)O(0) \right].

Cyclicity of the trace gives the Kubo–Martin–Schwinger relation. For bosonic operators, the imaginary-time-ordered correlator is periodic in βℏ\beta\hbar; a fermionic two-point function is antiperiodic. These are statements about the ordered thermal correlator, not merely about solving a differential equation on a circle.

At zero temperature, Euclidean vacuum correlation functions are often called Schwinger functions:

Sn(xE1,…,xEn)=⟨0∣Tτ{ϕ(xE1)⋯ϕ(xEn)}∣0⟩.S_n(x_{E1},\ldots,x_{En}) = \langle0\rvert \mathcal T_\tau \left\{ \phi(x_{E1})\cdots\phi(x_{En}) \right\} \lvert0\rangle.

Formally, they have the normalized path-integral representation

Sn=∫Dϕ ϕ(xE1)⋯ϕ(xEn)e−SE[ϕ]/ℏ∫Dϕ e−SE[ϕ]/ℏ.S_n = \frac{ \displaystyle \int\mathcal D\phi\, \phi(x_{E1})\cdots\phi(x_{En}) e^{-S_E[\phi]/\hbar} }{ \displaystyle \int\mathcal D\phi\, e^{-S_E[\phi]/\hbar} }.

The denominator removes vacuum normalization factors. The regulator, boundary conditions, and renormalization prescription remain part of the definition.

From Correlation Functions to QFT Observables explains how Euclidean decay, spectral reconstruction, and continuum extrapolation convert these functions into physical information.

Worked Example: Free Scalar Thermal Propagator

Section titled “Worked Example: Free Scalar Thermal Propagator”

For this example set ℏ=c=kB=1\hbar=c=k_B=1. A free scalar field on the thermal circle has Euclidean action

SE=12∫0βdτ∫ddx [(∂τϕ)2+(∇ϕ)2+m2ϕ2].S_E = \frac12 \int_0^\beta d\tau \int d^d x\, \left[ (\partial_\tau\phi)^2 + (\boldsymbol\nabla\phi)^2 + m^2\phi^2 \right].

Let

Ek:=k2+m2.E_{\mathbf k} := \sqrt{\mathbf k^2+m^2}.

The momentum-space Euclidean propagator is

GE(iωn,k)=1ωn2+Ek2,ωn=2πnβ.G_E(i\omega_n,\mathbf k) = \frac{1} {\omega_n^2+E_{\mathbf k}^2}, \qquad \omega_n = \frac{2\pi n}{\beta}.

For 0≤τ≤β0\le\tau\le\beta, summing the Matsubara modes gives

GE(τ,k)=12Ek[(1+nB(Ek))e−Ekτ+nB(Ek)eEkτ],\begin{aligned} G_E(\tau,\mathbf k) &= \frac{1}{2E_{\mathbf k}} \big[ \big(1+n_B(E_{\mathbf k})\big) e^{-E_{\mathbf k}\tau} \\ &\qquad\qquad + n_B(E_{\mathbf k}) e^{E_{\mathbf k}\tau} \big], \end{aligned}

where

nB(E):=1eβE−1.n_B(E) := \frac{1}{e^{\beta E}-1}.

The two terms describe propagation around the thermal circle in the two possible orientations. The expression agrees at τ=0\tau=0 and τ=β\tau=\beta, as periodicity requires. In the zero-temperature limit, nB(E)→0n_B(E)\to0 and

GE(τ,k)⟶e−Ekτ2Ek(τ≥0).G_E(\tau,\mathbf k) \longrightarrow \frac{e^{-E_{\mathbf k}\tau}} {2E_{\mathbf k}} \qquad (\tau\ge0).

Exponential Euclidean decay therefore reveals the excitation energy. In an interacting theory, long-distance decay similarly carries spectral information, although multiparticle thresholds, anomalous dimensions, finite volume, and operator overlap complicate the extraction.

The Euclidean weight resembles a Gibbs weight:

e−SE[ϕ]/ℏ↔e−βclHcl[ϕ].e^{-S_E[\phi]/\hbar} \quad\leftrightarrow\quad e^{-\beta_{\rm cl}\mathcal H_{\rm cl}[\phi]}.

This is more than a visual analogy. Euclidean QFT methods and classical statistical-field-theory methods share partition functions, correlation functions, coarse graining, critical points, and renormalization-group flows. A dd-dimensional quantum system at finite temperature becomes a Euclidean problem with dd spatial dimensions and one compact imaginary-time direction. A zero-temperature relativistic theory has an unbounded Euclidean-time direction and, when continuation is valid, Euclidean rotational symmetry.

The analogy has limits. The formal weight need not define a positive probability measure. Fermion determinants, finite density, topological terms, and some gauge-theory formulations can make the effective weight negative or complex. Even when the regulated measure is positive, taking a continuum limit and reconstructing a unitary Lorentzian theory are additional problems.

Statistical Field Theory Preview develops the classical field-measure side of this correspondence, including coarse-field definitions, sources, fluctuations, dimensional reduction, and the distinction between static measures and dynamics.

Euclidean correlators encode real-time information through analyticity and spectral representations. In one common convention,

GE(iωn,k)=∫−∞∞dω′2π ρ(ω′,k)iωn−ω′,G_E(i\omega_n,\mathbf k) = \int_{-\infty}^{\infty} \frac{d\omega'}{2\pi}\, \frac{\rho(\omega',\mathbf k)} {i\omega_n-\omega'},

while the retarded correlator is

GR(ω,k)=∫−∞∞dω′2π ρ(ω′,k)ω−ω′+i0.G_R(\omega,\mathbf k) = \int_{-\infty}^{\infty} \frac{d\omega'}{2\pi}\, \frac{\rho(\omega',\mathbf k)} {\omega-\omega'+i0}.

This motivates the continuation

iωn⟶ω+i0,i\omega_n \longrightarrow \omega+i0,

after the Euclidean function has been identified as the boundary value of the appropriate analytic function. Signs and factors of ii vary with Green-function conventions, so the spectral definition must accompany the continuation.

Exact analytic continuation is a structural relation. Numerical continuation from finitely many noisy Euclidean data points is an ill-conditioned inverse problem and requires additional information or regularization. Euclidean Monte Carlo data therefore does not automatically provide precise real-time spectra or transport coefficients.

A collection of Euclidean-invariant functions is not automatically the analytic continuation of a unitary Lorentzian QFT. One central requirement is reflection positivity. Let Θ\Theta reflect Euclidean time, τ↦−τ\tau\mapsto-\tau, together with complex conjugation. For suitable functionals FF supported at positive Euclidean time, reflection positivity requires

⟨(ΘF)F⟩E≥0.\langle (\Theta F)F \rangle_E \ge 0.

This condition is the Euclidean shadow of a positive Hilbert-space inner product. Together with regularity, symmetry, Euclidean invariance, and clustering assumptions, the Osterwalder–Schrader framework gives conditions under which Euclidean Schwinger functions reconstruct a Lorentzian quantum field theory.

Reflection positivity is not a routine consequence of writing e−SEe^{-S_E}. It can fail for an arbitrary discretization, an effective nonlocal action, or a complex weight. Gauge-fixed and fermionic theories also require formulations adapted to their constraints and grading. The reconstruction theorem is therefore a boundary on informal Wick-rotation arguments, not merely an optional mathematical refinement.

This page owns the transition from imaginary-time quantum mechanics to Euclidean field theory:

  • field wavefunctionals and vacuum projection;
  • the thermal circle;
  • periodic and antiperiodic temporal boundary conditions;
  • the field-theory realization of bosonic and fermionic Matsubara modes;
  • the relation between Euclidean correlators, statistical field theory, and Lorentzian reconstruction.

It does not rederive the quantum-mechanical Euclidean kernel, define the regulated field measure, or develop full thermal QFT. Use Bosonic and Fermionic Matsubara Frequencies for the QM-level grids, units, and index conventions; Matsubara Formalism Preview for the compact-time transform, frequency sums, and convergence checks; Finite-Temperature QFT Bridge for thermal propagators, loop sums, screening, zero-mode matching, and the Euclidean-to-real-time boundary; Euclidean and Imaginary-Time Path Integrals for the canonical QM derivation; From Path Integrals in QM to Field Path Integrals for the field configuration space and measure; and From Propagators in QM to Propagators in QFT for the distinctions among Feynman, retarded, and Euclidean two-point functions.

  • Treating Wick rotation as the substitution t=−iτt=-i\tau without checking singularities or boundary data.
  • Confusing a fixed-boundary Euclidean kernel, a vacuum path integral, and a thermal trace.
  • Assigning periodic boundary conditions to thermal fermion fields.
  • Forgetting the factor of ℏ\hbar in the thermal circumference βℏ\beta\hbar when not using natural units.
  • Assuming e−SEe^{-S_E} always defines a positive probability measure.
  • Reading a Matsubara correlator as a retarded correlator without analytic continuation and a convention check.
  • Assuming numerical Euclidean data determines real-time spectral functions uniquely.
  • Believing Euclidean regularization removes the need for renormalization or a continuum-limit analysis.
  • J. Schwinger, “Euclidean Quantum Electrodynamics,” Physical Review 115, 721–731, 1959, doi:10.1103/PhysRev.115.721.
  • T. Matsubara, “A New Approach to Quantum-Statistical Mechanics,” Progress of Theoretical Physics 14, 351–378, 1955, doi:10.1143/PTP.14.351.
  • P. C. Martin and J. Schwinger, “Theory of Many-Particle Systems. I,” Physical Review 115, 1342–1373, 1959, doi:10.1103/PhysRev.115.1342.
  • K. Osterwalder and R. Schrader, “Axioms for Euclidean Green’s Functions,” Communications in Mathematical Physics 31, 83–112, 1973, doi:10.1007/BF01645738.
  • K. Osterwalder and R. Schrader, “Axioms for Euclidean Green’s Functions II,” Communications in Mathematical Physics 42, 281–305, 1975, doi:10.1007/BF01608978.
  • M. Le Bellac, Thermal Field Theory, Cambridge University Press, 1996.
  • J. I. Kapusta and C. Gale, Finite-Temperature Field Theory: Principles and Applications, 2nd ed., Cambridge University Press, 2006.
  • J. Glimm and A. Jaffe, Quantum Physics: A Functional Integral Point of View, 2nd ed., Springer, 1987.

1. Vacuum wavefunctional from a long cylinder

Section titled “1. Vacuum wavefunctional from a long cylinder”

Starting from the spectral representation of KE[φf,T;φi,0]K_E[\varphi_f,\mathcal T;\varphi_i,0], show that a long Euclidean cylinder isolates the vacuum wavefunctional when the vacuum is nondegenerate and both boundaries have nonzero vacuum overlap. State what changes for a gg-fold degenerate ground space.

Solution

Factor out the lowest exponential:

KE=e−E0T/ℏ[Ψ0[φf]Ψ0∗[φi]+∑n>0Ψn[φf]Ψn∗[φi]e−(En−E0)T/ℏ].\begin{aligned} K_E &= e^{-E_0\mathcal T/\hbar} \bigg[ \Psi_0[\varphi_f]\Psi_0^*[\varphi_i] \\ &\qquad + \sum_{n\gt0} \Psi_n[\varphi_f]\Psi_n^*[\varphi_i] e^{-(E_n-E_0)\mathcal T/\hbar} \bigg]. \end{aligned}

If En−E0>0E_n-E_0\gt0 for the excited states and the vacuum overlaps are nonzero, multiplying by eE0T/ℏe^{E_0\mathcal T/\hbar} and taking T→∞\mathcal T\to\infty leaves

Ψ0[φf]Ψ0∗[φi].\Psi_0[\varphi_f]\Psi_0^*[\varphi_i].

For a gg-fold degenerate ground space, the limit is the ground-space projector kernel

∑a=1gΨ0a[φf]Ψ0a∗[φi].\sum_{a=1}^{g} \Psi_{0a}[\varphi_f] \Psi_{0a}^*[\varphi_i].

Boundary conditions or small symmetry-breaking sources may select a particular linear combination only after the relevant limits are specified.

Insert a mode e−iωτe^{-i\omega\tau} into periodic and antiperiodic boundary conditions on a circle of length βℏ\beta\hbar. Derive the two allowed frequency sets.

Solution

Periodicity requires

e−iω(τ+βℏ)=e−iωτ,e^{-i\omega(\tau+\beta\hbar)} = e^{-i\omega\tau},

so

e−iωβℏ=1.e^{-i\omega\beta\hbar}=1.

Therefore

ω=2πnβℏ,n∈Z.\omega = \frac{2\pi n}{\beta\hbar}, \qquad n\in\mathbb Z.

Antiperiodicity instead requires

e−iνβℏ=−1,e^{-i\nu\beta\hbar}=-1,

which gives

ν=(2n+1)πβℏ.\nu = \frac{(2n+1)\pi}{\beta\hbar}.

The periodic set contains n=0n=0; the antiperiodic set does not.

In natural units, evaluate

GE(τ,k)=∫−∞∞dω2π eiωτω2+Ek2.G_E(\tau,\mathbf k) = \int_{-\infty}^{\infty} \frac{d\omega}{2\pi}\, \frac{e^{i\omega\tau}} {\omega^2+E_{\mathbf k}^2}.

Interpret the large-∣τ∣\lvert\tau\rvert behavior.

Solution

For τ>0\tau\gt0, close the contour in the upper half-plane. The pole at ω=iEk\omega=iE_{\mathbf k} gives

GE(τ,k)=e−Ekτ2Ek.G_E(\tau,\mathbf k) = \frac{e^{-E_{\mathbf k}\tau}} {2E_{\mathbf k}}.

For τ<0\tau\lt0, closing in the lower half-plane gives the same result with τ\tau replaced by ∣τ∣\lvert\tau\rvert:

GE(τ,k)=e−Ek∣τ∣2Ek.G_E(\tau,\mathbf k) = \frac{e^{-E_{\mathbf k}\lvert\tau\rvert}} {2E_{\mathbf k}}.

The slowest exponential at large Euclidean separation identifies the lowest energy coupled to the chosen field operator at momentum k\mathbf k. In an interacting theory, the coefficient measures operator overlap and additional exponentials or continua represent higher states.

For the free scalar result

GE(τ,k)=(1+nB)e−Eτ+nBeEτ2E,0≤τ≤β,G_E(\tau,\mathbf k) = \frac{ \big(1+n_B\big)e^{-E\tau} + n_B e^{E\tau} }{2E}, \qquad 0\le\tau\le\beta,

use nB=(eβE−1)−1n_B=(e^{\beta E}-1)^{-1} to verify GE(β,k)=GE(0,k)G_E(\beta,\mathbf k)=G_E(0,\mathbf k). Then take β→∞\beta\to\infty.

Solution

The Bose factor obeys

(1+nB)e−βE=nB,nBeβE=1+nB.\big(1+n_B\big)e^{-\beta E} = n_B, \qquad n_B e^{\beta E} = 1+n_B.

Hence

GE(β,k)=1+2nB2E=GE(0,k).G_E(\beta,\mathbf k) = \frac{1+2n_B}{2E} = G_E(0,\mathbf k).

As β→∞\beta\to\infty, nB→0n_B\to0. For fixed τ≥0\tau\ge0,

GE(τ,k)→e−Eτ2E,G_E(\tau,\mathbf k) \to \frac{e^{-E\tau}}{2E},

the vacuum Euclidean propagator.

Explain why a bosonic field can have a static Matsubara mode while a thermal fermion cannot. What does this suggest about their long-distance roles at temperatures large compared with the spatial momentum scale?

Solution

Bosonic frequencies are ωn=2πn/β\omega_n=2\pi n/\beta, so ω0=0\omega_0=0. Fermionic frequencies are νn=(2n+1)π/β\nu_n=(2n+1)\pi/\beta, whose smallest magnitude is π/β\pi/\beta. A bosonic zero mode can therefore remain light when all nonzero temporal modes cost energies of order TT.

At spatial momenta much smaller than TT, nonzero Matsubara modes can often be integrated out, leaving an effective theory for static bosonic fields in one fewer dimension. Fermions have no static mode and generally contribute through matching coefficients rather than as long-distance thermal degrees of freedom. This argument identifies a possible effective description; it does not by itself prove that perturbative matching or infrared dynamics is under control.

6. Reflection positivity for a free scalar

Section titled “6. Reflection positivity for a free scalar”

Let FF be linear in a free scalar field and supported at τ>0\tau\gt0. In spatial momentum space, show that the reflected quadratic form has the structure

⟨(ΘF)F⟩E=∫ddk(2π)d12Ek∣∫0∞dτ e−Ekτf(τ,k)∣2.\langle(\Theta F)F\rangle_E = \int\frac{d^d k}{(2\pi)^d} \frac{1}{2E_{\mathbf k}} \left| \int_0^\infty d\tau\, e^{-E_{\mathbf k}\tau} f(\tau,\mathbf k) \right|^2.

Why is this useful?

Solution

The free mixed-representation propagator is

GE(τ−τ′,k)=e−Ek∣τ−τ′∣2Ek.G_E(\tau-\tau',\mathbf k) = \frac{ e^{-E_{\mathbf k}\lvert\tau-\tau'\rvert} }{2E_{\mathbf k}}.

Reflection sends the first time argument to −τ-\tau. For τ,τ′>0\tau,\tau'\gt0,

∣−τ−τ′∣=τ+τ′,\lvert-\tau-\tau'\rvert = \tau+\tau',

so the kernel factorizes:

GE(−τ−τ′,k)=e−Ekτe−Ekτ′2Ek.G_E(-\tau-\tau',\mathbf k) = \frac{ e^{-E_{\mathbf k}\tau} e^{-E_{\mathbf k}\tau'} }{2E_{\mathbf k}}.

Inserting this into the smeared two-point function gives the stated absolute square, which is nonnegative. The factorization exhibits how Euclidean reflection defines a positive inner product and prepares the Hilbert-space interpretation needed for Lorentzian reconstruction.