Beta decay


Beta decay is a mode of decay in which a nucleus corrects an excess of protons or neutrons by converting one into the other. The conversion can occur in both directions, known as β+\beta^{+} and β\beta^{-}:

β+ decay:pn+e++νeβ decay:np+e+νˉe\begin{align} \beta^{+}\text{ decay:}\quad &p\to n+e^{+}+\nu_{e}\\ \beta^{-}\text{ decay:}\quad &n\to p+e^{-}+\bar\nu_{e} \end{align}

where ee^{-} is an electron, e+e^{+} a positron and νe\nu_{e} and νˉe\bar{\nu}_{e} are electron neutrinos and antineutrinos respectively.

There exists a third process associated with β\beta decay called electron capture (EC for short):

EC:p+en+νe\text{EC:}\quad p+e^{-}\to n+\nu_{e}

Electron capture does exactly what the name suggests: it captures an electron, typically from an inner electron shell, and consumes it to convert a proton into a neutron.

β\beta decay is governed by the weak interaction. Specifically, it consists of the emission of a W boson from an up or down quark, converting it to the other flavor and therefore converting a proton into a neutron or viceversa. The W boson then decays into an electron/antineutrino or positron/neutrino pair depending on its charge.

Since β\beta decay is a conversion process and does not change the number of nucleons, the decay products are guaranteed to be isobars.

Examples of β\beta decays include:

  • β\beta^{-}: \ce53131I78  \ce54131Xe77\ce{_{53}^{131}I_{78}}\ \to\ \ce{_{54}^{131}Xe_{77}} with a half-life of t1/2=8.0 dt_{1/2}=8.0\text{ d}.
  • β+\beta^{+}: \ce1325Al12  \ce1225Mg13\ce{_{13}^{25}Al_{12}}\ \to\ \ce{_{12}^{25}Mg_{13}} with a half-life of t1/2=7.2 dt_{1/2}=7.2\text{ d}.
  • EC: \ce2554Mn29  \ce2454Cr30\ce{_{25}^{54}Mn_{29}}\ \to\ \ce{_{24}^{54}Cr_{30}} with a half-life of t1/2=312 dt_{1/2}=312\text{ d}.

Mechanisms

There are three decay processes that need to analyzed: β+\beta^{+}, β\beta^{-} and EC. β+\beta^{+} and β\beta^{-} largely work the same way and can be treated in unison. EC require its own section.

β+\beta^{+} and β\beta^{-}

Start by writing the semi-empirical mass formula as

M(Z,A)=αAβZ+γZ2+δM(Z,A)=\alpha A-\beta Z+\gamma Z^{2}+\delta

where

α=mnav+asA1/3+aa4β=aa+(mnmpme)γ=aaA+acA1/3\begin{align} \alpha&=m_{n}-a_{v}+a_{s}A^{-1/3}+\frac{a_{a}}{4}\\[2pt] \beta&=a_{a}+(m_{n}-m_{p}-m_{e})\\[2pt] \gamma&=\frac{a_{a}}{A}+\frac{a_{c}}{A^{1/3}} \end{align}

This allows us to factor terms based on their proton count and this is useful because β\beta decay occurs among isobars, so we can take AA to be constant. The only variable in β\beta decay is ZZ (since NN is given by N=AZN=A-Z). Then, α,β,γ\alpha,\beta,\gamma and αA\alpha A are also all constants. Two different behaviors can be identified based on the parity of AA:

  • If AA is odd, then M(Z)M(Z) draws a single parabola.
  • If AA is even, then M(Z)M(Z) draws two distinct parabolas, one for even-even nuclei and one for odd-odd nuclei. This is due to the behavior of the δ\delta term which inverts sign depending on whether a nucleus is even-even or odd-odd. The odd-odd parabola lies 2δ2\delta above the even-even one.

In both cases, the parabola has a minimum at Z=β/2γZ=\beta/2\gamma.

Note that the mass MM is the atom mass, not the nucleus mass. This is an important distinction because atomic mass automatically takes the created electron/positron into account. Neutrino masses are negligible.

Plot Beta decay odd A.svg|90%

Decay behavior of A=101A=101 isobars. The parabola minimum corresponds with \ce44Ru\ce{_{44}Ru}, to which all the other isobars converge to. Isobars with too few protons create new ones with β\beta^{-} at the expense of protons. Conversely, ones with too many convert them into neutrons with β+\beta^{+}.

Plot Beta decay even A.svg|90%

