The Journal of Chemical Physics

A semiclassical approach to dissipation in quantum mechanics

Frank Großmann

Published Sep 1, 1995• 27 citations

Abstract

Employing the path integral formalism, we study a quantum (test) system coupled to an environment consisting of infinitely many harmonic oscillators. A simple semiclassical approxi- mation, in which only real classical trajectories are required, is used to derive the propagator of the reduced density matrix for Ohmic dissipation and high temperatures. The inclusion of a summation over trajectories with the correct Maslov phases in the final expression for the propagator permits investigation of nonharmonic test systems. In a numerical example for a Morse oscillator interacting with a high temperature heat bath, the semiclassical result, obtained by the method described here, correctly displays the suppression of quantum behavior as observed in the exact calculation.

DOI: 10.1063/1.470046

Related papers

Decoherence and dissipation in a molecular system coupled to an environment: An application of semiclassical hybrid dynamics

Christoph-Marian Goletz, Frank Großmann

Applying the recently developed semiclassical hybrid dynamics [Grossmann, J. Chem. Phys. 125, 014111 (2006)], we study the decay of interference patterns in the reduced density as well as of the purity in a Morse oscillator test system due to the interaction with a finite harmonic bath at zero temperature. In the case that the bath mimics a continuous Ohmic spectral density, in addition to the quantum classical transition induced by the interaction with the environment, we corroborate the existence of a blueshift due to the bath coupling, predicted by Pollak [Phys. Rev. A 33, 4244 (1986)]. Furthermore, the decoherence dynamics of cat states is confirmed to be faster than that of single coherent states and we show that for a resonant bath the dissipation leads to an increase in the decoherence rate as compared to the low frequency bath.

The Journal of Chemical PhysicsJun 25, 200929 citations

Path integral formulation of dissipative quantum dynamics

Alexey Novikov

In this thesis the path integral formalism is applied to the calculation of the dynamics of dissipative quantum systems. The time evolution of a system of bilinearly coupled bosonic modes is treated using the real-time path integral technique in coherent-state representation. This method is applied to a damped harmonic oscillator within the Caldeira-Leggett model. In order to get the stationary trajectories the corresponding Lagrangian function is diagonalized and then the path integrals are evaluated by means of the stationary-phase method. The time evolution of the reduced density matrix in the basis of coherent states is given in simple analytic form for weak system-bath coupling, i.e. the so-called rotating-wave terms can be evaluated exactly but the non-rotating-wave terms only in a perturbative manner. The validity range of the rotating-wave approximation is discussed from the viewpoint of spectral equations. In addition, it is shown that systems without initial system-bath correlations can exhibit initial jumps in the population dynamics even for rather weak dissipation. Only with initial correlations the classical trajectories for the system coordinate can be recovered. The path integral formalism in a combined phase-space and coherent-state representation is applied to the problem of curve-crossing dynamics. The system of interest is described by two coupled one-dimensional harmonic potential energy surfaces interacting with a heat bath. The mapping approach is used to rewrite the Lagrangian function of the electronic part of the system. Using the Feynman-Vernon influence-functional method the bath is eliminated whereas the non-Gaussian part of the path integral is treated using the perturbation theory in the small coordinate shift between potential energy surfaces. The vibrational and the population dynamics is considered in a lowest order of the perturbation. The dynamics of a Gaussian wave packet is analyzed along a one-dimensional reaction coordinate. Also the damping rate of coherence in the electronic part of the relevant system is evaluated within the ordinary and variational perturbation theory. The analytic expressions for the rate functions are obtained in the low and high temperature regimes.

Jan 20, 20050 citations

Dissipative quantum dynamics: Driven molecular vibrations

Mingwhei Tung, Jian‐Min Yuan

We have studied quantum-mechanical behavior of a driven Morse oscillator coupled to a bath of harmonic oscillators. The purpose is to compare the quantum behavior of such a system with the classical solutions of a driven damped Morse oscillator. We start with the quantum Liouville equation, in which the Hamiltonian of the Morse oscillator is expressed in terms of generators of an su(2) Lie algebra. This algebra and the Markovian approximation allow us to derive the generalized master equation (for the reduced density matrix of the Morse oscillator), which contains level-dependent energy and phase relaxation terms. We have numerically integrated the differential equations of the matrix elements to obtain the time evolution of the reduced density matrix and the energy of the Morse oscillator. We show that the energy of the Morse oscillator in general varies with time and eventually reaches an asymptotic oscillatory state. The mean energy value of the asymptotic oscillatory state is studied as a function of the relaxation rates and laser frequency and amplitude. Vibrational distributions have also been found as functions of laser frequency and amplitude. The shape of the distribution changes gradually as laser amplitude increases; it first peaks at the ground state at a small field amplitude, then it peaks at some excited state at a large field amplitude, and for an amplitude in between it can have two maxima. This bimodal vibrational distribution reflects the bistability observed in the classical and semiclassical models. Thus bistability exists in quantum results, not as a hysteresis loop, but as a bimodal distribution. Finally, we show that as laser intensity increases, the time series of the oscillator energy evolves from regular to quasiperiodic behavior, and eventually a chaotic-looking series appears.

