TY - JOUR

T1 - Relativistic Cholesky-decomposed density matrix MP2

AU - Helmich-Paris, Benjamin

AU - Repisky, Michal

AU - Visscher, Lucas

PY - 2019/2/1

Y1 - 2019/2/1

N2 - In the present article, we introduce the relativistic Cholesky-decomposed density (CDD) matrix second-order Møller-Plesset perturbation theory (MP2) energies. The working equations are formulated in terms of the usual intermediates of MP2 when employing the resolution-of-the-identity approximation (RI) for two-electron integrals. Those intermediates are obtained by substituting the occupied and virtual quaternion pseudo-density matrices of our previously proposed two-component (2C) atomic orbital-based MP2 (Helmich-Paris et al., 2016) by the corresponding pivoted quaternion Cholesky factors. While working within the Kramers-restricted formalism, we obtain a formal spin-orbit overhead of 16 and 28 for the Coulomb and exchange contribution to the 2C MP2 correlation energy, respectively, compared to a non-relativistic (NR) spin-free CDD-MP2 implementation. This compact quaternion formulation could also be easily explored in any other algorithm to compute the 2C MP2 energy. The quaternion Cholesky factors become sparse for large molecules and, with a block-wise screening, block sparse-matrix multiplication algorithm, we observed an effective quadratic scaling of the total wall time for heavy-element containing linear molecules with increasing system size. The total run time for both NR and 2C calculations was dominated by the contraction to the exchange energy. We have also investigated a bulky Te-containing supramolecular complex. For such bulky, three-dimensionally extended molecules the present screening scheme has a much larger prefactor and is less effective.

AB - In the present article, we introduce the relativistic Cholesky-decomposed density (CDD) matrix second-order Møller-Plesset perturbation theory (MP2) energies. The working equations are formulated in terms of the usual intermediates of MP2 when employing the resolution-of-the-identity approximation (RI) for two-electron integrals. Those intermediates are obtained by substituting the occupied and virtual quaternion pseudo-density matrices of our previously proposed two-component (2C) atomic orbital-based MP2 (Helmich-Paris et al., 2016) by the corresponding pivoted quaternion Cholesky factors. While working within the Kramers-restricted formalism, we obtain a formal spin-orbit overhead of 16 and 28 for the Coulomb and exchange contribution to the 2C MP2 correlation energy, respectively, compared to a non-relativistic (NR) spin-free CDD-MP2 implementation. This compact quaternion formulation could also be easily explored in any other algorithm to compute the 2C MP2 energy. The quaternion Cholesky factors become sparse for large molecules and, with a block-wise screening, block sparse-matrix multiplication algorithm, we observed an effective quadratic scaling of the total wall time for heavy-element containing linear molecules with increasing system size. The total run time for both NR and 2C calculations was dominated by the contraction to the exchange energy. We have also investigated a bulky Te-containing supramolecular complex. For such bulky, three-dimensionally extended molecules the present screening scheme has a much larger prefactor and is less effective.

UR - http://www.scopus.com/inward/record.url?scp=85057141820&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85057141820&partnerID=8YFLogxK

U2 - 10.1016/j.chemphys.2018.11.009

DO - 10.1016/j.chemphys.2018.11.009

M3 - Article

AN - SCOPUS:85057141820

VL - 518

SP - 38

EP - 46

JO - Chemical Physics

JF - Chemical Physics

SN - 0301-0104

ER -