Decay behavior of A=106A=106 isobars. The lower parabola shows even-even nuclei, the upper one odd-odd nuclei. The minimum corresponds with \ce46Pd\ce{_{46}Pd}. Each β\beta decay jumps from one parabola to the other. Notice how \ce48Cd\ce{_{48}Cd} also happens to be a local minimum to \ce47Ag\ce{_{47}Ag} and \ce49In\ce{_{49}In} and ends up being stable too. This is allowed by the 2δ2\delta energy difference between the two parabolas, making it possible for there to be local minima on top of a global one. It is possible, albeit extremely rare, for \ce48Cd\ce{_{48}Cd} to decay into \ce46Pd\ce{_{46}Pd} through double β+\beta^{+} decay.

Electron capture

Because of the uncertainty principle, an electron always has a non-zero Probability of being found inside the nucleus. As such, there is always a non-zero probability of an interaction between the electron and a proton. Since the mass of the proton is slightly smaller than that of the neutron, it's possible for a proton and electron to "merge" to form a neutron, emitting a neutrino as a byproduct.

Electron capture is more probable for heavy nuclei, as bigger nuclei get closer and closer to the innermost electron shell. As one might expect, captured electrons come almost exclusively from inner shells (1s1s, 2p2p, etc.), as they are the closest to the nucleus. The hole left in the shell due to the capture causes a cascade of state transitions as higher-shell electrons "fall" to lower shells to fill the hole. This results in a cascade of X-ray photons being released jointly with the EC event.

You might notice that β+\beta^{+} and EC result in the same outcome: one proton is converted to one neutron. This is quite relevant, as the two modes compete energetically, meaning that one of the two (the most energetically efficient one) is going to occur more often than the other despite having the same effect. The condition for EC to happen is simply

M(A,Z)>M(A,Z1)BeM(A,Z)>M(A,Z-1)-B_{e^{-}}

under the assumption that MM is the atomic mass and not the nuclear mass. BeB_{e^{-}} is the binding energy of the electron that's captured. We'll pretend it doesn't exist for now, but see > Electron capture in the nucleus for more. Meanwhile, the condition for β+\beta^{+} is

M(A,Z)>M(A,Z1)+2mec2M(A,Z)>M(A,Z-1)+2m_{e}c^{2}

There's an additional 2mec22m_{e}c^{2} that comes due to the creation of the positron (one mec2m_{e}c^{2}) and due to one electron being lost from a shell since one proton was removed and the charge is imbalanced (the other mec2m_{e}c^{2}). Thus, EC is about 2mec21 MeV2m_{e}c^{2}\simeq 1\text{ MeV} more energy-efficient than β+\beta^{+}. When the difference between M(A,Z)M(A,Z) and M(A,Z1)M(A,Z-1) is less than 1 MeV1\text{ MeV}, β+\beta^{+} is completely suppressed, whereas EC is allowed: in these nuclides, EC is the only way that protons can be converted into neutrons. In all others, both processes can occur, but the branching ratio of EC is pretty much always higher than β+\beta^{+}.

Lifetimes

β\beta decay half-lives span a massive range: from a few milliseconds to 101610^{16} years. It all depends on available energy and nuclear properties.

Probably the most important example of a decently fast β\beta^{-} decay is the free-neutron decay:

np+e+νˉen\to p+e^{-}+\bar\nu_{e}

This occurs for free neutrons floating in the vacuum, has a Q value of about 0.78 MeV0.78\text{ MeV} and has mean lifetime of around τ=889\tau=889 seconds, or about 15 minutes. This is why we don't see free neutrons floating through space: they all decay into protons far before they can get anywhere. Note that this only happens in the vacuum: in nuclei neutrons are generally kept stable by the strong force. When they do decay, it's due to circumstance.

An example of some very slow processes is \ce1940K\ce{^{40}_{19}K}. This nuclide is β\beta-unstable in every way: it can β+\beta^{+}, β\beta^{-} and EC. The branching ratios greatly favor β\beta^{-} and EC, with β+\beta^{+} being a distant possibility. The potassium-40 half-life is about t1/2=1.27×109t_{1/2}=1.27\times 10^{9} years.

Diagram Beta decay potassium-40.svg|80%|center

Possible decay paths of \ce40K\ce{^{40}K} and their branching ratios. Notice how suppressed β+\beta^{+} is compared to EC. The numbers near the energy levels are the JPJ^{P} nuclear state (total angular momentum + parity). γ\gamma refers to gamma decay. The bend in the β+\beta^{+} arrow indicates the extra 2mec22m_{e}c^{2} energy cost over EC used up in handling electrons.

