Electron Repulsion Integrals Research Articles

Tensor-Hypercontracted MP2 First Derivatives: Runtime and Memory Efficient Computation of Hyperfine Coupling Constants.

We employ our recently introduced tensor-hypercontracted (THC) second-order Møller–Plesset perturbation theory (MP2) method [Bangerter, F. H., Glasbrenner, M., Ochsenfeld, C. J. Chem. Theory Comput.2021, 17, 211–221] for the computation of hyperfine coupling constants (HFCCs). The implementation leverages the tensor structure of the THC factorized electron repulsion integrals for an efficient formation of the integral-based intermediates. The computational complexity of the most expensive and formally quintic scaling exchange-like contribution is reduced to effectively subquadratic, by making use of the intrinsic, exponentially decaying coupling between tensor indices through screening based on natural blocking. Overall, this yields an effective subquadratic scaling with a low prefactor for the presented THC-based AO-MP2 method for the computation of isotropic HFCCs on DNA fragments with up to 500 atoms and 5000 basis functions. Furthermore, the implementation achieves considerable speedups with up to a factor of roughly 600–1000 compared to previous implementations [Vogler, S., Ludwig, M., Maurer, M., Ochsenfeld, C. J. Chem. Phys.2017, 147, 024101] for medium-sized organic radicals, while also significantly reducing storage requirements.

Periodic Density Matrix Embedding for CO Adsorption on the MgO(001) Surface.

The adsorption of simple gas molecules to metal oxide surfaces is a primary step in many heterogeneous catalysis applications. Quantum chemical modeling of these reactions is a challenge in terms of both cost and accuracy, and quantum-embedding methods are promising, especially for localized chemical phenomena. In this work, we employ density matrix embedding theory (DMET) for periodic systems to calculate the adsorption energy of CO to the MgO(001) surface. Using coupled-cluster theory with single and double excitations and second-order Møller–Plesset perturbation theory as quantum chemical solvers, we perform calculations with embedding clusters up to 266 electrons in 306 orbitals, with the largest embedding models agreeing to within 1.2 kcal/mol of the non-embedding references. Moreover, we present a memory-efficient procedure of storing and manipulating electron repulsion integrals in the embedding space within the framework of periodic DMET.

The numerical evaluation of Slater integrals on graphics processing units.

This article presents SlaterGPU, a graphics processing unit (GPU) accelerated library that uses OpenACC to numerically compute Slater-type orbital (STO) integrals. The electron repulsion integrals (ERI) are computed under the RI approximation using the Coulomb potential of the Slater basis function. To fully realize the performance capabilities of modern GPUs, the Slater integrals are evaluated in mixed-precision, resulting in speedups for the ERIs of over 80×. Parallelization on multiple GPUs allows for integral throughput of over 3 million integrals per second. This places STO integral throughput within reach of single-threaded, conventional Gaussian integration schemes. To test the quality of the integrals, the fluorine exchange reaction barrier in fluoromethane was computed using heat-bath configuration interaction (HBCI). In addition, the singlet-triplet gap of cyclobutadiene was examined using HBCI in a triple- , polarized basis set. These benchmarks demonstrate the library's ability to generate the full set of integrals necessary for configuration interaction with up to functions in the auxiliary basis.

Implementation of Occupied and Virtual Edmiston-Ruedenberg Orbitals Using Cholesky Decomposed Integrals.

We present a trust-region optimization of the Edmiston–Ruedenberg orbital localization function. The approach is used to localize both the occupied and the virtual orbitals and is the first demonstration of general virtual orbital localization using the Edmiston–Ruedenberg localization function. In the Edmiston–Ruedenberg approach, the sum of the orbital self-repulsion energies is maximized to obtain the localized orbitals. The Cholesky decomposition reduces the cost of transforming the electron repulsion integrals, and the overall scaling of our implementation is . The optimization is performed with all quantities in the molecular orbital basis, and the localization of the occupied orbitals is often less expensive than the corresponding self-consistent field (SCF) optimization. Furthermore, the occupied orbital localization scales linearly with the basis set. For the virtual space, the cost is significantly higher than the SCF optimization. The orbital spreads of the resulting virtual Edmiston–Ruedenberg orbitals are larger than for other, less expensive, orbital localization functions. This indicates that other localization procedures are more suitable for applications such as local post-Hartree–Fock calculations.

