Exact solution of the isotropic and anisotropic Hamiltonian of two coupled harmonic oscillators

Abstract

We study the isotropic and anisotropic Hamiltonian of two coupled harmonic oscillators from an algebraic approach of the $SU(1,1)$ and $SU(2)$ groups. In order to obtain the energy spectrum and eigenfunctions of this problem, we write its Hamiltonian in terms of the boson generators of the $SU(1,1)$ and $SU(2)$ groups. We use the one boson and two boson realizations of the $su(1,1)$ Lie algebra, and the one boson realization of the $su(2)$ Lie algebra to apply three tilting transformations to diagonalize the original Hamiltonian. These transformations let us to obtain the exact solutions of the isotropic and the anisotropic cases, from which the particular expected results are obtained for the cases where the coupling is neglected.

1 Introduction

Of all the problems that have contributed to our understanding of the universe around us, perhaps the simple harmonic oscillator is the most important. It has been a fundamental piece in the development of all branches of physics and its applications have extended to many different branches of science [1]. The mathematical simplicity of harmonic oscillators is contrasted by the enormous variety of methods and techniques that can be used to describe them. However, harmonic oscillators are not always disconnected from each other, and this is where the importance of studying coupled harmonic oscillators lies.

The study of coupled oscillators has found its application in nonlinear systems, as can be seen in Refs. [2, 3, 4, 5, 6, 7, 8, 9, 10]. One of the most studied characteristics of coupled oscillators is their quantum entanglement [11, 12, 13, 14], which has been applied to biophysics, quantum cryptography, quantum-coding, quantum computing algorithms, quantum state teleportation, among others [15, 16, 17, 18, 19, 20, 21]. In particular, in Ref. [22] the quantum entanglement of oscillators with position-position coupling is studied and the solution to the nonstationary Schrödinger equation is found.

The aim of this work is to study and exactly solve isotropic and anisotropic Hamiltonian of two coupled harmonic oscillators from an algebraic approach. We obtain its eigenfunctions and the energy spectrum by using similarity transformations in terms of the $su(1,1)$ and $su(2)$ Lie algebras.

This work is organized as it follows. In Sec. $2$ we write the Hamiltonian of the two isotropic coupled oscillators in terms of the bosonic operators $a$ and $b$. Then, we introduce Schwinger realizations in terms of one and two bosons of the $su(1,1)$ and $su(2)$ Lie algebras and apply three tilting transformations in terms of these realizations to the Hamiltonian to diagonalize it. These similarity transformations let us exactly solve the problem obtaining the eigenfunctions and the energy spectrum. In. Section $3$ we apply the three tilting transformations to the anisotropic case of our Hamiltonian in order to obtain the most general solution to our problem. Finally, we give some concluding remarks.

2 Hamiltonian of two coupled harmonic oscillators: The isotropic case

A Hamiltonian which describes the behavior of a system of two coupled harmonic oscillators with identical free frequencies can be seen as a Hamiltonian which coupling is of the position-position type. It is defined as

$$H_{xy}=\frac{1}{2m}\left(p_{x}^{2}+p_{y}^{2}\right)+\frac{m\omega^{2}}{2}\left(x^{2}+y^{2}\right)+2\kappa m\omega xy,$$ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 end_ARG ( italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) + 2 italic_κかっぱ italic_m italic_ωおめが italic_x italic_y , (1)

where $\kappa$κかっぱ is a constant that describes the strength of the coupling. This Hamiltonian can be written in terms of the bosonic operators as

$$H_{xy}=\hbar\omega\left(a^{{\dagger}}a+\frac{1}{2}\right)+\hbar\omega\left(b^{{\dagger}}b+\frac{1}{2}\right)+\hbar\kappa(a^{{\dagger}}+a)(b^{{\dagger}}+b).$$ωおめが ( italic_a start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_a + divide start_ARG 1 end_ARG start_ARG 2 end_ARG ) + roman_ℏ italic_ωおめが ( italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_b + divide start_ARG 1 end_ARG start_ARG 2 end_ARG ) + roman_ℏ italic_κかっぱ ( italic_a start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT + italic_a ) ( italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT + italic_b ) . (2)

In order to study the system of two coupled harmonic oscillators with identical free frequencies we propose the following more general Hamiltonian (with $\hbar=1$)

$$H=\omega(a^{\dagger}a+b^{\dagger}b+1)+\lambda e^{-i\psi}(a^{\dagger}b^{\dagger}+a^{\dagger}b)+\lambda e^{i\psi}(b^{\dagger}a+ba).$$ωおめが ( italic_a start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_a + italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_b + 1 ) + italic_λらむだ italic_e start_POSTSUPERSCRIPT - italic_i italic_ψぷさい end_POSTSUPERSCRIPT ( italic_a start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT + italic_a start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_b ) + italic_λらむだ italic_e start_POSTSUPERSCRIPT italic_i italic_ψぷさい end_POSTSUPERSCRIPT ( italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_a + italic_b italic_a ) . (3)

Using the two-bosons Schwinger realizations of the $SU(2)$ and $SU(1,1)$ Lie algebras of equations (A.13) and (A.14), we can write this Hamiltonian as

$$H=2\omega K_{0}+\lambda e^{-i\psi}(K_{+}+J_{+})+\lambda e^{i\psi}(K_{-}+J_{-}).$$ωおめが italic_K start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_λらむだ italic_e start_POSTSUPERSCRIPT - italic_i italic_ψぷさい end_POSTSUPERSCRIPT ( italic_K start_POSTSUBSCRIPT + end_POSTSUBSCRIPT + italic_J start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) + italic_λらむだ italic_e start_POSTSUPERSCRIPT italic_i italic_ψぷさい end_POSTSUPERSCRIPT ( italic_K start_POSTSUBSCRIPT - end_POSTSUBSCRIPT + italic_J start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) . (4)

From now on we will diagonalize the Hamiltonian $H$ by applying three consecutive tilting transformations, expressed in terms of the $SU(1,1)$ and $SU(2)$ displacement operators (see Appendix). Then, we first apply a tilting transformation to the stationary Schrödinger equation $H\mathrm{<ruby><rb>\Psi </rb><rt>ぷさい</rt></ruby>}=E\mathrm{<ruby><rb>\Psi </rb><rt>ぷさい</rt></ruby>}$Ψぷさい𝐸ΨぷさいH\Psi=E\Psiitalic_H roman_Ψぷさい = italic_E roman_Ψぷさい as it follows [23, 24, 25]

$$D^{†}(\mathrm{<ruby><rb>\xi </rb><rt>くしー</rt></ruby>})HD(\mathrm{<ruby><rb>\xi </rb><rt>くしー</rt></ruby>})D^{†}(\mathrm{<ruby><rb>\xi </rb><rt>くしー</rt></ruby>})\mathrm{<ruby><rb>\Psi </rb><rt>ぷさい</rt></ruby>}=E{D}^{†}(\mathrm{<ruby><rb>\xi </rb><rt>くしー</rt></ruby>})\mathrm{<ruby><rb>\Psi </rb><rt>ぷさい</rt></ruby>}.$$Ψぷさい𝐸superscript𝐷†𝜉ΨぷさいD^{\dagger}(\xi)HD(\xi)D^{\dagger}(\xi)\Psi=ED^{\dagger}(\xi)\Psi.italic_D start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_ξくしー ) italic_H italic_D ( italic_ξくしー ) italic_D start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_ξくしー ) roman_Ψぷさい = italic_E italic_D start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_ξくしー ) roman_Ψぷさい . (5)

