6. Radioactive Decays#

6.1. Introduction#

As we saw in the previous unit, nuclei in excited states can emit radiation in the form of gamma rays when nucleons move between different shell levels. This explains one of our three dominant sources of radiation. What about the others, β and α decay?

These become a bit trickier to solve, as both decays lead to a fundamental change in the nucleon content inside the nucleus itself. As a reminder, in all cases of radioactive decay, we characterize the process in terms of three important numbers:

  • The probability of decay per unit time λ

  • The average time a nucleus lasts before decay τ=1/λ

  • The time for half the sample to decay (half-life)

T1/2=τln(2)=0.693τ

Radioactive decay itself is a stochastic process. That is why we define it in terms of the average times it takes bulk samples of material to decay. The main thing that governs the rate of this decay in a nucleus is actually whether the process is energetically favorable or not. In gamma radiation we saw it was almost always energetically favorable for an excited nuclear state to drop down some energy levels by emitting radiation.

What about other forms of radiation? In alpha and beta radiation the nucleus physically changes from one type to another. It turns out that a nucleus is likely to be unstable and decay through these channels if the net energy of the final products Ef is less than that of the starting nucleus Ei. This energy deficit is defined as:

(6.1)#Q=M(A,Z)2c2[M(AX,ZY)c2+M(X,Y)c2]>0

Decay is possible only if the Q value is positive, and the decay is not disallowed by any other selection rule or conservation law. The energy released in the decay will come out in the kinetic energy of the release fragments in the decay. Let’s consider both beta decay and alpha decay in turn to see how theoretical concepts of the nucleus help us understand their behavior.

6.2. Beta Decay#

In beta decay we typically see experimentally an emission of an Electron from the nucleus combined with a transmutation of a nucleus from one type into another with a lower neutron number. This is referred to as beta-minus decay and typically occurs in neutron rich nuclei which are trying to lower their energy by rebalancing their proton/neutron ratios.

In the early years of Beta Decay studies the process observed was highly mysterious and showed very different a behavior to Alpha Decay which was being studied at a similar time. This is illustrated by the Cloud Chamber image shown below.

Bubble Chamber example of Beta Decay

Fig. 6.1 Bubble chamber photo of He-6 decay into Li-6, showing what appears to be non-conservation of momentum.#

We see the decay of 6He into 6Li and a charged electron. We see two tracks, curving because there is a magnetic field applied. But from analysis of the plot it looks like energy, momentum and even Total Angular Momentum I are not conserved if we assume the interaction:

6He  6Li+e

Of course the reason is that there is a Neutrino emitted as well, but this escapes the detector and we don’t see it. There are actually four things we need to consider in beta-minus(-plus) decay not just the neutrino.

  1. Beta Decay is a Three-body Process, unlike alpha decay (two body). We have to account for the energies and behavior of 3 objects - the nucleus, the beta and neutrino.

  2. The electron (positron) and neutrino produced do not exist before the decay, we need to account for their formation.

  3. The electron and neutrino are relativistic, unlike the alpha in alpha decay. This makes calculations of the energies and kinematics of the decay more complicated.

  4. The energy of the beta emitted by the nucleus in question is not the same in every decay. We get a continuum of energies put to a maximum. In contrast to alpha decay and gamma decays where we get one “Mono energetic” particle depending on the energy levels or Q value considered. The quark level interaction is given in the figure below. Typically the ‘Matrix Element’ is calculated based on the interacting particles considering the d quark in the initial state and u quark and electron and anti-neutrino in the final state.

Alt text

An example to illustrate this behavior is:

83210Bi84210Po+e+ν¯e

Note how with an electron emitted in beta decay we get an Anti-neutrino (this is Lepton Number conservation in particle physics)

Beta decay is a common type of radioactive decay that involves the transformation of nucleons within the nucleus, resulting in the emission of beta particles and a change in the nuclear charge. There are three main types of beta decay: electron emission (β decay), positron emission (β+ decay), and electron capture (EC). Each of these processes alters the nucleus in distinct ways, driven by the weak nuclear force. Here’s a breakdown of their differences:

6.2.1. Electron Emission (β Decay)#

  • Process: In β decay, a neutron inside the nucleus is transformed into a proton, an electron (e), and an electron antineutrino (ν¯e).

  • Equation: np+e+ν¯e

  • Result: The atomic number (Z) of the nucleus increases by one, while the mass number (A) remains unchanged. This process leads to the element changing to its next higher neighbor in the periodic table.

  • Driving Force: β decay occurs in neutron-rich nuclei, where the neutron-to-proton ratio is higher than that for a stable isotope of the element.

6.2.2. Positron Emission (β+ Decay)#

  • Process: In β+ decay, a proton is converted into a neutron, a positron (e+), and an electron neutrino (νe).

  • Equation: pn+e++νe

  • Result: The atomic number (Z) decreases by one, but the mass number (A) remains the same. The element changes to its previous neighbor in the periodic table.

  • Driving Force: β+ decay is favored in proton-rich nuclei, where the proton-to-neutron ratio is higher than the stable ratio for that element. It requires more energy than β decay, as the mass of the nucleus must compensate for the mass of the positron and the difference in mass between protons and neutrons.

6.2.3. Electron Capture (EC)#

  • Process: In electron capture, an inner orbital electron is captured by a proton in the nucleus, resulting in the formation of a neutron and the emission of an electron neutrino.

  • Equation: p+en+νe

  • Result: Similar to β+ decay, electron capture decreases the atomic number by one without changing the mass number, shifting the element to its previous neighbor in the periodic table.

  • Driving Force: EC also occurs in proton-rich nuclei and can be thought of as an alternative to β+ decay, especially in cases where the energy difference between the initial and final states does not favor positron emission. It strongly competes with e+ decay and is the only decay available to nuclei with a Proton Excess if the mass difference between the atoms is less than 1 MeV/c2.

6.3. Energy Considerations and Occurrence#

Note we have to be careful in Beta Decay to consider the mass of everything including the Atomic Electrons, not just the nucleons. Hence the mass used are atomic masses, not the nuclear masses. We were less careful earlier when dealing with Binding Energy calculations and the SEMF.

  • β Decay: This process is energetically favored in neutron-rich nuclei and does not require additional energy to occur since the mass of the neutron is greater than that of the proton. The Q value for this interaction is calculated based on the atomic masses as

(6.2)#Qβ=[M(A,Z)M(A,Z1)]c2
  • β+ Decay and EC: Both require the nucleus to have excess energy because the mass of a neutron is less than that of a proton (when considering the mass of the positron for β+ decay or the binding energy of the captured electron for EC). These processes are less common than β decay due to the energy requirements and occur in proton-rich nuclei.

The Q value for β+ decay based on atomic masses is calculated as

(6.3)#Qβ+=[M(A,Z)M(A,Z1)2me]c2

In the case for electron capture, the Q value needs to be calculated based on the energy of the electron that is captured based on which atomic orbital shell it is in (typically an inner shell).

6.4. Beta Decay Matrix Element#

Typically we find that the average Beta decay rate is proportional to 1/Q5 of the interaction itself. This is the Sergeant rule that derives from the theory of weak interactions. Note that this Q value is the total energy available for the decay based on the mass differences, but we expect a continuous spectrum when looking at just the electron energy for example.

As discussed in Krane, the probability of a beta decay occurring is proportional to the square of the matrix element, M, which depends on the overlap of the initial and final nuclear wave functions and the operators involved in the transition.

(6.4)#λ=2π |Mif|2 dnfdEf

where λ is the rate Mif is the nuclear matrix element, and dnfdEf is the ‘density of states’. We expect that as the density of possible states increases so to does the rate of interaction.

The nuclear matrix element is commonly written based on some interaction Hamiltonian, H^ that describes the transition,

(6.5)#|Mfi|=ψH^ψdV.

Weak interactions typically follow a VA form for the Hamiltonian, where V corresponds to vector couplings (typically the electromagnetic part in electroweak interactions), and the A corresponds to axial couplings (which carry the weak force).

If we consider that the neutrino and electron in β decay are being produced some distance r from the centre of the nucleus with momenta p (electron) and q (neutrino), then we can estimate the density of states at that momentum based on the volume of a spherical shell of thickness dp or dq as

(6.6)#dne=4πp2dpVh3
(6.7)#dnv=4πq2dqVh3