The reaction formulas are

β: 1940K2040Ca+e+νˉeβ+: 1940K1840Ar+e++νeEC: 1940K+e1840Ar+νe\begin{align} \beta^{-}:\ &^{40}_{19}\text{K}\to\,^{40}_{20}\text{Ca}+e^{-}+\bar\nu_{e}\\[2pt] \beta^{+}:\ &^{40}_{19}\text{K}\to\,^{40}_{18}\text{Ar}+e^{+}+\nu_{e}\\[2pt] \text{EC}:\ &^{40}_{19}\text{K}+e^{-}\to\,^{40}_{18}\text{Ar}+\nu_{e} \end{align}

Energy conservation

Energy conservation arguments provide some useful insight in the different modes of β\beta decay. As a reference, β+\beta^{+} and β\beta^{-} have Q values in the order of 1 MeV\sim 1\text{ MeV}.

β\beta^{-} in the vacuum

A free neutron in the vacuum may β\beta^{-} decay as

np+e+νˉen\to p+e^{-}+\bar\nu_{e}

The Q value is

Q=(mnmpmemν)c20.8 MeV\boxed{Q=(m_{n}-m_{p}-m_{e}-m_{\nu})c^{2}\simeq 0.8\ \text{MeV}}

(neutrino mass is basically zero). Assuming the neutron is at rest

Q=Tp+Te+TνˉQ=T_{p}+T_{e}+T_{\bar\nu}

The proton recoil energy is known to be tiny: Tp0.3 keVT_{p}\simeq 0.3\text{ keV}, so it's negligible. The remaining 0.8 MeV0.8\text{ MeV} is shared between the electron and the antineutrino. Since the neutrino is essentially massless, it's guaranteed to be relativistic. Similarly, the electron is also surely relativistic since QQ is comparable to its rest energy of mec20.5 MeVm_{e}c^{2}\simeq 0.5\text{ MeV}. So in short we have

  • small, non-relativistic proton recoil: Ep=Tp0.3 keVE_{p}=T_{p}\simeq 0.3\text{ keV}.
  • high electron relativistic energy: Ee=Te+mec21 MeVE_{e}=T_{e}+m_{e}c^{2}\sim 1\text{ MeV}.
  • high neutrino relativistic energy: Eν=TνE_{\nu}=T_{\nu}.

β+\beta^{+} in the vacuum

A free proton in the vacuum cannot β+\beta^{+} decay. This is because in the vacuum there is no external energy input, so the proton is solely responsible for respecting conservation of energy. However, the mass (and thus rest energy) of the proton is smaller than that of the neutron, so without energy coming in from the environment, decaying in the vacuum would require creating energy from nothing, and is therefore strictly impossible.

β\beta^{-} in the nucleus

