In the book, Introduction to Nuclear Reaction, by C. A. Bertulani, Chapter 2, there is a section that shows an example of phase shift from a square well. I follow the receipt and able to reproduce the result.
The actual method I did is little different. I use Mathematica to solve the radial Schrodinger equation in , the differential equation I solve is
where is the Woods-Saxon function, set the diffusiveness parameter to 0.001, the Woods-Saxon shape becomes square well.
In solving the differential equation, I set the right-hand side to be zero when . That solved to problem when , it has to deal with the . In fact, since , the right-hand side is always zero when .
For the solving range, I set
This is because the “first” peak of the Riccati-Bessel function appear after . This will make sure the range has 5 oscillations in the solution.
Next, I table the solution for the last 2 oscillations, find the maximum, normalize the solution, and then fit with
and extract , and .
The phase shift is
.
The negative is due to the is negatively defined in Mathematica.
The elastics cross section is
From the calculation, we found that, at very small energy, only s-wave scattering cross section is non-zero, which is agree with theory.
[20220531] update
The Mathematical notebook can be downloaded here.