Quantum computers have the potential to revolutionise our understanding of the microscopic behaviour of materials and chemical processes by enabling high-accuracy electronic structure calculations to scale more efficiently than is possible using classical computers. Current quantum computing hardware devices suffer from the dual challenges of noise and cost, which raises the question of what practical value these devices might offer before full fault tolerance is achieved and economies of scale enable cheaper access. Here we examine the practical value of noisy quantum computers as tools for high-accuracy electronic structure, by using a Quantinuum ion-trap quantum computer to predict the ionisation potential of helium. By combining a series of techniques suited for use with current hardware including qubit-efficient encoding coupled with chemical insight, low-cost variational optimisation with hardware-adapted quantum circuits, and moments-based corrections, we obtain an ionisation potential of 24.5536 (+0.0011, -0.0005) eV, which agrees with the experimentally measured value to within true chemical accuracy, and with high statistical confidence. The methods employed here can be generalised to predict other properties and expand our understanding of the value that might be provided by near-term quantum computers.
Abstract Quantum computers have the potential to perform accurate and efficient electronic structure calculations, enabling the simulation of properties of materials. However, today’s noisy, intermediate-scale quantum (NISQ) devices have a limited number of qubits and gate operations due to the presence of errors. Here, we propose a systematically improvable end-to-end pipeline to alleviate these limitations. Our proposed resource-efficient pipeline combines problem decomposition techniques for compact molecular representations, circuit optimization methods for compilation, solving the eigenvalue problem on advanced quantum hardware, and error-mitigation techniques in post-processing the results. Using the density matrix embedding theory for compact representation, and an ion-trap quantum computer, we simulate a ring of 10 hydrogen atoms taking into account all electrons equally and explicitly in the electronic structure calculation. In our experiment, we simulated the largest molecular system on a quantum computer within chemical accuracy with respect to total molecular energy calculated by the full CI method. Our methods reduce the number of qubits required for high-accuracy quantum simulations by an order of magnitude in the present work, enabling the simulation of larger, more industrially relevant molecules using NISQ devices. They are further systematically improvable as devices’ computational capacity continues to grow.
One of the most important application areas of molecular quantum chemistry is the study and prediction of chemical reactivity. Large-scale, fully error-tolerant quantum computers could provide exact or near-exact solutions to the underlying electronic structure problem with exponentially less effort than a classical computer thus enabling highly accurate predictions for comparably large molecular systems. In the nearer future, however, only “noisy” devices with a limited number of qubits that are subject to decoherence will be available. For such near-term quantum computers the hybrid quantum-classical variational quantum eigensolver algorithm in combination with the unitary coupled-cluster ansatz (UCCSD-VQE) [Peruzzo et al. Nat. Commun. 2014, 5, 4213 and McClean et al. New J. Phys. 2016, 18, 023023] has become an intensively discussed approach that could provide accurate results before the dawn of error-tolerant quantum computing. In this work we present an implementation of UCCSD-VQE that allows for the first time to treat both open- and closed-shell molecules. We study the accuracy of the obtained energies for nine small molecular systems as well as for four exemplary chemical reactions by comparing to well-established electronic structure methods like (nonunitary) coupled-cluster and density functional theory. Finally, we roughly estimate the required quantum hardware resources to obtain “useful” results for practical purposes.
Emergent quantum computing technologies are widely expected to provide novel approaches in the simulation of quantum chemistry. Despite rapid improvements in the scale and fidelity of quantum computers, high resource requirements make the execution of quantum chemistry experiments challenging. Typical experiments are limited in the number of qubits used, incur a substantial shot cost, or require complex architecture-specific optimization and error mitigation techniques. In this paper, we propose a conceptually simple benchmarking approach involving the use of multi-ancilla quantum phase estimation. Our approach is restricted to very small chemical systems, and does not scale favorably beyond molecular systems that can be described with
2 qubits; however, this restriction allows us to generate circuits that scale quadratically in gate count with the number of qubits in the readout register. This enables the execution of quantum chemistry circuits that act on many qubits, while producing meaningful results with limited shot counts. We use this technique (with
200 shots per experiment) to calculate the ground state energy of molecular hydrogen to
50 bits of precision (
8.9×10−16 hartree) on a
56-qubit trapped-ion quantum computer, negating Trotter error. Including Trotter error, we obtain between
32 and
36 bits of precision (
1.5∗10−10 and
6.0∗10−11 hartree respectively), vastly exceeding chemical accuracy (
1.6∗10−3 hartree) against Full Configuration Interaction. We consider application of the approach to deeper circuits, and discuss potential as a benchmark task for near-term quantum devices.
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.
Quantum computing is emerging as a new computational paradigm with the potential to transform several research fields including quantum chemistry. However, current hardware limitations (including limited coherence times, gate infidelities, and connectivity) hamper the implementation of most quantum algorithms and call for more noise-resilient solutions. We propose an explicitly correlated Ansatz based on the transcorrelated (TC) approach to target these major roadblocks directly. This method transfers, without any approximation, correlations from the wave function directly into the Hamiltonian, thus reducing the resources needed to achieve accurate results with noisy quantum devices. We show that the TC approach allows for shallower circuits and improves the convergence toward the complete basis set limit, providing energies within chemical accuracy to experiment with smaller basis sets and, thus, fewer qubits. We demonstrate our method by computing bond lengths, dissociation energies, and vibrational frequencies close to experimental results for the hydrogen dimer and lithium hydride using two and four qubits, respectively. To demonstrate our approach's current and near-term potential, we perform hardware experiments, where our results confirm that the TC method paves the way toward accurate quantum chemistry calculations already on today's quantum hardware.