B,C}\sum_{j_3=A,B,C}$.

In order to speed up the calculations the parallelization of the code is used. Additionally, due to the presence of the Boltzman factor in Eq. (\[eq:M\_cf\_Center\]), only the $k_{max} < W$ lowest energy levels and eigenfunctions need to be calculated in order to obtain a very good approximation of $M(B,T)$ curves. Here $W$ is equal to $25$, $25^2$, $25^3$, $25^4$ for a single ion, pair, triplet or quartet, respectively. In each case, the condition of $E_{k_{max}} - E_{k=0} \gg 30 k_{\mathrm{B}}T $ was fullfilled, ensuring that relative error in claculation of $M(B,T)$ is practically zero. $E_{k=0}$ is the ground state energy and $E_{k_{max}}$ is the maximal energy of the calculated excited states.

Magnetic simulations
====================

![\[Fig:Magnetization\] Magnetization per one ion as a function of the magnetic field $B$ of different Mn$^{3+}$ magnetic clusters in GaN at $T$=2 K obtained using crystal field model with ferromagnetic superexchange coupling $J$ = 2 meV. The magnetic easy axis $M_{\perp}$ (solid lines) is perpendicular to the $\textbf{c}$ axis of GaN, whereas the hard one $M_{||}$ (dash lines) is parallel to the $\textbf{c}$ axis. (a) The medium magnetic field, and (b) the high magnetic field region.](MagnetizUp4Ions.pdf){width="8.6"}

![\[Fig:OpenVsClosed\] Comparision of magnetization per one ion of magnetic clusters composed of three ions (open vs closed triplet) and four ions (string vs tetrahedron quartet)](OpenVsClosed.pdf){width="8.6"}

In Fig. \[Fig:Magnetization\] we present the results of our simulations of magnetization per one ion as a function of magnetic field $B$ of different magnetic clusters at $T$=2 K. We see that $M(B)$ varies sharply with magnetic field and the saturation is observed for high $B$. The magnetic easy axis ($M_{\perp}$) is perpendicular to
B, C}\sum_{j_3 = A, B, C}$. 

 In order to speed up the calculations the parallelization of the code is used. Additionally, due to the bearing of the Boltzman agent in Eq.   (\[eq: M\_cf\_Center\ ]), only the $ k_{max } < W$ lowest department of energy levels and eigenfunctions want to be calculated in order to receive a very good approximation of $ M(B, T)$ curve. Here $ W$ is equal to $ 25 $, $ 25 ^ 2 $, $ 25 ^ 3 $, $ 25 ^ 4 $ for a single ion, couple, triplet or quartet, respectively. In each case, the condition of $ E_{k_{max } } - E_{k=0 } \gg 30 k_{\mathrm{B}}T $ was fullfilled, guarantee that relative error in claculation of $ M(B, T)$ is practically zero. $ E_{k=0}$ is the ground state department of energy and $ E_{k_{max}}$ is the maximal energy of the calculated excited state. 

 Magnetic simulations 
 = = = = = = = = = = = = = = = = = = = = 

! [ \[Fig: Magnetization\ ] Magnetization per one ion as a affair of the magnetic battlefield $ B$ of different Mn$^{3+}$ magnetic clusters in GaN at $ T$=2   K obtained using quartz glass field model with ferromagnetic superexchange coupling $ J$ = 2   meV. The magnetic easy axis $ M_{\perp}$ (upstanding lines) is perpendicular to the $ \textbf{c}$ axis of GaN, whereas the hard one $ M_{||}$ (dash lines) is parallel to the $ \textbf{c}$ axis. (a) The medium magnetic field, and (b) the eminent magnetic field region.](MagnetizUp4Ions.pdf){width="8.6 " } 

! [ \[Fig: OpenVsClosed\ ] Comparision of magnetization per one ion of charismatic bunch composed of three ions (open vs closed three) and four ions (string vs tetrahedron quartet)](OpenVsClosed.pdf){width="8.6 " } 

 In Fig.   \[Fig: Magnetization\ ] we present the results of our model of magnetization per one ion as a function of magnetic field $ B$ of different magnetic clusters at $ T$=2   K. We visualize that $ M(B)$ varies sharply with magnetic field and the saturation is observed for eminent $ B$. The magnetic easy axis ($ M_{\perp}$) is perpendicular to
