Communication: Curing basis set overcompleteness with pivoted Cholesky decompositions

Abstract

The description of weakly bound electronic states is especially difficult with atomic orbital basis sets. The diffuse atomic basis functions that are necessary to describe the extended electronic state generate significant linear dependencies in the molecular basis set, which may make the electronic structure calculations ill-convergent. We propose a method where the over-complete molecular basis set is pruned by a pivoted Cholesky decomposition of the overlap matrix, yielding an optimal low-rank approximation that is numerically stable; the pivot indices determining a reduced basis set that is complete enough to describe all the basis functions in the original over-complete basis. The method can be implemented either by a simple modification to the usual canonical orthogonalization procedure, which hides the excess functions and yields fewer efficiency benefits, or by generating custom basis sets for all the atoms in the system, yielding significant cost reductions in electronic structure calculations. The pruned basis sets from the latter choice allow accurate calculations to be performed at a lower cost even at the self-consistent field level, as illustrated on a solvated (H2O)24- anion. Our results indicate that the Cholesky procedure allows one to perform calculations with accuracies close to standard augmented basis sets with cost savings which increase with the size of the basis set, ranging from 9% fewer functions in single-ζ basis sets to 28% fewer functions in triple-ζ basis sets.

6 pages

Cited by in corpus (13)

• DQC: a Python program package for Differentiable Quantum Chemistry
• Towards Pair Atomic Density Fitting for Correlation Energies with Benchmark Accuracy
• Reproducibility of density functional approximations: how new functionals should be reported
• Automatic Generation of Accurate and Cost-efficient Auxiliary Basis Sets
• Repulsive interatomic potentials calculated at three levels of theory
• Challenges with relativistic GW calculations in solids and molecules
• A variational approach for linearly dependent moving bases in quantum dynamics: application to Gaussian functions
• OpenOrbitalOptimizer -- a reusable open source library for self-consistent field calculations
• Shortcut to Chemically Accurate Quantum Computing via Density-based Basis-set Correction
• Simultaneous Spin and Point-Group Adaptation in Exact Diagonalization of Spin Clusters
• Projective spin adaptation for the exact diagonalization of isotropic spin clusters
• Importance profiles. Visualization of atomic basis set requirements
• Encoding electronic ground-state information with variational even-tempered basis sets