Continuum Wave-Function Solver For Grasp

Relativistic V(N-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 radial distances is automatically determined to provide at least 10 grid points per half cycle of the wavefunction. Lagrange multipliers are automatically determined.