Skip to content

Lippmann–Schwinger Equation

The Lippmann–Schwinger equation rewrites the stationary scattering problem as a free incoming state plus the wave generated by the potential. It is the formal starting point for scattering Green functions and the Born expansion. Scattering States and Boundary Conditions owns the incoming and outgoing prescriptions; S-Matrix owns the asymptotic operator; T-Matrix owns the transition operator; Green Function for Scattering owns the outgoing-kernel and far-field extraction; Born Series owns repeated iteration and convergence; and Lippmann–Schwinger Equation Preview owns the general resolvent map.

Let

H=H0+V,H=H_0+V,

where H0H_0 is the free Hamiltonian and VV is a short-range scattering potential. A scattering state at energy EE satisfies

H∣ψ(±)⟩=E∣ψ(±)⟩.H\lvert\psi^{(\pm)}\rangle = E\lvert\psi^{(\pm)}\rangle.

The Lippmann–Schwinger equation is

∣ψ(±)⟩=∣ϕ⟩+G0(±)(E)V∣ψ(±)⟩,\lvert\psi^{(\pm)}\rangle = \lvert\phi\rangle + G_0^{(\pm)}(E) V \lvert\psi^{(\pm)}\rangle,

where ∣ϕ⟩\lvert\phi\rangle is a free solution at the same energy and

G0(±)(E)=1E−H0±i0G_0^{(\pm)}(E) = \frac{1}{E-H_0\pm i0}

is the free resolvent with outgoing or incoming boundary conditions. Its coordinate kernel and normalization are derived in Energy Green Function.

Start from

(E−H0)∣ψ⟩=V∣ψ⟩.(E-H_0)\lvert\psi\rangle = V\lvert\psi\rangle.

If E−H0E-H_0 were invertible in the ordinary sense, one could multiply by (E−H0)−1(E-H_0)^{-1}. But EE lies in the continuous spectrum of H0H_0, so the inverse needs a boundary-condition prescription. This is why the resolvents are written as

G0(+)(E)=lim⁡ϵ→0+1E−H0+iϵ,G_0^{(+)}(E) = \lim_{\epsilon\to0^+} \frac{1}{E-H_0+i\epsilon},

and

G0(−)(E)=lim⁡ϵ→0+1E−H0−iϵ.G_0^{(-)}(E) = \lim_{\epsilon\to0^+} \frac{1}{E-H_0-i\epsilon}.

The sign selects outgoing or incoming spherical waves at infinity.

The state ∣ψ(+)⟩\lvert\psi^{(+)}\rangle is the physical scattering state built from a specified incoming free state and outgoing scattered radiation:

∣ψ(+)⟩=∣ϕ⟩+G0(+)V∣ψ(+)⟩.\lvert\psi^{(+)}\rangle = \lvert\phi\rangle + G_0^{(+)}V\lvert\psi^{(+)}\rangle.

The state ∣ψ(−)⟩\lvert\psi^{(-)}\rangle has the opposite boundary condition and is useful for formal scattering theory and matrix elements:

∣ψ(−)⟩=∣ϕ⟩+G0(−)V∣ψ(−)⟩.\lvert\psi^{(-)}\rangle = \lvert\phi\rangle + G_0^{(-)}V\lvert\psi^{(-)}\rangle.

The superscripts do not mean positive and negative energy. They label boundary conditions.

For a particle in three dimensions with incoming plane wave eik⋅re^{i\mathbf k\cdot\mathbf r},

E=ℏ2k22m.E = \frac{\hbar^2k^2}{2m}.

The outgoing free Green function gives

