Cholesky decomposition of complex two-electron integrals over GIAOs: Efficient MP2 computations for large molecules in strong magnetic fields

Abstract

In large-scale quantum-chemical calculations the electron-repulsion integral (ERI) tensor rapidly becomes the bottleneck in terms of memory and disk space. When an external finite magnetic field is employed, this problem becomes even more pronounced because of the reduced permutational symmetry and the need to work with complex integrals and wave-function parameters. One way to alleviate the problem is to employ a Cholesky decomposition (CD) to the complex ERIs over gauge-including atomic orbitals. The CD scheme establishes favourable compression rates by selectively discarding linearly dependent product densities from the chosen basis set while maintaining a rigorous and robust error control. This error control constitutes the main advantage over conceptually similar methods such as density fitting which rely on employing pre-defined auxiliary basis sets. We implemented the use of the CD in the framework of finite-field (ff) Hartree-Fock and ff second-order Møller Plesset perturbation theory. Our work demonstrates that the CD compression rates are particularly beneficial in calculations in the presence of a finite magnetic field. The ff-CD-MP2 scheme enables the correlated treatment of systems with more than 2000 basis functions in strong magnetic fields within a reasonable time span.

References in corpus (6)

• Straightforward and accurate automatic auxiliary basis set generation for molecular calculations with atomic orbital basis sets
• Ab-Initio Molecular Dynamics with Screened Lorentz Forces. Part I. Calculation and Atomic Charge Interpretation of Berry Curvature
• Complex Ground-State and Excitation Energies in Coupled-Cluster Theory
• Ab-Initio Molecular Dynamics with Screened Lorentz Forces. Part II. Efficient Propagators and Rovibrational Spectra in Strong Magnetic Fields
• NMR Chemical Shift Computations at Second-Order Møller-Plesset Perturbation Theory Using Gauge-Including Atomic Orbitals and Cholesky-Decomposed Two-Electron Integrals
• A black-box, general purpose quadratic self-consistent field code with and without Cholesky Decomposition of the two-electron integrals

Cited by in corpus (7)

• Cholesky decomposition of two-electron integrals in quantum-chemical calculations with perturbative or finite magnetic fields using gauge-including atomic orbitals
• Magnetic-Translational Sum Rule and Approximate Models of the Molecular Berry Curvature
• The approximate Coupled-Cluster methods CC2 and CC3 in a finite magnetic field
• Improved modularity and new features in ipie: Toward even larger AFQMC calculations on CPUs and GPUs at zero and finite temperatures
• Insight on Gaussian basis set truncation errors in weak to intermediate magnetic fields with an approximate Hamiltonian
• Computation of NMR shieldings at the CASSCF level using gauge-including atomic orbitals and Cholesky decomposition
• Equation-of-motion coupled-cluster variants in combination with perturbative triples corrections in strong magnetic fields