right)},\label{eq:lax_vector}$$ the frequency spectrum is related to the roots of the equation $\sqrt{L^{2}\left(u\right)}=0$. We analyze them numerically in the limit $N\gg1$, and find that, in the high coupling regime $\Omega\gg\Delta$, all roots merge into a continuous line except two complex conjugated pairs, see example in Figure \[fig:L\_roots\]. Note that all coefficients of the polynomial $L^{2}\left(u\right)$ are real thus every complex root has a complex conjugated partner. The dynamical variables that correspond to the continuous band form a decay part of the solution and the two discreet frequencies give an oscillating part that we will be interested in. The discreet roots can be found by turning the summation over $j$ in Eq. (\[eq:lax\_vector\]) into an integral, $\sum_{j}\rightarrow\frac{N}{\Delta}\int_{\omega-\Delta/2}^{\omega+\Delta/2}d\epsilon$. Then, the equation $\sqrt{L^{2}\left(u\right)}=0$ turns into $$\pm2iga\left(0\right)=u-\frac{V^{2}N}{\Delta}\ln\left(\frac{u+\Delta/2}{u-\Delta/2}\right),\label{eq:lax2}$$ where the total width of splittings $\Delta=\epsilon_{N}-\epsilon_1$ is finite.

In the limit $M=N$, opposite to $M\ll N$, the roots of Eq. (\[eq:lax2\]) have zero real part. Parametrizing the roots as $u=iu_{0}\Delta/2$ we obtain

$$\frac{\pm4Va\left(0\right)}{\Delta}=u_{0}+\frac{2V^{2}N}{\Delta^{2}}\left(\pi-2\tan^{-1}u_{0}\right).$$

Then a $1/u_{0}$-expansion gives the imaginary parts as $u_{1,2}=u_{0}\Delta/2=\pm\left(\frac{V\sqrt{N}}{2}\pm\frac{\Delta}{2\sqrt{3}}\right)$.

In the case of two discreet roots the hyperelliptic function of many variables reduces to an elliptic function of only one variable [@relaxation] which corresponds to an effective model of a single
right)},\label{eq: lax_vector}$$ the frequency spectrum is related to the roots of the equation $ \sqrt{L^{2}\left(u\right)}=0$. We analyze them numerically in the limit $ N\gg1 $, and witness that, in the eminent coupling regime $ \Omega\gg\Delta$, all roots blend into a continuous line except two complex conjugated pair, experience example in Figure \[fig: L\_roots\ ]. notice that all coefficient of the polynomial $ L^{2}\left(u\right)$ are real thus every complex etymon has a complex conjugated partner. The dynamic variables that correspond to the continuous isthmus form a decay part of the solution and the two discerning frequencies give an oscillating part that we will be concerned in. The discreet roots can be find by turn the summation over $ j$ in Eq. (\[eq: lax\_vector\ ]) into an built-in, $ \sum_{j}\rightarrow\frac{N}{\Delta}\int_{\omega-\Delta/2}^{\omega+\Delta/2}d\epsilon$. Then, the equation $ \sqrt{L^{2}\left(u\right)}=0 $ turns into $ $ \pm2iga\left(0\right)=u-\frac{V^{2}N}{\Delta}\ln\left(\frac{u+\Delta/2}{u-\Delta/2}\right),\label{eq: lax2}$$ where the total width of splittings $ \Delta=\epsilon_{N}-\epsilon_1 $ is finite. 

 In the limit $ molarity = N$, opposite to $ M\ll N$, the roots of Eq. (\[eq: lax2\ ]) have zero real part. Parametrizing the etymon as $ u = iu_{0}\Delta/2 $ we obtain 

 $ $ \frac{\pm4Va\left(0\right)}{\Delta}=u_{0}+\frac{2V^{2}N}{\Delta^{2}}\left(\pi-2\tan^{-1}u_{0}\right).$$ 

 Then a $ 1 / u_{0}$-expansion gives the imaginary parts as $ u_{1,2}=u_{0}\Delta/2=\pm\left(\frac{V\sqrt{N}}{2}\pm\frac{\Delta}{2\sqrt{3}}\right)$. 

 In the case of two discreet roots the hyperelliptic function of many variables shrink to an elliptic function of only one variable star [ @relaxation ] which represent to an effective model of a single
