Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2002.07700 · Physical review. B./Physical review. B

Exactly thermalized quantum dynamics of the spin-boson model coupled to a dissipative environment

M. A. Lane, D. Matos, Ian J. Ford, Lev N Kantorovich

Published Jun 22, 2020• 5 citations
arXivHTMLPDF

Abstract

We present an application of the extended stochastic Liouville--von Neumann equation (ESLN) method introduced earlier [G. M. G. McCaul, C. D. Lorenz, and L. Kantorovich, Phys. Rev. B 95, 125124 (2017); G. M. G. McCaul, C. D. Lorenz, and L. Kantorovich, Phys. Rev. B 97, 224310 (2018)], which describes the dynamics of an exactly thermalized open quantum system reduced density matrix coupled to a non-Markovian harmonic environment. Critically, the combined system of the open system fully coupled to its environment is thermalized at finite temperature using an imaginary-time evolution procedure before the application of real-time evolution. This initializes the combined system in the correct canonical equilibrium state rather than being initially decoupled. Here we apply our theory to the spin-boson Hamiltonian and develop a number of competing ESLN variants designed to reduce the numerical divergence of the trace of the open-system density matrix. We find that a careful choice of the driving noises is essential for improving numerical stability. We also investigate the effect of applying higher-order numerical schemes for solving stochastic differential equations, such as the Stratonovich-Heun scheme, and conclude that stochastic sampling dominates convergence with the improvement associated with the numerical scheme being less important for short times but required for late times. To verify the method and its numerical implementation, we first consider evolution under a fixed Hamiltonian and show that the system either remains in, or approaches, the correct canonical equilibrium state at long times. Additionally, evolution of the open system under nonequilibrium Landau-Zener (LZ) driving is considered and the asymptotic convergence to the LZ limit was observed for vanishing system-environment coupling and temperature. When coupling and temperature are nonzero, initially thermalizing the combined system at a finite time in the past was found to be a better approximation of the true LZ initial state than starting in a pure state.

DOI
10.1103/physrevb.101.224306

Related papers

Driving spin-boson models from equilibrium using exact quantum dynamics

Gerard McCaul, Christian D. Lorenz, Lev N Kantorovich

We present an application of the extended stochastic Liouville equation [G. M. G. McCaul et al., Phys. Rev. B 95, 125124 (2017)], which gives an exact solution for the reduced density matrix of an open system surrounded by a harmonic heat bath. This method considers the extended system (the open system and the bath) being thermally equilibrated prior to the action of a time-dependent perturbation, as opposed to the usual assumption that the system and the bath are initially partitioned. Our method presents an exact technique capable of accounting for arbitrary parameter regimes of the model. Here, we present our first numerical implementation of the method in the simplest case of a Caldeira-Leggett representation of the bath Hamiltonian, and apply it to a spin-boson system driven from coupled equilibrium. We observe significant behaviors in both the transient dynamics and asymptotic states of the reduced density matrix not present in the usual approximation.

arXiv:1712.06397•Physical review. B./Physical review. B•Jun 25, 2018•8 citations

Nonperturbative stochastic method for driven spin-boson model

Peter P. Orth, Adilet Imambekov, Karyn Le Hur

We introduce and apply a numerically exact method for investigating the real-time dissipative dynamics of quantum impurities embedded in a macroscopic environment beyond the weak-coupling limit. We focus on the spin-boson Hamiltonian that describes a two-level system interacting with a bosonic bath of harmonic oscillators. This model is archetypal for investigating dissipation in quantum systems, and tunable experimental realizations exist in mesoscopic and cold-atom systems. It finds abundant applications in physics ranging from the study of decoherence in quantum computing and quantum optics to extended dynamical mean-field theory. Starting from the real-time Feynman-Vernon path integral, we derive an exact stochastic Schr\"odinger equation that allows us to compute the full spin density matrix and spin-spin correlation functions beyond weak coupling. We greatly extend our earlier work [P. P. Orth, A. Imambekov, and K. Le Hur, Phys. Rev. A 82, 032118 (2010)] by fleshing out the core concepts of the method and by presenting a number of interesting applications. Methodologically, we present an analogy between the dissipative dynamics of a quantum spin and that of a classical spin in a random magnetic field. This analogy is used to recover the well-known noninteracting-blip approximation in the weak-coupling limit. We explain in detail how to compute spin-spin autocorrelation functions. As interesting applications of our method, we explore the non-Markovian effects of the initial spin-bath preparation on the dynamics of the coherence ${\ensuremath{\sigma}}^{x}(t)$ and of ${\ensuremath{\sigma}}^{z}(t)$ under a Landau-Zener sweep of the bias field. We also compute to a high precision the asymptotic long-time dynamics of ${\ensuremath{\sigma}}^{z}(t)$ without bias and demonstrate the wide applicability of our approach by calculating the spin dynamics at nonzero bias and different temperatures.

