 the subsystem Hamiltonian, $H_b=\sum_{j=1}^NP_j^2/2+1/2\omega_j^2R_j^2=
\sum_{j=1}^NP_j^2/2+V_b(R)$ is the Hamiltonian of a classical bath of $N$ harmonic oscillators, and $\hat{V}_c(R)=-\sum_{j=1}^Nc_j\hat{\sigma}_zR_j=\gamma(R)\hat{\sigma}_z$ is the interaction between the subsystem and the bath. An Ohmic spectral density is assumed for the bath. Hence, denoting the Kondo parameter as $\xi_K$ and the cut-off frequency as $\omega_{\rm max}$, the frequencies of the oscillators are defined by $\omega_j=-\ln(1-j\omega_0)$, where $\omega_0=N^{-1}(1-\exp(-\omega_{\rm max}))$, and the constants entering the coupling by and $c_j=\sqrt{\xi_K\omega_0}~\omega_j$. The adiabatic eigenvalues and eigenvectors, respectively, are $$E_{1,2}=V_b\mp\sqrt{\Omega^2+\gamma^2(R)}\;,$$ $$\begin{aligned}
\vert 1;R\rangle&=& \frac{1}{\sqrt{2(1+G^2}}\left(\begin{array}{c} 1+G\\
1-G\end{array} \right)
\nonumber\\
\vert 2;R\rangle&=& \frac{1}{\sqrt{2(1+G^2}}\left(\begin{array}{c}-1+G\\
1+G\end{array} \right)\;,\end{aligned}$$ where $$G(R)=\gamma^{-1}(R)\left[-\Omega+\sqrt{\Omega^2+\gamma^2(R)}\right]\;.$$ The coupling vector $d_{\alpha\alpha'}=\langle\alpha;R\vert\overrightarrow{\partial}
/\partial R\vert\alpha';R\rangle$ is $$d_{12}=-d_{21}=(1+G^2)^{-1}\partial G/\partial R\;.$$ Assuming an initially uncorrelated density matrix, where the bath is in thermal equilibrium and the subsystem is in state $\vert\uparrow\rangle$, the initial quantum-classical density matrix
the subsystem Hamiltonian, $ H_b=\sum_{j=1}^NP_j^2/2 + 1/2\omega_j^2R_j^2= 
 \sum_{j=1}^NP_j^2/2+V_b(R)$ is the Hamiltonian of a classical bath of $ N$ harmonic oscillators, and $ \hat{V}_c(R)=-\sum_{j=1}^Nc_j\hat{\sigma}_zR_j=\gamma(R)\hat{\sigma}_z$ is the interaction between the subsystem and the bathtub. An Ohmic apparitional density is assumed for the bathtub. therefore, denoting the Kondo parameter as $ \xi_K$ and the cut - off frequency as $ \omega_{\rm max}$, the frequency of the oscillators are defined by $ \omega_j=-\ln(1 - j\omega_0)$, where $ \omega_0 = N^{-1}(1-\exp(-\omega_{\rm max}))$, and the constant entering the coupling by and $ c_j=\sqrt{\xi_K\omega_0}~\omega_j$. The adiabatic eigenvalue and eigenvectors, respectively, are $ $ E_{1,2}=V_b\mp\sqrt{\Omega^2+\gamma^2(R)}\;,$$ $ $ \begin{aligned } 
 \vert 1;R\rangle&= & \frac{1}{\sqrt{2(1+G^2}}\left(\begin{array}{c } 1+G\\ 
 1 - G\end{array } \right) 
 \nonumber\\ 
 \vert 2;R\rangle&= & \frac{1}{\sqrt{2(1+G^2}}\left(\begin{array}{c}-1+G\\ 
 1+G\end{array } \right)\;,\end{aligned}$$ where $ $ G(R)=\gamma^{-1}(R)\left[-\Omega+\sqrt{\Omega^2+\gamma^2(R)}\right]\;.$$ The coupling vector $ d_{\alpha\alpha'}=\langle\alpha;R\vert\overrightarrow{\partial } 
 /\partial R\vert\alpha';R\rangle$ is $ $ d_{12}=-d_{21}=(1+G^2)^{-1}\partial G/\partial R\;.$$ Assuming an initially uncorrelated density matrix, where the bath is in thermal chemical equilibrium and the subsystem is in state $ \vert\uparrow\rangle$, the initial quantum - classical density matrix
 thf subsystem Hamiltonian, $M_b=\sum_{j=1}^NP_j^2/2+1/2\omega_j^2R_l^2=
\wum_{j=1}^NP_o^2/2+V_b(R)$ is the Hamkltonian of a classical bath od $N$ hqrmonic oscillators, ana $\hat{V}_c(R)=-\slm_{j=1}^Nc_j\hat{\wigme}_zR_j=\gamma(R)\hat{\sigma}_z$ is tmz intsvactimi between the sobsystem and the bath. An Ovmkc spectral density is assumed for thq bath. Nejce, denoting tre Kpgdo larameter as $\xi_K$ and the cut-off frsquency as $\omega_{\rm msx}$, the frequencies of the lsciplators are definef by $\omega_j=-\on(1-j\oiwga_0)$, where $\omdga_0=N^{-1}(1-\exp(-\omega_{\rm max}))$, and the constants entering the couoling by and $c_j=\wqet{\xl_N\omega_0}~\omega_o$. The wdiabatic einvnvaluev and eogenvectors, rexpertivwly, are $$E_{1,2}=V_b\mp\sqrt{\Omege^2+\gamma^2(R)}\;,$$ $$\begin{aligned}
\dert 1;R\ranclz&=& \frac{1}{\sqrt{2(1+G^2}}\left(\begin{aerqy}{c} 1+G\\
1-C\end{drrah} \rieht)
\hoiumger\\
\verh 2;R\cangle&=& \frac{1}{\aqrt{2(1+G^2}}\left(\betin{array}{c}-1+G\\
1+G\end{array} \rytht)\;,\end{aligned}$$ where $$G(W)=\gamma^{-1}(R)\left[-\Omega+\sqrt{\Omega^2+\gamma^2(R)}\right]\;.$$ Tht coulling vector $d_{\alpha\alphq'}=\langle\alpha;R\vert\overtightarrow{\[artial}
/\partial R\vert\alpha';R\rangle$ is $$d_{12}=-d_{21}=(1+G^2)^{-1}\partial G/\pdrtiam R\;.$$ Cwsumivt wn initially uncorrelated density matrix, wherq tne bath is in thcrmal equilibrium snf yre subsystem ks in stafe $\vert\uparrow\rangpe$, the ynitiql quantui-claxsical density matrix
the subsystem Hamiltonian, $H_b=\sum_{j=1}^NP_j^2/2+1/2\omega_j^2R_j^2= \sum_{j=1}^NP_j^2/2+V_b(R)$ is the a bath of harmonic oscillators, and the and the bath. Ohmic spectral density assumed for the bath. Hence, denoting Kondo parameter as $\xi_K$ and the cut-off frequency as $\omega_{\rm max}$, the frequencies the oscillators are defined by $\omega_j=-\ln(1-j\omega_0)$, where $\omega_0=N^{-1}(1-\exp(-\omega_{\rm max}))$, and the constants entering coupling and The eigenvalues and eigenvectors, respectively, are $$E_{1,2}=V_b\mp\sqrt{\Omega^2+\gamma^2(R)}\;,$$ $$\begin{aligned} \vert 1;R\rangle&=& \frac{1}{\sqrt{2(1+G^2}}\left(\begin{array}{c} 1+G\\ 1-G\end{array} \right) \nonumber\\ \vert 2;R\rangle&=& 1+G\end{array} \right)\;,\end{aligned}$$ where $$G(R)=\gamma^{-1}(R)\left[-\Omega+\sqrt{\Omega^2+\gamma^2(R)}\right]\;.$$ The coupling vector $d_{\alpha\alpha'}=\langle\alpha;R\vert\overrightarrow{\partial} R\vert\alpha';R\rangle$ is $$d_{12}=-d_{21}=(1+G^2)^{-1}\partial G/\partial Assuming an initially uncorrelated density where bath is thermal and subsystem is in $\vert\uparrow\rangle$, the initial quantum-classical density matrix
 the subsystem Hamiltonian, $H_b=\Sum_{j=1}^NP_j^2/2+1/2\omeGa_j^2R_j^2=
\Sum_{J=1}^NP_J^2/2+V_B(R)$ is The HAmiltonian of a cLAssiCal bath of $N$ harmonic osciLlatoRs, ANd $\haT{v}_c(r)=-\sum_{j=1}^nc_j\hat{\sIGmA}_Zr_j=\gAmMa(r)\haT{\sIGmA}_z$ is tHe iNteractIon between The SuBsystem and thE BaTh. An Ohmic sPecTral density iS asSumed fOr The BAth. HeNce, DenotIng the kOndo paRameter as $\Xi_k$ And the CUt-off frEQUeNcy aS $\omega_{\rm max}$, the freQUeNCies of the oscilLators ArE DeFINed By $\oMega_j=-\ln(1-j\omEgA_0)$, wherE $\Omega_0=N^{-1}(1-\eXP(-\oMEGA_{\rm MAx}))$, and the constAnts enterinG The CoupliNg By aND $c_j=\sqrT{\xi_K\oMeGA_0}~\omEga_j$. The adiaBatiC eigenvalUes and EIgenvecTOrs, respEctiveLy, aRe $$E_{1,2}=v_b\mp\SQrT{\OMegA^2+\gAMma^2(r)}\;,$$ $$\BeGin{ALigNed}
\vert 1;R\RaNgLe&=& \fraC{1}{\sqrT{2(1+g^2}}\LEFt(\beGin{ArraY}{c} 1+G\\
1-G\eNd{array} \right)
\nOnuMber\\
\VErt 2;r\rangLe&=& \fraC{1}{\sqrT{2(1+G^2}}\Left(\bEgin{arRay}{c}-1+G\\
1+g\eNd{array} \right)\;,\end{AligNed}$$ where $$G(r)=\gaMmA^{-1}(R)\lEfT[-\OmegA+\Sqrt{\OmEga^2+\GamMa^2(R)}\righT]\;.$$ The couPLinG vECTOr $D_{\alpha\alpha'}=\langle\aLpHA;r\vErt\overrIghtarROw{\PaRTial}
/\partIaL R\vErt\aLPHa';R\raNgle$ IS $$d_{12}=-D_{21}=(1+G^2)^{-1}\partiaL G/\partIAl r\;.$$ ASsuming An InitiaLlY unCorRelatED denSity maTrix, wherE the bATh is in thermal eQUilibrium and tHE sUBSySTem iS in State $\vert\upArroW\RangLe$, thE InItiAL quanTum-clAsSIcAL density matrix
 the subsystem Hamiltonian , $H_b=\su m_{j= 1}^ NP_ j^ 2/2+ 1/2\ omega_j^2R_j^2 = 
\su m_{j=1}^NP_j^2/2+V_b(R )$ is  t h e Ha m il tonia n of ac la s s ica lba thof $N $ har mon ic osci llators, a nd$\ hat{V}_c(R)= - \s um_{j=1}^N c_j \hat{\sigma} _zR _j=\ga mm a(R ) \hat{ \si gma}_ z$ ist he int eractionbe t ween t h e subsy s t em  and  the bath. An Ohm i cs pectral densit y is a ss u me d for  th e bath. He nc e, de n oting t h eK o n dop arameter as $ \xi_K$ andt hecut-of ffre q uencyas $\ om e ga_ {\rm max}$,  the  frequenc ies of the osc i llators  are d efi ned  by$ \o me ga_ j= - \ln ( 1- j\o m ega _0)$, wh er e$\ome ga_0 = N ^ { -1}( 1-\ exp( -\ome ga_{\rm max}) )$,  and the  cons tants  ent er ing t he cou pling  b y and $c_j=\sqr t{\x i_K\omega _0} ~\ ome ga _j$.T he adi aba tic  eigenv alues a n d e ig e n v ec tors, respectively ,a r e$$E_{1,2 }=V_b\ m p\ sq r t{\Omega ^2 +\g amma ^ 2 (R)}\ ;,$$ $$ \begin{a ligned } 
\ ve rt 1;R\ ra ngle&= &\fr ac{ 1}{\s q rt{2 (1+G^2 }}\left( \begi n {array}{c} 1+G \ \
1-G\end{arr a y} \ ri g ht)\no number\\
\v ert2 ;R\r angl e &= & \ f rac{1 }{\sq rt { 2( 1 +G^2}}\left(\begin{ ar ray}{c }-1+G \\
1+G\end{ar ray} \righ t ) \ ;,\end{a lign e d} $ $ where $$G(R) =\gam ma^{-1}(R) \ left[-\O mega+ \sqrt{\O mega^2+\g a m ma^2(R)} \ri ght ]\; .$$ T he  coupling vec t o r $d _{ \alpha\ alp ha'}=\l ang le\ alp ha; R\ vert\over rightarr ow {\ pa rt ial }
/\p a rtial R\ ve rt\ al pha ';R\r a ngle$is $$ d_{1 2} =- d _{2 1}=(1+G ^ 2) ^ { -1}\ pa rt ialG/\ pa rtial  R\; . $$Assumin g an init ial l y un co rr elateddensity matri x,  where the  b ath  is in t hermal e quilibrium and the subs y stem is  in  stat e $\ vert\upar row \rangl e$, the in itialquant um -cl a s sical d en sit ymatrix
 the_subsystem Hamiltonian,_$H_b=\sum_{j=1}^NP_j^2/2+1/2\omega_j^2R_j^2=
\sum_{j=1}^NP_j^2/2+V_b(R)$ is the Hamiltonian_of a_classical_bath of_$N$_harmonic oscillators, and_$\hat{V}_c(R)=-\sum_{j=1}^Nc_j\hat{\sigma}_zR_j=\gamma(R)\hat{\sigma}_z$ is the_interaction between the subsystem_and the bath._An_Ohmic spectral density is assumed for the bath. Hence, denoting the Kondo parameter as_$\xi_K$_and the_cut-off_frequency_as $\omega_{\rm max}$, the frequencies_of the oscillators are defined_by $\omega_j=-\ln(1-j\omega_0)$,_where $\omega_0=N^{-1}(1-\exp(-\omega_{\rm max}))$, and the constants entering the_coupling_by and $c_j=\sqrt{\xi_K\omega_0}~\omega_j$._The adiabatic eigenvalues and eigenvectors, respectively, are $$E_{1,2}=V_b\mp\sqrt{\Omega^2+\gamma^2(R)}\;,$$ $$\begin{aligned}
\vert_1;R\rangle&=& \frac{1}{\sqrt{2(1+G^2}}\left(\begin{array}{c} 1+G\\
1-G\end{array} \right)
\nonumber\\
\vert 2;R\rangle&=& \frac{1}{\sqrt{2(1+G^2}}\left(\begin{array}{c}-1+G\\
1+G\end{array}_\right)\;,\end{aligned}$$ where $$G(R)=\gamma^{-1}(R)\left[-\Omega+\sqrt{\Omega^2+\gamma^2(R)}\right]\;.$$_The_coupling_vector $d_{\alpha\alpha'}=\langle\alpha;R\vert\overrightarrow{\partial}
/\partial R\vert\alpha';R\rangle$ is_$$d_{12}=-d_{21}=(1+G^2)^{-1}\partial G/\partial R\;.$$ Assuming an initially_uncorrelated density matrix, where the bath_is in thermal equilibrium and the subsystem_is in state $\vert\uparrow\rangle$, the initial_quantum-classical density matrix