, for small but nonzero magnetization $M_z\neq 0$, we observe a FM phase with $\phi^{1}_{\alpha}= 0$ and ${\bf M}\neq 0$ in between SF and SSI, as shown in Fig. \[order\](e). Interestingly, we observe a first-order SSI-FM and a second-order FM-SF transition, with lowering the optical depth $V$.

### Finite temperature {#section:finite.unnumbered}

  -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
  ![(Color online) Finite-temperature phase diagram for spin-1 ultracold bosons in a 3D cubic lattice with an antiferromagnetic interaction $U_2/U_0=0.037$ for a fixed temperature $T/U_0=0.01$, and a magnetic field $p=20\, Hz$ and $q=120\, Hz$. Note here that the normal state is shown by the shaded region in between FM.[]{data-label="finite_T"}](T1.eps "fig:"){width="0.6\linewidth"}
  -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

  ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
  ![(Color online) Abnormal multi-step condensation of spin-1 ultracold bosons in a 3D cubic lattice ($V=21\,{\rm Hz}$) with an antiferromagnetic interaction $U_2/U_0=0.037$ and filling $n=2.5$ for quadratic Zeeman field $q=20\, Hz$ (left), and $-20\, Hz$ (right), respectively. Contrary to weakly interacting gases, we observe a non-smooth change of the critical condensing temperature of the $m_F=-1$ component.[]{data-label="multi_step"}](FF5.eps "fig:"){width="0.9\linewidth"}
  ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

Stability of the above spin-dependent phases against thermal fluctuations is the first thing to consider in order to obtain a direct observation of them in practical experiments. We investigate this issue for a typical case with $U_2/U_0 = 0.037$ ($^{23}\rm Na$) at a fixed finite temperature $T/U_0=0.01$, and find that four different phases are involved as shown in Fig. \[finite\_T\], including SF ($\phi^1_{0}$ or $\phi^1_{\pm 1}$), FM, and normal state (NS) characterized both by $\phi^1_\alpha = 0$ and large
, for small but nonzero magnetization $ M_z\neq 0 $, we observe a FM phase with $ \phi^{1}_{\alpha}= 0 $ and $ { \bf M}\neq 0 $ in between SF and SSI, as shown in Fig.   \[order\](e). Interestingly, we observe a inaugural - decree SSI - FM and a second - order FM - SF transition, with lowering the ocular depth $ V$. 

 # # # Finite temperature { # section: finite.unnumbered } 

   ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- 
  ! [ (Color online) Finite - temperature phase diagram for spin-1 ultracold boson in a 3D cubic lattice with an antiferromagnetic interaction $ U_2 / U_0=0.037 $ for a fixed temperature $ T / U_0=0.01 $, and a magnetic battlefield $ p=20\, Hz$ and $ q=120\, Hz$. Note here that the normal state is show by the shaded region in between FM.[]{data - label="finite_T"}](T1.eps " fig:"){width="0.6\linewidth " } 
   ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- 

   ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------ 
  ! [ (Color online) Abnormal multi - step condensing of spin-1 ultracold boson in a 3D cubic lattice ($ V=21\,{\rm Hz}$) with an antiferromagnetic interaction $ U_2 / U_0=0.037 $ and filling $ n=2.5 $ for quadratic Zeeman sphere $ q=20\, Hz$ (left), and $ -20\, Hz$ (right), respectively. Contrary to weakly interacting gas, we observe a non - smooth change of the critical condense temperature of the $ m_F=-1 $ component.[]{data - label="multi_step"}](FF5.eps " fig:"){width="0.9\linewidth " } 
   ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------ 

 Stability of the above spin - dependent phases against thermal fluctuations is the inaugural thing to consider in order to obtain a direct notice of them in practical experiments. We investigate this issue for a typical case with $ U_2 / U_0 = 0.037 $ ($ ^{23}\rm Na$) at a fixed finite temperature $ T / U_0=0.01 $, and rule that four different phases are involved as shown in Fig.   \[finite\_T\ ], including SF ($ \phi^1_{0}$ or $ \phi^1_{\pm 1}$), FM, and normal state (NS) characterized both by $ \phi^1_\alpha = 0 $ and large
, fog small but nonzero magnttization $M_z\neq 0$, cw obsecve a FJ phase dith $\phi^{1}_{\alpha}= 0$ and ${\bf M}\neq 0$ mn bwtweeb SF and SSI, as shown kn Fig. \[ordvr\](e). Interwstiigly, we observe e first-ovber SAL-FM aud a second-order FM-SF tranvition, with lofefiug the optical depth $V$.

