Skip to content

Pole Residue and Wavefunction Renormalization

An isolated real pole in an exact two-point function tells us that the chosen field couples to a stable particle. Its residue tells us how strongly the field overlaps with that particle. This is more specific than the pole language of potential scattering on the previous page: poles of a TT-matrix can also describe bound states or, after continuation to another sheet, unstable resonances. LSZ uses the stable one-particle pole of a two-point function.

That distinction is easy to miss. In a free scalar theory the field ϕ\phi is normalized so that its propagator has residue one at p2=m2p^2=m^2. In an interacting theory, the physical one-particle state is no longer the bare oscillator quantum. It is a dressed object: the field can create it, but it can also create multiparticle states with the same quantum numbers. The exact two-point function therefore contains an isolated one-particle pole plus a continuum. The coefficient of the isolated pole is the wavefunction renormalization factor ZZ.

This page develops the formula

G~(p)iZp2M2+i0near the stable one-particle pole,\widetilde G(p) \simeq {iZ\over p^2-M^2+i0} \qquad \text{near the stable one-particle pole},

and explains why

Z=11Π(M2)Z={1\over 1-\Pi'(M^2)}

when the exact propagator is written in terms of a Lorentz-invariant self-energy Π(p2)\Pi(p^2). The same ZZ is the factor that appears on external legs in LSZ reduction. The next page turns this into the full scattering formula.

Three related quantities will appear, and keeping them separate prevents many normalization errors:

QuantityMeaning
iZiZresidue of the covariant propagator as a function of p2p^2
iZ/(2Ep)iZ/(2E_{\mathbf p})residue of the positive-energy pole as a function of p0p^0 at fixed p\mathbf p
Z\sqrt Zmatrix element $\langle0

They are the same physical overlap viewed in three different variables.

For a free real scalar field,

G~0(p)=ip2m2+i0.\widetilde G_0(p)={i\over p^2-m^2+i0}.

At fixed spatial momentum p\mathbf p, define

ωp=p2+m2.\omega_{\mathbf p}=\sqrt{\mathbf p^2+m^2}.

Then

p2m2=(p0)2ωp2,p^2-m^2=(p^0)^2-\omega_{\mathbf p}^2,

and the propagator decomposes as

G~0(p0,p)=i(p0)2ωp2+i0=i2ωp(1p0ωp+i01p0+ωpi0).\widetilde G_0(p^0,\mathbf p) ={i\over (p^0)^2-\omega_{\mathbf p}^2+i0} ={i\over 2\omega_{\mathbf p}} \left( {1\over p^0-\omega_{\mathbf p}+i0} - {1\over p^0+\omega_{\mathbf p}-i0} \right).

The positive-energy pole is therefore at

p0=ωpi0,p^0=\omega_{\mathbf p}-i0,

with residue

Resp0=ωpG~0(p0,p)=i2ωp.\operatorname*{Res}_{p^0=\omega_{\mathbf p}}\widetilde G_0(p^0,\mathbf p) ={i\over 2\omega_{\mathbf p}}.

Equivalently, in the Lorentz-invariant variable s=p2s=p^2, the pole has residue ii:

G~0(p)=ism2+i0.\widetilde G_0(p)={i\over s-m^2+i0}.

Both statements are correct. The p0p^0-residue contains an additional Jacobian because

sm2=(p0)2ωp22ωp(p0ωp)s-m^2=(p^0)^2-\omega_{\mathbf p}^2 \simeq 2\omega_{\mathbf p}(p^0-\omega_{\mathbf p})

near the positive-energy pole.

The time-domain propagator makes the physical meaning especially clear. For t>0t>0,

G0(t,p)=dp02πeip0ti(p0)2ωp2+i0=12ωpeiωpt.G_0(t,\mathbf p) =\int{dp^0\over 2\pi}\, e^{-ip^0t}\,{i\over (p^0)^2-\omega_{\mathbf p}^2+i0} ={1\over 2\omega_{\mathbf p}}e^{-i\omega_{\mathbf p}t}.

The pole gives a long-lived oscillation. In an interacting theory, the exact one-particle pole gives the same kind of oscillation, but with shifted frequency and changed amplitude.

The exact propagator is obtained by summing all 1PI self-energy insertions:

G~(p)=iD(p2)+i0,D(s)=sm02Π(s),s=p2.\widetilde G(p)={i\over D(p^2)+i0}, \qquad D(s)=s-m_0^2-\Pi(s), \qquad s=p^2.

A stable scalar particle of physical mass MM appears as an isolated zero of the exact inverse propagator:

D(M2)=0.D(M^2)=0.

Thus

M2=m02+Π(M2).\boxed{ M^2=m_0^2+\Pi(M^2). }

This equation is implicit. The self-energy must be evaluated at the physical pole, not at the bare mass unless one is working only to the lowest nontrivial order.

To find the residue, expand the inverse propagator around s=M2s=M^2:

D(s)=D(M2)+D(M2)(sM2)+O((sM2)2).D(s)=D(M^2)+D'(M^2)(s-M^2)+O((s-M^2)^2).

Since D(M2)=0D(M^2)=0,

D(s)=(1Π(M2))(sM2)+O((sM2)2).D(s)=(1-\Pi'(M^2))(s-M^2)+O((s-M^2)^2).

Therefore

G~(p)i(1Π(M2))(p2M2)+i0.\widetilde G(p) \simeq {i\over (1-\Pi'(M^2))(p^2-M^2)+i0}.

The pole residue is

Z=11Π(M2).\boxed{ Z={1\over 1-\Pi'(M^2)}. }

So near the physical pole,

G~(p)iZp2M2+i0.\boxed{ \widetilde G(p) \simeq {iZ\over p^2-M^2+i0}. }

The exact inverse propagator has a zero at the physical mass and slope inverse to the residue

The physical mass MM is the zero of the exact inverse propagator D(s)=sm02Π(s)D(s)=s-m_0^2-\Pi(s). The slope at the zero is D(M2)=1Π(M2)=Z1D'(M^2)=1-\Pi'(M^2)=Z^{-1}, so the propagator has residue iZiZ.

The mass shift and the residue are different effects. The value of Π\Pi at the pole moves the pole; the derivative of Π\Pi at the pole changes the normalization of the pole. Constant self-energy insertions renormalize the mass but do not produce wavefunction renormalization. Momentum-dependent self-energy insertions produce both.

In perturbation theory one often expands around the bare mass, but the clean definition is the pole definition above. The residue must be evaluated at s=M2s=M^2, not merely at s=m02s=m_0^2, unless the difference is beyond the order being computed.

The exact propagator also has a nonperturbative interpretation. Insert a complete set of exact energy-momentum eigenstates between the two fields. For t>0t>0,

G(t,p)=d3xeipx0ϕ(t,x)ϕ(0)0G(t,\mathbf p) =\int d^3x\,e^{-i\mathbf p\cdot\mathbf x} \langle 0|\phi(t,\mathbf x)\phi(0)|0\rangle

contains terms of the form

0ϕ(0)n,p22En,peiEn,pt,{|\langle 0|\phi(0)|n,\mathbf p\rangle|^2\over 2E_{n,\mathbf p}} e^{-iE_{n,\mathbf p}t},

with the appropriate normalization of states. If there is a stable one-particle state p|\mathbf p\rangle of mass MM, its energy is

Ep=p2+M2.E_{\mathbf p}=\sqrt{\mathbf p^2+M^2}.

We define ZZ by

0ϕ(0)p=Z\boxed{ \langle 0|\phi(0)|\mathbf p\rangle=\sqrt Z }

in the standard relativistic normalization

pq=(2π)32Epδ(3)(pq).\langle \mathbf p|\mathbf q\rangle=(2\pi)^3\,2E_{\mathbf p}\,\delta^{(3)}(\mathbf p-\mathbf q).

Then the one-particle part of the exact two-point function is

G(t,p)Z2EpeiEpt(t>0),G(t,\mathbf p)\supset {Z\over 2E_{\mathbf p}}e^{-iE_{\mathbf p}t} \qquad (t>0),

which Fourier transforms to

G~(p)iZp2M2+i0.\widetilde G(p)\supset {iZ\over p^2-M^2+i0}.

Thus ZZ measures the amplitude squared for the field ϕ\phi to create the physical one-particle state from the vacuum. A badly chosen field may have a small overlap with the particle even though the particle is present in the spectrum. A field with the wrong quantum numbers has Z=0Z=0 for that particle.

A local field creates a physical one-particle state plus multiparticle continuum states

A local field acting on the vacuum does not create only a single bare quantum. In an interacting theory it creates the physical one-particle sector with amplitude Z\sqrt Z, together with continuum states carrying the same quantum numbers. In the diagram, dΠp=d3p/[(2π)32Ep]d\Pi_{\mathbf p}=d^3\mathbf p/[(2\pi)^3 2E_{\mathbf p}] is the invariant one-particle measure and Pcont\mathcal P_{\rm cont} projects onto the multiparticle continuum.

This explains the historical phrase “wavefunction renormalization.” The word is slightly misleading in relativistic QFT: we are not renormalizing a Schrödinger wavefunction. We are renormalizing the normalization of the field operator relative to the physical one-particle state.

Spectral representation and the bound on Z

Section titled “Spectral representation and the bound on Z”

The exact two-point function admits the Källén–Lehmann representation

G~(p)=0dμ2ρ(μ2)ip2μ2+i0,\widetilde G(p) = \int_0^\infty d\mu^2\, {\rho(\mu^2)\,i\over p^2-\mu^2+i0},

where ρ(μ2)\rho(\mu^2) is a positive spectral density in a unitary theory with positive-norm states. If the lightest state created by ϕ\phi is a stable particle of mass MM, the spectral density contains a delta function. For clarity, assume here that it is the only isolated stable state below the first continuum threshold:

ρ(μ2)=Zδ(μ2M2)+ρcont(μ2).\rho(\mu^2) = Z\,\delta(\mu^2-M^2)+\rho_{\rm cont}(\mu^2).

Therefore

G~(p)=iZp2M2+i0+μth2dμ2ρcont(μ2)ip2μ2+i0.\boxed{ \widetilde G(p) = {iZ\over p^2-M^2+i0} + \int_{\mu_{\rm th}^2}^{\infty}d\mu^2\, {\rho_{\rm cont}(\mu^2)\,i\over p^2-\mu^2+i0}. }

Here μth\mu_{\rm th} is the multiparticle threshold in the same quantum-number channel. The first term is the isolated particle pole. The integral is the branch cut associated with multiparticle states.

If several stable particles with the same quantum numbers lie below the continuum, the first term is replaced by aZaδ(μ2Ma2)\sum_a Z_a\delta(\mu^2-M_a^2) in ρ\rho, and the propagator contains one isolated pole for each such state. Nothing in the LSZ argument changes: each asymptotic species is reduced at its own pole.

Spectral density with an isolated particle pole and a multiparticle continuum

In the single-pole case shown, the exact spectral density splits into an isolated one-particle contribution Zδ(μ2M2)Z\delta(\mu^2-M^2) and a continuum beginning at the first multiparticle threshold μth2\mu_{\rm th}^2. The pole and the cut are two pieces of the same exact two-point function.

The pole and continuum also have sharply different late-time behavior. The isolated pole gives a persistent oscillation. A continuum contribution dephases and is controlled asymptotically by its threshold singularity. For example, suppose the physical discontinuity has the near-threshold behavior

DiscG(E)θ(EEth)(EEth)1/2.\operatorname{Disc}G(E) \propto \theta(E-E_{\rm th})(E-E_{\rm th})^{-1/2}.

Its large-time contribution then has the characteristic form

G(t)t1/2eiEthtG(t)\sim t^{-1/2}e^{-iE_{\rm th}t}

up to a convention-dependent complex coefficient. Thus a branch point produces a power-law tail, whereas an isolated stable pole produces an undamped one-particle oscillation.

For a canonically normalized scalar field, the equal-time commutation relation implies the spectral sum rule

0dμ2ρ(μ2)=1.\int_0^\infty d\mu^2\,\rho(\mu^2)=1.

Since ρcont0\rho_{\rm cont}\ge 0,

0Z1.\boxed{ 0\le Z\le 1. }

This inequality is a statement about the spectral weight of the bare canonical field. It is not a universal statement about every possible rescaled field. If we define a renormalized field

ϕR=Z1/2ϕ,\phi_R=Z^{-1/2}\phi,

then ϕR\phi_R has unit one-particle overlap by construction, and its propagator has unit residue at the pole.

There is also an important limitation: if the would-be particle is unstable, the pole moves off the physical real axis onto a different sheet. Then there is no normalizable asymptotic one-particle state and no positive real ZZ with the same direct probability interpretation. Resonance residues are still meaningful, but they are analytic-continuation data rather than probabilities.

This distinction is important in phenomenology. Stable external particles belong in the LSZ formula. Unstable resonances appear as internal structures in amplitudes and are described by pole positions and residues after analytic continuation.

The Lorentz-invariant formula Z=1/(1Π(M2))Z=1/(1-\Pi'(M^2)) assumes that the vacuum is Lorentz invariant and that the self-energy depends only on s=p2s=p^2. In some settings, especially nonrelativistic systems or finite-density media, the denominator is more naturally written as

D(ω,p)=ω2p2m02Σ(ω,p).D(\omega,\mathbf p) = \omega^2-\mathbf p^2-m_0^2-\Sigma(\omega,\mathbf p).

The dispersion relation is determined by

D(Ωp,p)=0.D(\Omega_{\mathbf p},\mathbf p)=0.

Near the positive-energy pole,

D(ω,p)(ωΩp)D(ω,p)ωω=Ωp.D(\omega,\mathbf p) \simeq (\omega-\Omega_{\mathbf p}) \left. {\partial D(\omega,\mathbf p)\over \partial\omega} \right|_{\omega=\Omega_{\mathbf p}}.

Therefore

G~(ω,p)i(ωΩp)ωDΩp.\widetilde G(\omega,\mathbf p) \simeq {i\over (\omega-\Omega_{\mathbf p}) \left.\partial_\omega D\right|_{\Omega_{\mathbf p}} }.

The residue in the ω\omega plane is

Resω=ΩpG~(ω,p)=iωDΩp.\operatorname*{Res}_{\omega=\Omega_{\mathbf p}}\widetilde G(\omega,\mathbf p) = {i\over \left.\partial_\omega D\right|_{\Omega_{\mathbf p}} }.

If Σ(ω,p)=Π(ω2p2)\Sigma(\omega,\mathbf p)=\Pi(\omega^2-\mathbf p^2), then

ωDΩp=2Ωp[1Π(M2)],\left.\partial_\omega D\right|_{\Omega_{\mathbf p}} = 2\Omega_{\mathbf p}\,[1-\Pi'(M^2)],

and the fixed-p\mathbf p residue becomes

iZ2Ωp.{iZ\over 2\Omega_{\mathbf p}}.

This is the same ZZ as the covariant pole residue. The extra 2Ωp2\Omega_{\mathbf p} is just the conversion from the variable p2p^2 to the variable p0p^0.

The reason ZZ matters for scattering is that an nn-point Green function has one exact propagator attached to each external field insertion. Near the simultaneous external one-particle poles, the connected momentum-space Green function has the schematic form

G~c(n)(p1,,pn)[r=1niZpr2M2+i0]i(2π)4δ(4) ⁣(rpr)M(p1,,pn),\widetilde G_c^{(n)}(p_1,\ldots,p_n) \sim \left[ \prod_{r=1}^n {i\sqrt Z\over p_r^2-M^2+i0} \right] i(2\pi)^4\delta^{(4)}\!\left(\sum_r p_r\right) \mathcal M(p_1,\ldots,p_n),

up to sign and crossing conventions for which momenta are taken incoming or outgoing.

This formula packages two facts. First, each field insertion creates or annihilates a physical particle with amplitude Z\sqrt Z. Second, once the external propagation has been stripped away, the remaining object is the scattering amplitude M\mathcal M.

External pole residues factor out of an n-point Green function

Near the one-particle poles, an exact Green function factorizes into external pole factors and an amputated scattering amplitude. On each leg, LSZ cancels the pole and its factor Z1/2Z^{-1/2} compensates the field–particle overlap Z\sqrt Z.

This is the bridge to LSZ reduction. Acting on an external leg with the Klein–Gordon operator contributes a factor pr2M2p_r^2-M^2, which cancels the pole. What remains is a factor Z\sqrt Z from that external leg. If we use the field ϕ\phi whose propagator has residue ZZ, LSZ includes a compensating Z1/2Z^{-1/2} per external particle. If instead we use the renormalized field ϕR=Z1/2ϕ\phi_R=Z^{-1/2}\phi, the pole residue is one and those compensating factors are already built into the field normalization.

Suppose the self-energy near the pole can be approximated by

Π(s)=a+b(sm02)+O((sm02)2),\Pi(s)=a+b(s-m_0^2)+O((s-m_0^2)^2),

where aa and bb are perturbatively small. The exact inverse propagator is

D(s)=sm02ab(sm02)+.D(s)=s-m_0^2-a-b(s-m_0^2)+\cdots.

The pole equation D(M2)=0D(M^2)=0 gives

(1b)(M2m02)a=0,(1-b)(M^2-m_0^2)-a=0,

so

M2=m02+a1b=m02+a+ab+.M^2=m_0^2+{a\over 1-b} =m_0^2+a+ab+\cdots.

The residue is

Z=1D(M2)=11b=1+b+O(b2).Z={1\over D'(M^2)} ={1\over 1-b} =1+b+O(b^2).

The constant part aa shifts the mass. The coefficient bb changes the residue. This is why a momentum-independent tadpole in ϕ4\phi^4 theory contributes to mass renormalization but not to wavefunction renormalization at that order.

For the exact propagator of a canonically normalized scalar field, spectral positivity gives 0<Z10<Z\le1. In the present sign convention this implies b=Π(M2)0b=\Pi'(M^2)\le0 when the linear approximation is evaluated at the pole. The algebraic toy parametrization itself does not impose that inequality unless it is required to represent such a unitary canonical two-point function.

An interacting scalar two-point function contains more information than the position of its pole. The physical mass MM is defined by the zero of the exact inverse propagator,

M2m02Π(M2)=0,M^2-m_0^2-\Pi(M^2)=0,

while the pole residue is

Z=11Π(M2).Z={1\over 1-\Pi'(M^2)}.

The same ZZ has a spectral interpretation:

0ϕ(0)p=Z,\langle 0|\phi(0)|\mathbf p\rangle=\sqrt Z,

so it measures the overlap of the chosen field with the physical one-particle state. In the Källén–Lehmann representation, the exact propagator is an isolated pole plus a continuum. For a canonically normalized scalar field in a positive-norm Hilbert space, the spectral sum rule gives 0Z10\le Z\le 1.

For scattering, the important consequence is factorization near external poles. Each external field insertion supplies a pole and a factor Z\sqrt Z. LSZ reduction removes the external poles and converts Green functions into SS-matrix elements with the correct field-normalization factors.

  • The pole location and pole residue are not the same datum. Π(M2)\Pi(M^2) fixes the mass shift; Π(M2)\Pi'(M^2) fixes ZZ.
  • The inequality 0Z10\le Z\le 1 assumes a canonically normalized field and a positive Hilbert-space metric. It is not a statement about an arbitrarily rescaled field.
  • A momentum-independent self-energy can renormalize the mass without renormalizing the field.
  • The p0p^0-plane residue and the covariant p2p^2 residue differ by a factor 2Ep2E_{\mathbf p}.
  • An unstable resonance does not have the same LSZ external-particle interpretation as a stable particle. Its pole is off the physical real axis.
  • External legs in LSZ must use the physical mass MM, not the bare mass parameter m0m_0.
  • The symbol ZZ can mean either the pole residue or a field-strength counterterm depending on context. In this page it is the physical pole residue of the chosen field.

Let

G~(p)=ip2m02Π(p2)+i0.\widetilde G(p)={i\over p^2-m_0^2-\Pi(p^2)+i0}.

Assume a stable pole at p2=M2p^2=M^2. Derive the pole residue ZZ.

Solution

Write

D(s)=sm02Π(s),s=p2.D(s)=s-m_0^2-\Pi(s), \qquad s=p^2.

The pole mass obeys

D(M2)=0.D(M^2)=0.

Expanding around s=M2s=M^2,

D(s)=D(M2)+D(M2)(sM2)+O((sM2)2).D(s)=D(M^2)+D'(M^2)(s-M^2)+O((s-M^2)^2).

The first term vanishes, and

D(M2)=1Π(M2).D'(M^2)=1-\Pi'(M^2).

Therefore

G~(p)i[1Π(M2)](p2M2)+i0=iZp2M2+i0,\widetilde G(p) \simeq {i\over [1-\Pi'(M^2)](p^2-M^2)+i0} ={iZ\over p^2-M^2+i0},

with

Z=11Π(M2).Z={1\over 1-\Pi'(M^2)}.

Use the Källén–Lehmann form

G~(p)=0dμ2ρ(μ2)ip2μ2+i0\widetilde G(p)=\int_0^\infty d\mu^2\, {\rho(\mu^2)i\over p^2-\mu^2+i0}

and the canonical sum rule

0dμ2ρ(μ2)=1\int_0^\infty d\mu^2\,\rho(\mu^2)=1

to show that 0Z10\le Z\le 1 when the spectral density contains

ρ(μ2)=Zδ(μ2M2)+ρcont(μ2).\rho(\mu^2)=Z\delta(\mu^2-M^2)+\rho_{\rm cont}(\mu^2).
Solution

Substitute the decomposition into the sum rule:

1=0dμ2[Zδ(μ2M2)+ρcont(μ2)].1=\int_0^\infty d\mu^2\, \left[ Z\delta(\mu^2-M^2)+\rho_{\rm cont}(\mu^2) \right].

This gives

1=Z+0dμ2ρcont(μ2).1=Z+\int_0^\infty d\mu^2\,\rho_{\rm cont}(\mu^2).

In a positive-norm Hilbert space, the spectral density is nonnegative:

ρcont(μ2)0.\rho_{\rm cont}(\mu^2)\ge 0.

Hence

Z1.Z\le 1.

Also the isolated pole weight itself is nonnegative, so

Z0.Z\ge 0.

Therefore

0Z1.0\le Z\le 1.

At fixed p\mathbf p, suppose an exact propagator has a stable pole

G~(ω,p)iZω2Ep2+i0.\widetilde G(\omega,\mathbf p) \simeq {iZ\over \omega^2-E_{\mathbf p}^2+i0}.

Compute the large-positive-time contribution to

G(t,p)=dω2πeiωtG~(ω,p).G(t,\mathbf p)=\int{d\omega\over2\pi}e^{-i\omega t}\widetilde G(\omega,\mathbf p).
Solution

For t>0t>0, close the ω\omega contour in the lower half-plane. The positive-energy Feynman pole lies at

ω=Epi0.\omega=E_{\mathbf p}-i0.

Using

1ω2Ep2+i0=12Ep(1ωEp+i01ω+Epi0),{1\over \omega^2-E_{\mathbf p}^2+i0} = {1\over 2E_{\mathbf p}} \left( {1\over \omega-E_{\mathbf p}+i0} - {1\over \omega+E_{\mathbf p}-i0} \right),

only the first pole contributes for t>0t>0. Therefore

G(t,p)Z2EpeiEpt.G(t,\mathbf p) \simeq {Z\over 2E_{\mathbf p}}e^{-iE_{\mathbf p}t}.

The continuum contributes additional oscillatory integrals beginning at threshold. At very large tt, the isolated pole gives the clean one-particle oscillation, while the continuum gives threshold-suppressed or branch-cut contributions.

Let a four-point connected Green function near its external poles behave as

G~c(4)[r=14iZpr2M2+i0]i(2π)4δ(4) ⁣(rpr)M.\widetilde G_c^{(4)} \sim \left[ \prod_{r=1}^4 {i\sqrt Z\over p_r^2-M^2+i0} \right]i(2\pi)^4\delta^{(4)}\!\left(\sum_r p_r\right)\mathcal M.

Show what remains after applying one factor of Z1/2(pr2M2)Z^{-1/2}(p_r^2-M^2) to each external leg.

Solution

Apply

r=14Z1/2(pr2M2)\prod_{r=1}^4 Z^{-1/2}(p_r^2-M^2)

to the pole form. The factor on each leg gives

Z1/2(pr2M2)iZpr2M2+i0iZ^{-1/2}(p_r^2-M^2) {i\sqrt Z\over p_r^2-M^2+i0} \longrightarrow i

as the on-shell limit is taken. Therefore all four external pole factors are removed, and the overlap factors Z\sqrt Z are canceled by the four Z1/2Z^{-1/2} factors. Up to the standard overall powers of ii fixed by the LSZ convention, what remains is

(2π)4δ(4) ⁣(rpr)M.(2\pi)^4\delta^{(4)}\!\left(\sum_r p_r\right)\mathcal M.

This is the core mechanism of LSZ reduction: external propagation is stripped away, leaving the invariant scattering amplitude.

  • Sidney Coleman, Lectures on Quantum Field Theory, edited by Bryan Gin-ge Chen et al., World Scientific, 2019, chapters 14–15.
  • H. Lehmann, “On the Properties of Propagation Functions and Renormalization Constants of Quantized Fields,” Il Nuovo Cimento 11 (1954), 342–357.
  • H. Lehmann, K. Symanzik, and W. Zimmermann, “On the Formulation of Quantized Field Theories,” Il Nuovo Cimento 1 (1955), 205–225.
  • Mark Srednicki, Quantum Field Theory, Cambridge University Press, 2007, sections 13–14 and 27.
  • K. Symanzik, “On the Many-Particle Structure of Green’s Functions in Quantum Field Theory,” Journal of Mathematical Physics 1 (1960), 249–273.
  • Steven Weinberg, The Quantum Theory of Fields, Volume I: Foundations, Cambridge University Press, 1995, sections 10.3 and 10.7.
  • A. Zee, Quantum Field Theory in a Nutshell, 2nd edition, Princeton University Press, 2010, chapter III.3.