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

arXiv:2502.17963 · Wiley Interdisciplinary Reviews Computational Molecular Science

ByteQC : GPU ‐Accelerated Quantum Chemistry Package for Large‐Scale Systems

Zhen Guo, Zigeng Huang, Qiaorui Chen, Jiang Shao, Guangcheng Liu, Hung Q. Pham, Yifei Huang, Changsu Cao, Ji Chen, Dingshun Lv

Published May 1, 2025• 7 citations
arXivHTMLPDF

Abstract

ABSTRACT Applying quantum chemistry algorithms to large‐scale systems requires substantial computational resources scaling with the system size and the desired accuracy. To address this, ByteQC, a fully functional and efficient package for large‐scale quantum chemistry simulations, has been open sourced at https://github.com/bytedance/byteqc , leveraging recent advances in computational power and many‐body algorithms. Regarding computational power, several standard algorithms are efficiently implemented on modern GPUs, ranging from mean‐field calculations (Hartree‐Fock and density functional theory) to post‐Hartree‐Fock methods such as Møller‐Plesset perturbation theory and coupled cluster methods. For the algorithmic approach, we also employ a quantum embedding method, which significantly expands the tractable system size while preserving high accuracy at the gold‐standard level. All these features have been systematically benchmarked. For standalone algorithms, the benchmark results demonstrate up to a 60× speedup when compared to 100‐core CPUs. Additionally, the tractable system sizes have been significantly expanded: 1610 orbitals for coupled cluster with single and double excitations (1380 orbitals with perturbative triple excitations), 11,040 orbitals for Møller‐Plesset perturbation theory of second order, 37,120 orbitals for mean‐field calculations under open boundary conditions, and over 100,000 orbitals for periodic boundary conditions. For the advanced quantum embedding feature, two representative examples are demonstrated: the water cluster problem (2752 orbitals) and a water monomer adsorbed on a boron nitride surface (3929 orbitals), achieving the gold‐standard accuracy. With these efforts, ByteQC is expected to significantly advance research in quantum chemistry, particularly in large‐scale, high‐accuracy calculations.

DOI
10.1002/wcms.70034

Related papers

macroqc 2.0: An electronic structure theory software for large-scale applications

Uğur Bozkaya, Betül Ermiş, Yavuz Alagöz +2 more

macroqc is a quantum chemistry software package designed for high-accuracy electronic structure computations and large-scale molecular simulations. The macroqc framework provides efficient implementations of energies and analytic gradients for a wide range of many-body perturbation theory and coupled-cluster (CC) methods. In particular, the analytic gradient implementations of second-order Møller-Plesset perturbation theory (MP2), coupled-cluster singles and doubles (CCSD), and CCSD with perturbative triples [CCSD(T)] exhibit excellent computational performance, even in comparison with established commercial quantum chemistry programs. In addition to conventional correlated methods, macroqc incorporates several advanced capabilities, including molecular fragmentation techniques for ground and excited state molecular systems and density-fitted analytic gradients. The software also supports interoperability with selected external packages through a flexible plugin framework, allowing developers to construct customized interfaces and extensions. The notable features of macroqc are the linear-scaling systematic molecular fragmentation (LSSMF) and the cluster-in-molecule (CIM) approaches, which enable the application of high-level post-Hartree-Fock approaches to large molecular systems that would otherwise be computationally prohibitive. Furthermore, the new machine learning (ML) library enables the integration of data-driven approaches into quantum chemical workflows, providing efficient tools for developing and applying ML-based models for molecular property prediction and electronic structure computations. With its combination of efficient implementations, advanced theoretical methods, and extensible design, macroqc is expected to serve as a useful platform for a broad range of computational chemistry applications.

The Journal of Chemical Physics•Oct 22, 2026•0 citations

Introducing GPU Acceleration into the Python-Based Simulations of Chemistry Framework

Rui Li, Qiming Sun, Xing Zhang +1 more

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.

arXiv:2407.09700•The Journal of Physical Chemistry A•Jan 23, 2025•42 citations

GPU ‐Accelerated One‐Electron Integral Computation for Quantum Chemistry

Nobuya Yokogawa, Yasuaki Ito, Satoki Tsuji +5 more

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.

Concurrency and Computation Practice and Experience•Feb 20, 2026•0 citations

A GPU Implementation of the Second- and Third-Order Møller-Plesset Perturbation Theory

Satoki Tsuji, Yasuaki Ito, Nobuya Yokogawa +4 more

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.

International Symposium on Computing and Networking Workshops (CANDARW)•Nov 25, 2025•1 citation

Parallel iQCC Enables 200 Qubit Scale Quantum Chemistry on Accelerated Computing Platforms Surpassing Classical Benchmarks in Ruthenium Catalysts

Seyyed Mehdi Hosseini Jenab, Brandon Henderson, Scott N. Genin

