.  \label{Veff}$$ Here $\hat U_{ol} \equiv A \cos(2 x) $ is the OL and $\hat U_s$ is the potential associated with a given eigenstate of the quantum problem (\[schro\]). For a self-consistent solution, one starts with a trial wavefunction for $\psi_s$ (typically a gaussian waveform), calculates the effective potential and solves the corresponding eigenvalue problem (\[schro\]). Then, one selects a given eigenstate (for example the ground state but not necessarily) as new trial function and iterates the procedure until convergence is reached.

![ Panel [**(a)**]{} Energy spectrum for the effective potential (\[Veff\]) with $A=3$ and $\chi=0$ (Mathieu equation). Full lines represent exact values of the band edges of the Mathieu equation while dots are the eigenvalues obtained with the above procedure on a lattice of length $L=40 \pi$, with $N=512$ points. Panel [**(b)**]{} Lowest energy band for the effective potential in Eq. (\[Veff\]) with $\psi_s$ taken as the ground state of the system and for $\chi=-1$ (attractive case). Parameters are fixed as in panel (a). Panel [**(c)**]{} The same as in panel (b) but for $A=-3$. Panel [**(d)**]{} Transition from the metastable IS mode to the OS ground state corresponding to the lower level of panel (c). The optical lattice (scaled by a factor 3) is reported as an help to locate the symmetry center of the solutions. Parameters are fixed as in panel (c).[]{data-label="fig1"}](fig1a.eps "fig:"){width="3.8cm" height="3.3cm"}![ Panel [**(a)**]{} Energy spectrum for the effective potential (\[Veff\]) with $A=3$ and $\chi=0$ (Mathieu equation). Full lines represent exact values of the band edges of the Mathieu equation while dots are the eigenvalues obtained with the above procedure on a lattice of length $L=40 \pi$, with $N=512$ points. Panel [**(b)**]{} Lowest energy band for the effective potential in Eq. (\[Veff\]) with $\psi_s$ taken as the ground state of the system and for $\chi=-1$ (attractive case). Parameters are fixed as in panel (a). Panel [**(c
.   \label{Veff}$$ Here $ \hat U_{ol } \equiv A \cos(2 x) $ is the OL and $ \hat U_s$ is the potential associated with a given eigenstate of the quantum problem (\[schro\ ]). For a self - reproducible solution, one start with a trial wavefunction for $ \psi_s$ (typically a gaussian waveform), calculates the effective electric potential and solves the corresponding eigenvalue problem (\[schro\ ]). Then, one choose a impart eigenstate (for example the land state but not necessarily) as modern trial function and iterate the operation until convergence is reached. 

! [ Panel [ * * (a) * * ] { } Energy spectrum for the effective potential (\[Veff\ ]) with $ A=3 $ and $ \chi=0 $ (Mathieu equation). Full credit line represent exact values of the ring edges of the Mathieu equation while dots are the eigenvalues obtain with the above procedure on a lattice of length $ L=40 \pi$, with $ N=512 $ points. Panel [ * * (b) * * ] { } low energy band for the effective potential in Eq. (\[Veff\ ]) with $ \psi_s$ taken as the ground department of state of the system and for $ \chi=-1 $ (attractive case). Parameters are fixed as in panel (a). Panel [ * * (c) * * ] { } The same as in panel (b) but for $ A=-3$. Panel [ * * (d) * * ] { } Transition from the metastable IS manner to the OS ground state represent to the low level of panel (c). The ocular lattice (scaled by a factor 3) is reported as an aid to locate the symmetry center of the solutions. Parameters are fixed as in gore (c).[]{data - label="fig1"}](fig1a.eps " fig:"){width="3.8 cm " height="3.3 cm " }! [ Panel [ * * (a) * * ] { } Energy spectrum for the effective potential (\[Veff\ ]) with $ A=3 $ and $ \chi=0 $ (Mathieu equation). Full telephone line represent accurate value of the band edges of the Mathieu equation while dots are the eigenvalues obtained with the above procedure on a wicket of length $ L=40 \pi$, with $ N=512 $ point. Panel [ * * (b) * * ] { } Lowest department of energy band for the effective potential in Eq. (\[Veff\ ]) with $ \psi_s$ taken as the ground state of the arrangement and for $ \chi=-1 $ (attractive case). parameter are fixed as in jury (a). Panel [ * * (c
.  \lwbel{Veff}$$ Here $\hat U_{ol} \eqmiv A \cos(2 x) $ is jhw OL aid $\hat H_s$ is thd potential associated with e gicen eugenstate of the quantjm problel (\[schro\]). Dor e self-consistent solutiok, one dtarcs with a trial eavefunctimn for $\psi_s$ (ty[izaply a gaussian waveform), calculates ehe effrchive potential and folvss the corresponding eigenvalue prkblem (\[sbhro\]). Then, one selrcts a given eigenstate (fog exwmple the ground shate but noj neswssarily) as vew trial function and iterates the procedure until cunveryence is reqcyed.

![ Kanel [**(a)**]{} Enerjy spebtrum for the effectiee potemtial (\[Veff\]) witm $A=3$ aid $\cyi=0$ (Mathieu equation). Fnll lines represent gxact valuas of the band edgew if thg Matvieu wquxtikn wgile dlts are the ejgenvalues ibtained with the anodv procedure oh a laetyce of length $L=40 \pi$, with $N=512$ points. Panel [**(t)**]{} Lkwest energy band for tye effective potentiap in Eq. (\[Vqff\]) with $\psi_s$ taken as the ground state of the syvtem ena fir $\zyi=-1$ (attractive case). Parameters are fixed as in pwhek (s). Panel [**(c)**]{} The fame as in lajek (b) but for $A=-3$. Kanel [**(d)**]{} Trznsition from the letastafle IW mode to the OS ground state correspondung to the ljqer level of panel (c). The optieal lajtice (xcaled by a factor 3) is xeportsd as an hepp to loczge the symmetry zenner mf the solutions. Parameterf are fixxd as in pandl (c).[]{cata-lafel="fig1"}](fig1a.fps "flc:"){width="3.8cm" height="3.3cm"}![ Oanel [**(a)**]{} Energy spfctrum for the effective potentmel (\[Veff\]) with $S=3$ dnd $\chi=0$ (Matkieu eauation). Full jines represenj exact vclues uf the bans edges of the Matrieu equation ahile dots ace the eidenvqluew obtaivdd with the abpve procebbre on a oattice of length $K=40 \pk$, with $N=512$ points. 'cntl [**(b)**]{} Lowest enerby canq voc the affective podentkal on Eq. (\[Veff\]) with $\isi_r$ talen as the ground stdte kf the system and gov $\chi=-1$ (atttactive cwse). Parameterx are fixed as in kanel (e). Paneu [**(c
. \label{Veff}$$ Here $\hat U_{ol} \equiv A $ the OL $\hat U_s$ is given of the quantum (\[schro\]). For a solution, one starts with a trial for $\psi_s$ (typically a gaussian waveform), calculates the effective potential and solves the eigenvalue problem (\[schro\]). Then, one selects a given eigenstate (for example the ground but necessarily) new function and iterates the procedure until convergence is reached. ![ Panel [**(a)**]{} Energy spectrum for the potential (\[Veff\]) with $A=3$ and $\chi=0$ (Mathieu equation). lines represent exact values the band edges of the equation dots are eigenvalues with above procedure on lattice of length $L=40 \pi$, with $N=512$ points. Panel [**(b)**]{} Lowest energy band for the effective potential Eq. (\[Veff\]) taken as ground of system and for case). Parameters are fixed as in [**(c)**]{} The same as in panel (b) but $A=-3$. Panel Transition from the metastable IS mode the OS ground state corresponding to the lower of panel (c). The optical lattice (scaled by a factor 3) is reported as an locate the symmetry center the solutions. Parameters fixed in (c).[]{data-label="fig1"}](fig1a.eps height="3.3cm"}![ Panel Energy spectrum for the effective potential (\[Veff\]) with $A=3$ and $\chi=0$ equation). Full lines represent exact values of the band edges Mathieu while dots are eigenvalues obtained with the procedure a lattice of length with points. energy for effective potential in Eq. with $\psi_s$ taken as the state of the system Parameters are fixed as in panel (a). Panel
.  \label{Veff}$$ Here $\hat U_{ol} \equiv A \Cos(2 x) $ is the Ol and $\hAt U_S$ is ThE potEntiAl associated wiTH a giVen eigenstate of the quanTum prObLEm (\[scHRo\]). for a sElf-consIStENT soLuTiOn, oNe STaRts wiTh a Trial waVefunction For $\PsI_s$ (typically a GAuSsian wavefOrm), Calculates thE efFectivE pOteNTial aNd sOlves The corREspondIng eigenvAlUE problEM (\[schro\]). THEN, oNe seLects a given eigensTAtE (For example the gRound sTaTE bUT Not NecEssarily) as NeW triaL FunctioN AnD ITEraTEs the procedurE until conveRGenCe is reAcHed.

![ pAnel [**(a)**]{} ENergy SpECtrUm for the effEctiVe potentiAl (\[Veff\]) WIth $A=3$ and $\CHi=0$ (MathiEu equaTioN). FuLl liNEs RePreSeNT exACt ValUEs oF the band EdGeS of thE MatHIEU EquaTioN whiLe dotS are the eigenvAluEs obTAinEd witH the aBove PrOceduRe on a lAtticE oF length $L=40 \pi$, with $N=512$ PoinTs. Panel [**(b)**]{} LOweSt EneRgY band FOr the eFfeCtiVe potenTial in EQ. (\[vefF\]) wITH $\PsI_s$ taken as the ground StATE oF the systEm and fOR $\cHi=-1$ (ATtractivE cAse). paraMETers aRe fiXEd As in paneL (a). PaneL [**(C)**]{} THe Same as iN pAnel (b) bUt For $a=-3$. PaNel [**(d)**]{} TRAnsiTion frOm the metAstabLE IS mode to the OS GRound state corREsPONdINg to The Lower level oF panEL (c). ThE optICaL laTTice (sCaled By A FaCTor 3) is reported as an heLp To locaTe the Symmetry centeR of the soluTIONs. ParameTers ARe FIxed as in panel (c).[]{Data-lAbel="fig1"}](fig1A.Eps "fig:"){wiDth="3.8cm" Height="3.3cm"}![ panel [**(a)**]{} EneRGY spectruM foR thE efFecTIVe Potential (\[Veff\]) WITh $A=3$ aNd $\Chi=0$ (MathIeu EquatioN). FuLl lIneS rePrEsent exacT values oF tHe BaNd EdgEs of tHE Mathieu EqUatIoN whIle doTS are thE eigeNvalUeS oBTaiNed with THe ABOve pRoCeDure On a LaTtice Of leNGth $l=40 \pi$, with $n=512$ points. PaNel [**(B)**]{} loweSt EnErgy banD for the effectIvE potential In eq. (\[VEff\]) witH $\PSi_s$ taken As the ground state of the sySTem and fOr $\cHi=-1$ (attRactIve case). PaRamEters aRe fIXed as iN panel (A). PaneL [**(c
.  \label{Veff}$$ Here $\h at U_{ol}\equi v A  \c os (2 x ) $is the OL and$ \hat  U_s$ is the potential  asso ci a tedw it h a g iven ei g en s t ate  o fthe  q u an tum p rob lem (\[ schro\]).For  a  self-consis t en t solution , o ne starts wi tha tria lwav e funct ion  for$\psi_ s $ (typ ically aga u ssianw aveform ) ,  c alcu lates the effecti v ep otential and s olvesth e  c o r res pon ding eigen va lue p r oblem ( \ [s c h r o\] ) . Then, one s elects a gi v eneigens ta te( for ex ample  t h e g round state  but  not nece ssaril y ) as ne w  trialfuncti onand  ite r at es  th ep roc e du reu nti l conver ge nc e isreac h e d . 

![  Pa nel[**(a )**]{} Energy  sp ectr u m f or th e eff ecti ve  pote ntial(\[Ve ff \]) with $A=3$and$\chi=0$(Ma th ieu  e quati o n). Fu lllin es repr esent e x act  v a l u es  of the band edges  o f th e Mathie u equa t io nw hile dot sare  the e igenv alue s  o btainedwith t h eab ove pro ce dure o na l att ice o f  len gth $L =40 \pi$ , wit h  $N=512$ point s . Panel [**(b ) ** ] { }L owes t e nergy bandfort he e ffec t iv e p o tenti al in  E q .( \[Veff\]) with $\ps i_ s$ tak en as  the ground s tate of th e s ystem an d fo r  $ \ chi=-1$ (attra ctive  case). Pa r ametersare f ixed asin panel( a ). Panel  [* *(c )** ]{} T he  same as in p a n el ( b)  but fo r $ A=-3$.Pan el[** (d) ** ]{} Trans ition fr om  t he  m eta stabl e  IS mode  t o t he  OS  grou n d stat e cor resp on di n g t o the l o we r leve lof  pan el(c ). Th e op t ica l latti ce (scale d b y  a f ac to r 3) is  reported asan  help to l oc ate  the s y m metry ce nter of the solutions.P aramete rsare f ixed  as in pa nel  (c).[ ]{d a ta-lab el="fi g1"}] (f ig1 a . eps " f i g: "){ wi dth="3.8cm " hei ght=" 3. 3cm" }![ Pan el [**(a)**]{} Ene r gyspectrum forthe  eff e c ti vep ot e nti al (\[ V e ff\]) with $A=3 $ and $\ch i= 0 $(Mathieu e q uat io n). Ful l lines  repr e sent ex act value s of theba nd e d g esof the Mat hieu equ ation whi l e dot s  a re th e e igenva lu esobtai ned wi t h t he ab ove pr oc edureon ala ttice of  length $L=40 \pi$, wit h $N=5 12$ p oin ts. Panel  [* * (b) **]{} Low estenergy ban d f orthe e ffe c tivepote n ti ali n Eq.  (\[ V eff\]) wi t h$\p s i _s $ taken ast h e  gr oundsta t e of t he s ystem and for $\c h i=-1$ (attract ivec a se) . P a rame te rs are fixed a s i np a nel (a).  P anel [**(c
. _\label{Veff}$$ Here_$\hat U_{ol} \equiv A_\cos(2 x)_$_is the_OL_and $\hat U_s$_is the potential_associated with a given_eigenstate of the_quantum_problem (\[schro\]). For a self-consistent solution, one starts with a trial wavefunction for $\psi_s$_(typically_a gaussian_waveform),_calculates_the effective potential and solves_the corresponding eigenvalue problem (\[schro\])._Then, one_selects a given eigenstate (for example the ground_state_but not necessarily)_as new trial function and iterates the procedure until_convergence is reached.

![ Panel [**(a)**]{} Energy_spectrum for the_effective_potential_(\[Veff\]) with $A=3$ and_$\chi=0$ (Mathieu equation). Full lines represent_exact values of the band edges_of the Mathieu equation while dots are_the eigenvalues obtained with the above_procedure on a lattice of_length $L=40_\pi$, with $N=512$ points. Panel_[**(b)**]{} Lowest energy_band for_the effective potential_in Eq. (\[Veff\]) with $\psi_s$ taken_as the ground_state of the system and for_$\chi=-1$_(attractive case). Parameters_are_fixed_as in_panel (a). Panel_[**(c)**]{}_The same_as_in panel (b) but for $A=-3$._Panel_[**(d)**]{} Transition from the metastable IS mode_to the OS ground_state_corresponding to the lower_level of panel (c). The_optical lattice (scaled by a factor_3) is_reported as_an help to locate the symmetry center of the solutions. Parameters_are fixed as in panel (c).[]{data-label="fig1"}](fig1a.eps_"fig:"){width="3.8cm" height="3.3cm"}![ Panel [**(a)**]{}_Energy spectrum_for_the effective potential_(\[Veff\])_with $A=3$_and $\chi=0$ (Mathieu equation). Full lines represent_exact values_of the band edges of the_Mathieu equation while dots_are_the eigenvalues obtained with the above_procedure on a lattice of length_$L=40 \pi$, with $N=512$ points._Panel_[**(b)**]{}_Lowest energy band for the_effective potential in Eq. (\[Veff\]) with_$\psi_s$ taken as_the ground state of the system and_for_$\chi=-1$ (attractive case). Parameters are fixed_as_in panel (a). Panel [**(c