Abstract: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.

Quasi-direct Quantum Molecular Dynamics: The Time-Dependent Adaptive Density-Guided Approach for Potential Energy Surface Construction

An overview of nonadiabatic dynamics simulations methods, with focus on the direct approach versus the fitting of potential energy surfaces

Calculating Potential Energy Surfaces with Quantum Computers by Measuring Only the Density Along Adiabatic Transitions

Nonadiabatic Dynamics Via the Symmetrical Quasi-Classical Method in the Presence of Anharmonicity.

Quantum Dynamics from Classical Trajectories

Molecular Quantum Dynamics: A Quantum Computing Perspective

QDπ: A Quantum Deep Potential Interaction Model for Drug Discovery

Quasiclassical Trajectory Simulation as a Protocol to Build Locally Accurate Machine Learning Potentials

Improved on-the-Fly MCTDH Simulations with Many-Body-Potential Tensor Decomposition and Projection Diabatization

Advanced Potential Energy Surfaces for Molecular Simulation

The symmetric quasi-classical model using on-the-fly time-dependent density functional theory within the Tamm-Dancoff approximation

Full Quantum Dynamics Simulation of a Realistic Molecular System Using the Adaptive Time-Dependent Density Matrix Renormalization Group Method

Ab initio potential energy surfaces calculation via restricted Hartree–Fock for molecular dynamics simulation: a comprehensive review

Quantum Car-Parrinello Molecular Dynamics: A Cost-Efficient Molecular Simulation Method on Near-Term Quantum Computers

Deep Potential Molecular Dynamics: A Scalable Model with the Accuracy of Quantum Mechanics

Quantum Quality with Classical Cost: Ab Initio Nonadiabatic Dynamics Simulations Using the Mapping Approach to Surface Hopping

Refining Potential Energy Surface through Dynamical Properties via Differentiable Molecular Simulation

Quantum Reaction Dynamics Based on a New Generation Density Functional and Neural Network Potential Energy Surfaces

Nonadiabatic Molecular Dynamics with Subsystem Density Functional Theory: Application to Crystalline Pentacene