Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

The Journal of Chemical Physics

The updates in Libcint 6: More integrals, API refinements, and SIMD optimization techniques

Qiming Sun

Published May 6, 2024• 15 citations
Publisher

Abstract

Libcint is a library designed for the evaluation of analytical integrals for Gaussian type orbitals. It prioritizes simplicity, ease of use, and efficiency for the development of quantum chemistry programs. In the release of version 6.0, Libcint supports the computation of integrals for various operators, such as overlap, Coulomb, Gaunt, Breit, attenuated Coulomb, Slater-type geminals, and Yukawa potential, as well as arbitrary orders of derivatives for these operators. To enhance the usability of the library, Libcint provides a uniform function signature for all integral functions. A code generator is included to automate the implementation of new integrals. To achieve better performance on modern central processing unit architectures, the library employs explicit single instruction multiple data parallelization in the code implementation.

DOI
10.1063/5.0200293

Related papers

Implementation of McMurchie–Davidson Algorithm for Gaussian AO Integrals Suited for SIMD Processors

Andrey Asadchev, Edward F. Valeev

High Resolution Image Download MS PowerPoint Slide We report an implementation of the McMurchie–Davidson evaluation scheme for 1- and 2-particle Gaussian AO integrals designed for processors with Single Instruction Multiple Data (SIMD) instruction sets. Like in our recent MD implementation for graphical processing units (GPUs) [ Asadchev, A.; Valeev, E. F. . J. Chem. Phys. 2024, 160, 244109 .], variable-sized batches of shellsets of integrals are evaluated at a time. By optimizing for the floating point instruction throughput rather than minimizing the number of operations, this approach achieves up to 50% of the theoretical hardware peak FP64 performance for many common SIMD-equipped platforms (AVX2, AVX512, NEON), which translates to speedups of up to 30 over the state-of-the-art one-shellset-at-a-time implementation of Obara–Saika-type schemes in Libint for a variety of primitive and contracted integrals. As with our previous work, we rely on the standard C++ programming language─such as the std::simd standard library feature to be included in the 2026 ISO C++ standard─without any explicit code generation to keep the code base small and portable. The implementation is part of the open source LibintX library freely available at https://github.com/ValeevGroup/libintx .

arXiv:2506.12501•The Journal of Physical Chemistry A•Oct 13, 2025•1 citation

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 Computation•Jan 25, 2008•606 citations

The numerical evaluation of Slater integrals on graphics processing units

Duy‐Khoi Dang, Leighton Wilson, Paul Martin Zimmerman

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.

Journal of Computational Chemistry•Jul 21, 2022•12 citations

Accelerating quantum chemistry with vectorized and batched integrals

Hua Huang, Edmond Chow

This paper presents the first quantum chemistry calculations using a recently developed vectorized library for computing electron repulsion integrals. To lengthen the SIMD loop and thus improve SIMD utilization, the approach used in this paper is to batch together multiple integrals that can be computed using the same code path. The standard approach is to compute integrals one at a time, and thus a batching procedure had to be developed. This paper shows proof-of-concept and demonstrates the performance gains possible when the batched approach is used. Batching also enables certain optimizations when the integrals are used to compute the Fock matrix. We further describe several other optimizations that were needed to obtain up to a 270% speedup over the non-batching version of the code, making a compelling case for adopting the presented techniques in quantum chemistry software.

IEEE International Conference on High Performance Computing, Data, and Analytics•Nov 11, 2018•6 citations

Accelerating Quantum Chemistry with Vectorized and Batched Integrals

Hua Huang, Edmond Chow

This paper presents the first quantum chemistry calculations using a recently developed vectorized library for computing electron repulsion integrals. To lengthen the SIMD loop and thus improve SIMD utilization, the approach used in this paper is to batch together multiple integrals that can be computed using the same code path. The standard approach is to compute integrals one at a time, and thus a batching procedure had to be developed. This paper shows proof-of-concept and demonstrates the performance gains possible when the batched approach is used. Batching also enables certain optimizations when the integrals are used to compute the Fock matrix. We further describe several other optimizations that were needed to obtain up to a 270% speedup over the non-batching version of the code, making a compelling case for adopting the presented techniques in quantum chemistry software.

International Conference for High Performance Computing, Networking, Storage and Analysis (SC)•Nov 1, 2018•13 citations

Designing quantum chemistry algorithms with just-in-time compilation

Xiaojie Wu, Qiming Sun, Yuanheng Wang

