Skip to content

Scattering Amplitudes and Dyson Resummation

The previous pages taught us to read physics from singularities of Green functions: poles describe stable particles, cuts describe continua, and imaginary parts appear when intermediate states can go on shell. This page turns that analytic information into a scattering amplitude. The cleanest setting is nonrelativistic potential scattering, because it contains the same architecture as relativistic QFT but without spin, antiparticles, or complicated phase space.

The central idea is that a Green function contains two kinds of information. Its external propagators describe how free particles travel to and from the interaction region. Its middle part describes the actual scattering. To get the scattering amplitude, we amputate the external propagators and put the remaining object on shell. The result is the TT-matrix.

The same organization also explains Dyson resummation. Instead of summing every diagram from scratch, one isolates a kernel that cannot be separated by cutting a chosen propagator, and then sums repeated insertions of that kernel. For a one-particle propagator the kernel is the self-energy. For potential scattering the kernel is the potential, or more generally an irreducible scattering kernel.

The key distinction for this page is:

ObjectOff shell?Contains external propagation?Directly observable?
G(ω;p,p)G(\omega;\mathbf p',\mathbf p)yesyesno
Γ(ω;p,p)\Gamma(\omega;\mathbf p',\mathbf p)yesnono
T(p,p)T(\mathbf p',\mathbf p)nonoyes, through SS and cross-sections

The off-shell objects are essential for calculation, but the physical scattering amplitude is the on-shell amputated object.

Free propagation in a potential background

Section titled “Free propagation in a potential background”

Consider a nonrelativistic field moving in a time-independent external potential U(x)U(\mathbf x). A convenient quadratic action is

S=dtd3xψ(it+22mU(x))ψ.S=\int dt\,d^3x\, \psi^\dagger\left(i\partial_t+{\nabla^2\over2m}-U(\mathbf x)\right)\psi.

The equation of motion is the Schrödinger equation

itψ(t,x)=(22m+U(x))ψ(t,x).i\partial_t\psi(t,\mathbf x) =\left(-{\nabla^2\over2m}+U(\mathbf x)\right)\psi(t,\mathbf x).

In the absence of UU, the time-ordered Green function is

G0(tt,p)=dω2πeiω(tt)ωϵp+i0.G_0(t-t',\mathbf p) =\int {d\omega\over2\pi}\, {e^{-i\omega(t-t')}\over \omega-\epsilon_{\mathbf p}+i0}.

The i0i0 prescription says that a particle propagates forward in time. For t>tt>t', closing the contour in the lower half-plane gives

G0(tt,p)=ieiϵp(tt),G_0(t-t',\mathbf p)=-i\,e^{-i\epsilon_{\mathbf p}(t-t')},

while for t<tt<t' the contour encloses no pole. This is the nonrelativistic version of the causal boundary condition already used for the Feynman propagator.

The potential has momentum-space matrix element

U(p,p)=d3xei(pp)xU(x).U(\mathbf p',\mathbf p) =\int d^3x\,e^{-i(\mathbf p'-\mathbf p)\cdot\mathbf x}U(\mathbf x).

Because the potential is time-independent, it conserves energy but not necessarily momentum. Each insertion of UU changes the spatial momentum by the amount carried by the background.

The Born series as a Green-function expansion

Section titled “The Born series as a Green-function expansion”

The exact Green function in the potential background obeys the operator identity

G=G0+G0UG.G=G_0+G_0UG.

Iterating it gives the Born series

G=G0+G0UG0+G0UG0UG0+.\boxed{ G=G_0+G_0UG_0+G_0UG_0UG_0+\cdots. }

In momentum space,

