Numerically Exact Configuration Interaction at Quadrillion-Determinant Scale
arXiv:2505.20375 · doi:10.1038/s41467-025-65967-7
Abstract
The combinatorial scaling of configuration interaction (CI) has long restricted its applicability to only the simplest molecular systems. Here, we report the first numerically exact CI calculation exceeding one quadrillion () determinants, enabled by categorical compression within the small-tensor-product distributed active space (STP-DAS) framework. As a demonstration, we converged the relativistic complete active space CI (CASCI) ground state of HBrTe involving over complex-valued 2-spinor determinants in under 34.5 hours (time-to-completion) using 1000 nodes, representing the largest CASCI calculation reported to date. Additionally, we achieved -build times of just 5 minutes for systems with approximately 150 billion complex-valued 2-spinor determinants using only a few compute nodes. Extensive benchmarks confirm that the method retains numerical exactness with drastically reduced resource demands. Compared to previous state-of-the-art CI calculations, this work represents a 3-orders-of-magnitude increase in CI space, a 6-orders-of-magnitude increase in FLOP count, and a 6-orders-of-magnitude improvement in computational speed. By introducing a numerically exact, categorically compressed representation of the CI expansion vectors and reformulating the -build accordingly, we eliminate memory bottlenecks associated with storing excitation lists and CI vectors while significantly reducing computational cost. A compression-compatible preconditioner further enhances performance by generating compressed CI expansion vectors throughout Davidson iterations. This work establishes a new computational frontier for numerically exact CI methods, enabling chemically and physically accurate simulations of strongly correlated, spin-orbit coupled systems previously thought to be beyond reach.
References in corpus (22)
- Heat-bath Configuration Interaction: An efficient selected CI algorithm inspired by heat-bath sampling
- Semistochastic Heat-bath Configuration Interaction method: selected configuration interaction with semistochastic perturbation theory
- A deterministic alternative to the full configuration interaction quantum Monte Carlo method
- An adaptive configuration interaction approach for strongly correlated electrons with tunable accuracy
- An efficient implementation of two-component relativistic exact-decoupling methods for large molecules
- Pushing Configuration-Interaction to the Limit: Towards Massively Parallel MCSCF Calculations
- Projected Quasiparticle Theory for Molecular Electronic Structure
- Excited states using semistochastic heat-bath configuration interaction
- Nuclear energy gradients for internally contracted complete active space second-order perturbation theory: Multistate extensions
- Fast Semistochastic Heat-Bath Configuration Interaction
- Range-separated density-functional theory with random phase approximation applied to noncovalent intermolecular interactions
- The Chronus Quantum (ChronusQ) Software Package
- Four-Component Density Matrix Renormalization Group
- Relativistic Internally Contracted Multireference Electron Correlation Methods
- Large-scale relativistic complete active space self-consistent field with robust convergence
- Assessing the Accuracy of the Jastrow Antisymmetrized Geminal Power in the H4 Model System
- Convergence to the fixed-node limit in deep variational Monte Carlo
- Reduced density matrix functional theory from an ab initio seniority-zero wave function: Exact and approximate formulations along adiabatic connection paths
- Bounds for the Rayleigh quotient and the spectrum of self-adjoint operators
- A greedy algorithm for computing eigenvalues of a symmetric matrix
- A truncated Davidson method for the efficient "chemically accurate" calculation of full configuration interaction wavefunctions without any large matrix diagonalization
- Assessing many-body methods on the potential energy surface of the (H) hydrogen dimer