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

TL;DR

ByteQC is a GPU-accelerated quantum chemistry software that enables large-scale simulations with significant speedups and expanded system sizes, utilizing advanced algorithms and quantum embedding techniques.

Contribution

This paper introduces ByteQC, a GPU-accelerated quantum chemistry package that significantly enhances computational efficiency and system size capabilities for large-scale simulations.

Findings

Up to 60× speedup over 100-core CPUs.

Expanded system sizes up to over 100,000 orbitals.

Successful demonstration of quantum embedding for large systems.

Abstract

Applying quantum chemistry algorithms to large-scale systems requires substantial computational resources scaled 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{\o}ller-Plesset perturbation theory, random phase approximation, coupled cluster methods, and quantum Monte Carlo 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$\times$ speedup when compared to 100-core CPUs. Additionally, the tractable system sizes have been significantly expanded: 1,610 orbitals for coupled cluster with single and double excitations (1,380 orbitals with perturbative triple excitations), 11,040 orbitals for M{\o}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 (2,752 orbitals) and a water monomer adsorbed on a boron nitride surface (3,929 orbitals), achieving the gold-standard accuracy.