Here, $D(\xi)$ξくしー ) is the $SU(1,1)$ displacement operator defined as $D(\xi)=\exp(\xi K_{+}-\xi^{*}K_{-})$ξくしー ) = roman_exp ( start_ARG italic_ξくしー italic_K start_POSTSUBSCRIPT + end_POSTSUBSCRIPT - italic_ξくしー start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT italic_K start_POSTSUBSCRIPT - end_POSTSUBSCRIPT end_ARG ), with $\xi=-{}^{1}\!/_{2}\tau e^{-i\phi_{\xi}}$ξくしー = - start_FLOATSUPERSCRIPT 1 end_FLOATSUPERSCRIPT / start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_τたう italic_e start_POSTSUPERSCRIPT - italic_i italic_ϕ start_POSTSUBSCRIPT italic_ξくしー end_POSTSUBSCRIPT end_POSTSUPERSCRIPT. Thus, if we define the tilted Hamiltonian $H^{\prime}$ as $H^{\prime}=D^{\dagger}(\xi)HD(\xi)$ξくしー ) italic_H italic_D ( italic_ξくしー ) and its wavefunctions ${\mathrm{<ruby><rb>\Psi </rb><rt>ぷさい</rt></ruby>}}^{\prime }$Ψぷさい′\Psi^{\prime}roman_Ψぷさい start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT as ${\mathrm{<ruby><rb>\Psi </rb><rt>ぷさい</rt></ruby>}}^{\prime }=D^{†}(\mathrm{<ruby><rb>\xi </rb><rt>くしー</rt></ruby>})\mathrm{<ruby><rb>\Psi </rb><rt>ぷさい</rt></ruby>}$Ψぷさい′superscript𝐷†𝜉Ψぷさい\Psi^{\prime}=D^{\dagger}(\xi)\Psiroman_Ψぷさい start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = italic_D start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_ξくしー ) roman_Ψぷさい, from the similarity transformations (B.2) and (B.3) of Appendix B we obtain

$$H^{\prime} = 2\omega\left[(2\beta_{\xi}+1)K_{0}+\frac{\alpha_{\xi}\xi}{2\absolutevalue{\xi}}K_{+}+\frac{\alpha_{\xi}\xi^{*}}{2\absolutevalue{\xi}}K_{-}\right]$$ ωおめが [ ( 2 italic_βべーた start_POSTSUBSCRIPT italic_ξくしー end_POSTSUBSCRIPT + 1 ) italic_K start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + divide start_ARG italic_αあるふぁ start_POSTSUBSCRIPT italic_ξくしー end_POSTSUBSCRIPT italic_ξくしー end_ARG start_ARG 2 | start_ARG italic_ξくしー end_ARG | end_ARG italic_K start_POSTSUBSCRIPT + end_POSTSUBSCRIPT + divide start_ARG italic_αあるふぁ start_POSTSUBSCRIPT italic_ξくしー end_POSTSUBSCRIPT italic_ξくしー start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT end_ARG start_ARG 2 | start_ARG italic_ξくしー end_ARG | end_ARG italic_K start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ] (6)

$$+\lambda e^{-i\psi}\left[\left[\frac{\xi^{*}}{\absolutevalue{\xi}}\alpha_{\xi}K_{0}+\beta_{\xi}\left(K_{+}+\frac{\xi^{*}}{\xi}K_{-}\right)+K_{+}\right]+\left[\frac{\xi^{*}}{\absolutevalue{\xi}}\alpha_{\xi}K_{-}^{(b)}+\frac{\xi}{\absolutevalue{\xi}}\alpha_{\xi}K_{+}^{(a)}+(2\beta_{\xi}+1)J_{+}\right]\right]$$λらむだ italic_e start_POSTSUPERSCRIPT - italic_i italic_ψぷさい end_POSTSUPERSCRIPT [ [ divide start_ARG italic_ξくしー start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT end_ARG start_ARG | start_ARG italic_ξくしー end_ARG | end_ARG italic_αあるふぁ start_POSTSUBSCRIPT italic_ξくしー end_POSTSUBSCRIPT italic_K start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_βべーた start_POSTSUBSCRIPT italic_ξくしー end_POSTSUBSCRIPT ( italic_K start_POSTSUBSCRIPT + end_POSTSUBSCRIPT + divide start_ARG italic_ξくしー start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT end_ARG start_ARG italic_ξくしー end_ARG italic_K start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) + italic_K start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ] + [ divide start_ARG italic_ξくしー start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT end_ARG start_ARG | start_ARG italic_ξくしー end_ARG | end_ARG italic_αあるふぁ start_POSTSUBSCRIPT italic_ξくしー end_POSTSUBSCRIPT italic_K start_POSTSUBSCRIPT - end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT + divide start_ARG italic_ξくしー end_ARG start_ARG | start_ARG italic_ξくしー end_ARG | end_ARG italic_αあるふぁ start_POSTSUBSCRIPT italic_ξくしー end_POSTSUBSCRIPT italic_K start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT + ( 2 italic_βべーた start_POSTSUBSCRIPT italic_ξくしー end_POSTSUBSCRIPT + 1 ) italic_J start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ] ]

$$+\lambda e^{i\psi}\left[\left[\frac{\xi}{\absolutevalue{\xi}}\alpha_{\xi}K_{0}+\beta_{\xi}\left(K_{-}+\frac{\xi}{\xi^{*}}K_{+}\right)+K_{-}\right]+\left[\frac{\xi^{*}}{\absolutevalue{\xi}}\alpha_{\xi}K_{-}^{(a)}+\frac{\xi}{\absolutevalue{\xi}}\alpha_{\xi}K_{+}^{(b)}+(2\beta_{\xi}+1)J_{-}\right]\right],$$λらむだ italic_e start_POSTSUPERSCRIPT italic_i italic_ψぷさい end_POSTSUPERSCRIPT [ [ divide start_ARG italic_ξくしー end_ARG start_ARG | start_ARG italic_ξくしー end_ARG | end_ARG italic_αあるふぁ start_POSTSUBSCRIPT italic_ξくしー end_POSTSUBSCRIPT italic_K start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_βべーた start_POSTSUBSCRIPT italic_ξくしー end_POSTSUBSCRIPT ( italic_K start_POSTSUBSCRIPT - end_POSTSUBSCRIPT + divide start_ARG italic_ξくしー end_ARG start_ARG italic_ξくしー start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT end_ARG italic_K start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) + italic_K start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ] + [ divide start_ARG italic_ξくしー start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT end_ARG start_ARG | start_ARG italic_ξくしー end_ARG | end_ARG italic_αあるふぁ start_POSTSUBSCRIPT italic_ξくしー end_POSTSUBSCRIPT italic_K start_POSTSUBSCRIPT - end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT + divide start_ARG italic_ξくしー end_ARG start_ARG | start_ARG italic_ξくしー end_ARG | end_ARG italic_αあるふぁ start_POSTSUBSCRIPT italic_ξくしー end_POSTSUBSCRIPT italic_K start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT + ( 2 italic_βべーた start_POSTSUBSCRIPT italic_ξくしー end_POSTSUBSCRIPT + 1 ) italic_J start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ] ] ,

where $\alpha_{\xi}=\sinh{(2\absolutevalue{\xi})}$αあるふぁ start_POSTSUBSCRIPT italic_ξくしー end_POSTSUBSCRIPT = roman_sinh ( 2 | start_ARG italic_ξくしー end_ARG | ) and $\beta_{\xi}={}^{1}\!/_{2}[\cosh{(2\absolutevalue{\xi})}-1]$βべーた start_POSTSUBSCRIPT italic_ξくしー end_POSTSUBSCRIPT = start_FLOATSUPERSCRIPT 1 end_FLOATSUPERSCRIPT / start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT [ roman_cosh ( start_ARG ( 2 | start_ARG italic_ξくしー end_ARG | ) end_ARG ) - 1 ]. Notice that in this expression we have introduced $\{K_{\pm}^{a},K_{\pm}^{b}\}$, which are the one-boson ladder operators for the $SU(1,1)$ realization (see Eq. (A.15) of Appendix A). This tilted Hamiltonian $H^{\prime}$ can be written as

$$H^{\prime} = \left[2\omega(2\beta_{\xi}+1)+\frac{\xi^{*}}{\absolutevalue{\xi}}\alpha_{\xi}\lambda e^{-i\psi}+\frac{\xi}{\absolutevalue{\xi}}\alpha_{\xi}\lambda e^{i\psi}\right]K_{0}$$ ωおめが ( 2 italic_βべーた start_POSTSUBSCRIPT italic_ξくしー end_POSTSUBSCRIPT + 1 ) + divide start_ARG italic_ξくしー start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT end_ARG start_ARG | start_ARG italic_ξくしー end_ARG | end_ARG italic_αあるふぁ start_POSTSUBSCRIPT italic_ξくしー end_POSTSUBSCRIPT italic_λらむだ italic_e start_POSTSUPERSCRIPT - italic_i italic_ψぷさい end_POSTSUPERSCRIPT + divide start_ARG italic_ξくしー end_ARG start_ARG | start_ARG italic_ξくしー end_ARG | end_ARG italic_αあるふぁ start_POSTSUBSCRIPT italic_ξくしー end_POSTSUBSCRIPT italic_λらむだ italic_e start_POSTSUPERSCRIPT italic_i italic_ψぷさい end_POSTSUPERSCRIPT ] italic_K start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT (7)

