 The constraint is imposed on the transverse magnetization $M$ in the form $$\label{Mconst}
M \equiv \sum_i S_i^N = const.,$$ where the Hamiltonian of the constrained model is given by $$\label{Hconst} 
{\cal H}_M = -J \sum_{\langle i j \rangle} 
\sum_{x=1}^{N-1} S_i^x S_j^x = {\cal H} + hM.$$ The critical temperature of the constrained model is a symmetric and monotonically decreasing function of the prescribed transverse magnetization $M$. The transverse field $h$ here plays the part of the chemical potential $\mu_i$ of impurities (see Sect. 1) and the transverse magnetization $M$ accordingly plays the part of the impurity concentration $n_i$. It is also possible to implement $O(N)$ spin models with impurities or vacancies with diffusion in order to mimic the situation discussed in [@Fisher68]. However, the fact that the Hamiltonians given by (\[H\]) and (\[Hconst\]) only require a single ’species’ with a single coupling constant leads to some simplifications in the algorithms. Note that the symmetric constrained model $(M = 0)$ becomes equivalent to the symmetric unconstrained model $(h = 0)$ for sufficiently large lattices. In particular, both versions of the symmetric model have the same $T_c$.

The Monte-Carlo algorithm is chosen as a hybrid scheme, where each hybrid Monte-Carlo step consists of 10 updates each of which can be one of the following: one Metropolis sweep of the whole lattice, one single cluster Wolff update [@Wolff89], or one overrelaxation update of the whole lattice [@CL94], where the latter can only be applied for $N \geq 3$. The Metropolis algorithm updates the lattice sequentially and works in the standard way for the unconstrained model. For the constrained model the constraint $M = const.$ is observed locally by applying a Kawasaki update dynamics for the $N^{th}$ components of the spins. For each lattice site $i$ a nearest neighbor site $j$ is chosen randomly and a random amount of the $N^{th}$ spin component is proposed for exchange such that $S_i^N + S_j^N$ remains constant. Then new spin components $(S_i^1,S_i^2,\dots,S_i^{N-1})$ are proposed and
The constraint is imposed on the transverse magnetization $ M$ in the form $ $ \label{Mconst } 
 M \equiv \sum_i S_i^N = const. ,$$ where the Hamiltonian of the constrained exemplar is render by $ $ \label{Hconst } 
 { \cal H}_M = -J \sum_{\langle i j \rangle } 
 \sum_{x=1}^{N-1 } S_i^x S_j^x = { \cal H } + hM.$$ The critical temperature of the constrained model is a symmetric and monotonically decreasing routine of the prescribed transverse magnetization $ M$. The transverse playing field $ h$ here play the part of the chemical electric potential $ \mu_i$ of impurity (see Sect.   1) and the cross magnetization $ M$ accordingly plays the part of the impurity concentration $ n_i$. It is besides possible to implement $ O(N)$ tailspin models with impurity or vacancy with diffusion in order to mimic the situation hash out in [ @Fisher68 ]. However, the fact that the Hamiltonians given by (\[H\ ]) and (\[Hconst\ ]) only require a single ’ coinage ’ with a single yoke constant leads to some simplifications in the algorithms. notice that the symmetric constrained model $ (M = 0)$ becomes equivalent to the symmetric unconstrained model $ (h = 0)$ for sufficiently large lattices. In particular, both versions of the symmetric model have the same $ T_c$. 

 The Monte - Carlo algorithm is chosen as a hybrid dodge, where each hybrid Monte - Carlo step consists of 10 update each of which can be one of the pursuit: one Metropolis sweep of the whole lattice, one single bunch Wolff update [ @Wolff89 ], or one overrelaxation update of the whole wicket [ @CL94 ], where the latter can only be applied for $ N \geq 3$. The Metropolis algorithm updates the wicket sequentially and works in the standard way for the unconstrained model. For the constrained model the constraint $ thousand = const.$ is observed locally by applying a Kawasaki update dynamics for the $ N^{th}$ components of the spins. For each lattice web site $ i$ a near neighbor site $ j$ is chosen randomly and a random amount of the $ N^{th}$ spin component is proposed for exchange such that $ S_i^N + S_j^N$ remain constant. Then new spin components $ (S_i^1,S_i^2,\dots, S_i^{N-1})$ are nominate and
 Thf constraint is imposed un the transverse magnxtizatikn $M$ in ghe form $$\label{Mconst}