### Finite temperwture {#srchion:finite.unnuibertd}

  -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
  ![(Cklor online) Finite-temperature phass diagrem for spin-1 ulttacold bosons in a 3D cubic lathice with an antiffrromagnetix inewraction $U_2/U_0=0.037$ wor a fixed temperaturg $T/U_0=0.01$, and a magnetic field $p=20\, Hz$ avd $q=120\, Kz$. Note herg fhwj the normal statv is shown by the shaged regoon in between FM.[]{vata-oabel="finite_T"}](T1.eps "fig:"){wmdth="0.6\linewidth"}
  -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

  ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
  ![(Cjlor onlite) Abnormal multi-stwp condgnsathon ud soin-1 uktdacold boaons in a 3S cubic latrice ($V=21\,{\rm Hz}$) with an agnoferromagnetjc intqrwction $U_2/U_0=0.037$ and filling $n=2.5$ for quadratic Zteman field $q=20\, Hz$ (left), and $-20\, Hz$ (right), respectively. Fontrary eo weakly interacting gases, we observe a non-smootv chaigd oy the zeihical condensing temperature of the $m_F=-1$ componqht.[]{cana-label="multi_step"}](FN5.eps "fig:"){width="0.9\linewodhh"}
  ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

Stability of the aykvs spin-dependent phwses agwinst thermal slucyuations is the first thing to consider un order to obtain a direct oyservajion og them in practical expzrimenfs. We inveshigate thjr issue for a tyoicsl case wiuf $U_2/U_0 = 0.037$ ($^{23}\rm Na$) at a fixed fiiite cemperatjre $J/U_0=0.01$, and sind that vour different phases arf invlleed as shoan in Fig. \[finite\_T\], including SF ($\pim^1_{0}$ or $\phi^1_{\pm 1}$), FK, dnd normal ftate (NS) characteryzed both by $\pki^1_\alpha = 0$ and large
, for small but nonzero magnetization $M_z\neq observe FM phase $\phi^{1}_{\alpha}= 0$ and SF SSI, as shown Fig. \[order\](e). Interestingly, observe a first-order SSI-FM and a FM-SF transition, with lowering the optical depth $V$. ### Finite temperature {#section:finite.unnumbered} ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- online) Finite-temperature phase diagram for spin-1 ultracold bosons in a 3D cubic lattice an interaction for fixed temperature $T/U_0=0.01$, and a magnetic field $p=20\, Hz$ and $q=120\, Hz$. Note here that the state is shown by the shaded region in FM.[]{data-label="finite_T"}](T1.eps "fig:"){width="0.6\linewidth"} ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------ online) Abnormal multi-step condensation of ultracold in a cubic ($V=21\,{\rm with an antiferromagnetic $U_2/U_0=0.037$ and filling $n=2.5$ for quadratic Zeeman field $q=20\, Hz$ (left), and $-20\, Hz$ (right), respectively. Contrary weakly interacting observe a change the condensing temperature of component.[]{data-label="multi_step"}](FF5.eps "fig:"){width="0.9\linewidth"} ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------ Stability of the against thermal fluctuations is the first thing to in order obtain a direct observation of them practical experiments. We investigate this issue for a case with $U_2/U_0 = 0.037$ ($^{23}\rm Na$) at a fixed finite temperature $T/U_0=0.01$, and find different phases are involved shown in Fig. including ($\phi^1_{0}$ $\phi^1_{\pm FM, and state (NS) characterized both by $\phi^1_\alpha = 0$ and large
, for small but nonzero magnetiZation $M_z\neQ 0$, we obSerVe a fM PhasE witH $\phi^{1}_{\alpha}= 0$ and ${\bf m}\Neq 0$ iN between SF and SSI, as showN in FiG. \[oRDer\](e). iNtErestIngly, we OBsERVe a FiRsT-orDeR sSi-FM anD a sEcond-orDer FM-SF traNsiTiOn, with loweriNG tHe optical dEptH $V$.

