H=H0+HB+HI
H0=∑i,jJx(f+ij+fij)(f+i+1,j+1−fi+1,j+1)+Jy(f+ij+fij)(f+i+1,j−fi+1,j)+σJz(2f+ijfij−1)
(flux=1:σ=1,flux=−1:for i=0,j>0,σ=−1.)
HB=Jx(f+0,0+f0,0)(f+1,1−f1,1)(2χ+xχx−1)+Jy(f+0,0+f0,0)(f+1,0−f1,0)(2χ+yχy−1)+Jz(2f+0,0f0,0−1)(2χ+zχz−1)
HI=Sx(f+0+f0)(χ+x+χx)+Sy(f+0+f0)(χ+y+χy)+Sz(f+0+f0)(χ+z+χz)