A time-reversible integrator for the time-dependent Schrödinger equation on an adaptive grid.

One of the most accurate methods for solving the time-dependent Schrödinger equation uses a combination of the dynamic Fourier method with the split-operator algorithm on a tensor-product grid. To reduce the number of required grid points, we let the grid move together with the wavepacket but find that the naïve algorithm based on an alternate evolution of the wavefunction and grid destroys the time reversibility of the exact evolution. Yet, we show that the time reversibility is recovered if the wavefunction and grid are evolved simultaneously during each kinetic or potential step; this is achieved by using the Ehrenfest theorem together with the splitting method. The proposed algorithm is conditionally stable, symmetric, and time-reversible and conserves the norm of the wavefunction. The preservation of these geometric properties is shown analytically and demonstrated numerically on a three-dimensional harmonic model and collinear model of He-H2 scattering. We also show that the proposed algorithm can be symmetrically composed to obtain time-reversible integrators of an arbitrary even order. We observed 10 000-fold speedup by using the tenth-order instead of the second-order method to obtain a solution with a time discretization error below 10-9. Moreover, using the adaptive grid instead of the fixed grid resulted in a 64-fold reduction in the required number of grid points in the harmonic system and made it possible to simulate the He-H2 scattering for six times longer while maintaining reasonable accuracy. Applicability of the algorithm to high-dimensional quantum dynamics is demonstrated using the strongly anharmonic eight-dimensional Hénon-Heiles model.

[1]  P. Ehrenfest.  Bemerkung über die angenäherte Gültigkeit der klassischen Mechanik innerhalb der Quantenmechanik , 1927 .

[2]  J. Crank,et al.  A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction type , 1947 .

[3]  C. Lanczos.  An iteration method for the solution of the eigenvalue problem of linear differential and integral operators , 1950 .

[4]  H. Trotter.  On the product of semi-groups of operators , 1959 .

[5]  D. Secrest,et al.  Exact Quantum‐Mechanical Calculation of a Collinear Collision of a Particle with a Harmonic Oscillator , 1966 .

[6]  L. Verlet.  Computer "Experiments" on Classical Fluids. I. Thermodynamical Properties of Lennard-Jones Molecules , 1967 .

[7]  G. Strang.  On the Construction and Comparison of Difference Schemes , 1968 .

[8]  R. Wyatt,et al.  Dynamics of the Collinear H+H2 Reaction. I. Probability Density and Flux , 1971 .

[9]  E. Heller.  Time‐dependent approach to semiclassical dynamics , 1975 .

[10]  Michael J Davis,et al.  Semiclassical Gaussian basis set method for molecular vibrational wave functions , 1979 .

[11]  M. Feit,et al.  Solution of the Schrödinger equation by a spectral method , 1982 .

[12]  G. A. Parker,et al.  Discrete variable representations and sudden models in quantum scattering theory , 1982 .

[13]  R. Kosloff,et al.  A fourier method solution for the time dependent Schrödinger equation as a tool in molecular dynamics , 1983 .

[14]  H. Tal-Ezer,et al.  An accurate and efficient scheme for propagating the time dependent Schrödinger equation , 1984 .

[15]  T. Park,et al.  Unitary quantum time evolution by iterative Lanczos reduction , 1986 .

[16]  A. Zewail,et al.  Molecular state evolution after excitation with an ultra-short laser pulse: A quantum analysis of NaI and NaBr dissociation , 1988 .

[17]  R. Kosloff.  Time-dependent quantum-mechanical methods for molecular dynamics , 1988 .

[18]  M. Suzuki,et al.  Fractal decomposition of exponential operators with applications to many-body theories and Monte Carlo simulations , 1990 .

[19]  H. Yoshida.  Construction of higher order symplectic integrators , 1990 .

[20]  U. Manthe,et al.  The multi-configurational time-dependent Hartree approach , 1990 .

[21]  R. Coalson,et al.  Adding configuration interaction to the time‐dependent Hartree grid approximation , 1990 .

[22]  U. Manthe,et al.  Wave‐packet dynamics within the multiconfiguration Hartree framework: General aspects and application to NOCl , 1992 .

[23]  T. Carrington,et al.  A general discrete variable method to calculate vibrational energy levels of three‐ and four‐atom molecules , 1993 .

[24]  N. Makri,et al.  Time-dependent discrete variable representations for quantum wave packet propagation , 1995 .

[25]  G. Stock,et al.  Resonance Raman spectroscopy of the S1 and S2 states of pyrazine: Experiment and first principles calculation of spectra , 1995 .

[26]  William Kahan,et al.  Composition constants for raising the orders of unconventional schemes for ordinary differential equations , 1997, Math. Comput..

[27]  R. Wyatt,et al.  Quantum Wave Packet Dynamics with Trajectories , 1999 .

[28]  R. Wyatt.  Quantum wave packet dynamics with trajectories: Application to reactive scattering , 1999 .

[29]  R. Wyatt,et al.  Quantum wave packet dynamics with trajectories: Implementation with adaptive Lagrangian grids , 2000 .

[30]  S. Adhikari,et al.  A time-dependent discrete variable representation method , 2000 .

[31]  T. Martínez,et al.  Ab Initio Multiple Spawning: Photochemistry from First Principles Quantum Molecular Dynamics , 2000 .

[32]  A. Bandrauk,et al.  Moving adaptive grid methods for numerical solution of the time-dependent molecular Schrödinger equation in laser fields , 2001 .