ψk(+)(r)=eik⋅r−2mℏ2∫d3r′ eik∣r−r′∣4π∣r−r′∣V(r′)ψk(+)(r′).\psi_{\mathbf k}^{(+)}(\mathbf r) = e^{i\mathbf k\cdot\mathbf r} - \frac{2m}{\hbar^2} \int d^3r'\, \frac{ e^{ik|\mathbf r-\mathbf r'|} }{ 4\pi|\mathbf r-\mathbf r'| } V(\mathbf r') \psi_{\mathbf k}^{(+)}(\mathbf r').

This is an integral equation: the unknown wavefunction appears on both sides.

The sign, source normalization, current direction, and controlled far-field expansion of this kernel are developed in Green Function for Scattering.

At large rr with r′\mathbf r' in the finite range of the potential,

∣r−r′∣≈r−r^⋅r′.|\mathbf r-\mathbf r'| \approx r-\hat{\mathbf r}\cdot\mathbf r'.

Then

eik∣r−r′∣∣r−r′∣≈eikrre−ik′⋅r′,k′=kr^.\frac{ e^{ik|\mathbf r-\mathbf r'|} }{ |\mathbf r-\mathbf r'| } \approx \frac{e^{ikr}}{r} e^{-i\mathbf k'\cdot\mathbf r'}, \qquad \mathbf k'=k\hat{\mathbf r}.

Comparing with

ψk(+)(r)∼eik⋅r+f(θ,ϕ)eikrr,\psi_{\mathbf k}^{(+)}(\mathbf r) \sim e^{i\mathbf k\cdot\mathbf r} + f(\theta,\phi) \frac{e^{ikr}}{r},

one obtains

f(θ,ϕ)=−m2πℏ2∫d3r′ e−ik′⋅r′V(r′)ψk(+)(r′).f(\theta,\phi) = - \frac{m}{2\pi\hbar^2} \int d^3r'\, e^{-i\mathbf k'\cdot\mathbf r'} V(\mathbf r') \psi_{\mathbf k}^{(+)}(\mathbf r').

This formula is exact within the stated potential-scattering assumptions. It becomes a practical approximation when ψk(+)\psi_{\mathbf k}^{(+)} under the integral is approximated.

Repeated substitution formally gives

∣ψ(+)⟩=∣ϕ⟩+G0(+)V∣ϕ⟩+G0(+)VG0(+)V∣ϕ⟩+⋯ .\lvert\psi^{(+)}\rangle = \lvert\phi\rangle + G_0^{(+)}V\lvert\phi\rangle + G_0^{(+)}VG_0^{(+)}V\lvert\phi\rangle + \cdots.

This equation records where the expansion comes from. Born Series owns its interaction-order interpretation, second-order momentum and coordinate forms, convergence conditions, unitarity bookkeeping, and failure modes.

Plane waves are not normalizable Hilbert-space vectors, and scattering states live naturally in a continuum-normalized or rigged-Hilbert-space setting. The i0i0 prescription is also a distributional boundary condition, not an ordinary finite matrix inverse.

For practical quantum mechanics, the equation is still extremely useful: it packages boundary conditions, Green functions, and scattering amplitudes into one compact formalism. For mathematical scattering theory, existence and completeness of wave operators require additional hypotheses on VV.

  • Interpreting the +i0+i0 sign as a small physical absorption rather than an outgoing boundary condition.
  • Forgetting that ∣ψ(+)⟩\lvert\psi^{(+)}\rangle appears on both sides of the equation.
  • Replacing ψ(+)\psi^{(+)} by the incoming plane wave without recognizing that this is the first Born approximation.
  • Applying the short-range asymptotic form without checking long-range Coulomb behavior.
  • Treating plane-wave states as square-normalizable bound states.
  1. Explain why G0(+)G_0^{(+)} and G0(−)G_0^{(-)} are different even though both contain the same free Hamiltonian.
Solution

The difference is the boundary-condition prescription:

G0(±)(E)=1E−H0±i0.G_0^{(\pm)}(E) = \frac{1}{E-H_0\pm i0}.

Because EE lies in the continuous spectrum, the limiting prescription matters. The plus sign selects outgoing-wave behavior, while the minus sign selects incoming-wave behavior.

  1. Starting from the coordinate-space Lippmann–Schwinger equation, identify the term that produces the outgoing spherical wave.
Solution

The integral term contains

eik∣r−r′∣∣r−r′∣.\frac{ e^{ik|\mathbf r-\mathbf r'|} }{ |\mathbf r-\mathbf r'| }.

At large rr this behaves as

eikrre−ik′⋅r′,\frac{e^{ikr}}{r} e^{-i\mathbf k'\cdot\mathbf r'},

which has the outgoing spherical factor eikr/re^{ikr}/r. Its coefficient is the scattering amplitude.

  1. What approximation is made when ψk(+)(r′)\psi_{\mathbf k}^{(+)}(\mathbf r') is replaced by eik⋅r′e^{i\mathbf k\cdot\mathbf r'} inside the amplitude integral?
Solution

That replacement keeps only the incident free wave inside the potential region and neglects rescattering. It is the first Born approximation. It is valid only when the scattered wave generated by the potential is small enough compared with the incident wave in the region that contributes to the integral.

  • J. R. Taylor, Scattering Theory: The Quantum Theory of Nonrelativistic Collisions, Dover, 2006.
  • R. G. Newton, Scattering Theory of Waves and Particles, 2nd ed., Springer, 1982.
  • 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.