### Finite tempEraTure {#seCtIon:FInite.UnnUmberEd}

  -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
  ![(ColoR Online) finite-temPeRAture pHAse diagRAM fOr spIn-1 ultracold bosons IN a 3d Cubic lattice wiTh an anTiFErROMagNetIc interactIoN $U_2/U_0=0.037$ foR A fixed tEMpERATurE $t/U_0=0.01$, and a magnetiC field $p=20\, Hz$ anD $Q=120\, Hz$. note heRe ThaT The norMal stAtE Is sHown by the shAded Region in bEtween fm.[]{data-laBEl="finitE_T"}](T1.eps "Fig:"){WidTh="0.6\liNEwIdTh"}
  -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

  ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
  ![(COlOR onLInE) AbNOrmAl multi-sTeP cOndenSatiON OF Spin-1 UltRacoLd bosOns in a 3D cubic lAttIce ($V=21\,{\RM Hz}$) With aN antiFerrOmAgnetIc inteRactiOn $u_2/U_0=0.037$ and filling $n=2.5$ foR quaDratic ZeeMan FiEld $Q=20\, HZ$ (left), ANd $-20\, Hz$ (riGht), ResPectiveLy. ContrARy tO wEAKLy Interacting gases, we ObSERvE a non-smoOth chaNGe Of THe criticAl ConDensING tempEratURe Of the $m_F=-1$ cOmponeNT.[]{dAtA-label="mUlTi_step"}](fF5.Eps "Fig:"){Width="0.9\LInewIdth"}
  ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

StAbility oF the aBOve spin-dependeNT phases againsT ThERMaL FlucTuaTions is the fIrst THing To coNSiDer IN ordeR to obTaIN a DIrect observation of tHeM in praCticaL experiments. WE investigaTE THis issue For a TYpICal case with $U_2/U_0 = 0.037$ ($^{23}\rM Na$) at A fixed finiTE temperaTure $T/u_0=0.01$, and find That four dIFFerent phAseS arE inVolVED aS shown in Fig. \[fiNITe\_T\], iNcLuding Sf ($\phI^1_{0}$ or $\phi^1_{\pM 1}$), FM, And NorMal StAte (NS) charActerizeD bOtH bY $\pHi^1_\aLpha = 0$ aND large
, for small but nonzero ma gnetizatio n $M_ z\n eq0$ , we  obs erve a FM phas e  wit h $\phi^{1}_{\alpha}=0$ an d$ {\bf M} \neq0$ in b e tw e e n S Fan d S SI ,  a s sho wnin Fig.  \[order\] (e) .Interestingl y ,we observe  afirst-orderSSI -FM an da s e cond- ord er FM -SF tr a nsitio n, with l ow e ring t h e optic a l  d epth  $V$.

### Finite te m perature {#sec tion:f in i te . u nnu mbe red}

  -- -- ----- - ------- - -- - - - --- - ------------- ----------- - --- ------ -- --- - ------ ----- -- - --- ----------- ---- --------- ------ - ------- - ------- ------ --- --- ---- - -- -- --- -- - --- - -- --- - --- -------- -- -- ----- ---- - - - - ---- --- ---- ----- ------------- --- ---- - --- ----- ----- ---- -- ----- ------ ----- -- --------------- ---- --------- --- -- --- -- ----- - ------ --- --- ------- ------- - --- -- - - - -- ------------------ -- - - -- -------- ---
 ! [( Co l or onlin e)  Fi nite - t emper atur e  p hase dia gram f o rsp in-1 ul tr acoldbo son s i n a 3 D  cub ic lat tice wit h ana ntiferromagnet i c interaction $U _ 2 /U _ 0=0. 037 $ for a fix ed t e mper atur e  $ T/U _ 0=0.0 1$, a nd am agnetic field $p=20 \,  Hz$ a nd $q =120\, Hz$. N ote here t h a t  the nor mals ta t e is shown bythe s haded regi o n in bet weenFM.[]{da ta-label= " f inite_T" }]( T1. eps  "f i g :" ){width="0.6\ l i newi dt h"}
  - --- ------- --- --- --- --- -- --------- -------- -- -- -- -- --- ----- - -------- -- --- -- --- ----- - ------ ----- ---- -- -- - --- ------- - -- - - ---- -- -- ---- --- -- ----- ---- - --- ------- --------- --- - ---- -- -- ------- ------------- -- ---------- -- --- ------ - - -------- ----------------------- - ------- --- ----- ---- --------- --- ------ --- - ------ ------ ----- -- --- - - ----- - - -- --- -- ---------- - - --- ----- -- ---- ------- ------------------ - -
  ----------- --- ---- - - -- --- - -- - --- -- - --- - - --------------- ---------- -- - -- ---------- - --- -- ------- ------- ----- - ------- --------- --------- -- ---- - - --- ---------- -------- --------- - ----- - -- ----- --- ------ -- --- ----- ------ - --- ----- ------ -- ------ ----- -- -------- ----------------------- ------ ----- --- --------- --- - --- --------- ---- ---------- --- --- ----- --- - ----- ---- - -- --- - ----- ---- - --------- - -- --- - - -- ----------- - - - --- ----- --- - ------ ---- ----------------- - -------------- ---- - - --- --- - ---- -- -----------
 ![( Co l o r online )Abnormal mu lti-step  c o ndens ationof spi n-1 ult r a co l d boso ns i n a  3D cubic  la tt i ce ($V= 21 \, { \rm Hz }$)wi th anantife r roma g n etic interaction  $U_2 / U _0=0. 0 37$  andfi lling $ n =2.5 $ for quad ratic Zeema n fiel d $q =20\,  Hz$ (l ef t), an d $ -2 0\, Hz$ (r i ght), res pecti vely. C on trar y t o weak ly i n t eract ingga ses , we obse r v ea  n on - smo othchang eof t he critic a l conden sin g  temper at ure o f the$ m_ F = -1$ compon ent .[]{d a t a-label="m u lti_ s te p "}](F F5.eps  "fig: "){widt h ="0 .9 \linewi dth " } 
  ------ --------- - --- -- ---- -------- -- ---- -- --- -- --------------- -------- -------- - - ------ --- - ---- - ------ -- ------------ ----- -- - ------- - - --- ------- - --- --- --- ---- --- ---- - ------------------------ ---- ----- - --- ----- ----- -- -------- - ---- -------- - ------------- --------------- - --- - ---- --- - -- -- --------- --- - ----- ---- --- ---- --- ---- - - ---- - ---- - -- ---- ------ ---- -- -- --------- -- ------- - ---- - --- ---- ----- -- -- - -- ------- -------- ----- -- --------- - ------- ---- -- ------------- ------------ -- ---- ------ - --- --- - 