The combined density of states can then be estimated from

(6.8)#dnednv=(4π)2 V2 q2dq p2dph6

If we consider the possible wave functions for our electron and neutrino we can consider free particle solutions in the first case

ψe(r)eipr       ψν(r)eiqr

If we consider the energy scales of beta decay, the electron is expected to have an energy in the range of around 1 MeV. This corresponds to around p/0.007 fm1 much smaller than the nucleus itself. As a result we can approximate the wavefunction as

(6.9)#eipr1+ipr+....1

This last approximation where we assume the wavefunction is constant means the terms are constant in the matrix element calculations and only way that the rate depends on the momentum is in the density of states.

When we account for this our density of states can be used to arrive at a form given in Krane for the number of beta particles at a given momentum of

(6.10)#N(p)=Cc2p2(QTe)2

where C is a constant, c is the speed of light, p is the momentum of the emitted beta particle, Q is the mass deficit for the decay, and Te is the kinetic energy of the beta particle. You can see how this has a zero point both at p=0 and Q=Te corresponding to our end points of beta decay.

As shown in the figures below for Cu-64 we expect a continuous distribution as a result. In fact what we find is that additional corrections are needed for β vs β+ decay since positrons typically observe additional coulomb repulsion inside the nucleus.

Alt text

One way to understand allowed vs forbidden transitions is by producing what are called “Fermi-Kurie” plots. Based on the form shown above, if we plot the observed number of particles at a fixed momentum divided by p2 and an additional correction due to coulomb repulsion, what we find is a roughly linear fit to the data for many transitions. The point where this line crosses the x-axis determines the end point energy (maximum Te of the transition). Linearity in this fit ends up only being true when our “allowed transition” assumption are valid. Any deviation from a straight line suggest a degree of “forbidden-ness” in the interaction. Importantly a Fermi-Kurie plot can be used to estimate the neutrino mass as the end point is expected to deviate from a straight line if the neutrino mass is non-zero as shown in the Figure below.

End point of the Fermi-Kurie Plot

Fig. 6.2 End point of the Fermi-Kurie plot which plots beta particle kinetic energy vs a scaled estimate of the number of particles observed at that energy. Deviations away from a straight line suggest forbidden transitions. The end point T0 here corresponds to the maximum Q value for the decay.#

For allowed transitions, the matrix element involves simpler operators (e.g., Fermi and Gamow-Teller operators), while for forbidden transitions, more complex operators that include higher angular momentum components are required in the calculation. These higher order corrections that include the second term in the exponential expansion above are referred to as ‘forbidden transitions’.

6.5. Transitions#

In beta decay, transitions between nuclear states are classified as “allowed” or “forbidden” based on certain selection rules derived from conservation laws and quantum mechanics. These rules are related to the spin (I) and parity (P) of the initial and final nuclear states, as well as the angular momentum and parity of the emitted beta particle and neutrino system.

Allowed transitions are those that require the least change in angular momentum and have no change in parity between the initial and final states. The simplest case is the “superallowed” transition, which involves no change in nuclear spin and parity (ΔIπ=0+) and is the most favorable in terms of transition probability.

We expect that the spin of the electron and anti-neutrino (both spin half particles) can come out either parrallel or anti-parallel with one another, carrying either S=0 or S=1 respectively. The transitions where they are parallel are referred to as Gamow-Teller transitions, whilst those where they are anti-parallel are referred to as Fermi transitions.

Superallowed Transition:

ΔI=0,Δπ=none

Here, ΔI represents the change in total angular momentum quantum number, and Δπ represents the change in parity; 0+ indicates no change in either. We typically find that transitions with the shortest corrected half-life (corrections account for differences in the energy distribution and density of states) all correspond to this class of allowed transitions. This is why they are called “super-allowed”.

Allowed Transition: For “allowed” transitions (not superallowed), the changes are still minimal:

ΔI=0 or ±1,Δπ=none

except for 00 transitions, which are forbidden.

Superallowed and allowed transitions are characterized by a high transition probability because they involve minimal changes in the nuclear state, making them the most common type of beta decay. They correspond to the case where we produced the particle in our density of states at r=0 such that the wavefunction does not depend in any way on the momentum for the transition.