arXiv:1211.1201•Physical Review B•Jan 25, 2013•83 citations

Non-perturbative stochastic method for driven quantum impurity systems

Peter P. Orth, Adilet Imambekov, Karyn Le Hur

We introduce and apply a numerically exact method for investigating the real-time dissipative dynamics of quantum impurities embedded in a macroscopic environment beyond the weak-coupling limit. We focus on the spin-boson Hamiltonian that describes a two-level system interacting with a bosonic bath of harmonic oscillators. This model is archetypal for investigating dissipation in quantum systems and tunable experimental realizations exist in mesoscopic and cold-atom systems. It finds abundant applications in physics ranging from the study of decoherence in quantum computing and quantum optics to extended dynamical mean-field theory. Starting from the real-time Feynman-Vernon path integral, we derive an exact stochastic Schr\odinger equation that allows to compute the full spin density matrix and spin-spin correlation functions beyond weak coupling. We greatly extend our earlier work (P. P. Orth, A. Imambekov, K. Le Hur, Phys. Rev. A {\bf 82}, 032118 (2010)) by fleshing out the core concepts of the method and by presenting a number of interesting applications. Methodologically, we present an analogy between the dissipative dynamics of a quantum spin and that of a classical spin in a random magnetic field. This analogy is used to recover the well-known Non-Interacting-Blip-Approximation (NIBA) in the weak-coupling limit. We explain in detail how to compute spin-spin autocorrelation functions. As interesting applications of our method, we explore the non-Markovian effects of the initial spin-bath preparation on the dynamics of the coherence $\sigma^x(t)$ and of $\sigma^z(t)$ under a Landau-Zener sweep of the bias field. We also compute to a high precision the asymptotic long-time dynamics of $\sigma^z(t)$ without bias and demonstrate the wide applicability of our approach by calculating the spin dynamics at non-zero bias and different temperatures.

•Nov 6, 2012•0 citations

Stochastic simulation of dissipation and non-Markovian effects in open quantum systems

Denis Lacroix A

The exact dynamics of a system coupled to an environment can be described by an integro-differential stochastic equation for the reduced density. The influence of the environment is incorporated through a mean field which is both stochastic and nonlocal in time and where the standard two-time correlation functions of the environment appear naturally. Since no approximation is made, the presented theory incorporates exactly dissipative and non-Markovian effects. Applications to the spin-boson model coupled to a heat bath with various coupling regimes and temperature show that the presented stochastic theory can be a valuable tool to simulate exactly the dynamics of open quantum systems. Links with the stochastic Schrödinger equation method and possible extensions to "imaginary time" propagation are discussed.

arXiv:0802.1981•Physical Review E•Apr 29, 2008•22 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, 2010•1 citation

Non-Markovian stochastic Schrödinger equations in different temperature regimes: A study of the spin-boson model

Inés de Vega, Daniel Alonso, Pierre Gaspard +1 more

