\mu_2&=&E,\\ 
\mu_3&=&\frac{1}{2}(D-E-h'), \nonumber
\label{case3}\end{aligned}$$ where $h'= \sqrt{(D+E)^2+4g_{\perp}^2\mu_B^2H_{a}^2}$ with $a=x$ or $y$ [@foot1]. Subsequently, we have $$\begin{aligned}
C_0&=&   B_{(D-E)/2,D},\nonumber\\
A_+&=&B_{(D+E)/2,D-E}/B_{(D-E)/2,D}^2, \label{SP1c3}\\
A_-&=&\exp(D/T)/B_{(D-E)/2,D}^3, \nonumber\end{aligned}$$ with now $B_{x,y}=2\exp(x/T)\cosh(h'/T)+\exp(y/T)$.

Eq. (\[fe\]) for the free energy $f(T,H)$ is our key result. Physical properties such as the susceptibility, magnetization and the specific heat follow in the usual way by differentiation. We also use $f(T,H)$ to calculate the phase diagram for both $T \simeq 0$ K and finite temperatures (see figure \[PD-NENC\]). We find that considering up to 3rd order in $J/T$ is sufficient as higher orders are negligibly small. This is in stark contrast to other series expansions which need many orders to accurately describe physical properties. This is mainly because here the coefficents are not just constants, but functions of the external model parameters, e.g., the magnetic field and the coupling strength.

Spin-$1$ compounds
==================

The compound NENC
-----------------

It is known that antiferromagnetic spin-$1$ chains [@SP1C1; @SP1C2] with weak planar anisotropy can exhibit a non-magnetic gapped phase. The large $D$ gapped phase has been observed in the compounds NENC, NDPK and NBYC [@NENC; @NBYC]. In these compounds the in-plane anisotropy $x^2-y^2$ breaks the $z^2$ symmetry and weakens the planar anisotropy. From experimental analysis, it was inferred that the in-plane anisotropy $E$ in NENC [@NENC] is negligible in comparison with
\mu_2&=&E,\\ 
 \mu_3&=&\frac{1}{2}(D - E - h'), \nonumber 
 \label{case3}\end{aligned}$$ where $ h'= \sqrt{(D+E)^2 + 4g_{\perp}^2\mu_B^2H_{a}^2}$ with $ a = x$ or $ y$ [ @foot1 ]. Subsequently, we have $ $ \begin{aligned } 
 C_0&= &    B_{(D - E)/2,D},\nonumber\\ 
 A_+&=&B_{(D+E)/2,D - E}/B_{(D - E)/2,D}^2, \label{SP1c3}\\ 
 A_-&=&\exp(D / T)/B_{(D - E)/2,D}^3, \nonumber\end{aligned}$$ with now $ B_{x, y}=2\exp(x / T)\cosh(h'/T)+\exp(y / T)$. 

 Eq.   (\[fe\ ]) for the spare department of energy $ f(T, H)$ is our key result. forcible properties such as the susceptibility, magnetization and the specific estrus follow in the usual manner by differentiation. We also practice $ f(T, H)$ to count the phase diagram for both $ T \simeq 0 $ kilobyte and finite temperatures (see design \[PD - NENC\ ]). We rule that considering up to 3rd order in $ J / T$ is sufficient as higher order are negligibly small. This is in stark contrast to early series expansions which necessitate many orders to accurately report physical properties. This is mainly because here the coefficents are not just constants, but routine of the external model parameters, e.g., the magnetic field and the coupling strength. 

 Spin-$1 $ compounds 
 = = = = = = = = = = = = = = = = = = 

 The compound NENC 
 ----------------- 

 It is known that antiferromagnetic spin-$1 $ chains [ @SP1C1; @SP1C2 ] with weak planar anisotropy can exhibit a non - magnetic gapped phase. The large $ D$ gapped phase has been observed in the compounds NENC, NDPK and NBYC [ @NENC; @NBYC ]. In these compound the in - plane anisotropy $ x^2 - y^2 $ breaks the $ z^2 $ isotropy and sabotage the planar anisotropy. From experimental psychoanalysis, it was inferred that the in - plane anisotropy $ E$ in NENC [ @NENC ] is negligible in comparison with
\mu_2&=&E,\\ 
\mu_3&=&\frac{1}{2}(D-E-h'), \nonumber
\label{gase3}\end{aligned}$$ wkwre $h'= \vqrt{(D+E)^2+4f_{\perp}^2\mu_B^2F_{a}^2}$ with $a=x$ or $y$ [@foot1]. Subsequxntlt, we yave $$\begin{aligned}
C_0&=&   B_{(A-E)/2,D},\nonumbvr\\
A_+&=&B_{(D+E)/2,D-E}/B_{(E-E)/2,D}^2, \oqbel{SP1c3}\\
A_-&=&\ex'(S/T)/B_{(D-E)/2,D}^3, \kjnumgcr\end{clmgned}$$ with now $N_{x,y}=2\exp(x/T)\cosv(h'/T)+\exp(y/T)$.

Eq. (\[fe\]) fmr tke free energy $f(T,H)$ is our key result. Physicsl properties sush ax the susceptibility, magnetization and fhe spebific heat follow in the usual way by diffegentlation. We also use $f(T,H)$ to calxulaew the phase aiagram fog both $T \sijeq 0$ K and finite temperatures (ree fngure \[PD-NENX\]). Qe vhnd that coisiderpng up to 3rd order in $J/T$ is xufficient as mighec oreers are negligibly skall. This is in stwrk contrdsc to other series expqnwions whiwh ndwd oang prsers tl arcurately dsscribe phywical properties. Thos px mainly beczuse hqrq the coefficents are not just constantv, bht functions of the extwrnal model parameterd, e.g., the iagnetic field and the coupling strength.

Spin-$1$ compmunds
==================

Vhd cimiouna NFNC
-----------------

It is known that antiferromagnetic spin-$1$ chwjnx [@XP1C1; @SP1C2] with wcak planar anisotrppj vwn exhibit a von-magustjc gapped phase. Thf large $D$ gapped phast has been observed in the compoynds NENC, NDIK abd NBYC [@NENC; @NBYC]. Nn these com'ounds the on-plane anisotropy $x^2-y^2$ bxeaks fhe $z^2$ symmehry and wsxkens the planar anpsotsopy. From experimental anajysis, it xas iuferred ghat the ig-plane anidotroif $E$ in NENC [@NENC] id negpicible in clmparison with
\mu_2&=&E,\\ \mu_3&=&\frac{1}{2}(D-E-h'), \nonumber \label{case3}\end{aligned}$$ where $h'= \sqrt{(D+E)^2+4g_{\perp}^2\mu_B^2H_{a}^2}$ or [@foot1]. Subsequently, have $$\begin{aligned} C_0&=& with $B_{x,y}=2\exp(x/T)\cosh(h'/T)+\exp(y/T)$. Eq. (\[fe\]) the free energy is our key result. Physical properties as the susceptibility, magnetization and the specific heat follow in the usual way differentiation. We also use $f(T,H)$ to calculate the phase diagram for both $T 0$ and temperatures figure \[PD-NENC\]). We find that considering up to 3rd order in $J/T$ is sufficient as higher are negligibly small. This is in stark contrast other series expansions which many orders to accurately describe properties. is mainly here coefficents not just constants, functions of the external model parameters, e.g., the magnetic field and the coupling strength. Spin-$1$ compounds ================== compound NENC is known antiferromagnetic chains @SP1C2] with weak can exhibit a non-magnetic gapped phase. gapped phase has been observed in the compounds NDPK and [@NENC; @NBYC]. In these compounds the anisotropy $x^2-y^2$ breaks the $z^2$ symmetry and weakens planar anisotropy. From experimental analysis, it was inferred that the in-plane anisotropy $E$ in NENC negligible in comparison with
\mu_2&=&E,\\ 
\mu_3&=&\frac{1}{2}(D-E-h'), \nonumber
\label{Case3}\end{aliGned}$$ wHerE $h'= \sQrT{(D+E)^2+4g_{\Perp}^2\Mu_B^2H_{a}^2}$ with $a=x$ or $y$ [@FOot1]. SUbsequently, we have $$\begin{AlignEd}
c_0&=&   b_{(D-E)/2,D},\NOnUmber\\
a_+&=&B_{(D+E)/2,D-E}/B_{(d-e)/2,D}^2, \LABel{sP1C3}\\
A_-&=&\Exp(d/T)/b_{(d-E)/2,d}^3, \nonuMbeR\end{aliGned}$$ with noW $B_{x,Y}=2\eXp(x/T)\cosh(h'/T)+\exP(Y/T)$.

eq. (\[fe\]) for the FreE energy $f(T,H)$ is Our Key resUlT. PhYSical ProPertiEs such AS the suSceptibilItY, MagnetIZation aND ThE speCific heat follow in THe USual way by diffeRentiaTiON. WE ALso Use $F(T,H)$ to calcuLaTe the PHase diaGRaM FOR boTH $T \simeq 0$ K and fiNite temperaTUreS (see fiGuRe \[Pd-nENC\]). We Find tHaT ConSidering up tO 3rd oRder in $J/T$ iS suffiCIent as hIGher ordErs are NegLigIbly SMaLl. thiS iS In sTArK coNTraSt to otheR sErIes exPansIONS WhicH neEd maNy ordErs to accurateLy dEscrIBe pHysicAl proPertIeS. This Is mainLy becAuSe here the coeffiCentS are not juSt cOnStaNtS, but fUNctionS of The ExternaL model pARamEtERS, E.g., The magnetic field anD tHE CoUpling stRength.

sPiN-$1$ cOMpounds
==================

THe ComPounD neNC
-----------------

It Is knOWn That antiFerromAGnEtIc spin-$1$ cHaIns [@SP1C1; @sP1c2] wiTh wEak plANar aNisotrOpy can exHibit A Non-magnetic gapPEd phase. The larGE $D$ GAPpED phaSe hAs been obserVed iN The cOmpoUNdS NEnc, NDPK And NByC [@neNc; @nBYC]. In these compoundS tHe in-plAne anIsotropy $x^2-y^2$ breAks the $z^2$ symMETRy and weaKens THe PLanar anisotropY. From ExperimentAL analysiS, it waS inferreD that the iN-PLane anisOtrOpy $e$ in nENc [@neNc] is negligible IN CompArIson witH
\mu_2&=&E,\\ 
\mu_3&=&\fra c{1}{2}(D- E-h') , \ non um ber\lab el{case3}\end{ a lign ed}$$ where $h'= \sqrt {(D+E )^ 2 +4g_ { \p erp}^ 2\mu_B^ 2 H_ { a }^2 }$  w ith  $ a =x $ or$y$  [@foot 1]. Subseq uen tl y, we have $ $ \b egin{align ed} 
C_0&=&   B_ {(D -E)/2, D} ,\n o numbe r\\ 
A_+& =&B_{( D +E)/2, D-E}/B_{( D- E )/2,D} ^ 2, \lab e l {S P1c3 }\\
A_-&=&\exp(D/ T )/ B _{(D-E)/2,D}^3 , \non um b er \ e nd{ ali gned}$$ wi th  now$ B_{x,y} = 2\ e x p (x/ T )\cosh(h'/T)+ \exp(y/T)$. 
Eq . (\[f e\ ])f or the  free  e n erg y $f(T,H)$is o ur key re sult.P hysical propert ies su chasthes us ce pti bi l ity ,  m agn e tiz ation an dth e spe cifi c h e at f oll ow i n the  usual way by  di ffer e nti ation . Wealso  u se $f (T,H)$  to c al culate the phas e di agram for  bo th  $T  \ simeq 0$ K a ndfin ite tem peratur e s ( se e f ig ure \[PD-NENC\]).We f in d that c onside r in gu p to 3rd  o rde r in $ J/T$is s u ff icient a s high e ror ders ar eneglig ib lysma ll. T h is i s in s tark con trast to other serie s  expansions w h ic h ne e d ma nyorders to a ccur a tely  des c ri bep hysic al pr op e rt i es. This is mainlybe causeherethe coefficen ts are not j u st const ants ,  b u t functions of  theexternal m o del para meter s, e.g.,  the magn e t ic field  an d t hecou p l in g strength.
S p in-$ 1$  compou nds 
====== === === === === 
The compo und NENC 
- -- -- -- --- ----- - -

It is  k now ntha t ant i ferrom agnet ic s pi n- $ 1$chains[ @S P 1 C1;@S P1 C2]wit hweakplan a r a nisotro py can ex hib i t ano n- magneti c gapped phas e.  The large  $ D$gapped p hase has  been observed in the c o mpounds  NE NC, N DPKand NBYC[@N ENC; @ NBY C ]. Inthesecompo un dst h e in- p l an e a ni sotropy $x ^ 2 -y^ 2$ br ea ks t he $z^2 $ symmetry and wea k ens  the planar a nis otro p y .Fro m  e x per im e nta l analysis, it wa s inferred  t h at  the in-pl a nean isotrop y $E$ i n NEN C  [@NENC ] is negl igible in  c ompa r i son  with
\mu_2&=&E,\\ 
\mu_3&=&\frac{1}{2}(D-E-h'),_\nonumber
\label{case3}\end{aligned}$$ where_$h'= \sqrt{(D+E)^2+4g_{\perp}^2\mu_B^2H_{a}^2}$ with $a=x$_or $y$_[@foot1]._Subsequently, we_have_$$\begin{aligned}
C_0&=&  _B_{(D-E)/2,D},\nonumber\\
A_+&=&B_{(D+E)/2,D-E}/B_{(D-E)/2,D}^2, \label{SP1c3}\\
A_-&=&\exp(D/T)/B_{(D-E)/2,D}^3, \nonumber\end{aligned}$$_with now $B_{x,y}=2\exp(x/T)\cosh(h'/T)+\exp(y/T)$.

Eq. (\[fe\]) for_the free energy_$f(T,H)$_is our key result. Physical properties such as the susceptibility, magnetization and the specific_heat_follow in_the_usual_way by differentiation. We also_use $f(T,H)$ to calculate the_phase diagram_for both $T \simeq 0$ K and finite_temperatures_(see figure \[PD-NENC\])._We find that considering up to 3rd order in_$J/T$ is sufficient as higher orders_are negligibly small._This_is_in stark contrast to_other series expansions which need many_orders to accurately describe physical properties._This is mainly because here the coefficents_are not just constants, but functions_of the external model parameters,_e.g., the_magnetic field and the coupling_strength.

Spin-$1$ compounds
==================

The compound_NENC
-----------------

It is_known that antiferromagnetic_spin-$1$ chains [@SP1C1; @SP1C2] with weak_planar anisotropy can_exhibit a non-magnetic gapped phase. The_large_$D$ gapped phase_has_been_observed in_the compounds NENC,_NDPK_and NBYC_[@NENC;_@NBYC]. In these compounds the in-plane_anisotropy_$x^2-y^2$ breaks the $z^2$ symmetry and weakens_the planar anisotropy. From_experimental_analysis, it was inferred_that the in-plane anisotropy $E$_in NENC [@NENC] is negligible in_comparison with