Skip to content

Lippmann–Schwinger Equation Preview

The Lippmann–Schwinger equation turns a stationary scattering problem into a free incident state plus a wave emitted by an effective source. For

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

an outgoing scattering state at energy EE is written formally as

∣ψ(+)⟩=∣ϕ⟩+1E−H0+i0V∣ψ(+)⟩.\lvert\psi^{(+)}\rangle = \lvert\phi\rangle + \frac{1}{E-H_0+i0} V\lvert\psi^{(+)}\rangle.

The three terms have distinct roles:

  • ∣ϕ⟩\lvert\phi\rangle is a prescribed free solution at energy EE;
  • V∣ψ(+)⟩V\lvert\psi^{(+)}\rangle acts as the source generated inside the interaction region;
  • (E−H0+i0)−1(E-H_0+i0)^{-1} propagates that source with outgoing boundary conditions.

This page is the canonical operator-level bridge between resolvents and scattering states. The exact coordinate state equation belongs to Lippmann–Schwinger Equation, while repeated iteration and convergence belong to Born Series in the scattering volume.

Split the Hamiltonian into a reference problem and an interaction:

H=H0+V.H=H_0+V.

Usually H0H_0 describes free relative motion and VV is a sufficiently short-range potential. A stationary interacting state obeys

(E−H)∣ψ⟩=0,(E-H)\lvert\psi\rangle=0,

or equivalently

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

This has the form of an inhomogeneous free equation. The quantity V∣ψ⟩V\lvert\psi\rangle is not an externally prescribed source because it contains the unknown state, but the Green-function logic is the same: invert the free operator with a boundary prescription and impose self-consistency.

The incident state satisfies the homogeneous free equation

(E−H0)∣ϕ⟩=0.(E-H_0)\lvert\phi\rangle=0.

It cannot be discarded. At a continuum energy, E−H0E-H_0 is not ordinarily invertible, and the general solution of the inhomogeneous equation equals one particular solution plus a homogeneous solution. In scattering, ∣ϕ⟩\lvert\phi\rangle specifies the beam or channel incident on the interaction region.

The choice of H0H_0 is part of the formulation. For an unscreened Coulomb force, ordinary plane waves and free spherical asymptotics are not the correct long-distance reference; Coulomb Scattering treats that long-range modification.

For complex zz outside the spectra, define the free and full resolvents

R0(z)=(z−H0)−1,R(z)=(z−H)−1.R_0(z)=(z-H_0)^{-1}, \qquad R(z)=(z-H)^{-1}.

Since

z−H=(z−H0)−V,z-H = (z-H_0)-V,

the second resolvent identity gives

R(z)=R0(z)+R0(z)VR(z).R(z) = R_0(z) + R_0(z)VR(z).

An equally valid ordering is

R(z)=R0(z)+R(z)VR0(z).R(z) = R_0(z) + R(z)VR_0(z).

These are exact operator identities wherever the resolvents exist. Iteration produces the formal resolvent series

R=R0+R0VR0+R0VR0VR0+⋯ ,R = R_0 + R_0VR_0 + R_0VR_0VR_0 + \cdots,

but convergence requires more than writing the series. The Resolvent Operator page owns the abstract inverse and analytic identities; the scattering volume owns the approximation theory built from iteration.

At a real continuum energy, take boundary values

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

The stationary equation then becomes the Lippmann–Schwinger equation

∣ψ(±)⟩=∣ϕ⟩+R0(±)(E)V∣ψ(±)⟩.\lvert\psi^{(\pm)}\rangle = \lvert\phi\rangle + R_0^{(\pm)}(E) V\lvert\psi^{(\pm)}\rangle.

Applying E−H0E-H_0 recovers

(E−H0)∣ψ(±)⟩=V∣ψ(±)⟩(E-H_0)\lvert\psi^{(\pm)}\rangle = V\lvert\psi^{(\pm)}\rangle

in the appropriate generalized sense. The homogeneous incident state is precisely what would be lost by treating R0(E)R_0(E) as an ordinary inverse and multiplying naively.

For the time convention e−iEt/ℏe^{-iEt/\hbar}, the free three-dimensional coordinate kernels at