Forbidden Transitions

Forbidden transitions occur when the change in angular momentum and/or parity does not meet the criteria for allowed transitions. They are less likely to occur due to the need for higher changes in angular momentum between the initial and final states. Forbidden transitions are categorized by the degree of “forbiddenness,” which corresponds to the change in angular momentum of the emitted particles.

ΔI=0,±1,±2,Δπ=1

except for 00 transitions, which are still completely forbidden.

In general, the “n-th forbidden” transition is characterized by:

ΔI=n,n±1,n±2,...,Δπ=(1)n

The transition probability decreases significantly with increasing order of forbidden-ness, making higher-order forbidden transitions increasingly rare.

6.6. Selection Rules Summary#

This results in the following complete set of selection rules

  • Allowed ΔI=0,1 and Δπ=no

  • 1st Forbidden ΔI=0,1,2 and Δπ=yes

  • 2st Forbidden ΔI=2,3 and Δπ=no

  • 3rd Forbidden ΔI=3,4 and Δπ=yes

Here the 2nd forbidden transitions do actually allow for ΔI=0,1 but since these are already covered by the allowed transitions these changes in the nuclear spin are far less likely.

6.7. Example: Beta Decay of Nucleus X to Y#

Suppose we have a beta decay process where nucleus X decays to nucleus Y. Let’s say:

  • The initial state of nucleus X has a total angular momentum quantum number Ii=2 and positive parity (πi=+1).

  • The final state of nucleus Y has a total angular momentum quantum number If=1 and negative parity (πf=1).

Step 1: Calculate the Change in Angular Momentum (ΔI)

ΔI=|IfIi|=|12|=1

Step 2: Determine the Change in Parity (Δπ)

Since the parity changes from positive to negative, there is a change in parity (Δπ=+1 or 1).

Step 3: Classify the Transition

Using the calculated ΔI and Δπ, we can classify the transition. For allowed transitions, we expect ΔI=0 or ±1 (excluding 00 without a change in parity), and Δπ=0.

  • In our example, ΔI=1 fits the criterion for ΔI in allowed transitions, but the change in parity (Δπ=yes) indicates this is a forbidden transition.

  • Specifically, since there’s a change in parity and the ΔI value fits within the first forbidden category (considering ΔI=0,±1,±2 with a parity change), this would be a first forbidden transition.

Example 1 For another example of an allowed transition in beta decay, let’s consider the beta-minus decay of tritium (3H) to helium-3 (3He):

  • Tritium (3H) has one proton and two neutrons. Its nuclear spin and parity (IP) are 1/2+ in its ground state.

  • Helium-3 (3He) has two protons and one neutron. Its ground state also has a nuclear spin and parity of 1/2+.

In this decay, one of the neutrons in tritium is converted into a proton, emitting an electron (e) and an electron antineutrino (ν¯e):

3H 3He+e+ν¯e

Given that both the initial and final states have spins of 1/2, the change in angular momentum is:

ΔI=|IfIi|=|1/21/2|=0

In this case we know that the electron and anti-neutrino needed to have taken away either S=0 or S=1, so we expect this to be Fermi decay. Both the initial and final states have positive parity (+), so there is no change in parity:

Δπ=no

This means for the transition to be allowed we expect the orbital AM l to be even. Our transition meets these criteria (ΔI=0 and Δπ=no), indicating it is an allowed transition.

Example 2 Now consider the example of beta decay in 6He 6Li. In this case we expect a transition from 0+ to 1+ so we expect a change in one for the nuclear spin, but no change in parity.

We therefore know that that orbital angular momentum has to be even, requiring that the combined spin is S=1, so that we get an odd number in L+S.

This is therefore a Gamow Teller Transition, which is classified as allowed.

Example 3 Finally let’s consider the beta decay from 17N 17O where we expect a transition between nuclear spin states 12+52+. Here we are expecting ΔI=2 and Δπ=yes.

We know that the orbital angular momentum has to be odd, and therefore we expect again that the spin state must be S=1 to give an even total change in spin.

This is therefore a 1st forbidden Gamow Teller decay and is significantly less likely than the decays in the previous examples.

6.8. Double Beta Decay#