Efficient implementation of molecular CCSD gradients with Cholesky-decomposed electron repulsion integrals.

We present an efficient implementation of ground and excited state coupled cluster singles and doubles (CCSD) gradients based on Cholesky-decomposed electron repulsion integrals. Cholesky decomposition and density fitting are both inner projection methods, and, thus, similar implementation schemes can be applied for both methods. One well-known advantage of inner projection methods, which we exploit in our implementation, is that one can avoid storing large V3O and V4 arrays by instead considering three-index intermediates. Furthermore, our implementation does not require the formation and storage of Cholesky vector derivatives. The new implementation is shown to perform well, with less than 10% of the time spent calculating the gradients in geometry optimizations. Furthermore, the computational time per optimization cycle is significantly lower compared to other implementations based on an inner projection method. We showcase the capabilities of the implementation by optimizing the geometry of the retinal molecule (C20H28O) at the CCSD/aug-cc-pVDZ level of theory.

Simplified GW/BSE Approach for Charged and Neutral Excitation Energies of Large Molecules and Nanomaterials.

Inspired by Grimme's simplified Tamm-Dancoff density functional theory approach [Grimme, S. J. Chem. Phys. 2013, 138, 244104], we describe a simplified approach to excited-state calculations within the GW approximation to the self-energy and the Bethe-Salpeter equation (BSE), which we call sGW/sBSE. The primary simplification to the electron repulsion integrals yields the same structure as with tensor hypercontraction, such that our method has a storage requirement that grows quadratically with system size and computational timing that grows cubically with system size. The performance of sGW is tested on the ionization potential of the molecules in the GW100 test set, for which it differs from ab initio GW calculations by only 0.2 eV. The performance of sBSE (based on the sGW input) is tested on the excitation energies of molecules in Thiel's set, for which it differs from ab initio GW/BSE calculations by about 0.5 eV. As examples of the systems that can be routinely studied with sGW/sBSE, we calculate the band gap and excitation energy of hydrogen-passivated silicon nanocrystals with up to 2650 electrons in 4678 spatial orbitals and the absorption spectra of two large organic dye molecules with hundreds of atoms.

Rank-reduced coupled-cluster. III. Tensor hypercontraction of the doubles amplitudes.

We develop a quartic-scaling implementation of coupled-cluster singles and doubles (CCSD) based on low-rank tensor hypercontraction (THC) factorizations of both the electron repulsion integrals (ERIs) and the doubles amplitudes. This extends our rank-reduced (RR) coupled-cluster method to incorporate higher-order tensor factorizations. The THC factorization of the doubles amplitudes accounts for most of the gain in computational efficiency as it is sufficient, in conjunction with a Cholesky decomposition of the ERIs, to reduce the computational complexity of most contributions to the CCSD amplitude equations. Further THC factorization of the ERIs reduces the complexity of certain terms arising from nested commutators between the doubles excitation operator and the two-electron operator. We implement this new algorithm using graphical processing units and demonstrate that it enables CCSD calculations for molecules with 250 atoms and 2500 basis functions using a single computer node. Furthermore, we show that the new method computes correlation energies with comparable accuracy to the underlying RR-CCSD method.

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

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 favorable 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 relies 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 (MP2). 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.

Full-frequency dynamical Bethe-Salpeter equation without frequency and a study of double excitations.

