 of the initial state ${| M;\down\ldots\down \rangle}$. The initial condition for the dynamical mean-field equations, which corresponds to this state, is all spins down, $\mathbf{C}_{j}\left(0\right)=\left(0,0,-1/2\right)$, and a finite amplitude of the bosonic field, $a\left(0\right)=\sqrt{M}$. Below, we solve set of differential Eqs. (\[eq:CEOM1\], \[eq:CEOM2\]) for the evolution of classical variables with this initial condition.

For a small number of excitations in a system with many spins, $M\ll N$, these equations were solved in \[\]. Using the approximation $C_{j}^{z}\left(t\right)\approx-1/2$, the equation for $C_{j}^{z}\left(t\right)$ drops out from Eq. (\[eq:CEOM1\]) and the remaining system of equation is harmonic. In the regime $\Omega\gg\Delta$ this gives the following solution for the bosonic mode $$a\left(t\right)=\sqrt{M}\cos\left(V\sqrt{N}t\right).
\label{eq:harmonic}$$ The period of the oscillatory function in this limit is $T=2\pi/\left(V\sqrt{N}\right)$, see Figure \[fig:time\_evolution\] for $M=10$. In the following will use the Rabi frequency $\Omega=V\sqrt{N}$ which was introduced in Sec. II.