In rare cases we find that the Q value of a decay or the selection rules limit our typical β decay processes. For example in the case of 128-Te, we find that the next nuclei on the beta decay chain would be 128-I, which would give a negative Q value. However a direct transition to 128-Xe would give a positive Q value, and is enough to allow the production of two β particles at once. This process whilst extremely rare (order of 1024 years) can occur naturally.

Alt text

A similar case also happens for the decay of 48-Ca, which we would expect to decay into 48-Sc. Whilst the Q can be positive in some cases, we find that all the possible transitions are high-order forbidden transitions and are extremely unlikely to occur. Instead it is far more likely that a double beta decay transition will occur (albeit with a very long half-life).

Alt text

Double beta decay is of particular importance in some high energy particle physics experiments as there is a long standing debate on whether neutrinos may in fact be their own anti-particles (Majorana neutrinos). Since they have no charge this is allowed in some models, and if so would allow a possible neutrino-less double beta decay transition to also take place in which the two neutrino’s emitted in double-beta decay annihilate with one another, resulting in both beta particles carrying the entire Q energy for the interaction. Observation of a clear excess in number of particles right at the end of the beta decay spectrum for such nuclei therefore would confirm the true nature of the neutrino.

6.9. Alpha Decay#

Finally let’s turn our attention to alpha decay. In the context of the nuclear shell model, alpha decay can be understood in terms of the rearrangement of the nucleons in the daughter nucleus and how this rearrangement affects the nuclear stability. The shell model predicts that nuclei near magic numbers (where shells are either fully filled or empty) are more stable. Thus, alpha decay often moves a nucleus closer to these magic numbers, thereby increasing its stability.

While the selection rules for alpha decay are not as strictly defined as for beta decay in terms of angular momentum and parity changes, the energy and stability considerations play a crucial role. The key factors influencing alpha decay include:

  1. Q-Value: The energy release (Q-value) in alpha decay is a critical factor. The decay can only occur if the Q-value is positive, meaning the mass of the parent nucleus is greater than the combined mass of the daughter nucleus and the alpha particle.

  2. Nuclear Shell Effects: The shell model predicts that nuclei with a full shell of protons or neutrons (magic numbers) are particularly stable. Alpha decay is found for mostly higher A neutron rich nuclei and results in a daughter nucleus that is closer to these magic numbers.

  3. Overlap and Tunneling: The probability of alpha decay is also determined by the quantum mechanical tunneling effect. As we will see alpha particles need to tunnel out of a potential barrier if they are to escape the nucleus.

Based on our understanding of the nucleus there are some important questions we need to ask for the alpha decay process:

  • Why are the α particles emitted and not some other combination of nucleons?

  • Why is there such a range of lifetimes and energies? What is the relationship between the two?

  • What determines the frequency of α decay?

  • How are the α particles formed in the nucleus and emitted?

  • Do alphas carry information on nuclear structure?

One plot that begins to answer some of these questions is a comparison of the Q energy for different alpha-emitting particles separated by different radioactive decay series. This is called the Geiger-Nuttall relation, and it demonstrates that actually for many different series, there is a strong correlation between the radioactive half-life and the energy of the emitted alpha.

Alt text

Notice also how steep the relation is in the Geiger-Nuttall plot (it’s a log scale). An alpha emitter that produces 8 MeV alpha particles is likely to have a half-life over 22 orders of magnitude shorter than one which produces alphas of 4 MeV kinetic energy. The Geiger-Nuttall relation is a fit to the observed data of the form:

log10t1/2=b1ZQ1/2+b2

where Z is the number of protons, Q is the energy released, b1 and b2 are constants. An equivalent form is:

lnλ=a1ZQ1/2+a2

It turns out that trying to understand why the data follows this shape leads us once again to the need for a quantum mechanical treatment of Alpha Decay. To explain the relation, it turns out we need to consider that the alpha forms inside the nucleus and then escapes by Quantum Mechanical Barrier Penetration.