$$+\left[\frac{\omega\alpha_{\xi}\xi}{\absolutevalue{\xi}}+(\beta_{\xi}+1)\lambda e^{-i\psi}+\frac{\beta_{\xi}\xi}{\xi^{*}}\lambda e^{i\psi}\right]K_{+}+\left[\frac{\omega\alpha_{\xi}\xi^{*}}{\absolutevalue{\xi}}+\frac{\beta_{\xi}\xi^{*}}{\xi}\lambda e^{-i\psi}+(\beta_{\xi}+1)\lambda e^{i\psi}\right]K_{-}$$ωおめが italic_αあるふぁ start_POSTSUBSCRIPT italic_ξくしー end_POSTSUBSCRIPT italic_ξくしー end_ARG start_ARG | start_ARG italic_ξくしー end_ARG | end_ARG + ( italic_βべーた start_POSTSUBSCRIPT italic_ξくしー end_POSTSUBSCRIPT + 1 ) italic_λらむだ italic_e start_POSTSUPERSCRIPT - italic_i italic_ψぷさい end_POSTSUPERSCRIPT + divide start_ARG italic_βべーた start_POSTSUBSCRIPT italic_ξくしー end_POSTSUBSCRIPT italic_ξくしー end_ARG start_ARG italic_ξくしー start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT end_ARG italic_λらむだ italic_e start_POSTSUPERSCRIPT italic_i italic_ψぷさい end_POSTSUPERSCRIPT ] italic_K start_POSTSUBSCRIPT + end_POSTSUBSCRIPT + [ divide start_ARG italic_ωおめが italic_αあるふぁ start_POSTSUBSCRIPT italic_ξくしー end_POSTSUBSCRIPT italic_ξくしー start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT end_ARG start_ARG | start_ARG italic_ξくしー end_ARG | end_ARG + divide start_ARG italic_βべーた start_POSTSUBSCRIPT italic_ξくしー end_POSTSUBSCRIPT italic_ξくしー start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT end_ARG start_ARG italic_ξくしー end_ARG italic_λらむだ italic_e start_POSTSUPERSCRIPT - italic_i italic_ψぷさい end_POSTSUPERSCRIPT + ( italic_βべーた start_POSTSUBSCRIPT italic_ξくしー end_POSTSUBSCRIPT + 1 ) italic_λらむだ italic_e start_POSTSUPERSCRIPT italic_i italic_ψぷさい end_POSTSUPERSCRIPT ] italic_K start_POSTSUBSCRIPT - end_POSTSUBSCRIPT

$$+\frac{\xi^{*}}{\absolutevalue{\xi}}\alpha_{\xi}\lambda(K_{-}^{(b)}e^{-i\psi}+K_{-}^{(a)}e^{i\psi})+\frac{\xi}{\absolutevalue{\xi}}\alpha_{\xi}\lambda(K_{+}^{(a)}e^{-i\psi}+K_{+}^{(b)}e^{i\psi})+(2\beta_{\xi}+1)\lambda(J_{+}e^{-i\psi}+J_{-}e^{i\psi}).$$ξくしー start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT end_ARG start_ARG | start_ARG italic_ξくしー end_ARG | end_ARG italic_αあるふぁ start_POSTSUBSCRIPT italic_ξくしー end_POSTSUBSCRIPT italic_λらむだ ( italic_K start_POSTSUBSCRIPT - end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT - italic_i italic_ψぷさい end_POSTSUPERSCRIPT + italic_K start_POSTSUBSCRIPT - end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT italic_i italic_ψぷさい end_POSTSUPERSCRIPT ) + divide start_ARG italic_ξくしー end_ARG start_ARG | start_ARG italic_ξくしー end_ARG | end_ARG italic_αあるふぁ start_POSTSUBSCRIPT italic_ξくしー end_POSTSUBSCRIPT italic_λらむだ ( italic_K start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT - italic_i italic_ψぷさい end_POSTSUPERSCRIPT + italic_K start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT italic_i italic_ψぷさい end_POSTSUPERSCRIPT ) + ( 2 italic_βべーた start_POSTSUBSCRIPT italic_ξくしー end_POSTSUBSCRIPT + 1 ) italic_λらむだ ( italic_J start_POSTSUBSCRIPT + end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT - italic_i italic_ψぷさい end_POSTSUPERSCRIPT + italic_J start_POSTSUBSCRIPT - end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_i italic_ψぷさい end_POSTSUPERSCRIPT ) .

If we choose the values of the coherent state parameters $\tau$τたう and $\phi_{\xi}$ξくしー end_POSTSUBSCRIPT such that the coefficients of the operators $K_{\pm}$ vanish, we find

$$\tau=\tanh^{-1}{\left(\frac{\lambda}{\omega}\right)},\hskip 8.5359pt\phi_{\xi}=\psi.$$τたう = roman_tanh start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ( divide start_ARG italic_λらむだ end_ARG start_ARG italic_ωおめが end_ARG ) , italic_ϕ start_POSTSUBSCRIPT italic_ξくしー end_POSTSUBSCRIPT = italic_ψぷさい . (8)

This let us obtain the following reduced expression of the Hamiltonian $H^{\prime}$

$$H^{\prime} = 2[\omega\cosh{(\tau)}-\lambda\sinh{(\tau)}]K_{0}+\lambda\cosh{(\tau)}(J_{+}e^{-i\psi}+J_{-}e^{i\psi})$$ ωおめが roman_cosh ( start_ARG ( italic_τたう ) end_ARG ) - italic_λらむだ roman_sinh ( italic_τたう ) ] italic_K start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_λらむだ roman_cosh ( start_ARG ( italic_τたう ) end_ARG ) ( italic_J start_POSTSUBSCRIPT + end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT - italic_i italic_ψぷさい end_POSTSUPERSCRIPT + italic_J start_POSTSUBSCRIPT - end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_i italic_ψぷさい end_POSTSUPERSCRIPT ) (9)

$$-\lambda\sinh{(\tau)}(K_{+}^{(b)}+K_{-}^{(b)})-\lambda\sinh{(\tau)}(K_{+}^{(a)}e^{-2i\psi}+K_{-}^{(a)}e^{2i\psi}).$$λらむだ roman_sinh ( italic_τたう ) ( italic_K start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT + italic_K start_POSTSUBSCRIPT - end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT ) - italic_λらむだ roman_sinh ( italic_τたう ) ( italic_K start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT - 2 italic_i italic_ψぷさい end_POSTSUPERSCRIPT + italic_K start_POSTSUBSCRIPT - end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT 2 italic_i italic_ψぷさい end_POSTSUPERSCRIPT ) .

Using the expressions $\cosh{(\tau)}=\frac{\omega}{\sqrt{\omega^{2}-\lambda^{2}}}$τたう ) end_ARG ) = divide start_ARG italic_ωおめが end_ARG start_ARG square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG and $\sinh{(\tau)}=\frac{\lambda}{\sqrt{\omega^{2}-\lambda^{2}}}$τたう ) = divide start_ARG italic_λらむだ end_ARG start_ARG square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG we arrive to

$$H^{\prime} = 2\sqrt{\omega^{2}-\lambda^{2}}K_{0}+\frac{\omega\lambda}{\sqrt{\omega^{2}-\lambda^{2}}}(J_{+}e^{-i\psi}+J_{-}e^{i\psi})$$ ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG italic_K start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + divide start_ARG italic_ωおめが italic_λらむだ end_ARG start_ARG square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG ( italic_J start_POSTSUBSCRIPT + end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT - italic_i italic_ψぷさい end_POSTSUPERSCRIPT + italic_J start_POSTSUBSCRIPT - end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_i italic_ψぷさい end_POSTSUPERSCRIPT ) (10)

$$-\frac{\lambda^{2}}{\sqrt{\omega^{2}-\lambda^{2}}}(K_{+}^{(b)}+K_{-}^{(b)})-\frac{\lambda^{2}}{\sqrt{\omega^{2}-\lambda^{2}}}(K_{+}^{(a)}e^{-2i\psi}+K_{-}^{(a)}e^{2i\psi}).$$λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG ( italic_K start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT + italic_K start_POSTSUBSCRIPT - end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT ) - divide start_ARG italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG ( italic_K start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT - 2 italic_i italic_ψぷさい end_POSTSUPERSCRIPT + italic_K start_POSTSUBSCRIPT - end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT 2 italic_i italic_ψぷさい end_POSTSUPERSCRIPT ) .

