b}$. Compared with $H_{1}^{a}$ in the single-qubit ASRM of Sec. II, the main difference is the presence of the $J_{z}^{2}$ operator term in $H_{2}^{a}$, but in the single-qubit ASRM the corresponding term $\sigma_z^2=1$ is just a constant. Therefore, $H_{2}^{a}$ here represents a renormalized three-level system in which we need to diagonalize $H_{2}^{a}$ to remove counterrotating-wave terms.

The eigenvalues $\nu_{k}$ ($k=1,2,3$) and eigenstates $|\varphi_{k}\rangle_{A}$ of the Hamiltonian $H_{2}^{''}=H_{2}^{a}-w_{b}b^{\dagger}b$ are: $$\begin{aligned}
\label{25}
\nu_{1}&=&\frac{A}{2}-\frac{1}{2}\sqrt{A^{2}+8B^2}, \cr
|\varphi_{1}\rangle_{A}&=&\frac{1}{N_{1}}\bigg\{|-1\rangle_{A}-\frac{(A+\sqrt{A^2+8B^2})}{2B}|0\rangle_{A}+|1\rangle_{A}\bigg\},\cr\cr
\nu_{2}&=&A, \cr\cr
|\varphi_{2}\rangle_{A}&=&\frac{1}{N_{2}}\bigg\{-|-1\rangle_{A}+|1\rangle_{A}\bigg\},\cr\cr
\nu_{3}&=&\frac{A}{2}+\frac{1}{2}\sqrt{A^{2}+8B^2}, \cr\cr
|\varphi_{3}\rangle_{A}&=&\frac{1}{N_{3}}\bigg\{|-1\rangle_{A}-\frac{(A-\sqrt{A^2+8B^2})}{2B}|0\rangle_{A}+|1\rangle_{A}\bigg\},\cr&&\end{aligned}$$ with $$\begin{aligned}
\label{26}
A&=&w_{b}\xi_{2}^2-2\xi_{2}(g_{1}+g_{2}), \cr
B&=&\frac{1}{\sqrt{2}}[w_{a}\eta_{2}-\eta_{2}\xi_{2}(g_{1}-g_{
b}$. Compared with $ H_{1}^{a}$ in the single - qubit ASRM of Sec. II, the main difference is the bearing of the $ J_{z}^{2}$ hustler terminus in $ H_{2}^{a}$, but in the single - qubit ASRM the comparable term $ \sigma_z^2=1 $ is precisely a constant. Therefore, $ H_{2}^{a}$ here represent a renormalized three - level system in which we necessitate to diagonalize $ H_{2}^{a}$ to remove counterrotating - wave terms. 

 The eigenvalues $ \nu_{k}$ ($ k=1,2,3 $) and eigenstates $ |\varphi_{k}\rangle_{A}$ of the Hamiltonian $ H_{2}^{''}=H_{2}^{a}-w_{b}b^{\dagger}b$ are: $ $ \begin{aligned } 
 \label{25 } 
 \nu_{1}&=&\frac{A}{2}-\frac{1}{2}\sqrt{A^{2}+8B^2 }, \cr 
 |\varphi_{1}\rangle_{A}&=&\frac{1}{N_{1}}\bigg\{|-1\rangle_{A}-\frac{(A+\sqrt{A^2 + 8B^2})}{2B}|0\rangle_{A}+|1\rangle_{A}\bigg\},\cr\cr 
 \nu_{2}&=&A, \cr\cr 
 |\varphi_{2}\rangle_{A}&=&\frac{1}{N_{2}}\bigg\{-|-1\rangle_{A}+|1\rangle_{A}\bigg\},\cr\cr 
 \nu_{3}&=&\frac{A}{2}+\frac{1}{2}\sqrt{A^{2}+8B^2 }, \cr\cr 
 |\varphi_{3}\rangle_{A}&=&\frac{1}{N_{3}}\bigg\{|-1\rangle_{A}-\frac{(A-\sqrt{A^2 + 8B^2})}{2B}|0\rangle_{A}+|1\rangle_{A}\bigg\},\cr&&\end{aligned}$$ with $ $ \begin{aligned } 
 \label{26 } 
 A&=&w_{b}\xi_{2}^2 - 2\xi_{2}(g_{1}+g_{2 }), \cr 
 B&=&\frac{1}{\sqrt{2}}[w_{a}\eta_{2}-\eta_{2}\xi_{2}(g_{1}-g _ {
b}$. Clmpared with $H_{1}^{a}$ in the slngle-qubit ASRM of Sec. MI, the jain difwerence is the presence of tie $J_{z}^{2}$ opeeator term in $H_{2}^{a}$, but iv the sinhle-qubit ASRN the corresponding term $\algma_z^2=1$ ms just a constsnt. Therefmre, $H_{2}^{a}$ here re[rdsznts a renormalized three-level systei in whocj we need to dyagomwlizs $H_{2}^{a}$ to remove counterrotating-wave terms.

Tie eigenvalues $\mu_{k}$ ($k=1,2,3$) and eigenstates $|\varpji_{k}\rwngle_{A}$ of the Hamiptonian $H_{2}^{''}=H_{2}^{a}-q_{b}b^{\dwtger}b$ are: $$\beein{aligned}
\label{25}
\nu_{1}&=&\frac{A}{2}-\rrac{1}{2}\sqrt{A^{2}+8B^2}, \cr
|\varphi_{1}\rangle_{A}&=&\frac{1}{N_{1}}\cigg\{|-1\rcngle_{A}-\frac{(A+\wqet{A^2+8H^2})}{2T}|0\rangle_{A}+|1\ranjle_{A}\bidg\},\cr\cr
\nu_{2}&=&A, \cr\gg
|\varphi_{2}\sangle_{A}&=&\grac{1}{N_{2}}\bigg\{-|-1\ranglc_{A}+|1\ranjle_{A}\vigg\},\cr\cr
\nu_{3}&=&\frac{A}{2}+\frac{1}{2}\sqct{A^{2}+8B^2}, \cr\cr
|\varphi_{3}\ranglg_{A}&=&\frac{1}{N_{3}}\bigc\{|-1\rcngle_{A}-\frac{(A-\sqrt{A^2+8B^2})}{2B}|0\rantlw_{A}+|1\rancle_{A}\tigg\},\ze&&\ena{aljgied}$$ with $$\hegmn{aligned}
\lagel{26}
A&=&w_{b}\xi_{2}^2-2\xi_{2}(g_{1}+t_{2}), \cr
B&=&\frac{1}{\sqrt{2}}[w_{a}\eta_{2}-\eua_{2}\xy_{2}(t_{1}-g_{
b}$. Compared with $H_{1}^{a}$ in the single-qubit Sec. the main is the presence in but in the ASRM the corresponding $\sigma_z^2=1$ is just a constant. Therefore, here represents a renormalized three-level system in which we need to diagonalize $H_{2}^{a}$ remove counterrotating-wave terms. The eigenvalues $\nu_{k}$ ($k=1,2,3$) and eigenstates $|\varphi_{k}\rangle_{A}$ of the Hamiltonian are: \label{25} \cr \nu_{2}&=&A, \cr\cr |\varphi_{2}\rangle_{A}&=&\frac{1}{N_{2}}\bigg\{-|-1\rangle_{A}+|1\rangle_{A}\bigg\},\cr\cr \nu_{3}&=&\frac{A}{2}+\frac{1}{2}\sqrt{A^{2}+8B^2}, \cr\cr |\varphi_{3}\rangle_{A}&=&\frac{1}{N_{3}}\bigg\{|-1\rangle_{A}-\frac{(A-\sqrt{A^2+8B^2})}{2B}|0\rangle_{A}+|1\rangle_{A}\bigg\},\cr&&\end{aligned}$$ with $$\begin{aligned} \label{26} A&=&w_{b}\xi_{2}^2-2\xi_{2}(g_{1}+g_{2}), \cr B&=&\frac{1}{\sqrt{2}}[w_{a}\eta_{2}-\eta_{2}\xi_{2}(g_{1}-g_{
b}$. Compared with $H_{1}^{a}$ in the singlE-qubit ASRM Of Sec. iI, tHe mAiN difFereNce is the presenCE of tHe $J_{z}^{2}$ operator term in $H_{2}^{a}$, buT in thE sINgle-QUbIt ASRm the corREsPONdiNg TeRm $\sIgMA_z^2=1$ Is jusT a cOnstant. therefore, $H_{2}^{A}$ heRe Represents a rENoRmalized thRee-Level system iN whIch we nEeD to DIagonAliZe $H_{2}^{a}$ tO removE CounteRrotating-WaVE terms.

tHe eigenVALuEs $\nu_{K}$ ($k=1,2,3$) and eigenstates $|\vARpHI_{k}\rangle_{A}$ of the hamiltOnIAn $h_{2}^{''}=h_{2}^{A}-w_{b}B^{\daGger}b$ are: $$\beGiN{aligNEd}
\label{25}
\NU_{1}&=&\fRAC{a}{2}-\frAC{1}{2}\sqrt{A^{2}+8B^2}, \cr
|\varpHi_{1}\rangle_{A}&=&\frAC{1}{N_{1}}\bIgg\{|-1\ranGlE_{A}-\fRAc{(A+\sqrT{A^2+8B^2})}{2B}|0\rAnGLe_{A}+|1\Rangle_{A}\bigg\},\Cr\cr
\Nu_{2}&=&A, \cr\cr
|\vaRphi_{2}\raNGle_{A}&=&\fraC{1}{n_{2}}\bigg\{-|-1\raNgle_{A}+|1\rAngLe_{A}\Bigg\},\CR\cR
\nU_{3}&=&\frAc{a}{2}+\FraC{1}{2}\SqRt{A^{2}+8b^2}, \Cr\cR
|\varphi_{3}\rAnGlE_{A}&=&\fraC{1}{N_{3}}\biGG\{|-1\RANgle_{a}-\frAc{(A-\sQrt{A^2+8B^2})}{2b}|0\rangle_{A}+|1\ranglE_{A}\bIgg\},\cR&&\End{AlignEd}$$ witH $$\begIn{AlignEd}
\labeL{26}
A&=&w_{b}\xI_{2}^2-2\xI_{2}(g_{1}+g_{2}), \cr
B&=&\frac{1}{\sqrt{2}}[w_{A}\eta_{2}-\Eta_{2}\xi_{2}(g_{1}-g_{
b}$. Compared with $H_{1}^ {a}$ in th e sin gle -qu bi t AS RM o f Sec. II, the main  difference is the pre sence  o f  the $J _{z}^ {2}$ op e ra t o r t er min$H _ {2 }^{a} $,but inthe single -qu bi t ASRM the c o rr espondingter m $\sigma_z^ 2=1 $ is j us t a const ant . The refore ,  $H_{2 }^{a}$ he re repres e nts a r e n or mali zed three-level s y st e m in which weneed t od ia g o nal ize  $H_{2}^{a }$  to r e move co u nt e r r ota t ing-wave term s.

The eig e nva lues $ \n u_{ k }$ ($k =1,2, 3$ )  an d eigenstat es $ |\varphi_ {k}\ra n gle_{A} $  of the  Hamil ton ian  $H_ { 2} ^{ ''} =H _ {2} ^ {a }-w _ {b} b^{\dagg er }b $ are : $$ \ b e g in{a lig ned} 
\lab el{25}
\nu_{1 }&= &\fr a c{A }{2}- \frac {1}{ 2} \sqrt {A^{2} +8B^2 },  \cr
|\varphi_{ 1}\r angle_{A} &=& \f rac {1 }{N_{ 1 }}\big g\{ |-1 \rangle _{A}-\f r ac{ (A + \ s qr t{A^2+8B^2})}{2B}| 0\ r a ng le_{A}+| 1\rang l e_ {A } \bigg\}, \c r\c r
\n u _ {2}&= &A,\ cr \cr
|\va rphi_{ 2 }\ ra ngle_{A }& =&\fra c{ 1}{ N_{ 2}}\b i gg\{ -|-1\r angle_{A }+|1\ r angle_{A}\bigg \ },\cr\cr
\nu_ { 3} & = &\ f rac{ A}{ 2}+\frac{1} {2}\ s qrt{ A^{2 } +8 B^2 } , \cr \cr
| \v a rp h i_{3}\rangle_{A}&=& \f rac{1} {N_{3 }}\bigg\{|-1\ rangle_{A} - \ f rac{(A-\ sqrt { A^ 2 +8B^2})}{2B}|0 \rang le_{A}+|1\ r angle_{A }\big g\},\cr& &\end{ali g n ed}$$ wi th$$\ beg in{ a l ig ned}
\label{2 6 } 
A&= &w _{b}\xi _{2 }^2-2\x i_{ 2}( g_{ 1}+ g_ {2}), \cr 
B&=&\fr ac {1 }{ \s qrt {2}}[ w _{a}\eta _{ 2}- \e ta_ {2}\x i _{2}(g _{1}- g_{
b}$. Compared_with $H_{1}^{a}$_in the single-qubit ASRM_of Sec._II,_the main_difference_is the presence_of the $J_{z}^{2}$_operator term in $H_{2}^{a}$,_but in the_single-qubit_ASRM the corresponding term $\sigma_z^2=1$ is just a constant. Therefore, $H_{2}^{a}$ here represents a_renormalized_three-level system_in_which_we need to diagonalize $H_{2}^{a}$_to remove counterrotating-wave terms.

The eigenvalues_$\nu_{k}$ ($k=1,2,3$)_and eigenstates $|\varphi_{k}\rangle_{A}$ of the Hamiltonian $H_{2}^{''}=H_{2}^{a}-w_{b}b^{\dagger}b$ are:_$$\begin{aligned}
\label{25}
\nu_{1}&=&\frac{A}{2}-\frac{1}{2}\sqrt{A^{2}+8B^2},_\cr
|\varphi_{1}\rangle_{A}&=&\frac{1}{N_{1}}\bigg\{|-1\rangle_{A}-\frac{(A+\sqrt{A^2+8B^2})}{2B}|0\rangle_{A}+|1\rangle_{A}\bigg\},\cr\cr
\nu_{2}&=&A, \cr\cr
|\varphi_{2}\rangle_{A}&=&\frac{1}{N_{2}}\bigg\{-|-1\rangle_{A}+|1\rangle_{A}\bigg\},\cr\cr
\nu_{3}&=&\frac{A}{2}+\frac{1}{2}\sqrt{A^{2}+8B^2}, \cr\cr
|\varphi_{3}\rangle_{A}&=&\frac{1}{N_{3}}\bigg\{|-1\rangle_{A}-\frac{(A-\sqrt{A^2+8B^2})}{2B}|0\rangle_{A}+|1\rangle_{A}\bigg\},\cr&&\end{aligned}$$_with $$\begin{aligned}
\label{26}
A&=&w_{b}\xi_{2}^2-2\xi_{2}(g_{1}+g_{2}), \cr
B&=&\frac{1}{\sqrt{2}}[w_{a}\eta_{2}-\eta_{2}\xi_{2}(g_{1}-g_{