I am trying to do some calculation using 2D dmrg code in Itensor.
2D dmrg code in ITENSOR includes nearest neighbor spin interactions.
for(auto bnd : lattice)
ampo += 0.5,"S+",bnd.s1,"S-",bnd.s2;
ampo += 0.5,"S-",bnd.s1,"S+",bnd.s2;
ampo += "Sz",bnd.s1,"Sz",bnd.s2;
Is there any way so that next nearest neighbor spin interaction can be included in 2D dmrg code for square lattice or triangular lattice?