The alpha is envisioned as a free-particle wavefunction subject to a combination of the Nuclear Potential binding it to the nucleus, and a coulomb potential that repels it from the other charged particles in the nucleus. Typically we expect a coulomb repulsion to push the charged alpha away from the rest of the nucleus, but in this case we are imagining that the alpha is inside the rest of the nucleus, therefore the remaining protons form a charged barrier that the alpha needs to tunnel through to escape. The Coulomb barrier arises because the electrostatic force is repulsive between two positively charged entities—the alpha particle and the rest of the nucleus.

  • As the alpha particle gets closer to the nuclear surface (from within), the electrostatic potential energy increases because the repulsive force increases. At the surface of the nucleus, this potential energy reaches a maximum, forming the Coulomb barrier.

  • Inside the Nucleus: While inside the nucleus, the alpha particle is bound by the strong nuclear force, which is attractive and much stronger than the Coulomb force at short distances. This force keeps the nucleons (including those in the alpha particle) together.

  • Escape Requirement: For the alpha particle to escape, it must have enough energy to overcome the maximum potential energy at the surface created by the Coulomb repulsion. Typically, the alpha particle’s kinetic energy (within the nucleus) is insufficient to overcome this barrier due to the energy binding it inside the nucleus.

Alt text

Of course, in a classical model, if the Potential Barrier is too high, the alpha never escapes. In Quantum Mechanics, we allow Tunneling with a small probability that the alpha can penetrate the barrier. Calculating the ratio of the wavefunctions inside and outside and squaring gives the probability. We can do a simple calculation of this by considering the diagram of the potential below. We assume the alpha is pre-formed in the nucleus, oscillates around, repeatedly hitting the edge inside the nucleus until at some point it escapes by Barrier Penetration. Once formed, the alpha particle finds itself in a potential well created by the nuclear forces. Outside this well, there’s a potential energy barrier formed by the Coulomb repulsion between the positively charged alpha particle and the rest of the nucleus. Classically, the alpha particle would need to have kinetic energy greater than this barrier to escape the nucleus. However, typically, the alpha particle’s kinetic energy is less than the height of this Coulomb barrier

Alt text

Notice how the Potential Barrier drops as 1/r. Rather than try to solve this, we can make a simple approximation that the barrier has a fixed height U0 and fixed length L, and that the alpha has energy E, as illustrated below. Also, the potential seen by the alpha particle is spherically symmetric, and so the quantum mechanical problem can be approached by first separating the variables into just a radial component.

Alt text As a particle approaches the barrier, it is described by a Free Particle Wavefunction. When it reaches the barrier, it must satisfy the simplified Schrödinger Equation in the form:

(6.11)#22md2dr2ψ(r)=[EU0]ψ(r)

Solving this QM equation of simple Barrier Penetration is beyond our scope here. What is important is the result. It gives us a value for X, the probability of barrier penetration as:

(6.12)#X=AeαL

where α=(2m(U0E))/2 is called the Gamow Factor, and A is a normalization constant we will call 1 for now to make the derivation simpler. Now, if we assume the alpha makes contact with the barrier once every second, then the average probability that the alpha leaves (a decay takes place) and the half-life is just related to the tunneling probability by:

(6.13)#λ=X
(6.14)#T1/2=ln(2)/λ=ln(2)/X

If the alpha hits the barrier faster than 1 Hz, then we expect the probability that it leaves the nucleus to scale upwards as well. If we can work out the frequency that the alpha hits the barrier, which will be determined by its velocity and the size of the nucleus, and we know the height of the barrier U0 and the width L and the energy of the alpha E, then we can find the half-life T1/2. We can then see if this agrees with the Geiger-Nuttall observation from data.

6.10. Alpha 212-Po Example#

The way to proceed is illustrated below, where we attempt to find the half-life of 212-Po, which emits an 8.78 MeV alpha. First, we need to find the width of the barrier. For this, we need to know the Nuclear Separation R, i.e., the distance between the center of the nucleus and the center of the alpha. Be careful that at this point the nucleus we are talking about is not A=212, but A=208, because we have removed the alpha.

Alt text

This results in an estimate of the nuclear seperation as

(6.15)#R=1.2(A)1/3+1.2(208)1/3=9.01 fm

We now find the height of the barrier by applying Coulomb’s Law at the nuclear separation as illustrated here

(6.16)#(k2e(Z2)e)r=(2(842)1.44 MeV fm)9.01 fm=26.21 MeV

