Skip to content

The Breit–Wigner form is the standard local description of an isolated, narrow resonance. In its simplest elastic partial-wave form,

Sℓ(E)=e2iδℓ,bg(E)E−ER−iΓ/2E−ER+iΓ/2.S_\ell(E) = e^{2i\delta_{\ell,\mathrm{bg}}(E)} \frac{ E-E_R-i\Gamma/2 }{ E-E_R+i\Gamma/2 }.

The factor has a pole at ER−iΓ/2E_R-i\Gamma/2, advances the resonant phase by π\pi, and produces a Lorentzian cross-section profile when the background is negligible. The parameters are the real pole energy ERE_R and pole width Γ\Gamma. The latter equals the observed full width at half maximum only in this elementary Lorentzian limit with nearly constant kinematic factors.

The formula is powerful because it ties together four descriptions of the same local structure:

  • a pole of the analytically continued amplitude;
  • rapid phase motion on the real energy axis;
  • a peak or peak–dip structure in an observable;
  • a lifetime scale ℏ/Γ\hbar/\Gamma.

It is not a universal fit function for every bump. Thresholds, energy-dependent phase space, nonresonant backgrounds, overlapping poles, and coupled channels can all deform the elementary profile.

Helpful background. Resonances introduces temporary trapping; scattering poles explains the analytic sheets used below.

Begin with distinguishable spinless particles in one elastic partial wave ℓ\ell. Assume:

  1. one resonance pole is much closer to the physical energy region than other poles;
  2. no threshold or branch point lies within a few widths of the resonance;
  3. nonresonant factors vary slowly across an interval of order Γ\Gamma;
  4. the resonance is narrow enough that k(E)k(E) and other kinematic prefactors can be evaluated at ERE_R when deriving the elementary line shape.

The useful control ratio is schematically

ΓΔ≪1,\frac{\Gamma}{\Delta}\ll1,

where Δ\Delta is the energy distance to the nearest omitted singularity or comparably rapid dynamical structure. A smooth background need not be small, but its variation across the resonance should be small:

Γ∣dδℓ,bgdE∣≪1.\Gamma \left| \frac{d\delta_{\ell,\mathrm{bg}}}{dE} \right| \ll1.

These conditions are more informative than saying only that the peak “looks narrow.”

Write

x=E−ER,γ=Γ2.x=E-E_R, \qquad \gamma=\frac{\Gamma}{2}.

The resonant factor is

Sres(E)=x−iγx+iγ.S_{\mathrm{res}}(E) = \frac{x-i\gamma}{x+i\gamma}.

For real EE,

∣Sres(E)∣=1,\left| S_{\mathrm{res}}(E) \right|=1,

so this factor is exactly elastic-unitary. Its analytic continuation has

pole:Ep=ER−iΓ2,zero:Ez=ER+iΓ2.\begin{aligned} \text{pole:}\quad E_p &= E_R-i\frac{\Gamma}{2}, \\ \text{zero:}\quad E_z &= E_R+i\frac{\Gamma}{2}. \end{aligned}

The lower-half-plane pole gives decay rather than exponential growth. In a continuum problem it lies on the sheet reached by analytically continuing through the appropriate threshold cut. The conjugate zero is required by the real-axis unitarity of this one-channel factor.

The full elastic element may contain a smooth phase:

Sℓ(E)=e2iδℓ,bg(E)Sres(E).S_\ell(E) = e^{2i\delta_{\ell,\mathrm{bg}}(E)} S_{\mathrm{res}}(E).

Multiplication is essential. Simply adding a pole term to an arbitrary background does not automatically preserve ∣Sℓ∣=1\lvert S_\ell\rvert=1.

Define the resonant phase through

Sres(E)=e2iδres(E).S_{\mathrm{res}}(E) = e^{2i\delta_{\mathrm{res}}(E)}.

Comparison with the pole factor gives

cot⁡δres(E)=ER−EΓ/2.\cot\delta_{\mathrm{res}}(E) = \frac{E_R-E}{\Gamma/2}.

A continuous branch can be chosen so that

δres(E)=π2+arctan⁡(2(E−ER)Γ).\delta_{\mathrm{res}}(E) = \frac{\pi}{2} + \arctan \left( \frac{2(E-E_R)}{\Gamma} \right).