Physical Review ANov 1, 198725 citations

Non-Markovian Dissipative Quantum Mechanics with Stochastic Trajectories

Werner Koch

All fields of physics - be it nuclear, atomic and molecular, solid state, or optical - offer examples of systems which are strongly influenced by the environment of the actual system under investigation. The scope of what is called "the environment" may vary, i.e., how far from the system of interest an interaction between the two does persist. Typically, however, it is much larger than the open system itself. Hence, a fully quantum mechanical treatment of the combined system without approximations and without limitations of the type of system is currently out of reach. With the single assumption of the environment to consist of an internally thermalized set of infinitely many harmonic oscillators, the seminal work of Stockburger and Grabert [Chem. Phys., 268:249-256, 2001] introduced an open system description that captures the environmental influence by means of a stochastic driving of the reduced system. The resulting stochastic Liouville-von Neumann equation describes the full non-Markovian dynamics without explicit memory but instead accounts for it implicitly through the correlations of the complex-valued noise forces. The present thesis provides a first application of the Stockburger-Grabert stochastic Liouville-von Neumann equation to the computation of the dynamics of anharmonic, continuous open systems. In particular, it is demonstrated that trajectory based propagators allow for the construction of a numerically stable propagation scheme. With this approach it becomes possible to achieve the tremendous increase of the noise sample count necessary to stochastically converge the results when investigating such systems with continuous variables. After a test against available analytic results for the dissipative harmonic oscillator, the approach is subsequently applied to the analysis of two different realistic, physical systems. As a first example, the dynamics of a dissipative molecular oscillator is investigated. Long time propagation - until thermalization is reached - is shown to be possible with the presented approach. The properties of the thermalized density are determined and they are ascertained to be independent of the system's initial state. Furthermore, the dependence on the bath's temperature and coupling strength is analyzed and it is demonstrated how a change of the bath parameters can be used to tune the system from the dissociative to the bound regime. A second investigation is conducted for a dissipative tunneling scenario in which a wave packet impinges on a barrier. The dependence of the transmission probability on the initial state's kinetic energy as well as the bath's temperature and coupling strength is computed. For both systems, a comparison with the high-temperature Markovian quantum Brownian limit is performed. The importance of a full non-Markovian treatment is demonstrated as deviations are shown to exist between the two descriptions both in the low temperature cases where they are expected and in some of the high temperature cases where their appearance might not be anticipated as easily.:1 Introduction. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9 2 Theory of Open Quantum Systems . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13 2.1 Influence Functional Formalism . . . . . . . . . . . . . . . . . . . . . . . . . 14 2.2 Quantum Brownian Limit . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18 2.3 Stochastic Unraveling of the Influence Functional . . . . . . . . . . . . . . . 20 2.4 Improved Statistics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23 2.4.1 Modified Dynamic Response . . . . . . . . . . . . . . . . . . . . . . . 23 2.4.2 Guide Trajectory Transformation . . . . . . . . . . . . . . . . . . . . 24 2.5 Obtaining Properly Correlated Stochastic Samples from Filtered White Noise 24 3 Unified Stochastic Trajectory Propagation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27 3.1 Semiclassical Brownian Motion . . . . . . . . . . . . . . . . . . . . . . . . . 28 3.1.1 Guide Trajectory . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30 3.1.2 Real Coherent State Center Coordinates . . . . . . . . . . . . . . . . 31 3.1.3 Propagation Scheme Including Stochastic Forces . . . . . . . . . . . 32 3.2 Stochastic Bohmian Mechanics with Complex Action . . . . . . . . . . . . . 33 3.2.1 Hydrodynamic Formulation of Bohmian Mechanics . . . . . . . . . . 33 3.2.2 Bohmian Mechanics with Complex Action . . . . . . . . . . . . . . . 34 3.2.3 Stochastic BOMCA Trajectories . . . . . . . . . . . . . . . . . . . . 38 3.3 Noise Distribution Preserving Removal of Adverse Samples . . . . . . . . . . 39 4 Dissipative Harmonic Oscillator. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41 4.1 Reservoir Specification . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41 4.2 Harmonic Oscillator Analytic Expectation Values . . . . . . . . . . . . . . . 42 4.2.1 Ohmic Bath . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42 4.2.2 Drude Regularized Bath . . . . . . . . . . . . . . . . . . . . . . . . . 43 4.3 Sampling Strategies and Analytic Comparison . . . . . . . . . . . . . . . . . 44 4.4 Limits of the Markovian Approximation . . . . . . . . . . . . . . . . . . . . 45 5 Dissipative Vibrational Dynamics of Diatomics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49 5.1 Molecular Morse Potential . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49 5.2 Anharmonic Phenomena . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51 5.3 Transient Non-Markovian Effects . . . . . . . . . . . . . . . . . . . . . . . . 53 5.4 Trapping by Dissipation and Thermalization . . . . . . . . . . . . . . . . . . 53 6 Tunneling with Dissipation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57 6.1 Eckart Barrier Scattering . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58 6.2 Dissipative Tunneling . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59 6.3 Investigation of Markovianity . . . . . . . . . . . . . . . . . . . . . . . . . . 61 7 Summary and Outlook . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65 Appendix A Conventions for Constants, Reservoir Kernels, and Influence Phases 69 Appendix B Stochastic Calculus. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71 B.1 Stochastic Differential Equations . . . . . . . . . . . . . . . . . . . . . . . . 71 B.2 Position Verlet Scheme . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72 B.3 Runge-Kutta Fourth Order Scheme . . . . . . . . . . . . . . . . . . . . . . . 73 Appendix CMorse Oscillator Expectation Values . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75 Appendix DPrerequisites of a Successful Stochastic Propagation . . . . . . . . . . . . . . 79 D.1 Hubbard-Stratonovich Transformation . . . . . . . . . . . . . . . . . . . . . 79 D.2 Kernels of the Reservoir . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82 D.2.1 Quadratic Cutoff . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82 D.2.2 Quartic Cutoff . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86 D.2.3 Strictly Ohmic Reservoir . . . . . . . . . . . . . . . . . . . . . . . . . 89 D.3 Guide Trajectory Integration . . . . . . . . . . . . . . . . . . . . . . . . . . 90 D.3.1 Quadratic Cutoff . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90 D.3.2 Quartic Cutoff . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91 D.3.3 Strictly Ohmic Reservoir . . . . . . . . . . . . . . . . . . . . . . . . . 92 D.4 Computation of Matrix Elements and Expectation Values . . . . . . . . . . 92 Bibliography . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97