Now, following a similar procedure, we can apply to $H^{\prime}$ a second tilting transformation. Thus, in terms of the $SU(2)$ displacement operator $D(\chi)=\exp(\chi J_{+}-\chi^{*}J_{-})$χかい ) = roman_exp ( start_ARG italic_χかい italic_J start_POSTSUBSCRIPT + end_POSTSUBSCRIPT - italic_χかい start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT italic_J start_POSTSUBSCRIPT - end_POSTSUBSCRIPT end_ARG ), we can define $H^{\prime\prime}$ and ${\mathrm{<ruby><rb>\Psi </rb><rt>ぷさい</rt></ruby>}}^{\prime \prime }$Ψぷさい′′\Psi^{\prime\prime}roman_Ψぷさい start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT as $H^{\prime\prime}=D^{\dagger}(\chi)H^{\prime}D(\chi)$χかい ) italic_H start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_D ( italic_χかい ) and ${\mathrm{<ruby><rb>\Psi </rb><rt>ぷさい</rt></ruby>}}^{\prime \prime }=D^{†}(\mathrm{<ruby><rb>\chi </rb><rt>かい</rt></ruby>}){\mathrm{<ruby><rb>\Psi </rb><rt>ぷさい</rt></ruby>}}^{\prime }$Ψぷさい′′superscript𝐷†𝜒superscriptΨぷさい′\Psi^{\prime\prime}=D^{\dagger}(\chi)\Psi^{\prime}roman_Ψぷさい start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT = italic_D start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_χかい ) roman_Ψぷさい start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT, respectively. Hence, the similarity transformations of Eqs. (B.4)-(B.6) of Appendix B let us to obtain

$$\begin{gathered}H^{\prime\prime}=2\sqrt{\omega^{2}-\lambda^{2}}K_{0}\\ +\frac{\omega\lambda}{\sqrt{\omega^{2}-\lambda^{2}}}\left[e^{-i\psi}\left[-\frac{\chi^{*}}{\absolutevalue{\chi}}\alpha_{\chi}J_{0}+\beta_{\chi}\left(J_{+}+\frac{\chi^{*}}{\chi}J_{-}\right)+J_{+}\right]+e^{i\psi}\left[-\frac{\chi}{\absolutevalue{\chi}}\alpha_{\chi}J_{0}+\beta_{\chi}\left(J_{-}+\frac{\chi}{\chi^{*}}J_{+}\right)+J_{-}\right]\right]\\ -\frac{\lambda^{2}}{\sqrt{\omega^{2}-\lambda^{2}}}\left[\left[\frac{K_{-}^{(b)}}{2}\left(\cos{(2\absolutevalue{\chi})}+1\right)-\frac{\chi^{*}}{2\absolutevalue{\chi}}\sin{(2\absolutevalue{\chi})}K_{-}-\frac{\chi^{*}}{2\chi}\left(\cos{(2\absolutevalue{\chi})}-1\right)K_{-}^{(a)}\right]\right.\\ \left.+\left[\frac{K_{+}^{(b)}}{2}\left(\cos{(2\absolutevalue{\chi})}+1\right)-\frac{\chi}{2\absolutevalue{\chi}}\sin{(2\absolutevalue{\chi})}K_{+}-\frac{\chi}{2\chi^{*}}\left(\cos{(2\absolutevalue{\chi})}-1\right)K_{+}^{(a)}\right]\right]\\ -\frac{\lambda^{2}}{\sqrt{\omega^{2}-\lambda^{2}}}\left[e^{-2i\psi}\left[\frac{K_{+}^{(a)}}{2}\left(\cos{(2\absolutevalue{\chi})}+1\right)+\frac{\chi^{*}}{2\absolutevalue{\chi}}\sin{(2\absolutevalue{\chi})}K_{+}-\frac{\chi^{*}}{2\chi}\left(\cos{(2\absolutevalue{\chi})}-1\right)K_{+}^{(b)}\right]\right.\\ \left.+e^{2i\psi}\left[\frac{K_{-}^{(a)}}{2}\left(\cos{(2\absolutevalue{\chi})}+1\right)+\frac{\chi}{2\absolutevalue{\chi}}\sin{(2\absolutevalue{\chi})}K_{-}-\frac{\chi}{2\chi^{*}}\left(\cos{(2\absolutevalue{\chi})}-1\right)K_{-}^{(b)}\right]\right].\end{gathered}$$ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG italic_K start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL + divide start_ARG italic_ωおめが italic_λらむだ end_ARG start_ARG square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG [ italic_e start_POSTSUPERSCRIPT - italic_i italic_ψぷさい end_POSTSUPERSCRIPT [ - divide start_ARG italic_χかい start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT end_ARG start_ARG | start_ARG italic_χかい end_ARG | end_ARG italic_αあるふぁ start_POSTSUBSCRIPT italic_χかい end_POSTSUBSCRIPT italic_J start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_βべーた start_POSTSUBSCRIPT italic_χかい end_POSTSUBSCRIPT ( italic_J start_POSTSUBSCRIPT + end_POSTSUBSCRIPT + divide start_ARG italic_χかい start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT end_ARG start_ARG italic_χかい end_ARG italic_J start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) + italic_J start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ] + italic_e start_POSTSUPERSCRIPT italic_i italic_ψぷさい end_POSTSUPERSCRIPT [ - divide start_ARG italic_χかい end_ARG start_ARG | start_ARG italic_χかい end_ARG | end_ARG italic_αあるふぁ start_POSTSUBSCRIPT italic_χかい end_POSTSUBSCRIPT italic_J start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_βべーた start_POSTSUBSCRIPT italic_χかい end_POSTSUBSCRIPT ( italic_J start_POSTSUBSCRIPT - end_POSTSUBSCRIPT + divide start_ARG italic_χかい end_ARG start_ARG italic_χかい start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT end_ARG italic_J start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) + italic_J start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ] ] end_CELL end_ROW start_ROW start_CELL - divide start_ARG italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG [ [ divide start_ARG italic_K start_POSTSUBSCRIPT - end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT end_ARG start_ARG 2 end_ARG ( roman_cos ( start_ARG ( 2 | start_ARG italic_χかい end_ARG | ) end_ARG ) + 1 ) - divide start_ARG italic_χかい start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT end_ARG start_ARG 2 | start_ARG italic_χかい end_ARG | end_ARG roman_sin ( start_ARG ( 2 | start_ARG italic_χかい end_ARG | ) end_ARG ) italic_K start_POSTSUBSCRIPT - end_POSTSUBSCRIPT - divide start_ARG italic_χかい start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_χかい end_ARG ( roman_cos ( start_ARG ( 2 | start_ARG italic_χかい end_ARG | ) end_ARG ) - 1 ) italic_K start_POSTSUBSCRIPT - end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT ] end_CELL end_ROW start_ROW start_CELL + [ divide start_ARG italic_K start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT end_ARG start_ARG 2 end_ARG ( roman_cos ( start_ARG ( 2 | start_ARG italic_χかい end_ARG | ) end_ARG ) + 1 ) - divide start_ARG italic_χかい end_ARG start_ARG 2 | start_ARG italic_χかい end_ARG | end_ARG roman_sin ( start_ARG ( 2 | start_ARG italic_χかい end_ARG | ) end_ARG ) italic_K start_POSTSUBSCRIPT + end_POSTSUBSCRIPT - divide start_ARG italic_χかい end_ARG start_ARG 2 italic_χかい start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT end_ARG ( roman_cos ( start_ARG ( 2 | start_ARG italic_χかい end_ARG | ) end_ARG ) - 1 ) italic_K start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT ] ] end_CELL end_ROW start_ROW start_CELL - divide start_ARG italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG [ italic_e start_POSTSUPERSCRIPT - 2 italic_i italic_ψぷさい end_POSTSUPERSCRIPT [ divide start_ARG italic_K start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT end_ARG start_ARG 2 end_ARG ( roman_cos ( start_ARG ( 2 | start_ARG italic_χかい end_ARG | ) end_ARG ) + 1 ) + divide start_ARG italic_χかい start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT end_ARG start_ARG 2 | start_ARG italic_χかい end_ARG | end_ARG roman_sin ( start_ARG ( 2 | start_ARG italic_χかい end_ARG | ) end_ARG ) italic_K start_POSTSUBSCRIPT + end_POSTSUBSCRIPT - divide start_ARG italic_χかい start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_χかい end_ARG ( roman_cos ( start_ARG ( 2 | start_ARG italic_χかい end_ARG | ) end_ARG ) - 1 ) italic_K start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT ] end_CELL end_ROW start_ROW start_CELL + italic_e start_POSTSUPERSCRIPT 2 italic_i italic_ψぷさい end_POSTSUPERSCRIPT [ divide start_ARG italic_K start_POSTSUBSCRIPT - end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT end_ARG start_ARG 2 end_ARG ( roman_cos ( start_ARG ( 2 | start_ARG italic_χかい end_ARG | ) end_ARG ) + 1 ) + divide start_ARG italic_χかい end_ARG start_ARG 2 | start_ARG italic_χかい end_ARG | end_ARG roman_sin ( start_ARG ( 2 | start_ARG italic_χかい end_ARG | ) end_ARG ) italic_K start_POSTSUBSCRIPT - end_POSTSUBSCRIPT - divide start_ARG italic_χかい end_ARG start_ARG 2 italic_χかい start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT end_ARG ( roman_cos ( start_ARG ( 2 | start_ARG italic_χかい end_ARG | ) end_ARG ) - 1 ) italic_K start_POSTSUBSCRIPT - end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT ] ] . end_CELL end_ROW (11)