where k is equal to coulombs constant k=1/(4πε0)=8.987552×109Nm2/C2. Now the distance at which the Coulomb Potential drops to the level of energy of the observed alpha is:

(6.17)#8.78MeV=2(82)1.44 MeV~fmrr=26.9fm

So the barrier width L and Gammow factor is:

(6.18)#L=26.90fm9.01 fm=17.9 fm
(6.19)#α=2m(26.28.78)2=1.879 fm1

When assuming a ‘top-hat’ potential barrier this results in a transmission probability of

(6.20)#X=eαL=e(1.879)(17.9)=2.47×1015

The real frequency of the alpha hitting the barrier can also be estimated by determining how many times an alpha of energy 8.78 MeV makes contact with the walls.

(6.21)#8.78MeV=mv22v=2.06×107 ms1
(6.22)#f=v2R1.14×1021s1

For a given alpha, the combined Tunnelling Probability per second for emission is therefore the product of the barrier frequency and the transmission

(6.23)#λ=fX=(2.47×1015)(1.14×1021 s1)=2.82×106 s1
(6.24)#T1/2=ln2λ=0.6932.82×106s1=2.5×107 s=0.25 μs

The result is quite good. But note we could do better by modelling the potential as a series of decreasing barriers. As shown in the figure below we consider the decaying barrier as a combination of many fixed height barriers. The transmission probability through the entire set of barriers is simply obtained by multiplying all their individual transmission probabilities together.

Alt text

This can be written as:

(6.25)#X=X1×X2×X3×...×Xn
(6.26)#X=Aexp(1(2m(UiQ))Δr)

In the limit as Δr0, this analysis turns into an integral from R up to r2 based on the known form of V(r)r1

(6.27)#X=exp(1Rr22m(V(r)Q)dr)

No we can calculate the potential difference in the Gamow factor as based on the potential

(6.28)#V(r)=2Ze24πϵ01r

We have defined rc as the point where

(6.29)#Q=V(rc)rc=2Ze24πϵ01Q

Now subbing this back into our potential tells us V(rc) based only on Q, r, and rc results in a potential dependent on Q

(6.30)#Vc(r)=Qrcr

Subbing this into our integral gives us a transmission of the form

(6.31)#X=exp(1Rr22m(QrcrQ)dr)

where we can pull the Q value out of the integral

(6.32)#X=exp(Q1/2Rr22m(rcr1)dr)=exp(G)

Rearranging this ends up in a final form where the Gamow factor in our transmission factor is given by

(6.33)#G=4Zα(2mc2Q)1/2(cos1(Rrc)Rrc(1Rrc))

In this solution

Last term goes to π2 if R<<rc

and we arrive at a relation for the Gamow factor’s correlation with Z and Q,

(6.34)#GZQ.

Based on our calculations before the half-life is still proportional to the frequency that the alpha contacts the barrier and our new more realistic integrated transmission

(6.35)#t1/2λ=feG=feDZQ

where D is a new constant we have introduced based on the full form of Gamow values above. Taking the log of the equation gives us that below which ultimately ends up matching our Geiger Nutall form,

(6.36)#log10t1/2=a+bQα1/2,

where a=ln(f) and b=DZ. This confirms why we see this behavior in the first place. The exponential drop in half life for alpha emitters is due to every alpha undergoing tunneling through a falling potential barrier at the edge of the nucleus.

The procedures shown here can be applied to many nuclei to gain an estimate on the half-life and potential barrier properties if a decay process’ Q value is known. Importantly the free parameters in the Geiger Nutall form are dependent on the specific starting isotope (indicated by the lines in the GN plot), therefore measurements of the half life and Q value of one specific isotope can be used to estimate the half-life of any other nuclei in the chain provided the Q value can be calculated to a reasonable estimate (which it can with the SEMF!).

6.11. Summary#

In this unit we’ve discussed the final dominant radioactive decay channels and the effect these can have on the nucleon content of a nucleus, and the various factors that need to be considered when trying to calculate their probability and daughter product kinetic energies.

In the next half of this course we’ll start to look at different possible reactions with nuclei in more detail before looking at how knowledge of nuclear structure and decays can be used in applied physics.