[33]  Keith H. Hughes,et al.  Wavepacket dynamics on dynamically adapting grids: application of the equidistribution principle , 2002 .

[34]  H. Meyer,et al.  Benchmark calculations on high-dimensional Henon–Heiles potentials with the multi-configuration time dependent Hartree (MCTDH) method , 2002 .

[35]  R. Wyatt.  Wave packet dynamics on adaptive moving grids , 2002 .

[36]  Haobin Wang,et al.  Multilayer formulation of the multiconfiguration time-dependent Hartree theory , 2003 .

[37]  Mark Sofroniou,et al.  Derivation of symmetric composition constants for symmetric integrators , 2005, Optim. Methods Softw..

[38]  S. Adhikari,et al.  Quantum dynamics of inelastic scattering with a moving grid , 2005 .

[39]  Steven G. Johnson,et al.  The Design and Implementation of FFTW3 , 2005, Proceedings of the IEEE.

[40]  Wave packet dynamics with adaptive grids: The moving boundary truncation method , 2006 .

[41]  T. Carrington,et al.  Discrete‐Variable Representations and their Utilization , 2007 .

[42]  Vasile Gradinaru,et al.  Strang Splitting for the Time-Dependent Schrödinger Equation on Sparse Grids , 2007, SIAM J. Numer. Anal..

[43]  Lorenz S. Cederbaum,et al.  Using the MCTDH wavepacket propagation method to describe multimode non-adiabatic dynamics , 2008 .

[44]  C. Lubich.  From Quantum to Classical Molecular Dynamics: Reduced Models and Numerical Analysis , 2008 .

[45]  A. Pounds.  Introduction to Quantum Mechanics: A Time-Dependent Perspective (David J. Tannor) , 2008 .

[46]  T. Carrington,et al.  A discrete variable representation method for studying the rovibrational quantum dynamics of molecules with more than three atoms. , 2009, The Journal of chemical physics.

[47]  T. Carrington,et al.  Nonproduct quadrature grids for solving the vibrational Schrödinger equation. , 2009, The Journal of chemical physics.

[48]  D. Tannor,et al.  Phase-space approach to solving the time-independent Schrödinger equation. , 2012, Physical review letters.

[49]  B. Poirier,et al.  Accurate quantum dynamics calculations using symmetrized Gaussians on a doubly dense Von Neumann lattice. , 2012, The Journal of chemical physics.

[50]  S. Adhikari,et al.  NEARLY LINEAR SCALABILITY OF TIME-DEPENDENT DISCRETE VARIABLE REPRESENTATION (TDDVR) METHOD FOR THE DYNAMICS OF MULTI-SURFACE MULTI-MODE HAMILTONIAN , 2013 .

[51]  S. Adhikari,et al.  Multisurface multimode molecular dynamical simulation of naphthalene and anthracene radical cations by using nearly linear scalable time-dependent discrete variable representation method. , 2014, The journal of physical chemistry. A.

[52]  Andr'e Nauts,et al.  Quantum dynamics with sparse grids: a combination of Smolyak scheme and cubature. Application to methanol in full dimensionality. , 2014, Spectrochimica acta. Part A, Molecular and biomolecular spectroscopy.

[53]  J. Vaníček,et al.  On-the-fly ab initio semiclassical dynamics: identifying degrees of freedom essential for emission spectra of oligothiophenes. , 2014, The Journal of chemical physics.

[54]  Tucker Carrington,et al.  Using a pruned, nondirect product basis in conjunction with the multi-configuration time-dependent Hartree (MCTDH) method. , 2016, The Journal of chemical physics.

[55]  Oliver T. Unke,et al.  Nonadiabatic effects in electronic and nuclear dynamics , 2017, Structural dynamics.

[56]  Tomislav Beguvsi'c,et al.  On-the-fly ab initio three thawed Gaussians approximation: A semiclassical approach to Herzberg-Teller spectra , 2018, Chemical Physics.

[57]  H. Kono,et al.  Automatic spatial extension of a time-dependent wavefunction expanded in terms of Gaussians: Application to multidimensional tunneling , 2018, Chemical Physics Letters.

[58]  T. Carrington,et al.  A new collocation-based multi-configuration time-dependent Hartree (MCTDH) approach for solving the Schrödinger equation with a general potential energy surface. , 2018, The Journal of chemical physics.

[59]  Chia-Chun Chou,et al.  Moving Boundary Truncated Grid Method for Wave Packet Dynamics. , 2018, The journal of physical chemistry. A.

[60]  Souvik Mandal,et al.  The TDDVR approach for molecular photoexcitation, molecule–surface and triatomic reactive scattering processes , 2018, International Reviews in Physical Chemistry.

[61]  J. Vaníček,et al.  Efficient geometric integrators for nonadiabatic quantum dynamics. I. The adiabatic representation. , 2019, The Journal of chemical physics.

[62]  J. Vaníček,et al.  Single-Hessian thawed Gaussian approximation. , 2019, The Journal of chemical physics.

[63]  T. Carrington,et al.  A pruned collocation-based multiconfiguration time-dependent Hartree approach using a Smolyak grid for solving the Schrödinger equation with a general potential energy surface. , 2019, The Journal of chemical physics.

[64]  J. Vaníček,et al.  Efficient geometric integrators for nonadiabatic quantum dynamics. II. The diabatic representation. , 2019, The Journal of chemical physics.