^+({\bf r}_{ij}){\sigma^{d}_{\hat\alpha}}^-({\bf r}_{jk})\end{aligned}$$ At this point, we are not interested for the possible orientations of the dimers with respect to the holes, as the hole density operator commutes with the dimer density operator. Now, we will repeatedly use that $$\begin{aligned}
\left[{\sigma^{h}}^{\pm},\sigma^{h}\right] = \mp 2{\sigma^{h}}^{\pm}\end{aligned}$$ at the same position in real space. At any distance different from zero, the commutator vanishes. For a given hole at a position ${\bf R}$, there are four possible positions ${\bf R}'={\bf R} \pm \hat x \pm\hat y$ where it may move through the available hopping term. By counting contributions from all these terms for every site of the lattice, we exactly take into account all the terms of the Hamiltonian including Hermitean conjugates.

If ${\bf R}$ is the initial hole position and ${\bf R}+{\bf r}_0$ the final one, then ${\bf r}_0$ has eight possible values: ${\bf r}_0=\pm
\hat x \pm \hat y $. The first commutator can now be computed for any site ${\bf R}$(The operator $\mathcal{T}_{\rm hole}({\bf R},{\bf r}_0)$ contains each of the above eight hopping terms). We have:

$$\begin{aligned}
\left[-t_{\rm hole}\mathcal{T}_{\rm hole}({\bf R},{\bf r}_0),{\sigma^{h}}({\bf k})\right] &=&
-t_{\rm hole}\mathcal{T}_{\rm hole}({\bf R},{\bf r}_0)\left(e^{-i{\bf k}\cdot{\bf R}}-e^{-i{\bf k}\cdot({\bf R}+{\bf r}_0)}\right)\\
\left[{\sigma^{h}}(-{\bf k}),-t_{\rm hole}\mathcal{T}_{\rm hole}({\bf R},{\bf r}_0)\right] &=&
 t_{\rm hole}\mathcal{T}_{\rm hole}({\bf R},{\bf r}_0)\left(e^{i{\bf k}\cdot{\bf R}}-e^{i{\bf k}\cdot({\bf R}+{\bf r}_0)}\right)\\
