Memory-Efficient Recursive Evaluation of 3-Center Gaussian Integrals

Abstract

To improve the efficiency of Gaussian integral evaluation on modern accelerated architectures FLOP-efficient Obara-Saika-based recursive evaluation schemes are optimized for the memory footprint. For the 3-center 2-particle integrals that are key for the evaluation of Coulomb and other 2-particle interactions in the density-fitting approximation the use of multi-quantal recurrences (in which multiple quanta are created or transferred at once) is shown to produce significant memory savings. Other innovation include leveraging register memory for reduced memory footprint and direct compile-time generation of optimized kernels (instead of custom code generation) with compile-time features of modern C++/CUDA. Performance of conventional and CUDA-based implementations of the proposed schemes is illustrated for both the individual batches of integrals involving up to Gaussians with low and high angular momenta (up to $L=6$) and contraction degrees, as well as for the density-fitting-based evaluation of the Coulomb potential. The computer implementation is available in the open-source LibintX library.

37 pages, 2 figures, 6 tables

References in corpus (1)

• Fast evaluation of solid harmonic Gaussian integrals for local resolution-of-the-identity methods and range-separated hybrid functionals

Cited by in corpus (6)

• A call to arms: making the case for more reusable libraries
• High-performance evaluation of high angular momentum 4-center Gaussian integrals on modern accelerated processors
• Distributed Memory, GPU Accelerated Fock Construction for Hybrid, Gaussian Basis Density Functional Theory
• 3-center and 4-center 2-particle Gaussian AO integrals on modern accelerated processors
• Benchmarking the accuracy of the separable resolution of the identity approach for correlated methods in the numeric atom-centered orbitals framework
• Implementation of McMurchie-Davidson algorithm for Gaussian AO integrals suited for SIMD processors