Journal of Computational Chemistry

The numerical evaluation of Slater integrals on graphics processing units

Duy‐Khoi Dang, Leighton Wilson, Paul Martin Zimmerman

Published Jul 21, 2022• 12 citations

Abstract

Abstract 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.

DOI: 10.1002/jcc.26968

Related papers

Quantum Chemistry on Graphical Processing Units. 1. Strategies for Two-Electron Integral Evaluation

Ivan S. Ufimtsev, Todd J. Martı́nez

Modern videogames place increasing demands on the computational and graphical hardware, leading to novel architectures that have great potential in the context of high performance computing and molecular simulation. We demonstrate that Graphical Processing Units (GPUs) can be used very efficiently to calculate two-electron repulsion integrals over Gaussian basis functions [Formula: see text] the first step in most quantum chemistry calculations. A benchmark test performed for the evaluation of approximately 10(6) (ss|ss) integrals over contracted s-orbitals showed that a naïve algorithm implemented on the GPU achieves up to 130-fold speedup over a traditional CPU implementation on an AMD Opteron. Subsequent calculations of the Coulomb operator for a 256-atom DNA strand show that the GPU advantage is maintained for basis sets including higher angular momentum functions.

Journal of Chemical Theory and ComputationJan 25, 2008606 citations

Efficient GPU Implementation of the McMurchie–Davidson Method for Shell-Based ERI Computations

Haruto Fujii, Yasuaki Ito, Nobuya Yokogawa +5 more

Quantum chemistry offers the formal machinery to derive molecular and physical properties arising from (sub)atomic interactions. However, as molecules of practical interest are largely polyatomic, contemporary approximation schemes such as the Hartree–Fock scheme are computationally expensive due to the large number of electron repulsion integrals (ERIs). Central to the Hartree–Fock method is the efficient computation of ERIs over Gaussian functions (GTO-ERIs). Here, the well-known McMurchie–Davidson method (MD) offers an elegant formalism by incrementally extending Hermite Gaussian functions and auxiliary tabulated functions. Although the MD method offers a high degree of versatility to acceleration schemes through Graphics Processing Units (GPUs), the current GPU implementations limit the practical use of supported values of the azimuthal quantum number. In this paper, we propose a generalized framework capable of computing GTO-ERIs for arbitrary azimuthal quantum numbers, provided that the intermediate terms of the MD method can be stored. Our approach benefits from extending the MD recurrence relations through shells, batches, and triple-buffering of the shared memory, and ordering similar ERIs, thus enabling the effective parallelization and use of GPU resources. Furthermore, our approach proposes four GPU implementation schemes considering the suitable mappings between Gaussian basis and CUDA blocks and threads. Our computational experiments involving the GTO-ERI computations of molecules of interest on an NVIDIA A100 Tensor Core GPU (NVIDIA, Santa Clara, CA, USA) have revealed the merits of the proposed acceleration schemes in terms of computation time, including up to a 72× improvement over our previous GPU implementation and up to a 4500× speedup compared to a naive CPU implementation, highlighting the effectiveness of our method in accelerating ERI computations for both monatomic and polyatomic molecules. Our work has the potential to explore new parallelization schemes of distinct and complex computation paths involved in ERI computation.

Applied SciencesFeb 27, 20256 citations

Dynamic Schwarz Screening for GPU ‐Accelerated Fock Matrix Computation

Satoki Tsuji, Yasuaki Ito, Haruto Fujii +5 more

