I tried this way,too:
for i in range(len(Kpoints)):
H_tot=H00+H01*exp(-1j*Kpoints*3*a)+H01.T*exp(1j*Kpoints*3*a)
Ev[:,i]=np.linalg.eigvals(H_tot)
Ev=np.sort(Ev)
and still the same result
Hello again, and thank you very much
that day I didn't really understand what you mean by kx=Kx. I understood yesterday and changed my code this way:
for kx in Kpoints:
H_tot=H00+H01*exp(-1j*kx*3*a)+H01.T*exp(1j*kx*3*a)
ii=0
Ev[:,ii]=np.linalg.eigvals(H_tot)
Ev=np.sort(Ev)...
Obviously not, they're not the same. I think the problem is in the calculating of H_tot. the final matrix is different from the same matrix in Matlab. but I don't know how to correct it.
The procedure is,
First, I must calculate total Hamiltonian, i.e., H_tot for each kx. Then find the...
Hi, I have a Matlab code for calculating the band structure of graphene nanoribbon which is working fine, but I wanted to convert it to python and I've done it. I guess I made a mistake somewhere because the plot is totally wrong. Could someone check it for me, please?