Next: About this document ...
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
As the typical binding energy per nucleon is 8 Mev, the deuteron is about
8 times more weakly bound than a typical nucleon.
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
(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 .
The finite spherical square well in three dimensions, as a model
of the deuteron.
When we transform the Schrödinger Wave Equation
from the cartesian co-ordinate system to the spherical polar
the radial part takes the form
We would like to investigate bound states where and
. The analytic
expressions for the solutions are :
The continuity condition at leads to
the transcendental equation
Because of the structure of and , this is effectively a family of relationships
between the depth of the potential well and the spatial extent of the well .
From the the electron scattering experiments, we know that the
size of the deuteron is fm. Inserting this value into the transcendental equation
allows us to (numerically) deduce that the depth of the potential well is MeV.
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,
Don't forget that is the reduced mass of the deuteron.
The first order coupled ODE's will be
Where the last one essentially handles the eigenvalue nature of the problem
The boundary conditions will be
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 .
You should now have developed confidence in your numerical solution.
- 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
is often used. It is simply a square well with rounded edges.
- 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.
Next: About this document ...