We introduce just-in-time (JIT) compilation to the integral kernels for Gaussian-type orbitals to enhance the efficiency of electron repulsion integral computations. For Coulomb and exchange (JK) matrices, JIT-based algorithms yield a 2× speedup for the small 6-31G* basis set over GPU4PySCF v1.4 on an NVIDIA A100-80G GPU. By incorporating a novel algorithm designed for orbitals with high angular momentum, the efficiency of JK evaluations with the large def2-TZVPP basis set is improved by up to 4×. The core CUDA implementation is compact, comprising only ∼1000 lines of code, including support for single-precision arithmetic. Furthermore, the single-precision implementation achieves a 3× speedup over the previous state-of-the-art.

arXiv:2507.09772•The Journal of Chemical Physics•May 15, 2026•2 citations
Libcint is a library designed for the evaluation of analytical integrals for Gaussian type orbitals. It prioritizes simplicity, ease of use, and efficiency for the development of quantum chemistry programs. In the release of version 6.0, Libcint supports the computation of integrals for various operators, such as overlap, Coulomb, Gaunt, Breit, attenuated Coulomb, Slater-type geminals, and Yukawa potential, as well as arbitrary orders of derivatives for these operators. To enhance the usability of the library, Libcint provides a uniform function signature for all integral functions. A code generator is included to automate the implementation of new integrals. To achieve better performance on modern central processing unit architectures, the library employs explicit single instruction multiple data parallelization in the code implementation.
High Resolution Image Download MS PowerPoint Slide We report an implementation of the McMurchie–Davidson evaluation scheme for 1- and 2-particle Gaussian AO integrals designed for processors with Single Instruction Multiple Data (SIMD) instruction sets. Like in our recent MD implementation for graphical processing units (GPUs) [ Asadchev, A.; Valeev, E. F. . J. Chem. Phys. 2024, 160, 244109 .], variable-sized batches of shellsets of integrals are evaluated at a time. By optimizing for the floating point instruction throughput rather than minimizing the number of operations, this approach achieves up to 50% of the theoretical hardware peak FP64 performance for many common SIMD-equipped platforms (AVX2, AVX512, NEON), which translates to speedups of up to 30 over the state-of-the-art one-shellset-at-a-time implementation of Obara–Saika-type schemes in Libint for a variety of primitive and contracted integrals. As with our previous work, we rely on the standard C++ programming language─such as the std::simd standard library feature to be included in the 2026 ISO C++ standard─without any explicit code generation to keep the code base small and portable. The implementation is part of the open source LibintX library freely available at https://github.com/ValeevGroup/libintx .
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.
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.
This paper presents the first quantum chemistry calculations using a recently developed vectorized library for computing electron repulsion integrals. To lengthen the SIMD loop and thus improve SIMD utilization, the approach used in this paper is to batch together multiple integrals that can be computed using the same code path. The standard approach is to compute integrals one at a time, and thus a batching procedure had to be developed. This paper shows proof-of-concept and demonstrates the performance gains possible when the batched approach is used. Batching also enables certain optimizations when the integrals are used to compute the Fock matrix. We further describe several other optimizations that were needed to obtain up to a 270% speedup over the non-batching version of the code, making a compelling case for adopting the presented techniques in quantum chemistry software.
This paper presents the first quantum chemistry calculations using a recently developed vectorized library for computing electron repulsion integrals. To lengthen the SIMD loop and thus improve SIMD utilization, the approach used in this paper is to batch together multiple integrals that can be computed using the same code path. The standard approach is to compute integrals one at a time, and thus a batching procedure had to be developed. This paper shows proof-of-concept and demonstrates the performance gains possible when the batched approach is used. Batching also enables certain optimizations when the integrals are used to compute the Fock matrix. We further describe several other optimizations that were needed to obtain up to a 270% speedup over the non-batching version of the code, making a compelling case for adopting the presented techniques in quantum chemistry software.
We introduce just-in-time (JIT) compilation to the integral kernels for Gaussian-type orbitals to enhance the efficiency of electron repulsion integral computations. For Coulomb and exchange (JK) matrices, JIT-based algorithms yield a 2× speedup for the small 6-31G* basis set over GPU4PySCF v1.4 on an NVIDIA A100-80G GPU. By incorporating a novel algorithm designed for orbitals with high angular momentum, the efficiency of JK evaluations with the large def2-TZVPP basis set is improved by up to 4×. The core CUDA implementation is compact, comprising only ∼1000 lines of code, including support for single-precision arithmetic. Furthermore, the single-precision implementation achieves a 3× speedup over the previous state-of-the-art.