We provide a systematic evaluation of the sample-based quantum diagonalization (SQD) method for electronic structure based on the W4-11 thermochemistry dataset, comprising 124 total atomization, 83 bond dissociation, 20 isomerization, 505 heavy-atom transfer, and 13 nucleophilic substitution processes, covering diverse bonding situations and reaction mechanisms. This is the largest study assessing the accuracy and precision of a quantum-hybrid algorithm on a digital quantum device across a variety of molecular systems and chemical reactions, using 16.85 hours on the superconducting quantum processor ibm_rensselaer and 724.22 node hours on the supercomputer AiMOS. To ensure a fair comparison, our study employs commensurate resource allocation for both classical and quantum simulations. Although SQD exhibits large statistical deviations from ground-state reference energies, energy extrapolations yield CCSD-level accuracy. While bond-breaking reactions show a systematic improvement as computational resources increase, nucleophilic substitution or heavy atom transfer reactions do not. The limitations quantified in this manuscript indicate opportunities for improvement in SQD-based algorithms. This work provides a benchmark and community resource for exploring new quantum algorithms and devices, supported by an online benchmark challenge and an open-source Python library for direct comparison.
We provide a systematic evaluation of the sample-based quantum diagonalization (SQD) method for electronic structure based on the W4-11 thermochemistry dataset, comprising 124 total atomization, 83 bond dissociation, 20 isomerization, 505 heavy-atom transfer, and 13 nucleophilic substitution processes, covering diverse bonding situations and reaction mechanisms. This is the largest study assessing the accuracy and precision of a quantum-hybrid algorithm on a digital quantum device across a variety of molecular systems and chemical reactions, using 16.85 hours on the superconducting quantum processor ibm_rensselaer and 724.22 node hours on the supercomputer AiMOS. To ensure a fair comparison, our study employs commensurate resource allocation for both classical and quantum simulations. Although SQD exhibits large statistical deviations from ground-state reference energies, energy extrapolations yield CCSD-level accuracy. While bond-breaking reactions show a systematic improvement as computational resources increase, nucleophilic substitution or heavy atom transfer reactions do not. The limitations quantified in this manuscript indicate opportunities for improvement in SQD-based algorithms. This work provides a benchmark and community resource for exploring new quantum algorithms and devices, supported by an online benchmark challenge and an open-source Python library for direct comparison.
Quantum computers offer the potential to simulate chemical systems beyond the capabilities of classical computers. However, current quantum computing approaches to ab initio molecular dynamics (AIMD) and vibrational analysis largely rely on the Variational Quantum Eigensolver (VQE), whose circuit depth, iterative optimization, and measurement overhead hinder the practical simulation of molecular systems. Here, we extend Sample-based Quantum Diagonalization (SQD), a quantum sampling-based electronic-structure method, to AIMD and vibrational analysis on noisy quantum hardware. The resulting framework enables nuclear forces and Hessians to be derived from quantum-computed electronic energies, facilitating molecular dynamics and characterization of potential energy surfaces. We demonstrate the approach using diverse molecular processes, including H2 bond vibration, the proton-exchange reaction (H+ + H2), the bimolecular nucleophilic substitution (SN2) reaction (Cl- + CH3Cl), and proton-transfer dynamics between two water molecules involving transient formation of the Zundel cation (H2O···H+···OH2). Across these systems, SQD-AIMD trajectories and vibrational frequencies obtained with noisy quantum hardware/noise-models show excellent agreement with classical reference calculations and substantially improved accuracy relative to previously reported VQE-based simulations. Importantly, the SQD-based approach eliminates costly variational optimization and significantly reduces measurement requirements relative to VQE, enabling AIMD simulations spanning approximately 500 timesteps (∼200 fs), substantially longer than previously reported AIMD simulations on quantum hardware. These results establish SQD-AIMD and SQD-vibrational analysis as practical and scalable frameworks for simulating chemical reaction dynamics and elucidating reaction mechanisms on near-term quantum hardware, extending quantum computing applications beyond static electronic structure calculations.
Quantum and hybrid quantum-classical algorithms have the potential to improve the accuracy and scale of quantum chemistry and impact a broad range of applications such as drug discovery and catalysis. In particular, subspace methods such as sample-based quantum diagonalization (SQD) can already tackle classically challenging systems. However, practical algorithms that go beyond the computation of electronic ground and excited states are still scarce. Here, we introduce the variable-geometry sample-based quantum diagonalization (VGSQD) algorithm, which extends SQD to variable molecular geometries by reusing quantum samples from a reference geometry at displaced ones, thus significantly reducing quantum computational cost. Using circuits of up to 124 qubits on ibm-quebec, an IBM Heron r2 quantum processor, we determine the equilibrium geometry and energy of the weakly bound, strongly correlated beryllium dimer, as well as the water molecule. Finally, we extend VGSQD to chemical reactions and compute the minimum energy path and energy barrier of a hydrogen exchange reaction. The results are comparable with state-of-the-art classical methods such as CCSD(T) and DMRG. This work broadens the scope of quantum algorithms for chemistry, from energies to geometry optimization and reaction pathways, both essential for practical applications.
We present a large-scale experimental study of quantum-computing-based molecular simulation carried out on IQM's Sirius 24-qubit superconducting processor, utilizing up to 16 operational qubits. The work employs Sample-based Quantum Diagonalization (SQD) together with the Local Unitary Cluster Jastrow (LUCJ) ansatz to estimate ground-state energies for a set of benchmark molecules, including H2, LiH, BeH2, H2O, and NH3. In addition, we introduce a Linear-CNOT variant of the Unitary Coupled-Cluster Singles and Doubles (LCNot-UCCSD) ansatz within the SQD workflow, trading higher circuit depth for reduced classical preprocessing. A comparison between these ansätze is provided, clarifying their respective strengths, limitations, and suitability for near-term quantum hardware. We further explore potential energy landscapes through 1D scans for H2 and HeH+ using both STO-3G and 6-31G basis sets, and for LiH and BeH2 in STO-3G. Extending beyond this, we demonstrate the experimental construction of a full 2D potential energy surface for the water molecule on quantum hardware, mapped over a 32 × 32 grid in bond length and bond angle. To move beyond small benchmark systems, we combine SQD(LUCJ) with Density Matrix Embedding Theory (DMET) to compute active-space energies for a set of ligand-like molecules, as well as the pharmacologically relevant amantadine system. Across all studies, the majority of quantum-computed energies agree with reference FCI results, as well as with DMET-CASCI energies for embedded systems, to within chemical accuracy for the chosen basis sets. These results demonstrate the reliability of sample-based diagonalization approaches and underscore the potential of hybrid embedding strategies for extending quantum simulations to increasingly complex molecular systems, while also highlighting their practicality on current IQM quantum hardware.
The accurate treatment of electron correlation in extended molecular systems remains computationally challenging using classical electronic structure methods. Hybrid quantum-classical algorithms offer a potential route to overcome these limitations; however, their practical deployment on existing quantum computers requires strategies that both reduce problem size and mitigate hardware noise. In this work, we investigate ground-state energy calculations of ligand-like molecules using Sample-based Quantum Diagonalization (SQD) within the Density Matrix Embedding Theory (DMET) framework, focusing on low-symmetry systems with diverse bonding motifs that exhibit subsystem-dependent variations in fragment-environment entanglement. These entanglement-based variations directly influence bath orbital construction, impurity sizes, and the structure of the embedded Hamiltonians, posing nontrivial challenges for both embedding and quantum sampling. By combining DMET fragmentation with SQD-based construction of reduced configuration spaces through quantum sampling and iterative configuration recovery, we perform quantum simulations on IBM's Eagle R3 (IBM Sherbrooke) and IBM's Heron R3 (IBM Boston) superconducting quantum hardware thereby, showing that the entanglement structure across embedding subsystems plays a central role in determining the efficiency and accuracy of the simulations. Despite these complexities, we show that the DMET-SQD framework yields ground-state energies in strong agreement with DMET-FCI benchmarks, achieving chemical accuracy (1 kcal/mol) across all systems studied. These results demonstrate that SQD-based quantum simulations can be robustly extended to low-symmetry, chemically realistic, industry relevant molecules, and highlight the importance of entanglement-aware embedding strategies for scalable quantum electronic structure calculations.
Abstract The simulation of molecular electronic structure is an important application of quantum devices. Recently, it has been shown that quantum devices can be effectively combined with classical supercomputing centers in the context of the sample-based quantum diagonalization (SQD) algorithm. This allowed the largest electronic structure quantum simulation to date (77 qubits) and opened near-term devices to practical use cases in chemistry toward the hundred-qubit mark. However, the description of many important physical and chemical properties of those systems, such as photo-absorption/-emission, requires a treatment that goes beyond the ground state alone. In this work, we extend the SQD algorithm to determine low-lying molecular excited states. The extended-SQD method improves over the original SQD method in accuracy, at the cost of an additional computational step. It also improves over quantum subspace expansion based on single and double electronic excitations, a widespread approach to excited states on pre-fault-tolerant quantum devices, in both accuracy and efficiency. We employ the extended SQD method to compute the first singlet (S 1) and triplet (T 1) excited states of the nitrogen molecule with a correlation-consistent basis set, and the ground- and excited-state properties of the [2Fe-2S] cluster.
Accurately capturing electron correlation in large-scale molecular systems remains one of the foremost challenges in quantum chemistry and a primary driver for the development of quantum algorithms. Classical configuration-interaction methods, while rigorous, suffer from exponential scaling, rendering them impractical for large or strongly correlated systems. Overcoming this limitation is central to realizing the promise of quantum computing in chemistry. Here, we investigate the pyridine-Li ion complex using three quantum algorithms: the variational quantum eigensolver (VQE), the subspace quantum diagonalization (SQD) method, and the recently introduced handover iterative VQE (HI-VQE). Our results demonstrate how new generations of hybrid quantum-classical frameworks overcome the scalability and noise sensitivity that constrain conventional VQE approaches. SQD and HI-VQE achieve ground-state energy calculations for problem sizes inaccessible to classical computation, marking a clear advance toward quantum advantage. In particular, HI-VQE enables calculations within active spaces as large as (24e,22o), requiring 44 qubits-well beyond the reach of classical CASCI and VQE. This capability provides a systematic pathway for incorporating increasing numbers of electrons into quantum treatment, thereby approaching exact molecular energies. Importantly, both SQD and HI-VQE exhibit robustness against hardware noise, a critical improvement over earlier approaches. By enabling quantum simulations of molecular systems previously deemed intractable, SQD and HI-VQE offer a realistic route toward practical quantum advantage in computational chemistry. The comparison between HI-VQE and SQD shows that optimizing circuit parameters is crucial for accurate simulation.