Time dependent linear and nonlinear equations govern the evolution of an extensive set of physical systems and processes describing, to enumerate just a few, Bose-Einstein condensates; soliton propagation in optical and photonic band-gap fibers; quantum control of atomic and molecular collisions and reactions; highly-compressed liquids; and dense and ultracold plasmas. While the media vary substantially, the basic computational procedures have many common features. We focus on the nonlinear Schrodinger equation and discuss two powerful approaches to its propagation: the Arnoldi/Lanczos(AL)l and Real Space Product Formula(RSPF)2. Both provide efficient systematic approximations to the short-time exponential propagator that moves the solution between time steps. We implement the former in a discrete variable representation (DVR)3 both in spatial grid and finite element forms and the latter in a spatial mesh with a finite difference representation of the kinetic energy operator. Both approaches require O(N) operations to propagate the wavefunction between time steps and handle multidimensional systems. We shall also draw connections with Liouville formulations used in quantum molecular dynamics simulations of large collections of atoms and molecules. After briefly outlining these formulations, we shall discuss some of the varied applications.
Time dependent linear and nonlinear equations govern the evolution of an extensive set of physical systems and processes describing, to enumerate just a few, Bose-Einstein condensates; soliton propagation in optical and photonic band-gap fibers; quantum control of atomic and molecular collisions and reactions; highly-compressed liquids; and dense and ultracold plasmas. While the media vary substantially, the basic computational procedures have many common features. We focus on the nonlinear Schrodinger equation and discuss two powerful approaches to its propagation: the Arnoldi/Lanczos(AL)l and Real Space Product Formula(RSPF)2. Both provide efficient systematic approximations to the short-time exponential propagator that moves the solution between time steps. We implement the former in a discrete variable representation (DVR)3 both in spatial grid and finite element forms and the latter in a spatial mesh with a finite difference representation of the kinetic energy operator. Both approaches require O(N) operations to propagate the wavefunction between time steps and handle multidimensional systems. We shall also draw connections with Liouville formulations used in quantum molecular dynamics simulations of large collections of atoms and molecules. After briefly outlining these formulations, we shall discuss some of the varied applications.
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 time-dependent Schrödinger Equation (TDSE) is a parabolic partial differential equation (PDE) comparable to a diffusion equation but with imaginary time. Due to its first order time derivative, exponential integrators or propagators are natural methods to describe evolution in time of the TDSE, both for time-independent and time-dependent potentials. Two splitting methods based on Fer and/or Magnus expansions allow for developing unitary factorizations of exponentials with different accuracies in the time step △t. The unitary factorization of exponentials to high order accuracy depends on commutators of kinetic energy operators with potentials. Fourth-order accuracy propagators can involve negative or complex time steps, or real time steps only but with gradients of potentials, i.e. forces. Extending the propagators of TDSE's to imaginary time allows to also apply these methods to classical many-body dynamics, and quantum statistical mechanics of molecular systems.
Several theoretical methods for the computation of quantum dynamical quantities are formulated, implemented, and applied, in the overall context of overcoming the exponential scaling property of quantum mechanics. The convergence of the approximate forward-backward semiclassical dynamics (FBSD) method has been accelerated by working with only the imaginary part of the real-time position and momentum time correlation functions. Simple manipulations allow the computationally favorable limit of a single factorization (“bead”) of the Boltzmann operator to yield greatly improved results. Numerical examples for a one-dimensional anharmonic oscillator and Lennard-Jones fluids are presented. Bohm's trajectory formulation of quantum mechanics is an exact mapping of the time-dependent Schrödinger equation onto classical fluid dynamics, plus one extra term responsible for all non-classical effects, the dreaded quantum potential. A computational procedure has been introduced to stably solve Bohm's hydrodynamic equations of motion in the forward-backward quantum dynamics (FBQD) formulation of time correlation functions, and applied to one-dimensional bound systems at zero and low temperature. A key step in the formulation of this procedure was the use of Hamilton's law of varying action, a generalization of Hamilton's principle of stationary action applicable to initial-value problems. This has resulted in the first reported synthetic computation of non-stationary quantum trajectories for a fully anharmonic oscillator. Success in one dimension has laid the groundwork for future extension to a rigorous quantum-classical scheme, which is briefly discussed. Other unpublished methods for obtaining time correlation functions and propagators in quantum dynamics are also discussed.
We consider the numerical propagation of models that combine both quantum and classical degrees of freedom, usually, electrons and nuclei, respectively. We focus, in our computational examples, on the case in which the quantum electrons are modeled with time-dependent density-functional theory, although the methods discussed below can be used with any other level of theory. Often, for these so-called quantum-classical molecular dynamics models, one uses some propagation technique to deal with the quantum part and a different one for the classical equations. While the resulting procedure may, in principle, be consistent, it can however spoil some of the properties of the methods, such as the accuracy order with respect to the time step or the preservation of the geometrical structure of the equations. Few methods have been developed specifically for hybrid quantum-classical models. We propose using the same method for both the quantum and classical particles, in particular, one family of techniques that proves to be very efficient for the propagation of Schrödinger-like equations: the (quasi)-commutator free Magnus expansions. These have been developed, however, for linear systems, yet our problem is nonlinear: formally, the full quantum-classical system can be rewritten as a nonlinear Schrödinger equation, i.e., one in which the Hamiltonian depends on the system itself. The Magnus expansion algorithms for linear systems require the application of the Hamiltonian at intermediate points in a given propagating interval. For nonlinear systems, this poses a problem as this Hamiltonian is unknown due to its dependence on the state. We approximate it by employing a higher order extrapolation using previous steps as input. The resulting technique can then be regarded as a multistep technique or, alternatively, as a predictor corrector formula.
The time-dependent Schrödinger equation (TDSE) in real space is fundamental to understanding the dynamics of many-electron quantum systems, with applications ranging from quantum chemistry to condensed matter physics and materials science. However, solving the TDSE for complex fermionic systems remains a significant challenge, particularly due to the need to capture the time-evolving many-body correlations, while the antisymmetric nature of fermionic wavefunctions complicates the function space in which these solutions must be represented. We propose a general-purpose neural network framework for solving the real-space TDSE, Fermionic Antisymmetric Spatio-Temporal Network, which treats time as an explicit input alongside spatial coordinates, enabling a unified spatiotemporal representation of complex, antisymmetric wavefunctions for fermionic systems. This approach formulates the TDSE as a global optimization problem, avoiding step-by-step propagation and supporting highly parallelizable training. The method is demonstrated on five benchmark problems, achieving excellent agreement with reference solutions across all cases. These results demonstrate the method's accuracy and flexibility within the bound-state manifold across various dimensions and interaction regimes. While the current localized Ansatz inherently restricts the description of extensive ionization and continuum states, the method demonstrates the capability to stably simulate coherent multi-electron dynamics over extended time windows. Our framework offers a highly expressive alternative to traditional basis-dependent or mean-field methods, opening new possibilities for ab initio simulations of time-dependent quantum systems, with applications in quantum dynamics, molecular control, and ultrafast spectroscopy.
Various numerical integration schemes to calculate the propagation of a state following the time-dependent Schrodinger equation in the one dimensional position representation are presented and compared to each other. Three potentials have been used: a harmonic, a double-well and a zero potential. Eigenstates and a coherent state have been chosen as initial states. Special attention has been given to the long-time stability of the algorithms. These are: kinetic referenced split operator (KRSO), kinetic referenced Cayley (KRC), distributed approximating functions (DAF), Chebysheff expansion (CH), residuum minimization (RES), second-order differencing (SOD), an eigenstate expansion (EE) and a corrected kinetic referenced split operator (CKRSO). In addition, a speedup of the KRC and KRSO methods is presented which is specially suited when very few grid points are used. Numerical results are compared to analytically calculated values. Mixed classical/quantum mechanical simulations require a representation of the quantum state on a limited number of grid points, classical integration time steps of about one femtosecond and compatibility with methods to solve the time-ordering problem. For the considered potentials which differ quite essentially from the potentials used for scattering problems in particle physics, the EE method has been found to be faster, more accurate and more stable than the other methods if only a few grid points are required. Otherwise, good results have been obtained with KRC, KRSO, CH, DAF and RES. SOD has been found to be too slow, and CKRSO is not stable enough for long simulation times.