The simulation of ion-atom collisions remains a formidable challenge due to the complex interplay between electronic and nuclear degrees of freedom. We present a hybrid quantum-classical computing framework for simulating time-dependent ion-atom collision dynamics, within which two variational quantum time evolution algorithms are implemented. To validate our framework, we simulate the charge-transfer dynamics and compute the corresponding cross sections for the proton-hydrogen collision system across an energy range 1--25 keV. Our results accurately reproduce the charge-transfer dynamics with high fidelity and exhibit very good agreement with available experimental and theoretical cross-section data across the entire energy range. These results highlight the accuracy and applicability of our hybrid quantum-classical framework for scattering cross-section calculations. Our work demonstrates an effective approach for mapping time-dependent many-body collision problems onto near-term quantum computing devices, and also provides promising directions for practical applications of universal quantum computing in the noisy intermediate-scale quantum era.
The simulation of ion-atom collisions remains a formidable challenge due to the complex interplay between electronic and nuclear degrees of freedom. We present a hybrid quantum-classical computing framework for simulating time-dependent ion-atom collision dynamics, within which two variational quantum time evolution algorithms are implemented. To validate our framework, we simulate the charge-transfer dynamics and compute the corresponding cross sections for the proton-hydrogen collision system across an energy range 1--25 keV. Our results accurately reproduce the charge-transfer dynamics with high fidelity and exhibit very good agreement with available experimental and theoretical cross-section data across the entire energy range. These results highlight the accuracy and applicability of our hybrid quantum-classical framework for scattering cross-section calculations. Our work demonstrates an effective approach for mapping time-dependent many-body collision problems onto near-term quantum computing devices, and also provides promising directions for practical applications of universal quantum computing in the noisy intermediate-scale quantum era.
Nuclear quantum effects are critical for describing proton transfer and hydrogen bonding, but their incorporation into quantum chemistry calculations is often computationally prohibitive on classical hardware. A promising alternative are quantum computers due to their linear scaling in the space requirements with system size. We introduce a novel hybrid quantum-classical algorithm for ab-initio molecular dynamics that incorporates nuclear quantum effects via the nuclear-electronic orbital method. The proposed algorithm evaluates ground state energies and forces on the quantum computer using a variational quantum eigensolver, while the molecular geometries are updated classically. We validate our approach through simulations of H2, H2O and the Zundel ion H5O2+, comparing the simulated vibrational spectra with experimental data. Upon inclusion of nuclear quantum effects, the proton shuttling movement in the Zundel ion becomes effectively barrierless, and errors in the simulated frequencies improve significantly. Employing compact hardware-efficient ansatz circuits we achieve results comparable to the more accurate UCCSD ansatzes, which hints towards the feasibility of executing our algorithm on near-term quantum devices.
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.
The most widely used approach for simulating the dynamics of time-dependent Hamiltonians via quantum computation depends on the quantum-classical hybrid variational quantum time evolution algorithm, in which ordinary differential equations of the variational coefficients for determining time evolution are solved via classical simulations with a time discretization method. We here present a full-quantum approach, in which ordinary differential equations of the variational coefficients are transformed into static linear equations via the Chebyshev spectral discretization method and then solved via the quantum singular value transformation algorithm. Our full quantum algorithm avoids classical feedback, achieves exponential convergence for smooth Hamiltonians, and yields a quantum circuit depth that is independent of the number of time steps. We demonstrate two implementation strategies, with a global formulation designed for fault-tolerant architectures and a sequential formulation tailored to near-term devices, and validate the approach through numerical simulations of proton-hydrogen charge-transfer dynamics, a prototypical time-dependent quantum chemistry problem. This work establishes a systematic pathway from quantum-classical hybrid variational quantum algorithms to full-quantum solvers for general time-dependent Hamiltonians, particularly those whose dynamics admit compact variational descriptions, opening a route toward full quantum computational advantages in time-dependent simulations.
We propose a hybrid quantum-classical algorithm for solving the time-independent Schrödinger equation for atomic and molecular collisions. The algorithm is based on the S -matrix version of the Kohn variational principle, which computes the fundamental scattering S -matrix by inverting the Hamiltonian matrix expressed in the basis of square-integrable functions. The computational bottleneck of the classical algorithm─symmetric matrix inversion─is addressed here using the variational quantum linear solver (VQLS), a recently developed noisy intermediate-scale quantum (NISQ) algorithm for solving systems of linear equations. We apply our algorithm to single- and multichannel quantum scattering problems, obtaining accurate vibrational relaxation probabilities in collinear atom–molecule collisions. We also show how the algorithm could be scaled up to simulate collisions of large polyatomic molecules. Our results demonstrate that it is possible to calculate scattering cross sections and rates for complex molecular collisions on NISQ quantum processors, opening up the possibility of scalable digital quantum computation of gas-phase bimolecular collisions and reactions of relevance to astrochemistry and ultracold chemistry.
Accurate simulation of dynamic processes in molecules and reactions is among the most challenging problems in quantum chemistry. Quantum computers promise efficient chemical simulation, but the existing quantum algorithms require many logical qubits and gates, placing practical applications beyond existing technology. Here, we carry out the first quantum simulations of chemical dynamics by employing a more hardware-efficient encoding scheme that uses both qubits and bosonic degrees of freedom. Our trapped-ion device accurately simulates the dynamics of nonadiabatic chemical processes, which are among the most difficult problems in computational chemistry because they involve strong coupling between electronic and nuclear motions. We demonstrate the programmability and versatility of our approach by simulating the dynamics of three different molecules, as well as open-system dynamics in the condensed phase, all with the same quantum resources. Our approach requires orders of magnitude fewer resources than equivalent qubit-only quantum simulations, demonstrating the potential of using hybrid encoding schemes to accelerate quantum simulations of complex chemical processes, which could have applications in fields ranging from energy conversion and storage to biology and drug design.
Coupled quantum electron-nuclear dynamics is often associated with the Born-Huang expansion of the molecular wave function and the appearance of nonadiabatic effects as a perturbation. On the other hand, native multicomponent representations of electrons and nuclei also exist, which do not rely on any a priori approximation. However, their implementation is hampered by prohibitive scaling. Consequently, quantum computers offer a unique opportunity for extending their use to larger systems. Here, we propose a quantum algorithm for simulating the time-evolution of molecular systems and apply it to proton transfer dynamics in malonaldehyde, described as a rigid scaffold. The proposed quantum algorithm can be easily generalized to include the explicit dynamics of the classically described molecular scaffold. We show how entanglement between electronic and nuclear degrees of freedom can persist over long times if electrons do not follow the nuclear displacement adiabatically. The proposed quantum algorithm may become a valid candidate for the study of such phenomena when sufficiently powerful quantum computers become available.