E=ℏ2k22m,k>0,E=\frac{\hbar^2k^2}{2m}, \qquad k\gt0,

are

G0(±)(r,r′;E)=−m2πℏ2e±ik∣r−r′∣∣r−r′∣.G_0^{(\pm)}(\mathbf r,\mathbf r';E) = -\frac{m}{2\pi\hbar^2} \frac{ e^{\pm ik\lvert\mathbf r-\mathbf r'\rvert} }{ \lvert\mathbf r-\mathbf r'\rvert }.

Thus +i0+i0 selects a wave moving outward from a localized source, while −i0-i0 selects the incoming counterpart. The normalization and source equation are derived in Energy Green Function.

The notation is easy to misread:

StateFree componentScattered boundary condition
∣ψ(+)⟩\lvert\psi^{(+)}\rangleprescribed incident stateoutgoing from the interaction region
∣ψ(−)⟩\lvert\psi^{(-)}\rangleprescribed asymptotic stateincoming toward the interaction region

The superscript does not label positive or negative energy. Nor does it mean that the whole ∣ψ(+)⟩\lvert\psi^{(+)}\rangle state is outgoing: it contains an incident free component and an outgoing scattered component.

Energy-domain outgoing behavior and time-domain retardation are related by analytic boundary values, but they answer different formulations of a problem. In the stationary scattering equation, the immediate role of i0i0 is to select the spatial radiation condition at infinity.

Continuum scattering states are generally not normalizable Hilbert-space vectors. They are generalized eigenstates with delta-function normalization, or distributions in a rigged-Hilbert-space description. Accordingly, the limit

R0(E±i0)R_0(E\pm i0)

may exist only between suitable weighted spaces or inside matrix elements. The notation is compact, but it encodes a limiting-absorption statement.

When wave operators exist, one may write

∣ψ(±)⟩=Ω(±)∣ϕ⟩,\lvert\psi^{(\pm)}\rangle = \Omega^{(\pm)}\lvert\phi\rangle,

with

Ω(±)=s-lim⁡t→∓∞eiHt/ℏe−iH0t/ℏ.\Omega^{(\pm)} = \underset{t\to\mp\infty}{\operatorname{s-lim}} e^{iHt/\hbar}e^{-iH_0t/\hbar}.

The ++ state is matched to free data in the distant past and develops outgoing scattered radiation; the −- state is matched in the distant future. Existence, isometry, and asymptotic completeness require hypotheses on VV. They are not automatic consequences of writing the formal limit.

Bound states lie outside the scattering subspace generated by these wave operators. Resonances are also not ordinary normalizable scattering eigenstates; they appear through analytic continuation. Their pole structure is discussed in Bound States and Scattering Poles.

Taking a coordinate matrix element gives

ψ(+)(r)=ϕ(r)+∫d3r′ G0(+)(r,r′;E)×V(r′)ψ(+)(r′).\begin{aligned} \psi^{(+)}(\mathbf r) &= \phi(\mathbf r) + \int d^3r'\, G_0^{(+)}(\mathbf r,\mathbf r';E)\\ &\qquad\times V(\mathbf r')\psi^{(+)}(\mathbf r'). \end{aligned}

This makes the source picture concrete: V(r′)ψ(+)(r′)V(\mathbf r')\psi^{(+)}(\mathbf r') emits an outgoing wave from every point in the interaction region, and the exact state is the self-consistent sum of those contributions plus the incident wave.

The next steps have separate canonical homes:

