Keith H. Hughes

A new method for wave packet dynamics: Derivative propagation along quantum trajectories

A new method is proposed for computing the time evolution of quantum mechanical wave packets. Equations of motion for the real-valued functions C and S in the complex action S=C(r,t)+iS(r,t)/ℏ, with ψ(r,t)=exp(S), involve gradients and curvatures of C and S. In previous implementations of the hydrodynamic formulation, various time-consuming fitting techniques of limited accuracy were used to evaluate these derivatives around each fluid element in an evolving ensemble. In this study, equations of motion are developed for the spatial derivatives themselves and a small set of these are integrated along quantum trajectories concurrently with the equations for C and S. Significantly, quantum effects can be included at various orders of approximation, no spatial fitting is involved, there are no basis set expansions, and single quantum trajectories (rather than correlated ensembles) may be propagated, one at a time. Excellent results are obtained when the derivative propagation method is applied to anharmonic p...

Effective-mode representation of non-Markovian dynamics: A hierarchical approximation of the spectral density. I. Application to single surface dynamics

An approach to non-Markovian system-environment dynamics is described which is based on the construction of a hierarchy of coupled effective environmental modes that is terminated by coupling the final member of the hierarchy to a Markovian bath. For an arbitrary environment, which is linearly coupled to the subsystem, the discretized spectral density is replaced by a series of approximate spectral densities involving an increasing number of effective modes. This series of approximants, which are constructed analytically in this paper, guarantees the accurate representation of the overall system-plus-bath dynamics up to increasing times. The hierarchical structure is manifested in the approximate spectral densities in the form of the imaginary part of a continued fraction similar to Mori theory. The results are described for cases where the hierarchy is truncated at the first-, second-, and third-order level. It is demonstrated that the results generated from a reduced density matrix equation of motion and large dimensional system-plus-bath wavepacket calculations are in excellent agreement. For the reduced density matrix calculations, the system and hierarchy of effective modes are treated explicitly and the effects of the bath on the final member of the hierarchy are described by the Caldeira-Leggett equation and its generalization to zero temperature.

Effective-mode representation of non-Markovian dynamics: A hierarchical approximation of the spectral density. II. Application to environment-induced nonadiabatic dynamics

The non-Markovian approach developed in the companion paper [Hughes et al., J. Chem. Phys. 131, 024109 (2009)], which employs a hierarchical series of approximate spectral densities, is extended to the treatment of nonadiabatic dynamics of coupled electronic states. We focus on a spin-boson-type Hamiltonian including a subset of system vibrational modes which are treated without any approximation, while a set of bath modes is transformed to a chain of effective modes and treated in a reduced-dimensional space. Only the first member of the chain is coupled to the electronic subsystem. The chain construction can be truncated at successive orders and is terminated by a Markovian closure acting on the end of the chain. From this Mori-type construction, a hierarchy of approximate spectral densities is obtained which approach the true bath spectral density with increasing accuracy. Applications are presented for the dynamics of a vibronic subsystem comprising a high-frequency mode and interacting with a low-frequency bath. The bath is shown to have a striking effect on the nonadiabatic dynamics, which can be rationalized in the effective-mode picture. A reduced two-dimensional subspace is constructed which accounts for the essential features of the nonadiabatic process induced by the effective environmental mode. Electronic coherence is found to be preserved on the shortest time scale determined by the effective mode, while decoherence sets in on a longer time scale. Numerical simulations are carried out using either an explicit wave function representation of the system and overall bath or else an explicit representation of the system and effective-mode part in conjunction with a Caldeira-Leggett master equation.

Wavepacket dynamics on dynamically adapting grids: application of the equidistribution principle

A moving grid approach to wavepacket dynamics is described that enables grid points to be used efficiently in regions where high resolution of the wavepacket is required. The grid movement is based on the principle of equidistribution and by using a grid smoothing technique the grid points trace a path that continuously adapt to reflect the dynamics of the wavepacket. The technique is robust and allows accurate computations to be obtained for long wavepacket propagation times. Results are presented for two systems: tunnelling dynamics in a double well potential and scattering of a wavepacket from a repulsive Eckart barrier.

Visualization of molecular quantum dynamics: a molecular visualization tool with integrated Web3D and haptics

The Department of Chemistry and the School of Informatics at the University of Wales, Bangor are working together to create tools for the visualization of molecular quantum dynamics. This paper presents the results of our initial work. A prototype Molecular Visualiser (MV) application has been developed based on Web3D standards, plus extensions for support of haptic interaction. MV provides the user with visualizations of molecular systems, potential energy surfaces, and wavepacket dynamics. These can be displayed in a web browser using VRML, or be delivered to a virtual environment in which haptic properties have been assigned based on the molecular dynamics of the system. The use of MV for both research and teaching is discussed.

