A full-pivoting algorithm for the Cholesky decomposition of two-electron repulsion and spin-orbit coupling integrals.

A significant reduction in the computational effort for the evaluation of the electronic repulsion integrals (ERI) in ab initio quantum chemistry calculations is obtained by using Cholesky decomposition (CD), a numerical procedure that can remove the zero or small eigenvalues of the ERI positive (semi)definite matrix, while avoiding the calculation of the entire matrix. Conversely, due to its antisymmetric character, CD cannot be directly applied to the matrix representation of the spatial part of the two-electron spin-orbit coupling (2e-SOC) integrals. Here, we present a computational strategy to achieve a Cholesky representation of the spatial part of the 2e-SOC integrals, and propose a new efficient CD algorithm for both ERI and 2e-SOC integrals. The proposed algorithm differs from previous CD implementations by the extensive use of a full-pivoting design, which allows a univocal definition of the Cholesky basis, once the CD δ threshold is made explicit. We show that 2δ is the upper limit for the errors affecting the reconstructed 2e-SOC integrals. The proposed strategy was implemented in the ab initio program Computational Emulator of Rare Earth Systems (CERES), and tested for computational performance on both the ERI and 2e-SOC integrals evaluation. © 2017 Wiley Periodicals, Inc.