The Bethe-Salpeter equation (BSE) that results from the GW approximation to the self-energy is a frequency-dependent (nonlinear) eigenvalue problem due to the dynamically screened Coulomb interaction between electrons and holes. The computational time required for a numerically exact treatment of this frequency dependence is O(N6), where N is the system size. To avoid the common static screening approximation, we show that the full-frequency dynamical BSE can be exactly reformulated as a frequency-independent eigenvalue problem in an expanded space of single and double excitations. When combined with an iterative eigensolver and the density fitting approximation to the electron repulsion integrals, this reformulation yields a dynamical BSE algorithm whose computational time is O(N5), which we verify numerically. Furthermore, the reformulation provides direct access to excited states with dominant double excitation character, which are completely absent in the spectrum of the statically screened BSE. We study the 21Ag state of butadiene, hexatriene, and octatetraene and find that GW/BSE overestimates the excitation energy by about 1.5-2eV and significantly underestimates the double excitation character.

Scalable Molecular GW Calculations: Valence and Core Spectra.

We present a scalable implementation of the GW approximation using Gaussian atomic orbitals to study the valence and core ionization spectroscopies of molecules. The implementation of the standard spectral decomposition approach to the screened-Coulomb interaction, as well as a contour-deformation method, is described. We have implemented both of these approaches using the robust variational fitting approximation to the four-center electron repulsion integrals. We have utilized the MINRES solver with the contour-deformation approach to reduce the computational scaling by 1 order of magnitude. A complex heuristic in the quasiparticle equation solver further allows a speed-up of the computation of core and semicore ionization energies. Benchmark tests using the GW100 and CORE65 data sets and the carbon 1s binding energy of the well-studied ethyl trifluoroacetate, or ESCA molecule, were performed to validate the accuracy of our implementation. We also demonstrate and discuss the parallel performance and computational scaling of our implementation using a range of water clusters of increasing size.

Faster Self-Consistent Field (SCF) Calculations on GPU Clusters.

A novel implementation of the self-consistent field (SCF) procedure specifically designed for high-performance execution on multiple graphics processing units (GPUs) is presented. The algorithm offloads to GPUs the three major computational stages of the SCF, namely, the calculation of one-electron integrals, the calculation and digestion of electron repulsion integrals, and the diagonalization of the Fock matrix, including SCF acceleration via DIIS. Performance results for a variety of test molecules and basis sets show remarkable speedups with respect to the state-of-the-art parallel GAMESS CPU code and relative to other widely used GPU codes for both single and multi-GPU execution. The new code outperforms all existing multi-GPU implementations when using eight V100 GPUs, with speedups relative to Terachem ranging from 1.2× to 3.3× and speedups of up to 28× over QUICK on one GPU and 15× using eight GPUs. Strong scaling calculations show nearly ideal scalability up to 8 GPUs while retaining high parallel efficiency for up to 18 GPUs.

GPU acceleration of rank-reduced coupled-cluster singles and doubles.

We have developed a graphical processing unit (GPU) accelerated implementation of our recently introduced rank-reduced coupled-cluster singles and doubles (RR-CCSD) method. RR-CCSDintroduces a low-rank approximation of the doubles amplitudes. This is combined with a low-rank approximation of the electron repulsion integrals via Cholesky decomposition. The result of these two low-rank approximations is the replacement of the usual fourth-order CCSDtensors with products of second- and third-order tensors. In our implementation, only a single fourth-order tensor must be constructed as an intermediate during the solution of the amplitude equations. Owing in large part to the compression of the doubles amplitudes, the GPU-accelerated implementation shows excellent parallel efficiency (95% on eight GPUs). Our implementation can solve the RR-CCSDequations for up to 400 electrons and 1550 basis functions-roughly 50% larger than the largest canonical CCSDcomputations that have been performed on any hardware. In addition to increased scalability, the RR-CCSDcomputations are faster than the corresponding CCSDcomputations for all but the smallest molecules. We test the accuracy of RR-CCSDfor a variety of chemical systems including up to 1000 basis functions and determine that accuracy to better than 0.1% error in the correlation energy can be achieved with roughly 95% compression of the ov space for the largest systems considered. We also demonstrate that conformational energies can be predicted to be within 0.1kcal mol-1 with efficient compression applied to the wavefunction. Finally, we find that low-rank approximations of the CCSDdoubles amplitudes used in the similarity transformation of the Hamiltonian prior to a conventional equation-of-motion CCSDcomputation will not introduce significant errors (on the order of a few hundredths of an electronvolt) into the resulting excitation energies.

