A B-spline Galerkin method for the Dirac equation Original Research Article

The B-spline Galerkin method is first investigated for the simple eigenvalue problem, image, that can also be written as a pair of first-order equations image, image. Expanding both image and image in the image basis results in many spurious solutions such as those observed for the Dirac equation. However, when image is expanded in the image basis and image in the image basis, solutions of the well-behaved second-order differential equation are obtained. From this analysis, we propose a stable method (image) basis for the Dirac equation and evaluate its accuracy by comparing the computed and exact R-matrix for a wide range of nuclear charges Z and angular quantum numbers κ. When splines of the same order are used, many spurious solutions are found whereas none are found for splines of different order. Excellent agreement is obtained for the R-matrix and energies for bound states for low values of Z. For high Z, accuracy requires the use of a grid with many points near the nucleus. We demonstrate the accuracy of the bound-state wavefunctions by comparing integrals arising in hyperfine interaction matrix elements with exact analytic expressions. We also show that the Thomas–Reiche–Kuhn sum rule is not a good measure of the quality of the solutions obtained by the B-spline Galerkin method whereas the R-matrix is very sensitive to the appearance of pseudo-states.