The quantum dynamics of a low-dimensional system in contact with a large but finite harmonic bath is theoretically investigated by coarse-graining the bath into a reduced set of effective energy states. In this model, the couplings between the system and the bath are obtained from statistically averaging over the discrete, degenerate effective states. Our model is aimed at intermediate bath sizes in which non-Markovian processes and energy transfer between the bath and the main system are important. The method is applied to a model system of a Morse oscillator coupled to 40 harmonic modes. The results are found to be in excellent agreement with the direct quantum dynamics simulations presented in the work of Bouakline et al. [J. Phys. Chem. A 116, 11118-11127 (2012)], but at a much lower computational cost. Extension to larger baths is discussed in comparison to the time-convolutionless method. We also extend this study to the case of a microcanonical bath with finite initial internal energies. The computational efficiency and convergence properties of the effective bath states model with respect to relevant parameters are also discussed.
The quantum dynamics of a low-dimensional system in contact with a large but finite harmonic bath is theoretically investigated by coarse-graining the bath into a reduced set of effective energy states. In this model, the couplings between the system and the bath are obtained from statistically averaging over the discrete, degenerate effective states. Our model is aimed at intermediate bath sizes in which non-Markovian processes and energy transfer between the bath and the main system are important. The method is applied to a model system of a Morse oscillator coupled to 40 harmonic modes. The results are found to be in excellent agreement with the direct quantum dynamics simulations presented in the work of Bouakline et al. [J. Phys. Chem. A 116, 11118-11127 (2012)], but at a much lower computational cost. Extension to larger baths is discussed in comparison to the time-convolutionless method. We also extend this study to the case of a microcanonical bath with finite initial internal energies. The computational efficiency and convergence properties of the effective bath states model with respect to relevant parameters are also discussed.
The dissipative quantum dynamics of an anharmonic oscillator coupled to a bath is studied with the purpose of elucidating the differences between the relaxation to a spin bath and to a harmonic bath. Converged results are obtained for the spin bath by the surrogate Hamiltonian approach. This method is based on constructing a system-bath Hamiltonian, with a finite but large number of spin bath modes, that mimics exactly a bath with an infinite number of modes for a finite time interval. Convergence with respect to the number of simultaneous excitations of bath modes can be checked. The results are compared to calculations that include a finite number of harmonic modes carried out by using the multiconfiguration time-dependent Hartree method of Nest and Meyer [J. Chem. Phys. 119, 24 (2003)]. In the weak coupling regime, at zero temperature and for small excitations of the primary system, both methods converge to the Markovian limit. When initially the primary system is significantly excited, the spin bath can saturate restricting the energy acceptance. An interaction term between bath modes that spreads the excitation eliminates the saturation. The loss of phase between two cat states has been analyzed and the results for the spin and harmonic baths are almost identical. For stronger couplings, the dynamics induced by the two types of baths deviate. The accumulation and degree of entanglement between the bath modes have been characterized. Only in the spin bath the dynamics generate entanglement between the bath modes.
We review an iterative path integral method that allows the computation of exact, long-time dynamics of small systems interacting with a dissipative bath. The method takes advantage of the finite memory of large environments at finite temperature to allow an iterative evaluation of the dynamics, thereby replacing an exponential scaling in simulation time with a linear scaling. This method is applied to calculate the dynamics of two model systems. The first consists of two-qubits interacting with a common bath. In this case we observe a variety of entanglement effects. If the qubits are initially separable, through an indirect coupling, the bath can create steady state entanglement between the qubits. This effect is due to the existence of a decoherence-free subspace as a result of the form of the qubit Hamiltonian and system-bath coupling. Entanglement created by the bath is shown to decrease with increasing temperature and system-bath coupling strength. However large system-bath coupling causes a faster increase in the entanglement. Initially entangled qubits lose their entanglement as a result of interactions with the bath, an effect that is heightened at higher temperatures. Direct coupling between the qubits is shown to slow the decay of entanglement and preserve some entanglement at long times; however at high temperatures this steady state entanglement becomes negligible. The second system we consider is vibrational degree of freedom coupled to a bath of harmonic oscillators or two-level systems. We compute four-time correlation functions which are used to calculate response functions relevant to third order infrared or seventh order Raman experiments for harmonic, Morse, and quadratic-quartic potentials. Our calculations reveal the role of potential features (anharmonicity and eigenvalue spectrum), both on short and long time scales, on the response function. Further, thermal excitation causes dramatic changes in the appearance of the response function, introducing symmetry with respect to the main diagonal. Finally, coupling to harmonic dissipative baths leads to decay of the response function (primarily along the t3 direction) and a broadening of the peaks in its Fourier transform. At high temperatures two-level-system baths are less efficient in destroying coherence than harmonic baths of similar parameters.
We study the dynamics of relaxation and thermalization in an exactly solvable model of a particle interacting with a harmonic oscillator bath. Our goal is to understand the effects of non-Markovian processes on the relaxational dynamics and to compare the exact evolution of the distribution function with approximate Markovian and non-Markovian quantum kinetics. There are two different cases that are studied in detail: (i) a quasiparticle (resonance) when the renormalized frequency of the particle is above the frequency threshold of the bath and (ii) a stable renormalized "particle" state below this threshold. The time evolution of the occupation number for the particle is evaluated exactly using different approaches that yield to complementary insights. The exact solution allows us to investigate the concept of the formation time of a quasiparticle and to study the difference between the relaxation of the distribution of bare particles and that of quasiparticles. For the case of quasiparticles, the exact occupation number asymptotically tends to a statistical equilibrium distribution that differs from a simple Bose-Einstein form as a result of off-shell processes whereas in the stable particle case, the distribution of particles does not thermalize with the bath. We derive a non-Markovian quantum kinetic equation which resums the perturbative series and includes off-shell effects. A Markovian approximation that includes off-shell contributions and the usual Boltzmann equation (energy conserving) are obtained from the quantum kinetic equation in the limit of wide separation of time scales upon different coarse-graining assumptions. The relaxational dynamics predicted by the non-Markovian, Markovian, and Boltzmann approximations are compared to the exact result. The Boltzmann approach is seen to fail in the case of wide resonances and when threshold and renormalization effects are important.
For many open quantum systems, a master equation approach employing the Markov approximation cannot reliably describe the dynamical behaviour. This is the case, for example, in a number of solid state or biological systems, and it has motivated a line of research aimed at quantifying the amount of non-Markovian behaviour in a given model. Within this framework, we investigate the dynamics of a quantum harmonic oscillator linearly coupled to a bosonic bath. We focus on Gaussian states, which are suitably treated using a covariance matrix approach. Concentrating on an entanglement based non-Markovian behaviour quantifier (NMBQ) proposed by Rivas et. al. [1], we consider the role that near resonant and off-resonant modes play in affecting the NMBQ. By using a large but finite bath of oscillators for both Ohmic and super Ohmic spectral densities we find, by systematically increasing the coupling strength, initially the near resonant modes provide the most significant non-Markovian effects, while after a certain threshold of coupling strength the off-resonant modes play the dominant role. We also consider the NMBQ for two other models where we add a single strongly coupled oscillator to the model in extra bath mode and 'buffer' configurations, which affects the modes that determine non-Markovian behaviour.
We investigate the effect of inter-mode coupling on the vibrational relaxation dynamics of molecules in weak dissipative environments. The simulations are performed within the reduced density matrix formalism in the Markovian regime, assuming a Lindblad form for the system-bath interaction. The prototypical two-dimensional model system representing two CO molecules approaching a Cu(100) surface is adapted from an ab initio potential, while the diatom-diatom vibrational coupling strength is systematically varied. In the weak system-bath coupling limit and at low temperatures, only first order non-adiabatic uni-modal coupling terms contribute to surface-mediated vibrational relaxation. Since dissipative dynamics is non-unitary, the choice of representation will affect the evolution of the reduced density matrix. Two alternative representations for computing the relaxation rates and the associated operators are thus compared: the fully coupled spectral basis, and a factorizable ansatz. The former is well-established and serves as a benchmark for the solution of Liouville-von Neumann equation. In the latter, a contracted grid basis of potential-optimized discrete variable representation is tailored to incorporate most of the inter-mode coupling, while the Lindblad operators are represented as tensor products of one-dimensional operators, for consistency. This procedure results in a marked reduction of the grid size and in a much more advantageous scaling of the computational cost with respect to the increase of the dimensionality of the system. The factorizable method is found to provide an accurate description of the dissipative quantum dynamics of the model system, specifically of the time evolution of the state populations and of the probability density distribution of the molecular wave packet. The influence of intra-molecular vibrational energy redistribution appears to be properly taken into account by the new model on the whole range of coupling strengths. It demontrates that most of the mode mixing during relaxation is due to the potential part of the Hamiltonian and not to the coupling among relaxation operators.
Non-Markovian effects are important in modeling the behavior of open quantum systems arising in solid-state physics, quantum optics as well as in study of biological and chemical systems. The non-Markovian environment is often approximated by discrete bosonic modes, thus mapping it to a Lindbladian or Hamiltonian simulation problem. While systematic constructions of such modes have been previously proposed, the resulting approximation lacks rigorous and general convergence guarantees. In this Letter, we show that under some physically motivated assumptions on the system-environment interaction, the finite-time dynamics of the non-Markovian open quantum system computed with a sufficiently large number of modes is guaranteed to converge to the true result. Furthermore, we show that this approximation error typically falls off polynomially with the number of modes. Our results lend rigor to classical and quantum algorithms for approximating non-Markovian dynamics.