Faster, more accurate single-particle energy and wave function calculations for nuclear structure
This reference describes a Fortran-based program designed to compute eigenvalues and wave functions for single-particle states in nuclear structure models. It improves speed, accuracy, and usability for bound and continuum states, using a real Woods-Saxon potential with a spin-orbit term and a refined integration approach.
The program is structured to run as a subroutine within larger nuclear-structure code, with three calculation modes and an automatic search feature that homes in on eigenvalues. It provides normalized wave functions at up to 201 radial points and includes updates that boost performance by roughly five times over earlier versions, along with checks against harmonic oscillator potentials for validation.
What you’ll gain or observe:
- A practical method for solving the radial Schrödinger equation with realistic nuclear potentials
- An automatic eigenvalue search that compares inside and outside derivatives to find matches
- Flexible input options and multiple modes to compute, normalize, or test wave functions
- Clear handling of bound and scattering states, including a 90-degree phase-shift criterion for continuum cases
- Documentation of parameters and typical values used to configure the potential, mass, and geometry
Ideal for readers of practical computational nuclear physics and researchers needing a robust tool for shell-model calculations.