\left[{\sigma^{h}}(-{\bf k}),\left[-t_{\rm hole}\mathcal
^+({\bf r}_{ij}){\sigma^{d}_{\hat\alpha}}^-({\bf r}_{jk})\end{aligned}$$ At this point, we are not interested for the possible orientations of the dimer with obedience to the holes, as the hole concentration operator commute with the dimer density operator. immediately, we will repeatedly practice that $ $ \begin{aligned } 
 \left[{\sigma^{h}}^{\pm},\sigma^{h}\right ] = \mp 2{\sigma^{h}}^{\pm}\end{aligned}$$ at the same position in real space. At any distance different from zero, the commutator vanishes. For a give hole at a position $ { \bf R}$, there be four possible positions $ { \bf R}'={\bf R } \pm \hat x \pm\hat y$ where it may move through the available hopping term. By reckon contributions from all these terms for every site of the wicket, we exactly take into report all the terms of the Hamiltonian include Hermitean conjugates. 

 If $ { \bf R}$ is the initial hole position and $ { \bf R}+{\bf r}_0 $ the final one, then $ { \bf r}_0 $ has eight potential values: $ { \bf r}_0=\pm 
 \hat x \pm \hat y $. The first commutator can now be computed for any site $ { \bf R}$(The operator $ \mathcal{T}_{\rm hole}({\bf R},{\bf r}_0)$ contain each of the above eight hop terms). We have: 

 $ $ \begin{aligned } 
 \left[-t_{\rm hole}\mathcal{T}_{\rm hole}({\bf R},{\bf r}_0),{\sigma^{h}}({\bf k})\right ] & = & 
 -t_{\rm hole}\mathcal{T}_{\rm hole}({\bf R},{\bf r}_0)\left(e^{-i{\bf k}\cdot{\bf R}}-e^{-i{\bf k}\cdot({\bf R}+{\bf r}_0)}\right)\\ 
 \left[{\sigma^{h}}(-{\bf k}),-t_{\rm hole}\mathcal{T}_{\rm hole}({\bf R},{\bf r}_0)\right ] & = & 
  t_{\rm hole}\mathcal{T}_{\rm hole}({\bf R},{\bf r}_0)\left(e^{i{\bf k}\cdot{\bf R}}-e^{i{\bf k}\cdot({\bf R}+{\bf r}_0)}\right)\\ 
 \left[{\sigma^{h}}(-{\bf k}),\left[-t_{\rm hole}\mathcal
^+({\bf g}_{ij}){\sigma^{d}_{\hat\alpha}}^-({\bf r}_{jk})\ekd{aligned}$$ At this point, we ars not ingerested for the possible ormentqtionw of the dimers with rdspect to the holws, aw the hole vsnsity operatkv comkntes with the dlmer densitf operator. Now, wd cill repeatedly use that $$\begin{aligneq}
\left[{\sibmw^{h}}^{\pm},\sigma^{h}\righj] = \mp 2{\figmz^{h}}^{\pm}\end{aligned}$$ at the same positioh in reel space. At any distance different from zfro, hhe commutator vanlshes. For a givqb hole at a oosition ${\bf R}$, there arg four possible positions ${\bf R}'={\bf R} \pm \hat x \pm\hqt y$ avere it may move through the availabne hoppong term. By comntinj cobtributions from all vhese terms for everi site of dhz lattice, we exactly raje injo acwoung alu tge tsrms ov tie Hamiltonjan includibg Hermitean conjugstqw.

If ${\bf R}$ is tge inieiwl hole position and ${\bf R}+{\bf r}_0$ the final ons, then ${\bf r}_0$ has eight pissible values: ${\bf r}_0=\pm
\jat x \pm \rat y $. The first commutator can now be computed fmr anb ritt ${\nn R}$(Tfw lperator $\mathcal{T}_{\rm hole}({\bf R},{\bf r}_0)$ contains eacr og nhe above eight hjpping termx). Ae rave:

$$\begin{aligved}
\lefc[-f_{\rj hole}\mathcal{T}_{\rm hlle}({\bf R},{\ff r}_0),{\sugma^{h}}({\bf k})\wighy] &=&
-t_{\rm hole}\mathcal{T}_{\rm hole}({\bf R},{\bf r}_0)\left(e^{-i{\yf j}\cdot{\bf R}}-e^{-i{\bf k}\cdoc({\bf R}+{\bf r}_0)}\rigkt)\\
\left[{\xigma^{n}}(-{\bf k}),-t_{\rm hole}\mathcal{T}_{\rm hols}({\bf R},{\bf r}_0)\rihht] &=&
 t_{\rm gule}\mathcal{T}_{\rm houe}({\bn R},{\tf r}_0)\left(t^{k{\bf k}\cdot{\bf R}}-e^{i{\bf k}\cdot({\bf C}+{\bf r}_0)}\xight)\\
\lefg[{\sigka^{h}}(-{\bf h}),\left[-t_{\rm hlle}\mabvcal
^+({\bf r}_{ij}){\sigma^{d}_{\hat\alpha}}^-({\bf r}_{jk})\end{aligned}$$ At this point, we interested the possible of the dimers as hole density operator with the dimer operator. Now, we will repeatedly use $$\begin{aligned} \left[{\sigma^{h}}^{\pm},\sigma^{h}\right] = \mp 2{\sigma^{h}}^{\pm}\end{aligned}$$ at the same position in real space. At distance different from zero, the commutator vanishes. For a given hole at a ${\bf there four positions ${\bf R}'={\bf R} \pm \hat x \pm\hat y$ where it may move through the available term. By counting contributions from all these terms every site of the we exactly take into account the of the including conjugates. ${\bf R}$ is initial hole position and ${\bf R}+{\bf r}_0$ the final one, then ${\bf r}_0$ has eight possible values: r}_0=\pm \hat \hat y The commutator now be computed site ${\bf R}$(The operator $\mathcal{T}_{\rm hole}({\bf each of the above eight hopping terms). We $$\begin{aligned} \left[-t_{\rm hole}({\bf R},{\bf r}_0),{\sigma^{h}}({\bf k})\right] &=& -t_{\rm hole}({\bf R},{\bf r}_0)\left(e^{-i{\bf k}\cdot{\bf R}}-e^{-i{\bf k}\cdot({\bf R}+{\bf r}_0)}\right)\\ k}),-t_{\rm hole}\mathcal{T}_{\rm hole}({\bf R},{\bf r}_0)\right] &=& t_{\rm hole}\mathcal{T}_{\rm hole}({\bf R},{\bf r}_0)\left(e^{i{\bf k}\cdot{\bf R}}-e^{i{\bf k}\cdot({\bf R}+{\bf k}),\left[-t_{\rm hole}\mathcal
^+({\bf r}_{ij}){\sigma^{d}_{\hat\alpha}}^-({\bf r}_{jk})\eNd{aligned}$$ AT this PoiNt, wE aRe noT intErested for the pOSsibLe orientations of the dimErs wiTh REspeCT tO the hOles, as tHE hOLE deNsItY opErAToR commUteS with thE dimer densIty OpErator. Now, we wILl Repeatedly Use That $$\begin{aliGneD}
\left[{\sIgMa^{h}}^{\PM},\sigmA^{h}\rIght] = \mP 2{\sigma^{H}}^{\Pm}\end{aLigned}$$ at tHe SAme posITion in rEAL sPace. at any distance diffEReNT from zero, the coMmutatOr VAnISHes. for A given hole At A posiTIon ${\bf R}$, tHErE ARE foUR possible posiTions ${\bf R}'={\bf R} \PM \haT x \pm\haT y$ WheRE it may Move tHrOUgh The availablE hopPing term. BY countINg contrIButions From alL thEse TermS FoR eVerY sITe oF ThE laTTicE, we exactLy TaKe intO accOUNT All tHe tErms Of the hamiltonian inCluDing hErmItean ConjuGateS.

IF ${\bf R}$ iS the inItial HoLe position and ${\bf r}+{\bf r}_0$ The final oNe, tHeN ${\bf R}_0$ hAs eigHT possiBle ValUes: ${\bf r}_0=\pM
\hat x \pm \HAt y $. thE FIRsT commutator can now bE cOMPuTed for anY site ${\bF r}$(THe OPerator $\mAtHcaL{T}_{\rm HOLe}({\bf R},{\Bf r}_0)$ cONtAins each Of the aBOvE eIght hopPiNg termS). WE haVe:

$$\bEgin{aLIgneD}
\left[-t_{\Rm hole}\maThcal{t}_{\Rm hole}({\bf R},{\bf r}_0),{\siGMa^{h}}({\bf k})\right] &=&
-t_{\rM HoLE}\MaTHcal{t}_{\rm Hole}({\bf R},{\bf r}_0)\lEft(e^{-I{\Bf k}\cDot{\bF r}}-e^{-I{\bf K}\Cdot({\bF R}+{\bf r}_0)}\RiGHt)\\
\LEft[{\sigma^{h}}(-{\bf k}),-t_{\rm hole}\MaThcal{T}_{\Rm holE}({\bf R},{\bf r}_0)\right] &=&
 t_{\Rm hole}\mathCAL{t}_{\rm hole}({\bF R},{\bf R}_0)\LeFT(e^{i{\bf k}\cdot{\bf R}}-e^{I{\bf k}\cDot({\bf R}+{\bf r}_0)}\rIGht)\\
\left[{\sIgma^{h}}(-{\Bf k}),\left[-t_{\Rm hole}\matHCAl
^+({\bf r}_{ij}){\sigma^{d }_{\hat\al pha}} ^-( {\b fr}_{ jk}) \end{aligned}$ $  Atthis point, we are not  inte re s tedf or  thepossibl e  o r i ent at io nsof th e dim ers  with r espect tothe  h oles, as the ho le density  op erator commu tes  withth e d i mer d ens ity o perato r . Now,  we willre p eatedl y  use th a t  $ $\be gin{aligned}
\lef t [{ \ sigma^{h}}^{\p m},\si gm a ^{ h } \ri ght ] = \mp 2{ \s igma^ { h}}^{\p m }\ e n d {al i gned}$$ at th e same posi t ion  in re al  sp a ce. At  anydi s tan ce differen t fr om zero,the co m mutator vanishe s. For  agiv en h o le  a t a  p o sit i on  ${ \ bfR}$, the re  a re fo ur p o s s i blepos itio ns ${ \bf R}'={\bfR}\pm\ hat  x \p m\hat  y$wh ere i t maymoveth rough the avail able  hoppingter m.  By  c ounti n g cont rib uti ons fro m all t h ese  t e r m sfor every site ofth e la ttice, w e exac t ly  t a ke intoac cou nt a l l  theterm s  o f the Ha milton i an  i ncludin gHermit ea n c onj ugate s .

I f ${\b f R}$ is  thei nitial hole po s ition and ${\ b fR } +{ \ bf r }_0 $ the final  one ,  the n ${ \ bf  r} _ 0$ ha s eig ht po s sible values: ${\bf  r }_0=\p m
\ha t x \pm \haty $. The f i r s t commut ator ca n  now be comput ed fo r any site ${\bf R} $(The  operato r $\mathc a l {T}_{\rm  ho le} ({\ bfR } ,{ \bf r}_0)$ co n t ains  e ach ofthe  aboveeig hthop pin gterms). W e have:
$ $\ be gi n{a ligne d }
\left[ -t _{\ rm  ho le}\m a thcal{ T}_{\ rm h ol e} ( {\b f R},{\ b fr } _0), {\ si gma^ {h} }( {\bfk})\ r igh t] &=&-t_{\rm h ole } \mat hc al {T}_{\r m hole}({\bfR} ,{\bf r}_0 )\ lef t(e^{- i { \bf k}\c dot{\bf R}}-e^{-i{\bf k } \cdot({ \bf  R}+{ \bfr}_0)}\ri ght )\\
\l eft [ {\sigm a^{h}} (-{\b fk}) , - t_{\r m ho le} \m athcal{T}_ { \ rmhole} ({ \bfR},{\bf  r}_0)\right] &=&t_{ \rm hole}\mat hca l{T} _ { \r m h o le } ({\ bf R}, { \ bf r}_0)\left(e ^{i{\bf k} \c d ot {\bf R}}-e ^ {i{ \b f k}\cd ot({\bf  R}+{ \ bf r}_0 )}\right) \\
\left[ {\ sigm a ^ {h} }(-{\bf k} ),\left[ -t_{\rm h o le}\m a th cal
^+({\bf r}_{ij}){\sigma^{d}_{\hat\alpha}}^-({\bf_r}_{jk})\end{aligned}$$ At_this point, we are_not interested_for_the possible_orientations_of the dimers_with respect to_the holes, as the_hole density operator_commutes_with the dimer density operator. Now, we will repeatedly use that $$\begin{aligned}
\left[{\sigma^{h}}^{\pm},\sigma^{h}\right] = \mp_2{\sigma^{h}}^{\pm}\end{aligned}$$_at the_same_position_in real space. At any_distance different from zero, the_commutator vanishes._For a given hole at a position ${\bf_R}$,_there are four_possible positions ${\bf R}'={\bf R} \pm \hat x \pm\hat_y$ where it may move through_the available hopping_term._By_counting contributions from all_these terms for every site of_the lattice, we exactly take into_account all the terms of the Hamiltonian_including Hermitean conjugates.

If ${\bf R}$ is_the initial hole position and_${\bf R}+{\bf_r}_0$ the final one, then_${\bf r}_0$ has_eight possible_values: ${\bf r}_0=\pm
\hat_x \pm \hat y $. The_first commutator can_now be computed for any site_${\bf_R}$(The operator $\mathcal{T}_{\rm_hole}({\bf_R},{\bf_r}_0)$ contains_each of the_above_eight hopping_terms)._We have:

$$\begin{aligned}
\left[-t_{\rm hole}\mathcal{T}_{\rm hole}({\bf R},{\bf r}_0),{\sigma^{h}}({\bf_k})\right]_&=&
-t_{\rm hole}\mathcal{T}_{\rm hole}({\bf R},{\bf r}_0)\left(e^{-i{\bf k}\cdot{\bf R}}-e^{-i{\bf_k}\cdot({\bf R}+{\bf r}_0)}\right)\\
\left[{\sigma^{h}}(-{\bf k}),-t_{\rm_hole}\mathcal{T}_{\rm_hole}({\bf R},{\bf r}_0)\right] &=&
_t_{\rm hole}\mathcal{T}_{\rm hole}({\bf R},{\bf r}_0)\left(e^{i{\bf_k}\cdot{\bf R}}-e^{i{\bf k}\cdot({\bf R}+{\bf r}_0)}\right)\\
\left[{\sigma^{h}}(-{\bf k}),\left[-t_{\rm_hole}\mathcal