 \[Ne\_1\], where two phase diagrams are shown in the $p$-$q$ plane with an antiferromagnetic interaction $U_2/U_0\approx0.037\;(\rm ^{23}Na)$ and filling $n=2$ (left) and $3$ (right), respectively. The depth of the optical lattice is chosen to be $V=20\, E_R$ (t/$U_{0} \approx 0.012$), where the system is in a typical Mott insulating regime, with $E_R$ being the recoil energy of atoms. Three distinct phases, namely ferromagnetic phase (FM), nematic insulator (NI) and spin-singlet insulator (SSI), are observed and we employ the values of the condensate order parameter $\phi^1_\alpha\equiv\langle b_\alpha\rangle$, the nematic order parameter $\phi^2_{\alpha\beta} \equiv \langle S_\alpha S_\beta \rangle -\delta_{\alpha\beta}/3 \langle S^2 \rangle$, and the local magnetization $\bf {M}\equiv\langle {\bf S} \rangle$ to characterize them respectively.

For filling $n=2$ and $V=20\, E_R$, our BDMFT calculations predict that different parameters will result in different Mott insulating phases with different types of spin order including FM featured by $\phi^1_\alpha =0 $, ${\bf M} \neq 0$, NI by $\phi^1_\alpha =0 $, $\phi^2_{\alpha\beta}>0$ and ${\bf M}=0$ and SSI by $\phi^1_\alpha =0$, $\phi^2_{\alpha\beta}=0$ and $\langle {\bf S}^2\rangle=0$ with $S$ being the local total spin (due to the symmetry of the spin wave-function on each site, $n+S=({\rm even})$ [@Imambekov]). It is worth mention that the formation of singlet pairs will give rise to an even number of atoms per site which characterize SSI phase. Actually, for the small linear Zeeman energy $p$ (equivalent to magnetization $M_z$) and quadratic Zeeman energy $q$, the system is in a regime dominated by spin-dependent interactions, and favors the SSI phase by forming a singlet pair and lowering the spin-dependent energy. For larger quadratic Zeeman energy $q$, the degeneracy of
\[Ne\_1\ ], where two phase diagrams are shown in the $ p$-$q$ plane with an antiferromagnetic interaction $ U_2 / U_0\approx0.037\;(\rm ^{23}Na)$ and filling $ n=2 $ (leftover) and $ 3 $ (right field), respectively. The astuteness of the optical lattice is chosen to be $ V=20\, E_R$ (t/$U_{0 } \approx 0.012 $), where the organization is in a typical Mott insulating regimen, with $ E_R$ being the recoil department of energy of atoms. Three distinct phases, namely ferromagnetic phase (FM), nematic insulator (NI) and spin - singlet insulator (SSI), are observed and we hire the values of the condensate order argument $ \phi^1_\alpha\equiv\langle b_\alpha\rangle$, the nematic order parameter $ \phi^2_{\alpha\beta } \equiv \langle S_\alpha S_\beta \rangle -\delta_{\alpha\beta}/3 \langle S^2 \rangle$, and the local magnetization $ \bf { M}\equiv\langle { \bf S } \rangle$ to characterize them respectively. 

 For filling $ n=2 $ and $ V=20\, E_R$, our BDMFT calculations predict that different parameter will result in different Mott insulating phase with different types of spin order include FM featured by $ \phi^1_\alpha = 0 $, $ { \bf M } \neq 0 $, NI by $ \phi^1_\alpha = 0 $, $ \phi^2_{\alpha\beta}>0 $ and $ { \bf M}=0 $ and SSI by $ \phi^1_\alpha = 0 $, $ \phi^2_{\alpha\beta}=0 $ and $ \langle { \bf S}^2\rangle=0 $ with $ S$ being the local total spin (due to the symmetry of the spin wave - affair on each site, $ n+S=({\rm even})$   [ @Imambekov ]). It is worth mention that the formation of singlet pairs will give rise to an even number of atoms per site which characterize SSI phase. Actually, for the small linear Zeeman energy $ p$ (equivalent to magnetization $ M_z$) and quadratic Zeeman energy $ q$, the system is in a regime predominate by spin - dependent interaction, and favor the SSI phase by forming a singlet pair and lowering the tailspin - dependent department of energy. For larger quadratic Zeeman energy $ q$, the degeneracy of
 \[Ne\_1\], ahere two phase diagrams are shown in tkw $p$-$q$ pnane wjth an avtiferromagnetic interaction $U_2/Y_0\apprix0.037\;(\rm ^{23}Na)$ and filling $n=2$ (ueft) and $3$ (right), rwspertively. The depti of the opticzp lactmce is chosen tp be $V=20\, E_R$ (d/$U_{0} \approx 0.012$), whese tke system is in a typical Mott insulwting rrglme, with $E_R$ beyng uhe reckpl energy of atoms. Three distindt phasts, namely ferromagmetic phase (FM), nematic insklatlr (NI) and spin-singpet insulatir (SFU), are observdd and we tm'loy the vamues of the condensate order pafametzr $\phi^1_\alpha\gqhig\nangle b_\alpia\rangje$, the nematlb order parameyer $\phi^2_{\alpha\beba} \eqniv \oangle S_\alpha S_\beta \rengle -\delta_{\alpha\beta}/3 \langle S^2 \rcngle$, and the local mqgbetizdtiot $\bf {N}\eqjiv\maigls {\bf S} \raigle$ to chadacterize tyem respectively.

