Particle decay
Particle decay is the spontaneous process of one elementary particle transforming into other elementary particles. During this process, an elementary particle becomes a different particle with less mass and an intermediate particle such as W boson in muon decay. The intermediate particle then transforms into other particles. If the particles created are not stable, the decay process can continue.
Particle decay is also used to refer to the decay of hadrons. However, the term is not typically used to describe radioactive decay, in which an unstable atomic nucleus is transformed into a lighter nucleus accompanied by the emission of particles or radiation, although the two are conceptually similar.
Note that this article uses natural units, where c=\hbar=1. \,
Contents
Probability of survival and particle lifetime
Particle decay is a Poisson process, and hence the probability that a particle survives for time t before decaying is given by an exponential distribution whose time constant depends on the particle's velocity:

 P(t) = e^{t/(\gamma \tau)} \,
 where

 \tau is the mean lifetime of the particle (when at rest), and
 \gamma = \frac{1}{\sqrt{1v^2/c^2}} is the Lorentz factor of the particle.
Table of some elementary and composite particle lifetimes
All data is from the Particle Data Group.

Type Name Symbol Energy (MeV) Mean lifetime Lepton Electron / Positron e^ \, / \, e^+ 0.511 > 4.6 \times 10^{26} \ \mathrm{years} \, Muon / Antimuon \mu^ \, / \, \mu^+ 105.7 2.2\times 10^{6} \ \mathrm{seconds} \, Tau lepton / Antitau \tau^ \, / \, \tau^+ 1777 2.9 \times 10^{13} \ \mathrm{seconds} \, Meson Neutral Pion \pi^0\, 135 8.4 \times 10^{17} \ \mathrm{seconds} \, Charged Pion \pi^+ \, / \, \pi^ 139.6 2.6 \times 10^{8} \ \mathrm{seconds} \, Baryon Proton / Antiproton p^+ \, / \, p^ 938.2 > 10^{29} \ \mathrm{years} \, Neutron / Antineutron n \, / \, \bar{n} 939.6 885.7 \ \mathrm{seconds} \, Boson W boson W^+ \, / \, W^ 80,400 10^{25} \ \mathrm{seconds} \, Z boson Z^0 \, 91,000 10^{25} \ \mathrm{seconds} \,
Decay rate
The lifetime of a particle is given by the inverse of its decay rate, \Gamma, the probability per unit time that the particle will decay. For a particle of a mass M and fourmomentum P decaying into particles with momenta p_i, the differential decay rate is given by the general formula (see also Fermi's golden rule)

 d \Gamma_n = \frac{S \left\mathcal{M} \right^2}{2M} d \Phi_n (P; p_1, p_2,\dots, p_n) \,
 where

 n is the number of particles created by the decay of the original,
 S is a combinatorial factor to account for indistinguishable final states (see below),
 \mathcal{M}\, is the invariant matrix element or amplitude connecting the initial state to the final state (usually calculated using Feynman diagrams),
 d\Phi_n \, is an element of the phase space, and
 p_i \, is the fourmomentum of particle i.
The factor S is given by

 S = \prod_{j=1}^m \frac{1}{k_j!}\,

where
 m is the number of sets of indistinguishable particles in the final state, and
 k_j \, is the number of particles of type j, so that \sum_{j=1}^m k_j = n \,.
The phase space can be determined from

 d \Phi_n (P; p_1, p_2,\dots, p_n) = (2\pi)^4 \delta^4\left(P  \sum_{i=1}^n p_i\right) \prod_{i=1}^n \frac{d^3 \vec{p}_i}{2(2\pi)^3 E_i}

where
 \delta^4 \, is a fourdimensional Dirac delta function,
 \vec{p}_i \, is the (three)momentum of particle i, and
 E_i \, is the energy of particle i.
One may integrate over the phase space to obtain the total decay rate for the specified final state.
If a particle has multiple decay branches or modes with different final states, its full decay rate is obtained by summing the decay rates for all branches. The branching ratio for each mode is given by its decay rate divided by the full decay rate.
Twobody decay
Decay rate
Say a parent particle of mass M decays into two particles, labeled 1 and 2. In the rest frame of the parent particle,
 \vec{p}_1 = \vec{p_2} = \frac{2M}, \,
which is obtained by requiring that fourmomentum be conserved in the decay, i.e.
 (M, \vec{0}) = (E_1, \vec{p}_1) + (E_2, \vec{p}_2).\,
Also, in spherical coordinates,
 d^3 \vec{p} = \vec{p}\,^2\, d\vec{p}\,\, d\phi\, d\left(\cos \theta \right). \,
Using the delta function to perform the d^3 \vec{p}_2 and d\vec{p}_1\, integrals in the phasespace for a twobody final state, one finds that the decay rate in the rest frame of the parent particle is
 d\Gamma = \frac{ \left \mathcal{M} \right^2}{32 \pi^2} \frac{\vec{p}_1}{M^2}\, d\phi_1\, d\left( \cos \theta_1 \right). \,
From two different frames
The angle of an emitted particle in the lab frame is related to the angle it has emitted in the center of momentum frame by the equation

 \tan{\theta'} = \frac{\sin{\theta}}{\gamma \left(\beta / \beta' + \cos{\theta} \right)}
3body decay
The phase space element of one particle decaying into three is

 d\Phi_3 = \frac{1}{(2\pi)^5} \delta^4(P  p_1  p_2  p_3) \frac{d^3 \vec{p}_1}{2 E_1} \frac{d^3 \vec{p}_2}{2 E_2} \frac{d^3 \vec{p}_3}{2 E_3} \,
Complex mass and decay rate
The mass of an unstable particle is formally a complex number, with the real part being its mass in the usual sense, and the imaginary part being its decay rate in natural units. When the imaginary part is large compared to the real part, the particle is usually thought of as a resonance more than a particle. This is because in quantum field theory a particle of mass M (a real number) is often exchanged between two other particles when there is not enough energy to create it, if the time to travel between these other particles is short enough, of order 1/M, according to the uncertainty principle. For a particle of mass \scriptstyle M+i\Gamma, the particle can travel for time 1/M, but decays after time of order of \scriptstyle 1/\Gamma. If \scriptstyle \Gamma > M then the particle usually decays before it completes its travel.^{[1]}
See also
Notes
 ^ "The Particle Adventures"
External references
  See page 2.
 Particle Data Group.
 "The Particle Adventure" Particle Data Group, Lawrence Berkeley National Laboratory.