G(ω;p,p)=(2π)3δ(3)(pp)G0(ω,p)+G0(ω,p)U(p,p)G0(ω,p)+d3k(2π)3G0(ω,p)U(p,k)G0(ω,k)U(k,p)G0(ω,p)+.\begin{aligned} G(\omega;\mathbf p',\mathbf p) &=(2\pi)^3\delta^{(3)}(\mathbf p'-\mathbf p)G_0(\omega,\mathbf p)\\ &\quad+G_0(\omega,\mathbf p')U(\mathbf p',\mathbf p)G_0(\omega,\mathbf p)\\ &\quad+\int{d^3k\over(2\pi)^3} G_0(\omega,\mathbf p')U(\mathbf p',\mathbf k) G_0(\omega,\mathbf k)U(\mathbf k,\mathbf p)G_0(\omega,\mathbf p)+\cdots. \end{aligned}

The first term is free propagation. The second term says that the particle propagates freely, scatters once, then propagates freely again. The third term has one intermediate momentum integral and two scatterings. Higher terms repeat the same pattern.

Born series for potential scattering

The Born series for the Green function in a background potential. Thin lines denote G0G_0, crosses denote insertions of UU, and the intermediate momenta are integrated over. The external free propagators will be removed when extracting the scattering amplitude.

This expansion is the nonrelativistic ancestor of perturbation theory in QFT. The only difference is the meaning of the vertex. Here UU is a fixed external potential; in QFT, vertices come from interaction terms in the action and internal lines can represent all allowed virtual particles.

It is inefficient to keep writing the external propagators. Define the amputated scattering kernel Γ\Gamma by

G(ω;p,p)=(2π)3δ(3)(pp)G0(ω,p)+G0(ω,p)Γ(ω;p,p)G0(ω,p).\boxed{ G(\omega;\mathbf p',\mathbf p) =(2\pi)^3\delta^{(3)}(\mathbf p'-\mathbf p)G_0(\omega,\mathbf p) +G_0(\omega,\mathbf p')\, \Gamma(\omega;\mathbf p',\mathbf p)\, G_0(\omega,\mathbf p). }

Then the Born series for Γ\Gamma is

Γ=U+UG0U+UG0UG0U+.\boxed{ \Gamma=U+UG_0U+UG_0UG_0U+\cdots. }

Equivalently, Γ\Gamma satisfies the Lippmann–Schwinger equation

Γ(ω)=U+UG0(ω)Γ(ω).\boxed{ \Gamma(\omega)=U+UG_0(\omega)\Gamma(\omega). }

In momentum variables this means

