Help Privacy Policy Disclaimer
  Advanced SearchBrowse




Journal Article

Relativistic Cholesky-decomposed density matrix MP2


Helmich-Paris,  Benjamin
Research Group Helmich-Paris, Max-Planck-Institut für Kohlenforschung, Max Planck Society;
Division of Theoretical Chemistry, Vrije Universiteit Amsterdam;

External Resource
No external resources are shared
Fulltext (restricted access)
There are currently no full texts shared for your IP range.
Fulltext (public)
There are no public fulltexts stored in PuRe
Supplementary Material (public)
There is no public supplementary material available

Helmich-Paris, B., Repisky, M., & Visscher, L. (2019). Relativistic Cholesky-decomposed density matrix MP2. Chemical Physics, 518(2), 38-46. doi:10.1016/j.chemphys.2018.11.009.

Cite as: https://hdl.handle.net/21.11116/0000-0002-E55E-3
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.