Calculating the molecular ground-state energy is a central challenge in computational chemistry. Conventional methods such as the Complete Active Space Configuration Interaction scale exponentially with molecular size, limiting their applicability to large molecules. Quantum computing offers a promising alternative by mapping molecular Hamiltonians by qubits, enabling cheaper computational scaling. Previous studies have shown that it is possible to formulate molecular ground state calculations as discrete optimization problems, addressable by quantum annealing. However, these efforts have been limited by previous generations of hardware and suboptimal annealing techniques. Here, the
H2O ground-state problem is mapped to an Ising Hamiltonian using the Xian-Bias-Kas (XBK) method. By taking advantage of enhanced qubit connectivity and shorter embedding chains, it is solved with a more than doubled probability of achieving Hartree-Fock-level solutions with respect to the most advanced predecessor. Advanced annealing strategies extend Hartree-Fock-level accuracy to significantly larger problem instances, enabling solutions that use nearly 2.5 times more physically embedded qubits than the largest cases previously reported and allowing to improve annealing results by two orders of magnitude, reaching an energy difference of 0.120~Hartree relative to Hartree-Fock. These results show tangible progress toward practical quantum annealing applications in NISQ era.
Quantum annealers are an alternative approach to quantum computing, which make use of the adiabatic theorem to efficiently find the ground state of a physically realizable Hamiltonian. Such devices are currently commercially available and have been successfully applied to several combinatorial and discrete optimization problems. However, the application of quantum annealers to problems in chemistry remains a relatively sparse area of research due to the difficulty in mapping molecular systems to the Ising model Hamiltonian. In this paper, we review two different methods for finding the ground state of molecular Hamiltonians using Ising model-based quantum annealers. In addition, we compare the relative effectiveness of each method by calculating the binding energies, bond lengths, and bond angles of the H3+ and H2O molecules and mapping their potential energy curves. We also assess the resource requirements of each method by determining the number of qubits and computation time required to simulate each molecule using various parameter values. While each of these methods is capable of accurately predicting the ground state properties of small molecules, we find that they are still outperformed by modern classical algorithms and that the scaling of the resource requirements remains a challenge.
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 chemistry calculations for small molecules on quantum hardware have been demonstrated to date only on universal-gate quantum computers, not quantum annealers. The latter devices are limited to finding the lowest eigenstate of the Ising Hamiltonian whereas the electronic Hamiltonian could not be mapped to the Ising form without exponential growth of the Ising Hamiltonian with the size of the system [J. Phys. Chem. B 122, 3384 (2018)]. Here we propose a novel mixed discrete-continuous optimization algorithm, which finds the lowest eigenstate of the qubit coupled cluster (QCC) method using a quantum annealer for solving a discrete part of the problem. The QCC method is a potentially exact approach for constructing the electronic wave function in the qubit space. Therefore, our methodology allows for systematically improvable quantum chemistry calculations using quantum annealears. We illustrate capabilities of our approach by calculating QCC ground electronic states for the LiH, H
2O, and C
6H
6 molecules. C
6H
6 calculations involve 36 qubits and are the largest quantum chemistry calculations made on a quantum annealer (the D-Wave 2000Q system) to date. Our findings opens up a new perspective for use quantum annealers in high-throughput material discovery.
Quantum computation has the potential to revolutionize quantum chemistry through major speedups in computation times and an exponential reduction in computational resources. Here, we combine the symmetry-adapted Jordan-Wigner encoding based on the full Boolean symmetry group Z2k with our new implementation of the Xia-Bian-Kais (XBK) method for improving the efficiency of electronic structure theory calculations on quantum annealers, particularly by reducing the number of qubits needed to achieve the same accuracy. By providing a more extensive symmetry-adapted encoding (SAE) than previous work, we are able to simulate molecules larger than those previously reported that have been studied using methods developed for quantum annealers and without using an active space. We calculated the potential energy surfaces of H2, LiH, He2, H2O, O2, N2, Li2, F2, CO, BH3, NH3, and CH4, with the largest molecule in the STO-6G basis set requiring 16 qubits with our SAE, and compared them with full configuration interaction results. The application of SAE to the XBK method provides an exponential reduction in the size of the Hilbert space and scales well with the size of the problem. It does not introduce significant additional errors for even or large values of a key variational parameter that determines the number of ancilla qubits used in the XBK method's Hamiltonian embedding, or for certain molecules such as He2 and H2O. We provide an explanation for this behavior and a recommendation on the usage of our method. In addition, we briefly discuss the potential of extracting electronic excited states from our method.
Abstract To understand the properties and interactions of materials, and determining the ground state energies is one of the important challenges in quantum chemistry, materials science, and quantum mechanics, where quantum computing can play an important role for studying the properties of materials. In this study, we have explored the quantum processor application to compute the He molecule ground state energy which utilizes the Variational Quantum Eigensolver (VQE) algorithm. In here, we have implemented VQE on a state-of-the-art quantum processor, optimizing a parameterized quantum circuit to minimize the energy expectation value of the He molecule's Hamiltonian on the four qubits processor. The obtained results of this work show a significant improvement in accuracy compared to classical computational methods, such as Hartree-Fock and density functional theory, which demonstrate the compute potential of quantum algorithms in quantum many-body problems. Thus, these results demonstrate the advantages of quantum computing in achieving high accuracy in simulations of molecular and material properties, and pave the way for future applications in more complex systems. This work highlights the potential of quantum processors in the fields of quantum chemistry, computational physics, and data science.