We introduce a parallel, GPU-accelerated implementation of the iterative qubit coupled cluster (iQCC) method that overcomes the exponential growth of the transformed Hamiltonian -- the principal bottleneck for classical emulation of quantum chemistry circuits. By distributing Hamiltonian terms across compute nodes via bit-wise partitioning and offloading Pauli contractions to GPUs, we achieve speedups exceeding two orders of magnitude over the serial CPU approach. Crucially, iQCC confines the variational evolution to a classically simulable operator subspace by selecting entanglers exclusively from the Direct Interaction Space, which guarantees non-vanishing energy gradients at every iteration and thereby naturally avoids the barren-plateau phenomenon that renders highly expressive quantum circuits untrainable. Leveraging these algorithmic and hardware advances, we simulate electronic-structure Hamiltonians for industrially relevant ruthenium catalysts in the 100--124 qubit regime, completing full ground-state calculations on NVIDIA GPUs in the ranges of 1.2 - 45 hrs and surpassing the accuracy of Density Matrix Renormalization Group. These results effectively de-quantize a significant portion of the NISQ roadmap: quantum advantage for chemistry is often assumed to emerge beyond ${\sim}50$ qubits, yet our work demonstrates that this frontier lies significantly further -- potentially past 200 qubits -- reshaping expectations for where genuine quantum advantage may first appear.

quant-pharXiv:2603.08883•Mar 11, 2026•0 citations

GANSU: A GPU-Native Quantum Chemistry Framework for Efficient Hartree–Fock and Post-HF Calculations

Yasuaki Ito, Satoki Tsuji, Koji Nakano +1 more

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.

Eng—Advances in Engineering•Apr 28, 2026•1 citation
ABSTRACT Applying quantum chemistry algorithms to large‐scale systems requires substantial computational resources scaling with the system size and the desired accuracy. To address this, ByteQC, a fully functional and efficient package for large‐scale quantum chemistry simulations, has been open sourced at https://github.com/bytedance/byteqc , leveraging recent advances in computational power and many‐body algorithms. Regarding computational power, several standard algorithms are efficiently implemented on modern GPUs, ranging from mean‐field calculations (Hartree‐Fock and density functional theory) to post‐Hartree‐Fock methods such as Møller‐Plesset perturbation theory and coupled cluster methods. For the algorithmic approach, we also employ a quantum embedding method, which significantly expands the tractable system size while preserving high accuracy at the gold‐standard level. All these features have been systematically benchmarked. For standalone algorithms, the benchmark results demonstrate up to a 60× speedup when compared to 100‐core CPUs. Additionally, the tractable system sizes have been significantly expanded: 1610 orbitals for coupled cluster with single and double excitations (1380 orbitals with perturbative triple excitations), 11,040 orbitals for Møller‐Plesset perturbation theory of second order, 37,120 orbitals for mean‐field calculations under open boundary conditions, and over 100,000 orbitals for periodic boundary conditions. For the advanced quantum embedding feature, two representative examples are demonstrated: the water cluster problem (2752 orbitals) and a water monomer adsorbed on a boron nitride surface (3929 orbitals), achieving the gold‐standard accuracy. With these efforts, ByteQC is expected to significantly advance research in quantum chemistry, particularly in large‐scale, high‐accuracy calculations.
macroqc is a quantum chemistry software package designed for high-accuracy electronic structure computations and large-scale molecular simulations. The macroqc framework provides efficient implementations of energies and analytic gradients for a wide range of many-body perturbation theory and coupled-cluster (CC) methods. In particular, the analytic gradient implementations of second-order Møller-Plesset perturbation theory (MP2), coupled-cluster singles and doubles (CCSD), and CCSD with perturbative triples [CCSD(T)] exhibit excellent computational performance, even in comparison with established commercial quantum chemistry programs. In addition to conventional correlated methods, macroqc incorporates several advanced capabilities, including molecular fragmentation techniques for ground and excited state molecular systems and density-fitted analytic gradients. The software also supports interoperability with selected external packages through a flexible plugin framework, allowing developers to construct customized interfaces and extensions. The notable features of macroqc are the linear-scaling systematic molecular fragmentation (LSSMF) and the cluster-in-molecule (CIM) approaches, which enable the application of high-level post-Hartree-Fock approaches to large molecular systems that would otherwise be computationally prohibitive. Furthermore, the new machine learning (ML) library enables the integration of data-driven approaches into quantum chemical workflows, providing efficient tools for developing and applying ML-based models for molecular property prediction and electronic structure computations. With its combination of efficient implementations, advanced theoretical methods, and extensible design, macroqc is expected to serve as a useful platform for a broad range of computational chemistry applications.
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.
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.
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.
We introduce a parallel, GPU-accelerated implementation of the iterative qubit coupled cluster (iQCC) method that overcomes the exponential growth of the transformed Hamiltonian -- the principal bottleneck for classical emulation of quantum chemistry circuits. By distributing Hamiltonian terms across compute nodes via bit-wise partitioning and offloading Pauli contractions to GPUs, we achieve speedups exceeding two orders of magnitude over the serial CPU approach. Crucially, iQCC confines the variational evolution to a classically simulable operator subspace by selecting entanglers exclusively from the Direct Interaction Space, which guarantees non-vanishing energy gradients at every iteration and thereby naturally avoids the barren-plateau phenomenon that renders highly expressive quantum circuits untrainable. Leveraging these algorithmic and hardware advances, we simulate electronic-structure Hamiltonians for industrially relevant ruthenium catalysts in the 100--124 qubit regime, completing full ground-state calculations on NVIDIA GPUs in the ranges of 1.2 - 45 hrs and surpassing the accuracy of Density Matrix Renormalization Group. These results effectively de-quantize a significant portion of the NISQ roadmap: quantum advantage for chemistry is often assumed to emerge beyond ∼50{\sim}50∼50 qubits, yet our work demonstrates that this frontier lies significantly further -- potentially past 200 qubits -- reshaping expectations for where genuine quantum advantage may first appear.
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.