Hello everyone, I have been very busy lately, working on many things at the same time. Writing experiment proposal, testing a new detector array, flighting with the noise. Coding a program to read digitizer data for real time. Update the analysis code for coming experiments. Update/making some calculation codes, one of them is the calculation of Woods-Saxon energy levels.

Using RK4 method, a 2ndary 1-D differential potential can be solved easily (with some effort).

The Woods-Saxon potential with spin-orbital coupling is:

There are 6 parameters, 3 for the central term, 3 for the spin-orbital term.

And the equation is

Since the radial wavefunction has to be zero when it is far away, we can search for energy, such that the wavefunction is close to zero at a certain range. ( this is the only method I know how to practically implemented. I can imagine that, there may be some other method, transform the whole equation with a fixed r and directly solve the energy. Something like least-square method for linear system can be transformed into Matrix formalism and be directly solved. But I don’t know how)

The code is available at github (currently only for neutron).

I also have an excel to solve a single energy. Feel free to leave in comment and I will email it to you.