For fyoling $n=2$ and $V=20\, E_R$, ouw FDMFT calculations predict that differett larameters will result un different Mott insolating phwses with different types of spin order including FM fxagurtd by $\pfu^1_\appha =0 $, ${\bf M} \neq 0$, NI by $\phi^1_\alpha =0 $, $\phi^2_{\alpha\betw}>0$ amd ${\bf M}=0$ and SSI ny $\phi^1_\alpha =0$, $\phi^2_{\allhw\brja}=0$ and $\langle {\cf S}^2\raufls=0$ with $S$ being the local jotal wpin (due uo thr symmetry of the spin wave-dunction on vach site, $n+S=({\rm even})$ [@Imalbekov]). It ir wotth memtion that the formatiou of sjnglet paird will gibd rise to an evev nlmbes of atoms per site which sharactermze SXI phasd. Acjually, sor the smwll lltear Zeeman energy $p$ (eqoivalett to magnftization $M_z$) and quadratic Zeemei energy $q$, thg sfstvm is in c regike dominated fy spin-dependeut intercctionr, and favogs the SSM phase by fjrming a singngt pair and lmwering ehe wpin-eependevg energy. For lsrger quabxatic Zeenan energy $q$, the dcgenetady of
\[Ne\_1\], where two phase diagrams are shown $p$-$q$ with an interaction $U_2/U_0\approx0.037\;(\rm ^{23}Na)$ $3$ respectively. The depth the optical lattice chosen to be $V=20\, E_R$ (t/$U_{0} 0.012$), where the system is in a typical Mott insulating regime, with $E_R$ the recoil energy of atoms. Three distinct phases, namely ferromagnetic phase (FM), nematic (NI) spin-singlet (SSI), observed and we employ the values of the condensate order parameter $\phi^1_\alpha\equiv\langle b_\alpha\rangle$, the nematic order $\phi^2_{\alpha\beta} \equiv \langle S_\alpha S_\beta \rangle -\delta_{\alpha\beta}/3 \langle \rangle$, and the local $\bf {M}\equiv\langle {\bf S} \rangle$ characterize respectively. For $n=2$ $V=20\, our BDMFT calculations that different parameters will result in different Mott insulating phases with different types of spin order including featured by $, ${\bf \neq NI $\phi^1_\alpha =0 $, ${\bf M}=0$ and SSI by $\phi^1_\alpha $\langle {\bf S}^2\rangle=0$ with $S$ being the local spin (due the symmetry of the spin wave-function each site, $n+S=({\rm even})$ [@Imambekov]). It is worth that the formation of singlet pairs will give rise to an even number of atoms which characterize SSI phase. for the small Zeeman $p$ to $M_z$) and Zeeman energy $q$, the system is in a regime dominated by interactions, and favors the SSI phase by forming a singlet lowering spin-dependent energy. For quadratic Zeeman energy $q$, degeneracy
 \[Ne\_1\], where two phase diagrams arE shown in thE $p$-$q$ plAne WitH aN antIferRomagnetic inteRActiOn $U_2/U_0\approx0.037\;(\rm ^{23}Na)$ and filliNg $n=2$ (leFt) ANd $3$ (riGHt), RespeCtively. tHe DEPth Of ThE opTiCAl LattiCe iS chosen To be $V=20\, E_R$ (t/$U_{0} \aPprOx 0.012$), Where the systEM iS in a typicaL MoTt insulating RegIme, witH $E_r$ beINg the RecOil enErgy of AToms. ThRee distinCt PHases, nAMely ferROMaGnetIc phase (FM), nematic iNSuLAtor (NI) and spin-sInglet InSUlATOr (SsI), aRe observed AnD we emPLoy the vALuES OF thE Condensate ordEr parameter $\PHi^1_\aLpha\eqUiV\laNGle b_\alPha\raNgLE$, thE nematic ordEr paRameter $\phI^2_{\alpha\BEta} \equiV \Langle S_\Alpha S_\BetA \raNgle -\DElTa_{\AlpHa\BEta}/3 \LAnGle s^2 \RanGle$, and thE lOcAl magNetiZATIOn $\bf {m}\eqUiv\lAngle {\Bf S} \rangle$ to chAraCterIZe tHem reSpectIvelY.

