Straightforward and accurate automatic auxiliary basis set generation for molecular calculations with atomic orbital basis sets

Abstract

Density fitting (DF), also known as the resolution of the identity (RI), is a widely used technique in quantum chemical calculations with various types of atomic basis sets - Gaussian-type orbitals, Slater-type orbitals, as well as numerical atomic orbitals - to speed up density functional, Hartree-Fock, and post-Hartree-Fock calculations. Traditionally, custom auxiliary basis sets are hand-optimized for each orbital basis set; however, some automatic schemes have also been suggested. In this work, we propose a simple yet numerically stable automated scheme for forming auxiliary basis sets with the help of a pivoted Cholesky decomposition, which is applicable to any type of atomic basis function. We exemplify the scheme with proof-of-concept calculations with Gaussian basis sets and show that the proposed approach leads to negligible DF/RI errors in Hartree-Fock (HF) and second-order Møller-Plesset (MP2) total energies of the non-multireference part of the W4-17 test set when used with orbital basis sets of at least polarized triple-$\zeta$ quality. The results are promising for future applications employing Slater-type orbitals or numerical atomic orbitals, as well as schemes based on the use of local fitting approaches. Global fitting approaches can also be used, in which case the high-angular-momentum functions produced by the present scheme can be truncated to minimize the computational cost.

41 pages, 3 figures

References in corpus (2)

• BAGEL: Brilliantly Advanced General Electronic-structure Library
• Communication: Curing basis set overcompleteness with pivoted Cholesky decompositions

Cited by in corpus (8)

• Cholesky decomposition of complex two-electron integrals over GIAOs: Efficient MP2 computations for large molecules in strong magnetic fields
• Towards Pair Atomic Density Fitting for Correlation Energies with Benchmark Accuracy
• Automatic Generation of Accurate and Cost-efficient Auxiliary Basis Sets
• Accurate Hellmann-Feynman forces from density functional calculations with augmented Gaussian basis sets
• Revisiting Gauge-Independent Kinetic Energy Densities in Meta-GGAs and Local Hybrid Calculations of Magnetizabilities
• Automatic generation of complementary auxiliary basis sets (CABS) for explicitly correlated methods
• OpenOrbitalOptimizer -- a reusable open source library for self-consistent field calculations
• General embedded cluster protocol for accurate modeling of oxygen vacancies in metal-oxides