Woods Saxon Potential Solver

Notes: The program above finds eigenstates of a particle in the potential given by the Hamiltonian

$~H=\frac{{\bf p}^{2}}{2\mu}+V(r)+V_{c}(r)+\frac{1}{r}\,\frac{\partial}{\partial r}V_{\rm SO}(r) \,{\bf l}\cdot {\bf s}$

The potential consists of three components, central potential, coulomb and spin orbit. The radial dependence for all of them (the charge density for the Coulomb part) is given by the Woods-Saxon form factor

$~f(r,R,a)=\left[1+\exp\left(\frac{r-R}{a}\right)\right]^{-1}$

thus depth (scaling factor), radius and diffuseness fully determine each of these three components. In case of Coulomb the scaling factor is given by the charge and code is limited to zero diffuseness.

The parameters best suited for a particular nucleus can be generated automatically by selecting number of protons and neutrons on the top of the page; note that mean-field potential for protons and neutrons is slightly different therefore select an appropriate option.

If you have question/comments fell free to drop me a line Alexander Volya

Reduced mass	[MeV]	Protons	Neutrons
Central Potential Radius Diffuseness Depth	[fm] [fm] [MeV]
Spin Orbit Potential Radius Diffuseness Depth	[fm] [fm] [MeV]
Coulomb Potential Radius Diffuseness Charge (mean field)	[fm] [fm]* [integer]