from NanoTCAD_ViDES import * Lx=1e-9; Ly=2e-9; a=0.5e-9; nx=int(Lx/a)+1; ny=int(Ly/a)+1; me=(0.98*0.19*0.19)**(1/3)*m0; t=hbar**2/(2*me*a**2)/q; print nx,ny # The grid is ordered as follows # # 1--4--7--10--13 # | | | | | # 2--5--8--11--14 # | | | | | # 3--6--9--12--15 # I define the Hamiltonian H=Hamiltonian(nx,ny); #Initialization of the list h=[]; #I consider only one orbital h.append([1,0,0]); #I work on the on-site energies for j in range(0,ny): for i in range(0,nx): ix=i+j*nx; h.append([ix+1,ix+1,4*t]); #I work on the hopping parameters #First and second off-diagonal for j in range(0,ny): for i in range(0,nx): ix=i+j*nx; if (i