It then obeys

δres⟶0(E≪ER),δres(ER)=π2,δres⟶π(E≫ER).\begin{aligned} \delta_{\mathrm{res}} &\longrightarrow 0 && (E\ll E_R), \\ \delta_{\mathrm{res}}(E_R) &= \frac{\pi}{2}, \\ \delta_{\mathrm{res}} &\longrightarrow \pi && (E\gg E_R). \end{aligned}

The phase therefore advances by π\pi across an isolated resonance. A displayed principal-value arctangent may jump by π\pi even though the physical phase is continuous; phase unwrapping is part of the analysis, not an optional plotting choice.

The complete phase shift is

δℓ(E)=δℓ,bg(E)+δres(E).\delta_\ell(E) = \delta_{\ell,\mathrm{bg}}(E) + \delta_{\mathrm{res}}(E).

A resonance is diagnosed more reliably by this rapid phase motion or by a nearby pole than by a peak alone.

Use the dimensionless partial-wave amplitude

aℓ(E)=Sℓ(E)−12i.a_\ell(E) = \frac{S_\ell(E)-1}{2i}.

With no background phase,

aℓres(E)=12i(x−iγx+iγ−1)=−γx+iγ.\begin{aligned} a_\ell^{\mathrm{res}}(E) &= \frac{1}{2i} \left( \frac{x-i\gamma}{x+i\gamma} -1 \right) \\ &= -\frac{\gamma}{x+i\gamma}. \end{aligned}

At E=ERE=E_R,

aℓres(ER)=i,a_\ell^{\mathrm{res}}(E_R)=i,

which lies at the top of the elastic unitarity circle. Its magnitude is

∣aℓres(E)∣2=Γ2/4(E−ER)2+Γ2/4.\left| a_\ell^{\mathrm{res}}(E) \right|^2 = \frac{\Gamma^2/4}{ (E-E_R)^2+\Gamma^2/4 }.

For distinguishable spinless particles, the elastic contribution of this partial wave is therefore

σℓel(E)=4πk2(2ℓ+1)Γ2/4(E−ER)2+Γ2/4.\sigma_\ell^{\mathrm{el}}(E) = \frac{4\pi}{k^2} (2\ell+1) \frac{\Gamma^2/4}{ (E-E_R)^2+\Gamma^2/4 }.

If kk is nearly constant across the resonance, the peak reaches

σℓmax⁡=4πkR2(2ℓ+1),kR=k(ER).\sigma_\ell^{\max} = \frac{4\pi}{k_R^2} (2\ell+1), \qquad k_R=k(E_R).

This is the one-channel elastic unitarity limit. Partial-Wave Cross Sections derives the general channel sums and their bounds.

Introduce the dimensionless detuning

ε=2(E−ER)Γ.\varepsilon = \frac{2(E-E_R)}{\Gamma}.

The normalized line shape is

L(ε)=σℓel(E)σℓmax⁡=11+ε2.L(\varepsilon) = \frac{ \sigma_\ell^{\mathrm{el}}(E) }{ \sigma_\ell^{\max} } = \frac{1}{1+\varepsilon^2}.

Half maximum occurs at

ε=±1,\varepsilon=\pm1,

or

E=ER±Γ2.E=E_R\pm\frac{\Gamma}{2}.

The distance between these two energies is Γ\Gamma, which is why Γ\Gamma is the full width at half maximum. The half width at half maximum is Γ/2\Gamma/2.

Breit–Wigner line shape and resonant phase shift plotted against detuning

For an isolated elastic resonance with negligible background, the same detuning controls the Lorentzian cross section and the continuous π\pi phase advance. The half-maximum points lie at ER±Γ/2E_R\pm\Gamma/2.

The elementary profile is symmetric about ERE_R. An asymmetric profile can still come from a resonance: other amplitudes may interfere with the pole term, or kinematic factors may change appreciably across the width.

The pole energy

Ep=ER−iΓ2E_p = E_R-i\frac{\Gamma}{2}

produces the time factor

exp⁡(−iEptℏ)=exp⁡(−iERtℏ)×exp⁡(−Γt2ℏ).\exp \left( -\frac{iE_p t}{\hbar} \right) = \exp \left( -\frac{iE_R t}{\hbar} \right) \times \exp \left( -\frac{\Gamma t}{2\hbar} \right).

