University of the Witwatersrand Advanced Techniques in Physics : Part 2 - 2004 Project 2 :
Study of the Deuteron |

**The Deuteron** is the simplest nucleus of all. It is for nuclear physics
what the hydrogen atom is for atomic physics.
The binding energy of the deuteron is

(1) |

We would like a quantum mechanical model of the deuteron.
As usual, we think of the simplest potential first.
We choose a finite spherical square well in three dimensions
(figure 1),

(2) |

(Here is the well depth which is an approximation to the strength of the nucleon-nucleon potential and may be taken as a size parameter for the deuteron.) We consider the lowest energy state of the deuteron which has .

When we transform the Schrödinger Wave Equation
from the cartesian co-ordinate system to the spherical polar
co-ordinate system
,
the radial part takes the form

(3) |

We would like to investigate bound states where and
. The analytic
expressions for the solutions are :

(4) |

(5) |

It is clear therefore that the deuteron is only barely bound.

- Assuming a well depth of MeV, solve the Schrödinger equation (eigenstates and wave functions)
for the deuteron in the case of the square well, using the shooting method (to a fitting point).
The method of shooting to a fitting point should be used.
You should also as usual work in dimensionless variables,

(6)

(7)

(8) The first order coupled ODE's will be

(9)

Where the last one essentially handles the eigenvalue nature of the problem

(10) The boundary conditions will be

(11)

and

(12)

so that there is one unknown boundary condition at the origin and two unknown boundary conditions for large . The homogeneous nature of the problem implies that is correct to a normalisation factor. The eigenvalue spectrum is determined sequentially by guessing start values of the shooting algorithm that begin at small numbers and gradually increasing subsequent guesses, thus probing the space of eigenvalues. Clearly the first guess at will be the same as .- Verify that there is only one eigenvalue.
- Check wave functions of your numerical procedure against the analytic solutions. (Shape is sufficient, although it is possible to normalise both the analytic and the numerical solutions as usual.)

- Now consider a more general potential. The phenomenological form

(13) - Show that this solution evolves from that of the square well. This is another check of your programme.
- Is deuteron physics very sensitive to the shape of the potential ? Why do you think this is the case.

- You may make any other observations regarding other possible tests of the numerical routine or
explorations of the physics relevant to this model.