During the noisy intermediate-scale quantum (NISQ) era, quantum computational approaches refined to overcome the challenge of limited quantum resources are highly valuable. A comprehensive benchmark for a quantum computational approach in this spirit could provide insights toward further improvements. On the other hand, the accuracy of the molecular properties predicted by most of the quantum computations nowadays is still far off (not within chemical accuracy) compared to their corresponding experimental data. In this work, we propose a promising qubit-efficient quantum computational approach and present a comprehensive investigation by benchmarking quantum computation of the harmonic vibrational frequencies of a large set of neutral closed-shell diatomic molecules with results in great agreement with their experimental data. To this end, we construct the accurate Hamiltonian using molecular orbitals, derived from density functional theory to account for the electron correlation and expanded in the Daubechies wavelet basis set to allow an accurate representation in real space grid points, where an optimized compact active space is further selected so that only a reduced small number of qubits is sufficient to yield an accurate result. Typically, calculations achieved with 2 to 12 qubits using our approach would need 20 to 60 qubits using a traditional cc-pVDZ basis set with frozen core approximation to achieve similar accuracy. To justify the approach, we benchmark the performance of the Hamiltonians spanned by the selected molecular orbitals by first transforming the molecular Hamiltonians into qubit Hamiltonians and then using the exact diagonalization method to calculate the results, regarded as the best results achievable by quantum computation to compare to the experimental data. Furthermore, using the variational quantum eigensolver algorithm with the constructed qubit Hamiltonians, we show that the variational quantum circuit with the chemistry-inspired UCCSD ansatz can achieve the same accuracy as the exact diagonalization method except for systems whose Mayer bond order indices are larger than 2. For those systems, we then demonstrate that the heuristic hardware-efficient RealAmplitudes ansatz, even with a substantially shorter circuit depth, can provide a significant improvement over the UCCSD ansatz, verifying that the harmonic vibrational frequencies could be calculated accurately by quantum computation in the NISQ era.
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.
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.
High Resolution Image Download MS PowerPoint Slide The variational quantum eigensolver (VQE) algorithm has positioned quantum chemistry calculation methods as promising applications for noisy intermediate-scale quantum (NISQ) devices. In NISQ devices, where the quantum noise is considerable, it is desirable to have an ansatz for the quantum circuit that is practical to implement and capable of achieving high chemical accuracy. Hardware-efficient Ansätze (HEA), such as the RyRz linear ansatz (RLA) and its modified form that preserves physical quantities─the symmetry preserving ansatz (SPA)─are constructed using quantum gates that can be implemented easily on a quantum computer. However, the extent to which HEA can accurately yield the electronic state energies of actual molecules and the number of quantum gates required to achieve chemical accuracy─which, in practical cases, should be within 1 kcal/mol of the exact energy value─are not well understood. In this study, we aimed to gain a detailed understanding of these aspects by performing noiseless simulations to obtain the ground and low-lying excited state using high-depth HEA quantum circuits for several molecules, including LiH, H 2 O, BeH 2, CH 4, and N 2 . From the results, we demonstrated that symmetry preserving HEA, such as SPA, can achieve accurate computational results that maintain CCSD-level chemical accuracy by increasing the number of layers. We quantitatively analyze how the expressibility and entangling capability of RLA and SPA quantum circuits evolve with circuit depth. Our results reveal that while increased layers generally expand accessible Hilbert space and entanglement, differences in their limitations highlight the importance of considering physically allowed entanglement for performance of a parametric quantum circuit in VQE. We also studied the potential energy surface of the dissociation of these molecules and found that quantum circuits using SPA can capture static electron correlation, which is challenging to address with classical single-reference quantum chemistry methods such as CCSD. These results from SPA demonstrate its potential to represent electronic states, as it achieves highly accurate results with fewer gate operations compared to physically inspired quantum circuits such as the Unitary Coupled Cluster (UCC) method. This study establishes a possible path for the application of quantum devices in solving quantum chemistry problems.
CO2 capture is critical to solving global warming. Amine-based solvents are extensively used to chemically absorb CO2. Thus, it is crucial to study the chemical absorption of CO2 by amine-based solvents to better understand and optimize CO2 capture processes. Here, we use quantum computing algorithms to quantify molecular vibrational energies and reaction pathways between CO2 and a simplified amine-based solvent model—NH3. Molecular vibrational properties are important to understanding kinetics of reactions. However, the molecule size correlates with the strength of anharmonicity effect on vibrational properties, which can be challenging to address using classical computing. Quantum computing can help enhance molecular vibrational calculations by including anharmonicity. We implement a variational quantum eigensolver (VQE) algorithm in a quantum simulator to calculate ground state vibrational energies of reactants and products of the CO2 and NH3 reaction. The VQE calculations yield ground vibrational energies of CO2 and NH3 with similar accuracy to classical computing. In the presence of hardware noise, Compact Heuristic for Chemistry (CHC) ansatz with shallower circuit depth performs better than Unitary Vibrational Coupled Cluster. The “Zero Noise Extrapolation” error-mitigation approach in combination with CHC ansatz improves the vibrational calculation accuracy. Excited vibrational states are accessed with quantum equation of motion method for CO2 and NH3. Using quantum Hartree–Fock (HF) embedding algorithm to calculate electronic energies, the corresponding reaction profile compares favorably with Coupled Cluster Singles and Doubles while being more accurate than HF. Our research showcases quantum computing applications in the study of CO2 capture reactions.
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.