A projection operator method in Heisenberg picture is developed. Corresponding to the equations for the reduced density matrix, there are two types of Langevin-like equations; time convolution (TC) and time convolutionless (TCL) equations. When the whole system is composed of relevant system and irrelevant system, unnecessary information is eliminated by projection operator. Explicit expansion formulae are obtained for the both types of equations. Relaxation operators and fluctuating forces are expressed in terms of certain cumulants. The method is illustrated by simple examples.
We systematically derive the quantum generalized nonlinear Langevin equation using Morozov's projection operator method. This approach extends the linear Mori-Zwanzig projection operator technique, allowing for the inclusion of nonlinear interactions among macroscopic modes. Additionally, we obtain the quantum generalized Fokker-Planck equation within the Heisenberg picture, which is consistent with Morozov's original formulation. These equations are fundamentally significant in nonequilibrium statistical physics, particularly in scenarios characterized by enhanced fluctuations, such as anomalous transport phenomena near critical points. The quantum nature of the derived generalized Langevin and Fokker-Planck equations is anticipated to provide a more detailed description than their classical equivalents. Specifically, the noise kernel in the quantum generalized Langevin equation is multiplicative, which broadens the applicability beyond Gaussian approximations. Given specific interactions, these equations are expected to be instrumental in investigating critical transport phenomena.
Building upon the work of Hu, Paz, and Zhang [1,2] on open quantum systems we consider the quantum Brownian motion (QBM) model with one oscillator (position variable
x) as the system, {\it nonlinearly} coupled to an environment of
N harmonic oscillators (with mass
mn, natural frequency
ωn, position
qn and momentum
pn variables) in the form
∑n(vn1(x)qnk+vn2(x)pnl) where
k,l are integers (the present work only considers the
k=l=2 cases). The vertex functions $v_{n1}, v_{n2} $ are of the form
vn1=λCn1f(x),vn2(x)=−λCn2mn−2ωn−2f(x) where
Cn1,2 are the coupling constants with the
nth oscillator,
f(x) is any arbitrary function of
x, and
λ is a dimensionless constant. Employing the closed-time-path formalism the influence action
SIF is calculated using a perturbative expansion in
λ. It is possible to identify the terms in
SIF quadratic or higher in
Δ(s)≡f(x+(s))−f(x−(s)) to constitute the noise kernel, while terms linear in
Δ to that of the dissipation kernel. The non-Gaussian noise kernel gives rise to non-zero three-point correlation function of the corresponding stochastic force. The pathway presented here should be useful for the exploration of \textit{non-Gaussian properties of systems nonlinearly coupled with their environments}; examples in early universe cosmology and in quantum optomechanics (QOM) are mentioned. A modified fluctuation-dissipation relation (FDR) is also established, which ensures the consistency of the model and the accuracy of results even at higher perturbative orders. Another result of significance is the derivation of a nonlinear Langevin equation which is expected to be useful for many open quantum system applications.
Overview Projection operator technique is simple in purpose. It allows us to separate the dynamical evolution of a quantity into two parts: known and unknown, the latter to be evaluated. However, this technique serves multiple purposes in TDSM. It is used to derive a reduced equation of motion for the relaxation of a dynamical variable. It is also used to derive expression for an experimentally observed TCF in terms of theoretically known TCFs. This technique was developed by Zwanzig, who was inspired by Kubo’s approach of “reduction in degrees of freedom” in an interacting many-body system. In the projection operator method, we divide the time evolution of a dynamical variable A ( t ) into two parts: (i) one through a relevant subspace that is closely related to A and varies on a similar time scale, and (ii) an orthogonal subspace, which could be fast and therefore can be averaged out to obtain a time-dependent friction kernel. The idea then is that the first contribution is known while the second may be described accurately by projecting it onto the relevant subspace [that is, (i)]. In the derivation of the reduced equation of motion, one starts with the Liouville equation and projects out irrelevant degrees of freedom to obtain, after specific algebraic manipulations, an equation of motion of the dynamical variable. The final relaxation equation has the form of the Langevin equation with microscopic expressions for the terms involved (frequency term, the friction, and the random force term). This final equation of motion is often termed relaxation equation or memory function equation. It finds wide use in the study of time-dependent statistical mechanics. In the present chapter we derive the memory function equation and discuss a few applications.
In the previous chapter, we developed the equation of motion for a system as it evolved under the influence of an unobserved (reservoir) system. We used the density matrix approach and worked in the interaction picture. In this chapter, we consider the same problem of the system-reservoir interaction using a quantum operator approach. We again eliminate the reservoir variables. The resulting equations for the system operators include, in addition to the damping terms, the noise operators which produce fluctuations. These equations have the form of classical Langevin equations, which describe, for example, the Brownian motion of a particle suspended in a liquid. The Heisenberg–Langevin approach discussed in this chapter is particularly suitable for the calculation of two-time correlation functions of the system operator as is, for example, required for the determination of the natural linewidth of a laser. We first consider the damping of the harmonic oscillator by an interaction with a reservoir consisting of many other simple harmonic oscillators. This system describes, for example, the damping of a single-mode field inside a cavity with lossy mirrors. The reservoir, in this case, consists of a large number of phonon-like modes in the mirrors. We also consider the decay of the field due to its interaction with an atomic reservoir. An interesting application of the theory of the system–reservoir interaction is the evolution of an atom inside a damped cavity. It is shown that the spontaneous transition rate of the atom can be substantially enhanced if it is placed in a resonant cavity.