The corresponding probability decays as

P(t)∝e−Γt/ℏ,P(t)\propto e^{-\Gamma t/\hbar},

so the exponential lifetime is

τ=ℏΓ.\tau=\frac{\hbar}{\Gamma}.

This relation concerns the pole decay law. It should not be identified term by term with every collision-time definition.

For example, the one-channel Wigner–Smith delay is

τW(E)=2ℏdδℓdE.\tau_{\mathrm W}(E) = 2\hbar \frac{d\delta_\ell}{dE}.

With a constant width and negligible background,

dδresdE=Γ/2(E−ER)2+Γ2/4,\frac{d\delta_{\mathrm{res}}}{dE} = \frac{\Gamma/2}{ (E-E_R)^2+\Gamma^2/4 },

and hence

τW(ER)=4ℏΓ.\tau_{\mathrm W}(E_R) = \frac{4\hbar}{\Gamma}.

The factor of four relative to ℏ/Γ\hbar/\Gamma is not a contradiction. The Wigner–Smith quantity compares scattering phases and wave-packet arrival times, whereas ℏ/Γ\hbar/\Gamma is the exponential survival lifetime of the pole state.

Let the background phase be approximately constant over the resonance. The elastic partial cross section remains

σℓel(E)=4πk2(2ℓ+1)sin⁡2(δℓ,bg+δres).\sigma_\ell^{\mathrm{el}}(E) = \frac{4\pi}{k^2} (2\ell+1) \sin^2 \left( \delta_{\ell,\mathrm{bg}} + \delta_{\mathrm{res}} \right).

For the detuning ε\varepsilon defined above,

sin⁡δres=11+ε2,cos⁡δres=−ε1+ε2.\begin{aligned} \sin\delta_{\mathrm{res}} &= \frac{1}{\sqrt{1+\varepsilon^2}}, \\ \cos\delta_{\mathrm{res}} &= -\frac{\varepsilon}{ \sqrt{1+\varepsilon^2} }. \end{aligned}

The interference factor is

sin⁡2(δℓ,bg+δres)=(cos⁡δℓ,bg−εsin⁡δℓ,bg)21+ε2.\sin^2 \left( \delta_{\ell,\mathrm{bg}} + \delta_{\mathrm{res}} \right) = \frac{ \left( \cos\delta_{\ell,\mathrm{bg}} - \varepsilon\sin\delta_{\ell,\mathrm{bg}} \right)^2 }{ 1+\varepsilon^2 }.

Three consequences follow immediately:

  • at E=ERE=E_R, the cross-section factor is cos⁡2δℓ,bg\cos^2\delta_{\ell,\mathrm{bg}}, not necessarily one;
  • elastic saturation occurs at
Emax⁡=ER−Γ2tan⁡δℓ,bg;E_{\max} = E_R - \frac{\Gamma}{2} \tan\delta_{\ell,\mathrm{bg}};
  • a zero can occur at
Ezero=ER+Γ2cot⁡δℓ,bg.E_{\mathrm{zero}} = E_R + \frac{\Gamma}{2} \cot\delta_{\ell,\mathrm{bg}}.

These extrema refer to the angular-momentum factor with kk frozen across the resonance. A predicted extremum far outside the local fit interval should not be trusted. The same pole can therefore produce a displaced peak, an asymmetric peak–dip profile, or primarily a dip. This is amplitude interference. The Fano profile is a standard broader framework for such asymmetric resonance shapes.

Suppose an isolated state couples to open channels a,b,…a,b,\ldots. Assign a partial width Γc\Gamma_c to each open decay channel and define

Γ=∑cΓc.\Gamma = \sum_c\Gamma_c.

Ignoring smooth channel phases for the moment, a unitary resonant channel matrix can be written locally as

Sba(E)=δba−iΓbΓaE−ER+iΓ/2.S_{ba}(E) = \delta_{ba} - \frac{ i\sqrt{\Gamma_b\Gamma_a} }{ E-E_R+i\Gamma/2 }.

For one channel, Γa=Γ\Gamma_a=\Gamma, and this reduces to the elastic pole factor. The factorized numerator expresses the fact that the state must be formed from channel aa and then decay into channel bb.

For spinless scattering in a definite partial wave, the resonant transition cross section has the form