FOr filLing $n=2$ aNd $V=20\, E_R$, OuR BDMFT calculatiOns pRedict thaT diFfEreNt ParamETers wiLl rEsuLt in difFerent MOTt iNsULATiNg phases with differEnT TYpEs of spin Order iNClUdINg FM featUrEd bY $\phi^1_\ALPha =0 $, ${\bf m} \neq 0$, ni bY $\phi^1_\alphA =0 $, $\phi^2_{\alPHa\BeTa}>0$ and ${\bf m}=0$ aNd SSI bY $\pHi^1_\aLphA =0$, $\phi^2_{\aLPha\bEta}=0$ and $\Langle {\bf s}^2\rangLE=0$ with $S$ being the LOcal total spin (DUe TO ThE SymmEtrY of the spin wAve-fUNctiOn on EAcH siTE, $n+S=({\rm Even})$ [@IMaMBeKOv]). It is worth mention tHaT the foRmatiOn of singlet paIrs will givE RISe to an evEn nuMBeR Of atoms per site Which CharacteriZE SSI phasE. ActuAlly, for tHe small liNEAr Zeeman EneRgy $P$ (eqUivALEnT to magnetizatION $M_z$) aNd QuadratIc ZEeman enErgY $q$, tHe sYstEm Is in a regiMe dominaTeD bY sPiN-dePendeNT interacTiOns, AnD faVors tHE SSI phAse by FormInG a SIngLet pair ANd LOWeriNg ThE spiN-dePeNdent EnerGY. FoR larger Quadratic zeeMAn enErGy $Q$, the degEneracy of
 \[Ne\_1\], where two phas e diagrams  aresho wnin  the  $p$ -$q$ plane wit h  anantiferromagnetic inte racti on $U_2 / U_ 0\app rox0.03 7 \; ( \ rm^{ 23 }Na )$ an d fil lin g $n=2$  (left) an d $ 3$  (right), re s pe ctively. T hedepth of the  op ticalla tti c e ischo sen t o be $ V =20\,E_R$ (t/$ U_ { 0} \ap p rox 0.0 1 2 $) , wh ere the system is in a typical Mott  insul at i ng r egi me,  with $E_R $being the rec o il e n erg y  of atoms. Th ree distinc t  ph ases,na mel y  ferro magne ti c  ph ase (FM), n emat ic insula tor (N I ) and s p in-sing let in sul ato r (S S I) ,are  o b ser v ed  an d  we  employth evalue s of t h e  con den sate  orde r parameter $ \ph i^1_ \ alp ha\eq uiv\l angl eb_\al pha\ra ngle$ ,the nematic ord er p arameter$\p hi ^2_ {\ alpha \ beta}\eq uiv  \langl e S_\al p haS_ \ b e ta  \rangle -\delta_{ \a l p ha \beta}/3  \lang l eS^ 2  \rangle $,  an d th e local  mag n et ization$\bf { M }\ eq uiv\lan gl e {\bf  S } \ ran gle$t o ch aracte rize the m res p ectively.

