鋸齒石墨烯磁性hubbard matlab計算程序

主程序python Nx=3; Ny=50; n=Nx*Ny; [x,y]=zigzagP(Ny,Nx); t1=-2.7; H=Hamiltonian_NN_graphene(x,y,t1); HD=H(n/3+1:2/3*n,n/3+1:2/3*n); HDL=H(n/3+1:2/3*n,1:n/3); HDR=H(n/3+1:2/3*n,2*n/3+1:n); x=x(n/3+1:2/3*n

>>阅读原文<<