Continuum wavefunction solver for GRASP

Document Type


Publication Date



Relativistic VN-1 continuum orbitals are solved numerically within the framework of the GRASP code [K. Dyall et al., Comput. Phys. Commun. 55 (1989) 425]. Normalization is accomplished with a relativistic WKB method or curve fitting, depending on whether the core is an ion or is neutral. In the normalization, the phase shift (relative to hydrogen) is also calculated. The spacing of the radial grid at large distances is automatically determined to provide at least 10 grid points per half cycle of the wavefunction. Lagrange multipliers are automatically determined. © 1993.

Publication Title

Computer Physics Communications