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.
This article presents a full operator analytical method for studying the quadratic nonlinear interactions in quantum optomechanics. The method is based on the application of higher-order operators, using a six-dimensional basis of second order operators which constitute an exactly closed commutators. We consider both types of standard position-field and the recently predicted non-standard momentum-field quadratic interactions, which is significant when the ratio of mechanical frequency to optical frequency is not negligible. This unexplored regime of large mechanical frequency can be investigated in few platforms including the superconducting electromechanics and simulating quantum cavity electrodynamic circuits. It has been shown that the existence of non-standard quadratic interaction could be observable under appropriate conditions.
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 success of quantum physics in description of various physical interaction phenomena relies primarily on the accuracy of analytical methods used. In quantum mechanics, many of such interactions such as those found in quantum optomechanics and quantum computing have a highly nonlinear nature, which makes their analysis extraordinarily difficult using classical schemes. Typically, modern quantum systems of interest nowadays come with four basic properties: (i) quantumness, (ii) openness, (iii) randomness, and (iv) nonlinearity. The newly introduced method of higher-order operators targets analytical solutions to such systems, and while providing at least mathematically approximate expressions with improved accuracy over the fully linearized schemes, some cases admit exact solutions. Many different applications of this method in quantum and classically nonlinear systems are demonstrated throughout. This review is purposed to provide the reader with ease of access to this recent and well-established operator algebra, while going over a moderate amount of literature review. The reader with basic knowledge of quantum mechanics and quantum noise theory should be able to start using this scheme to his or her own problem of interest.
With the help of second-order nonlinear interactions and hence Hamiltonian, we construct the equations of motion corresponding to pump and signal (second harmonic) fields in a second harmonic generation. The corresponding coupled differential equations involving the non-commuting field operators are not solvable in closed analytical forms. With the help of a perturbation method, we obtain analytical solutions of these field operators up to cubic orders in the interaction constant. In an appropriate limit (truncating the solution up to the second order in the interaction constant), these solutions lead to the existing solutions under the short-time approximation. The present analytical solutions are exploited to investigate the squeezing and the antibunching of photons of the input coherent light coupled to the said second-order nonlinear medium. We report the squeezing and antibunching effects of the pump field even for solutions up to the second order and hence the present results are consistent with earlier results. However, for the second harmonic mode, we report the squeezing involving the leading cubic term in the interaction constant.
We generalise the quantum mechanical toy model for gravitationally induced decoherence presented in Xu, Blencowe (2022) and Domi et al. (2024). In contrast to earlier formulations, in which the Hamiltonian of the system of interest is linearly coupled to the position operators of the oscillators in the environment, we consider an interaction formulated in terms of Weyl elements of the environment's position operators. This extension is motivated by polymer quantum mechanics, in which Weyl elements are fundamental operators, as well as by the possibility of generating non-linear interactions through suitable truncations of the exponential Weyl elements. Here we focus on a sinus-like coupling that is still quantised using the Schrödinger representation and, in the limit of a small Weyl parameter, reproduces the conventional linear interaction. To derive the corresponding master equation, we developed two complementary methods for the analytical calculation of the environmental correlation functions. The first utilises Wick's theorem for thermal expectation values in conjunction with annihilation and creation operators, while the second is based on the short-time Fourier transform and completely avoids the use of annihilation and creation operators, making it more readily transferable and generalisable to a polymer quantisation. Both approaches yield identical results. We further generalise the spectral density required for the exponential coupling structure. A numerical analysis shows that the environmental correlation functions decay rapidly with time, which supports the validity of the Markov approximation. Using a Taylor expansion in the Weyl parameter, we show that the first-order term reproduces the decoherence model of Xu, Blencowe (2022) and Domi et al. (2024). Finally, we derive the solution to the renormalised master equation.
The out-of-time-order correlator (OTOC) of simple harmonic oscillator with extra anharmonic (quartic) interaction are calculated by the second quantization method. We obtain the analytic formulas of spectrum, Fock space states and matrix elements of coordinate to the second order of anharmonic interaction. These relations clearly reveal the property that the correction of the interaction is proportional to the quantum number of the energy level, and shows the enhancement for some physical quantities. The analytic results explicitly show that the OTOC is raising in the quadratic power law at early times. Then, we use the formulas to do numerical summation to calculate the OTOC, which shows that, at late times, while the first-order perturbative OTOC is oscillating as that in a simple harmonic oscillator, the second-order perturbative OTOC saturates to a constant value. We compare it with
2⟨x2⟩T⟨p2⟩T, which is associated with quantum chaotic behavior in systems that exhibit chaos, and discuss the validity of the second-order perturbation.