The theoretical foundations of relativistic electronic structure theory within quantum electrodynamics (QED) and the computational basis of the atomic structure code GRASP are briefly surveyed. A class of four-component basis set is introduced, which we denote the CKG-spinor set, that enforces the charge-conjugation symmetry of the Dirac equation. This formalism has been implemented using the Gaussian function technology that is routinely used in computational quantum chemistry, including in our relativistic molecular structure code, BERTHA. We demonstrate that, unlike the kinetically matched two-component basis sets that are widely employed in relativistic quantum chemistry, the CKG-spinor basis is able to reproduce the well-known eigenvalue spectrum of point-nuclear hydrogenic systems to high accuracy for all atomic symmetry types. Calculations are reported of third- and higher-order vacuum polarization effects in hydrogenic systems using the CKG-spinor set. These results reveal that Gaussian basis set expansions are able to calculate accurately these QED effects without recourse to the apparatus of regularization and in agreement with existing methods. An approach to the evaluation of the electron self-energy is outlined that extends our earlier work using partial-wave expansions in QED. Combined with the treatment of vacuum polarization effects described in this article, these basis set methods suggest the development of a comprehensive ab initio approach to the calculation of radiative and QED effects in future versions of the GRASP code.
Read full abstract