M \equiv \dun_i S_i^B = const.,$$ where the Hamkltonian lf the cinstcained model is jjven by $$\label{Ggonst} 
{\ral H}_M = -J \sum_{\lakgle i j \ratgle} 
\sum_{x=1}^{N-1} S_i^x S_g^x = {\cal H} + hM.$$ The critical temperaturq of thr fonstrained moqel ps a sgmmetric and monotonically decreasjng funbtion of the presvribed transverse magnetizwtioj $M$. The transverse field $h$ hets pjqys the part of the chtmncal potentjal $\mu_i$ of impurities (see Sect. 1) xnd tke transverwe mahtetization $N$ accjrdingly planx the [art of the impurity gonceitrarion $n_i$. It is also povsible to implemenj $O(N)$ spin kobels with impurities ir vacatciev wigy dkffhsmon in orfer to mimic fhe situatiin discussed in [@Fisnew68]. However, the ract trae the Hamiltonians given by (\[H\]) and (\[Hconsu\]) onlg require a single ’specues’ with a single coukling conseant leads to some simplifications in the algoritvms. Nkge uhqt thd sjmmetric constrained model $(M = 0)$ becomes equivajsnu tp the symmetrig unconstrained mocep $(n = 0)$ for suffickently ladge lattices. In pagticulat, both versions of yhe symmetric model have thw same $T_c$.

The Nonte-Carlo algoritkm is chosen as s hybtid scheme, where each hvbrid Jonte-Carlo dtep consjrts of 10 updates dacm ox which can be one of the sollowing: one Metropulis sweep of the whlle lattice, one single cpustet Wolfx update [@Wllff89], or one overrelaxation updavx of the wholg ldttpce [@CL94], whzre thc latter can onjy be applied yor $N \gez 3$. Thd Metropolps algorivhm updates ehe lattice saeuentially aid works yn tye srandard day for the unvonstrainvd modwl. For the constralned okdel the constrcnnu $N = const.$ is obxerxed llcelly ff applying a Kawxsayo updxte dynamicw for the $N^{th}$ components mf tge spins. For each kabtice sitg $i$ a neawest neighbor site $j$ is chosen gandokly and s rwndom amount of the $N^{th}$ spin ckmponent ls iroposed for qxchqnge such thct $S_i^N + S_j^N$ remains constant. Then new spmn components $(S_i^1,S_i^2,\dots,S_u^{N-1})$ are proposed and
The constraint is imposed on the transverse in form $$\label{Mconst} \equiv \sum_i S_i^N of constrained model is by $$\label{Hconst} {\cal = -J \sum_{\langle i j \rangle} S_i^x S_j^x = {\cal H} + hM.$$ The critical temperature of the constrained is a symmetric and monotonically decreasing function of the prescribed transverse magnetization $M$. transverse $h$ plays part of the chemical potential $\mu_i$ of impurities (see Sect. 1) and the transverse magnetization $M$ plays the part of the impurity concentration $n_i$. is also possible to $O(N)$ spin models with impurities vacancies diffusion in to the discussed in [@Fisher68]. the fact that the Hamiltonians given by (\[H\]) and (\[Hconst\]) only require a single ’species’ with a coupling constant some simplifications the Note the symmetric constrained = 0)$ becomes equivalent to the $(h = 0)$ for sufficiently large lattices. In both versions the symmetric model have the same The Monte-Carlo algorithm is chosen as a hybrid where each hybrid Monte-Carlo step consists of 10 updates each of which can be one following: one Metropolis sweep the whole lattice, single Wolff [@Wolff89], one overrelaxation of the whole lattice [@CL94], where the latter can only be for $N \geq 3$. The Metropolis algorithm updates the lattice works the standard way the unconstrained model. For constrained the constraint $M = observed by update for $N^{th}$ components of the For each lattice site $i$ nearest neighbor site $j$ random amount of the $N^{th}$ spin component is for exchange such that $S_i^N + S_j^N$ constant. Then new spin components $(S_i^1,S_i^2,\dots,S_i^{N-1})$ are proposed and
 The constraint is imposed on tHe transverSe magNetIzaTiOn $M$ iN the Form $$\label{MconsT}
m \equIv \sum_i S_i^N = const.,$$ where the hamilToNIan oF ThE consTrained MOdEL Is gIvEn By $$\lAbEL{HConst} 
{\Cal h}_M = -J \sum_{\lAngle i j \ranGle} 
\SuM_{x=1}^{N-1} S_i^x S_j^x = {\cal h} + HM.$$ the criticaL teMperature of tHe cOnstraInEd mODel is A syMmetrIc and mONotoniCally decrEaSIng funCTion of tHE PrEscrIbed transverse magNEtIZation $M$. The tranSverse FiELd $H$ HEre PlaYs the part oF tHe cheMIcal potENtIAL $\Mu_i$ OF impurities (seE Sect. 1) and the TRanSverse MaGneTIzatioN $M$ accOrDIngLy plays the pArt oF the impurIty conCEntratiON $n_i$. It is Also poSsiBle To imPLeMeNt $O(n)$ sPIn mODeLs wITh iMpuritieS oR vAcancIes wITH DIffuSioN in oRder tO mimic the situAtiOn diSCusSed in [@fisheR68]. HowEvEr, the Fact thAt the haMiltonians given By (\[H\]) aNd (\[Hconst\]) oNly ReQuiRe A singLE ’speciEs’ wIth A single CouplinG ConStANT LeAds to some simplificAtIONs In the algOrithmS. noTe THat the syMmEtrIc coNSTrainEd moDEl $(m = 0)$ becomes EquivaLEnT tO the symMeTric unCoNstRaiNed moDEl $(h = 0)$ fOr suffIciently Large LAttices. In partiCUlar, both versiONs OF ThE SymmEtrIc model have The sAMe $T_c$.

the MONtE-CaRLo algOrithM iS ChOSen as a hybrid scheme, wHeRe each HybriD Monte-Carlo stEp consists OF 10 UPdates eaCh of WHiCH can be one of the FolloWing: one MetROpolis swEep of The whole Lattice, onE SIngle cluSteR WoLff UpdATE [@WOlff89], or one overRELaxaTiOn updatE of The wholE laTtiCe [@Cl94], whErE the latteR can only Be ApPlIeD foR $N \geq 3$. tHe MetropOlIs aLgOriThm upDAtes thE lattIce sEqUeNTiaLly and wORkS IN the StAnDard Way FoR the uNconSTraIned modEl. For the cOnsTRainEd MoDel the cOnstraint $M = conSt.$ Is observed LoCalLy by apPLYing a KawAsaki update dynamics for tHE $N^{th}$ comPonEnts oF the Spins. For eAch LatticE siTE $i$ a neaRest neIghboR sIte $J$ IS chosEN RaNdoMlY and a randoM AMouNt of tHe $n^{th}$ sPin compOnent is proposed for EXchAnge such that $S_I^N + S_J^N$ reMAInS coNStANt. THeN New SPIn components $(S_i^1,S_I^2,\dots,S_i^{N-1})$ arE pROpOsed and
 The constraint is imposed  on the tr ansve rse  ma gn etiz atio n $M$ in the f o rm $ $\label{Mconst}
M \equ iv \s um _ i S_ i ^N  = co nst.,$$ wh e r e t he  H ami lt o ni an of  th e const rained mod elis  given by $$ \ la bel{Hconst } {\cal H}_M =  -J  \sum_ {\ lan g le ij \ rangl e} 
\s u m_{x=1 }^{N-1} S _i ^ x S_j^ x  = {\ca l H}  + h M.$$ The critical te m perature of th e cons tr a in e d  mo del  is a symm et ric a n d monot o ni c a l lyd ecreasing fun ction of th e  pr escrib ed  tr a nsvers e mag ne t iza tion $M$. T he t ransverse  field $h$ her e  playsthe pa rtofthec he mi cal  p o ten t ia l $ \ mu_ i$ of im pu ri ties(see S e c t. 1 ) a nd t he tr ansverse magn eti zati o n $ M$ ac cordi ngly  p laysthe pa rt of  t he impurity con cent ration $n _i$ .Itis  also possib letoimpleme nt $O(N ) $ s pi n m od els with impuritie so r  v acancies  withd if fu s ion in o rd erto m i m ic th e si t ua tion dis cussed in  [ @Fisher 68 ]. How ev er,  th e fac t  tha t theHamilton iansg iven by (\[H\] )  and (\[Hcons t \] ) on l y re qui re a single  ’sp e cies ’ wi t ha s i nglecoupl in g  c o nstant leads to som esimpli ficat ions in the a lgorithms. N o te thatthes ym m etric constrai ned m odel $(M = 0)$ beco mes e quivalen t to thes y mmetricunc ons tra ine d mo del $(h = 0)$ f or s uf ficient lylarge l att ice s.Inpa rticular,  both ve rs io ns  o f t he sy m metric m od elha vethe s a me $T_ c$.
TheMo nt e -Ca rlo alg o ri t h m is  c ho senasahybri d sc h eme , where  each hyb rid Mont e- Ca rlo ste p consists of  1 0 updatesea chof whi c h  can beone of the following: o n e Metro pol is sw eepof the wh ole  latti ce, one si ngle c luste rWol f f  upda t e  [ @Wo lf f89], or o n e  ov errel ax atio n updat e of the whole lat t ice  [@CL94], whe rethel a tt erc an onl yb e a p p lied for $N \ge q 3$. TheMe t ro polis algo r ith mupdates  the la ttice sequent ially and  works in  t he s t a nda rd way for  the unc onstraine d  mode l .For t heconstr ai ned  mode l thec ons train t $M =  c onst.$  is o bs erved lo cally by applying a Kaw asakiupdat e d ynamics f ort he$N^{th}$comp onents ofthe  sp ins.For eachlatt i ce  si t e $i$  a n e arest nei g hb ors i te  $j$ is cho s e n  ra ndoml y a n d a ra ndom  amount of the $N ^ {th}$ spin com pone n t  is  pr o pose dfor exchange s uch  t h a t $S_i^N  +  S_j^N$ rem ains con st a nt. T hen ne w spin  compon e n ts $(S_i^ 1,S_ i^2 ,\dots,S_ i^{ N- 1 })$ are  p ro p osed a nd
 The_constraint is_imposed on the transverse_magnetization $M$_in_the form_$$\label{Mconst}
M_\equiv \sum_i S_i^N_= const.,$$ where_the Hamiltonian of the_constrained model is_given_by $$\label{Hconst} 
{\cal H}_M = -J \sum_{\langle i j \rangle} 
\sum_{x=1}^{N-1} S_i^x S_j^x =_{\cal_H} +_hM.$$_The_critical temperature of the constrained_model is a symmetric and_monotonically decreasing_function of the prescribed transverse magnetization $M$. The_transverse_field $h$ here_plays the part of the chemical potential $\mu_i$ of_impurities (see Sect. 1) and the transverse_magnetization $M$ accordingly_plays_the_part of the impurity_concentration $n_i$. It is also possible_to implement $O(N)$ spin models with_impurities or vacancies with diffusion in order_to mimic the situation discussed in_[@Fisher68]. However, the fact that_the Hamiltonians_given by (\[H\]) and (\[Hconst\])_only require a_single ’species’_with a single_coupling constant leads to some simplifications_in the algorithms._Note that the symmetric constrained model_$(M_= 0)$ becomes_equivalent_to_the symmetric_unconstrained model $(h_=_0)$ for_sufficiently_large lattices. In particular, both versions_of_the symmetric model have the same $T_c$.

The_Monte-Carlo algorithm is chosen_as_a hybrid scheme, where_each hybrid Monte-Carlo step consists_of 10 updates each of which_can be_one of_the following: one Metropolis sweep of the whole lattice, one single_cluster Wolff update [@Wolff89], or one_overrelaxation update of the_whole lattice_[@CL94],_where the latter_can_only be_applied for $N \geq 3$. The Metropolis_algorithm updates_the lattice sequentially and works in_the standard way for_the_unconstrained model. For the constrained model_the constraint $M = const.$ is_observed locally by applying a_Kawasaki_update_dynamics for the $N^{th}$ components_of the spins. For each lattice_site $i$ a_nearest neighbor site $j$ is chosen randomly_and_a random amount of the $N^{th}$_spin_component is proposed for exchange such_that_$S_i^N_+ S_j^N$ remains constant. Then_new spin components $(S_i^1,S_i^2,\dots,S_i^{N-1})$ are proposed_and