Skip to content

Phase-Space Formulation from QM to QFT

For a mechanical system, phase space is coordinatized by canonical pairs (qi,pi)(q^i,p_i). For a classical field, the corresponding data at one time are a field configuration and its canonical momentum:

ϕ(x),π(x).\phi(\mathbf x), \qquad \pi(\mathbf x).

This replacement is the shortest phase-space bridge from quantum mechanics to field theory:

(qi,pi)⟶(ϕ(x),π(x)).(q^i,p_i) \quad \longrightarrow \quad \bigl( \phi(\mathbf x), \pi(\mathbf x) \bigr).

The longer bridge requires a regulator, infinitely many degrees of freedom, local operator products, constraints, and renormalization. Phase-space functions become functionals, ordinary derivatives become functional derivatives, and finite-dimensional measures become regulated functional measures.

This page is an orientation. The canonical mechanical tools remain Weyl Transform, Wigner Function, and Phase-Space Dynamics. The detailed construction of field path integrals belongs beyond this quantum-mechanical bridge.

Classical Phase Space in Mechanics and Fields

Section titled “Classical Phase Space in Mechanics and Fields”

For NN canonical coordinates, a classical state is a point

z=(q1,…,qN,p1,…,pN).z = \left( q^1,\ldots,q^N, p_1,\ldots,p_N \right).

The Poisson bracket is

{F,G}PB=∑i=1N(∂F∂qi∂G∂pi−∂F∂pi∂G∂qi).\{F,G\}_{\mathrm{PB}} = \sum_{i=1}^{N} \left( \frac{\partial F}{\partial q^i} \frac{\partial G}{\partial p_i} - \frac{\partial F}{\partial p_i} \frac{\partial G}{\partial q^i} \right).

A scalar field is the continuum limit of a many-coordinate system. On a spatial lattice with NN sites, the field values ϕn\phi_n are NN coordinates and their conjugate momenta πn\pi_n supply the other NN phase-space coordinates. The regulated phase space is therefore 2N2N-dimensional.

In continuum notation, a functional F[ϕ,π]F[\phi,\pi] has field Poisson bracket

{F,G}PB=∫ddx[δFδϕ(x)δGδπ(x)−δFδπ(x)δGδϕ(x)].\begin{aligned} \{F,G\}_{\mathrm{PB}} &= \int d^d x \left[ \frac{\delta F}{\delta\phi(\mathbf x)} \frac{\delta G}{\delta\pi(\mathbf x)} \right. \\ &\qquad\left. - \frac{\delta F}{\delta\pi(\mathbf x)} \frac{\delta G}{\delta\phi(\mathbf x)} \right]. \end{aligned}

Hamilton’s field equations are

ϕ˙(x)=δHδπ(x),π˙(x)=−δHδϕ(x).\dot\phi(\mathbf x) = \frac{\delta H}{\delta\pi(\mathbf x)}, \qquad \dot\pi(\mathbf x) = - \frac{\delta H}{\delta\phi(\mathbf x)}.

The continuum formulas summarize a regulated limit. Boundary conditions, spatial topology, lattice spacing or momentum cutoff, and the definition of the canonical momentum are part of the phase-space structure.

A field configuration at one time is not a spacetime history. The canonical phase space labels instantaneous data; a path integral integrates over trajectories through that phase space.

Canonical quantization replaces the field Poisson bracket with equal-time commutators. For a real bosonic scalar field,

[ϕ^(t,x),π^(t,y)]=iℏδ(d)(x−y),\left[ \hat\phi(t,\mathbf x), \hat\pi(t,\mathbf y) \right] = i\hbar \delta^{(d)}(\mathbf x-\mathbf y),

with

[ϕ^(t,x),ϕ^(t,y)]=0,[π^(t,x),π^(t,y)]=0.\left[ \hat\phi(t,\mathbf x), \hat\phi(t,\mathbf y) \right] = 0, \qquad \left[ \hat\pi(t,\mathbf x), \hat\pi(t,\mathbf y) \right] = 0.