rigjt)},\label{eq:lax_vector}$$ the fvequency spectrum is renated fo the ruots of the equation $\sqrt{L^{2}\leht(u\rught)}=0$. Qe analyze them numerizally in nhe limit $N\gg1$, qnd find tizt, in tmz higg cou'lmng regime $\Omegs\gg\Delta$, anl roots merge ivtl a continuous line except two compjex conkuhated pairs, seg exak[le jn Figure \[fig:L\_roots\]. Note that all doefficpents of the polymomial $L^{2}\left(u\right)$ are reap thks every complex rlot has a cimplqz conjugated partner. The dynamical variables that correspond to tfe coutinuous babd fogk a decay pert of the solutiok and tve two ciscreet frequcnciev guve an oscillating pact that we will be igterested iu. The discreet roots xab be xoung by rurvinf vhe summahioi over $j$ in Eq. (\[eq:lax\_vextor\]) into an integrsl, $\wum_{j}\rightarros\frac{N}{\Qejta}\int_{\omega-\Delta/2}^{\omega+\Delta/2}d\epsilon$. Then, ths equation $\sqrt{L^{2}\left(u\ritht)}=0$ turns into $$\pm2iga\lgft(0\right)=u-\fwac{V^{2}N}{\Delta}\ln\left(\frac{u+\Delta/2}{u-\Delta/2}\right),\label{eq:lax2}$$ wvere vhd titwu wldth of splittings $\Delta=\epsilon_{N}-\epsilon_1$ is figjtt.

In the limit $M=N$, jpposite to $M\pl G$, the roots ow Eq. (\[eq:las2\]) have zero real pwrt. Parwmetruzing the rooys as $u=iu_{0}\Delta/2$ we obtain

$$\frax{\pm4Va\left(0\rigkt)}{\Dwlta}=u_{0}+\frac{2V^{2}N}{\Delta^{2}}\leyt(\pi-2\tan^{-1}u_{0}\righc).$$

Then s $1/u_{0}$-exlansion gives the imagiuary pzrts as $u_{1,2}=u_{0}\Dflta/2=\pm\leff(\wrac{V\sqrt{N}}{2}\pm\frac{\Aelna}{2\sqst{3}}\right)$.

