In the field of ab initio polaritonic chemistry, the use of quantum computers may offer, in the long term, a promising computational pathway for simulating and exploring novel cavity-induced light-matter phenomena. Addressing such problems on emerging platforms, however, already presents significant challenges that need to be considered for near-term applications. Within the framework of the variational quantum eigensolver (VQE), a key question concerns the development of quantum circuit Ans\"atze capable of encoding strongly correlated electron-photon states at a reasonable resource cost. Motivated by this question, in this work we develop different Ansatz strategies covering not only conventional qubit-based approaches but also higher-dimensional information units. Specifically, we introduce a set of circuits tailored to three different quantum computing architectures: (i) qubits, (ii) qudits, and (iii) hybrid qubit-qumode units of information. All circuits are designed to be both compact and physically motivated, with a central element being the development of spin-adapted electron-photon entangling fragment circuits tailored to the native capabilities of each platform. Integrating the different Ans\"atze within the state-averaged VQE, we numerically benchmark the three strategies on a cavity-embedded ${\mathrm{H}}_{2}$ molecule using state-vector simulations. Our results show that each proposed Ansatz achieves high accuracy in encoding the ground, first, and second polaritonic eigenstates of the system, even reproducing characteristic light-matter phenomena such as light-induced avoided crossings and conical intersections. This prospective work, beyond introducing multiple computational strategies adapted to different platforms, also highlights the potential of higher-dimensional quantum architectures to address hybrid fermion-boson systems.
In the field of ab initio polaritonic chemistry, the use of quantum computers may offer, in the long term, a promising computational pathway for simulating and exploring novel cavity-induced light-matter phenomena. Addressing such problems on emerging platforms, however, already presents significant challenges that need to be considered for near-term applications. Within the framework of the variational quantum eigensolver (VQE), a key question concerns the development of quantum circuit Ans"atze capable of encoding strongly correlated electron-photon states at a reasonable resource cost. Motivated by this question, in this work we develop different Ansatz strategies covering not only conventional qubit-based approaches but also higher-dimensional information units. Specifically, we introduce a set of circuits tailored to three different quantum computing architectures: (i) qubits, (ii) qudits, and (iii) hybrid qubit-qumode units of information. All circuits are designed to be both compact and physically motivated, with a central element being the development of spin-adapted electron-photon entangling fragment circuits tailored to the native capabilities of each platform. Integrating the different Ans"atze within the state-averaged VQE, we numerically benchmark the three strategies on a cavity-embedded H2 molecule using state-vector simulations. Our results show that each proposed Ansatz achieves high accuracy in encoding the ground, first, and second polaritonic eigenstates of the system, even reproducing characteristic light-matter phenomena such as light-induced avoided crossings and conical intersections. This prospective work, beyond introducing multiple computational strategies adapted to different platforms, also highlights the potential of higher-dimensional quantum architectures to address hybrid fermion-boson systems.
One of the first practical applications of quantum computers is expected to be molecular modelling. Performing this task would profoundly affect areas such as chemistry, materials science and drug synthesis. Modelling of molecules, which are classically intractable, can be achieved with just over 30 qubits, whereas state of the art quantum computers already have more than 50 qubits. The Variational Quantum Eigensolver (VQE) algorithm and VQE based protocols, are promising candidates to enable this task on emerging Noisy Intermediate-Scale Quantum (NISQ) computers. These protocols require short quantum circuits and short coherence times, and are particularly resilient to quantum errors. Nevertheless, there is still a significant gap between the accuracy and the coherence times of current NISQ computers, and the hardware requirements of VQE protocols to simulate practically interesting molecules. In this thesis, I present my contribution to narrowing this gap by developing VQE protocols for molecular modelling that are less demanding on quantum hardware. The VQE relies on the Rayleigh-Ritz variational principle to estimate the eigenvalues of a Hamiltonian operator, by minimizing its expectation value with respect to a trial quantum state, prepared by an ansatz. A major challenge for the practical realisation of VQE protocols on NISQ computers is to construct an ansatz that: (1) can accurately approximate the eigenstates of the Hamiltonian; (2) is easy to optimize; and (3) can be implemented by a shallow circuit, within the capabilities of a NISQ computer. The most widely used, unitary coupled cluster (UCC), type of ans\"atze mathematically correspond to a product of unitary evolutions of fermionic excitation operators. Owing to their fermionic structure, UCC ans\"atze preserve the symmetries of electronic wavefunctions, and thus are accurate and easy to optimize. Nevertheless, UCC ans\"atze are implemented by high depth circuits, which severely limit the size of the molecules that can be reliably simulated on NISQ computers. In this thesis, I begin by constructing efficient quantum circuits to perform evolutions of fermionic excitation operators. The circuits are optimized in the number of two-qubit entangling gates, which are the current bottleneck of NISQ computers. Compared to the standard circuits used to implement evolutions of fermionic excitation operators, the circuits derived in this thesis reduce the number of two-qubit entangling gates by more than 70% on average. As an intermediate result, I also derive efficient circuits to perform evolutions of qubit excitation operators (excitation operators that account for qubit, rather than fermionic commutation relations). Even with the fermionic-excitation-evolution circuits derived here, UCC ans\"atze still require very long circuits, with a particularly large number of two-qubit entangling gates. In this thesis, I consider the use of alternative VQE ans\"atze, based on evolutions of qubit excitation operators. Due to not accounting for fermionic anticommutation, evolutions of qubit excitation operators can be performed by circuits that require asymptotically fewer two-qubit entangling gates. Furthermore, qubit excitation operators preserve many of the physical properties of fermionic excitation operators. Performing a number of classical numerical VQE simulations for small molecules, I show that qubit-excitation-based ans\"atze can approximate molecular electronic wavefunctions almost as accurately as fermionic-excitation-based ans\"atze. Hence, I argue that evolutions of qubit excitation operators are more suitable to construct molecular ans\"atze than evolutions of fermionic excitation operators, especially in the era of NISQ computers. Motivated by the advantage of qubit-excitation-based ans\"atze, I introduce the qubit-excitation-based adaptive variational quantum eigensolver (QEB-ADAPT-VQE). The QEB-ADAPT-VQE belongs to a family of ADAPT-VQE protocols for molecular modelling that grow a problem-tailored ansatz by iteratively appending unitary operators sampled from a predefined finite-size pool of operators. The operator at each iteration is sampled based on an ansatz-growing strategy, which aims to achieve the lowest estimate for the Hamiltonian expectation value at each iteration. In this way, ADAPT-VQE protocols construct shallow-circuit, few-parameter ans\"atze tailored specifically to the molecular systems of interest. In the case of the QEB-ADAPT-VQE, the operator pool is defined by a set of evolutions of single and double qubit excitation operators. I benchmark the performance of the QEB-ADAPT-VQE, by performing classical numerical simulations. I demonstrate that it can construct ans\"atze that are several orders of magnitude more accurate, and require significantly shallower circuits, than standard UCC ans\"atze. I also compare the QEB-ADAPT-VQE against the original fermionic-ADAPT-VQE, which utilizes a pool of fermionic excitation evolutions, and the qubit-ADAPT-VQE, which utilizes a pool of Pauli-string evolutions. I demonstrate that, in terms of circuit efficiency and convergence speed, the QEB-ADAPT-VQE systematically outperforms the qubit-ADAPT-VQE, which to my knowledge was the previous most circuit-efficient, scalable VQE protocol for molecular modeling. The QEB-ADAPT-VQE protocol, therefore represents a significant improvement in the field of VQE protocols for molecular modelling and brings us closer to achieving practical quantum advantage. Lastly, I outline a modified version of the QEB-ADAPT-VQE, the excited-QEB-ADAPT-VQE, designed to estimate energies of excited molecular states. The excited-QEB-ADAPT-VQE is more robust to initial simulation conditions, at the expense of increased computational complexity.
Quantum harmonic oscillators, or qumodes, provide a promising and versatile framework for quantum computing. Unlike qubits, which are limited to two discrete levels, qumodes have an infinite-dimensional Hilbert space, making them well-suited for a wide range of quantum simulations. In this work, we focus on the molecular electronic structure problem. We propose an approach to map the electronic Hamiltonian into a qumode bosonic problem that can be solved on bosonic quantum devices using the variational quantum eigensolver (VQE). Our approach is demonstrated through the computation of ground potential energy surfaces for benchmark model systems, including H 2 and the linear H 4 molecule. The preparation of trial qumode states and the computation of expectation values leverage universal ansatzes based on the echoed conditional displacement (ECD), or the selective number-dependent arbitrary phase (SNAP) operations. These techniques are compatible with circuit quantum electrodynamics (cQED) platforms, where microwave resonators coupled to superconducting transmon qubits can offer an efficient hardware realization. This work establishes a new pathway for simulating many-fermion systems, highlighting the potential of hybrid qubit-qumode quantum devices in advancing quantum computational chemistry.
Recent years have witnessed a surge of experimental and theoretical interest in controlling the properties of matter, such as its chemical reactivity, by confining it in optical cavities, where the enhancement of the light-matter coupling strength leads to the creation of hybrid light-matter states known as polaritons. However, ab initio calculations that account for the quantum nature of both the electromagnetic field and matter are challenging and have only started to be developed in recent years. We introduce a deep learning variational quantum Monte Carlo method to solve the electronic and photonic Schrödinger equations of molecules trapped in optical cavities. We extend typical electronic neural network wave function ansätze to describe joint fermionic and bosonic systems, i.e., electron-photon systems, in a quantum Monte Carlo framework. We apply our method to hydrogen molecules in a cavity, computing both ground and excited states. We assess their energy, dipole moment, charge density shift due to the cavity, the state of the photonic field, and the entanglement developed between the electrons and photons. Where possible, we compare our results with more conventional quantum chemistry methods proposed in the literature, finding good qualitative agreement, thus extending the range of scientific problems that can be tackled using machine learning techniques.
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.
One of the most significant challenges in quantum chemistry simulations is the accurate and efficient treatment of strongly correlated systems using classical computers. Quantum computational chemistry is an emerging field that has grown alongside the development of quantum computers, aiming to tackle challenges beyond the reach of traditional computational chemistry methods. Given the current limitations of quantum hardware, hybrid quantum-classical algorithms have become the dominant approach in quantum computational chemistry, with the Variational Quantum Eigensolver (VQE) serving as a representative example. In the realm of variational quantum algorithms, the design of the ansatz plays a pivotal role, as it directly determines both the accuracy of the results and the optimization efficiency of the algorithm, making it a central component in achieving success in quantum computational chemistry. However, existing ansätze often fails to strike a balance between precision and hardware efficiency, posing a major bottleneck for practical applications. In our previous work, we proposed a novel strategy for designing hardware-efficient ansätze with built-in physical constraints (PCHEA). This approach demonstrated outstanding performance in simulating one-dimensional Heisenberg models and several typical molecular systems. By embedding physical insights into the design, PCHEA not only improved computational accuracy but also significantly enhanced hardware compatibility, making it a promising candidate for near-term quantum devices. Nevertheless, challenges remain. When applying the one-dimensional PCHEA to two-dimensional systems, the circuit costs increase significantly. To address this challenge, in this work, we simplify the implementation of the physics-constrained hardware-efficient ansatz, using fSim gates arranged in hardware-efficient linear and brickwall layouts to construct one-dimensional fSim linear and fSim brickwall ansätze. This modification substantially reduces both the number of two-qubit gates and the circuit depth. These ansätze ensure universality, systematic improvement, and size consistency within a particle-number-conserving Hilbert space, with significant improvements observed in the simulation of the one-dimensional Heisenberg model. Furthermore, building on the advancements in quantum hardware connectivity, we extend the one-dimensional PCHEA to a two-dimensional version by arranging it in a grid pattern. The two-dimensional PCHEA not only retains the physical constraints of the one-dimensional PCHEA but also exploits hardware connectivity to further reduce circuit costs. Simulations on two-dimensional Heisenberg models show that the two-dimensional PCHEA achieves chemical accuracy with considerably lower circuit depth, fewer parameters, and fewer two-qubit gates compared to the one-dimensional PCHEA. This work reinforces the effectiveness of designing hardware-efficient ansätze with embedded physical constraints, presenting a promising approach to advancing the practical application of quantum computers for simulating two-dimensional systems. Such progress brings us closer to achieving more efficient and scalable quantum simulations while addressing key challenges in quantum computational chemistry.
The recent advent of quantum algorithms for noisy quantum devices offers a new route toward simulating strong light–matter interactions of molecules in optical cavities for polaritonic chemistry. In this work, we introduce a general framework for simulating electron–photon-coupled systems on small, noisy quantum devices. This method is based on the variational quantum eigensolver (VQE) with the polaritonic unitary coupled cluster (PUCC) ansatz. To achieve chemical accuracy, we exploit various symmetries in qubit reduction methods, such as electron–photon parity, and use recently developed error mitigation schemes, such as the reference zero-noise extrapolation method. We explore the robustness of the VQE-PUCC approach across a diverse set of regimes for the bond length, cavity frequency, and coupling strength of the H 2 molecule in an optical cavity. To quantify the performance, we measure two properties: ground-state energy, fundamentally relevant to chemical reactivity, and photon number, an experimentally accessible general indicator of electron–photon correlation.