Sep 9, 20101 citation

Quantum dissipation and decay in a medium

I. Joichi, Sh. Matsumoto, M. Yoshimura

Quantum dissipation in a thermal environment is investigated, using the path-integral approach. The reduced density matrix of the harmonic oscillator system coupled to a thermal bath of oscillators is derived for an arbitrary spectrum of bath oscillators. Time evolution and the end point of the two-body decay of unstable particles are then elucidated: After early transient times unstable particles undergo the exponential decay, followed by the power-law decay and finally ending in a mixed state of residual particles containing contributions from both on and off the mass shell, whose abundance does not suffer from the Boltzmann suppression.

arXiv:hep-th/9609223Physical Review AFeb 1, 199829 citations

Dissipative quantum and classical dynamics: driven molecular vibration

Jian-Min Yuan, Mingwhei Tung

We study quantum mechanical behavior of a Morse oscillator coupled to a bath of harmonic oscillators and driven by an infrared laser. The purpose is to compare quantum behavior of such a system with the classical and semiclassical solutions of a driven damped Morse oscillator. The general picture coming out of classical and semiclassical studies is that the total energy of the Morse oscillator in steady states or steady oscillations show jump and hysteretic behavior, as expected in a cusp catastrophe, when laser intensity or frequency is varied. Along with the bistable behavior we have also found period-doubling bifurcations leading to chaos on one or both branches. In the quantum treatment a generalized master equation is derived from the Liouville equation, in which the Hamiltonian of the Morse oscillator is expressed in terms of generators of an SU(2) algebra. This algebra and the Markoffian approximation allows us to derive the equation of motion of the reduced density matrix, which is solved numerically. Bistable behavior does not exist in this quantum treatment, but solutions do show jump behavior and bifurcation phenomenon. Experimental implications of these studies are also discussed.

International Laser Science ConferenceJan 1, 19860 citations