A systematic procedure for deriving the master equation of a dissipative system is reported in the framework of stochastic description. For the Caldeira-Leggett model of the harmonic-oscillator bath, a detailed and elementary derivation of the bath-induced stochastic field is presented. The dynamics of the system is thereby fully described by a stochastic differential equation, and the desired master equation would be acquired with statistical averaging. It is shown that the existence of a closed-form master equation depends on the specificity of the system as well as the feature of the dissipation characterized by the spectral density function. For a dissipative harmonic oscillator it is observed that the correlation between the stochastic field due to the bath and the system can be decoupled, and the master equation naturally results. Such an equation possesses the Lindblad form in which time-dependent coefficients are determined by a set of integral equations. It is proved that the obtained master equation is equivalent to the well-known Hu-Paz-Zhang equation based on the path-integral technique. The procedure is also used to obtain the master equation of a dissipative harmonic oscillator in time-dependent fields.
A systematic procedure for deriving the master equation of a dissipative system is reported in the framework of stochastic description. For the Caldeira-Leggett model of the harmonic-oscillator bath, a detailed and elementary derivation of the bath-induced stochastic field is presented. The dynamics of the system is thereby fully described by a stochastic differential equation, and the desired master equation would be acquired with statistical averaging. It is shown that the existence of a closed-form master equation depends on the specificity of the system as well as the feature of the dissipation characterized by the spectral density function. For a dissipative harmonic oscillator it is observed that the correlation between the stochastic field due to the bath and the system can be decoupled, and the master equation naturally results. Such an equation possesses the Lindblad form in which time-dependent coefficients are determined by a set of integral equations. It is proved that the obtained master equation is equivalent to the well-known Hu-Paz-Zhang equation based on the path-integral technique. The procedure is also used to obtain the master equation of a dissipative harmonic oscillator in time-dependent fields.
The methodology of stochastic description for dissipation, a generic scheme to decouple the interaction between two subsystems, is applied to the study of dissipative dynamics in quantum optics. It is shown that the influence of the coupled thermal or vacuum field on the quantum mode can be exactly represented by the induced stochastic fields. The quantum mode thereby satisfies a stochastic differential equation and dissipation effect due to the coupling with the environment is obtained through statistical averaging. Within the framework of stochastic description, it is demonstrated how to derive the master equation for a single optical mode interacting with the bosonic bath. A numerical algorithm for solving the master equation in which the coefficients are determined by a set of integral equations is discussed and a comparison with the known results is displayed. The derivation of the master equation for the spontaneous decay of two-state atoms in the vacuum is also presented.
Time evolution of a harmonic oscillator linearly coupled to a heat bath is compared for three classes of initial states for the bath modes-grand canonical ensemble, number states, and coherent states. It is shown that for a wide class of number states the behavior of the oscillator is similar to the case of the equilibrium bath. If the bath modes are initially in coherent states, then the variances of the oscillator coordinate and momentum, as well as its entanglement to the bath, asymptotically approach the same values as for the oscillator at zero temperature and the average coordinate and momentum show a Brownian-like behavior. We derive an exact master equation for the characteristic function of the oscillator valid for arbitrary factorized initial conditions. In the case of the equilibrium bath this equation reduces to an equation of the Hu-Paz-Zhang type, while for the coherent states bath it leads to an exact stochastic master equation with a multiplicative noise.
Driven nonlinear quantum oscillators are a central platform for quantum technologies, yet their dissipative dynamics are typically described using Lindblad or Caldeira-Leggett master equations derived under assumptions that exclude nonlinearities and driving. Here, we derive a generalized Caldeira-Leggett master equation for driven nonlinear oscillators by retaining the full nonlinear and time-dependent system dynamics in the construction of the dissipator. For position- and momentum-dependent system-bath coupling, the dissipator itself becomes dynamically dressed, generating nonlinear and drive-dependent dissipative channels beyond conventional fixed-dissipator approaches. This produces nonlinear damping together with dissipation-induced corrections to the effective drive. The resulting dissipative dynamics suppress large-amplitude excitations and reduce phase-space fluctuations. For a driven Kerr oscillator, this leads to the suppression of bistability, asymmetric resonance responses, and strongly modified fluctuation distributions. More broadly, our results establish a microscopic framework in which nonlinear dynamics and driving directly reshape the dissipative sector of driven open quantum systems.
Classical Brownian motion has well been investigated since the pioneering work of Einstein, which inspired mathematicians to lay the theoretical foundation of stochastic processes. A stochastic formulation for quantum dynamics of dissipative systems described by the system-plus-bath model has been developed and found many applications in chemical dynamics, spectroscopy, quantum transport, and other fields. This article provides a tutorial review of the stochastic formulation for quantum dissipative dynamics. The key idea is to decouple the interaction between the system and the bath by virtue of the Hubbard-Stratonovich transformation or Itô calculus so that the system and the bath are not directly entangled during evolution, rather they are correlated due to the complex white noises introduced. The influence of the bath on the system is thereby defined by an induced stochastic field, which leads to the stochastic Liouville equation for the system. The exact reduced density matrix can be calculated as the stochastic average in the presence of bath-induced fields. In general, the plain implementation of the stochastic formulation is only useful for short-time dynamics, but not efficient for long-time dynamics as the statistical errors go very fast. For linear and other specific systems, the stochastic Liouville equation is a good starting point to derive the master equation. For general systems with decomposable bath-induced processes, the hierarchical approach in the form of a set of deterministic equations of motion is derived based on the stochastic formulation and provides an effective means for simulating the dissipative dynamics. A combination of the stochastic simulation and the hierarchical approach is suggested to solve the zero-temperature dynamics of the spin-boson model. This scheme correctly describes the coherent-incoherent transition (Toulouse limit) at moderate dissipation and predicts a rate dynamics in the overdamped regime. Challenging problems such as the dynamical description of quantum phase transition (local- ization) and the numerical stability of the trace-conserving, nonlinear stochastic Liouville equation are outlined.
In this paper we derive an exact master equation for two coupled quantum harmonic oscillators interacting via bilinear coupling with a common environment at arbitrary temperature made up of many harmonic oscillators with a general spectral density function. We first show a simple derivation based on the observation that the two harmonic oscillator model can be effectively mapped into that of a single harmonic oscillator in a general environment plus a free harmonic oscillator. Since the exact one harmonic oscillator master equation is available [B. L. Hu, J. P. Paz, and Y. Zhang, Phys. Rev. D 45, 2843 (1992)], the exact master equation with all its coefficients for this two harmonic oscillator model can be easily deduced from the known results of the single harmonic oscillator case. In the second part we give an influence functional treatment of this model and provide explicit expressions for the evolutionary operator of the reduced density matrix which are useful for the study of decoherence and disentanglement issues. We show three applications of this master equation: on the decoherence and disentanglement of two harmonic oscillators due to their interaction with a common environment under Markovian approximation, and a derivation of the uncertainty principle at finite temperature for a composite object, modeled by two interacting harmonic oscillators. The exact master equation for two, and its generalization to N, harmonic oscillators interacting with a general environment are expected to be useful for the analysis of quantum coherence, entanglement, fluctuations, and dissipation of mesoscopic objects toward the construction of a theoretical framework for macroscopic quantum phenomena.
The master equation for a linear open quantum system in a general environment is derived using a stochastic approach. This is an alternative derivation to that of Hu, Paz and Zhang, which was based on the direct computation of path integrals, or to that of Halliwell and Yu, based on the evolution of the Wigner function for a linear closed quantum system. We first show by using the influence functinal formalism that the reduced Wigner function for the open system coincides with a distribution function resulting from averaging both over the initial conditions and the stochastic source of a formal Langevin equation. The master equation for the reduced Wigner function can then be deduced as a Fokker-Planck equation obtained from the formal Langevin equation.