石墨烯纳米带超胞hubbard计算程序

子函数与前面程序相同，使用了超胞计算的技术python Nx=30; Ny=32; 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:

>>阅读原文<<