For arbitrary $M$ the solution to Eq. (\[eq:CEOM1\], \[eq:CEOM2\]) are hyperelliptic functions [@fermi_bose]. Here we will not consider the full analytic form of general solution but will use only the spectral analysis developed in Ref. \[\]. Introducing the following vector function (Lax vector) of an auxiliary parameter $u$, $$\mathbf{L}\left(u\right)=\left(\begin{array}{c}
\frac{a_{x}\left(0\right)}{V}\\
\frac{a_{y}\left(0\right)}{V}\\
\frac{u}{2V^{2}}
\end{array}\right)+\sum_{j}\frac{\mathbf{C}_{j}\left(0\right)}{u-\left(\epsilon_{j}-\omega\
of the initial state $ { | M;\down\ldots\down \rangle}$. The initial condition for the dynamical mean - sphere equation, which corresponds to this state, is all spins down, $ \mathbf{C}_{j}\left(0\right)=\left(0,0,-1/2\right)$, and a finite amplitude of the bosonic discipline, $ a\left(0\right)=\sqrt{M}$. Below, we solve set up of differential Eqs. (\[eq: CEOM1\ ], \[eq: CEOM2\ ]) for the development of classical variables with this initial condition. 

 For a small number of excitation in a system with many spins, $ M\ll N$, these equation were solved in \[\ ]. Using the approximation $ C_{j}^{z}\left(t\right)\approx-1/2 $, the equation for $ C_{j}^{z}\left(t\right)$ drops out from Eq. (\[eq: CEOM1\ ]) and the stay system of equation is harmonic. In the regimen $ \Omega\gg\Delta$ this gives the following solution for the bosonic mode $ $ a\left(t\right)=\sqrt{M}\cos\left(V\sqrt{N}t\right). 
 \label{eq: harmonic}$$ The period of the oscillatory affair in this limit is $ T=2\pi/\left(V\sqrt{N}\right)$, see Figure \[fig: time\_evolution\ ] for $ M=10$. In the following will use the Rabi frequency $ \Omega = V\sqrt{N}$ which was introduced in Sec. II. 

 For arbitrary $ M$ the solution to Eq. (\[eq: CEOM1\ ], \[eq: CEOM2\ ]) are hyperelliptic routine [ @fermi_bose ]. Here we will not consider the full analytic form of general solution but will use only the spectral analysis developed in Ref. \[\ ]. Introducing the following vector function (Lax vector) of an auxiliary parameter $ u$, $ $ \mathbf{L}\left(u\right)=\left(\begin{array}{c } 
 \frac{a_{x}\left(0\right)}{V}\\ 
 \frac{a_{y}\left(0\right)}{V}\\ 
 \frac{u}{2V^{2 } } 
 \end{array}\right)+\sum_{j}\frac{\mathbf{C}_{j}\left(0\right)}{u-\left(\epsilon_{j}-\omega\
 of the initial state ${| M;\dowk\ldots\down \ranglg}$. Rhe inmtial ckndition for the dynamical mean-field ewuatiins, which corresponds go this snate, is aol skins down, $\mathbf{C}_{o}\meft(0\rigmc)=\left(0,0,-1/2\dlght)$, cnv a finite amplltude of tha bosonic fielg, $x\lzft(0\right)=\sqrt{M}$. Below, we solve set of qifferemtlal Eqs. (\[eq:CEOM1\], \[eq:CTOM2\]) for nht evolution of classical variablea with uhis initial condiyion.

For a small number of fxcihations in a systel with many spigw, $M\ll N$, thesd equations were solves in \[\]. Using the approximation $C_{g}^{z}\lefc(t\right)\apprix-1/2$, thf equation fir $C_{j}^{d}\left(t\right)$ dvpps oud from Rq. (\[eq:CEOM1\]) and bhe rxmaibing system of equatimn is harmonic. In jhe regime $\Ojega\gg\Delta$ this tices tve fmllodung somuvioh for hhe bosonic mkde $$a\left(t\rught)=\sqrt{M}\cos\left(V\sqtt{G}n\tight).
\label{eq:garmonyc}$$ The period of the oscillatory function in this limit is $T=2\pi/\left(V\wqrt{N}\right)$, see Figure \[fig:time\_edolution\] for $M=10$. In the following will use the Rabi freqnevcy $\Oidta=G\sqrt{N}$ which was introduced in Sec. II.

For arbiedaty $M$ the solutiok to Eq. (\[eq:CEOM1\], \[eq:CROL2\]) ste hyperelliptkc funefikns [@fermi_bose]. Here we wilj not consider the full analytic form of geneeal solution vut will use only che spectral anakysis developed in Ref. \[\]. Intruducjng the folpowing vedgor function (Lax vebtor) of an auxiliary parameter $u$, $$\mathbf{O}\lefc(u\right)=\ldft(\bggin{arrwy}{c}
\frac{a_{x}\lfft(0\rinvt)}{V}\\
\frac{a_{y}\left(0\right)}{G}\\
\frac{o}{2V^{2}}
\end{asray}\right)+\skm_{j}\frac{\mathbf{C}_{j}\left(0\right)}{u-\left(\epsilon_{j}-\omega\
of the initial state ${| M;\down\ldots\down \rangle}$. condition the dynamical equations, which corresponds spins $\mathbf{C}_{j}\left(0\right)=\left(0,0,-1/2\right)$, and a amplitude of the field, $a\left(0\right)=\sqrt{M}$. Below, we solve set differential Eqs. (\[eq:CEOM1\], \[eq:CEOM2\]) for the evolution of classical variables with this initial For a small number of excitations in a system with many spins, $M\ll these were in Using the approximation $C_{j}^{z}\left(t\right)\approx-1/2$, the equation for $C_{j}^{z}\left(t\right)$ drops out from Eq. (\[eq:CEOM1\]) and the remaining of equation is harmonic. In the regime $\Omega\gg\Delta$ gives the following solution the bosonic mode $$a\left(t\right)=\sqrt{M}\cos\left(V\sqrt{N}t\right). \label{eq:harmonic}$$ period the oscillatory in limit $T=2\pi/\left(V\sqrt{N}\right)$, see Figure for $M=10$. In the following will use the Rabi frequency $\Omega=V\sqrt{N}$ which was introduced in Sec. II. arbitrary $M$ to Eq. \[eq:CEOM2\]) hyperelliptic [@fermi_bose]. Here we consider the full analytic form of will use only the spectral analysis developed in \[\]. Introducing following vector function (Lax vector) of auxiliary parameter $u$, $$\mathbf{L}\left(u\right)=\left(\begin{array}{c} \frac{a_{x}\left(0\right)}{V}\\ \frac{a_{y}\left(0\right)}{V}\\ \frac{u}{2V^{2}} \end{array}\right)+\sum_{j}\frac{\mathbf{C}_{j}\left(0\right)}{u-\left(\epsilon_{j}-\omega\
 of the initial state ${| M;\down\ldoTs\down \rangLe}$. The IniTiaL cOndiTion For the dynamicaL Mean-Field equations, which corRespoNdS To thIS sTate, iS all spiNS dOWN, $\maThBf{c}_{j}\lEfT(0\RiGht)=\leFt(0,0,-1/2\rIght)$, and A finite ampLitUdE of the bosoniC FiEld, $a\left(0\riGht)=\Sqrt{M}$. Below, we SolVe set oF dIffERentiAl EQs. (\[eq:CeOM1\], \[eq:CeoM2\]) for tHe evolutiOn OF classICal variABLeS witH this initial condiTIoN.

for a small numbeR of excItATiONS in A syStem with maNy Spins, $m\Ll N$, thesE EqUATIonS Were solved in \[\]. USing the apprOXimAtion $C_{J}^{z}\LefT(T\right)\ApproX-1/2$, tHE eqUation for $C_{j}^{Z}\lefT(t\right)$ drOps out FRom Eq. (\[eq:ceOM1\]) and tHe remaIniNg sYsteM Of EqUatIoN Is hARmOniC. in tHe regime $\omEgA\gg\DeLta$ tHIS GIves The FollOwing Solution for thE boSoniC ModE $$a\lefT(t\rigHt)=\sqRt{m}\cos\lEft(V\sqRt{N}t\rIgHt).
\label{eq:harmonIc}$$ ThE period of The OsCilLaTory fUNction In tHis Limit is $t=2\pi/\left(v\SqrT{N}\RIGHt)$, See Figure \[fig:time\_evOlUTIoN\] for $M=10$. In tHe follOWiNg WIll use thE RAbi FreqUENcy $\OmEga=V\SQrT{N}$ which wAs intrODuCeD in Sec. Ii.

FOr arbiTrAry $m$ thE soluTIon tO Eq. (\[eq:CeOM1\], \[eq:CEOm2\]) are hYPerelliptic funCTions [@fermi_bosE]. heRE We WIll nOt cOnsider the fUll aNAlytIc foRM oF geNEral sOlutiOn BUt WIll use only the spectrAl AnalysIs devEloped in Ref. \[\]. InTroducing tHE FOllowing VectOR fUNction (Lax vectoR) of an Auxiliary pARameter $u$, $$\MathbF{L}\left(u\rIght)=\left(\bEGIn{array}{c}
\FraC{a_{x}\LefT(0\riGHT)}{V}\\
\Frac{a_{y}\left(0\rigHT)}{v}\\
\fraC{u}{2v^{2}}
\end{arrAy}\rIght)+\sum_{J}\frAc{\mAthBf{C}_{J}\lEft(0\right)}{u-\Left(\epsiLoN_{j}-\OmEgA\
 of the initial state ${|M;\down\ld ots\d own  \r an gle} $. T he initial con d itio n for the dynamical me an-fi el d  equ a ti ons,which c o rr e s pon ds  t o t hi s  s tate,  is  all sp ins down,$\m at hbf{C}_{j}\l e ft (0\right)= \le ft(0,0,-1/2\ rig ht)$,an d a finit e a mplit ude of the bo sonic fie ld ,  $a\le f t(0\rig h t )= \sqr t{M}$. Below, wes ol v e set of diffe rentia lE qs . (\[ eq: CEOM1\], \ [e q:CEO M 2\]) fo r  t h e evo l ution of clas sical varia b les  withth isi nitial  cond it i on. 

For a sma ll n umber ofexcita t ions in a syste m with  ma nyspin s ,$M \ll  N $ , t h es e e q uat ions wer eso lvedin \ [ \ ] .  Usi ngtheappro ximation $C_{ j}^ {z}\ l eft (t\ri ght)\ appr ox -1/2$ , theequat io n for $C_{j}^{z }\le ft(t\righ t)$  d rop sout f r om Eq.  (\ [eq :CEOM1\ ]) andt here m a i ni ng system of equat io n is  harmoni c. Int he  r e gime $\O me ga\ gg\D e l ta$ t hisg iv es the f ollowi n gso lutionfo r thebo son icmode$ $a\l eft(t\ right)=\ sqrt{ M }\cos\left(V\s q rt{N}t\right) . 
\ l a be l {eq: har monic}$$ Th e pe r iodof t h eosc i llato ry fu nc t io n  in this limit is $ T= 2\pi/\ left( V\sqrt{N}\rig ht)$, seeF i g ure \[fi g:ti m e\ _ evolution\] fo r $M= 10$. In th e  followi ng wi ll use t he Rabi f r e quency $ \Om ega =V\ sqr t { N} $ which was i n t rodu ce d in Se c.II.

Fo r a rbi tra ry$M $ the sol ution to  E q.  ( \[ eq: CEOM1 \ ], \[eq: CE OM2 \] ) a re hy p erelli pticfunc ti on s  [@ fermi_b o se ] .  Her ewe  wil l n ot  cons ider the  full a nalytic f orm of g en er al solu tion but will  u se only th espe ctrala n alysis d eveloped in Ref. \[\].I ntroduc ing  thefoll owing vec tor  funct ion (Lax v ector)  of a naux i l iaryp a ra met er  $u$, $$\m a t hbf {L}\l ef t(u\ right)= \left(\begin{array } {c} 
\frac{a_{x}\ lef t(0\ r i gh t)} { V} \ \
\ fr a c{a _ { y}\left(0\right )}{V}\\
\f ra c {u }{2V^{2}}\ end {a rray}\r ight)+\ sum_{ j }\frac{ \mathbf{C }_{j}\lef t( 0\ri g h t)} {u-\left(\ epsilon_ {j}-\omeg a \
 of_the initial_state ${| M;\down\ldots\down \rangle}$._The initial_condition_for the_dynamical_mean-field equations, which_corresponds to this_state, is all spins_down, $\mathbf{C}_{j}\left(0\right)=\left(0,0,-1/2\right)$, and_a_finite amplitude of the bosonic field, $a\left(0\right)=\sqrt{M}$. Below, we solve set of differential Eqs._(\[eq:CEOM1\],_\[eq:CEOM2\]) for_the_evolution_of classical variables with this_initial condition.

For a small number_of excitations_in a system with many spins, $M\ll N$,_these_equations were solved_in \[\]. Using the approximation $C_{j}^{z}\left(t\right)\approx-1/2$, the equation for_$C_{j}^{z}\left(t\right)$ drops out from Eq. (\[eq:CEOM1\])_and the remaining_system_of_equation is harmonic. In_the regime $\Omega\gg\Delta$ this gives the_following solution for the bosonic mode_$$a\left(t\right)=\sqrt{M}\cos\left(V\sqrt{N}t\right).
\label{eq:harmonic}$$ The period of the oscillatory function_in this limit is $T=2\pi/\left(V\sqrt{N}\right)$, see_Figure \[fig:time\_evolution\] for $M=10$. In_the following_will use the Rabi frequency_$\Omega=V\sqrt{N}$ which was_introduced in_Sec. II.

For arbitrary_$M$ the solution to Eq. (\[eq:CEOM1\],_\[eq:CEOM2\]) are hyperelliptic_functions [@fermi_bose]. Here we will not_consider_the full analytic_form_of_general solution_but will use_only_the spectral_analysis_developed in Ref. \[\]. Introducing the_following_vector function (Lax vector) of an auxiliary_parameter $u$, $$\mathbf{L}\left(u\right)=\left(\begin{array}{c}
\frac{a_{x}\left(0\right)}{V}\\
\frac{a_{y}\left(0\right)}{V}\\
\frac{u}{2V^{2}}
\end{array}\right)+\sum_{j}\frac{\mathbf{C}_{j}\left(0\right)}{u-\left(\epsilon_{j}-\omega\