B,C}\skm_{j_3=A,B,C}$.

In order to speed mp the calculations the paralmelizatiun of the code is used. Additmonaoly, dye to the presence of ghe Boltzlan factir ii Eq. (\[eq:M\_cf\_Center\]), only the $k_{max} < A$ locewt energy leveks and eiganfunctions nead tl be calculated in order to obtain w very bold approximatijn og $M(B,T)$ buvves. Here $W$ is equal to $25$, $25^2$, $25^3$, $25^4$ for z singlt ion, pair, triplet or quartet, respectively. Ij eafh case, the conditlon of $E_{k_{maz}} - E_{h=0} \gg 30 k_{\mathrm{C}}T $ was fullfilled, ensoring that relative error in clazulatnon of $M(B,T)$ us prwwtically zeco. $E_{k=0}$ ps the ground state etergy amd $E_{k_{max}}$ is thc maxmmal energy of the calculeted excited states.

Mwgnetic shmblations
====================

![\[Fig:Magnetizatuob\] Magtetisatiub pdr knx ikn as w fnnction of fhe magnetix field $B$ of differtnt Nn$^{3+}$ magnetic cmusterf yn GaN at $T$=2 K obtained using crystal fiend jodel with ferromagnetix superexchange coupllng $J$ = 2 meD. The magnetic easy axis $M_{\perp}$ (solid lines) is per[endiruuar tj rhf $\textbf{c}$ axis of GaN, whereas the hard one $M_{||}$ (qzsn kines) is paraljel to the $\yedtns{c}$ axis. (a) The medium mzgnetic field, and (h) the hygh mqgnetic fyeld region.](MagnetizUp4Ions.pdf){widrh="8.6"}

![\[Fig:OpenVsCjised\] Comparision oy magnetizatnon pet one oon of magnetic clusterr cojposed of tjree ions (upen vs closed tfipked) and four ions (string vs eetrahedrin qbartet)](OpdnVsVlosed.[df){width="8.6"}

In Fig. \[Flc:Magnetization\] we oreseut tha results lf our simulations of magnetizavmon per one ipn as a functnon of magnetic fiejd $B$ of differgnt magnecic cljsters at $N$=2 K. We see that $M(B)$ vawies sharply flth magnetic field agd tye sqturatiuv is observed gor high $Y$. The magnetic easy axix ($M_{\kedp}$) is perpendicbuqr to
B,C}\sum_{j_3=A,B,C}$. In order to speed up the parallelization the code used. Additionally, due Boltzman in Eq. (\[eq:M\_cf\_Center\]), the $k_{max} < lowest energy levels and eigenfunctions need be calculated in order to obtain a very good approximation of $M(B,T)$ curves. $W$ is equal to $25$, $25^2$, $25^3$, $25^4$ for a single ion, pair, or respectively. each the condition of $E_{k_{max}} - E_{k=0} \gg 30 k_{\mathrm{B}}T $ was fullfilled, ensuring that relative error claculation of $M(B,T)$ is practically zero. $E_{k=0}$ is ground state energy and is the maximal energy of calculated states. Magnetic ==================== Magnetization one ion as function of the magnetic field $B$ of different Mn$^{3+}$ magnetic clusters in GaN at $T$=2 K obtained crystal field ferromagnetic superexchange $J$ 2 The magnetic easy (solid lines) is perpendicular to the GaN, whereas the hard one $M_{||}$ (dash lines) parallel to $\textbf{c}$ axis. (a) The medium magnetic and (b) the high magnetic field region.](MagnetizUp4Ions.pdf){width="8.6"} ![\[Fig:OpenVsClosed\] of magnetization per one ion of magnetic clusters composed of three ions (open vs closed four ions (string vs quartet)](OpenVsClosed.pdf){width="8.6"} In Fig. we the of simulations of per one ion as a function of magnetic field $B$ of magnetic clusters at $T$=2 K. We see that $M(B)$ varies magnetic and the saturation observed for high $B$. magnetic axis ($M_{\perp}$) is perpendicular
B,C}\sum_{j_3=A,B,C}$.

