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.
The opportunities afforded by near-term quantum computers to calculate the ground-state properties of small molecules depend on the structure of the computational ansatz as well as the errors induced by device noise. Here we investigate the behavior of these noisy quantum circuits using numerical simulations to estimate the accuracy and fidelity of the prepared quantum states relative to the ground truth obtained by conventional means. We implement several different types of ansatz circuits derived from unitary coupled cluster theory for the purposes of estimating the ground-state energy of sodium hydride using the variational quantum eigensolver algorithm. We show how relative error in the energy and the fidelity scale with the levels of gate-based noise, the internuclear configuration, the ansatz circuit depth, and the parameter optimization methods.
In recent years, the Variational Quantum Eigensolver (VQE) has emerged as one of the most popular algorithms for solving the electronic structure problem on near-term quantum computers. The utility of VQE is often hindered by the limitations of current quantum hardware, including, but not limited to, short qubit coherence times and low gate fidelities. These limitations become particularly pronounced when VQE is used along with deep quantum circuits, such as those required by the ”Unitary Coupled Cluster Singles and Doubles” (UCCSD) ansatz, often resulting in significant errors. To address these issues, we propose a low-depth ansatz based on parallelized Givens rotations, which can recover substantial correlation energy while drastically reducing circuit depth and two-qubit gate counts for an arbitrary active space (AS). Also, considering the current hardware architectures with low qubit counts, we introduce a systematic way to select molecular orbitals to define active spaces (ASs) that retain significant electron correlation. We validate our approach by computing bond dissociation profiles of water and strongly correlated systems, such as molecular nitrogen and oxygen, across various ASs. Noiseless simulations using the new ansatz yield ground-state energies comparable to those from the UCCSD ansatz while reducing circuit depth by 50–70%. Moreover, in noisy simulations, our approach achieves energy error rates an order of magnitude lower than that of UCCSD. Considering the efficiency and practical usage of our ansatz, we hope that it becomes a potential choice for performing quantum chemistry calculations on near-term quantum devices.
While numerical simulations are presented in most papers introducing new methods to enhance the VQE performance, comprehensive, comparative, and applied studies remain relatively rare. We present a comprehensive, yet concise guide for the implementation of the VQE for molecular problems on NISQ devices, specifically applied to estimate the ground state energy of the BeH2 molecule using hardware-efficient and chemically informed ansätze. This work clarifies several under-documented aspects in the literature, such as the construction of the electronic Hamiltonian, the transformation of fermionic operators into qubit operators via second quantization, and the mathematical framework's details for the unitary coupled cluster single and double (UCCSD) ansatz. Our methodology, implemented using Qiskit 1.2, the latest release as of the date of this writing, is demonstrated on a noiseless simulator and further tested with noisy quantum circuits. The resilience of the VQE to quantum noise remains an open question. This study compares the computational accuracy of ground state energy estimations for molecules using the VQE across three different current quantum hardware noise models. Furthermore, our experiment on IBM's 156-qubit actual quantum computer revealed valuable insights on the real performance of the VQE on current quantum hardware.
Quantum computing has emerged as a transformative approach for tackling complex problems in quantum chemistry, particularly in simulating multielectron systems and electron-phonon interactions. However, the current noisy intermediate-scale quantum (NISQ) devices face significant challenges in implementing practical quantum algorithms due to error accumulation caused by increased circuit depth, qubit counts, and gate operations. To address these challenges, we present TenCirChem, an open-source Python package specifically designed for implementing variational quantum algorithms in quantum computational chemistry. TenCirChem demonstrates exceptional performance in simulating unitary coupled-cluster circuits through its innovative use of compact representations for quantum states and excitation operators. This package supports noisy circuit simulation and provides advanced algorithms for variational quantum dynamics, enabling researchers to explore complex chemical phenomena. Its capabilities are exemplified in various applications, including the calculation of potential energy curves and the investigation of quantum gate error impacts on molecular systems. Moreover, TenCirChem’s seamless integration with real quantum hardware makes it a versatile tool for both simulation and experimentation. A key innovation developed within the TenCirChem framework is the Clifford-based Hamiltonian engineering approach for molecules (CHEM). This algorithm addresses the critical challenge of achieving chemical accuracy with shallow quantum circuits, a fundamental requirement for practical applications on NISQ devices. CHEM employs a sophisticated Clifford-based Hamiltonian transformation that operates within the variational quantum eigensolver (VQE) framework using hardware-efficient ansatz. The method ensures four crucial advantages: (1) generation of initial circuit parameters corresponding to Hartree-Fock energy, (2) maximization of initial energy gradients with respect to circuit parameters, (3) minimal classical processing overhead without additional quantum resource requirements, and (4) compatibility with any circuit topology. Through quantum hardware emulator demonstrations, CHEM has achieved chemical accuracy for systems up to 12 qubits with fewer than 30 two-qubit gates, representing a significant advancement in practical quantum computational chemistry. To enhance the efficiency of variational quantum algorithms, we developed the sequential optimization with approximate parabola (SOAP) method, specifically designed for parameter optimization in unitary coupled-cluster ansatz. SOAP addresses the critical bottleneck of measurement requirements in VQE by implementing an innovative optimization strategy that approximates the energy landscape as quadratic functions. This approach minimizes the number of energy evaluations while incorporating parameter correlations through the integration of average directions from previous iterations. Benchmark studies demonstrate SOAP’s superior performance, showing faster convergence and enhanced noise robustness compared to traditional optimization methods. The method’s scalability has been validated through numerical simulations of up to 20 qubits, and its practical efficacy has been confirmed through experiments on superconducting quantum computers. For simulating electron-phonon systems, we introduce a variational basis state encoding algorithm that significantly reduces resource requirements compared to conventional unary and binary encoding schemes. Our approach achieves smaller scaling than traditional methods for qubits and gates for systems obeying the area law of entanglement entropy, this remarkable reduction in resource requirements comes at the cost of a constant amount of additional measurements. The algorithm’s effectiveness has been validated through both numerical simulations and quantum hardware experiments, demonstrating that using just one or two qubits per phonon mode can produce quantitatively accurate results across various coupling regimes. The integration of these innovations within the TenCirChem library represents a significant advancement in quantum computational chemistry. The software package provides researchers with a comprehensive toolkit for developing, testing, and implementing quantum algorithms, while the novel methods address fundamental challenges in circuit depth, parameter optimization, and resource efficiency. These developments collectively enhance the practicality of quantum computing for addressing real-world chemical problems in the NISQ era, offering improved accuracy, efficiency, and noise robustness.
molecule, emphasizing practical implementation and performance on current quantum hardware. Our research presents a comparative study of HEA and UCCSD ansätze on noiseless and noisy simulations and implements VQE on recent IBM quantum computer noise models and a real quantum computer, IBM Fez, providing a fully functional code employing Qiskit 1.2. Our experiments confirm UCCSD's reliability in ideal conditions, while the HEA demonstrates greater robustness to hardware noise, achieving chemical accuracy on state-vector simulation (SVS). The results reveal that achieving ground-state energy within chemical accuracy is feasible without error mitigation during VQE convergence. We demonstrate that current quantum devices effectively optimize circuit parameters despite misestimating simulated energies. The SVS-evaluated energies provide a more accurate representation of the solution quality compared to QPU-estimated energy values, indicating that VQE converges to the correct ground state despite quantum noise. Our study also applies noise mitigation as a postprocessing technique, using zero-noise extrapolation (ZNE) on a real quantum computer. The detailed methodologies presented in this study, including Hamiltonian construction and Fermionic-to-qubit transformations, facilitate flexible adaptation of the VQE approach for various algorithm variants and across different levels of algorithmic implementation.