Γ(ω;p,p)=U(p,p)+d3k(2π)3U(p,k)1ωϵk+i0Γ(ω;k,p).\Gamma(\omega;\mathbf p',\mathbf p) =U(\mathbf p',\mathbf p) +\int{d^3k\over(2\pi)^3} U(\mathbf p',\mathbf k) {1\over \omega-\epsilon_{\mathbf k}+i0} \Gamma(\omega;\mathbf k,\mathbf p).

The equation is not merely a mnemonic for perturbation theory. It is also a resummation: if the integral equation can be solved, it contains infinitely many potential insertions.

Lippmann–Schwinger equation for the amputated scattering kernel

The amputated kernel satisfies Γ=U+UG0Γ\Gamma=U+UG_0\Gamma. This is the potential-scattering version of the Dyson idea: isolate a kernel and sum its repeated insertions connected by free propagation.

The word “amputated” means exactly what it sounds like. The full connected Green function contains external factors G0(ω,p)G_0(\omega,\mathbf p') and G0(ω,p)G_0(\omega,\mathbf p). Removing them leaves Γ\Gamma. This step is algebraic and can be done off shell.

The next step is physical: scattering states live on the mass shell, so the physical TT-matrix is the on-shell value

T(p,p)=Γ(ω;p,p)ω=ϵp=ϵp.\boxed{ T(\mathbf p',\mathbf p)= \Gamma(\omega;\mathbf p',\mathbf p)\big|_{\omega=\epsilon_{\mathbf p}=\epsilon_{\mathbf p'}}. }

This formula is the simplest version of LSZ logic: amputate the external propagators and put the external particles on shell.

The on-shell delta function arises from the long-time limit. Look at the connected part of the two-point Green function,

Gconn(ω;p,p)=G0(ω,p)Γ(ω;p,p)G0(ω,p).G_{\rm conn}(\omega;\mathbf p',\mathbf p) =G_0(\omega,\mathbf p')\Gamma(\omega;\mathbf p',\mathbf p)G_0(\omega,\mathbf p).

Fourier transform it over a long time interval. If Γ(ω)\Gamma(\omega) is smooth near the external poles, the singular part is controlled by

I(T)=dω2πeiωTΓ(ω;p,p)(ωϵp+i0)(ωϵp+i0).I(T)=\int {d\omega\over2\pi}\, {e^{-i\omega T}\Gamma(\omega;\mathbf p',\mathbf p) \over(\omega-\epsilon_{\mathbf p'}+i0)(\omega-\epsilon_{\mathbf p}+i0)}.

For T>0T>0, the poles lie just below the real axis. Evaluating their residues gives a difference of two oscillatory exponentials divided by the energy difference:

I(T)ieiϵpTΓ(ϵp;p,p)eiϵpTΓ(ϵp;p,p)ϵpϵp.I(T)\simeq -i\,{e^{-i\epsilon_{\mathbf p'}T}\Gamma(\epsilon_{\mathbf p'};\mathbf p',\mathbf p) -e^{-i\epsilon_{\mathbf p}T}\Gamma(\epsilon_{\mathbf p};\mathbf p',\mathbf p) \over \epsilon_{\mathbf p'}-\epsilon_{\mathbf p}}.

Because G(ω)=(ωH+i0)1G(\omega)=(\omega-H+i0)^{-1} has no overall factor of ii in our convention, its Fourier transform is i-i times the evolution kernel. Thus the connected Schrödinger-picture evolution amplitude is iI(T)iI(T). To obtain the interaction-picture SS-matrix between symmetric endpoints T/2-T/2 and T/2T/2, strip the free phases from both external states:

Sppconn(T)=e+iϵpT/2[iI(T)]e+iϵpT/2i2sin[(ϵpϵp)T/2]ϵpϵpΓ(E;p,p),\begin{aligned} S^{\rm conn}_{\mathbf p'\mathbf p}(T) &=e^{+i\epsilon_{\mathbf p'}T/2}\,[iI(T)]\, e^{+i\epsilon_{\mathbf p}T/2}\\ &\simeq -i\,{2\sin[(\epsilon_{\mathbf p'}-\epsilon_{\mathbf p})T/2] \over \epsilon_{\mathbf p'}-\epsilon_{\mathbf p}}\, \Gamma(E;\mathbf p',\mathbf p), \end{aligned}

where the smooth kernel may be evaluated at the common energy EE in the distributional on-shell limit. Now

limTT/2T/2dtei(ϵpϵp)t=2πδ(ϵpϵp).\lim_{T\to\infty}\int_{-T/2}^{T/2}dt\, e^{i(\epsilon_{\mathbf p'}-\epsilon_{\mathbf p})t} =2\pi\delta(\epsilon_{\mathbf p'}-\epsilon_{\mathbf p}).

and therefore the transition part of the SS-matrix has the form

Sppconn=2πiδ(ϵpϵp)T(p,p),S_{\mathbf p'\mathbf p}^{\rm conn} =-2\pi i\, \delta(\epsilon_{\mathbf p'}-\epsilon_{\mathbf p}) \,T(\mathbf p',\mathbf p),

where TT is the on-shell amputated kernel. The explicit phase stripping is essential: the resolvent transform by itself is a Schrödinger-picture propagation amplitude, not yet the SS-matrix. After that step, the external propagator poles select the on-shell value of Γ\Gamma and supply the energy-conserving delta function.

Amputation of external propagators and on-shell projection

The connected Green function factorizes as G0ΓG0G_0\Gamma G_0. The scattering amplitude is obtained by removing the external propagators and evaluating the remaining kernel at ω=ϵp=ϵp\omega=\epsilon_{\mathbf p}=\epsilon_{\mathbf p'}.

Discontinuities and on-shell intermediate states

Section titled “Discontinuities and on-shell intermediate states”

Suppose a contribution is built from two subkernels AA and BB joined by one free intermediate propagator. Before any on-shell limit, their composition is

C(E;p,p)=d3k(2π)3A(E;p,k)1Eϵk+i0B(E;k,p).C(E;\mathbf p',\mathbf p) =\int {d^3k\over(2\pi)^3} A(E;\mathbf p',\mathbf k) {1\over E-\epsilon_{\mathbf k}+i0} B(E;\mathbf k,\mathbf p).

The intermediate momentum is off shell in this ordinary product. It is placed on shell only when we take the discontinuity of the propagator. If AA and BB have no discontinuity in the channel under discussion, then

DiscC(E;p,p)=2πid3k(2π)3A(E;p,k)δ(Eϵk)B(E;k,p).\boxed{ \operatorname{Disc}C(E;\mathbf p',\mathbf p) =-2\pi i\int {d^3k\over(2\pi)^3} A(E;\mathbf p',\mathbf k) \delta(E-\epsilon_{\mathbf k}) B(E;\mathbf k,\mathbf p). }

This is the one-particle version of a much broader principle: discontinuities factorize through the phase space of intermediate on-shell states. In relativistic QFT the one-particle measure becomes

d3k(2π)32Ek,{d^3k\over(2\pi)^3 2E_{\mathbf k}},

and the delta functions impose full four-momentum conservation. The optical theorem and cutting rules are sophisticated versions of the same statement.

Off-shell kernel composition and the on-shell discontinuity of its intermediate propagator

The ordinary composition C=AG0BC=AG_0B contains an off-shell propagator. Its discontinuity replaces that propagator by 2πiδ(Eϵk)-2\pi i\delta(E-\epsilon_{\mathbf k}), exposing the phase space of a real intermediate particle.

This is also where reducible and irreducible diagrams become useful language. A diagram is reducible with respect to a chosen propagator if cutting that propagator separates the diagram into two pieces. The repeated terms in the Born series are reducible in this sense. The potential UU is the irreducible kernel for the simple problem above. In an interacting QFT, the corresponding kernel is often a sum of many diagrams that are irreducible with respect to the relevant cut.

The self-energy resummation from earlier pages had the form

G=G0+G0ΣG0+G0ΣG0ΣG0+.G=G_0+G_0\Sigma G_0+G_0\Sigma G_0\Sigma G_0+\cdots.

This geometric series gives

G=G0+G0ΣG,G1=G01Σ.G=G_0+G_0\Sigma G, \qquad G^{-1}=G_0^{-1}-\Sigma.

Potential scattering has the parallel structure

Γ=U+UG0Γ,Γ=(1UG0)1U,\Gamma=U+UG_0\Gamma, \qquad \Gamma=(1-UG_0)^{-1}U,

and therefore

T=Γon shell.T=\Gamma\big|_{\rm on\ shell}.

The exact same logic appears later in many forms:

problemirreducible kernelresummed objectone-particle propagationΣGpotential scatteringUΓ off shell; T=Γon shelltwo-particle scatteringKfour-point function; M on shell\begin{array}{c|c|c} \text{problem} & \text{irreducible kernel} & \text{resummed object} \\ \hline \text{one-particle propagation} & \Sigma & G \\ \text{potential scattering} & U & \Gamma\ \text{off shell};\ T=\Gamma|_{\rm on\ shell} \\ \text{two-particle scattering} & K & \text{four-point function};\ \mathcal M\ \text{on shell} \end{array}

The common lesson is that the object one should compute is rarely the full answer diagram by diagram. It is usually better to compute an irreducible building block and then let an integral equation perform the resummation.

At leading order,

T(p,p)=U(p,p)+O(U2).T(\mathbf p',\mathbf p)=U(\mathbf p',\mathbf p)+O(U^2).

For elastic scattering, p=p=p|\mathbf p'|=|\mathbf p|=p, and the momentum transfer is

q=pp,q=2psinθ2.\mathbf q=\mathbf p'-\mathbf p, \qquad q=2p\sin{\theta\over2}.

With the normalization used here, the nonrelativistic scattering amplitude f(θ)f(\theta) is related to the TT-matrix by

f(θ)=m2πT(p,p),f(\theta)=-{m\over2\pi}T(\mathbf p',\mathbf p),

so

dσdΩ=f(θ)2=m24π2T(p,p)2.{d\sigma\over d\Omega}=|f(\theta)|^2 ={m^2\over4\pi^2}|T(\mathbf p',\mathbf p)|^2.

For a weak Gaussian potential

U(x)=U0ex2/(2a2),U(\mathbf x)=U_0 e^{-\mathbf x^2/(2a^2)},

we have

U(p,p)=U0(2π)3/2a3ea2(pp)2/2.U(\mathbf p',\mathbf p) =U_0(2\pi)^{3/2}a^3 e^{-a^2(\mathbf p'-\mathbf p)^2/2}.

Therefore the first Born cross-section is

dσdΩ=m24π2U0(2π)3/2a3e2a2p2sin2(θ/2)2.{d\sigma\over d\Omega} ={m^2\over4\pi^2} \left|U_0(2\pi)^{3/2}a^3 e^{-2a^2p^2\sin^2(\theta/2)}\right|^2.

Forward scattering is largest because the Fourier transform of a smooth potential is largest at small momentum transfer. A short-range potential gives a broad angular distribution; a long-range smooth potential strongly favors small angles.

The scattering amplitude is the on-shell, amputated part of a Green function. In the nonrelativistic potential problem,

Gconn=G0ΓG0,T(p,p)=Γ(ϵp;p,p)withϵp=ϵp.G_{\rm conn}=G_0\Gamma G_0, \qquad T(\mathbf p',\mathbf p)=\Gamma(\epsilon_{\mathbf p};\mathbf p',\mathbf p) \quad \text{with}\quad \epsilon_{\mathbf p'}=\epsilon_{\mathbf p}.

The external propagator poles provide the energy-conserving delta function in the SS-matrix, while the amputated kernel contains the real scattering dynamics.

The Born series is a perturbative expansion in the potential:

Γ=U+UG0U+UG0UG0U+.\Gamma=U+UG_0U+UG_0UG_0U+\cdots.

Written as

Γ=U+UG0Γ,\Gamma=U+UG_0\Gamma,

it becomes the Lippmann–Schwinger equation. This is the scattering version of Dyson resummation: isolate an irreducible kernel, then sum its repeated insertions. The same strategy organizes self-energies, scattering amplitudes, Bethe–Salpeter equations, and eventually the effective action.

Confusing the Green function with the scattering amplitude. The full Green function includes external propagation before and after the collision. The scattering amplitude is the amputated, on-shell kernel.

Putting internal lines on shell too early. The internal line in the Lippmann–Schwinger equation is G0(ω,k)=1/(ωϵk+i0)G_0(\omega,\mathbf k)=1/(\omega-\epsilon_{\mathbf k}+i0) and is integrated over all k\mathbf k. An on-shell delta function appears only when taking a discontinuity or an appropriate long-time/asymptotic projection.

Losing the i0i0 prescription. The sign of i0i0 decides whether the solution is outgoing or incoming. Scattering theory is not just algebraic inversion; it is inversion with a boundary condition.

Mixing normalization conventions. The relation between TT and dσ/dΩd\sigma/d\Omega changes if one uses box-normalized states or relativistic normalization. The invariant content is the SS-matrix element with its stated state normalization.

Calling every resummation “Dyson” without naming the kernel. A resummation becomes meaningful only after specifying what is irreducible. For the one-particle propagator it is the 1PI self-energy Σ\Sigma; for potential scattering it is UU or a more general irreducible scattering kernel KK.

Forgetting that off-shell kernels are not observables. Γ(ω;p,p)\Gamma(\omega;\mathbf p',\mathbf p) depends on how the calculation is organized. The on-shell TT-matrix, with the stated normalization, is what enters scattering probabilities.

Starting from the Born series

Γ=U+UG0U+UG0UG0U+,\Gamma=U+UG_0U+UG_0UG_0U+\cdots,

show that Γ\Gamma obeys both

Γ=U+UG0Γ\Gamma=U+UG_0\Gamma

and

Γ=U+ΓG0U.\Gamma=U+\Gamma G_0U.
Solution

Factor the Born series after the first UU:

Γ=U+UG0(U+UG0U+UG0UG0U+).\Gamma =U+UG_0\left(U+UG_0U+UG_0UG_0U+\cdots\right).

The expression in parentheses is Γ\Gamma, so

Γ=U+UG0Γ.\Gamma=U+UG_0\Gamma.

Alternatively, factor the last UU:

Γ=U+(U+UG0U+UG0UG0U+)G0U.\Gamma =U+\left(U+UG_0U+UG_0UG_0U+\cdots\right)G_0U.

Again the parenthesis is Γ\Gamma, giving

Γ=U+ΓG0U.\Gamma=U+\Gamma G_0U.

The two equations are equivalent for the same Born series. Written in momentum space, they differ in whether the potential next to the external final leg or the external initial leg is singled out.

Use the distribution identity

1x+i0=PV1xiπδ(x){1\over x+i0}=\operatorname{PV}{1\over x}-i\pi\delta(x)

to find the imaginary part of the second Born term

T(2)(E;p,p)=d3k(2π)3U(p,k)U(k,p)Eϵk+i0.T^{(2)}(E;\mathbf p',\mathbf p) =\int {d^3k\over(2\pi)^3} {U(\mathbf p',\mathbf k)U(\mathbf k,\mathbf p) \over E-\epsilon_{\mathbf k}+i0}.

Specialize to forward scattering, p=p\mathbf p'=\mathbf p, with E=ϵpE=\epsilon_{\mathbf p}, and take UU to be Hermitian.

Solution

Apply the identity with x=Eϵkx=E-\epsilon_{\mathbf k}:

T(2)(E;p,p)=d3k(2π)3U(k,p)2PV1Eϵkiπd3k(2π)3U(k,p)2δ(Eϵk).T^{(2)}(E;\mathbf p,\mathbf p) =\int {d^3k\over(2\pi)^3} |U(\mathbf k,\mathbf p)|^2 \operatorname{PV}{1\over E-\epsilon_{\mathbf k}} -i\pi\int {d^3k\over(2\pi)^3} |U(\mathbf k,\mathbf p)|^2 \delta(E-\epsilon_{\mathbf k}).

Therefore

ImT(2)(E;p,p)=πd3k(2π)3U(k,p)2δ(Eϵk).\operatorname{Im}T^{(2)}(E;\mathbf p,\mathbf p) =-\pi\int {d^3k\over(2\pi)^3} |U(\mathbf k,\mathbf p)|^2 \delta(E-\epsilon_{\mathbf k}).

The imaginary part is produced only by intermediate momenta satisfying the on-shell condition ϵk=E\epsilon_{\mathbf k}=E. Its negative sign is consistent with this page’s convention S=12πiδ(EfEi)TS=1-2\pi i\delta(E_f-E_i)T (and with f=mT/(2π)f=-mT/(2\pi)). This is the second-order, nonrelativistic shadow of the forward optical theorem and the cutting rule.

For the Gaussian potential

U(x)=U0ex2/(2a2),U(\mathbf x)=U_0e^{-\mathbf x^2/(2a^2)},

compute the first Born approximation to dσ/dΩd\sigma/d\Omega using

f(θ)=m2πT(p,p),T(1)(p,p)=U(p,p).f(\theta)=-{m\over2\pi}T(\mathbf p',\mathbf p), \qquad T^{(1)}(\mathbf p',\mathbf p)=U(\mathbf p',\mathbf p).
Solution

The Fourier transform is

U(p,p)=d3xei(pp)xU0ex2/(2a2).U(\mathbf p',\mathbf p) =\int d^3x\,e^{-i(\mathbf p'-\mathbf p)\cdot\mathbf x} U_0e^{-\mathbf x^2/(2a^2)}.

Using the standard Gaussian integral gives

U(p,p)=U0(2π)3/2a3ea2q2/2,q=pp.U(\mathbf p',\mathbf p) =U_0(2\pi)^{3/2}a^3e^{-a^2q^2/2}, \qquad \mathbf q=\mathbf p'-\mathbf p.

For elastic scattering, q=2psin(θ/2)q=2p\sin(\theta/2). Hence

f(θ)=m2πU0(2π)3/2a3exp[2a2p2sin2θ2].f(\theta)=-{m\over2\pi}U_0(2\pi)^{3/2}a^3 \exp\left[-2a^2p^2\sin^2{\theta\over2}\right].

Therefore

dσdΩ=m24π2U02(2π)3a6exp[4a2p2sin2θ2].{d\sigma\over d\Omega} ={m^2\over4\pi^2} |U_0|^2(2\pi)^3a^6 \exp\left[-4a^2p^2\sin^2{\theta\over2}\right].

The angular distribution narrows as aa grows, because a wider potential has a narrower Fourier transform.

Consider a separable potential

U(p,p)=gu(p)u(p),U(\mathbf p',\mathbf p)=g\,u(\mathbf p')u^*(\mathbf p),

where gg is real. Solve the Lippmann–Schwinger equation for Γ(E;p,p)\Gamma(E;\mathbf p',\mathbf p).

Solution

Use the ansatz

Γ(E;p,p)=τ(E)u(p)u(p).\Gamma(E;\mathbf p',\mathbf p)=\tau(E)u(\mathbf p')u^*(\mathbf p).

Substitute into

Γ=U+UG0Γ.\Gamma=U+UG_0\Gamma.

The right-hand side is

gu(p)u(p)+d3k(2π)3gu(p)u(k)1Eϵk+i0τ(E)u(k)u(p).g u(\mathbf p')u^*(\mathbf p) +\int{d^3k\over(2\pi)^3} g u(\mathbf p')u^*(\mathbf k) {1\over E-\epsilon_{\mathbf k}+i0} \tau(E)u(\mathbf k)u^*(\mathbf p).

Thus

τ(E)=g+gI(E)τ(E),\tau(E)=g+g I(E)\tau(E),

where

I(E)=d3k(2π)3u(k)2Eϵk+i0.I(E)=\int{d^3k\over(2\pi)^3} {|u(\mathbf k)|^2\over E-\epsilon_{\mathbf k}+i0}.

Solving,

τ(E)=g1gI(E)\boxed{ \tau(E)={g\over1-gI(E)} }

and therefore

Γ(E;p,p)=gu(p)u(p)1gI(E).\boxed{ \Gamma(E;\mathbf p',\mathbf p)= {g\,u(\mathbf p')u^*(\mathbf p) \over1-gI(E)}. }

A pole occurs when 1gI(E)=01-gI(E)=0. Depending on its location, this pole describes a bound state, a resonance, or an instability of the perturbative expansion.

  • B. A. Lippmann and J. Schwinger, “Variational Principles for Scattering Processes,” Physical Review 79 (1950), for the original integral-equation formulation of scattering states.
  • Steven Weinberg, The Quantum Theory of Fields, vol. I, sections 3.1–3.6 and 10.3, for scattering states, the SS-matrix, perturbation theory, and pole extraction.
  • Sidney Coleman, Lectures on Quantum Field Theory, chapters 7, 10–14, for perturbation theory, scattering, phase space, and the LSZ viewpoint.
  • Mark Srednicki, Quantum Field Theory, sections 5, 10, and 13–14, for LSZ reduction, scattering amplitudes, and exact propagator structure.
  • M. L. Goldberger and K. M. Watson, Collision Theory, for a classic systematic treatment of potential scattering, the TT-matrix, and the Lippmann–Schwinger equation.