arXiv:2310.01302 · Journal of Chemical Theory and Computation
Toward Accurate Post-Born–Oppenheimer Molecular Simulations on Quantum Computers: An Adaptive Variational Eigensolver with Nuclear-Electronic Frozen Natural Orbitals
Anton Nykänen, Aaron Miller, N. Walter Talarico, Stefan Knecht, Arseny Kovyrshin, Mårten Skogh, Lars Tornberg, Anders Broo, Stefano Mensa, Benjamin C. B. Symons, Emre Sahin, Jason N. Crain, Ivano Tavernelli, Fabijan Pavošević
Nuclear quantum effects such as zero-point energy and hydrogen tunneling play a central role in many biological and chemical processes. The nuclear-electronic orbital (NEO) approach captures these effects by treating selected nuclei quantum mechanically on the same footing as electrons. On classical computers, the resources required for an exact solution of NEO-based models grow exponentially with system size. By contrast, quantum computers offer a means of solving this problem with polynomial scaling. However, due to the limitations of current quantum devices, NEO simulations are confined to the smallest systems described by minimal basis sets, whereas realistic simulations beyond the Born–Oppenheimer approximation require more sophisticated basis sets. For this purpose, we herein extend a hardware-efficient ADAPT-VQE method to the NEO framework in the frozen natural orbital (FNO) basis. We demonstrate on H 2 and D 2 molecules that the NEO-FNO-ADAPT-VQE method reduces the CNOT count by several orders of magnitude relative to the NEO unitary coupled cluster method with singles and doubles while maintaining the desired accuracy. This extreme reduction in the CNOT gate count is sufficient to permit practical computations employing the NEO method─an important step toward accurate simulations involving nonclassical nuclei and non-Born–Oppenheimer effects on near-term quantum devices. We further show that the method can capture isotope effects, and we demonstrate that inclusion of correlation energy systematically improves the prediction of difference in the zero-point energy (ΔZPE) between isotopes.
Nuclear quantum effects such as zero-point energy and hydrogen tunneling play a central role in many biological and chemical processes. The nuclear-electronic orbital (NEO) approach captures these effects by treating selected nuclei quantum mechanically on the same footing as electrons. On classical computers, the resources required for an exact solution of NEO-based models grow exponentially with system size. By contrast, quantum computers offer a means of solving this problem with polynomial scaling. However, due to the limitations of current quantum devices, NEO simulations are confined to the smallest systems described by minimal basis sets, whereas realistic simulations beyond the Born–Oppenheimer approximation require more sophisticated basis sets. For this purpose, we herein extend a hardware-efficient ADAPT-VQE method to the NEO framework in the frozen natural orbital (FNO) basis. We demonstrate on H 2 and D 2 molecules that the NEO-FNO-ADAPT-VQE method reduces the CNOT count by several orders of magnitude relative to the NEO unitary coupled cluster method with singles and doubles while maintaining the desired accuracy. This extreme reduction in the CNOT gate count is sufficient to permit practical computations employing the NEO method─an important step toward accurate simulations involving nonclassical nuclei and non-Born–Oppenheimer effects on near-term quantum devices. We further show that the method can capture isotope effects, and we demonstrate that inclusion of correlation energy systematically improves the prediction of difference in the zero-point energy (ΔZPE) between isotopes.
Nuclear quantum phenomena beyond the Born-Oppenheimer approximation are known to play an important role in a growing number of chemical and biological processes. While there exists no unique consensus on a rigorous and efficient implementation of coupled electron-nuclear quantum dynamics, it is recognized that these problems scale exponentially with system size on classical processors and, therefore, may benefit from quantum computing implementations. Here, we introduce a methodology for the efficient quantum treatment of the electron-nuclear problem on near-term quantum computers, based upon the Nuclear-Electronic Orbital (NEO) approach. We generalize the electronic two-qubit tapering scheme to include nuclei by exploiting symmetries inherent in the NEO framework, thereby reducing the Hamiltonian dimension, number of qubits, gates, and measurements needed for calculations. We also develop parameter transfer and initialization techniques, which improve convergence behavior relative to conventional initialization. These techniques are applied to H2 and malonaldehyde for which results agree with NEO full configuration interaction and NEO complete active space configuration interaction benchmarks for ground state energy to within 10-6 hartree and entanglement entropy to within 10-4. These implementations therefore significantly reduce resource requirements for full quantum simulations of molecules on near-term quantum devices while maintaining high accuracy.
The variational quantum eigensolver (VQE) algorithm combined with the unitary coupled cluster (UCC) ansatz has been developed for the quantum computation of molecular energies and wave functions within the Born–Oppenheimer approximation. Herein, this approach is extended to multicomponent systems to enable the quantum mechanical treatment of more than one type of particle, such as electrons and positrons or electrons and nuclei, without invoking the Born–Oppenheimer approximation. Specifically, we introduce the multicomponent unitary coupled cluster (mcUCC) method combined with the VQE algorithm for the calculation of ground-state energies and wave functions as well as the multicomponent equation-of-motion (mcEOM) method for the calculation of excitation energies. These methods are developed within the nuclear-electronic orbital (NEO) framework and are formulated in the qubit basis to enable implementations on quantum computers. Moreover, these methods are used to calculate the ground-state energy and excitation energies of positronium hydride, where both electrons and the positron are treated quantum mechanically, as well as the H 2 molecule, where both electrons and one proton are treated quantum mechanically. These applications validate the implementation and provide benchmark data for future calculations. The errors due to Trotterization of the mcUCC ansatz are also analyzed. This formalism, as well as the accompanying computer code, will serve as the basis for applications to more complex multicomponent systems, such as simulations of photoinduced nonadiabatic molecular processes, on both classical and quantum computers.
The method of increments and frozen natural orbital (MI-FNO) framework is introduced to help expedite the application of noisy, intermediate-scale quantum (NISQ) devices for quantum chemistry simulations. The MI-FNO framework provides a systematic reduction of the occupied and virtual orbital spaces for quantum chemistry simulations. The correlation energies of the resulting increments from the MI-FNO reduction can then be solved by various algorithms, including quantum algorithms such as the phase estimation algorithm and the variational quantum eigensolver (VQE). The unitary coupled-cluster singles and doubles VQE framework is used to obtain correlation energies for the case of small molecules (i.e., BeH2, CH4, NH3, H2O, and HF) using the cc-pVDZ basis set. The quantum resource requirements are estimated for a constrained geometry complex catalyst that is utilized in industrial settings for the polymerization of α-olefins. We show that the MI-FNO approach provides a significant reduction in the quantum bit (qubit) requirements relative to the full system simulations. We propose that the MI-FNO framework can create scalable examples of quantum chemistry problems that are appropriate for assessing the progress of NISQ devices.
The variational quantum eigensolver (VQE) has emerged as one of the leading quantum algorithms for solving electronic structure problems on near-term noisy intermediate-scale quantum devices. However, its practical application to quantum chemistry remains challenging due to the limited coherence time, imperfect quantum gate fidelity, and the large number of measurements required, which together confine current electronic structure simulations to relatively small active spaces. In this work, we present an orbital-optimized VQE framework based on orbital compression, designed to improve the accuracy of electronic structure calculations while maintaining relatively small active spaces. Frozen natural orbitals (FNO) and split virtual orbitals (SVO) are first employed to construct compact active spaces for VQE simulations, leading to the FNO/SVO-VQE approach. Orbital optimization is then incorporated to further recover electron correlation effects, resulting in the FNO/SVO-OO-VQE methods. We apply the proposed method to simulate potential energy surfaces for molecular dissociation and the activation energy of formaldehyde decomposition. Numerical results demonstrate that both FNO-OO-VQE and SVO-OO-VQE improve the variational accuracy while substantially reducing measurement cost.
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.
Quantum computing offers promising advantages for computational chemistry, particularly through algorithms that efficiently model correlated methods. In this work, we apply quantum computing techniques to the constrained nuclear-electronic orbital (CNEO) framework, which enables the inclusion of nuclear quantum effects in chemical simulations while preserving a well-defined molecular structure. We present the development and implementation of two correlated wave function methods within this framework: CNEO full configuration interaction (CNEO-FCI) and CNEO unitary coupled-cluster with singles and doubles (CNEO-UCCSD), with the latter solved using the variational quantum eigensolver algorithm. These methods were applied to the hydrogen isotopologues H 2, HD, and D 2 and used to calculate potential energy surfaces, equilibrium geometries, harmonic vibrational frequencies, and von Neumann entropies. The CNEO-UCCSD results are in excellent agreement with CNEO-FCI, recovering over 99% of the correlation energy and accurately capturing geometries and vibrational frequencies. Additionally, we observe a strong connection between CNEO-FCI and CNEO-UCCSD energy and entropy differences as a function of bond length, highlighting the role of quantum entanglement in molecular dissociation. These results demonstrate the viability of CNEO-based quantum algorithms for capturing nuclear quantum effects and lay the groundwork for future quantum simulations and entropy analysis of multicomponent systems within the CNEO paradigm.