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.
Abstract We propose a Langevin equation to describe the quantum Brownian motion of bounded particles based on a distinctive formulation concerning both the fluctuation and dissipation forces. The fluctuation force is similar to that employed in the classical case. It is a white noise with a variance proportional to the temperature. The dissipation force is not restricted to be proportional to the velocity and is determined in a way as to guarantee that the stationary state is given by a density operator of the Gibbs canonical type. To this end we derived an equation that gives the time evolution of the density operator, which turns out to be a quantum Fokker–Planck–Kramers equation. The approach is applied to the harmonic oscillator in which case the dissipation force is found to be non Hermitian and proportional to the velocity and position.