A bound neutron in some nucleus \ceZAXN\ce{^{A}_{Z}X_{N}} may β\beta^{-} decay, leading to a new daughter nucleus \ceZ+1AXN1\ce{^{A}_{Z+1}X'_{N-1}} according to

ZAXNZ+1AXN1+e+νˉe^{A}_{Z}\text{X}_{N}\to\,^{A}_{Z+1}\text{X}'_{N-1}+e^{-}+\bar\nu_{e}

The Q value is, using nuclear masses MNM_{N}:

Q=[MN(ZAX)MN(Z+1AX)me]c2.Q=[M_{N}(^{A}_{Z}\text{X})-M_{N}(^{A}_{Z+1}\text{X}')-m_{e}]c^{2}.

(ignoring neutrino mass). It's convenient to use atomic masses instead:

M(\ceAX)c2=[MN(\ceZAX)+Zme]c2i=1ZBiM(\ce{^{A}X})c^{2}=[M_{N}(\ce{^{A}_{Z}X})+Zm_{e}]c^{2}-\sum_{i=1}^{Z}B_{i}

where BiB_{i} is the binding energy of the ii-th electron. So, extracting MNM_{N},

MN(\ceZAX)c2=[M(\ceAX)Zme]c2+i=1ZBiM_{N}(\ce{^{A}_{Z}X})c^{2}=[M(\ce{^{A}X})-Zm_{e}]c^{2}+\sum_{i=1}^{Z} B_{i}

Then by substitution

Q=[M(\ceAX)Zme(M(\ceAX)(Z+1)me)me]c2+(i=1ZBii=1Z+1Bi)Q=[M(\ce{^{A}X})-Zm_{e}-(M(\ce{^{A}X'})-(Z+1)m_{e})-m_{e}]c^{2}+\left( \sum_{i=1}^{Z} B_{i}-\sum_{i=1}^{Z+1} B_{i} \right)

Notice that the electron masses cancel out. The binding energies also cancel out with the exception of a single BZ+1-B_{Z+1} term. Outer shell electronic binding energy is very small (10100 eV\sim 10-100\text{ eV}), so we might as well drop it. We're left with

Q=[M(\ceAX)M(\ceAX)]c2\boxed{Q=[M(\ce{^{A}X})-M(\ce{^{A}X}')]c^{2}}

For example, \ce210Bi\ce{^{210}Bi} decays into \ce210Po\ce{^{210}Po}. The Q value is

Q=[M(\ce210Bi)M(\ce210Po)]c2=1.161 MeVQ=[M(\ce{^{210}Bi})-M(\ce{^{210}Po})]c^{2}=1.161\text{ MeV}

β+\beta^{+} in the nucleus

A bound proton in some nucleus \ceZAXN\ce{^{A}_{Z}X_{N}} may β+\beta^{+} decay, leading to a new daughter nucleus \ceZ1AXN+1\ce{^{A}_{Z-1}X'_{N+1}} according to

ZAXNZ1AXN+1+e++νe^{A}_{Z}\text{X}_{N}\to\,^{A}_{Z-1}\text{X}'_{N+1}+e^{+}+\nu_{e}

Running the same procedure as the β\beta^{-} cases gives

Q=[m(AX)m(AX)2me]c2\boxed{Q=[m(^{A}\text{X})-m(^{A}\text{X}')-2m_{e}]c^{2}}

Just like we saw when discussing EC, there is a 2mec22m_{e}c^{2} term that distinguishes β+\beta^{+}. This comes up because, unlike in β\beta^{-}, there isn't a perfect cancellation of electron masses. This makes β+\beta^{+} QQ values universally smaller by about 1 MeV\sim 1\text{ MeV} compared to β\beta^{-} decays. Since β\beta^{-} QQ values are already about 1 MeV1\text{ MeV} in scale, this explains the rarity of β+\beta^{+}1 and why sometimes it's just not possible at all, as opposed to EC. For QQ to be positive and hence β+\beta^{+} to be allowed, the mass difference between atoms must be greater than 2mec2=1.022 MeV2m_{e}c^{2}=1.022\text{ MeV}.

Electron capture in the nucleus

A bound proton in some nucleus \ceZAXN\ce{^{A}_{Z}X_{N}} may capture an electron from a shell, typically the innermost or second innermost shells, leading to a new daughter nucleus \ceZ1AXN+1\ce{^{A}_{Z-1}X'_{N+1}} according to

ZAXN+eZ1AXN+1+νe^{A}_{Z}\text{X}_{N}+e^{-}\to\,^{A}_{Z-1}\text{X}'_{N+1}+\nu_{e}

The energy of the neutrino is insignificant, so the QQ value ends up being almost identical to the β\beta^{-} decay

Q=[m(AX)m(AX)]c2Be\boxed{Q=[m(^{A}\text{X})-m(^{A}\text{X}')]c^{2}-B_{e^{-}}}

BeB_{e^{-}} is the binding energy of the electron that's being captured. Normally you'd ignore electronic binding energies, but since we're talking about innermost shells, the most tightly bound ones, it's possible in some atoms for BeB_{e^{-}} to climb into the 100 keV=0.1 MeV\sim100\text{ keV}=0.1\text{ MeV} range, which is no longer negligible. BeB_{e^{-}} becomes higher the more shells are filled in the atom, which means that it is most important in very heavy atoms2. In smaller atoms, it can safely be neglected, in which case EC has the same QQ value as β\beta^{-}. As such, EC also often has QQ values in the range of 1 MeV\sim 1\text{ MeV}, though take care of BeB_{e^{-}}.

As we've mentioned, the capture of an electron causes a cascade of high-energy photons. The first photon has energy BeB_{e^{-}}, the second has whatever energy is emitted due to the second transition, and so on and so forth. Since you're climbing into higher shells with each photon, the energy gets progressively lower and lower.

Footnotes

  1. Remember that the likelihood of a decay is proportional to the QQ value. Bigger QQ values make the decay more energetically efficient and therefore more likely to occur. Conversely, such as in this case, a tiny QQ value makes the decay technically possible but very unlikely.

  2. For instance, the highest BeB_{e^{-}} I could find is in the 1s1s electron in uranium, which is 116 keV116\text{ keV}. See here.