The delta distribution is the continuum counterpart of the Kronecker delta for finitely many canonical pairs. Its normalization depends on the Fourier and finite-volume conventions.

For a free scalar field, the Hamiltonian is schematically

H=12∫ddx [π2+(∇ϕ)2+m2ϕ2].H = \frac12 \int d^d x\, \left[ \pi^2 + (\nabla\phi)^2 + m^2\phi^2 \right].

After a mode decomposition and suitable reality conditions, each independent mode behaves like a harmonic oscillator:

H∼∑k[12pk2+12ωk2qk2].H \sim \sum_{\mathbf k} \left[ \frac12p_{\mathbf k}^2 + \frac12\omega_{\mathbf k}^2q_{\mathbf k}^2 \right].

This is the canonical content of Harmonic Oscillator to Fields. The oscillator ladder operators become mode creation and annihilation operators. An interacting field is not a collection of independent oscillators: interaction terms couple modes and modify the vacuum and excitations.

Canonical commutation relations are a starting point, not a complete field theory. Operator-valued distributions cannot generally be multiplied at the same point without regularization, and different continuum settings can admit representation issues absent for finitely many canonical pairs.

A mechanical Hamiltonian path integral has the formal form

K=∫Dp Dq exp⁡[iℏ∫dt (piq˙i−H(p,q))].K = \int \mathcal Dp\,\mathcal Dq\, \exp\left[ \frac{i}{\hbar} \int dt\, \left( p_i\dot q^i-H(p,q) \right) \right].

For a bosonic field, the analogous expression is

Z=∫Dπ Dϕ ×exp⁡{iℏ∫dt ddx [πϕ˙−H(ϕ,π)]}.\begin{aligned} Z &= \int \mathcal D\pi\,\mathcal D\phi\, \\ &\quad\times \exp\left\{ \frac{i}{\hbar} \int dt\,d^d x\, \left[ \pi\dot\phi - \mathcal H(\phi,\pi) \right] \right\}. \end{aligned}

At finite spatial lattice and finite time slicing, this is a many-coordinate phase-space integral. The functional notation is its continuum shorthand.

For a standard quadratic momentum dependence,

H=12π2+V(ϕ,∇ϕ),\mathcal H = \frac12\pi^2 + \mathcal V(\phi,\nabla\phi),

the π\pi integral is Gaussian. Completing the square,

πϕ˙−12π2=−12(π−ϕ˙)2+12ϕ˙2,\pi\dot\phi-\frac12\pi^2 = - \frac12 \left( \pi-\dot\phi \right)^2 + \frac12\dot\phi^2,

formally converts the Hamiltonian integral into a configuration-space integral with Lagrangian density

L=12ϕ˙2−V.\mathcal L = \frac12\dot\phi^2 - \mathcal V.

This equivalence is not automatic for arbitrary systems. Position-dependent kinetic terms, constraints, nontrivial measures, first-order actions, and gauge redundancy can leave determinants or constraint factors. Gauge fields require constraint and gauge-fixing machinery. Fermionic fields use Grassmann variables rather than ordinary commuting phase-space coordinates.

From Path Integrals in QM to Field Path Integrals owns the regulated path-to-field construction. The phase-space point here is that the first-order term πϕ˙\pi\dot\phi records the canonical symplectic structure before momenta are integrated out.

Bosonic coherent states are labeled by complex amplitudes αj\alpha_j and satisfy

∫d2αjπ∣αj⟩⟨αj∣=Ij\int \frac{d^2\alpha_j}{\pi} \lvert\alpha_j\rangle \langle\alpha_j\rvert = I_j

for each mode. Repeatedly inserting these overcomplete identities produces a path integral over complex phase-space labels.

For finitely many modes, the action has the schematic form

