LIBGRPP: A Library for the Evaluation of Molecular Integrals of the Generalized Relativistic Pseudopotential Operator over Gaussian Functions

Generalized relativistic pseudopotentials (GRPP) of atomic cores implying the use of different potentials for atomic electronic shells with different principal quantum numbers give rise to accurate and reliable relativistic electronic structure models of atoms, molecules, clusters, and solids. These models readily incorporate the effects of Breit electron-electron interactions and one-loop quantum electrodynamics effects. Here, we report the computational procedure for evaluating one-electron integrals of GRPP over contracted Gaussian functions. This procedure was implemented in a library of routines named LIBGRPP, which can be integrated into existing quantum chemistry software, thus enabling the application of various methods to solve the many-electron problem with GRPPs. Pilot applications to electronic transitions in the ThO and UO2 molecules using the new library and intermediate-Hamiltonian Fock space relativistic coupled cluster method are presented. Deviations of excitation energies obtained within the GRPP approach from their all-electron Dirac-Coulomb-Gaunt counterparts do not exceed 50 cm(-1) for the 31 lowest-energy states of ThO and 110 cm(-1) for the 79 states of UO2. The results clearly demonstrate that rather economical tiny-core GRPP models can exceed in accuracy relativistic all-electron models defined by Dirac-Coulomb and Dirac-Coulomb-Gaunt Hamiltonians.

Automated summary

Generalized relativistic pseudopotentials (GRPP) are effective in accurately modeling the relativistic electronic structure of atoms, molecules, clusters, and solids. This study presents a computational procedure for evaluating one-electron integrals of GRPP, implemented in a library named LIBGRPP. Pilot applications using the new library and a relativistic coupled cluster method demonstrate the accuracy of GRPP models compared to all-electron Dirac-Coulomb-Gaunt models.