Full configuration interaction (FCI) provides exact electronic structure within a given atomic basis, but its computational cost grows exponentially with the number of spin orbitals. Selected configuration interaction (SCI) methods alleviate this limitation by retaining only the most important Slater determinants. However, the repeated identification of important determinants remains a major computational bottleneck. We present a quantum assisted selected configuration interaction (QASCI) method that combines SCI with graph based block diagonalization (GBBD) of FCI Hamiltonian. The GBBD method partitions FCI Hamiltonian into independent blocks, within which determinant selection problem is formulated as a quadratic unconstrained binary optimization (QUBO) problem. The QUBO problems for selecting determinants to construct SCI space are iteratively solved using a fast annealing approach. We benchmark method on H8-H18 hydrogen chains and Li2S in STO3G basis. For Hn chains, chemical accuracy is achieved while retaining only a small fraction of Slater determinants, and this fraction decreases with increasing n, despite the exponential growth of the FCI Hilbert space. For Li2S, QASCI results remain within chemical accuracy while retaining substantially fewer determinants than the full FCI space. We apply QASCI to N2 using the 631G basis, considering both active orbital and full orbital treatments. The full orbital QASCI calculation, using 50000 determinants, yields a lower ground state energy than an FCI calculation within an active space comprising 12 spin orbitals and 12 electrons. These results demonstrate that the combination of QASCI and the GBBD approach can substantially reduce computational cost of determinant selection while maintaining the accuracy of FCI based electronic structure calculations, thereby enabling accurate calculations in larger orbital spaces.
Full configuration interaction (FCI) provides exact electronic structure within a given atomic basis, but its computational cost grows exponentially with the number of spin orbitals. Selected configuration interaction (SCI) methods alleviate this limitation by retaining only the most important Slater determinants. However, the repeated identification of important determinants remains a major computational bottleneck. We present a quantum assisted selected configuration interaction (QASCI) method that combines SCI with graph based block diagonalization (GBBD) of FCI Hamiltonian. The GBBD method partitions FCI Hamiltonian into independent blocks, within which determinant selection problem is formulated as a quadratic unconstrained binary optimization (QUBO) problem. The QUBO problems for selecting determinants to construct SCI space are iteratively solved using a fast annealing approach. We benchmark method on H8-H18 hydrogen chains and Li2S in STO3G basis. For Hn chains, chemical accuracy is achieved while retaining only a small fraction of Slater determinants, and this fraction decreases with increasing n, despite the exponential growth of the FCI Hilbert space. For Li2S, QASCI results remain within chemical accuracy while retaining substantially fewer determinants than the full FCI space. We apply QASCI to N2 using the 631G basis, considering both active orbital and full orbital treatments. The full orbital QASCI calculation, using 50000 determinants, yields a lower ground state energy than an FCI calculation within an active space comprising 12 spin orbitals and 12 electrons. These results demonstrate that the combination of QASCI and the GBBD approach can substantially reduce computational cost of determinant selection while maintaining the accuracy of FCI based electronic structure calculations, thereby enabling accurate calculations in larger orbital spaces.
Accurate solutions to the electronic Schrödinger equation remain a central challenge in quantum chemistry due to the exponential growth of the many-electron Hilbert space. Selected configuration interaction (SCI) methods mitigate this challenge by restricting the wavefunction expansion to a subset of important determinants, but their efficiency is limited by the quality/cost of determinant selection and Hamiltonian diagonalization. In this work, we propose QUBO-SCI, a novel approach that reformulates the selection step in SCI as a quadratic unconstrained binary optimization (QUBO) problem. This formulation enables the sampling of low-energy electronic configurations as bitstrings using both quantum-inspired approaches and quantum annealing hardware. The sampled configurations are then incorporated into an iterative SCI framework using a perturbative importance criterion. We demonstrate the method on small benchmark systems, including H 2 , C 2 , H 2 O, and N 2 , showing that QUBO-SCI achieves near-FCI accuracy with a significantly reduced number of determinants. The results highlight the potential of annealing-based sampling as an efficient alternative to both conventional determinant selection strategies and circuit-based quantum SCI approaches in electronic structure calculations.
Recently, a new distributed implementation of the full configuration interaction (FCI) method has been reported [Gao et al. J. Chem Theory Comput. 2024, 20, 1185]. Thanks to a hybrid parallelization scheme, the authors were able to compute the exact energy of propane (\ce{C3H8}) in the minimal basis STO-3G. This formidable task involves handling an active space of 26 electrons in 23 orbitals or a Hilbert space of \SI{1.3d12} determinants. This is, by far, the largest FCI calculation reported to date. Here, we illustrate how, from a general point of view, selected configuration interaction (SCI) can achieve microhartree accuracy at a fraction of the computational and memory cost, via a sparse exploration of the FCI space. The present SCI calculations are performed with the \textit{Configuration Interaction using a Perturbative Selection made Iteratively} (CIPSI) algorithm, as implemented in a determinant-driven way in the \textsc{quantum package} software. The present study reinforces the common wisdom that among the exponentially large number of determinants in the FCI space, only a tiny fraction of them significantly contribute to the energy. More importantly, it demonstrates the feasibility of achieving comparable accuracy using more reasonable and sustainable computational resources, hence reducing the ever-growing carbon footprint of computational chemistry.
Abstract Finding the ground-state energy of molecules is an important and challenging computational problem for which quantum computing can potentially find efficient solutions. The variational quantum eigensolver (VQE) is a quantum algorithm that tackles the molecular groundstate problem and is regarded as one of the flagships of quantum computing. Yet, to date, only very small molecules were computed via VQE, due to high noise levels in current quantum devices. Here we present an alternative variational quantum scheme that requires significantly less qubits than VQE. The reduction in the qubit number allows for shallower circuits to be sufficient, rendering the method more resistant to noise. The proposed algorithm, termed variational quantum selected-configuration-interaction (VQ-SCI), is based on: (a) representing the target groundstate as a superposition of Slater determinant configurations, encoded directly upon the quantum computational basis states; and (b) selecting a-priory only the most dominant configurations. This is demonstrated through a set of groundstate calculations of the H2, LiH, BeH2, H2O, NH3 and C2H4 molecules in the sto-3g basis set, performed on IBM quantum devices. We show that the VQ-SCI reaches the full configuration interaction energy within chemical accuracy using the lowest number of qubits reported to date. Moreover, when the SCI matrix is generated ‘on the fly’, the VQ-SCI requires exponentially less memory than classical SCI methods. This offers a potential remedy to a severe memory bottleneck problem in classical SCI calculations. Finally, the proposed scheme is general and can be straightforwardly applied for finding the groundstate of any Hermitian matrix, outside the chemical context.
Selected configuration interaction methods achieve near-exact electronic structure calculations by iteratively constructing compact variational spaces, but their efficiency depends critically on the heuristics used to identify important determinants. Here, we introduce a data-driven selection framework that recasts determinant importance as a binary classification task and integrates a neural-network classifier into the iterative CI workflow through an active-learning loop. At each iteration, a random subset of candidate determinants is labelled via temporary diagonalisation, and the trained classifier guides selection of the remaining configurations. We demonstrate the utility of this framework for both classical and quantum CI methods by calculating the ground-state energy of a diatomic molecule. Our method achieves result parity with traditional configuration interaction methods at substantially lower computational cost: roughly a ×5 reduction in memory and per-iteration cost for the classical cHCI variant, and convergence in markedly fewer iterations for the quantum-classical cSQD variant. These results establish classifier-assisted determinant selection as a lightweight, method-agnostic tool for compressing variational spaces and accelerating both classical and hybrid quantum-classical configuration interaction algorithms.
Selected configuration interaction (SCI) methods are currently enjoying a resurgence due to several recent developments which improve either the overall computational efficiency or the compactness of the resulting SCI vector. These recent advances have made it possible to get full CI (FCI) quality results for much larger orbital active spaces compared to conventional approaches. However, due to the starting assumption that the FCI vector has only a small number of significant Slater determinants, SCI becomes intractable for systems with strong correlation. This paper introduces a method for developing SCI algorithms in a way which exploits local molecular structure to significantly reduce the number of SCI variables. The proposed method is defined by first grouping the orbitals into clusters over which we can define many-particle cluster states. We then directly perform the SCI algorithm in a basis of tensor products of cluster states instead of Slater determinants. While the approach is general for arbitrarily defined cluster states, we find significantly improved performance by defining cluster states through a Tucker decomposition of the global (and sparse) SCI vector. To demonstrate the potential of this method, called tensor product selected configuration interaction (TPSCI), we present numerical results for a diverse set of examples: (1) modified Hubbard model with different inter- and intracluster hopping terms, (2) less obviously clusterable cases of bond breaking in N 2 and F 2, and (3) ground state energies of large planar π-conjugated systems with active spaces of up to 42 electrons in 42 orbitals. These numerical results show that TPSCI can be used to significantly reduce the number of SCI variables in the variational space, thus paving a path for extending these deterministic and variational SCI approaches to a wider range of physical systems.
High Resolution Image Download MS PowerPoint Slide Strongly correlated molecules remain out of reach for most electronic-structure solvers because the exact wave function spans a determinant space that grows exponentially with the number of orbitals. State-of-the-art selected-CI (e.g., CIPSI and HCI) and ML-CI methods mitigate search by exploiting Hamiltonian sparsity and learned prescreening. Nonetheless, in strongly multireference regimes, the excitation-based candidate lists can grow rapidly, which may increase memory and screening costs. We present the Hamiltonian-guided autoregressive selected-configuration interaction (HAAR-SCI), a learn–sample–compress workflow that runs on a single GPU. A gated Transformer samples determinants autoregressively; Gumbel Top-K noise encourages exploration, and GPU min-heap kernels keep only configurations with the largest Hamiltonian couplings. The network is retrained after each expansion, and iterations stop when successive energies differ by ≤1 mHa. Across an 18-molecule benchmark set reaching 116 spin–orbitals, HAAR-SCI attains a mean absolute error of 0.51 mHa while using much fewer determinants than heat-bath CI. It traces the entire N 2 dissociation curve within 0.67 mHa and achieves HCI accuracy for the 40-spin-orbital [Fe 2 S 2 (SCH 3 ) 4 ] 2– cluster with a determinant reduction of 72%, demonstrating its power on systems considered intractable for conventional selected-CI solvers. A final probability-mass pruning compresses typical wave functions by a further 10–50×, retaining <0.01% of the Hilbert space yet still capturing >99.9% of the correlation energy. HAAR-SCI thus offers a compact and truly scalable route to chemical-accuracy quantum chemistry on commodity hardware.