where $\alpha_{\chi}=\sin{(2\absolutevalue{\chi})}$αあるふぁ start_POSTSUBSCRIPT italic_χかい end_POSTSUBSCRIPT = roman_sin ( start_ARG ( 2 | start_ARG italic_χかい end_ARG | ) end_ARG ) and $\beta_{\chi}={}^{1}\!/_{2}[\cos{(2\absolutevalue{\chi})}-1]$βべーた start_POSTSUBSCRIPT italic_χかい end_POSTSUBSCRIPT = start_FLOATSUPERSCRIPT 1 end_FLOATSUPERSCRIPT / start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT [ roman_cos ( start_ARG ( 2 | start_ARG italic_χかい end_ARG | ) end_ARG ) - 1 ].
We can remove the dependence of the Hamiltonian $H^{\prime\prime}$ on the operators $J_{\pm}$ if the coherent state parameters are chosen as

$$\theta=(2n+1)\frac{\pi}{2}\hskip 8.5359pt\hbox{with}\hskip 8.5359ptn\in\boldsymbol{Z},\hskip 8.5359pt\phi_{\theta}=\psi.$$θしーた = ( 2 italic_n + 1 ) divide start_ARG italic_πぱい end_ARG start_ARG 2 end_ARG with italic_n ∈ bold_italic_Z , italic_ϕ start_POSTSUBSCRIPT italic_θしーた end_POSTSUBSCRIPT = italic_ψぷさい . (12)

In this case, since $\cos{(\theta)}=0$θしーた ) end_ARG ) = 0 and $\sin{(\theta)}=1$θしーた ) end_ARG ) = 1, we get the expression

$$\begin{gathered}H^{\prime\prime}=2\sqrt{\omega^{2}-\lambda^{2}}K_{0}+\frac{2\omega\lambda}{\sqrt{\omega^{2}-\lambda^{2}}}J_{0}-\frac{\lambda^{2}}{\sqrt{\omega^{2}-\lambda^{2}}}(K_{+}^{(b)}+K_{-}^{(b)})-\frac{\lambda^{2}}{\sqrt{\omega^{2}-\lambda^{2}}}(K_{+}^{(a)}e^{-2i\psi}+K_{-}^{(a)}e^{2i\psi}).\end{gathered}$$ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG italic_K start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + divide start_ARG 2 italic_ωおめが italic_λらむだ end_ARG start_ARG square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG italic_J start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - divide start_ARG italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG ( italic_K start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT + italic_K start_POSTSUBSCRIPT - end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT ) - divide start_ARG italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG ( italic_K start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT - 2 italic_i italic_ψぷさい end_POSTSUPERSCRIPT + italic_K start_POSTSUBSCRIPT - end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT 2 italic_i italic_ψぷさい end_POSTSUPERSCRIPT ) . end_CELL end_ROW (13)

Here it is important to note that to diagonalize $H^{\prime\prime}$, it is necessary to apply a new tilting transformation in terms of the $SU(1,1)$ displacement operator of each boson $a$ and $b$, which it is defined as

$$D(\xi)_{a,b}=\exp(\xi_{a}K^{(a)}_{+}-\xi^{*}_{a}K^{(a)}_{-})\exp(\xi_{b}K^{(b)}_{+}-\xi^{*}_{b}K^{(b)}_{-})=D(\xi_{a})D(\xi_{b}),$$ξくしー ) start_POSTSUBSCRIPT italic_a , italic_b end_POSTSUBSCRIPT = roman_exp ( start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_K start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT + end_POSTSUBSCRIPT - italic_ξくしー start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_K start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - end_POSTSUBSCRIPT end_ARG ) roman_exp ( start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_K start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT + end_POSTSUBSCRIPT - italic_ξくしー start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_K start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - end_POSTSUBSCRIPT end_ARG ) = italic_D ( italic_ξくしー start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ) italic_D ( italic_ξくしー start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ) , (14)

where $\xi_{a}=-\frac{\theta_{a}}{2}e^{-i\phi_{a}}$ξくしー start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = - divide start_ARG italic_θしーた start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG italic_e start_POSTSUPERSCRIPT - italic_i italic_ϕ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_POSTSUPERSCRIPT and $\xi_{b}=-\frac{\theta_{b}}{2}e^{-i\phi_{b}}$ξくしー start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT = - divide start_ARG italic_θしーた start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG italic_e start_POSTSUPERSCRIPT - italic_i italic_ϕ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_POSTSUPERSCRIPT. Therefore, with $H^{\prime\prime\prime}=D^{\dagger}(\xi)_{a,b}H^{\prime\prime}D(\xi)_{a,b}$ξくしー ) start_POSTSUBSCRIPT italic_a , italic_b end_POSTSUBSCRIPT italic_H start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT italic_D ( italic_ξくしー ) start_POSTSUBSCRIPT italic_a , italic_b end_POSTSUBSCRIPT and ${\mathrm{<ruby><rb>\Psi </rb><rt>ぷさい</rt></ruby>}}^{\prime \prime \prime }=D^{†}{(\mathrm{<ruby><rb>\xi </rb><rt>くしー</rt></ruby>})}_{a,b}{\mathrm{<ruby><rb>\Psi </rb><rt>ぷさい</rt></ruby>}}^{\prime \prime }$Ψぷさい′′′superscript𝐷†subscript𝜉𝑎𝑏superscriptΨぷさい′′\Psi^{\prime\prime\prime}=D^{\dagger}(\xi)_{a,b}\Psi^{\prime\prime}roman_Ψぷさい start_POSTSUPERSCRIPT ′ ′ ′ end_POSTSUPERSCRIPT = italic_D start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_ξくしー ) start_POSTSUBSCRIPT italic_a , italic_b end_POSTSUBSCRIPT roman_Ψぷさい start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT, and using the results of Eqs. (B.7) of Appendix B we obtain

