The variational quantum eigensolver (VQE) is a widely employed method to solve electronic structure problems in the current noisy intermediate-scale quantum (NISQ) devices. However, due to inherent noise in the NISQ devices, VQE results on NISQ devices often deviate significantly from the results obtained on noiseless statevector simulators or traditional classical computers. The iterative nature of VQE further amplifies the errors in each loop. Recent works have explored ways to integrate deep neural networks (DNN) with VQE to mitigate iterative errors, albeit primarily limited to the noiseless statevector simulators. In this work, we trained DNN models across various quantum circuits and examined the potential of two DNN-VQE approaches, DNN1 and DNNF, for predicting the ground state energies of small molecules in the presence of device noise. We carefully examined the accuracy of the DNN1, DNNF, and VQE methods on both noisy simulators and real quantum devices by considering different ansatzes of varying qubit counts and circuit depths. Our results illustrate the advantages and limitations of both VQE and DNN-VQE approaches. Notably, both DNN1 and DNNF methods consistently outperform the standard VQE method in providing more accurate ground state energies in noisy environments. However, despite being more accurate than VQE, the energies predicted using these methods on real quantum hardware remain meaningful only at reasonable circuit depths (depth = 15, gates = 21). At higher depths (depth = 83, gates = 112), they deviate significantly from the exact results. Additionally, we find that DNNF does not offer any notable advantage over VQE in terms of speed. Consequently, our study recommends DNN1 as the preferred method for obtaining quick and accurate ground state energies of molecules on current quantum hardware, particularly for quantum circuits with lower depth and fewer qubits.
Quantum computing offers a potential for algorithmic speedups for applications, such as large-scale simulations in chemistry and physics. However, these speedups must yield results that are sufficiently accurate to predict realistic outcomes of experiments precisely. Delivering on the promise of high accuracy and precision requires methods to evaluate the computational accuracy of the quantum computing devices. We develop a framework to estimate the computational accuracy of near-term noisy, intermediate scale quantum (NISQ) computing devices using a quantum chemistry application. Application benchmarks that run on NISQ devices require techniques for mitigating errors to improve accuracy and precision. We use device agnostic error-mitigation schemes, quantum error detection and readout error detection, with post-selection to mitigate the dominant sources of noise. We evaluate the framework by simulating the ground state of molecular hydrogen with the variational quantum eigensolver (VQE) algorithm, estimating the energy and calculating the precision of the estimate using numerical simulations with realistic noise models. We first quantify the improvement in the logical error rate and state fidelity of the VQE application when encoded with the [[4,2,2]] quantum error detection code. When additionally encoded with readout error detection, we show that compared to the unencoded simulation, the encoded simulation yields a more accurate estimate by more than 1 mHa (0.027 eV) with comparable precision and higher state fidelity. Additionally, unlike the best estimate from the unencoded simulations, the results from the encoded simulation fall within the chemical accuracy threshold of 1.6 mHa of the exact energy. The estimated accuracy and precision indicate that current quantum computers can achieve error rates that yield useful outcomes for chemical applications.
Quantum chemistry is one of the most promising applications of quantum computers in the near future. For noisy intermediate-scale quantum devices, the quantum-classical hybrid framework based on the variational quantum eigensolver (VQE) has become the method of choice. In the literature, there are many different variants of VQE, but it is not known which one is optimal for a given molecule. For this purpose, we perform a thorough benchmarking on more than ten different kinds of VQE ansatzes (in systems up to 30 qubits), based on their performance on the energy accuracy, runtime until convergence, and number of parameters. Our results show that the ADAPT ansatz can be used to obtain more accurate energy for small systems (below 14 qubits), but it costs much more computational resources. For larger molecules, UCCSD0 has better performance. However, all the tested ansatzes can hardly reach chemical accuracy at stretched bond lengths. Our results were obtained using MindSpore Quantum, where the codes and the benchmarking toolkit are publicly available at Gitee.
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.
We present a self-consistent field (SCF) approach within the adaptive derivative-assembled problem-tailored ansatz variational quantum eigensolver (ADAPT-VQE) framework for efficient quantum simulations of chemical systems on near-term quantum computers. To this end, our ADAPT-VQE-SCF approach combines the idea of generating an ansatz with a small number of parameters, resulting in shallow-depth quantum circuits with a direct minimization of an energy expression that is correct to second order with respect to changes in the molecular orbital basis. Our numerical analysis, including calculations for the transition-metal complex ferrocene [Fe (C 5 H 5 ) 2 ], indicates that convergence in the self-consistent orbital optimization loop can be reached without a considerable increase in the number of two-qubit gates in the quantum circuit by comparison to a VQE optimization in the initial molecular orbital basis. Moreover, the orbital optimization can be carried out simultaneously within each iteration of the ADAPT-VQE cycle. ADAPT-VQE-SCF thus allows us to implement a routine analogous to the complete active space SCF, a cornerstone of state-of-the-art computational chemistry, in a hardware-efficient manner on near-term quantum computers. Hence, ADAPT-VQE-SCF paves the way toward a paradigm shift for quantitative quantum-chemistry simulations on quantum computers by requiring fewer qubits and opening up for the use of large and flexible atomic orbital basis sets in contrast to earlier methods that are predominantly based on the idea of full active spaces with minimal basis sets.
The Variational Quantum Eigensolver (VQE) is widely regarded as a promising algorithm for calculating ground states of quantum systems that are intractable for classical computers. This promise is typically motivated by the hope of mitigating the exponential growth of Hilbert space with system size. Here we scrutinize how the computational cost of adaptive VQE scales with the size of the target system. We demonstrate that the Rényi entropy derived from classical simulations predicts the required number of adaptive iterations of VQE with high accuracy (
R2≈0.99). We validate this on a benchmarking set of more than 20 different molecules with active spaces ranging from four to ten orbitals. For these molecules, we find an exponential scaling of the number of adaptive iterations, and in turn, of the circuit depth with the system size. We therefore conclude that it is unlikely that VQE in its current form is able to simulate large molecular systems with high fidelity without exponential resource requirements.