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 J0, 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 J0 and defines a sequence {Jn} 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 J0, convergence is fast enough to be useful in practice for an effective Markovian reduction of quantum dissipative 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.

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.

Memory effects in quantum dynamical processes involving structured environments are presently difficult, if not impossible, to investigate using standard approaches. Progress can be made by transforming the environmental variables to a suitable chain representation which effectively performs a Markovian embedding of the dynamics. Here, we show that this effective-mode chain representation provides a unique way of unraveling the memory kernel κ(t) as a function of time. Truncated or Markov-closed chains with n effective modes exactly reproduce κ(t) to the 4nth order in time, up to an irrelevant constant of order κ(0)/n. These favorable convergence properties pave the way for efficient quantum simulations of fast (non-Markovian) processes by reduced dynamical models.

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 potentials involving barrier transmission.