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:2604.01983

Towards Chemically Accurate and Scalable Quantum Simulations on IQM Quantum Hardware: A Quantum-HPC Hybrid Approach

Anurag K. S. V., Ashish Kumar Patra, Manas Mukherjee, Alok Shukla, Sai Shankar P., Ruchika Bhat, Radhika T. S. L., Jaiganesh G

quant-phcs.ETphysics.chem-phphysics.comp-phSubmitted Apr 2, 2026• Announced Apr 3, 2026• 0 citations
v1 · Apr 2, 2026v2 · Jul 14, 2026
arXivHTMLPDF

Abstract

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 H$_2$, LiH, BeH$_2$, H$_2$O, and NH$_3$. 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 H$_2$ and HeH$^+$ using both STO-3G and 6-31G basis sets, and for LiH and BeH$_2$ 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 $\times$ 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.

Comments
86 pages, 41 figures

Related papers

Ground-State Energy Estimation of HeH$^{+}$, ArH$^{+}$, and H$_2$O via Sample-Based Quantum Diagonalization

Jubin Park, Chae-Hyun Yoon, Minkyu Lee +1 more

Accurate ground-state energies are essential for understanding molecular structure, chemical bonding, and reaction energetics in quantum chemistry. In this work, we investigate the ground-state properties of the molecular systems HeH$^+$, ArH$^+$, and H$_2$O using Sample-Based Quantum Diagonalization (SQD), a hybrid quantum-classical framework designed for near-term quantum devices. Unlike variational approaches such as VQE, which require deep parameterized circuits and repeated expectation-value measurements, SQD reconstructs a low-energy determinant subspace directly from measured bitstrings. For the present calculations, bitstrings were generated on IBM quantum hardware using shallow local unitary cluster Jastrow (LUCJ) circuits whose parameters were constructed from the $t_1$ and $t_2$ amplitudes of coupled-cluster singles and doubles (CCSD) calculations based on restricted Hartree--Fock (RHF) references. From these samples, we compute ground-state potential-energy curves of HeH$^+$, ArH$^+$, and H$_2$O with the 6-31G and cc-pVDZ basis sets. For all three systems, the SQD results obtained with the cc-pVDZ basis closely follow the corresponding same-basis CCSD energies and reproduce the equilibrium-region trends of the potential-energy curves. HeH$^+$ and ArH$^+$ were chosen as simple yet astrophysically important molecular-ion benchmarks, while H$_2$O was included as a representative polyatomic molecule to assess the applicability of SQD beyond diatomic ionic systems. At the adopted equilibrium geometries, the deviations from the same-active-space CASCI references are 0.00, 2.51, and 6.34 mHa for HeH$^+$, ArH$^+$, and H$_2$O, respectively. These results demonstrate the feasibility of hardware-assisted SQD for the present benchmark systems and motivate further studies of its accuracy and computational scaling for larger molecular active spaces.

arXiv:2608.06415•Aug 5, 2026•0 citations

Quantum Simulation of Ligand-like Molecules through Sample-based Quantum Diagonalization in Density Matrix Embedding Framework

Ashish Kumar Patra, Anurag K. S. V., Sai Shankar P. +4 more

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.

quant-pharXiv:2511.22158•The Journal of Chemical Physics•Nov 28, 2025•0 citations

Towards quantum-centric simulations of extended molecules: sample-based quantum diagonalization enhanced with density matrix embedding theory

Akhil Shajan, Danil S. Kaliakin, Abhishek Mitra +9 more

Computing ground-state properties of molecules is a promising application for quantum computers operating in concert with classical high-performance computing resources. Quantum embedding methods are a family of algorithms particularly suited to these computational platforms: they combine high-level calculations on active regions of a molecule with low-level calculations on the surrounding environment, thereby avoiding expensive high-level full-molecule calculations and allowing to distribute computational cost across multiple and heterogeneous computing units. Here, we present the first density matrix embedding theory (DMET) simulations performed in combination with the sample-based quantum diagonalization (SQD) method. We employ the DMET-SQD formalism to compute the ground-state energy of a ring of 18 hydrogen atoms, and the relative energies of the chair, half-chair, twist-boat, and boat conformers of cyclohexane. The full-molecule 41- and 89-qubit simulations are decomposed into 27- and 32-qubit active-region simulations, that we carry out on the ibm_cleveland device, obtaining results in agreement with reference classical methods. Our DMET-SQD calculations mark a tangible progress in the size of active regions that can be accurately tackled by near-term quantum computers, and are an early demonstration of the potential for quantum-centric simulations to accurately treat the electronic structure of large molecules, with the ultimate goal of tackling systems such as peptides and proteins.