For filling $n=2$ an d $V = 20\,  E_ R$, our BDM FT c a lcul atio n spre d ict t hat d if f er e nt parameters willre sult i n dif ferent Mott i nsulatingp h a ses with  dif f er e nt types of sp in or der includ i ng FM fe ature d by $\p hi^1_\alp h a  =0 $, $ {\b f M } \ neq 0 $,  NI by $\phi^ 1 _ \alp ha  =0 $,$\p hi^2_{\ alp ha\ bet a}> 0$  and ${\b f M}=0$an dSS Iby$\phi ^ 1_\alpha  = 0$,  $ \ph i^2_{ \ alpha\ beta} =0$an d$ \la ngle {\ b fS } ^2\r an gl e=0$  wi th  $S$bein g  th e local  total sp in( dueto  t he symm etry of the s pi n wave-fun ct ion  on ea c h  site, $ n+S=({\rm even})$ [@Ima m bekov]) . I t iswort h mention  th at the  fo r mation  of si nglet  p air s willg i ve  ri se  to an eve n num ber o fatom s per s ite which characte r ize  SSI phase. A ctu ally , fo r t h es mal ll ine a r  Zeeman energy$p$ (equiv al e nt  to magnet i zat io n $M_z$ ) and q uadra t ic Zeem an energy  $q$, the  s yste m isin a regim e domina ted by sp i n-dep e nd ent i nte ractio ns , a nd fa vors t h e S SI ph ase by  f orming  a si ng let pair  and lowering the spin- depend ent e ner gy. For l arg e r q uadraticZeem an energy$q$ , t he de gen e racyof
 \[Ne\_1\], where_two phase_diagrams are shown in_the $p$-$q$_plane_with an_antiferromagnetic_interaction $U_2/U_0\approx0.037\;(\rm ^{23}Na)$_and filling $n=2$_(left) and $3$ (right),_respectively. The depth_of_the optical lattice is chosen to be $V=20\, E_R$ (t/$U_{0} \approx 0.012$), where the_system_is in_a_typical_Mott insulating regime, with $E_R$_being the recoil energy of_atoms. Three_distinct phases, namely ferromagnetic phase (FM), nematic insulator_(NI)_and spin-singlet insulator_(SSI), are observed and we employ the values of_the condensate order parameter $\phi^1_\alpha\equiv\langle b_\alpha\rangle$,_the nematic order_parameter_$\phi^2_{\alpha\beta}_\equiv \langle S_\alpha S_\beta_\rangle -\delta_{\alpha\beta}/3 \langle S^2 \rangle$, and_the local magnetization $\bf {M}\equiv\langle {\bf_S} \rangle$ to characterize them respectively.

For filling_$n=2$ and $V=20\, E_R$, our BDMFT_calculations predict that different parameters_will result_in different Mott insulating phases_with different types_of spin_order including FM_featured by $\phi^1_\alpha =0 $, ${\bf_M} \neq 0$,_NI by $\phi^1_\alpha =0 $, $\phi^2_{\alpha\beta}>0$_and_${\bf M}=0$ and_SSI_by_$\phi^1_\alpha =0$,_$\phi^2_{\alpha\beta}=0$ and $\langle_{\bf_S}^2\rangle=0$ with_$S$_being the local total spin (due_to_the symmetry of the spin wave-function on_each site, $n+S=({\rm even})$ [@Imambekov])._It_is worth mention that_the formation of singlet pairs_will give rise to an even_number of_atoms per_site which characterize SSI phase. Actually, for the small linear Zeeman_energy $p$ (equivalent to magnetization $M_z$)_and quadratic Zeeman energy_$q$, the_system_is in a_regime_dominated by_spin-dependent interactions, and favors the SSI phase_by forming_a singlet pair and lowering the_spin-dependent energy. For larger_quadratic_Zeeman energy $q$, the degeneracy of