Skip to main content
Physics LibreTexts

7.A: Alpha Decay (Project)

  • Page ID
    10407
  • \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \) \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)\(\newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\) \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\) \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\) \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\) \( \newcommand{\Span}{\mathrm{span}}\) \(\newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\) \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\) \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\) \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\) \( \newcommand{\Span}{\mathrm{span}}\)\(\newcommand{\AA}{\unicode[.8,0]{x212B}}\)

    The radioactive process known as alpha decay involves the tunneling of an alpha particle (a bound state of two protons and two neutrons) through the Coulomb barrier to escape the nuclear potential.


    In this activity, you will construct a spreadsheet that calculates the tunneling probability for the alpha particle and, using this probability, calculates the half-life of the radioactive nucleus.

    I. Constructing the Spreadsheet

    Construct a spreadsheet that has the following general form. (The template AlphaDecay is available in the PHY 262 course folder.)

    Alpha Decay
    your name:
    your name:

    alpha energy = MeV T =
    A = f =
    Z =  =
    half-life = s
    nuclear depth = MeV
    nuclear radius = fm
    exit radius = fm

    Position Potential integrand
    (fm) (MeV)

    The spreadsheet should allow you to enter the mass number (A) and atomic number (Z) of a radioactive nucleus, the kinetic energy of the emitted alpha particle, and then calculate the half-life for the transition.

    One of the several approximations you will make is to imagine the alpha particle as a point particle moving in the potential well created by the rest of the nucleus. Since the alpha particle consists of 4 nucleons (two of which are protons), the effective mass number of the rest of the nucleus is (A – 4) and the effective charge of the rest of the nucleus is (Z – 2).

    For heavy nuclei, the depth of the nuclear potential (UNuclear) is approximately

    \[U_{Nuclear} = -50\text{ MeV}\]

    In addition, the radius of the nuclear potential (R) is
    \[R = (1.2\text{ fm})(A-4)^{1/3}\]

    The electrical potential energy between the alpha particle and the rest of the nucleus (UCoulomb) is given by

    \[U_{Coulomb} = qV\]

    \[U_{Coulomb} = (2e)\bigg(\frac{k((Z-2)e)}{r}\bigg)\]
    \[U_{Coulomb} = \frac{2(Z-2)ke^2}{r}\]

    (Note that in “nuclear” units, ke2 = 1.44 MeV fm.)

    Therefore, the exit radius (RExit) can be determined by equating the alpha energy (EAlpha) to the electrical energy
    \[E_{Alpha} = U_{Coulomb}\]

    \[E_{Alpha} = \frac{2(Z-2)ke^2}{R_{Exit}}\]

    \[R_{Exit} = \frac{2(Z-2)ke^2}{E_{Alpha}}\]

    Using these definitions, complete the Position and Potential columns of the spreadsheet for 500 equally spaced points between the nuclear radius and the exit radius.


    A. Tunneling Approximation Scheme

    The probability of tunneling through a square potential barrier of height U and width x is given by:

    \[T \cong e^{-x/\delta}\]
    where
    \[ \delta = \frac{\hbar c}{\sqrt{8mc^2 (U-E)}}\]

    If the barrier is not square, we can approximate the barrier by a series of extremely thin barriers which are approximately square. In the limit of infinitesimally thin barriers, this becomes

    \[T \cong e^{-\int dx/\delta}\]

    \[T \cong exp\bigg(-\frac{\sqrt{8mc^2}}{\hbar c}\int \sqrt{U-E}dx\bigg)\]

    Thus, to find the probability of the alpha particle tunneling through the barrier, you should calculate, at each of the 500 equally spaced points between the nuclear radius and the exit radius, the value of the integrand,

    \[\sqrt{U-E}(\Delta x)\]

    The sum of these contributions can then be used to find T. (Note that the rest energy of the alpha particle is 3728 MeV.)

    The number of collisions per second the alpha particle makes with the nuclear “wall”, f, can be determined from the speed of the alpha particle inside the nucleus and the size of the nucleus. The decay rate, , is the probability of decay per second, and can be determined from T and f. Finally, the half-life is directed related to . The specific relationships are left for you to determine.


    II. Using the Spreadsheet

    1. Calculate the half-life for \(_{92}^{222}U\) emitting an alpha particle of energy 9.50 MeV and record it below. (Don’t worry if you don’t get exactly the answer shown in the data table. Due to the approximations we’ve made your answer should be a couple hundred times larger than the experimentally measured value.) Print the first page of your spreadsheet and attach it to the end of this activity.

    2. Although your value for the half-life is larger than the measured value, your spreadsheet can still be used to explore how half-life depends on alpha energy. For each of the even-A uranium isotopes listed in the data table, determine the theoretical half-life using your spreadsheet. Record these results in a suitably labeled column in the data table.

    3. To compare the theoretical dependence between half-life and alpha energy to the experimentally measured dependence, create a graph of the log of the half-life vs. alpha energy. Show both the theoretical and experimental data on this graph. Print this graph and attach it to the end of the activity.

    4. Does your spreadsheet accurately reflect the dependence of half-life on alpha energy? Clearly explain why or why not.

    5. Clearly explain why a reduction in alpha energy by a factor of approximately two (9.50 MeV to 4.27 MeV) can result in a change in half-life by a factor of approximately 1023 (!).


    This page titled 7.A: Alpha Decay (Project) is shared under a CC BY-NC-SA 4.0 license and was authored, remixed, and/or curated by Paul D'Alessandris.

    • Was this article helpful?