arXiv:2411.09861•Nov 15, 2024•4 citations

Quantum Computing Enabled ab initio Molecular Dynamics Simulations

Susanta Das, Subhamoy Bhowmik, Zhen Li +4 more

We demonstrate a quantum-classical workflow for ab initio molecular dynamics (AIMD) in which quantum measurements from a chemistry-inspired LUCJ ansatz are post-processed using Sample-based Quantum Diagonalization (SQD) to recover determinant subspaces and deliver energies and analytical nuclear gradients for dynamics. As an exact benchmark, we use full configuration interaction (FCI) in the STO-3G basis, enabling a direct assessment of the accuracy of SQD. In gas-phase benchmarks, SQD reproduces FCI energies and gradients to within 1 kcal mol$^{-1}$ of the FCI reference and yields stable AIMD trajectories. In explicit-solvent QM/MM simulations, SQD retains this agreement, matching FCI energy fluctuations and RMS gradient profiles and reproducing solute-solvent structure as quantified by radial distribution functions. Overall, these benchmarks establish LUCJ+SQD as a practical route for integrating current quantum hardware into QM/MM molecular dynamics and provide an early demonstration of condensed-phase QM/MM dynamics driven by a quantum electronic-structure engine.

physics.chem-phCross-listedarXiv:2607.28548•Jul 31, 2026•0 citations

Toward Quantum-Centric Simulations of Extended Molecules: Sample-Based Quantum Diagonalization Enhanced with Density Matrix Embedding Theory

Akhil Shajan, Danil S. Kaliakin, Abhishek Mitra +9 more

Computing ground-state properties of molecules is a promising application for quantum computers operating in concert with classical high-performance computing resources. Quantum embedding methods are a family of algorithms particularly suited to these computational platforms: they combine high-level calculations on active regions of a molecule with low-level calculations on the surrounding environment, thereby avoiding expensive high-level full-molecule calculations and allowing the distribution of computational cost across multiple and heterogeneous computing units. Here, we present the first density matrix embedding theory (DMET) simulations performed in combination with a sample-based quantum diagonalization (SQD) method. We employ the DMET-SQD formalism to compute the ground-state energy of a ring of 18 hydrogen atoms and the relative energies of the chair, half-chair, twist-boat, and boat conformers of cyclohexane. The full-molecule 41- and 89-qubit simulations are decomposed into 27- and 32-qubit active-region simulations, which we carry out on the ibm_cleveland device, obtaining results in agreement with reference classical methods. Our DMET-SQD calculations mark tangible progress in the size of active regions that can be accurately tackled by near-term quantum computers and are an early demonstration of the potential for quantum-centric simulations to accurately treat the electronic structure of large molecules, with the ultimate goal of tackling systems such as peptides and proteins.

Journal of Chemical Theory and Computation•Jul 8, 2025•21 citations

From Promise to Practice: Benchmarking Quantum Chemistry on Quantum Hardware

Osama M. Raisuddin, Haimeng Zhang, Mario Motta +1 more

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.

