A key requirement to perform simulations of large quantum systems on near-term quantum hardware is the design of quantum algorithms with short circuit depth that finish within the available coherence time. A way to stay within the limits of coherence is to reduce the number of gates by implementing a gate set that matches the requirements of the specific algorithm of interest directly in hardware. Here, we show that exchange-type gates are a promising choice for simulating molecular eigenstates on near-term quantum devices since these gates preserve the number of excitations in the system. Complementing the theoretical work by Barkoutsos et al. [PRA 98, 022322 (2018)], we report on the experimental implementation of a variational algorithm on a superconducting qubit platform to compute the eigenstate energies of molecular hydrogen. We utilize a parametrically driven tunable coupler to realize exchange-type gates that are configurable in amplitude and phase on two fixed-frequency superconducting qubits. With gate fidelities around 95% we are able to compute the eigenstates within an accuracy of 50 mHartree on average, a limit set by the coherence time of the tunable coupler.
High Resolution Image Download MS PowerPoint Slide The recent developments of quantum computing present novel potential pathways for quantum chemistry as the scaling of the computational power of quantum computers could be harnessed to naturally encode and solve electronic structure problems. Theoretically exact quantum algorithms for chemistry have been proposed (e.g., quantum phase estimation), but the limited capabilities of current noisy intermediate-scale quantum devices motivated the development of less demanding hybrid algorithms. In this context, the variational quantum eigensolver (VQE) algorithm was successfully introduced as an effective method to compute the ground-state energies of small molecules. This study investigates the folded spectrum (FS) method as an extension of the VQE algorithm for the computation of molecular excited states. It provides the possibility of directly computing excited states around a selected target energy using the same or quantum circuit as for the ground-state calculation. Inspired by the variance-based methods from the quantum Monte Carlo literature, the FS method minimizes the energy variance, thus, in principle, requiring a computationally expensive squared Hamiltonian to be applied. We alleviate this potentially poor scaling by employing a Pauli grouping procedure to identify sets of commuting Pauli strings that can be evaluated simultaneously. This allows for a significant reduction in the computational cost. We applied the FS-VQE method to small molecules (H 2, LiH), obtaining all electronic excited states with chemical accuracy on ideal quantum simulators. Furthermore, we explore the application of quantum error mitigation techniques, demonstrating improved energy accuracy on noisy simulators compared with simulations without mitigation.
Variational quantum algorithms are emerging as promising candidates for near-term practical applications of quantum information processors, in the field of quantum chemistry. We implement the variational quantum eigensolver algorithm to calculate the molecular ground-state energy of the H2 molecule and experimentally demonstrated it on an NMR quantum processor. Further, we simulate the excited states of the H2 molecule using the variational quantum deflation algorithm and experimentally demonstrate it on the same NMR quantum processor. We also develop the first simulation of the energy calculation of the H2 molecule using only a single qubit, and verify the results on an NMR quantum computer. Our experimental results demonstrate that only a single NMR qubit suffices to calculate the molecular energies of the H2 molecule to the desired accuracy.
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.
New approaches into computational quantum chemistry can be developed through the use of quantum computing. While universal, fault-tolerant quantum computers are still not available, and we want to utilize today's noisy quantum processors. One of their flagship applications is the variational quantum eigensolver (VQE)-an algorithm for calculating the minimum energy of a physical Hamiltonian. In this study, we investigate how various types of errors affect the VQE and how to efficiently use the available resources to produce precise computational results. We utilize a simulator of a noisy quantum device, an exact statevector simulator, and physical quantum hardware to study the VQE algorithm for molecular hydrogen. We find that the optimal method of running the hybrid classical-quantum optimization is to: (i) allow some noise in intermediate energy evaluations, using fewer shots per step and fewer optimization iterations, but ensure a high final readout precision; (ii) emphasize efficient problem encoding and ansatz parametrization; and (iii) run all experiments within a short time-frame, avoiding parameter drift with time. Nevertheless, current publicly available quantum resources are still very noisy and scarce/expensive, and even when using them efficiently, it is quite difficult to perform trustworthy calculations of molecular energies.
Quantum mechanics has introduced a new theoretical framework for the study of molecules, enabling the prediction of properties and dynamics through the solution of the Schrödinger equation applied to these systems. However, solving this equation is computationally expensive, which has led to the development of various mathematical frameworks and computational methods designed to balance the available resources with the desired level of accuracy. In particular, quantum computers have emerged as a promising technology with the potential to address these problems more efficiently in the coming decades, whether through reductions in memory, time, and energy consumption,
i. e., reductions in computational complexity or by enhancing precision. This research field is known as Quantum Simulation. Given the current technological limitations of quantum computers, Variational Quantum Algorithms (VQAs), especially the Variational Quantum Eigensolver (VQE), stand out as a feasible approach to demonstrating advantages over classical methods in the near term. This feasibility arises from their lower demand for quantum gates and the reduced depth of the circuits required for their implementation. In this work, we present the application of quantum mechanics to the study of molecules, provide an introduction to the fundamentals of quantum computing, and explore the integration of these fields by employing the VQE in molecular simulations. Finally, we discuss the spatial and temporal complexity associated with the algorithm, highlighting its implications and challenges.