We present a GPU-accelerated software package for the evaluation of analytic two-electron energy and gradient contributions in Hartree-Fock (HF) and Density Functional theory (DFT) calculations. The implementation is provided as a Python library with a C++ backend, enabling straightforward integration into modern computational chemistry and drug-discovery workflows. The code supports single-point energy and nuclear gradient evaluations on both single- and multi-GPU systems, and employs MPI-based parallelization with dynamic load balancing in multi-node environments. We report comprehensive benchmarks demonstrating favorable scaling with respect to system size, as well as high throughput for batched evaluations relevant to molecular dynamics, geometry optimization, and large-scale virtual screening. Parallel execution of a single system was carried out on up to 24 A100 GPUs. The implementation builds on optimized GPU-enabled variants of the LibintX and GauXC libraries to efficiently compute density-fitted Coulomb, semi-numerical exchange, and exchange--correlation contributions.
We present a GPU-accelerated software package for the evaluation of analytic two-electron energy and gradient contributions in Hartree-Fock (HF) and Density Functional theory (DFT) calculations. The implementation is provided as a Python library with a C++ backend, enabling straightforward integration into modern computational chemistry and drug-discovery workflows. The code supports single-point energy and nuclear gradient evaluations on both single- and multi-GPU systems, and employs MPI-based parallelization with dynamic load balancing in multi-node environments. We report comprehensive benchmarks demonstrating favorable scaling with respect to system size, as well as high throughput for batched evaluations relevant to molecular dynamics, geometry optimization, and large-scale virtual screening. Parallel execution of a single system was carried out on up to 24 A100 GPUs. The implementation builds on optimized GPU-enabled variants of the LibintX and GauXC libraries to efficiently compute density-fitted Coulomb, semi-numerical exchange, and exchange--correlation contributions.
We present a GPU-accelerated software package for the evaluation of analytic two-electron energy and gradient contributions in Hartree-Fock (HF) and Density Functional theory (DFT) calculations. The implementation is provided as a Python library with a C++ backend, enabling straightforward integration into modern computational chemistry and drug-discovery workflows. The code supports single-point energy and nuclear gradient evaluations on both single- and multi-GPU systems, and employs MPI-based parallelization with dynamic load balancing in multi-node environments. We report comprehensive benchmarks demonstrating favorable scaling with respect to system size, as well as high throughput for batched evaluations relevant to molecular dynamics, geometry optimization, and large-scale virtual screening. Parallel execution of a single system was carried out on up to 24 A100 GPUs. The implementation builds on optimized GPU-enabled variants of the LibintX and GauXC libraries to efficiently compute density-fitted Coulomb, semi-numerical exchange, and exchange--correlation contributions.
High Resolution Image Download MS PowerPoint Slide We introduce the first version of GPU4PySCF, a module that provides GPU acceleration of methods in PySCF . As a core functionality, this provides a GPU implementation of two-electron repulsion integrals (ERIs) for contracted basis sets comprising up to g functions using the Rys quadrature. As an illustration of how this can accelerate a quantum chemistry workflow, we describe how to use the ERIs efficiently in the integral-direct Hartree–Fock build and nuclear gradient construction. Benchmark calculations show a significant speedup of 2 orders of magnitude with respect to the multithreaded CPU Hartree–Fock code of PySCF and the performance comparable to other open-source GPU-accelerated quantum chemical packages, including GAMESS and QUICK, on a single NVIDIA A100 GPU.
We introduce a GPU-accelerated multigrid Gaussian-Plane-Wave density fitting (FFTDF) approach for efficient Fock builds and nuclear gradient evaluations within Kohn-Sham density functional theory, as implemented in the GPU4PySCF module of PySCF. Our CUDA kernels employ a grid-based parallelization strategy for contracting Gaussian basis function pairs and achieve up to 80% of the FP64 peak performance on NVIDIA GPUs, with no loss of efficiency for high angular momentum (up to f-shell) functions. Benchmark calculations on molecules and solids with up to 1536 atoms and 20480 basis functions show up to 25x speedup on an H100 GPU relative to the CPU implementation on a 28-core shared memory node. For a 256-water cluster, the ground-state energy and nuclear gradients can be computed in ~30 seconds on a single H100 GPU. This implementation serves as an open-source foundation for many applications, such as ab initio molecular dynamics and high-throughput calculations.
ABSTRACT In Quantum chemical computation, numerical schemes such as the Hartree–Fock (HF) and density functional theory (DFT) are widely used to solve the Schrödinger equation numerically, to realize experiment‐free prediction and analysis of key molecular properties such as structure and energy. Computing one‐electron integrals, such as kinetic energy integrals and nuclear attraction integrals, is essential in both HF and DFT to characterize the molecular electronic states. However, as molecules of practical interest grow in size and angular momentum, computing one‐electron orbitals becomes computationally expensive in most cases. Although computing kinetic energy integrals on CPUs is straightforward, bottlenecks in CPU‐GPU data transfer have often been overlooked. In this study, we propose an efficient method to compute both the kinetic‐energy and nuclear‐attractive integrals on GPUs. First, we explicitly and symbolically expand recurrence relations based on the Obara–Saika and McMurchie–Davidson methods to eliminate redundant operations, thus improving computational efficiency. Second, we implemented a hybrid method that selects the best/fastest of both methods depending on the integration task. Third, we achieved further speedups by using CUDA streams to parallelize the execution of multiple kernels and efficiently utilize multiprocessor resources on the GPU. Computational experiments using NVIDIA A100 GPUs and Intel Xeon Gold 6338 CPU on relevant molecules of interest demonstrated the superiority of our one‐electron integral GPU implementations, achieving a speedup of 20.2 times over PySCF, and a speedup of 132.6 times over GPU4PySCF.
GPU-accelerated quantum chemistry programs can dramatically reduce the time required for electronic structure calculations, yet most existing implementations either retrofit GPU kernels onto legacy CPU codebases or optimize individual kernels without addressing workflow-level integration overhead. We present GANSU (GPU Accelerated Numerical Simulation Utility), an open-source quantum chemistry framework written entirely in CUDA/C++ that integrates GPU-accelerated kernels for electron repulsion integrals, Fock matrix construction, and post-Hartree–Fock (post-HF) methods into a unified, GPU-resident execution pipeline. The key design principle is to eliminate host–device data transfers between computational stages by keeping all intermediate data, including density matrices, integral buffers, and Fock matrix replicas, on the GPU throughout the self-consistent field (SCF) iteration, combined with runtime-selectable integral strategies (stored ERI, resolution-of-the-identity, and Direct-SCF) that adapt to system size and available memory. On an NVIDIA H200 GPU, GANSU achieves end-to-end speedups of up to 52× over PySCF for SCF, 45× for MP2 on molecules with up to 470 basis functions, and 44× for FCI, while outperforming GPU4PySCF by up to 34× for FCI, across a range of molecular systems with up to 650 basis functions. The framework further provides analytical energy gradients and geometry optimization with nine algorithms, all operating within the same GPU-resident data flow. These results demonstrate that workflow-aware kernel integration, not just kernel-level optimization, is essential for realizing the full potential of GPU acceleration in scientific computing. GANSU is publicly available under the BSD-3-Clause license.
Post-Hartree-Fock methods are wavefunction-based quantum chemistry approaches that explicitly account for the electron correlation beyond the mean-field approximation, allowing accurate determination of electron orbitals and energies in a given system. However, their computational cost scales as the fourth power or higher with respect to the number of atoms, severely limiting their practical applications. This paper proposes an efficient GPU implementation of the Møller-Plesset method, a post-Hartree-Fock approach based on perturbation theory. We focus on the computation of second- and thirdorder perturbation energies (MP2 and MP3) for closed-shell systems, which incorporate electron correlation on top of the Hartree-Fock energy. To accelerate the extensive contraction operations arising in the Møller-Plesset method, our implementation employs a formulation that maximizes the use of highly optimized matrix multiplication routines. Furthermore, we effectively utilize computational capabilities of GPU architecture through a hierarchical reduction using warp-level primitives and shared memory, as well as contiguous coalesced access to matrix elements. Performance evaluation on an NVIDIA A100 GPU shows up to 255.5 and 71.1 times speedups for the MP2 correlation energy computation, compared with the widely used quantum chemistry software package PySCF and its GPU plugin GPU4PySCF, respectively. In the MP3 calculations, our DGEMM-based implementation achieves a maximum speedup of 96.2 times over PySCF.