We present the first step toward the quantum computing (QC) formulation of the electron nuclear dynamics (END) method within the variational quantum simulator (VQS) scheme: END/QC/VQS. END is a time-dependent, variational, on-the-flight, and non-adiabatic method to simulate chemical reactions. END represents nuclei with frozen Gaussian wave packets and electrons with a single-determinantal state in the Thouless non-unitary representation. Within the hybrid quantum/classical VQS, END/QC/VQS currently evaluates the metric matrix M and gradient vector V of the symplectic END/QC equations on the QC software development kit QISKIT, and calculates basis function integrals and time evolution on a classical computer. To adapt END to QC, we substitute the Thouless non-unitary representation with Fukutome unitary representation. We derive the first END/QC/VQS version for pure electronic dynamics in multielectron chemical models consisting of two-electron units with fixed nuclei. Therein, Fukutome unitary matrices factorize into triads of one-qubit rotational matrices, which leads to a QC encoding of one electron per qubit. We design QC circuits to evaluate M and V in one-electron diatomic molecules. In log2-log2 plots, errors and deviations of those evaluations decrease linearly with the number of shots and with slopes = −1/2. We illustrate an END/QC/VQS simulation with the pure electronic dynamics of H2+ We discuss the present results and future END/QC/QVS extensions.
We present the first step toward the quantum computing (QC) formulation of the electron nuclear dynamics (END) method within the variational quantum simulator (VQS) scheme: END/QC/VQS. END is a time-dependent, variational, on-the-flight, and non-adiabatic method to simulate chemical reactions. END represents nuclei with frozen Gaussian wave packets and electrons with a single-determinantal state in the Thouless non-unitary representation. Within the hybrid quantum/classical VQS, END/QC/VQS currently evaluates the metric matrix M and gradient vector V of the symplectic END/QC equations on the QC software development kit QISKIT, and calculates basis function integrals and time evolution on a classical computer. To adapt END to QC, we substitute the Thouless non-unitary representation with Fukutome unitary representation. We derive the first END/QC/VQS version for pure electronic dynamics in multielectron chemical models consisting of two-electron units with fixed nuclei. Therein, Fukutome unitary matrices factorize into triads of one-qubit rotational matrices, which leads to a QC encoding of one electron per qubit. We design QC circuits to evaluate M and V in one-electron diatomic molecules. In log2-log2 plots, errors and deviations of those evaluations decrease linearly with the number of shots and with slopes = −1/2. We illustrate an END/QC/VQS simulation with the pure electronic dynamics of H2+ We discuss the present results and future END/QC/QVS extensions.
Simulations of chemical dynamics are a powerful means for understanding chemistry. However, classical computers struggle to simulate many chemical processes, especially non-adiabatic ones, where the Born-Oppenheimer approximation breaks down. Quantum computers could simulate quantum-chemical dynamics more efficiently than classical computers, but there is currently no complete quantum algorithm for calculating dynamical observables to within a known error. Here, we develop an efficient, end-to-end quantum algorithm for simulating chemical dynamics that avoids all uncontrolled approximations (including the Born-Oppenheimer approximation) and whose error is bounded subject to mild assumptions. To do so, we treat the nuclei and the electrons on an equal footing and simulate the full molecular wavefunction on a momentum-space grid in first quantization, including all algorithmic steps: initial-state preparation, time evolution using qubitization, and measurement of chemical observables such as reaction yields and rates. Our work gives the first algorithm for quantum simulation of chemistry whose end-to-end complexity achieves sublinear scaling in the size of the grid. We achieve this by developing an exponentially faster method for initial-state-preparation. Photochemistry is a likely early application of our algorithm and we estimate resources required for end-to-end simulations of non-adiabatic dynamics of atmospherically important molecules. Classically intractable photochemical computations could be performed using resources comparable to those required for other chemical applications of quantum computing.
The electronic energies of molecules have been successfully evaluated on quantum computers. However, more attention is paid to the dynamics simulation of molecules in practical applications. Based on the variational quantum eigensolver (VQE) algorithm, Fedorov et al. proposed a correlated sampling (CS) method and demonstrated the vibrational dynamics of H2 molecules [J. Chem. Phys. 154, 164103 (2021)]. In this study, we have developed a quantum approach by extending the CS method based on the VQE algorithm (labeled eCS-VQE) for simulating chemical reaction dynamics. First, the CS method is extended to the three-dimensional cases for calculation of first-order energy gradients, and then, it is further generalized to calculate the second-order gradients of energies. By calculating atomic forces and vibrational frequencies for H2, LiH, H+ + H2, and Cl- + CH3Cl systems, we have seen that the approach has achieved the CCSD level of accuracy. Thus, we have simulated dynamics processes for two typical chemical reactions, hydrogen exchange and chlorine substitution, and obtained high-precision reaction dynamics trajectories consistent with the classical methods. Our eCS-VQE approach, as measurement expectations and ground-state wave functions can be reused, is less demanding in quantum computing resources and is, therefore, a feasible means for the dynamics simulation of chemical reactions on the current noisy intermediate-scale quantum-era quantum devices.
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.
ConspectusSimulating molecular dynamics (MD) within a comprehensive quantum framework has been a long-standing challenge in computational chemistry. An exponential scaling of computational cost renders solving the time dependent Schrödinger equation (TDSE) of a molecular Hamiltonian, including both electronic and nuclear degrees of freedom (DOFs), as well as their couplings, infeasible for more than a few DOFs. In the Born-Oppenheimer (BO), or adiabatic, picture, electronic and nuclear parts of the wave function are decoupled and treated separately. Within this framework, the nuclear wave function evolves along potential energy surfaces (PESs) computed as solutions to the electronic Schrödinger equation parametrized in the nuclear DOFs. This approximation, together with increasingly elaborate numerical approaches to solve the nuclear time dependent Schrödinger equation (TDSE), enabled the treatment of up to a few dozens of degrees of freedom (DOFs). However, for particular applications, such as photochemistry, the BO approximation breaks down. In this regime of non-adiabatic dynamics, solving the full molecular problem including electron-nuclear couplings becomes essential, further increasing the complexity of the numerical solution. Although valuable methods such as multiconfigurational time-dependent Hartree (MCTDH) have been proposed for the solution of the coupled electron-nuclear dynamics, they remain hampered by an exponential scaling in the number of nuclear DOFs and by the difficulty of finding universal variational forms.In this Account, we present a perspective on novel quantum computational algorithms, aiming to alleviate the exponential scaling inherent to the simulation of many-body quantum dynamics. In particular, we focus on the derivation and application of quantum algorithms for adiabatic and non-adiabatic quantum dynamics, which include efficient approaches for the calculation of the BO potential energy surfaces (PESs). Thereafter, we study the time-evolution of a model system consisting of two coupled PESs in first and second quantization. In a first application, we discuss a recently introduced quantum algorithm for the evolution of a wavepacket in first quantization and exploit the potential quantum advantage of mapping its spatial grid representation to logarithmically many qubits. For the second demonstration, we move to the second quantization framework and review the scaling properties of two alternative time-evolution algorithms, namely, a variational quantum algorithm (VQA) (based on the McLachlan variational principle) and conventional Trotter-type evolution (based on a Lie-Trotter-Suzuki formula). Both methods clearly demonstrate the potential of quantum algorithms and their favorable scaling compared to the available classical approaches. However, a clear demonstration of quantum advantage in the context of molecular quantum dynamics may require the implementation of these algorithms in fault-tolerant quantum computers, while their application in near-term, noisy quantum devices is still unclear and deserves further investigation.
Gaussian wavepacket (GWP) methods are prevalent means of solving the nuclear time-dependent Schrodinger equation (TDSE) on classical computers. They consist of representing the nuclear wave function as a superposition of Gaussians. Here, we present a variational hybrid continuous-variable discrete-variable (CV-DV) algorithm for simulating adiabatic nuclear dynamics on qubit-oscillator quantum hardware. It encodes a superposition of frozen Gaussians (FGs) onto the qubit-oscillator register and evolves it according to the time-dependent variational principle (VP). We test two variants of the approach, one that evolves a single FG and another that evolves a superposition of FGs, on a set of prototypical 1D harmonic and anharmonic systems. We simulate the photoexcited state dynamics of SO2 using the single FG variant of the algorithm and compute its autocorrelation function on physical quantum hardware. For harmonic and Morse potentials, the variant evolving the superposition of FGs converges to the numerically exact solution as the number of Gaussians increases. In the case of the double-well potential, the algorithm simulates wavepacket bifurcation and recurrence correctly, and the ensuing irregular dynamics moderately well. Overall, this work establishes a pathway for simulating molecular dynamics on noisy intermediate-scale quantum devices.
Dynamical simulation of quantum many-body systems is a central task in quantum chemistry and requires efficient wavefunction representations. Although tensor-network and neural-network quantum states have achieved considerable success in ground-state calculations, entanglement growth hinders their application to quantum dynamics. Quantum computing may offer a route to quantum advantage, but algorithms such as Trotterization generally require fault-tolerant quantum computers, whereas near-term hardware supports only limited circuit sizes. Variational quantum simulation (VQS) instead represents the evolving state with a single shallow, fixed-size parameterized quantum circuit (PQC). Here, we introduce nonorthogonal variational quantum simulation (NOVQS), which applies linear combinations of parameterized quantum states to real- and imaginary-time evolutions. We design a shallow, hardware-friendly ansatz tailored to second-quantized electronic-structure Hamiltonians, together with resource-efficient protocols for measuring the matrices and vectors in the parameter equations of motion. Error analysis and resource estimation are also provided. Numerical simulations of hydrogen chains and the nitrogen molecule demonstrate that a collection of shallow, or even single-layer, PQCs can match or outperform a much deeper PQC in VQS. In particular, we identify a trade-off between circuit number and depth. NOVQS therefore provides a flexible route to enhancing wavefunction expressivity under circuit-depth constraints, making it a promising framework for quantum dynamics simulations on near-term quantum processors.