A method is described to solve the nonlinear Langevin equations arising from quadratic interactions in quantum mechanics. While the zeroth order linearization approximation to the operators is normally used, here, first and second order truncation perturbation schemes are proposed. These schemes employ higher-order system operators, and then approximate number operators with their corresponding mean boson numbers only where needed. Spectral densities of higher-order operators are derived, and an expression for the second-order correlation function at zero time-delay has been found, which reveals that the cavity photon occupation of an ideal laser at threshold reaches 6 − 2 , in good agreement with extensive numerical calculations. As further applications, analysis of the quantum anharmonic oscillator, calculation of Q-functions, analysis of quantum limited amplifiers, and nondemoliton measurements are provided.
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.
The quantum mechanical Langevin equation for a damped harmonic oscillator with operator-valued-noise, da t /dt=(-i omega - gamma )A t + lambda w(t), is reinvestigated. By assuming the canonical commutation relation at initial time, (A 0 , A 0 *)=1, and a Hermitian operator for the Gaussian quantum noise w(t), it is shown that Streater's procedure to satisfy (A t ,A t *)=1 for all later time yields an improved result: the noise operator w(t) can be determined universally and normalised such that, by choosing lambda 2 =2 gamma ( beta h(cross) omega ) -1 , the (symmetrised) power spectrum of w(t) is equal to the physically significant form 1 / 2 beta h(cross) omega coth( 1 / 2 beta h(cross) omega ), where beta is the inverse temperature of the heat bath, thus ensuring the beta -KMS correlation functions for any pair of operators.
We derive a
c-number generalized Langevin equation (GLE) describing the evolution of the expectation values ${\ensuremath{\langle}{x}_{i}\ensuremath{\rangle}}_{t}$ of the atomic position operators
xi of an open system. The latter is coupled linearly to a harmonic bath kept at a fixed temperature. The equations of motion contain a non-Markovian friction term with the classical kernel [L. Kantorovich, Phys. Rev. B 78, 094304 (2008)] and a zero mean non-Gaussian random force with correlation functions that depend on the initial preparation of the open system. We used a density operator formalism without assuming that initially the combined system was decoupled. The only approximation made in deriving quantum GLE consists of assuming that the Hamiltonian of the open system at time
t can be expanded up to the second order with respect to operators of atomic displacements ${u}_{i}={x}_{i}\ensuremath{-}{\ensuremath{\langle}{x}_{i}\ensuremath{\rangle}}_{t}$ in the open system around their exact atomic positions ${\ensuremath{\langle}{x}_{i}\ensuremath{\rangle}}_{t}$ (the ``harmonization'' approximation). The noise is introduced to ensure that sampling many quantum GLE trajectories yields exactly the average one. An explicit expression for the pair correlation function of the noise, consistent with the classical limit, is also proposed. Unlike the usually considered quantum operator GLE, the proposed
c-number quantum GLE can be used in direct molecular dynamic simulations of open systems under general equilibrium or nonequilibrium conditions.