quant-pharXiv:2512.01012•Dec 2, 2025•0 citations
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_22​, LiH, BeH2_22​, H2_22​O, and NH3_33​. 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_22​ and HeH+^++ using both STO-3G and 6-31G basis sets, and for LiH and BeH2_22​ 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 ×\times× 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.
Ground-State Energy Estimation of HeH+^{+}+, ArH+^{+}+, and H2_22​O via Sample-Based Quantum Diagonalization
Accurate ground-state energies are essential for understanding molecular structure, chemical bonding, and reaction energetics in quantum chemistry. In this work, we investigate the ground-state properties of the molecular systems HeH+^++, ArH+^++, and H2_22​O using Sample-Based Quantum Diagonalization (SQD), a hybrid quantum-classical framework designed for near-term quantum devices. Unlike variational approaches such as VQE, which require deep parameterized circuits and repeated expectation-value measurements, SQD reconstructs a low-energy determinant subspace directly from measured bitstrings. For the present calculations, bitstrings were generated on IBM quantum hardware using shallow local unitary cluster Jastrow (LUCJ) circuits whose parameters were constructed from the t1t_1t1​ and t2t_2t2​ amplitudes of coupled-cluster singles and doubles (CCSD) calculations based on restricted Hartree--Fock (RHF) references. From these samples, we compute ground-state potential-energy curves of HeH+^++, ArH+^++, and H2_22​O with the 6-31G and cc-pVDZ basis sets. For all three systems, the SQD results obtained with the cc-pVDZ basis closely follow the corresponding same-basis CCSD energies and reproduce the equilibrium-region trends of the potential-energy curves. HeH+^++ and ArH+^++ were chosen as simple yet astrophysically important molecular-ion benchmarks, while H2_22​O was included as a representative polyatomic molecule to assess the applicability of SQD beyond diatomic ionic systems. At the adopted equilibrium geometries, the deviations from the same-active-space CASCI references are 0.00, 2.51, and 6.34 mHa for HeH+^++, ArH+^++, and H2_22​O, respectively. These results demonstrate the feasibility of hardware-assisted SQD for the present benchmark systems and motivate further studies of its accuracy and computational scaling for larger molecular active spaces.
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.
Computing ground-state properties of molecules is a promising application for quantum computers operating in concert with classical high-performance computing resources. Quantum embedding methods are a family of algorithms particularly suited to these computational platforms: they combine high-level calculations on active regions of a molecule with low-level calculations on the surrounding environment, thereby avoiding expensive high-level full-molecule calculations and allowing to distribute computational cost across multiple and heterogeneous computing units. Here, we present the first density matrix embedding theory (DMET) simulations performed in combination with the sample-based quantum diagonalization (SQD) method. We employ the DMET-SQD formalism to compute the ground-state energy of a ring of 18 hydrogen atoms, and the relative energies of the chair, half-chair, twist-boat, and boat conformers of cyclohexane. The full-molecule 41- and 89-qubit simulations are decomposed into 27- and 32-qubit active-region simulations, that we carry out on the ibm_cleveland device, obtaining results in agreement with reference classical methods. Our DMET-SQD calculations mark a tangible progress in the size of active regions that can be accurately tackled by near-term quantum computers, and are an early demonstration of the potential for quantum-centric simulations to accurately treat the electronic structure of large molecules, with the ultimate goal of tackling systems such as peptides and proteins.
We demonstrate a quantum-classical workflow for ab initio molecular dynamics (AIMD) in which quantum measurements from a chemistry-inspired LUCJ ansatz are post-processed using Sample-based Quantum Diagonalization (SQD) to recover determinant subspaces and deliver energies and analytical nuclear gradients for dynamics. As an exact benchmark, we use full configuration interaction (FCI) in the STO-3G basis, enabling a direct assessment of the accuracy of SQD. In gas-phase benchmarks, SQD reproduces FCI energies and gradients to within 1 kcal mol−1^{-1}−1 of the FCI reference and yields stable AIMD trajectories. In explicit-solvent QM/MM simulations, SQD retains this agreement, matching FCI energy fluctuations and RMS gradient profiles and reproducing solute-solvent structure as quantified by radial distribution functions. Overall, these benchmarks establish LUCJ+SQD as a practical route for integrating current quantum hardware into QM/MM molecular dynamics and provide an early demonstration of condensed-phase QM/MM dynamics driven by a quantum electronic-structure engine.
Computing ground-state properties of molecules is a promising application for quantum computers operating in concert with classical high-performance computing resources. Quantum embedding methods are a family of algorithms particularly suited to these computational platforms: they combine high-level calculations on active regions of a molecule with low-level calculations on the surrounding environment, thereby avoiding expensive high-level full-molecule calculations and allowing the distribution of computational cost across multiple and heterogeneous computing units. Here, we present the first density matrix embedding theory (DMET) simulations performed in combination with a sample-based quantum diagonalization (SQD) method. We employ the DMET-SQD formalism to compute the ground-state energy of a ring of 18 hydrogen atoms and the relative energies of the chair, half-chair, twist-boat, and boat conformers of cyclohexane. The full-molecule 41- and 89-qubit simulations are decomposed into 27- and 32-qubit active-region simulations, which we carry out on the ibm_cleveland device, obtaining results in agreement with reference classical methods. Our DMET-SQD calculations mark tangible progress in the size of active regions that can be accurately tackled by near-term quantum computers and are an early demonstration of the potential for quantum-centric simulations to accurately treat the electronic structure of large molecules, with the ultimate goal of tackling systems such as peptides and proteins.
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.