σa→b(ℓ)(E)=πka2(2ℓ+1)×ΓaΓb(E−ER)2+Γ2/4.\sigma_{a\to b}^{(\ell)}(E) = \frac{\pi}{k_a^2} (2\ell+1) \times \frac{ \Gamma_a\Gamma_b }{ (E-E_R)^2+\Gamma^2/4 }.

Introduce branching fractions

Bc=ΓcΓ.B_c=\frac{\Gamma_c}{\Gamma}.

Then

σa→b(ℓ)(E)=4πka2(2ℓ+1)×Γ2/4(E−ER)2+Γ2/4BaBb.\sigma_{a\to b}^{(\ell)}(E) = \frac{4\pi}{k_a^2} (2\ell+1) \times \frac{ \Gamma^2/4 }{ (E-E_R)^2+\Gamma^2/4 } B_aB_b.

The resonance may be narrow yet weak in a chosen reaction because either its entrance or exit branching fraction is small.

For a resonance of total angular momentum JJ formed by incident spins s1s_1 and s2s_2, the commonly used unpolarized formation formula is

σa→b(E)=2J+1(2s1+1)(2s2+1)4πka2×Γ2/4(E−ER)2+Γ2/4BaBb.\sigma_{a\to b}(E) = \frac{2J+1}{ (2s_1+1)(2s_2+1) } \frac{4\pi}{k_a^2} \times \frac{ \Gamma^2/4 }{ (E-E_R)^2+\Gamma^2/4 } B_aB_b.

Additional symmetry, polarization, identical-particle, and channel-normalization factors must be restored for the process at hand. The scalar formula is a convention dictionary, not a substitute for deriving the measured observable.

A constant width is inappropriate when an open-channel momentum changes appreciably across the resonance. For a short-range two-body channel cc with orbital angular momentum ℓc\ell_c, threshold behavior gives schematically

Γc(E)∝kc(E)2ℓc+1\Gamma_c(E) \propto k_c(E)^{2\ell_c+1}

above threshold. The total width becomes

Γ(E)=∑cΓc(E).\Gamma(E) = \sum_c\Gamma_c(E).

A more realistic local denominator then has the structure

D(E)=E−E0−Δ(E)+i2Γ(E),D(E) = E-E_0-\Delta(E) + \frac{i}{2}\Gamma(E),

where Δ(E)\Delta(E) is an energy shift associated with the real part of the self-energy. Analyticity links Δ(E)\Delta(E) and Γ(E)\Gamma(E); changing only the imaginary part may be an expedient model but is not a complete analytic construction.

Near an ss-wave threshold, kc(E)k_c(E) has square-root behavior. The resulting cusp or rapid width variation can invalidate a symmetric Lorentzian even when a pole is present. A nearby closed or newly opened channel may require a coupled-channel parameterization rather than a constant-width Breit–Wigner form.

The practical criterion is not simply Γ/ER≪1\Gamma/E_R\ll1. One also needs

Γ≪Δthreshold,\Gamma \ll \Delta_{\mathrm{threshold}},

unless the threshold dependence is included explicitly and analytically.

In particle physics, an invariant-mass amplitude is often written schematically as

BW⁡(s)=1mBW2−s−imBWΓ(s).\operatorname{BW}(s) = \frac{1}{ m_{\mathrm{BW}}^2-s - i m_{\mathrm{BW}}\Gamma(s) }.

Here ss is the squared center-of-mass energy. For s=E\sqrt{s}=E close to mBWm_{\mathrm{BW}},

mBW2−s≃−2mBW(E−mBW),m_{\mathrm{BW}}^2-s \simeq -2m_{\mathrm{BW}} \left( E-m_{\mathrm{BW}} \right),

so the squared denominator reduces locally to the nonrelativistic Lorentzian after removal of a smooth overall factor.

Several distinctions matter:

  • mBWm_{\mathrm{BW}} and ΓBW\Gamma_{\mathrm{BW}} are parameters of a chosen real-axis parameterization;
  • a pole is defined by a complex solution sps_p of the analytically continued denominator;
  • pole mass and width may be introduced through