In the case of two qiscreet coots the hyoerekliptis function of many variables reducfs to at elliptic function of only one variable [@cxlaxation] whivh cogresponds to ak effective modql of a single
right)},\label{eq:lax_vector}$$ the frequency spectrum is related to of equation $\sqrt{L^{2}\left(u\right)}=0$. analyze them numerically find in the high regime $\Omega\gg\Delta$, all merge into a continuous line except complex conjugated pairs, see example in Figure \[fig:L\_roots\]. Note that all coefficients of polynomial $L^{2}\left(u\right)$ are real thus every complex root has a complex conjugated partner. dynamical that to continuous band form a decay part of the solution and the two discreet frequencies give an part that we will be interested in. The roots can be found turning the summation over $j$ Eq. into an $\sum_{j}\rightarrow\frac{N}{\Delta}\int_{\omega-\Delta/2}^{\omega+\Delta/2}d\epsilon$. the $\sqrt{L^{2}\left(u\right)}=0$ turns into where the total width of splittings $\Delta=\epsilon_{N}-\epsilon_1$ is finite. In the limit $M=N$, opposite to $M\ll N$, roots of have zero part. the as $u=iu_{0}\Delta/2$ we Then a $1/u_{0}$-expansion gives the imaginary In the case of two discreet roots the function of variables reduces to an elliptic function only one variable [@relaxation] which corresponds to an model of a single
right)},\label{eq:lax_vector}$$ the fRequency spEctruM is RelAtEd to The rOots of the equatIOn $\sqRt{L^{2}\left(u\right)}=0$. We analyze Them nUmERicaLLy In the Limit $N\gG1$, AnD FInd ThAt, In tHe HIgH coupLinG regime $\omega\gg\DelTa$, aLl Roots merge inTO a Continuous LinE except two coMplEx conjUgAteD Pairs, See ExampLe in FiGUre \[fig:l\_roots\]. NotE tHAt all cOEfficieNTS oF the Polynomial $L^{2}\left(u\rIGhT)$ Are real thus eveRy compLeX RoOT Has A coMplex conjuGaTed paRTner. The DYnAMICal VAriables that cOrrespond to THe cOntinuOuS baND form a Decay PaRT of The solution And tHe two discReet frEQuencieS Give an oScillaTinG paRt thAT wE wIll Be INteREsTed IN. ThE discreeT rOoTs can Be foUND BY turNinG the SummaTion over $j$ in Eq. (\[Eq:lAx\_veCTor\]) Into aN inteGral, $\SuM_{j}\rigHtarroW\frac{n}{\DElta}\int_{\omega-\DelTa/2}^{\omEga+\Delta/2}d\EpsIlOn$. THeN, the eQUation $\SqrT{L^{2}\lEft(u\rigHt)}=0$ turns INto $$\Pm2IGA\LeFt(0\right)=u-\frac{V^{2}N}{\DeltA}\lN\LEfT(\frac{u+\DeLta/2}{u-\DeLTa/2}\RiGHt),\label{eQ:lAx2}$$ wHere THE totaL widTH oF splittiNgs $\DelTA=\ePsIlon_{N}-\epSiLon_1$ is fInIte.

in tHe limIT $M=N$, oPpositE to $M\ll N$, tHe rooTS of Eq. (\[eq:lax2\]) have ZEro real part. PaRAmETRiZIng tHe rOots as $u=iu_{0}\DeLta/2$ wE ObtaIn

$$\frAC{\pM4Va\LEft(0\riGht)}{\DeLtA}=U_{0}+\fRAc{2V^{2}N}{\Delta^{2}}\left(\pi-2\tan^{-1}u_{0}\RiGht).$$

TheN a $1/u_{0}$-exPansion gives tHe imaginarY PARts as $u_{1,2}=u_{0}\DElta/2=\PM\lEFt(\frac{V\sqrt{N}}{2}\pm\Frac{\DElta}{2\sqrt{3}}\riGHt)$.

In the cAse of Two discrEet roots tHE HyperellIptIc fUncTioN OF mAny variables rEDUces To An ellipTic FunctioN of OnlY onE vaRiAble [@relaxAtion] whiCh CoRrEsPonDs to aN EffectivE mOdeL oF a sIngle
right)},\label{eq:lax_vect or}$$ thefrequ enc y s pe ctru m is  related to th e  roo ts of the equation $\s qrt{L ^{ 2 }\le f t( u\rig ht)}=0$ .  W e ana ly ze  th em nu meric all y in th e limit $N \gg 1$ , and find t h at , in the h igh  coupling re gim e $\Om eg a\g g \Delt a$,  allrootsm erge i nto a con ti n uous l i ne exce p t  t wo c omplex conjugated pa i rs, see exampl e in F ig u re \ [fi g:L \_roots\].  N ote t h at allc oe f f i cie n ts of the pol ynomial $L^ { 2}\ left(u \r igh t )$ are  real  t h usevery compl ex r oot has a  compl e x conju g ated pa rtner.  Th e d ynam i ca lvar ia b les th atc orr espond t oth e con tinu o u s band  fo rm a  deca y part of the  so luti o n a nd th e two  dis cr eet f requen ciesgi ve an oscillati ng p art thatwewi llbe  inte r estedin.  Th e discr eet roo t s c an b e  f ound by turning th es u mm ation ov er $j$ in  E q . (\[eq: la x\_ vect o r \]) i ntoa nintegral , $\su m _{ j} \righta rr ow\fra c{ N}{ \De lta}\ i nt_{ \omega -\Delta/ 2}^{\ o mega+\Delta/2} d \epsilon$. Th e n, t he equa tio n $\sqrt{L^ {2}\ l eft( u\ri g ht )}= 0 $ tur ns in to $$ \ pm2iga\left(0\right )= u-\fra c{V^{ 2}N}{\Delta}\ ln\left(\f r a c {u+\Delt a/2} { u- \ Delta/2}\right ),\la bel{eq:lax 2 }$$ wher e the  total w idth of s p l ittings$\D elt a=\ eps i l on _{N}-\epsilon _ 1 $ is  f inite.
In  the li mit  $M =N$ , o pp osite to$M\ll N$ ,th ero ots  of E q . (\[eq: la x2\ ])  ha ve ze r o real  part . Pa ra me t riz ing the ro o t s as  $ u= iu_{ 0}\ De lta/2 $ we obt ain

