# 420: A Rudimentary Model for Alpha Particle Decay

- Page ID
- 135910

Simple models for the potential energy experienced by an alpha particle in a nucleus have the form shown below. In the interest of mathematical simplicity we will not attempt to model any particular alpha emitter, but just try to capture the essentials of the quantum mechanical tunneling mechanism for alpha decay that was formulated by Gamov, Gurney and Condon in 1928.

\[ V(x) := if[(x \geq 2), 0, \frac{5}{x} ]\]

The attractive nuclear interaction (strong nuclear force) is represented by a well of depth 2.5 and range 2 in atomic units. The repulsive Coulomb interaction becomes dominant as the strong nuclear force fades for \(x\) values greater than 2. The trapped particle is assumed to have unit mass. A variational calculation will be carried out using the trial wavefunction given below. It is also possible to solve Schrödinger's equation for alpha decay by numerical integration.

Normalized trial wavefunction with variational parameter \(α\), a decay constant:

\[ \psi (x, \alpha ) := 2 \alpha^{ \frac{3}{2}} x e^{- \alpha x}\]

\[ \int_{0}^{ \infty} \psi (x, \alpha )^{2} dx |_{simplify}^{assume,~ \alpha > 0} \rightarrow 1\]

Set up the variational energy integral:

\[ E ( \alpha ) := \int_{0}^{ \infty} \psi (x, \alpha ) \frac{1}{2} \frac{d^{2}}{dx^{2}} \psi (x, \alpha ) dx + \int_{0}^{ \infty} V(x) \psi (x, \alpha )^{2} dx\]

Next the energy is minimized with respect to \(α\) numerically. This method requires a seed value for \(α\).

α := 3 α := Minimize(E, α) α:= 1.003 E(α) = 0.958

In what follows the results of the variational calculation will be displayed graphically and interpreted.

Display the results of the variational calculation:

Calculate the probability that the particle is in the classically forbidden region. This is the region where the particle's total energy is less than the potential energy.

Because the energy of the particle is 0.958, the classical forbidden region extends from \(x = 2\) to \(x = 5.22\).

\[ \frac{5}{x} = 0.958 |_{solve,~x}^{float,~3} \rightarrow 5.22\]

Probability in classically forbidden region:

\[ \int_{2}^{5.22} \psi (x, \alpha )^{2} dx = 0.234\]

Calculate the probability that the particle has tunneled beyond the classically forbidden region.

\[ \int_{5.22}^{ \infty} \psi (x, \alpha )^{2} dx \approx 1.879 \times 10^{-3}\]

Calculate the probability that the particle is still in the nucleus.

\[ \int_{0}^{2} \psi (x, \alpha )^{2} dx \approx 0.764\]