Extended hydrodynamic approach to quantum-classical nonequilibrium evolution. I. Theory

A mixed quantum-classical formulation is developed for a quantum subsystem in strong interaction with an N-particle environment, to be treated as classical in the framework of a hydrodynamic representation. Starting from the quantum Liouville equation for the N-particle distribution and the corresponding reduced single-particle distribution, exact quantum hydrodynamic equations are obtained for the momentum moments of the single-particle distribution coupled to a discretized quantum subsystem. The quantum-classical limit is subsequently taken and the resulting hierarchy of equations is further approximated by various closure schemes. These include, in particular, (i) a Grad-Hermite-type closure, (ii) a Gaussian closure at the level of a quantum-classical local Maxwellian distribution, and (iii) a dynamical density functional theory approximation by which the hydrodynamic pressure term is replaced by a free energy functional derivative. The latter limit yields a mixed quantum-classical formulation which has previously been introduced by I. Burghardt and B. Bagchi, Chem. Phys. 134, 343 (2006).

Communication: Universal Markovian reduction of Brownian particle dynamics

Non-Markovian processes can often be turned Markovian by enlarging the set of variables. Here we show, by an explicit construction, how this can be done for the dynamics of a Brownian particle obeying the generalized Langevin equation. Given an arbitrary bath spectral density J(0), we introduce an orthogonal transformation of the bath variables into effective modes, leading stepwise to a semi-infinite chain with nearest-neighbor interactions. The transformation is uniquely determined by J(0) and defines a sequence {J(n)}(n∈N) of residual spectral densities describing the interaction of the terminal chain mode, at each step, with the remaining bath. We derive a simple one-term recurrence relation for this sequence and show that its limit is the quasi-Ohmic expression provided by the Rubin model of dissipation. Numerical calculations show that, irrespective of the details of J(0), convergence is fast enough to be useful in practice for an effective Ohmic reduction of the dissipative dynamics.

Wavepacket dynamics on arbitrary Lagrangian–Eulerian grids: Application to an Eckart barrier

An adaptive grid approach to a computational study of the scattering of a wavepacket from a repulsive Eckart barrier is described. The grids move in an arbitrary Lagrangian–Eulerian (ALE) framework and a hybrid of the moving path transform of the Schrodinger equation and the hydrodynamic equations are used for the equations of motion. Boundary grid points follow Lagrangian trajectories and interior grid points follow non-Lagrangian paths. For the hydrodynamic equations the interior grid points are equally spaced between the evolving Lagrangian boundaries. For the moving path transform of the Schrodinger equation interior grid distribution is determined by the principle of equidistribution, and by using a grid smoothing technique these grid points trace a path that continuously adapts to reflect the dynamics of the wavepacket. The moving grid technique is robust and allows accurate computations to be obtained with a small number of grid points for wavepacket propagation times exceeding 5 ps.

Quantum dynamics of hydrogen atoms on graphene. I. System-bath modeling

An accurate system-bath model to investigate the quantum dynamics of hydrogen atoms chemisorbed on graphene is presented. The system comprises a hydrogen atom and the carbon atom from graphene that forms the covalent bond, and it is described by a previously developed 4D potential energy surface based on density functional theory ab initio data. The bath describes the rest of the carbon lattice and is obtained from an empirical force field through inversion of a classical equilibrium correlation function describing the hydrogen motion. By construction, model building easily accommodates improvements coming from the use of higher level electronic structure theory for the system. Further, it is well suited to a determination of the system-environment coupling by means of ab initio molecular dynamics. This paper details the system-bath modeling and shows its application to the quantum dynamics of vibrational relaxation of a chemisorbed hydrogen atom, which is here investigated at T = 0 K with the help of the multi-configuration time-dependent Hartree method. Paper II deals with the sticking dynamics.

Non-Markovian reduced dynamics based upon a hierarchical effective-mode representation.

A reduced dynamics representation is introduced which is tailored to a hierarchical, Mori-chain type representation of a bath of harmonic oscillators which are linearly coupled to a subsystem. We consider a spin-boson system where a single effective mode is constructed so as to absorb all system-environment interactions, while the residual bath modes are coupled bilinearly to the primary mode and among each other. Using a cumulant expansion of the memory kernel, correlation functions for the primary mode are obtained, which can be suitably approximated by truncated chains representing the primary-residual mode interactions. A series of reduced-dimensional bath correlation functions is thus obtained, which can be expressed as Fourier-Laplace transforms of spectral densities that are given in truncated continued-fraction form. For a master equation which is second order in the system-bath coupling, the memory kernel is re-expressed in terms of local-in-time equations involving auxiliary densities and auxiliary operators.