Skip to content

Green Functions and Resolvents

A Green function is an inverse kernel together with a boundary prescription. In quantum mechanics, Green functions solve sourced equations, encode spectral information, select incoming or outgoing behavior, and organize linear response. The same language later expands into many-body correlators and relativistic field propagators.

For a time-independent self-adjoint Hamiltonian, the central operator is the resolvent

R(z)=(zI−H)−1,R(z) = (zI-H)^{-1},

defined when zz lies in the resolvent set of HH. Its coordinate kernel is

G(x,x′;z)=⟨x∣R(z)∣x′⟩,G(x,x';z) = \langle x\rvert R(z)\lvert x'\rangle,

and it satisfies

(z−Hx)G(x,x′;z)=δ(x−x′).(z-H_x)G(x,x';z) = \delta(x-x').

This equation alone does not determine a unique physical Green function when zz approaches the spectrum. The choice of domain, asymptotic behavior, temporal support, and limiting prescription is part of the definition.

This chapter is the canonical home for

  • Green functions as inverse operators and source-response kernels;
  • the resolvent set, spectrum, and operator-valued resolvent;
  • coordinate-space energy Green functions;
  • incoming, outgoing, retarded, and advanced boundary prescriptions;
  • discrete, continuous, and mixed spectral representations;
  • poles, residues, branch cuts, and spectral densities;
  • total and local density of states from resolvent boundary values;
  • introductory bridges to linear response and the Kubo formula;
  • the resolvent structure behind the Lippmann–Schwinger equation;
  • the translation from one-particle quantum mechanics to many-body and QFT Green functions.

The chapter does not own the time-evolution kernel itself; that belongs to Propagators and Kernels. It previews response without replacing the canonical open-systems and response treatment. It previews scattering-state construction without replacing the full Lippmann–Schwinger Equation. It also introduces QFT terminology without importing field-theory locality, vacuum structure, or renormalization into a one-particle setting.

For a linear operator LL, a Green kernel formally satisfies

LxG(x,x′)=δ(x−x′).L_xG(x,x') = \delta(x-x').

If G1G_1 and G2G_2 solve the same source equation, their difference obeys the homogeneous equation:

Lx[G1(x,x′)−G2(x,x′)]=0.L_x \left[ G_1(x,x')-G_2(x,x') \right] = 0.

Boundary or support conditions select among these possibilities. For the Schrödinger operator, common choices include

  • square-integrable behavior for bound-state problems;
  • outgoing or incoming waves in stationary scattering;
  • retarded support after a source acts;
  • advanced support before a designated source time;
  • time ordering for Feynman correlators;
  • periodicity in imaginary time for thermal problems.

Calling two kernels “solutions of the same equation” does not make them the same Green function. Their global prescriptions carry physical information.

The resolvent set ρ(H)\rho(H) consists of complex zz for which zI−HzI-H has a bounded inverse. The spectrum is its complement:

σ(H)=C∖ρ(H).\sigma(H) = \mathbb C\setminus\rho(H).

For self-adjoint HH, every nonreal zz lies in the resolvent set. The spectral theorem gives

R(z)=∫R1z−λ dEH(λ),R(z) = \int_{\mathbb R} \frac{1}{z-\lambda} \,dE_H(\lambda),

where EHE_H is the projection-valued spectral measure. In a finite discrete system,

R(z)=∑nPnz−En.R(z) = \sum_n \frac{P_n}{z-E_n}.

An isolated eigenvalue produces a pole whose residue is the projector onto the corresponding eigenspace. Continuous spectrum appears through boundary values and discontinuities across the real axis rather than isolated normalizable-state poles.

The norm estimate

∥R(z)∥≤1∣Im⁡z∣\lVert R(z)\rVert \le \frac{1}{\lvert\operatorname{Im}z\rvert}

for nonreal zz shows why the resolvent becomes sensitive near the real spectrum. A finite imaginary part is both a regularization and, in numerical work, an energy broadening that must be reported.

At a real energy in the continuous spectrum, use limiting boundary values

R±(E)=lim⁡η→0+(E−H±iη)−1.R^\pm(E) = \lim_{\eta\to0^+} (E-H\pm i\eta)^{-1}.

For standard scattering conventions,

R+(E)selects outgoing behavior,R^+(E) \quad \text{selects outgoing behavior},

while

R−(E)selects incoming behavior.R^-(E) \quad \text{selects incoming behavior}.

The operator-valued distribution identity is

R+(E)−R−(E)=−2πi δ(E−H).R^+(E)-R^-(E) = -2\pi i\,\delta(E-H).

Consequently,

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

These identities connect analytic boundary values to spectral projectors and density of states. They are not ordinary pointwise operator equalities on the spectrum; interpret them through matrix elements, traces, or test-function smearing.

For a time-independent Hamiltonian, one common retarded normalization is

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

It satisfies

(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

and vanishes for t<t′t\lt t'. The advanced solution is

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

With a consistent Fourier convention, the retarded and advanced transforms yield the upper- and lower-half-plane boundary values of the resolvent. Factors of ii, ℏ\hbar, and 2π2\pi vary across conventions, so every comparison should begin by writing the defining source equation and transform pair.

For a finite system, the exact density of states is

ρ(E)=∑nδ(E−En)=Tr⁡δ(E−H).\rho(E) = \sum_n\delta(E-E_n) = \operatorname{Tr}\delta(E-H).

The retarded resolvent gives

ρ(E)=−1πIm⁡Tr⁡R+(E).\rho(E) = -\frac{1}{\pi} \operatorname{Im} \operatorname{Tr}R^+(E).

In a position basis, the local density of states is

ρ(x,E)=−1πIm⁡G+(x,x;E).\rho(x,E) = -\frac{1}{\pi} \operatorname{Im} G^+(x,x;E).

A finite broadening η\eta replaces exact delta peaks by broadened lines. That can model resolution or aid visualization, but it is not the exact finite-system spectrum. State whether ρ\rho is total, local, per volume, per cell, or resolved by spin or another internal label.

QuestionCanonical pageMain object
What does “Green function” mean across quantum mechanics?What Is a Green Function?inverse kernel plus prescription
Where is the inverse defined as an operator?Resolvent OperatorR(z)=(zI−H)−1R(z)=(zI-H)^{-1}
How does the resolvent become a coordinate-space source solution?Energy Green FunctionG(x,x′;z)G(x,x';z)
How do support and i0i0 signs select causal response?Retarded and Advanced Green FunctionsGRG^R and GAG^A
How do poles, continua, and spectral densities enter?Spectral Representation of Green Functionsspectral measures and discontinuities
How is state counting recovered from a resolvent?Green Functions and Density of States−π−1Im⁡Tr⁡R+-\pi^{-1}\operatorname{Im}\operatorname{Tr}R^+
How do retarded commutators encode weak-source response?Green Functions and Response Previewsusceptibility and Kubo bridge
How do free resolvents generate scattering states?Lippmann–Schwinger Equation Preview∣ψ±⟩\lvert\psi^\pm\rangle and R0±VR_0^\pm V
What changes for many-body systems and quantum fields?Green Functions from QM to QFTcorrelators, ordering, and fields

These nine articles form the planned chapter. The response and scattering pages are deliberately labeled previews because their detailed canonical homes require additional physical assumptions and machinery.

Read What Is a Green Function?, Resolvent Operator, and Energy Green Function in order. Then use Retarded and Advanced Green Functions to see how a local inverse equation acquires physical meaning through support and boundary values.

After the resolvent operator, read Spectral Representation and Green Functions and Density of States. Pair them with the Spectral Theorem to distinguish isolated eigenvalue poles, continuous spectral measure, and finite broadening.

Read the energy Green function and its outgoing free-particle examples, then continue to the Lippmann–Schwinger Preview. The full scattering chapter owns existence assumptions, asymptotic states, the TT matrix, cross sections, and practical solution methods.

Read retarded and advanced functions before Green Functions and Response Preview. Retarded support is necessary for causal linear response, but the response function also contains a commutator and a reference state. A one-particle evolution kernel is not automatically a susceptibility.

Finish with Green Functions from QM to QFT and Correlation Functions in Path Integrals. Keep operator ordering, vacuum or thermal state, field content, and normalization explicit. The QFT use of “propagator” is related to, but not identical with, the kernel of one-particle time evolution.

Before using a Green function, state

  1. the operator being inverted and its domain;
  2. the source equation and delta-function normalization;
  3. the coordinate measure or basis normalization;
  4. the complex spectral parameter and sign convention;
  5. the boundary, asymptotic, support, or ordering prescription;
  6. the Fourier or Laplace transform convention;
  7. whether the spectrum is discrete, continuous, or mixed;
  8. whether traces and coincident-point limits require regularization;
  9. whether a finite η\eta is numerical broadening or a limiting prescription;
  10. whether the object is a one-particle inverse, a response correlator, or a QFT propagator.
ObjectDefining structureMain warning
resolvent(zI−H)−1(zI-H)^{-1} off the spectrumsingular boundary values need a prescription
energy Green kernel⟨x∣R(z)∣x′⟩\langle x\rvert R(z)\lvert x'\rangledepends on measure, domain, and asymptotics
retarded Green operatorfuture-supported inversetemporal causality does not imply relativistic light-cone support
advanced Green operatorpast-supported inverseuseful analytically, not the usual initial-value response
propagator kernelmatrix element of U(t,t0)U(t,t_0)evolves states and is not itself a resolvent
response functionretarded commutator in a specified staterequires observables and an ensemble or state
QFT Feynman propagatortime-ordered field correlatordepends on field theory, state, and ordering conventions
  • Writing “the Green function” without naming the operator and prescription. The inverse equation is incomplete by itself.
  • Using (E−H)−1(E-H)^{-1} directly on the spectrum. Boundary values or a different generalized interpretation are required.
  • Dropping the sign of i0i0. Incoming and outgoing, retarded and advanced solutions are then conflated.
  • Treating a finite η\eta as harmless. It broadens spectral features and may hide nearby poles or gaps.
  • Reading every singularity as a normalizable bound state. Continua produce cuts and boundary-value discontinuities.
  • Comparing formulas before reconciling Fourier conventions. Factors of ii, ℏ\hbar, and 2π2\pi can differ legitimately.
  • Taking a trace or coincident-point limit without checking convergence. Ultraviolet or volume divergences may require regularization or normalization per volume.
  • Equating a one-particle retarded inverse with a many-body response function. The latter contains an operator commutator and a chosen state.
  • Using the Lippmann–Schwinger equation without its boundary label. The ++ and −- states encode different asymptotic conditions.
  • Importing relativistic causal claims into nonrelativistic kernels. The theories organize locality differently.
  • E. N. Economou, Green’s Functions in Quantum Physics, 3rd ed., Springer, 2006, doi:10.1007/3-540-28841-4.
  • G. D. Mahan, Many-Particle Physics, 3rd ed., Springer, 2000, doi:10.1007/978-1-4757-5714-9.
  • M. Reed and B. Simon, Methods of Modern Mathematical Physics, Vol. I: Functional Analysis, revised ed., Academic Press, 1980.
  • R. Kubo, “Statistical-Mechanical Theory of Irreversible Processes. I,” Journal of the Physical Society of Japan 12, 570–586 (1957), doi:10.1143/JPSJ.12.570.
  • J. J. Sakurai and J. Napolitano, Modern Quantum Mechanics, 3rd ed., Cambridge University Press, 2020.