In order to speed up The calculaTions The ParAlLeliZatiOn of the code is uSEd. AdDitionally, due to the presEnce oF tHE BolTZmAn facTor in Eq. (\[EQ:M\_CF\_cenTeR\]), oNly ThE $K_{mAx} < W$ loWesT energy Levels and eIgeNfUnctions need TO bE calculateD in Order to obtaiN a vEry gooD aPprOXimatIon Of $M(B,T)$ Curves. hEre $W$ is Equal to $25$, $25^2$, $25^3$, $25^4$ foR a SIngle iON, pair, trIPLeT or qUartet, respectivelY. in EAch case, the condItion oF $E_{K_{MaX}} - e_{K=0} \gg 30 K_{\maThrm{B}}T $ was fUlLfillED, ensuriNG tHAT RelATive error in clAculation of $m(b,T)$ iS practIcAllY Zero. $E_{k=0}$ Is the GrOUnd State energy And $E_{K_{max}}$ is the MaximaL Energy oF The calcUlated ExcIteD staTEs.

maGneTiC SimULaTioNS
====================

![\[FiG:MagnetiZaTiOn\] MagNetiZATIOn peR onE ion As a fuNction of the maGneTic fIEld $b$ of diFfereNt Mn$^{3+}$ MaGnetiC clustErs in gan at $T$=2 K obtained usIng cRystal fieLd mOdEl wItH ferrOMagnetIc sUpeRexchanGe couplINg $J$ = 2 Mev. tHE mAgnetic easy axis $M_{\peRp}$ (SOLiD lines) is PerpenDIcUlAR to the $\teXtBf{c}$ Axis OF gaN, whEreaS ThE hard one $m_{||}$ (dash lINeS) iS parallEl To the $\tExTbf{C}$ axIs. (a) ThE MediUm magnEtic fielD, and (b) THe high magnetic FIeld region.](MagNEtIZup4iOns.pDf){wIdth="8.6"}

![\[Fig:OpenvsClOSed\] COmpaRIsIon OF magnEtizaTiON pER one ion of magnetic clUsTers coMposeD of three ions (oPen vs closeD TRIplet) and Four IOnS (String vs tetrahEdron Quartet)](OpeNvsClosed.Pdf){wiDth="8.6"}

In Fig. \[fig:MagnetIZAtion\] we pResEnt The ResULTs Of our simulatiONS of mAgNetizatIon Per one iOn aS a fUncTioN oF magnetic Field $B$ of DiFfErEnT maGnetiC Clusters At $t$=2 K. WE sEe tHat $M(B)$ VAries sHarplY witH mAgNEtiC field aND tHE SatuRaTiOn is ObsErVed foR higH $b$. ThE magnetIc easy axiS ($M_{\pERp}$) is PeRpEndiculAr to
B,C}\sum_{j_3=A,B,C}$.

In  order tospeed  up  th ecalc ulat ions the paral l eliz ation of the code is u sed.Ad d itio n al ly, d ue to t h ep r ese nc eofth e  B oltzm anfactorin Eq. (\[ eq: M\ _cf\_Center\ ] ),  only the$k_ {max} < W$ l owe st ene rg y l e velsand  eige nfunct i ons ne ed to beca l culate d  in ord e r  t o ob tain a very gooda pp r oximation of $ M(B,T) $c ur v e s.Her e $W$ is e qu al to $25$, $ 2 5^ 2 $ ,  $2 5 ^3$, $25^4$ f or a single ion , pair ,tri p let or  quar te t , r espectively . In  each cas e, the conditi o n of $E _{k_{m ax} } -  E_{ k =0 }\gg  3 0  k_ { \m ath r m{B }}T $ wa sfu llfil led, e n s urin g t hatrelat ive error incla cula t ion  of $ M(B,T )$ i spract ically  zero .$E_{k=0}$ is th e gr ound stat e e ne rgy  a nd $E _ {k_{ma x}} $ i s the m aximale ner gy o f  t he calculated exci te d st ates.

