Large-scale Dirac-Fock-Breit method using density fitting and 2-spinor basis functions

We present an efficient theory and algorithm for computing four-component relativistic Dirac-Fock wave functions using the Coulomb, Gaunt, and full Breit interactions. Our implementation is based on density fitting, and is routinely applicable to systems with 100 atoms and a few heavy elements. The small components are expanded using 2-spinor basis functions. We show that the factorization of 3-index half-transformed integrals before building Coulomb and exchange matrices is essential for efficient evaluation of the Fock matrix. With the Coulomb interaction, the computational cost for evaluating the Fock operator has been found to be only 70-90 times that in the non-relativistic density-fitted Hartree-Fock method. The prefactors have been 170 and 350-450 for the Gaunt and Breit interactions, respectively. The largest molecule to which we have applied the Dirac-Fock-Coulomb method is an Ac(III) motexafin complex (130 atoms, 556 electrons, 1289 basis functions), for which one self-consistent iteration takes around 1100 s using 1024 CPU cores. In addition, we have found that, while the standard fitting basis sets are accurate for Dirac-Fock-Coulomb calculations, their accuracy is very poor for Dirac-Fock-Gaunt and Breit calculations. We report a prototype of accurate fitting basis sets for these cases.