$$ \frac{\pm 4Va \ left (0 \r ight)}{ \Delta}=u_{0} +\ frac{2V^{2 }N }{\ Delta^ { 2 }}\left( \pi-2\tan^{-1}u_{0}\rig h t).$$
The n a $ 1/u_ {0}$-expa nsi on giv est he ima ginary  part sas$ u _{1,2 } = u_ {0} \D elta/2=\pm \ l eft (\fra c{ V\sq rt{N}}{ 2}\pm\frac{\Delta} { 2\s qrt{3}}\right )$. 

In t he  ca s eo f t wo dis c r eet roots the h yperellipt ic fu nction ofm any  v ariable s reduc es to an elli ptic func tion of o nl y on e var iable [@re laxation ] which c o rresp o nd s toaneffect iv e m odelof a s i ngl e
right)},\label{eq:lax_vector}$$ the_frequency spectrum_is related to the_roots of_the_equation $\sqrt{L^{2}\left(u\right)}=0$._We_analyze them numerically_in the limit_$N\gg1$, and find that,_in the high_coupling_regime $\Omega\gg\Delta$, all roots merge into a continuous line except two complex conjugated pairs,_see_example in_Figure_\[fig:L\_roots\]._Note that all coefficients of_the polynomial $L^{2}\left(u\right)$ are real_thus every_complex root has a complex conjugated partner. The_dynamical_variables that correspond_to the continuous band form a decay part of_the solution and the two discreet_frequencies give an_oscillating_part_that we will be_interested in. The discreet roots can_be found by turning the summation_over $j$ in Eq. (\[eq:lax\_vector\]) into an_integral, $\sum_{j}\rightarrow\frac{N}{\Delta}\int_{\omega-\Delta/2}^{\omega+\Delta/2}d\epsilon$. Then, the equation $\sqrt{L^{2}\left(u\right)}=0$_turns into $$\pm2iga\left(0\right)=u-\frac{V^{2}N}{\Delta}\ln\left(\frac{u+\Delta/2}{u-\Delta/2}\right),\label{eq:lax2}$$ where the_total width_of splittings $\Delta=\epsilon_{N}-\epsilon_1$ is finite.

In_the limit $M=N$,_opposite to_$M\ll N$, the_roots of Eq. (\[eq:lax2\]) have zero_real part. Parametrizing_the roots as $u=iu_{0}\Delta/2$ we obtain

$$\frac{\pm4Va\left(0\right)}{\Delta}=u_{0}+\frac{2V^{2}N}{\Delta^{2}}\left(\pi-2\tan^{-1}u_{0}\right).$$

Then_a_$1/u_{0}$-expansion gives the_imaginary_parts_as $u_{1,2}=u_{0}\Delta/2=\pm\left(\frac{V\sqrt{N}}{2}\pm\frac{\Delta}{2\sqrt{3}}\right)$.

In_the case of_two_discreet roots_the_hyperelliptic function of many variables reduces_to_an elliptic function of only one variable_[@relaxation] which corresponds to_an_effective model of a_single