QuestionCanonical page
How is the large-rr coefficient extracted?Scattering Amplitude
How is the integral equation derived and iterated?Lippmann–Schwinger Equation
When may the exact state under the integral be replaced by the incident state?First Born Approximation
How does the amplitude become an observable rate?Cross Sections
How are rotationally invariant problems organized?Partial-Wave Expansion
  • Treating (E−H0)−1(E-H_0)^{-1} as an ordinary bounded inverse at a continuum energy.
  • Multiplying by a resolvent and accidentally dropping the incident homogeneous solution.
  • Reading the superscripts (±)(\pm) as signs of the energy.
  • Calling ∣ψ(+)⟩\lvert\psi^{(+)}\rangle a purely outgoing state and forgetting its incident component.
  • Replacing +i0+i0 with a finite absorption without stating that the physical problem has changed.
  • Iterating the equation and assuming the resulting Born series converges for every potential and energy.
  • Applying free short-range asymptotics to an unscreened Coulomb potential.
  • Treating plane waves as normalized vectors rather than generalized states.
  • Assuming formal wave operators exist and are complete without conditions on the interaction.
  • Confusing bound-state poles, resonance poles, and continuum boundary values.
  • 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.
  • M. Reed and B. Simon, Methods of Modern Mathematical Physics III: Scattering Theory, Academic Press, 1979.
  • D. R. Yafaev, Mathematical Scattering Theory: General Theory, American Mathematical Society, 1992.
  • J. J. Sakurai and J. Napolitano, Modern Quantum Mechanics, 3rd ed., Cambridge University Press, 2020.
  1. Derive R=R0+R0VRR=R_0+R_0VR from the definitions of the free and full resolvents.
Solution

Use

R−1=z−H0−V,R0−1=z−H0.R^{-1}=z-H_0-V, \qquad R_0^{-1}=z-H_0.

Then

R0−1−R−1=V.R_0^{-1}-R^{-1}=V.

Multiplying on the left by R0R_0 and on the right by RR gives

R−R0=R0VR.R-R_0=R_0VR.

Therefore

R=R0+R0VR.R=R_0+R_0VR.

Reversing the multiplication order gives the companion identity R=R0+RVR0R=R_0+RVR_0.

  1. Verify formally that a solution of the Lippmann–Schwinger equation satisfies the stationary Schrödinger equation.
Solution

Start from

∣ψ(±)⟩=∣ϕ⟩+R0(±)V∣ψ(±)⟩.\lvert\psi^{(\pm)}\rangle = \lvert\phi\rangle + R_0^{(\pm)}V \lvert\psi^{(\pm)}\rangle.

Apply E−H0E-H_0. Since (E−H0)∣ϕ⟩=0(E-H_0)\lvert\phi\rangle=0 and, in the generalized boundary-value sense,

(E−H0)R0(±)=I,(E-H_0)R_0^{(\pm)}=I,

one obtains

(E−H0)∣ψ(±)⟩=V∣ψ(±)⟩.(E-H_0)\lvert\psi^{(\pm)}\rangle = V\lvert\psi^{(\pm)}\rangle.

Moving the right side to the left gives (E−H)∣ψ(±)⟩=0(E-H)\lvert\psi^{(\pm)}\rangle=0.

  1. Why must the incident state be added even though a Green function is called an inverse?
Solution

At a continuum energy, E−H0E-H_0 has nontrivial generalized homogeneous solutions and no ordinary bounded inverse. A Green-function boundary value supplies one particular solution of

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

but the general solution also contains a homogeneous term. The chosen ∣ϕ⟩\lvert\phi\rangle specifies the incoming channel. Omitting it would erase the incident beam and change the boundary-value problem.

  1. For the convention e−iEt/ℏe^{-iEt/\hbar}, identify the spatial behavior selected by R0(E+i0)R_0(E+i0) in three dimensions.
Solution

Its coordinate kernel is proportional to

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

At large distance from a localized source, this behaves as eikr/re^{ikr}/r, an outgoing spherical wave for the stated time convention. Thus E+i0E+i0 selects outgoing scattered radiation. The opposite boundary value produces e−ikr/re^{-ikr}/r and incoming behavior.

  1. What is the first approximation produced by iterating the Lippmann–Schwinger equation once, and why can it fail near a resonance?
Solution

Insert the incident state for the exact state on the right:

∣ψ(+)⟩≈∣ϕ⟩+R0(+)V∣ϕ⟩.\lvert\psi^{(+)}\rangle \approx \lvert\phi\rangle + R_0^{(+)}V\lvert\phi\rangle.

This is the first Born approximation at the level of the state. Near a resonance, repeated propagation and interaction can be enhanced by small energy denominators. The omitted terms R0VR0V∣ϕ⟩R_0VR_0V\lvert\phi\rangle, and higher iterations, need not be small even when a naive potential-strength estimate looks modest.