Linearly Coupled Quantum Harmonic Oscillators and Their Quantum Entanglement

In many applications of quantum optics, nonlinear physics, molecular chemistry and biophysics, one can encounter models in which the coupled quantum harmonic oscillator provides an explanation for many physical phenomena and effects. In general, these are harmonic oscillators coupled via coordinates and momenta, which can be represented as H ^ = ∑ i = 1 2 p ^ i 2 2 m i + m i ω i 2 2 x i 2 + H ^ i n t , where the interaction of two oscillators H ^ i n t = i k 1 x 1 p ^ 2 + i k 2 x 2 p ^ 1 + k 3 x 1 x 2 − k 4 p ^ 1 p ^ 2 . Despite the importance of this system, there is currently no general solution to the Schrödinger equation that takes into account arbitrary initial states of the oscillators. Here, this problem is solved in analytical form, and it is shown that the probability of finding the system in any states and quantum entanglement depends only on one coefficient R ∈ ( 0 , 1 ) for the initial factorizable Fock states of the oscillator and depends on two parameters R ∈ ( 0 , 1 ) and ϕ for arbitrary initial states. These two parameters R ∈ ( 0 , 1 ) and ϕ include the entire set of variables of the system under consideration.