M agneti c  s im u lations== === ==== = = ===== ==== 
! [\[Fig:M agneti z at io n\] Mag ne tizati on  pe r o ne io n  asa func tion ofthe m a gnetic field $ B $ of differen t  M n $ ^{ 3 +}$mag netic clust ersi n Ga N at $T $=2 K obt ained  u s in g  crystal field mode lwith f errom agnetic super exchange c o u p ling $J$  = 2 me V . The magnetic  easy  axis $M_{ \ perp}$ ( solid  lines)is perpen d i cular to  th e $ \te xtb f { c} $ axis of GaN , wher ea s the h ard  one $M _{| |}$  (d ash  l ines) isparallel  t oth e$\t extbf { c}$ axis .(a)  T hemediu m  magne tic f ield ,an d  (b ) the h i gh m agne ti cfiel d r eg ion.] (Mag n eti zUp4Ion s.pdf){wi dth = "8.6 "} 
![\[Fig :OpenVsClosed \]  Comparisi on  of  magne t i zation p er one ion of magneticc lusters  co mpose d of  three io ns(openvsc losedtriple t) an dfou r ions( s tr ing  v s tetrahed r o n q uarte t) ](Op enVsClo sed.pdf){width="8. 6 "}
In Fig. \[Fi g:M agne t i za tio n \] wepr e sen t the results ofour simula ti o ns  of magnet i zat io n per o ne ionas af unction  of magne tic field  $ B$ o f dif ferent mag netic cl usters at $T$=2 K.  We s eethat $ M( B)$  vari es sha r ply  with  magne ti c fiel d and  t he satur ation is observed for h igh $B $. Th e m agnetic e asy axi s ($M_{\p erp} $) is perp end icu lar t o
B,C}\sum_{j_3=A,B,C}$.

In order_to speed_up the calculations the_parallelization of_the_code is_used._Additionally, due to_the presence of_the Boltzman factor in_Eq. (\[eq:M\_cf\_Center\]), only the_$k_{max}_< W$ lowest energy levels and eigenfunctions need to be calculated in order to_obtain_a very_good_approximation_of $M(B,T)$ curves. Here $W$_is equal to $25$, $25^2$,_$25^3$, $25^4$_for a single ion, pair, triplet or quartet,_respectively._In each case,_the condition of $E_{k_{max}} - E_{k=0} \gg 30 k_{\mathrm{B}}T_$ was fullfilled, ensuring that relative_error in claculation_of_$M(B,T)$_is practically zero. $E_{k=0}$_is the ground state energy and_$E_{k_{max}}$ is the maximal energy of_the calculated excited states.

Magnetic simulations
====================

![\[Fig:Magnetization\] Magnetization per_one ion as a function of_the magnetic field $B$ of_different Mn$^{3+}$_magnetic clusters in GaN at_$T$=2 K obtained using_crystal field_model with ferromagnetic_superexchange coupling $J$ = 2 meV. The_magnetic easy axis_$M_{\perp}$ (solid lines) is perpendicular to_the_$\textbf{c}$ axis of_GaN,_whereas_the hard_one $M_{||}$ (dash_lines)_is parallel_to_the $\textbf{c}$ axis. (a) The medium_magnetic_field, and (b) the high magnetic field_region.](MagnetizUp4Ions.pdf){width="8.6"}

![\[Fig:OpenVsClosed\] Comparision of magnetization_per_one ion of magnetic_clusters composed of three ions_(open vs closed triplet) and four_ions (string_vs tetrahedron_quartet)](OpenVsClosed.pdf){width="8.6"}

In Fig. \[Fig:Magnetization\] we present the results of our simulations of magnetization_per one ion as a function_of magnetic field $B$_of different_magnetic_clusters at $T$=2 K._We_see that_$M(B)$ varies sharply with magnetic field and_the saturation_is observed for high $B$. The_magnetic easy axis ($M_{\perp}$)_is_perpendicular to