Basis set generation for quantum dynamics simulations using simple trajectory-based methods

Methods for solving the time-dependent Schrödinger equation generally employ either a global static basis set, which is fixed at the outset, or a dynamic basis set, which evolves according to classical-like or variational equations of motion; the former approach results in the well-known exponential scaling with system size, while the latter can suffer from challenging numerical problems, such as singular matrices, as well as violation of energy conservation. Here, we suggest a middle road: building a basis set using trajectories to place time-independent basis functions in the regions of phase space relevant to wave function propagation. This simple approach, which potentially circumvents many of the problems traditionally associated with global or dynamic basis sets, is successfully demonstrated for two challenging benchmark problems in quantum dynamics, namely, relaxation dynamics following photoexcitation in pyrazine, and the spin Boson model

A quantum-classical treatment of non-adiabatic transitions

We have carried out molecular dynamics simulations of non-adiabatic processes with the help of a newly formulated potentially exact quantum-classical approach derived from a method proposed earlier [J. Chem. Phys. 118 (2003) 5302]. In this method, time-dependent Schroedinger equation is solved by representing Psi on a moving Gauss-Hermite DVR grid, the motion of grid-centre being handled classically, but self consistently with the quantum evolution of the wavefunction. Electronic transitions are allowed anywhere in the configuration space among any number of coupled states. We have tested the method on three model problems proposed by J.C. Tully [J. Chem. Phys. 93 (1990) 1061]. These models are relevant to a wide range of gas-phase and condensed-phase phenomena occurring even at low energies. Excellent agreement of computed transition probabilities with corresponding quantum mechanical (DVR/FFT) results even in the deep quantum regime and its cost-efficiency (computational) are encouraging. (C) 200

A quantum-classical treatment of non-adiabatic transitions

We have carried out molecular dynamics simulations of non-adiabatic processes with the help of a newly formulated potentially exact quantum-classical approach derived from a method proposed earlier [J. Chem. Phys. 118 (2003) 5302]. In this method, time-dependent Schroedinger equation is solved by representing Psi on a moving Gauss-Hermite DVR grid, the motion of grid-centre being handled classically, but self consistently with the quantum evolution of the wavefunction. Electronic transitions are allowed anywhere in the configuration space among any number of coupled states. We have tested the method on three model problems proposed by J.C. Tully [J. Chem. Phys. 93 (1990) 1061]. These models are relevant to a wide range of gas-phase and condensed-phase phenomena occurring even at low energies. Excellent agreement of computed transition probabilities with corresponding quantum mechanical (DVR/FFT) results even in the deep quantum regime and its cost-efficiency (computational) are encouraging. (C) 200