S=∫dt [iℏ∑jαj∗α˙j−Hsymb(α∗,α)].S = \int dt\, \left[ i\hbar \sum_j \alpha_j^*\dot\alpha_j - H_{\mathrm{symb}} \left( \alpha^*,\alpha \right) \right].

The Hamiltonian symbol and finite-slice rule depend on operator ordering. Normal, antinormal, and Weyl symbols are not interchangeable. Boundary terms also matter because coherent-state path integrals use an overcomplete, nonorthogonal basis.

For a nonrelativistic bosonic field, the mode amplitudes can be traded for a complex field ψ(x,t)\psi(\mathbf x,t):

S∼∫dt ddx [iℏψ∗∂tψ−H(ψ∗,ψ)].S \sim \int dt\,d^d x\, \left[ i\hbar\psi^*\partial_t\psi - \mathcal H \left( \psi^*,\psi \right) \right].

This language is central in many-body theory, quantum optics, condensate physics, and nonequilibrium methods. It does not mean ψ\psi is a one-particle wavefunction in every use; in a field path integral it is an integration variable associated with field-mode amplitudes.

Fermionic coherent states are labeled by Grassmann variables. Their measures, signs, and thermal boundary conditions differ fundamentally from the bosonic case. Replacing them by ordinary complex numbers destroys fermionic anticommutation.

Coherent States in Phase Space gives the one-mode geometry. Coherent-State Path Integrals Preview gives the finite-slice bosonic and fermionic thermal construction. Second Quantization gives the occupation-number bridge that makes the multimode notation natural.

For a lattice-regulated bosonic field, the density operator can be represented by a finite-dimensional Wigner function over all canonical site or mode variables. The continuum notation is a Wigner functional

WΛ[φ,π],W_\Lambda[\varphi,\pi],

where Λ\Lambda denotes the regulator data.

A formal field-basis definition is

WΛ[φ,π]=NΛ∫Dη exp⁡[−iℏ∫ddx π(x)η(x)]×⟨φ+η2|ρ|φ−η2⟩.\begin{aligned} W_\Lambda[\varphi,\pi] &= \mathcal N_\Lambda \int\mathcal D\eta\, \exp\left[ - \frac{i}{\hbar} \int d^d x\, \pi(\mathbf x)\eta(\mathbf x) \right] \\ &\quad\times \left\langle \varphi+\frac{\eta}{2} \middle| \rho \middle| \varphi-\frac{\eta}{2} \right\rangle. \end{aligned}

The normalization NΛ\mathcal N_\Lambda depends on the regulated measure and Fourier convention. It should not be guessed from a one-coordinate formula after the cutoff has been removed.

With matched conventions, integrating over all momentum fields gives the diagonal field-configuration probability functional:

∫Dπ WΛ[φ,π]=⟨φ∣ρ∣φ⟩.\int\mathcal D\pi\, W_\Lambda[\varphi,\pi] = \langle\varphi\rvert\rho\lvert\varphi\rangle.

As in ordinary quantum mechanics, the full Wigner functional is a quasiprobability and can be negative.

The exact evolution can be written with a functional Moyal bracket:

∂WΛ∂t={HΛ,WΛ}M,func.\frac{\partial W_\Lambda}{\partial t} = \{H_\Lambda,W_\Lambda\}_{M,\mathrm{func}}.

Its leading derivative term is the classical field Poisson bracket, with higher functional derivatives supplying quantum corrections. Dropping those derivatives leads to classical-statistical or truncated-Wigner-type approximations in suitable bosonic regimes.

Such a truncation is not universally controlled. Its accuracy depends on occupation, coupling, initial state, cutoff, observable, and evolution time. Vacuum fluctuations contribute roughly half a quantum per bosonic mode in common Wigner conventions, so increasing the ultraviolet cutoff can inject cutoff-dependent noise unless the approximation and renormalization are handled consistently.

The finite-mode dictionary is direct:

Quantum mechanicsRegulated field theory
canonical pairs (qi,pi)(q_i,p_i)lattice or mode pairs (ϕn,πn)(\phi_n,\pi_n)
[qi,pj]=iℏδij[q_i,p_j]=i\hbar\delta_{ij}equal-time field commutator
Wigner function W(q,p)W(q,p)multimode Wigner function or functional
coherent amplitude αi\alpha_imode amplitude or complex field
star productmultimode or functional star product
Moyal equationfunctional Wigner evolution

The continuum limit adds problems that have no finite-dimensional analogue:

  • infinitely many ultraviolet modes;
  • singular local products of fields;
  • regulator and renormalization dependence;
  • gauge constraints and redundant variables;
  • Grassmann integration for fermions;
  • possible infrared and infinite-volume limits;
  • dependence of particles and modes on the chosen background or basis.

The formal replacement of sums by integrals does not solve these issues. A continuum phase-space functional is trustworthy only after the regulated theory, observables, and limiting procedure are specified.

The most durable ideas carried from quantum mechanics are:

  • canonical coordinates and momenta encode a symplectic structure;
  • quantization replaces Poisson brackets by noncommutative operator algebra;
  • coherent states label bosonic mode phase space;
  • Wigner transforms turn density operators into quasiprobabilities;
  • Moyal brackets organize quantum corrections to classical flow.

The detailed field-theory continuation branches by purpose:

These bridges provide orientation, not a substitute for a regulated field-theory treatment.

  • Identifying a field configuration at one time with an entire spacetime history.
  • Writing continuum functional measures without first specifying a lattice, mode cutoff, or other regulator.
  • Treating each free-field oscillator mode as a distinct particle rather than a mode whose excitations define quanta.
  • Assuming interacting fields remain independent oscillators.
  • Integrating out canonical momenta without checking determinants, constraints, or boundary terms.
  • Using an ordinary complex coherent-state field for fermions instead of Grassmann variables.
  • Treating a Wigner functional as a positive probability distribution.
  • Calling truncated-Wigner evolution exact because it follows classical field equations.
  • Adding half a quantum to infinitely many modes without tracking cutoff dependence.
  • Ignoring gauge fixing, ultraviolet renormalization, or local-operator definitions in the continuum.
  • J. W. Negele and H. Orland, Quantum Many-Particle Systems, Addison-Wesley, 1988.
  • J. R. Klauder and B.-S. Skagerstam, Coherent States: Applications in Physics and Mathematical Physics, World Scientific, 1985.
  • A. Altland and B. Simons, Condensed Matter Field Theory, 2nd ed., Cambridge University Press, 2010.
  • A. Kamenev, Field Theory of Non-Equilibrium Systems, Cambridge University Press, 2011.
  • A. Polkovnikov, “Phase space representation of quantum dynamics,” Annals of Physics 325, 1790–1852 (2010).
  • J. Berges, “Introduction to nonequilibrium quantum field theory,” AIP Conference Proceedings 739, 3–62 (2004).
  1. A scalar field is regulated on a lattice with NN sites. What is the dimension of its canonical phase space, and why must the regulator be retained before taking a continuum limit?
Solution

Each site has one field coordinate ϕn\phi_n and one conjugate momentum πn\pi_n. The regulated phase space therefore has dimension

2N.2N.

At finite NN, Poisson brackets, commutators, measures, and Wigner transforms are finite-dimensional. Sending the lattice spacing to zero while holding physical volume fixed sends NN to infinity. Ultraviolet divergences and regulator-dependent normalization can then appear, so the limiting observables and renormalization conditions must be specified before the continuum notation is treated as defined.

  1. Show that the field Poisson bracket generates Hamilton’s field equations.
Solution

Take the functional F[ϕ,π]=ϕ(y)F[\phi,\pi]=\phi(\mathbf y). Its functional derivatives are