Self-Consistent Auxiliary Density Perturbation Theory.

The working equations for the extension of auxiliary density perturbation theory (ADPT) to hybrid functionals, employing the variational fitting of the Fock potential, are derived. The response equations in the resulting self-consistent ADPT (SC-ADPT) are solved iteratively with an adapted Eirola-Nevanlinna algorithm. As a result, a memory and CPU time efficient implementation of perturbation theory free of four-center electron repulsion integrals (ERIs) is obtained. Our validation calculations of SC-ADPT static and dynamic polarizabilities show quantitative agreement with corresponding coupled perturbed Hartree-Fock and Kohn-Sham results employing four-center ERIs. The comparison of SC-ADPT hybrid functional polarizabilities with coupled cluster reference calculations yield semiquantitative agreement. The presented systematic study of the dynamic polarizabilities of oligothiophenes shows that hybrid functionals can overcome the pathological misplacement of excitation poles by the local density and generalized gradient approximations. Good agreement with experimental dynamic polarizabilities for all studied oligothiophenes is achieved with range-separated hybrid functionals in the framework of SC-ADPT.

Tight distance-dependent estimators for screening two-center and three-center short-range Coulomb integrals over Gaussian basis functions.

We derive distance-dependent estimators for two-center and three-center electron repulsion integrals over a short-range Coulomb potential, erfc(ωr12)/r12. These estimators are much tighter than the ones based on the Schwarz inequality and can be viewed as a complement to the distance-dependent estimators for four-center short-range Coulomb integrals and for two-center and three-center full Coulomb integrals previously reported. Because the short-range Coulomb potential is commonly used in solid-state calculations, including those with the Heyd-Scuseria-Ernzerhof functional and with our recently introduced range-separated periodic Gaussian density fitting, we test our estimators on a diverse set of periodic systems using a wide range of the range-separation parameter ω. These tests demonstrate the robust tightness of our estimators, which are then used with integral screening to calculate periodic three-center short-range Coulomb integrals with linear scaling in system size.

Harnessing the Power of Multi-GPU Acceleration into the Quantum Interaction Computational Kernel Program.

We report a new multi-GPU capable ab initio Hartree-Fock/density functional theory implementation integrated into the open source QUantum Interaction Computational Kernel (QUICK) program. Details on the load balancing algorithms for electron repulsion integrals and exchange correlation quadrature across multiple GPUs are described. Benchmarking studies carried out on up to four GPU nodes, each containing four NVIDIA V100-SXM2 type GPUs demonstrate that our implementation is capable of achieving excellent load balancing and high parallel efficiency. For representative medium to large size protein/organic molecular systems, the observed parallel efficiencies remained above 82% for the Kohn-Sham matrix formation and above 90% for nuclear gradient calculations. The accelerations on NVIDIA A100, P100, and K80 platforms also have realized parallel efficiencies higher than 68% in all tested cases, paving the way for large-scale ab initio electronic structure calculations with QUICK.

Toward the Minimal Floating Operation Count Cholesky Decomposition of Electron Repulsion Integrals.

