Calculations of molecular spectral properties, like photodissociation rates and absorption bands, rely on knowledge of the excited state energies of the molecule of interest. Protocols based on the variational quantum eigensolver (VQE) are promising candidates to calculate such energies on emerging noisy intermediate scale quantum (NISQ) computers. The successful implementation of these protocols on NISQ computers, relies on ansätze that can accurately approximate the molecular states and that can be implemented by shallow quantum circuits. In this paper, we introduce the excited qubit-excitation-based adaptive (e-QEB-ADAPT)-VQE protocol to calculate molecular excited state energies. The e-QEB-ADAPT-VQE constructs efficient problem-tailored ansätze by iteratively appending evolutions of qubit excitation operators. The e-QEB-ADAPT-VQE is an adaptation of the QEB-ADAPT-VQE protocol, which is designed to be independent on the choice of an initial reference state. We perform classical numerical simulations for LiH and BeH
2 to benchmark the performance of the e-QEB-ADAPT-VQE. We demonstrate that the e-QEB-ADAPT-VQE can construct highly accurate ansätze that require at least an order of magnitude fewer
CNOTs than standard fixed UCC ansätze, such as the UCCSD and the GUCCSD.
Calculations of molecular spectral properties, like photodissociation rates and absorption bands, rely on knowledge of the excited state energies of the molecule of interest. Protocols based on the variational quantum eigensolver (VQE) are promising candidates to calculate such energies on emerging noisy intermediate-scale quantum (NISQ) computers. The successful implementation of these protocols on NISQ computers, relies on ans\"atze that can accurately approximate the molecular states and that can be implemented by shallow quantum circuits. We introduce the excited qubit-excitation-based adaptive (e-QEB-ADAPT)-VQE protocol to calculate molecular-excited-state energies. The e-QEB-ADAPT-VQE constructs efficient problem-tailored ans\"atze by iteratively appending evolutions of qubit excitation operators. The e-QEB-ADAPT-VQE also improves on previous ADAPT-VQE protocol in that it is designed to be independent on the choice of initial reference state. We perform classical numerical simulations for LiH and
BeH2 to benchmark the performance of the e-QEB-ADAPT-VQE. We demonstrate that the e-QEB-ADAPT-VQE can construct highly accurate ans\"atze that require at least an order of magnitude fewer cnot gates than standard fixed unitary coupled-cluster ans\"atze, such as the UCCSD and the GUCCSD. We also show that the e-QEB-ADAPT-VQE is more successful in constructing ans\"atze for excited molecular states than other ADAPT-VQE protocols.
Abstract Molecular simulations with the variational quantum eigensolver (VQE) are a promising application for emerging noisy intermediate-scale quantum computers. Constructing accurate molecular ansatze that are easy to optimize and implemented by shallow quantum circuits is crucial for the successful implementation of such simulations. Ansatze are, generally, constructed as series of fermionic-excitation evolutions. Instead, we demonstrate the usefulness of constructing ansatze with ``qubit-excitation evolutions', which, contrary to fermionic excitation evolutions, obey ``qubit commutation relations'. We show that qubit excitation evolutions, despite the lack of some of the physical features of fermionic excitation evolutions, accurately construct ansatze, while requiring asymptotically fewer gates. Utilizing qubit excitation evolutions, we introduce the iterative qubit excitation based VQE (IQEB-VQE) algorithm. The IQEB-VQE performs molecular simulations using a problem-tailored ansatz, grown iteratively by appending evolutions of single and double qubit excitation operators. By performing numerical simulations for small molecules, we benchmark the IQEB-VQE, and compare it against other competitive VQE algorithms. In terms of circuit efficiency and time complexity, we find that the IQEB-VQE systematically outperforms the previously most circuit-efficient, practically scalable VQE algorithms.
Classical simulation of molecular electronic structure is often limited by the rapid increase in computational cost associated with describing electron correlation. The Variational Quantum Eigensolver (VQE) is one hybrid quantum-classical approach ideally adapted to the Noisy Intermediate-Scale Quantum (NISQ) regime. This work applies the VQE algorithm to compute ground state energies of a set of molecules:
H2,H3+,LiH,HF,OH−,H2O,BH3,CH4, and
NH3. The approach utilizes the Unitary Coupled Cluster for Single and Double excitations (UCCSD) ansatz along with the Bravyi-Kitaev transformation for fermionic-to-qubit encoding. Quantum simulation outputs are systematically benchmarked against Full Configuration Interaction (FCI), a numerically exact traditional approach, and Unrestricted Hartree-Fock (UHF). The results establish a robust consistency of the VQE-obtained molecular ground state energies against standard FCI and UHF benchmarks. This research provides a critical benchmark, confirming that VQE has the capability for high-precision molecular simulations and playing the role of a pivotal stepping stone toward revolutionizing applications in drug discovery and materials science.
Simulation of electronic structure is one of the most promising applications on noisy intermediate-scale quantum (NISQ) era devices. However, NISQ devices suffer from a number of challenges like limited qubit connectivity, short coherence times, and sizable gate error rates. Thus, desired quantum algorithms should require shallow circuit depths and low qubit counts to take advantage of these devices. Here, we attempt to reduce quantum resource requirements for molecular simulations on a quantum computer while maintaining the desired accuracy with the help of classical quantum chemical theories of canonical transformation and explicit correlation. In this work, compact ab initio Hamiltonians are generated classically, in the second quantized form, through an approximate similarity transformation of the Hamiltonian with (a) an explicitly correlated two-body unitary operator with generalized pair excitations that remove the Coulombic electron-electron singularities from the Hamiltonian and (b) a unitary one-body operator to efficiently capture the orbital relaxation effects required for accurate description of the excited states. The resulting transcorrelated Hamiltonians are able to describe both the ground and the excited states of molecular systems in a balanced manner. Using the variational quantum eigensolver (VQE) method based on the unitary coupled cluster with singles and doubles (UCCSD) ansatz and only a minimal basis set (ANO-RCC-MB), we demonstrate that the transcorrelated Hamiltonians can produce ground state energies comparable to the reference CCSD energies with the much larger cc-pVTZ basis set. This leads to a reduction in the number of required CNOT gates by more than 3 orders of magnitude for the chemical species studied in this work. Furthermore, using the quantum equation of motion (qEOM) formalism in conjunction with the transcorrelated Hamiltonian, we are able to reduce the deviations in the excitation energies from the reference EOM-CCSD/cc-pVTZ values by an order of magnitude. The transcorrelated Hamiltonians developed here are Hermitian and contain only one- and two-body interaction terms and thus can be easily combined with any quantum algorithm for accurate electronic structure simulations.
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.