Skip to content

Retarded and Advanced Green Functions

A retarded Green function gives the response after a source acts. An advanced Green function gives the complementary boundary prescription in which the response is supported before the source. The local differential equation may be the same; the boundary condition changes the Green function. In nonrelativistic dynamics, temporal retardation does not by itself imply finite spatial propagation; see Causality, Support, and Interpretation in Nonrelativistic QM.

For a Schrödinger-type operator

L=iℏ∂∂t−H,L = i\hbar\frac{\partial}{\partial t} - H,

a time-domain Green function satisfies

LtG(t,t′)=δ(t−t′)I.L_tG(t,t') = \delta(t-t')I.

The retarded and advanced choices differ by support:

GR(t,t′)=0for t<t′,G^R(t,t')=0 \quad \text{for }t\lt t',

while

GA(t,t′)=0for t>t′.G^A(t,t')=0 \quad \text{for }t\gt t'.

The support condition is not a minor add-on. It is the part of the definition that states the physical question being answered.

For a time-independent Hamiltonian,

U(t,t′)=e−iH(t−t′)/ℏU(t,t') = e^{-iH(t-t')/\hbar}

is the unitary evolution operator from t′t' to tt. With a common normalization, the retarded Green operator is

GR(t,t′)=−iℏθ(t−t′)U(t,t′).G^R(t,t') = -\frac{i}{\hbar} \theta(t-t')U(t,t').

It is zero before the source time t′t'. To verify the Green-function equation, use

iℏ∂∂tU(t,t′)=HU(t,t′)i\hbar\frac{\partial}{\partial t}U(t,t') = HU(t,t')

and

ddtθ(t−t′)=δ(t−t′).\frac{d}{dt}\theta(t-t') = \delta(t-t').

The derivative of the step function supplies the delta function, while the derivative of UU cancels the explicit Hamiltonian term. Thus

(iℏ∂∂t−H)GR(t,t′)=δ(t−t′)I.\left( i\hbar\frac{\partial}{\partial t} - H \right) G^R(t,t') = \delta(t-t')I.

If a driven equation is written as

(iℏ∂∂t−H)ψ(t)=η(t),\left( i\hbar\frac{\partial}{\partial t} - H \right)\psi(t) = \eta(t),

then the retarded solution is schematically

ψR(t)=∫−∞∞dt′ GR(t,t′)η(t′).\psi_R(t) = \int_{-\infty}^{\infty}dt'\, G^R(t,t')\eta(t').

Because GR(t,t′)G^R(t,t') vanishes for t<t′t\lt t', only earlier source values contribute to the response at time tt.

The advanced Green operator uses the opposite support condition:

GA(t,t′)=iℏθ(t′−t)U(t,t′).G^A(t,t') = \frac{i}{\hbar} \theta(t'-t)U(t,t').

It also satisfies

(iℏ∂∂t−H)GA(t,t′)=δ(t−t′)I.\left( i\hbar\frac{\partial}{\partial t} - H \right) G^A(t,t') = \delta(t-t')I.

The sign differs from the retarded expression because

ddtθ(t′−t)=−δ(t−t′).\frac{d}{dt}\theta(t'-t) = -\delta(t-t').

Advanced Green functions are not the usual causal response functions for initial-value problems. They are nevertheless important as complementary boundary values, as ingredients in spectral identities, and as part of the analytic structure connecting retarded, advanced, time-ordered, and Euclidean objects.

Retarded and advanced Green functions solve the same local equation. The difference is global:

ObjectSupport conditionTypical use
GR(t,t′)G^R(t,t')zero for t<t′t\lt t'causal response to a source
GA(t,t′)G^A(t,t')zero for t>t′t\gt t'complementary boundary value, formal identities
GF(t,t′)G^F(t,t')time-ordered prescriptionvacuum perturbation theory preview

This is why the phrase “the Green function of LL” is incomplete. One must also specify the boundary prescription. The same differential expression can have retarded, advanced, outgoing, incoming, Feynman, Euclidean, or finite-temperature Green functions.

Assume time-translation invariance, so G(t,t′)=G(t−t′)G(t,t')=G(t-t'). Use the energy transform convention

G(E)=∫−∞∞dt eiEt/ℏG(t).G(E) = \int_{-\infty}^{\infty}dt\, e^{iEt/\hbar}G(t).

For the retarded Green function, the step function restricts the integral to t>0t\gt0. The convergence prescription is

ei(E−H)t/ℏ⟶ei(E−H)t/ℏe−ϵt/ℏ,ϵ>0.e^{i(E-H)t/\hbar} \quad\longrightarrow\quad e^{i(E-H)t/\hbar}e^{-\epsilon t/\hbar}, \qquad \epsilon\gt0.

This gives the boundary value

GR(E)=lim⁡ϵ→0+1E−H+iϵ.G^R(E) = \lim_{\epsilon\to0^+} \frac{1}{E-H+i\epsilon}.

Similarly,

GA(E)=lim⁡ϵ→0+1E−H−iϵ.G^A(E) = \lim_{\epsilon\to0^+} \frac{1}{E-H-i\epsilon}.

Thus the retarded prescription approaches the real energy axis from the upper half-plane, while the advanced prescription approaches from the lower half-plane:

GR(E)=G(E+i0),GA(E)=G(E−i0).G^R(E)=G(E+i0), \qquad G^A(E)=G(E-i0).

The sign of i0i0 is not decoration. It is the Fourier-domain record of the time-support condition. For coordinate-space source equations and outgoing free kernels, see Energy Green Function.

If the Hamiltonian has discrete eigenstates,

H∣n⟩=En∣n⟩,H|n\rangle=E_n|n\rangle,

then

GR(E)=∑n∣n⟩⟨n∣E−En+i0,G^R(E) = \sum_n \frac{|n\rangle\langle n|} {E-E_n+i0},

and

GA(E)=∑n∣n⟩⟨n∣E−En−i0.G^A(E) = \sum_n \frac{|n\rangle\langle n|} {E-E_n-i0}.

The distribution identity

1x+i0=PV⁡1x−iπδ(x)\frac{1}{x+i0} = \operatorname{PV}\frac{1}{x} - i\pi\delta(x)

implies

GR(E)−GA(E)=−2πi δ(E−H).G^R(E)-G^A(E) = -2\pi i\,\delta(E-H).

Equivalently,

−1πIm⁡GR(E)=δ(E−H),-\frac{1}{\pi} \operatorname{Im}G^R(E) = \delta(E-H),

as an operator-valued distribution, interpreted through matrix elements or traces. This is the seed of the density-of-states relation developed in Spectral Representation of Green Functions.

For a one-dimensional Hilbert space with Hamiltonian H=E0H=E_0, the retarded Green function is

GR(t)=−iℏθ(t)e−iE0t/ℏ.G^R(t) = -\frac{i}{\hbar} \theta(t) e^{-iE_0t/\hbar}.

It satisfies

(iℏddt−E0)GR(t)=δ(t).\left( i\hbar\frac{d}{dt} - E_0 \right)G^R(t) = \delta(t).

In energy space,

GR(E)=1E−E0+i0.G^R(E) = \frac{1}{E-E_0+i0}.

The pole prescription determines which inverse Fourier transform is allowed. With the retarded convention, the result vanishes for negative time. The advanced object has the opposite boundary value:

GA(E)=1E−E0−i0.G^A(E) = \frac{1}{E-E_0-i0}.

This example contains the whole idea: the local inverse is fixed by the pole, and the time boundary condition is fixed by how the pole is bypassed.

Retarded and Feynman Green functions are often confused because both are called propagators in some contexts. They answer different questions.

A retarded correlator in many-body or field-theory notation commonly has the form

GABR(t)=−iℏθ(t)⟨[A(t),B(0)]⟩,G^R_{AB}(t) = -\frac{i}{\hbar} \theta(t) \langle[A(t),B(0)]\rangle,

for bosonic operators. It measures causal linear response: a disturbance coupled to BB affects measurements of AA only at later times.

Green Functions in Many-Body QM owns the distinct normal fermionic propagator, whose retarded bracket is an anticommutator and whose intermediate states lie in N±1N\pm1 sectors. That object is not an ordinary Kubo susceptibility.

Retarded and Advanced Response owns the many-body observable commutator pair, its adjoint identities, response spectral density, dispersion relations, and stability diagnostics. The present page retains the inverse-kernel boundary prescription.

A Feynman or time-ordered correlator has a different ordering prescription, schematically

GABF(t)=−iℏ⟨TA(t)B(0)⟩.G^F_{AB}(t) = -\frac{i}{\hbar} \langle\mathcal T A(t)B(0)\rangle.

It is the natural object in vacuum perturbation theory. It is not the same as the retarded response function, and its analytic prescription is not obtained by replacing every denominator with the same +i0+i0 rule.

The QFT bridge page From Propagators in QM to Propagators in QFT explains this distinction from the field-theory side.

  • Dropping the i0i0 prescription after writing E−H±i0E-H\pm i0.
  • Calling GRG^R, GAG^A, and GFG^F interchangeable because they solve related equations.
  • Treating the advanced Green function as a causal response function for ordinary initial-value dynamics.
  • Forgetting that the sign convention depends on the Fourier transform convention and on the normalization of GG.
  • Inferring that a nonzero Feynman propagator at spacelike separation would by itself imply acausal signaling.
  • G. D. Mahan, Many-Particle Physics, 3rd ed., Kluwer Academic/Plenum, 2000.
  • A. L. Fetter and J. D. Walecka, Quantum Theory of Many-Particle Systems, Dover, 2003.
  • J. J. Sakurai and J. Napolitano, Modern Quantum Mechanics, 3rd ed., Cambridge University Press, 2020.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994.
  • M. E. Peskin and D. V. Schroeder, An Introduction to Quantum Field Theory, Addison-Wesley, 1995.
  • M. D. Schwartz, Quantum Field Theory and the Standard Model, Cambridge University Press, 2014.
  1. Verify that GR(t,t′)=−(i/ℏ)θ(t−t′)U(t,t′)G^R(t,t')=-(i/\hbar)\theta(t-t')U(t,t') satisfies the Green-function equation for a time-independent Hamiltonian.
Solution

Differentiate the product of the step function and the evolution operator:

∂GR∂t=−iℏδ(t−t′)U(t,t′)−iℏθ(t−t′)∂U(t,t′)∂t.\frac{\partial G^R}{\partial t} = -\frac{i}{\hbar} \delta(t-t')U(t,t') - \frac{i}{\hbar} \theta(t-t') \frac{\partial U(t,t')}{\partial t}.

The first term gives δ(t−t′)I\delta(t-t')I after multiplication by iℏi\hbar, because U(t′,t′)=IU(t',t')=I. The second term cancels HGRHG^R because iℏ ∂tU=HUi\hbar\,\partial_tU=HU. Therefore

(iℏ∂∂t−H)GR(t,t′)=δ(t−t′)I.\left( i\hbar\frac{\partial}{\partial t} - H \right)G^R(t,t') = \delta(t-t')I.
  1. Show why the retarded Fourier transform gives GR(E)=(E−H+i0)−1G^R(E)=(E-H+i0)^{-1}.
Solution

For t>0t\gt0,

GR(t)=−iℏe−iHt/ℏ.G^R(t) = -\frac{i}{\hbar}e^{-iHt/\hbar}.

Insert a convergence factor e−ϵt/ℏe^{-\epsilon t/\hbar}:

GR(E)=−iℏ∫0∞dt ei(E−H)t/ℏe−ϵt/ℏ.G^R(E) = -\frac{i}{\hbar} \int_0^\infty dt\, e^{i(E-H)t/\hbar} e^{-\epsilon t/\hbar}.

The integral gives

GR(E)=1E−H+iϵ.G^R(E) = \frac{1}{E-H+i\epsilon}.

Taking ϵ→0+\epsilon\to0^+ yields the retarded boundary value.

  1. Use the distribution identity for 1/(x±i0)1/(x\pm i0) to derive GR(E)−GA(E)=−2πi δ(E−H)G^R(E)-G^A(E)=-2\pi i\,\delta(E-H).
Solution

The scalar identities are

1x+i0=PV⁡1x−iπδ(x),\frac{1}{x+i0} = \operatorname{PV}\frac{1}{x} - i\pi\delta(x),

and

1x−i0=PV⁡1x+iπδ(x).\frac{1}{x-i0} = \operatorname{PV}\frac{1}{x} + i\pi\delta(x).

Subtracting gives

1x+i0−1x−i0=−2πi δ(x).\frac{1}{x+i0} - \frac{1}{x-i0} = -2\pi i\,\delta(x).

Applying this through the spectral theorem with x=E−Hx=E-H gives

GR(E)−GA(E)=−2πi δ(E−H).G^R(E)-G^A(E) = -2\pi i\,\delta(E-H).
  1. Explain why a time-ordered Green function is not the same object as a retarded Green function.
Solution

A retarded Green function is supported only after the source time and is tied to causal linear response. A time-ordered Green function orders operators according to time and is tied to vacuum or perturbative boundary prescriptions. In operator language, retarded functions involve commutators or anticommutators with a step function, while Feynman functions involve time-ordered products. They may be related by analytic continuation or spectral formulas in special settings, but they are not interchangeable definitions.