sp=Mp−iΓp2;\sqrt{s_p} = M_p-i\frac{\Gamma_p}{2};
  • Breit–Wigner and pole parameters agree closely only for a sufficiently narrow, isolated resonance with mild energy dependence;
  • numerators contain production couplings, spin tensors, barrier factors, and channel phase space, so the observed invariant-mass distribution is not just ∣BW⁡(s)∣2\lvert\operatorname{BW}(s)\rvert^2.

Bridge to QFT Scattering places resonance poles in relativistic amplitudes. QFT Bridge: Optical Theorem and Unitarity explains the corresponding relativistic unitarity constraints.

The normalized Lorentzian satisfies

lim⁡Γ→0+1πΓ/2(E−ER)2+Γ2/4=δ(E−ER)\lim_{\Gamma\to0^+} \frac{1}{\pi} \frac{\Gamma/2}{ (E-E_R)^2+\Gamma^2/4 } = \delta(E-E_R)

in the distributional sense. Consequently,

Γ2/4(E−ER)2+Γ2/4⟶πΓ2δ(E−ER).\frac{\Gamma^2/4}{ (E-E_R)^2+\Gamma^2/4 } \longrightarrow \frac{\pi\Gamma}{2} \delta(E-E_R).

This replacement is valid only inside an integral whose remaining factors are smooth across the width. It turns resonant production and decay into an on-shell factorization:

production×branching fraction.\text{production} \times \text{branching fraction}.

It fails when cuts, selection thresholds, detector response, parton luminosities, or interference vary significantly over the resonant region.

Take ℓ=1\ell=1, kR=2 fm−1k_R=2\,\mathrm{fm}^{-1}, and Γ=4 MeV\Gamma=4\,\mathrm{MeV}. The elastic peak is

σ1max⁡=4π(2 fm−1)2(3)=3π fm2≃9.42 fm2.\begin{aligned} \sigma_1^{\max} &= \frac{4\pi}{(2\,\mathrm{fm}^{-1})^2} (3) \\ &= 3\pi\,\mathrm{fm}^2 \simeq 9.42\,\mathrm{fm}^2. \end{aligned}

Since 1 fm2=10 mb1\,\mathrm{fm}^2=10\,\mathrm{mb}, this is about 94.2 mb94.2\,\mathrm{mb}. The half-maximum energies are ER±2 MeVE_R\pm2\,\mathrm{MeV}.

Let δℓ,bg=π/6\delta_{\ell,\mathrm{bg}}=\pi/6. At the pole-energy parameter,

σℓ(ER)σℓmax⁡=cos⁡2π6=34.\frac{ \sigma_\ell(E_R) }{ \sigma_\ell^{\max} } = \cos^2\frac{\pi}{6} = \frac{3}{4}.

The elastic maximum is shifted to

Emax⁡=ER−Γ23.E_{\max} = E_R - \frac{\Gamma}{2\sqrt{3}}.

The pole position has not moved in this constant-background model. The peak has moved because the observable depends on coherent phase addition.

If Ba=0.10B_a=0.10 and Bb=0.25B_b=0.25, then the peak cross section for a→ba\to b is only

BaBb=0.025B_aB_b=0.025

times the spin-weighted unitarity prefactor. A prominent pole can therefore be inconspicuous in one final state and clear in another.

  1. Verify that the elastic resonant factor
Sres(E)=E−ER−iΓ/2E−ER+iΓ/2S_{\mathrm{res}}(E) = \frac{E-E_R-i\Gamma/2}{ E-E_R+i\Gamma/2 }

has unit magnitude for real EE. Locate its pole and zero.

Solution

For real EE, the numerator is the complex conjugate of the denominator, so

∣Sres∣2=(E−ER)2+Γ2/4(E−ER)2+Γ2/4=1.\left| S_{\mathrm{res}} \right|^2 = \frac{ (E-E_R)^2+\Gamma^2/4 }{ (E-E_R)^2+\Gamma^2/4 } =1.

The denominator vanishes at

E=ER−iΓ2,E=E_R-i\frac{\Gamma}{2},

which is the pole. The numerator vanishes at

E=ER+iΓ2,E=E_R+i\frac{\Gamma}{2},

which is the conjugate zero.

  1. Starting from
L(E)=Γ2/4(E−ER)2+Γ2/4,L(E) = \frac{\Gamma^2/4}{ (E-E_R)^2+\Gamma^2/4 },