δFδϕ(x)=δ(d)(x−y),δFδπ(x)=0.\frac{\delta F}{\delta\phi(\mathbf x)} = \delta^{(d)}(\mathbf x-\mathbf y), \qquad \frac{\delta F}{\delta\pi(\mathbf x)} = 0.

Therefore

{ϕ(y),H}PB=∫ddx δ(d)(x−y)δHδπ(x)=δHδπ(y).\begin{aligned} \{\phi(\mathbf y),H\}_{\mathrm{PB}} &= \int d^d x\, \delta^{(d)}(\mathbf x-\mathbf y) \frac{\delta H}{\delta\pi(\mathbf x)} \\ &= \frac{\delta H}{\delta\pi(\mathbf y)}. \end{aligned}

Thus ϕ˙={ϕ,H}PB=δH/δπ\dot\phi=\{\phi,H\}_{\mathrm{PB}}=\delta H/\delta\pi. Repeating the calculation with F=π(y)F=\pi(\mathbf y) gives

π˙={π,H}PB=−δHδϕ.\dot\pi = \{\pi,H\}_{\mathrm{PB}} = - \frac{\delta H}{\delta\phi}.
  1. Integrate the momentum dependence of a scalar phase-space path integral formally.
Solution

For

H=12π2+V(ϕ,∇ϕ),\mathcal H = \frac12\pi^2 + \mathcal V(\phi,\nabla\phi),

the first-order integrand is

πϕ˙−H=−12(π−ϕ˙)2+12ϕ˙2−V.\pi\dot\phi-\mathcal H = - \frac12 \left( \pi-\dot\phi \right)^2 + \frac12\dot\phi^2 - \mathcal V.

The regulated integral over the shifted Gaussian variable π−ϕ˙\pi-\dot\phi contributes a normalization or determinant factor. What remains in the exponent is

L=12ϕ˙2−V.\mathcal L = \frac12\dot\phi^2 - \mathcal V.

The simple result depends on the quadratic, field-independent coefficient of π2\pi^2. More general kinetic terms or constraints require additional factors.

  1. Derive the field-configuration marginal of the regulated Wigner functional.
Solution

Start from

WΛ[φ,π]=NΛ∫Dη exp⁡[−iℏ∫ddx πη]×⟨φ+η2|ρ|φ−η2⟩.\begin{aligned} W_\Lambda[\varphi,\pi] &= \mathcal N_\Lambda \int\mathcal D\eta\, \exp\left[ - \frac{i}{\hbar} \int d^d x\,\pi\eta \right] \\ &\quad\times \left\langle \varphi+\frac{\eta}{2} \middle| \rho \middle| \varphi-\frac{\eta}{2} \right\rangle. \end{aligned}

With the normalization matched to the regulated Fourier transform, integration over every momentum variable produces a functional delta:

NΛ∫Dπ exp⁡[−iℏ∫ddx πη]=δΛ[η].\mathcal N_\Lambda \int\mathcal D\pi\, \exp\left[ - \frac{i}{\hbar} \int d^d x\,\pi\eta \right] = \delta_\Lambda[\eta].

The delta functional sets η=0\eta=0, leaving

∫Dπ WΛ[φ,π]=⟨φ∣ρ∣φ⟩.\int\mathcal D\pi\, W_\Lambda[\varphi,\pi] = \langle\varphi\rvert\rho\lvert\varphi\rangle.
  1. Why is truncated-Wigner evolution not automatically exact when the classical field equations are solved without numerical error?
Solution

The exact Wigner-functional equation contains a functional Moyal bracket. Classical field evolution keeps only its leading functional Poisson bracket and drops higher functional derivatives. Solving the retained equation exactly removes numerical integration error but does not restore the omitted quantum terms.

The truncation error depends on occupation, coupling, state gradients, cutoff, observable, and time. Classical flow can also generate fine structure that makes higher derivatives important. Numerical convergence of the truncated equation is therefore necessary but not sufficient for convergence to the exact quantum theory.