$$\begin{gathered}H^{\prime\prime\prime}=2\sqrt{\omega^{2}-\lambda^{2}}\left[\frac{1}{2}\left(\cosh(2\absolutevalue{\xi_{a}})+\cosh(2\absolutevalue{\xi_{b}})\right)K_{0}+\frac{1}{2}\left(\cosh(2\absolutevalue{\xi_{a}})-\cosh(2\absolutevalue{\xi_{b}})\right)J_{0}\right.\\ \left.+\frac{\sinh(2\absolutevalue{\xi_{a}})}{2\absolutevalue{\xi_{a}}}\left(\xi^{*}_{a}K^{(a)}_{-}+\xi_{a}K^{(a)}_{+}\right)+\frac{\sinh(2\absolutevalue{\xi_{b}})}{2\absolutevalue{\xi_{b}}}\left(\xi^{*}_{b}K^{(b)}_{-}+\xi_{b}K^{(b)}_{+}\right)\right]\\ +\frac{2\omega\lambda}{\sqrt{\omega^{2}-\lambda^{2}}}\left[\frac{\sinh(2\absolutevalue{\xi_{a}})}{2\absolutevalue{\xi_{a}}}\left(\xi^{*}_{a}K^{(a)}_{-}+\xi_{a}K^{(a)}_{+}\right)-\frac{\sinh(2\absolutevalue{\xi_{b}})}{2\absolutevalue{\xi_{b}}}\left(\xi^{*}_{b}K^{(b)}_{-}+\xi_{b}K^{(b)}_{+}\right)\right.\\ \left.+\frac{1}{2}\left(\cosh(2\absolutevalue{\xi_{a}})-\cosh(2\absolutevalue{\xi_{b}})\right)K_{0}+\frac{1}{2}\left(\cosh(2\absolutevalue{\xi_{a}})+\cosh(2\absolutevalue{\xi_{b}})\right)J_{0}\right]\end{gathered}$$ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG [ divide start_ARG 1 end_ARG start_ARG 2 end_ARG ( roman_cosh ( start_ARG 2 | start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG | end_ARG ) + roman_cosh ( start_ARG 2 | start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG | end_ARG ) ) italic_K start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + divide start_ARG 1 end_ARG start_ARG 2 end_ARG ( roman_cosh ( start_ARG 2 | start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG | end_ARG ) - roman_cosh ( start_ARG 2 | start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG | end_ARG ) ) italic_J start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL + divide start_ARG roman_sinh ( 2 | start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG | ) end_ARG start_ARG 2 | start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG | end_ARG ( italic_ξくしー start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_K start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - end_POSTSUBSCRIPT + italic_ξくしー start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_K start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) + divide start_ARG roman_sinh ( 2 | start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG | ) end_ARG start_ARG 2 | start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG | end_ARG ( italic_ξくしー start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_K start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - end_POSTSUBSCRIPT + italic_ξくしー start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_K start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) ] end_CELL end_ROW start_ROW start_CELL + divide start_ARG 2 italic_ωおめが italic_λらむだ end_ARG start_ARG square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG [ divide start_ARG roman_sinh ( 2 | start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG | ) end_ARG start_ARG 2 | start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG | end_ARG ( italic_ξくしー start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_K start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - end_POSTSUBSCRIPT + italic_ξくしー start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_K start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) - divide start_ARG roman_sinh ( 2 | start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG | ) end_ARG start_ARG 2 | start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG | end_ARG ( italic_ξくしー start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_K start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - end_POSTSUBSCRIPT + italic_ξくしー start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_K start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) end_CELL end_ROW start_ROW start_CELL + divide start_ARG 1 end_ARG start_ARG 2 end_ARG ( roman_cosh ( start_ARG 2 | start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG | end_ARG ) - roman_cosh ( start_ARG 2 | start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG | end_ARG ) ) italic_K start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + divide start_ARG 1 end_ARG start_ARG 2 end_ARG ( roman_cosh ( start_ARG 2 | start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG | end_ARG ) + roman_cosh ( start_ARG 2 | start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG | end_ARG ) ) italic_J start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ] end_CELL end_ROW

$$\begin{gathered}-\frac{\lambda^{2}}{\sqrt{\omega^{2}-\lambda^{2}}}\left[\frac{\xi_{b}}{\absolutevalue{\xi_{b}}}\alpha_{b}K^{(b)}_{0}+\beta_{b}\left(K_{-}^{(b)}+\frac{\xi_{b}}{\xi^{*}_{b}}K_{+}^{(b)}\right)+K^{(b)}_{-}\right]\\ -\frac{\lambda^{2}}{\sqrt{\omega^{2}-\lambda^{2}}}\left[\frac{\xi^{*}_{b}}{\absolutevalue{\xi_{b}}}\alpha_{b}K^{(b)}_{0}+\beta_{b}\left(K_{+}^{(b)}+\frac{\xi^{*}_{b}}{\xi_{b}}K_{-}^{(b)}\right)+K^{(b)}_{+}\right]\\ -\frac{\lambda^{2}}{\sqrt{\omega^{2}-\lambda^{2}}}e^{-2i\psi}\left[\frac{\xi^{*}_{a}}{\absolutevalue{\xi_{a}}}\alpha_{a}K^{(a)}_{0}+\beta_{a}\left(K_{+}^{(a)}+\frac{\xi^{*}_{a}}{\xi_{a}}K_{-}^{(a)}\right)+K^{(a)}_{+}\right]\\ -\frac{\lambda^{2}}{\sqrt{\omega^{2}-\lambda^{2}}}e^{2i\psi}\left[\frac{\xi_{a}}{\absolutevalue{\xi_{a}}}\alpha_{a}K^{(a)}_{0}+\beta_{a}\left(K_{-}^{(a)}+\frac{\xi_{a}}{\xi^{*}_{a}}K_{+}^{(a)}\right)+K^{(a)}_{-}\right].\end{gathered}$$λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG [ divide start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG start_ARG | start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG | end_ARG italic_αあるふぁ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_K start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_βべーた start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ( italic_K start_POSTSUBSCRIPT - end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT + divide start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG start_ARG italic_ξくしー start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG italic_K start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT ) + italic_K start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ] end_CELL end_ROW start_ROW start_CELL - divide start_ARG italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG [ divide start_ARG italic_ξくしー start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG start_ARG | start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG | end_ARG italic_αあるふぁ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_K start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_βべーた start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ( italic_K start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT + divide start_ARG italic_ξくしー start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG italic_K start_POSTSUBSCRIPT - end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT ) + italic_K start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ] end_CELL end_ROW start_ROW start_CELL - divide start_ARG italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG italic_e start_POSTSUPERSCRIPT - 2 italic_i italic_ψぷさい end_POSTSUPERSCRIPT [ divide start_ARG italic_ξくしー start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG start_ARG | start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG | end_ARG italic_αあるふぁ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_K start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_βべーた start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ( italic_K start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT + divide start_ARG italic_ξくしー start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG italic_K start_POSTSUBSCRIPT - end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT ) + italic_K start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ] end_CELL end_ROW start_ROW start_CELL - divide start_ARG italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG italic_e start_POSTSUPERSCRIPT 2 italic_i italic_ψぷさい end_POSTSUPERSCRIPT [ divide start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG start_ARG | start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG | end_ARG italic_αあるふぁ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_K start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_βべーた start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ( italic_K start_POSTSUBSCRIPT - end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT + divide start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG start_ARG italic_ξくしー start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG italic_K start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT ) + italic_K start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ] . end_CELL end_ROW (15)

In this expression $\alpha_{a}=\sinh(2\absolutevalue{\xi_{a}})$αあるふぁ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = roman_sinh ( 2 | start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG | ), $\beta_{a}=\frac{1}{2}[\cosh(2\absolutevalue{\xi_{a}})-1]$βべーた start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG [ roman_cosh ( start_ARG 2 | start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG | end_ARG ) - 1 ], $\alpha_{b}=\sinh(2\absolutevalue{\xi_{b}})$αあるふぁ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT = roman_sinh ( 2 | start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG | ) and $\beta_{b}=\frac{1}{2}[\cosh(2\absolutevalue{\xi_{b}})-1]$βべーた start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG [ roman_cosh ( start_ARG 2 | start_ARG italic_ξくしー start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG | end_ARG ) - 1 ]. The coefficients of the operators $K^{(a)}_{\pm}$ and $K^{(b)}_{\pm}$ will vanish if we set the values of the coherent state parameters given by

$$\begin{gathered}\theta_{a}=\tanh^{-1}{\left(\frac{-\lambda^{2}}{\omega^{2}-\lambda^{2}+\omega\lambda}\right)},\hskip 8.5359pt\phi_{a}=2\psi,\\ \theta_{b}=\tanh^{-1}{\left(\frac{-\lambda^{2}}{\omega^{2}-\lambda^{2}-\omega\lambda}\right)},\hskip 8.5359pt\phi_{b}=0.\end{gathered}$$θしーた start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = roman_tanh start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ( divide start_ARG - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_ωおめが italic_λらむだ end_ARG ) , italic_ϕ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = 2 italic_ψぷさい , end_CELL end_ROW start_ROW start_CELL italic_θしーた start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT = roman_tanh start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ( divide start_ARG - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ωおめが italic_λらむだ end_ARG ) , italic_ϕ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT = 0 . end_CELL end_ROW (16)

Therefore, we obtain that the Hamiltonian $H^{\prime\prime\prime}$ is diagonal and is reduced to

$$H^{\prime\prime\prime} = \left[\sqrt{\omega^{2}-\lambda^{2}}[\cosh(\theta_{a})+\cosh(\theta_{b})]+\frac{\omega\lambda}{\sqrt{\omega^{2}-\lambda^{2}}}[\cosh(\theta_{a})-\cosh(\theta_{b})]\right]K_{0}$$ ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG [ roman_cosh ( start_ARG italic_θしーた start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG ) + roman_cosh ( start_ARG italic_θしーた start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG ) ] + divide start_ARG italic_ωおめが italic_λらむだ end_ARG start_ARG square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG [ roman_cosh ( start_ARG italic_θしーた start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG ) - roman_cosh ( start_ARG italic_θしーた start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG ) ] ] italic_K start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT (17)