Stochastic Schrodinger equations are used to describe the dynamics of a quantum open system in contact with a large environment, as an alternative to the commonly used master equations. We present a study of the two main types of non-Markovian stochastic Schrodinger equations, linear and nonlinear ones. We compare them both analytically and numerically, the latter for the case of a spin-boson model. We show in this paper that two linear stochastic Schrodinger equations, derived from different perspectives by Diosi, Gisin, and Strunz [Phys. Rev. A 58, 1699 (1998)], and Gaspard and Nagaoka [J. Chem. Phys. 13, 5676 (1999)], respectively, are equivalent in the relevant order of perturbation theory. Nonlinear stochastic Schrodinger equations are in principle more efficient than linear ones, as they determine solutions with a higher weight in the ensemble average which recovers the reduced density matrix of the quantum open system. However, it will be shown in this paper that for the case of a spin-boson system and weak coupling, this improvement does only occur in the case of a bath at high temperature. For low temperatures, the sampling of realizations of the nonlinear equation is practically equivalent to the sampling of the linear ones. We study further this result by analyzing, for both temperature regimes, the driving noise of the linear equations in comparison to that of the nonlinear equations.

The Journal of Chemical Physics•Mar 22, 2005•40 citations
We present an application of the extended stochastic Liouville--von Neumann equation (ESLN) method introduced earlier [G. M. G. McCaul, C. D. Lorenz, and L. Kantorovich, Phys. Rev. B 95, 125124 (2017); G. M. G. McCaul, C. D. Lorenz, and L. Kantorovich, Phys. Rev. B 97, 224310 (2018)], which describes the dynamics of an exactly thermalized open quantum system reduced density matrix coupled to a non-Markovian harmonic environment. Critically, the combined system of the open system fully coupled to its environment is thermalized at finite temperature using an imaginary-time evolution procedure before the application of real-time evolution. This initializes the combined system in the correct canonical equilibrium state rather than being initially decoupled. Here we apply our theory to the spin-boson Hamiltonian and develop a number of competing ESLN variants designed to reduce the numerical divergence of the trace of the open-system density matrix. We find that a careful choice of the driving noises is essential for improving numerical stability. We also investigate the effect of applying higher-order numerical schemes for solving stochastic differential equations, such as the Stratonovich-Heun scheme, and conclude that stochastic sampling dominates convergence with the improvement associated with the numerical scheme being less important for short times but required for late times. To verify the method and its numerical implementation, we first consider evolution under a fixed Hamiltonian and show that the system either remains in, or approaches, the correct canonical equilibrium state at long times. Additionally, evolution of the open system under nonequilibrium Landau-Zener (LZ) driving is considered and the asymptotic convergence to the LZ limit was observed for vanishing system-environment coupling and temperature. When coupling and temperature are nonzero, initially thermalizing the combined system at a finite time in the past was found to be a better approximation of the true LZ initial state than starting in a pure state.
We present an application of the extended stochastic Liouville equation [G. M. G. McCaul et al., Phys. Rev. B 95, 125124 (2017)], which gives an exact solution for the reduced density matrix of an open system surrounded by a harmonic heat bath. This method considers the extended system (the open system and the bath) being thermally equilibrated prior to the action of a time-dependent perturbation, as opposed to the usual assumption that the system and the bath are initially partitioned. Our method presents an exact technique capable of accounting for arbitrary parameter regimes of the model. Here, we present our first numerical implementation of the method in the simplest case of a Caldeira-Leggett representation of the bath Hamiltonian, and apply it to a spin-boson system driven from coupled equilibrium. We observe significant behaviors in both the transient dynamics and asymptotic states of the reduced density matrix not present in the usual approximation.
We introduce and apply a numerically exact method for investigating the real-time dissipative dynamics of quantum impurities embedded in a macroscopic environment beyond the weak-coupling limit. We focus on the spin-boson Hamiltonian that describes a two-level system interacting with a bosonic bath of harmonic oscillators. This model is archetypal for investigating dissipation in quantum systems, and tunable experimental realizations exist in mesoscopic and cold-atom systems. It finds abundant applications in physics ranging from the study of decoherence in quantum computing and quantum optics to extended dynamical mean-field theory. Starting from the real-time Feynman-Vernon path integral, we derive an exact stochastic Schr\"odinger equation that allows us to compute the full spin density matrix and spin-spin correlation functions beyond weak coupling. We greatly extend our earlier work [P. P. Orth, A. Imambekov, and K. Le Hur, Phys. Rev. A 82, 032118 (2010)] by fleshing out the core concepts of the method and by presenting a number of interesting applications. Methodologically, we present an analogy between the dissipative dynamics of a quantum spin and that of a classical spin in a random magnetic field. This analogy is used to recover the well-known noninteracting-blip approximation in the weak-coupling limit. We explain in detail how to compute spin-spin autocorrelation functions. As interesting applications of our method, we explore the non-Markovian effects of the initial spin-bath preparation on the dynamics of the coherence ${\ensuremath{\sigma}}^{x}(t)$ and of ${\ensuremath{\sigma}}^{z}(t)$ under a Landau-Zener sweep of the bias field. We also compute to a high precision the asymptotic long-time dynamics of ${\ensuremath{\sigma}}^{z}(t)$ without bias and demonstrate the wide applicability of our approach by calculating the spin dynamics at nonzero bias and different temperatures.
We introduce and apply a numerically exact method for investigating the real-time dissipative dynamics of quantum impurities embedded in a macroscopic environment beyond the weak-coupling limit. We focus on the spin-boson Hamiltonian that describes a two-level system interacting with a bosonic bath of harmonic oscillators. This model is archetypal for investigating dissipation in quantum systems and tunable experimental realizations exist in mesoscopic and cold-atom systems. It finds abundant applications in physics ranging from the study of decoherence in quantum computing and quantum optics to extended dynamical mean-field theory. Starting from the real-time Feynman-Vernon path integral, we derive an exact stochastic Schr\odinger equation that allows to compute the full spin density matrix and spin-spin correlation functions beyond weak coupling. We greatly extend our earlier work (P. P. Orth, A. Imambekov, K. Le Hur, Phys. Rev. A {\bf 82}, 032118 (2010)) by fleshing out the core concepts of the method and by presenting a number of interesting applications. Methodologically, we present an analogy between the dissipative dynamics of a quantum spin and that of a classical spin in a random magnetic field. This analogy is used to recover the well-known Non-Interacting-Blip-Approximation (NIBA) in the weak-coupling limit. We explain in detail how to compute spin-spin autocorrelation functions. As interesting applications of our method, we explore the non-Markovian effects of the initial spin-bath preparation on the dynamics of the coherence σx(t)\sigma^x(t)σx(t) and of σz(t)\sigma^z(t)σz(t) under a Landau-Zener sweep of the bias field. We also compute to a high precision the asymptotic long-time dynamics of σz(t)\sigma^z(t)σz(t) without bias and demonstrate the wide applicability of our approach by calculating the spin dynamics at non-zero bias and different temperatures.
The exact dynamics of a system coupled to an environment can be described by an integro-differential stochastic equation for the reduced density. The influence of the environment is incorporated through a mean field which is both stochastic and nonlocal in time and where the standard two-time correlation functions of the environment appear naturally. Since no approximation is made, the presented theory incorporates exactly dissipative and non-Markovian effects. Applications to the spin-boson model coupled to a heat bath with various coupling regimes and temperature show that the presented stochastic theory can be a valuable tool to simulate exactly the dynamics of open quantum systems. Links with the stochastic Schrödinger equation method and possible extensions to "imaginary time" propagation are discussed.
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
Stochastic Schrodinger equations are used to describe the dynamics of a quantum open system in contact with a large environment, as an alternative to the commonly used master equations. We present a study of the two main types of non-Markovian stochastic Schrodinger equations, linear and nonlinear ones. We compare them both analytically and numerically, the latter for the case of a spin-boson model. We show in this paper that two linear stochastic Schrodinger equations, derived from different perspectives by Diosi, Gisin, and Strunz [Phys. Rev. A 58, 1699 (1998)], and Gaspard and Nagaoka [J. Chem. Phys. 13, 5676 (1999)], respectively, are equivalent in the relevant order of perturbation theory. Nonlinear stochastic Schrodinger equations are in principle more efficient than linear ones, as they determine solutions with a higher weight in the ensemble average which recovers the reduced density matrix of the quantum open system. However, it will be shown in this paper that for the case of a spin-boson system and weak coupling, this improvement does only occur in the case of a bath at high temperature. For low temperatures, the sampling of realizations of the nonlinear equation is practically equivalent to the sampling of the linear ones. We study further this result by analyzing, for both temperature regimes, the driving noise of the linear equations in comparison to that of the nonlinear equations.