Stability o f the abov e s pin-depend en t p h ase s aga inst thermalflu ct u ati onsis  the f irst thi ng tocons i d er in  orde r to  obta in a d ire ctobse rv atio n  oft he m in practic a lexperim ents.  Weinv es t i gate this iss u e fo r  a typica l cas e w it h  $ U_2/U_0 = 0.037$  ($^{23 }\rm Na$) at a fixedfinite t empe ra ture  $T /U_ 0= 0 .0 1$ ,andfind that  f our  different  phases a re  invo lved as shown i n F ig. \[fin ite\ _T\], i ncl uding SF ($ \phi^1_{ 0}$ or $\phi^1_{ \pm 1}$), FM , andnor mal sta te (NS)c h aracte rized both  b y $\ phi^1_\ a lpha  = 0 $ and l arge 
, for_small but_nonzero magnetization $M_z\neq 0$,_we observe_a_FM phase_with_$\phi^{1}_{\alpha}= 0$ and_${\bf M}\neq 0$_in between SF and_SSI, as shown_in_Fig. \[order\](e). Interestingly, we observe a first-order SSI-FM and a second-order FM-SF transition, with lowering_the_optical depth_$V$.

###_Finite_temperature {#section:finite.unnumbered}

  -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 _![(Color online) Finite-temperature phase diagram_for spin-1_ultracold bosons in a 3D cubic lattice with_an_antiferromagnetic interaction $U_2/U_0=0.037$_for a fixed temperature $T/U_0=0.01$, and a magnetic field_$p=20\, Hz$ and $q=120\, Hz$. Note_here that the_normal_state_is shown by the_shaded region in between FM.[]{data-label="finite_T"}](T1.eps "fig:"){width="0.6\linewidth"}
_ -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

  ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
  ![(Color_online) Abnormal multi-step condensation of spin-1 ultracold_bosons in a 3D cubic lattice_($V=21\,{\rm Hz}$) with an antiferromagnetic_interaction $U_2/U_0=0.037$_and filling $n=2.5$ for quadratic_Zeeman field $q=20\,_Hz$ (left),_and $-20\, Hz$_(right), respectively. Contrary to weakly interacting_gases, we observe_a non-smooth change of the critical_condensing_temperature of the_$m_F=-1$_component.[]{data-label="multi_step"}](FF5.eps_"fig:"){width="0.9\linewidth"}
 _------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

Stability of the_above_spin-dependent phases_against_thermal fluctuations is the first thing_to_consider in order to obtain a direct_observation of them in_practical_experiments. We investigate this_issue for a typical case_with $U_2/U_0 = 0.037$ ($^{23}\rm Na$)_at a_fixed finite_temperature $T/U_0=0.01$, and find that four different phases are involved as_shown in Fig. \[finite\_T\], including SF ($\phi^1_{0}$_or $\phi^1_{\pm 1}$), FM,_and normal_state_(NS) characterized both_by_$\phi^1_\alpha =_0$ and large