As quantum chemistry calculations deal with molecular systems of increasing size, the memory requirement to store electron-repulsion integrals (ERIs) greatly outpaces the physical memory available in computing hardware. The Cholesky decomposition of ERIs provides a convenient yet accurate technique to reduce the storage requirement of integrals. Recent developments of a two-step algorithm have drastically reduced the memory operation (MOP) count, leaving the floating operation (FLOP) count as the last frontier of cost reduction in the Cholesky ERI algorithm. In this report, we introduce a dynamic integral tracking, reusing, and compression/elimination protocol embedded in the two-step Cholesky ERI method. Benchmark studies suggest that this technique becomes particularly advantageous when the basis set consists of many computationally expensive high-angular-momentum basis functions. With this dynamic-ERI improvement, the Cholesky ERI approach proves to be a highly efficient algorithm with minimal FLOP and MOP count.

Development of Linear-Scaling Relativistic Quantum Chemistry Covering the Periodic Table

Abstract This Award Account focuses on the author’s studies on the theoretical developments of two-component (2c) relativistic quantum chemistry calculations for large systems with high efficiency and high accuracy, with a review of related studies as the background. The local unitary transformation scheme allows the linear-scaling computation cost to be applied to construct a 2c Hamiltonian, such as an infinite-order two-component version. The divide-and-conquer scheme can lead to linear-scaling computation costs to apply not only a Hartree-Fock (HF) method but also post-HF methods such as the second-order Møller-Plesset perturbation and couple cluster theory with singles and doubles for the 2c Hamiltonian in addition to a non-relativistic version. The frozen core potential scheme can naturally connect pseudopotential calculations with all-electron calculations. The accompanying coordinate expansion with a transfer recurrence relation scheme provides an efficient algorithm for the rapid evaluation of electron repulsion integrals for systems including heavy elements, the orbitals of which have long contractions and high angular momenta, such as f- and g-orbitals. Illustrative applications will help readers realize the advantages and usefulness of these schemes.

Smooth potential-energy surfaces in fragmentation-based local correlation methods for periodic systems

ABSTRACT Local approximations facilitate the application of post-Hartree–Fock methods in the condensed phase, but simultaneously introduce errors leading to discontinuous potential-energy surfaces. In this work, we explore how these discontinuities arise in periodic systems, their implications, and possible ways of controlling them. In addition, we present a fully periodic Divide-Expand-Consolidate second-order Møller–Plesset approach using an attenuated resolution-of-the-identity approximation for the electron repulsion integrals and a convenient class to handle translation-symmetric tensors in block-Toeplitz format.

Analytic gradients for multiconfiguration pair-density functional theory with density fitting: Development and application to geometry optimization in the ground and excited states.

Density fitting reduces the computational cost of both energy and gradient calculations by avoiding the computation and manipulation of four-index electron repulsion integrals. With this algorithm, one can efficiently optimize the geometries of large systems with an accurate multireference treatment. Here, we present the derivation of multiconfiguration pair-density functional theory for energies and analytic gradients with density fitting. Six systems are studied, and the results are compared to those obtained with no approximation to the electron repulsion integrals and to the results obtained by complete active space second-order perturbation theory. With the new approach, there is an increase in the speed of computation with a negligible loss in accuracy. Smaller grid sizes have also been used to reduce the computational cost of multiconfiguration pair-density functional theory with little effect on the optimized geometries and gradient values.

Resolution of the identity approximation applied to PNOF correlation calculations.

In this work, the required algebra to employ the resolution of the identity approximation within the Piris Natural Orbital Functional (PNOF) is developed, leading to an implementation named DoNOF-RI. The arithmetic scaling is reduced from fifth-order to fourth-order, and the memory scaling is reduced from fourth-order to third-order, allowing significant computational time savings. After the DoNOF-RI calculation has fully converged, a restart with four-center electron repulsion integrals can be performed to remove the effect of the auxiliary basis set incompleteness, quickly converging to the exact result. The proposed approach has been tested on cycloalkanes and other molecules of general interest to study the numerical results, as well as the speed-ups achieved by PNOF7-RI when compared with PNOF7.