show that the full width at half maximum is Γ\Gamma.

Solution

The maximum is L(ER)=1L(E_R)=1. Set L(E)=1/2L(E)=1/2:

Γ2/4(E−ER)2+Γ2/4=12.\frac{\Gamma^2/4}{ (E-E_R)^2+\Gamma^2/4 } = \frac12.

Cross-multiplication gives

(E−ER)2=Γ24.(E-E_R)^2 = \frac{\Gamma^2}{4}.

Thus the two half-maximum points are ER−Γ/2E_R-\Gamma/2 and ER+Γ/2E_R+\Gamma/2. Their separation is Γ\Gamma.

  1. Assume a constant background phase δℓ,bg\delta_{\ell,\mathrm{bg}}. Derive the energy at which the elastic channel saturates its unitarity bound.
Solution

Saturation requires

δℓ,bg+δres=π2(modπ).\delta_{\ell,\mathrm{bg}} + \delta_{\mathrm{res}} = \frac{\pi}{2} \pmod{\pi}.

For the branch through the resonance,

δres=π2−δℓ,bg.\delta_{\mathrm{res}} = \frac{\pi}{2} - \delta_{\ell,\mathrm{bg}}.

Using

cot⁡δres=ER−EΓ/2\cot\delta_{\mathrm{res}} = \frac{E_R-E}{\Gamma/2}

and

cot⁡(π2−δℓ,bg)=tan⁡δℓ,bg,\cot \left( \frac{\pi}{2} - \delta_{\ell,\mathrm{bg}} \right) = \tan\delta_{\ell,\mathrm{bg}},

one obtains

Emax⁡=ER−Γ2tan⁡δℓ,bg.E_{\max} = E_R - \frac{\Gamma}{2} \tan\delta_{\ell,\mathrm{bg}}.

The observed maximum can therefore differ from the pole-energy parameter even when the background is constant.

  1. A resonance has branching fractions Ba=1/4B_a=1/4 into its entrance channel and Bb=1/2B_b=1/2 into an observed exit channel. At E=ERE=E_R, what fraction of the spin-weighted unitarity prefactor appears in σa→b\sigma_{a\to b}?
Solution

At ERE_R, the Lorentzian factor is one. The remaining channel factor is

BaBb=1412=18.B_aB_b = \frac14\frac12 = \frac18.

Thus the resonant transition reaches one eighth of the relevant spin-weighted unitarity prefactor.

  1. Derive the peak Wigner–Smith delay for the constant-width resonance and compare it with the exponential lifetime.
Solution

Differentiate the continuous resonant phase:

dδresdE=Γ/2(E−ER)2+Γ2/4.\frac{d\delta_{\mathrm{res}}}{dE} = \frac{\Gamma/2}{ (E-E_R)^2+\Gamma^2/4 }.

At ERE_R,

dδresdE∣ER=2Γ.\left. \frac{d\delta_{\mathrm{res}}}{dE} \right|_{E_R} = \frac{2}{\Gamma}.

Therefore

τW(ER)=2ℏdδresdE∣ER=4ℏΓ.\tau_{\mathrm W}(E_R) = 2\hbar \left. \frac{d\delta_{\mathrm{res}}}{dE} \right|_{E_R} = \frac{4\hbar}{\Gamma}.

The pole survival lifetime is τ=ℏ/Γ\tau=\hbar/\Gamma. The two quantities answer different timing questions, so their numerical difference is expected.

  1. Why is a constant-width Breit–Wigner form especially suspect for a resonance lying close to a two-body pp-wave threshold?
Solution

For a short-range two-body channel,

Γc(E)∝kc(E)2ℓc+1.\Gamma_c(E) \propto k_c(E)^{2\ell_c+1}.

For a pp wave, ℓc=1\ell_c=1, so

Γc(E)∝kc(E)3.\Gamma_c(E)\propto k_c(E)^3.

The momentum begins at zero and varies rapidly just above threshold. The imaginary part of the denominator is therefore strongly energy dependent across the resonance. The associated real self-energy shift also varies because analyticity links the two. A constant width misses the threshold suppression, distorts the line shape, and may bias the inferred energy and width. Resonance from a Square Well supplies an exact pp-wave example in which a local pole factor succeeds only after a smooth background and a finite fit window are specified.