This work describes resource‑efficient quantum algorithms for molecular physics that exploit geometric structure and correspondences between physical and qubit symmetries. It begins with a review of classical approaches to the molecular electronic‑structure problem and of the variational quantum eigensolver in quantum computing. We then discuss the symmetry‑adapted encodings we introduced to reduce qubit requirements by eliminating redundant degrees of freedom arising from point‑group and other Boolean symmetries in the quantum simulation of molecular systems. Building on these encodings, we describe a qubit‑based symmetry‑adapted complete active‑space method that permits the simulation of larger molecular systems in a resource‑efficient manner on near‑term quantum devices. We present the QuantumSymmetry software package and illustrate its application to representative molecular systems. Next, we address electron‑molecule scattering by introducing a quantum computing algorithm to solve the inner‑region problem within the R‑matrix framework. In particular, we propose a novel simultaneous diagonalisation algorithm with subspace optimisation that requires only a limited number of Hamiltonian evaluations and could be generalised to a wide range of problems. The final chapter derives the explicit form of qubit operators for particle‑number projection. Altogether, these developments advance practical quantum simulations of complex molecular phenomena by leveraging symmetry and geometry in quantum‑algorithm design.
This work describes resource‑efficient quantum algorithms for molecular physics that exploit geometric structure and correspondences between physical and qubit symmetries. It begins with a review of classical approaches to the molecular electronic‑structure problem and of the variational quantum eigensolver in quantum computing. We then discuss the symmetry‑adapted encodings we introduced to reduce qubit requirements by eliminating redundant degrees of freedom arising from point‑group and other Boolean symmetries in the quantum simulation of molecular systems. Building on these encodings, we describe a qubit‑based symmetry‑adapted complete active‑space method that permits the simulation of larger molecular systems in a resource‑efficient manner on near‑term quantum devices. We present the QuantumSymmetry software package and illustrate its application to representative molecular systems. Next, we address electron‑molecule scattering by introducing a quantum computing algorithm to solve the inner‑region problem within the R‑matrix framework. In particular, we propose a novel simultaneous diagonalisation algorithm with subspace optimisation that requires only a limited number of Hamiltonian evaluations and could be generalised to a wide range of problems. The final chapter derives the explicit form of qubit operators for particle‑number projection. Altogether, these developments advance practical quantum simulations of complex molecular phenomena by leveraging symmetry and geometry in quantum‑algorithm design.
Computer simulation is central to understanding the behaviour of matter. Yet, more than a century after the development of quantum theory, first-principles simulation of chemistry remains a daunting challenge, owing to the exponential growth in memory required to represent interacting quantum particles such as electrons in molecules. Quantum computing offers, in principle, a fundamentally more efficient paradigm for simulating quantum systems, with quantum chemistry often cited as a flagship application. However, many quantum algorithms face intrinsic limitations, and any genuine quantum advantage must surpass the strong baseline set by decades of progress in classical algorithms and hardware. While classical methods often suffice for simulating weakly correlated systems, systems with many strongly correlated electrons remain a challenge for both classical and quantum approaches. In this thesis, we develop quantum algorithms for solving the time-independent Schrödinger equation, with a focus on molecular systems with strongly correlated electrons. We introduce a framework for identifying a compact set of states that parametrise the many-body wavefunction by capturing the correlation due to spin coupling between localised electrons. This representation requires very few variational parameters but incurs a very high memory cost on classical computers. To overcome this, we design quantum circuits that efficiently prepare such states by exploiting their symmetry, encoding the dominant entanglement structure of electronic wavefunctions into the initial qubit states. Through theoretical analysis and classical simulations of quantum computers, we show that this substantially reduces the quantum resources required to accurately compute ground- and excitedstate properties. Finally, we develop an approach for approximating the optimal ansatz within the variational quantum eigensolver, an algorithm designed for near-term quantum devices. Our method constructs wavefunctions from products of symmetrypreserving unitary operators and couples discrete and continuous optimisation to yield highly compact representations of strongly correlated electrons using shallow quantum circuits. Together, these results establish complementary directions for embedding problemspecific structure into quantum simulation algorithms, advancing the practical simulation of strongly correlated electrons.
Abstract A symmetry-adapted fermion-to-spin mapping or encoding that is able to store information about the occupancy of the n spin-orbitals of a molecular system into a lower number of n − k qubits in a quantum computer (where the number of reduced qubits k ranges from 2 to 5 depending on the symmetry of the system) is introduced. This mapping reduces the computational cost of a quantum computing simulation and at the same time enforces symmetry constraints. These symmetry-adapted encodings (SAEs) can be explicitly seen as a block-diagonalization of the Jordan–Wigner qubit Hamiltonian, followed by an orthogonal projection. We provide the form of the Clifford tableau for a general class of fermion-to-qubit encodings, and then use it to construct the map that block-diagonalizes the Hamiltonian in the SAEs. The algorithm proposed does not require any further computations to obtain this map, which is derived directly from the character table of the molecular point group. An implementation of the algorithm is presented as an open-source Python package, QuantumSymmetry, a user guide and code examples. QuantumSymmetry uses open-source quantum chemistry software PySCF for Hartree–Fock calculations, and is compatible with quantum computing toolsets OpenFermion and Qiskit. QuantumSymmetry takes arbitrary user input such as the molecular geometry and atomic basis set to construct the qubit operators that correspond in the appropriate SAE to fermionic operators on the molecular system, such as the second-quantized electronic structure Hamiltonian. QuantumSymmetry is used to produce numerical examples of variational quantum algorithm simulations to find the ground state energy for a number of example molecules, for both Unitary Coupled Clusters with Singles and Doubles and Adaptive Derivative Assembled Pseudo-Trotter Variational Quantum Eigensolver ansätze. We show that, beyond the advantage given by the lower qubit count, the proposed encodings consistently result in shallower and less complex circuits with a reduced number of variational parameters that are able to reach convergence faster and without any loss of computed accuracy.
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.
We present a symmetry-adapted qubit encoding with complete active space (SAE-CAS) for quantum chemistry on fault-tolerant and near-term quantum processors. Building on exact-symmetry encodings, we extend symmetry-adapted mappings to approximate Z-symmetries corresponding to frozen-core and virtual orbitals, thereby reducing qubit requirements without significant loss of accuracy. We derive the mapping from the second-quantised Hamiltonian to active-space qubit Hamiltonians, prove its equivalence to the canonical CAS Hamiltonian with frozen-core and virtual-orbital projection, and integrate it with point-group and spin-parity symmetry encodings via affine Clifford transformations to maximise qubit reduction while preserving the target symmetry sector. The same framework also accommodates the Bravyi--Kitaev mapping, yielding an SAE-CAS-BK variant that is unitarily equivalent to SAE-CAS. Numerical benchmarking on nine small molecules using UCCSD and a hardware-efficient shifted-circular-alternating (HE-SCA) ansatz shows that SAE-CAS reduces qubit counts and Pauli-operator weight, yields shallower circuits with fewer parameters, and often accelerates VQE convergence; with HE-SCA it consistently reaches CAS reference energies in cases where JW-CAS does not converge within the tested budgets. We provide an open-source implementation in the Python package QuantumSymmetry. SAE-CAS offers a route to resource-efficient molecular simulations on fault-tolerant and near-term quantum processors.
Accurate prediction of molecular structures and properties is vital for chemistry; materials science, and drug discovery, yet classical electronic-structure methods often fail for strongly correlated systems, large basis sets, and complex potential-energy landscapes. Quantum technology encompassing quantum computing, quantum machine learning, and hybrid quantum classical strategies offers a fundamentally new paradigm by encoding many-electron wave functions directly on qubits and exploiting superposition and entanglement to explore exponentially large Hilbert spaces. This review synthesizes recent algorithmic and hardware advances relevant to molecular modelling, including the Variational Quantum Eigensolver (VQE), Quantum Phase Estimation (QPE), quantum unitary coupled-cluster (q-UCC) families, equation-of-motion and subspace methods for excited states, tensor-network hybrids, and quantum kernel and variational QML approaches. We examine noise-aware hybrid workflows, error-mitigation techniques, symmetry-preserving ansätze, and operator-factorization methods that reduce measurement and gate overhead. Representative applications are discussed: ground- and excited-state energy prediction, potential-energy surface mapping, geometry and transition-state optimization, spectroscopic property estimation (IR, UV–Vis, NMR, EPR), and reaction-dynamics scenarios where non-adiabatic effects and conical intersections dominate. Resource estimates and scaling analyses clarify current NISQ limitations qubit counts, circuit depth, shot complexity and delineate the roadmap to fault-tolerant QPE for chemical accuracy. We compare quantum approaches with classical baselines (DFT, CCSD (T), multireference methods), identifying domains where quantum methods already show promise (strong correlation, multi-reference dissociation, spin-state ordering) and where classical methods remain competitive. Finally, we highlight near-term industrial opportunities in drug design, catalysis, CO2 capture, and energy materials, and outline critical research directions: algorithmic reductions in measurement/precision cost, hardware improvements in fidelity and connectivity, scalable ansatz design, and integrated software stacks for reproducible hybrid simulations. Together, these developments indicate that while practical, large-scale quantum advantage for general chemistry remains future work, quantum technologies are rapidly maturing into powerful tools for targeted molecular problems that are intractable with existing classical techniques.
The exact evaluation of the molecular ground state in quantum chemistry requires an exponentially increasing computational cost. Quantum computation is a promising way to overcome the exponential problem using polynomial-time quantum algorithms. A quantum-classical hybrid optimization scheme known as the variational quantum eigensolver is preferred for noisy intermediate-scale quantum devices. However, the circuit depth becomes one of the bottlenecks of its application to large molecules of more than 20 qubits. In this work, we employ point-group symmetry to reduce the number of operators in constructing ansatz so as to achieve a more compact quantum circuit. We illustrate this methodology with a series of molecules ranging from LiH (12 qubits) to C2H4 (28 qubits). A significant reduction of up to 82% of the operator numbers is reached on C2H4. This also sheds light onto further work in this direction to construct even shallower ansatz with enough expressive power and simulate even larger scale systems.