ABSTRACT In quantum chemistry, fundamental molecular orbital methods such as the Hartree–Fock theory are pragmatically constrained by the computational complexity behind two‐electron repulsion integrals (ERIs) encoding electron Coulomb interactions. Although the computational cost of ERIs roughly scales with the fourth power of the atom count in the target molecule, it can be reduced by screening techniques based on the Schwarz inequality. However, the Schwarz screening involves evaluating the necessity of computation for individual ERIs based on their upper bounds, which results in substantial overhead consuming most of the total computation time, especially for large‐sized molecules. To tackle this bottleneck, we propose a new dynamic screening algorithm and its efficient GPU implementation, which drastically reduces the number of upper bound evaluations. Our parallel algorithm can discard numerically insignificant ERIs without even evaluating their upper bounds by dynamic scheduling of thread blocks, allowing to traverse only the significant ERIs. The merits and performance of our GPU implementation has been evaluated through computational experiments using an NVIDIA A100/H200 GPU with a relevant set of low and medium molecular‐weight compounds. As a result, our dynamic screening reduces the Schwarz upper bound evaluations by up to and achieves a speedup of up to compared to examining the necessity of computation for all ERIs. Our work has the potential to further enhance the performance frontiers of GPU‐based large‐scale molecular quantum chemistry calculations.

Concurrency and Computation Practice and ExperienceApr 30, 20260 citations

Dynamic Screening of Two-Electron Repulsion Integrals in GPU Parallelization

Satoki Tsuji, Yasuaki Ito, Haruto Fujii +4 more

Two-electron repulsion integrals (ERIs) are among the most foundational quantities for numerically solving the Schrödinger equation, accounting for the Coulomb interactions between electrons in a molecule. Calculating ERIs is a computationally intensive task of determining integral values for every combination of four atomic orbitals, which dominates the execution time of fundamental frameworks in quantum chemistry such as the Hartree-Fock method. It is known that the numerous ERI calculations can be significantly reduced using the Schwarz inequality. However, the Schwarz screening requires evaluating the upper bound for every integral value, leading to the time complexity comparable to the entire ERI calculations themselves. This paper proposes a dynamic screening algorithm that minimizes the number of the upper bound evaluations using the asynchronous parallelism of GPUs. Our parallel algorithm performed by dynamically scheduled CUDA blocks can discard most ERI calculations without even evaluating their upper bounds. Furthermore, by screening ERIs at the level of the finest integral units, primitive integrals, we reduce more integral calculations over coarse-grained screening. Experimental results for various molecules using an NVIDIA A100 GPU demonstrate that the proposed dynamic screening achieves up to an 18.0-fold speedup compared to the non-screening ERI computation while keeping the energy error below 1.0 × 10−7hartree.

International Symposium on Computing and Networking Workshops (CANDARW)Nov 26, 20245 citations

Calculation of Quantum Chemical Two-Electron Integrals by Applying Compiler Technology on GPU

Gábor János Tornai, István Ladjánszki, Ádám Rák +2 more

In this article, we present an effective approach to calculate quantum chemical two-electron integrals over basis sets consisting of Gaussian-type basis functions on graphical processing unit (GPU). Our framework generates several different variants called routes to the same integral problem with different integral algorithms (McMurchie-Davidson, Head-Gordon-Pople, and Rys) and precision. Each route is benchmarked on more GPU architectures, and with this data, a model is fitted to select the best available route for an integral task given a GPU architecture. Moreover, this approach supports the computation of high angular momentum orbitals up to g effectively on GPU, tested up to cc-pVQZ-sized basis sets. Rigorous analysis is shown regarding the effectiveness of our method. Molecule simulations with several basis sets are measured using NVIDIA GTX 1080 Ti, NVIDIA P100, and NVIDIA V100 cards.

Journal of Chemical Theory and ComputationSep 10, 201955 citations

Generation of Kernels for Calculating Electron Repulsion Integrals of High Angular Momentum Functions on GPUs – Preliminary Results

A.V. Titov, Volodymyr Kindratenko, Ivan S. Ufimtsev +1 more

Abstract—Evaluation of electron repulsion integrals (ERIs) takes considerable time in modern quantum chemistry applications and also presents a certain difficulty to be efficiently computed on GPUs. Here, we describe a novel methodology for generating high-arithmetic-density kernels for ERI evaluation of d and higher angular momentum functions, as well as highlight challenges associated with the efficient use of GPUs for ERIs. Employing this approach allows to perform the ERI evaluations substantially faster than just using traditional CPU codes. Index Terms—GPUs; quantum chemistry; gaussian integrals; d-functions I.

Jan 1, 20116 citations