All elements of a quantum algorithm for calculations of rotationally inelastic molecule + atom scattering within the framework of a mixed quantum/classical theory are outlined. In this approach, the rotational motion of the molecule is described quantum mechanically using the time-dependent Schrödinger equation, while the scattering process of two collision partners is treated classically. The matrix of potential coupling is precomputed on a classical processor, whereas the quantum hardware is used to propagate the system of coupled equations for the rotational state-to-state transitions. All quantum circuits needed for practical implementation of the algorithm are presented. First, the quantum codes written in Qiskit are rigorously tested by running calculations for a N 2 + O collision on a classical emulator of quantum hardware using a realistic potential energy surface of this system and comparing these results against the results obtained by the MQCT code. Next, these codes are run on the actual quantum hardware, such as the publicly available IBM Brisbane, Kyiv, and Sherbrooke . A very good agreement with benchmark data was obtained. To the best of our knowledge, this is the first proof-of-principle calculation of inelastic scattering implemented successfully on a quantum computer using a case study within mixed quantum/classical framework.
All elements of a quantum algorithm for calculations of rotationally inelastic molecule + atom scattering within the framework of a mixed quantum/classical theory are outlined. In this approach, the rotational motion of the molecule is described quantum mechanically using the time-dependent Schrödinger equation, while the scattering process of two collision partners is treated classically. The matrix of potential coupling is precomputed on a classical processor, whereas the quantum hardware is used to propagate the system of coupled equations for the rotational state-to-state transitions. All quantum circuits needed for practical implementation of the algorithm are presented. First, the quantum codes written in Qiskit are rigorously tested by running calculations for a N 2 + O collision on a classical emulator of quantum hardware using a realistic potential energy surface of this system and comparing these results against the results obtained by the MQCT code. Next, these codes are run on the actual quantum hardware, such as the publicly available IBM Brisbane, Kyiv, and Sherbrooke . A very good agreement with benchmark data was obtained. To the best of our knowledge, this is the first proof-of-principle calculation of inelastic scattering implemented successfully on a quantum computer using a case study within mixed quantum/classical framework.
We explore how the fundamental problems in quantum molecular dynamics can be modeled using classical simulators (emulators) of quantum computers and the actual quantum hardware available to us today. The list of problems we tackle includes propagation of a free wave packet, vibration of a harmonic oscillator, and tunneling through a barrier. Each of these problems starts with the initial wave packet setup. Although Qiskit provides a general method for initializing wave functions, in most cases it generates deep quantum circuits. While these circuits perform well on noiseless simulators, they suffer from excessive noise on quantum hardware. To overcome this issue, we designed a shallower quantum circuit for preparing a Gaussian-like initial wave packet, which improves the performance of real hardware. Next, quantum circuits are implemented to apply the kinetic and potential energy operators for the evolution of a wave function over time. The results of our modeling on classical emulators of quantum hardware agree perfectly with the results obtained using the traditional (classical) methods. This serves as a benchmark and demonstrates that the quantum algorithms and Qiskit codes we developed are accurate. However, the results obtained on the actual quantum hardware available today, such as IBM's superconducting qubits and IonQ's trapped ions, indicate large discrepancies due to hardware limitations. This work highlights both the potential and challenges of using quantum computers to solve fundamental quantum molecular dynamics problems.
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.
Collisions between atoms and molecules govern how energy is exchanged and redistributed in the gas phase, underpinning phenomena from spectral line broadening and interstellar molecular excitation to the cooling and trapping of molecules. A quantitative account of these collisions requires connecting the theoretical description of the interaction between the collision partners to measurable quantities such as cross sections. However, scattering theory as presented in quantum mechanics textbooks does not readily translate into numerical methods. Here, we formulate the scattering problem in Dirac notation, as customary in present-day quantum mechanics courses, and show how to build a practical numerical framework from it. We derive the "coupled channel" representation of the scattering problem with correct asymptotic behaviour and introduce Numerov's method to obtain time-independent wavefunctions, from which we calculate cross sections, the observables measured in scattering experiments. Starting from single channel scattering, we build up to the realistic collision of a ground-state atom with a diatomic molecule.
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 coupled Schrödinger equations that describe nonadiabatic dynamics are recast using the Bohm formulation of quantum mechanics. The resulting coupled Bohm equations are solved numerically for two scattering models, giving results that are essentially identical to wave-packet solution of the original coupled Schrödinger equations. The classical limit of the set of coupled Bohm equations is then described, producing a mixed quantum-classical theory incorporating classical-like motion on each potential energy surface accompanied by quantum transitions between the quantum states. Numerical tests of the mixed quantum-classical method are in excellent agreement with the accurate full-quantum results for the model problems. The method is contrasted with the related surface-hopping approach. It is shown that computing the dynamics of a distribution of classical particles is more consistent and more accurate than computing the motion of independent point particles as with surface hopping.
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.