$$+\left[\sqrt{\omega^{2}-\lambda^{2}}[\cosh(\theta_{a})-\cosh(\theta_{b})]+\frac{\omega\lambda}{\sqrt{\omega^{2}-\lambda^{2}}}[\cosh(\theta_{a})+\cosh(\theta_{b})]\right]J_{0}$$ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG [ roman_cosh ( start_ARG italic_θしーた start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG ) - roman_cosh ( start_ARG italic_θしーた start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG ) ] + divide start_ARG italic_ωおめが italic_λらむだ end_ARG start_ARG square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG [ roman_cosh ( start_ARG italic_θしーた start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG ) + roman_cosh ( start_ARG italic_θしーた start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG ) ] ] italic_J start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT

$$+\frac{2\lambda^{2}}{\sqrt{\omega^{2}-\lambda^{2}}}\sinh(\theta_{a})K_{0}^{(a)}+\frac{2\lambda^{2}}{\sqrt{\omega^{2}-\lambda^{2}}}\sinh(\theta_{b})K_{0}^{(b)}.$$λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG roman_sinh ( italic_θしーた start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ) italic_K start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT + divide start_ARG 2 italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG roman_sinh ( italic_θしーた start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ) italic_K start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT .

We can write the Hamiltonian in Eq. (17) only in terms of the operators $K_{0}$ and $J_{0}$ considering that $K_{0}=K_{0}^{(a)}+K_{0}^{(b)}$ and $J_{0}=K_{0}^{(a)}-K_{0}^{(b)}$. Therefore, we obtain the following result for the Hamiltonian $H^{\prime\prime\prime}$

$$\begin{gathered}H^{\prime\prime\prime}=\left[\sqrt{\omega^{2}-\lambda^{2}}[\cosh(\theta_{a})+\cosh(\theta_{b})]+\frac{\omega\lambda}{\sqrt{\omega^{2}-\lambda^{2}}}[\cosh(\theta_{a})-\cosh(\theta_{b})]+\frac{\lambda^{2}}{\sqrt{\omega^{2}-\lambda^{2}}}[\sinh(\theta_{a})+\sinh(\theta_{b})]\right]K_{0}\\ +\left[\sqrt{\omega^{2}-\lambda^{2}}[\cosh(\theta_{a})-\cosh(\theta_{b})]+\frac{\omega\lambda}{\sqrt{\omega^{2}-\lambda^{2}}}[\cosh(\theta_{a})+\cosh(\theta_{b})]+\frac{\lambda^{2}}{\sqrt{\omega^{2}-\lambda^{2}}}[\sinh(\theta_{a})-\sinh(\theta_{b})]\right]J_{0},\end{gathered}$$ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG [ roman_cosh ( start_ARG italic_θしーた start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG ) + roman_cosh ( start_ARG italic_θしーた start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG ) ] + divide start_ARG italic_ωおめが italic_λらむだ end_ARG start_ARG square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG [ roman_cosh ( start_ARG italic_θしーた start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG ) - roman_cosh ( start_ARG italic_θしーた start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG ) ] + divide start_ARG italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG [ roman_sinh ( italic_θしーた start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ) + roman_sinh ( italic_θしーた start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ) ] ] italic_K start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL + [ square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG [ roman_cosh ( start_ARG italic_θしーた start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG ) - roman_cosh ( start_ARG italic_θしーた start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG ) ] + divide start_ARG italic_ωおめが italic_λらむだ end_ARG start_ARG square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG [ roman_cosh ( start_ARG italic_θしーた start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG ) + roman_cosh ( start_ARG italic_θしーた start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG ) ] + divide start_ARG italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG [ roman_sinh ( italic_θしーた start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ) - roman_sinh ( italic_θしーた start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ) ] ] italic_J start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , end_CELL end_ROW (18)

where $\cosh(\theta_{a})$θしーた start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG ), $\sinh(\theta_{a})$θしーた start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ), $\cosh(\theta_{b})$θしーた start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG ) and $\sinh(\theta_{b})$θしーた start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ) take the values

