Abstract
Quantum harmonic oscillators linearly coupled through coordinates and momenta, represented by the Hamiltonian H=Σ²_i=1( p²_i/2 m_i + m_i ω²_i/2 x²_i) +H_int , where the interaction of two oscillators H_int = i k₁ x₁ p ₂+ i k₂ x₂ p₁ + k₃ x₁ x₂-k₄ p₁ p₂, found in many applications of quantum optics, nonlinear physics, molecular chemistry and biophysics. Despite this, there is currently no general solution to the Schr\"{o}dinger equation for such a system. This is especially relevant for quantum entanglement of such a system in quantum optics applications. Here this problem is solved and it is shown that quantum entanglement depends on only one coefficient R ∈ (0,1), which includes all the parameters of the system under consideration. It has been shown that quantum entanglement can be very large at certain values of this coefficient. The results obtained have a fairly simple analytical form, which facilitates analysis.