Abstract
The irregular solutions of the stationary Schrödinger equation are important for the fundamental formal development of scattering theory. They are also necessary for the analytical properties of the Green function, which in practice can greatly speed up calculations. Nevertheless, they are seldom considered in numerical treatments because of their divergent behavior at origin. This divergence demands high numerical precision that is difficult to achieve, particularly for non-spherical potentials which lead to different divergence rates in the coupled angular momentum channels. Based on an unconventional treatment of boundary conditions, an integral-equation method is here developed which is capable of dealing with this problem. The available precision is illustrated by electron-density calculations for NiTi in its monoclinic B19’ structure.