research26 min studyLesson 81 in Physics

Relativistic decays and phase space

Derive two-body decay kinematics and separate an interaction amplitude from the states that conservation laws permit.

On this page

What you’ll understand

  • Compute daughter energies and threshold behavior in a parent's rest frame.
  • Use a consistent Lorentz-invariant phase-space normalization.
  • Distinguish partial width, branching fraction, and laboratory decay length.

A decay rate contains two different stories. Dynamics determines how strongly an initial state couples to a final state, while phase space counts the allowed final momenta. Separating them prevents a common error: assuming that an allowed interaction produces a large rate even when almost no momentum configurations are available. We use natural units with the speed of light and the reduced Planck constant equal to one, and the metric with positive time component.

Solve conservation before calculating a diagram

Place a parent with mass MM at rest. For a two-body decay, momentum conservation makes the daughter three-momenta equal in magnitude and opposite in direction. Combining energy conservation with each on-shell condition fixes both daughter energies. There remains an orientation to integrate over, but not an independently adjustable daughter energy. If M<m1+m2M<m_1+m_2, no real on-shell final momenta exist; a nonzero vertex cannot overcome this threshold.

The square root in the momentum formula is often written with the Källén polynomial. Its factorized form is useful because one factor vanishes at the physical threshold. For equal daughters, the momentum becomes M2/4−m2\sqrt{M^2/4-m^2}. As the masses approach M/2M/2 from below, the available momentum shrinks continuously. Additional angular-momentum selection rules can suppress the rate by further powers of this momentum, so threshold suppression is not always just one square root.

E1=M2+m12−m222M,p∗=[M2−(m1+m2)2][M2−(m1−m2)2]2ME_1=\frac{M^2+m_1^2-m_2^2}{2M},\qquad p_*=\frac{\sqrt{[M^2-(m_1+m_2)^2][M^2-(m_1-m_2)^2]}}{2M}
These are rest-frame values for on-shell daughters. The positive real root applies only in the kinematically allowed region.↔ Scroll wide equations horizontally.

Keep the normalization attached to the measure

Choose single-particle normalization ⟨p′∣p⟩=2E(2π)3δ3(p′−p)\langle p'|p\rangle=2E(2\pi)^3\delta^3(\mathbf p'-\mathbf p). In the convention below the four-dimensional delta function carries its own factor (2π)4(2\pi)^4. Some references place that factor outside their definition of phase space, so copying a formula without its definition can create an entire power of 2π2\pi error. Average over the parent spin when appropriate and sum over unobserved final spins.

For distinguishable daughters and an angle-independent spin-summed or spin-averaged amplitude, integration gives Γ=p∗∣M∣2‾/(8πM2)\Gamma=p_*\overline{|\mathcal M|^2}/(8\pi M^2). An identical-particle final state needs its symmetry factor, and an angularly dependent amplitude must actually be integrated. A width has dimensions of energy; the amplitude's mass dimension depends on the number of external fields. This dimensional check is especially helpful when comparing decay and scattering formulas with different flux factors.

dΓ=∣M∣2‾2M dΦn,dΦn=(2π)4δ(4) ⁣(P−∑ipi)∏id3pi(2π)32Eid\Gamma=\frac{\overline{|\mathcal M|^2}}{2M}\,d\Phi_n,\qquad d\Phi_n=(2\pi)^4\delta^{(4)}\!\left(P-\sum_i p_i\right)\prod_i\frac{d^3p_i}{(2\pi)^3 2E_i}
This definition includes the full delta-function normalization. Symmetry factors for identical daughters are additional.↔ Scroll wide equations horizontally.

Connect a width to a detected displacement

The total width sums exclusive partial widths, and a branching fraction divides one partial width by that total. In the exponential-decay approximation the mean proper lifetime is 1/Γ1/\Gamma. A boosted parent travels a mean laboratory distance βγ/Γ\beta\gamma/\Gamma in natural units. Restore a factor of ℏc\hbar c when converting an inverse energy to meters. Proper lifetime and laboratory lifetime therefore answer different measurement questions.

A measured branching fraction also depends on selecting and correcting event samples. Daughter acceptance can change with the parent boost or decay angle; observing twice as many candidates does not automatically mean twice the underlying partial width. Near broad resonances or overlapping channels, a naive production-times-branching-factor description can fail. State the narrow-width and exponential assumptions before interpreting an apparent reconstructed peak as an independently produced unstable particle.

A fully specified hypothetical decay

A spinless parent of mass 10 GeV decays to two distinguishable daughters of mass 3 GeV. Its constant squared amplitude is 25 GeV². Find the partial width.

  1. The daughters each carry energy 5 GeV, so their common momentum is √(25−9)=4 GeV.

  2. Insert the values into Γ=p∗∣M∣2/(8πM2)\Gamma=p_*|\mathcal M|^2/(8\pi M^2): the numerator is 4×25=100 GeV34\times25=100\ \mathrm{GeV}^3 and the denominator is 800π GeV2800\pi\ \mathrm{GeV}^2.

  3. The width is 1/(8π) GeV, approximately 0.03979 GeV. No identical-particle factor is needed because the daughter species are distinct.

The result

The partial width is about 39.79 MeV. Kinematics and a dimensionful amplitude were both needed; the masses alone do not predict a rate.

Explore the model

Interactive lab

Interactive lab · vectors

Resolve a daughter momentum

Set the magnitude to represent p* in a chosen momentum unit and rotate the one displayed vector. Inspect its x and y projections. The unshown second daughter has the opposite vector in the parent rest frame. This lab illustrates orientation and components; it does not calculate phase space, boost both daughters, or sample decays.

Resolve a daughter momentumPlayback traces a fixed vector and its projections. Angle and magnitude set the geometry; the moving marker does not change them. Current readings: x component: 1.393; y component: 0.975; Magnitude: 1.7. All parameters are adjustable with the labeled controls below.−7.46−3.7303.737.46−3.5−1.7501.753.5xyvvₓ = 1.39
Vectorx / y projectionsTrace marker
x component
1.393
y component
0.975
Magnitude
1.7
t = 0 s

Playback traces a fixed vector and its projections. Angle and magnitude set the geometry; the moving marker does not change them. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forA two-body rest-frame decay fixes the momentum magnitude while leaving an orientation to integrate over.

Check your understanding

For fixed finite constant amplitude and equal daughter masses, what happens to this two-body width as M approaches 2m from above?

Try it yourself

Keep M=10 GeV and the same amplitude, but take both distinguishable daughter masses to zero. Compute the width and its ratio to the worked example.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Massless daughters each have momentum M/2; only the phase-space momentum changes.

Compare with the solution

The momentum becomes 5 GeV. The width is 125/(800π) GeV, approximately 0.04974 GeV, which is 5/4 times the previous width. This comparison assumes the amplitude really stays constant as the masses change; a physical interaction can have additional mass dependence.

Take it with you

  • Solve on-shell conservation constraints before integrating over allowed final states.
  • Phase-space conventions and identical-particle factors must accompany rate formulas.
  • Widths determine proper lifetimes; boosts and acceptance determine observed displacements and yields.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.