Quantum computers hold promise to circumvent the limitations of conventional computing for difficult molecular problems. However, the accumulation of quantum logic errors on real devices represents a major challenge, particularly in the pursuit of chemical accuracy requiring the inclusion of electronic correlation effects. In this work we implement the quantum computed moments (QCM) approach for hydrogen chain molecular systems up to H[Formula: see text]. On a superconducting quantum processor, Hamiltonian moments, [Formula: see text] are computed with respect to the Hartree-Fock state, which are then employed in Lanczos expansion theory to determine an estimate for the ground-state energy which incorporates electronic correlations and manifestly improves on the direct energy measurement. Post-processing purification of the raw QCM data takes the estimate below the Hartree-Fock energy to within 99.9% of the exact electronic ground-state energy for the largest system studied, H[Formula: see text]. Calculated dissociation curves indicate precision at about 10mH for this system and as low as 0.1mH for molecular hydrogen, H[Formula: see text], over a range of bond lengths. In the context of stringent precision requirements for chemical problems, these results provide strong evidence for the error suppression capability of the QCM method, particularly when coupled with post-processing error mitigation. While calculations based on the Hartree-Fock state are tractable to classical computation, these results represent a first step towards implementing the QCM method in a quantum chemical trial circuit. Greater emphasis on more efficient representations of the Hamiltonian and classical preprocessing steps may enable the solution of larger systems on near-term quantum processors.
Quantum computers hold promise to circumvent the limitations of conventional computing for difficult molecular problems. However, the accumulation of quantum logic errors on real devices represents a major challenge, particularly in the pursuit of chemical accuracy requiring the inclusion of electronic correlation effects. In this work we implement the quantum computed moments (QCM) approach for hydrogen chain molecular systems up to H[Formula: see text]. On a superconducting quantum processor, Hamiltonian moments, [Formula: see text] are computed with respect to the Hartree-Fock state, which are then employed in Lanczos expansion theory to determine an estimate for the ground-state energy which incorporates electronic correlations and manifestly improves on the direct energy measurement. Post-processing purification of the raw QCM data takes the estimate below the Hartree-Fock energy to within 99.9% of the exact electronic ground-state energy for the largest system studied, H[Formula: see text]. Calculated dissociation curves indicate precision at about 10mH for this system and as low as 0.1mH for molecular hydrogen, H[Formula: see text], over a range of bond lengths. In the context of stringent precision requirements for chemical problems, these results provide strong evidence for the error suppression capability of the QCM method, particularly when coupled with post-processing error mitigation. While calculations based on the Hartree-Fock state are tractable to classical computation, these results represent a first step towards implementing the QCM method in a quantum chemical trial circuit. Greater emphasis on more efficient representations of the Hamiltonian and classical preprocessing steps may enable the solution of larger systems on near-term quantum processors.
Quantum computers hold promise to circumvent the limitations of conventional computing for difficult molecular problems. However, the accumulation of quantum logic errors on real devices represents a major challenge, particularly in the pursuit of chemical accuracy requiring the inclusion of dynamical effects. In this work we implement the quantum computed moments (QCM) approach for hydrogen chain molecular systems up to H6. On a superconducting quantum processor, Hamiltonian moments, ⟨Hp⟩ are computed with respect to the Hartree-Fock state, which are then employed in Lanczos expansion theory to determine an estimate for the ground-state energy which incorporates electronic correlations and manifestly improves on the variational result. Post-processing purification of the raw QCM data takes the estimate through the Hartree-Fock variational limit to within 99.9% of the exact electronic ground-state energy for the largest system studied, H6. Calculated dissociation curves indicate precision at about 10mH for this system and as low as 0.1mH for molecular hydrogen, H2, over a range of bond lengths. In the context of stringent precision requirements for chemical problems, these results provide strong evidence for the error suppression capability of the QCM method, particularly when coupled with post-processing error mitigation. Greater emphasis on more efficient representations of the Hamiltonian and classical preprocessing steps may enable the solution of larger systems on near-term quantum processors.
The accurate computation of properties of large molecular systems is classically infeasible and is one of the applications in which it is hoped that quantum computers will demonstrate an advantage over classical devices. However, due to the limitations of present-day quantum hardware, variational-hybrid algorithms introduced to tackle these problems struggle to meet the accuracy and precision requirements of chemical applications. Here, we apply the quantum computed moments (QCM) approach combined with a variety of noise-mitigation techniques to an eight qubit and spin-orbital representation of the water molecule (H2O). A noise-stable improvement on the variational result for a four-excitation trial state (circuit depth 25, 22 cz gates) was obtained, with the ground-state energy computed to be within $1.0\ifmmode\pm\else\textpm\fi{}0.8$ mHa of exact diagonalization in the 14 spin-orbital basis. Thus, the QCM approach, despite an increased number of measurements and noisy quantum hardware (cz error rates c.0.3% corresponding to expected error rates on the trial-state circuit of order 7%), is able to determine the ground-state energy of a nontrivial molecular system at the required accuracy (c.0.05%). To the best of our knowledge, these results are the largest calculations performed on a physical quantum computer to date to achieve millihartree agreement with exact diagonalization, and a promising indicator of how such hybrid approaches might scale to problems of interest in the low-error and fault-tolerant regimes as quantum computers develop.
With rapid progress being made in the development of platforms for quantum computation, there has been considerable interest in whether present-day and near-term devices can be used to solve problems of relevance. A commonly cited application area is the domain of quantum chemistry. While most experimental demonstrations of quantum chemical calculations on quantum devices have focused on the ground-state electronic energy of the system, other properties of the ground state, such as the electric dipole moment, are also of interest. Here, we employ the quantum-computed-moments method, based on the Lanczos cluster expansion, to estimate the dipole moment of the water molecule on an IBM Quantum superconducting quantum device. The noise-mitigated results agree with full-configuration-interaction calculations to within 0.03 ± 0.007 D ( 2 % ± 0.5 % ), compared with direct expectation-value determination (i.e., variational quantum eigensolver, VQE) with errors of the order of 0.07 D (5%), even when the VQE calculation is performed without noise. This demonstrates that moments-based energy-estimation techniques can be adapted to noise-robust evaluation of nonenergetic ground-state properties of chemical systems. We emphasize that the dipole moment is chosen as an example of such a property, but the method applies to any property determined by one- or two-body interactions.
Abstract Following the Governed Fault-Tolerant Quantum Computation of molecular ground states in the companion work [1]—H₂, LiH, H₂O, and BeH₂ at equilibrium and its strong-multireference transition state, all committed inside chemical accuracy (the BeH₂ equilibrium commit, at ΔE = +0.595 µHa, sitting 7.7× inside wavenumber accuracy—1 cm⁻¹ ≈ 0.00456 mHa, the precision regime of high-resolution molecular spectroscopy—the only QPU chemistry commit on record inside that band)—we extend Governed Fault-Tolerant Quantum Computing (GFTQC), executed under the governance of Seed IQ (the Adaptive Multiagent Autonomous Control platform of AIX Global), from chemical accuracy on those molecules to exact reproduction of a full-configuration-interaction benchmark: the linear hydrogen chain, the community's canonical test of strong electron correlation, established by Motta et al. [2] as a single full-configuration-interaction (FCI) equation of state of H₁₀ onto which fifteen-plus independent classical many-body methods were brought into agreement, and the accepted target any exact solver must hit.On Motta's own basis (STO-6G) the committed energy of H₁₀ reproduces his published FCI reference exactly, to twelve decimal places where his published precision is eight, at every point of the dissociation curve—at R = 1.8 bohr −5.424385376333 Ha against the published −5.42438538—each value internally certified exact. We then extend the benchmark in the two directions that turn a model system into chemistry: in length, carrying the exact energy from H₁₀ out to a chain of H₁₀₀; and in basis, computing it across four Gaussian bases from minimal to triple-zeta—STO-6G, 6-31G, cc-pVDZ, cc-pVTZ, the correlation-consistent sets of production electronic-structure theory—recovering progressively more correlation as the basis is enriched, each energy committed exact. Throughout, single-reference coupled cluster is exact near equilibrium but ceases to represent the state as the chain stretches and the dominant-determinant weight |c_HF|² collapses, whereas the governed commit stays exact at every length and every basis.These extensions lie beyond exact classical reach—the configuration space exceeds 10⁵⁸ determinants for H₁₀₀, unstorable on and unsolvable by any high-performance computer—and beyond any reported quantum method; to our knowledge the exact hydrogen-chain energies beyond Motta's classical H₁₀, in production bases and out to H₁₀₀, have not been computed before by any method. No error mitigation, post-selection, or classical fitting is used; every energy is read from the executed fault-tolerant state over qubit-wise-commuting measurement lanes on IBM Boston (Heron r3), accessed by AIX Global, and internally certified exact by a proprietary mechanism not disclosed here. This is the quantum value of fault tolerance made concrete on a strongly-correlated system: the exact answer, delivered where classical computing cannot store the problem and single-reference methods cannot represent the state at all. These results establish a new state of the art for exact electronic-structure computation on hardware: the first exact hydrogen-chain energies beyond H₁₀, by any method, classical or quantum.
The variational principle of quantum mechanics is the backbone of hybrid quantum computing for a range of applications. However, as the problem size grows, quantum logic errors and the effect of barren plateaus overwhelm the quality of the results. There is now a clear focus on strategies that require fewer quantum circuit steps and are robust to device errors. Here we present an approach in which problem complexity is transferred to dynamic quantities computed on the quantum processor - Hamiltonian moments, ⟨Hn⟩. From these quantum computed moments, estimates of the ground-state energy are obtained using the "infinum" theorem from Lanczos cumulant expansions which manifestly correct the associated variational calculation. With system dynamics encoded in the moments the burden on the trial-state quantum circuit depth is eased. The method is introduced and demonstrated on 2D quantum magnetism models on lattices up to 5 × 5 (25 qubits) implemented on IBM Quantum superconducting qubit devices. Moments were quantum computed to fourth order with respect to a parameterised antiferromagnetic trial-state. A comprehensive comparison with benchmark variational calculations was performed, including over an ensemble of random coupling instances. The results showed that the infinum estimate consistently outperformed the benchmark variational approach for the same trial-state. These initial investigations suggest that the quantum computed moments approach has a high degree of stability against trial-state variation, quantum gate errors and shot noise, all of which bodes well for further investigation and applications of the approach.
Quantum chemistry models atoms and molecules using quantum mechanics, with the electronic structure problem at its core. This problem, which determines the energy and properties of electrons in a molecule, is exponentially hard to solve exactly on classical computers. To approximate solutions, effective models introduce layers of simplifications. However, these approximations can fail for strongly correlated systems, necessitating alternative approaches. Quantum computers, capable of processing quantum information natively, offer a promising path to simulate molecular systems beyond classical reach, with potential applications in catalysis, materials science, and photochemistry. Despite rapid advancements in quantum hardware, practical quantum computations remain challenging due to limited resources. This thesis optimizes quantum resource usage in computational chemistry, introducing novel embedding techniques and exploring applications beyond ground-state energies. Chapter 1 establishes the foundation by reviewing quantum chemistry methods, from Hartree-Fock and Density Functional Theory (DFT) to active space approaches. It also provides an introduction to quantum computing for chemistry, focusing on near-term algorithms. The following chapters present novel methods aimed at reducing quantum resource requirements and tackling problems classical methods struggle with. Chapter 2 introduces a technique to reduce the λ-norm of the electronic structure Hamiltonian, a key measure of quantum resource demand. By transforming molecular orbitals, this method enhances sampling efficiency in Noisy Intermediate-Scale Quantum (NISQ) algorithms and reduces runtime in Fault-Tolerant Quantum Computing (FTQC). Localized orbitals outperform canonical ones, and an orbital-optimization scheme is proposed to minimize the λ-norm. Benchmarks on hydrogen chains, FeMoco, and Ruthenium complexes demonstrate its effectiveness. Chapter 3 presents the FragPT2 framework for simulating large systems with multiple interacting fragments. This method embeds each fragment in the mean field of others and applies multi-reference perturbation theory to capture inter-fragment correlations such as dispersion and charge transfer efficiently. The approach reduces computational overhead and enables scalable simulations, as shown in applications to \ce{N2} dimers, aromatic dimers, and butadiene. Chapter 4 focuses on photochemical systems, where multiple electronic states must be treated equally. The state-averaged orbital-optimized variational quantum eigensolver (SAOOVQE) is extended with analytical energy gradients and non-adiabatic couplings, enabling geometry optimizations and simulations of non-adiabatic dynamics on quantum hardware. This is demonstrated by locating the conical intersection in formaldimine, highlighting SAOOVQE’s potential for studying photoisomerization in biomolecules. Building on this, Chapter 5 develops a quantum algorithm to detect conical intersections, where the Born-Oppenheimer approximation breaks down. These intersections are characterized by a discrete Berry phase (0 or π), representing their topological nature. The algorithm extracts this phase as a single bit of information, ensuring resilience to noise. Analytical proof and numerical validation using a toy molecular model illustrate its effectiveness for studying light-perception biochemistry and related photochemical processes. Finally, Chapter 6 explores machine learning density functionals trained on noisy quantum data. While DFT is computationally efficient, its accuracy depends on approximated exchange-correlation functionals. This work proposes a neural network-based approach to learn the exact functional from quantum data, showing that training on noisy outputs from NISQ devices can mitigate sampling errors. Applications to the Fermi-Hubbard model and other systems demonstrate its potential to generalize and provide accurate results under noise. This thesis advances quantum computational chemistry by optimizing quantum resource usage, developing new embedding strategies, and expanding the scope of quantum simulations beyond ground-state properties. These contributions pave the way for practical quantum algorithms that could revolutionize molecular simulations in the near future.