$$\begin{gathered}\cosh(\theta_{a})=\frac{\omega^{2}-\lambda^{2}+\omega\lambda}{\sqrt{(\omega^{2}-\lambda^{2}+\omega\lambda)^{2}-\lambda^{4}}},\quad\quad\sinh(\theta_{a})=-\frac{\lambda^{2}}{\sqrt{(\omega^{2}-\lambda^{2}+\omega\lambda)^{2}-\lambda^{4}}},\\ \cosh(\theta_{b})=\frac{\omega^{2}-\lambda^{2}-\omega\lambda}{\sqrt{(\omega^{2}-\lambda^{2}-\omega\lambda)^{2}-\lambda^{4}}},\quad\quad\sinh(\theta_{a})=-\frac{\lambda^{2}}{\sqrt{(\omega^{2}-\lambda^{2}-\omega\lambda)^{2}-\lambda^{4}}}.\end{gathered}$$θしーた start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG ) = divide start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_ωおめが italic_λらむだ end_ARG start_ARG square-root start_ARG ( italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_ωおめが italic_λらむだ ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG end_ARG , roman_sinh ( italic_θしーた start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ) = - divide start_ARG italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG square-root start_ARG ( italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_ωおめが italic_λらむだ ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG end_ARG , end_CELL end_ROW start_ROW start_CELL roman_cosh ( start_ARG italic_θしーた start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG ) = divide start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ωおめが italic_λらむだ end_ARG start_ARG square-root start_ARG ( italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ωおめが italic_λらむだ ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG end_ARG , roman_sinh ( italic_θしーた start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ) = - divide start_ARG italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG square-root start_ARG ( italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ωおめが italic_λらむだ ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG end_ARG . end_CELL end_ROW (19)

As we know, the operator $K_{0}$ represents the Hamiltonian of the two-dimensional harmonic oscillator, which commutes with $J_{0}$, so the eigenfunctions of $H^{\prime\prime\prime}$ turn out to be

$${\mathrm{<ruby><rb>\Psi </rb><rt>ぷさい</rt></ruby>}}_{N,m}(r,\varphi )=\frac{1}{\sqrt{\mathrm{<ruby><rb>\pi </rb><rt>ぱい</rt></ruby>}}}{e}^{im\varphi }{(-1)}^{\frac{N-m}{2}}\sqrt{\frac{2\left(\frac{N-m}{2}\right)!}{\left(\frac{N+m}{2}\right)!}}{r}^m{L}_{\frac{1}{2}(N-m)}^m(r^2)e^{-\frac{1}{2}{r}^2},$$Ψぷさい𝑁𝑚𝑟italic-ϕ1𝜋superscript𝑒𝑖𝑚italic-ϕsuperscript1𝑁𝑚22𝑁𝑚2𝑁𝑚2superscript𝑟𝑚superscriptsubscript𝐿12𝑁𝑚𝑚superscript𝑟2superscript𝑒12superscript𝑟2\Psi_{N,m}(r,\phi)=\frac{1}{\sqrt{\pi}}e^{im\phi}(-1)^{\frac{N-m}{2}}\sqrt{% \frac{2\left(\frac{N-m}{2}\right)!}{\left(\frac{N+m}{2}\right)!}}r^{m}L_{\frac% {1}{2}(N-m)}^{m}(r^{2})e^{-\frac{1}{2}r^{2}},roman_Ψぷさい start_POSTSUBSCRIPT italic_N , italic_m end_POSTSUBSCRIPT ( italic_r , italic_ϕ ) = divide start_ARG 1 end_ARG start_ARG square-root start_ARG italic_πぱい end_ARG end_ARG italic_e start_POSTSUPERSCRIPT italic_i italic_m italic_ϕ end_POSTSUPERSCRIPT ( - 1 ) start_POSTSUPERSCRIPT divide start_ARG italic_N - italic_m end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT square-root start_ARG divide start_ARG 2 ( divide start_ARG italic_N - italic_m end_ARG start_ARG 2 end_ARG ) ! end_ARG start_ARG ( divide start_ARG italic_N + italic_m end_ARG start_ARG 2 end_ARG ) ! end_ARG end_ARG italic_r start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT italic_L start_POSTSUBSCRIPT divide start_ARG 1 end_ARG start_ARG 2 end_ARG ( italic_N - italic_m ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT ( italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) italic_e start_POSTSUPERSCRIPT - divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT , (20)

or

$${\mathrm{<ruby><rb>\Psi </rb><rt>ぷさい</rt></ruby>}}_{n,m}(r,\varphi )=\frac{1}{\sqrt{\mathrm{<ruby><rb>\pi </rb><rt>ぱい</rt></ruby>}}}{e}^{im\varphi }{(-1)}^{n_r}\sqrt{\frac{2\left(n_r\right)!}{\left(n_r+m\right)!}}{r}^m{L}_n^m(r^2)e^{-\frac{1}{2}{r}^2},$$Ψぷさい𝑛𝑚𝑟italic-ϕ1𝜋superscript𝑒𝑖𝑚italic-ϕsuperscript1subscript𝑛𝑟2subscript𝑛𝑟subscript𝑛𝑟𝑚superscript𝑟𝑚superscriptsubscript𝐿𝑛𝑚superscript𝑟2superscript𝑒12superscript𝑟2\Psi_{n,m}(r,\phi)=\frac{1}{\sqrt{\pi}}e^{im\phi}(-1)^{n_{r}}\sqrt{\frac{2% \left(n_{r}\right)!}{\left(n_{r}+m\right)!}}r^{m}L_{n}^{m}(r^{2})e^{-\frac{1}{% 2}r^{2}},roman_Ψぷさい start_POSTSUBSCRIPT italic_n , italic_m end_POSTSUBSCRIPT ( italic_r , italic_ϕ ) = divide start_ARG 1 end_ARG start_ARG square-root start_ARG italic_πぱい end_ARG end_ARG italic_e start_POSTSUPERSCRIPT italic_i italic_m italic_ϕ end_POSTSUPERSCRIPT ( - 1 ) start_POSTSUPERSCRIPT italic_n start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT end_POSTSUPERSCRIPT square-root start_ARG divide start_ARG 2 ( italic_n start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT ) ! end_ARG start_ARG ( italic_n start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT + italic_m ) ! end_ARG end_ARG italic_r start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT italic_L start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT ( italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) italic_e start_POSTSUPERSCRIPT - divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT , (21)

where $n=\frac{1}{2}(N-m)$ is the radial quantum number. Therefore, the eigenfunctions of the Hamiltonian of two coupled harmonic oscillators with identical free frequencies are given by

$$\mathrm{<ruby><rb>\Psi </rb><rt>ぷさい</rt></ruby>}=D(\mathrm{<ruby><rb>\xi </rb><rt>くしー</rt></ruby>})D(\mathrm{<ruby><rb>\chi </rb><rt>かい</rt></ruby>})D{(\mathrm{<ruby><rb>\xi </rb><rt>くしー</rt></ruby>})}_{a,b}{\mathrm{<ruby><rb>\Psi </rb><rt>ぷさい</rt></ruby>}}^{\prime \prime \prime }=D(\mathrm{<ruby><rb>\xi </rb><rt>くしー</rt></ruby>})D(\mathrm{<ruby><rb>\chi </rb><rt>かい</rt></ruby>})D{(\mathrm{<ruby><rb>\xi </rb><rt>くしー</rt></ruby>})}_{a,b}{\mathrm{<ruby><rb>\Psi </rb><rt>ぷさい</rt></ruby>}}_{N,m}(r,\varphi ).$$Ψぷさい𝐷𝜉𝐷𝜒𝐷subscript𝜉𝑎𝑏superscriptΨぷさい′′′𝐷𝜉𝐷𝜒𝐷subscript𝜉𝑎𝑏subscriptΨぷさい𝑁𝑚𝑟italic-ϕ\Psi=D(\xi)D(\chi)D(\xi)_{a,b}\Psi^{\prime\prime\prime}=D(\xi)D(\chi)D(\xi)_{a% ,b}\Psi_{N,m}(r,\phi).roman_Ψぷさい = italic_D ( italic_ξくしー ) italic_D ( italic_χかい ) italic_D ( italic_ξくしー ) start_POSTSUBSCRIPT italic_a , italic_b end_POSTSUBSCRIPT roman_Ψぷさい start_POSTSUPERSCRIPT ′ ′ ′ end_POSTSUPERSCRIPT = italic_D ( italic_ξくしー ) italic_D ( italic_χかい ) italic_D ( italic_ξくしー ) start_POSTSUBSCRIPT italic_a , italic_b end_POSTSUBSCRIPT roman_Ψぷさい start_POSTSUBSCRIPT italic_N , italic_m end_POSTSUBSCRIPT ( italic_r , italic_ϕ ) . (22)

It can be shown from the theory of irreducible representations of the $SU(1,1)$ and $SU(2)$ groups (see Appendix A) that the action of the operators $K_{0}$ and $J_{0}$ on the basis $|N,m\rangle$ are explicitly [24, 25]

$$K_{0}\ket{N,m}=\frac{1}{2}(\hat{a}^{{\dagger}}\hat{a}+\hat{b}^{{\dagger}}\hat{b}+1)\ket{N,m}=\frac{1}{2}(N+1)\ket{N,m},$$ (23)

$$J_{0}|N,m\rangle=\frac{1}{2}\left(a^{{\dagger}}a-b^{{\dagger}}b\right)|N,m\rangle=\frac{m}{2}|N,m\rangle.$$ (24)

Therefore, by substituting the Eqs. (23) and (24) into Eq. (18) we obtain that the energy eigenvalues in the $\ket{N,m}$ basis are

$$\begin{gathered}E=\left[\sqrt{\omega^{2}-\lambda^{2}}[\cosh(\theta_{a})+\cosh(\theta_{b})]+\frac{\omega\lambda}{\sqrt{\omega^{2}-\lambda^{2}}}[\cosh(\theta_{a})-\cosh(\theta_{b})]+\frac{\lambda^{2}}{\sqrt{\omega^{2}-\lambda^{2}}}[\sinh(\theta_{a})+\sinh(\theta_{b})]\right]\frac{(N+1)}{2}\\ +\left[\sqrt{\omega^{2}-\lambda^{2}}[\cosh(\theta_{a})-\cosh(\theta_{b})]+\frac{\omega\lambda}{\sqrt{\omega^{2}-\lambda^{2}}}[\cosh(\theta_{a})+\cosh(\theta_{b})]+\frac{\lambda^{2}}{\sqrt{\omega^{2}-\lambda^{2}}}[\sinh(\theta_{a})-\sinh(\theta_{b})]\right]\frac{m}{2}.\end{gathered}$$ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG [ roman_cosh ( start_ARG italic_θしーた start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG ) + roman_cosh ( start_ARG italic_θしーた start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG ) ] + divide start_ARG italic_ωおめが italic_λらむだ end_ARG start_ARG square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG [ roman_cosh ( start_ARG italic_θしーた start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG ) - roman_cosh ( start_ARG italic_θしーた start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG ) ] + divide start_ARG italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG [ roman_sinh ( italic_θしーた start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ) + roman_sinh ( italic_θしーた start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ) ] ] divide start_ARG ( italic_N + 1 ) end_ARG start_ARG 2 end_ARG end_CELL end_ROW start_ROW start_CELL + [ square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG [ roman_cosh ( start_ARG italic_θしーた start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG ) - roman_cosh ( start_ARG italic_θしーた start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG ) ] + divide start_ARG italic_ωおめが italic_λらむだ end_ARG start_ARG square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG [ roman_cosh ( start_ARG italic_θしーた start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG ) + roman_cosh ( start_ARG italic_θしーた start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG ) ] + divide start_ARG italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG square-root start_ARG italic_ωおめが start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG [ roman_sinh ( italic_θしーた start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ) - roman_sinh ( italic_θしーた start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ) ] ] divide start_ARG italic_m end_ARG start_ARG 2 end_ARG . end_CELL end_ROW (25)

Here is important to note that, if we set $\lambda=0$λらむだ = 0, we obtain that $\cosh(\theta_{a})=\cosh(\theta_{b})=1$θしーた start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG ) = roman_cosh ( start_ARG italic_θしーた start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG ) = 1 and $\sinh(\theta_{a})=\sinh(\theta_{b})=0$θしーた start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ) = roman_sinh ( italic_θしーた start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ) = 0. With this consideration we obtain that the energy spectrum is reduced to

$$E=\mathrm{<ruby><rb>\omega </rb><rt>おめが</rt></ruby>}(N+1)=\mathrm{<ruby><rb>\omega </rb><rt>おめが</rt></ruby>}(2n+m+1)=2\mathrm{<ruby><rb>\omega </rb><rt>おめが</rt></ruby>}($$