guanjihuan/guanjihuan.com
1
1
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Update calculation_of_Chern_number_by_Wilson_loop.py

Browse Source
This commit is contained in:
guanjihuan 2022-08-10 12:58:55 +08:00
parent 897d900c81
commit 794aad4a8e
1 changed files with 5 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
academic_codes/2021.12.28_calculation_of_Chern_number_by_Wilson_loop/calculation_of_Chern_number_by_Wilson_loop.py
View File

@ -32,7 +32,7 @@ def main():
for ky in np.arange(-pi, pi, delta): for ky in np.arange(-pi, pi, delta):
vector_array = [] vector_array = []
# line_1 # line_1
for i2 in range(n2+1): for i2 in range(n2):
H_delta = hamiltonian(kx+delta/n2*i2, ky) H_delta = hamiltonian(kx+delta/n2*i2, ky)
eigenvalue, eigenvector = np.linalg.eig(H_delta) eigenvalue, eigenvector = np.linalg.eig(H_delta)
vector_delta = eigenvector[:, np.argsort(np.real(eigenvalue))[0]] vector_delta = eigenvector[:, np.argsort(np.real(eigenvalue))[0]]
@ -40,19 +40,19 @@ def main():
vector_array.append(vector_delta) vector_array.append(vector_delta)
# line_2 # line_2
for i2 in range(n2): for i2 in range(n2):
H_delta = hamiltonian(kx+delta, ky+delta/n2*(i2+1)) H_delta = hamiltonian(kx+delta, ky+delta/n2*i2)
eigenvalue, eigenvector = np.linalg.eig(H_delta) eigenvalue, eigenvector = np.linalg.eig(H_delta)
vector_delta = eigenvector[:, np.argsort(np.real(eigenvalue))[0]] vector_delta = eigenvector[:, np.argsort(np.real(eigenvalue))[0]]
vector_array.append(vector_delta) vector_array.append(vector_delta)
# line_3 # line_3
for i2 in range(n2): for i2 in range(n2):
H_delta = hamiltonian(kx+delta-delta/n2*(i2+1), ky+delta) H_delta = hamiltonian(kx+delta-delta/n2*i2, ky+delta)
eigenvalue, eigenvector = np.linalg.eig(H_delta) eigenvalue, eigenvector = np.linalg.eig(H_delta)
vector_delta = eigenvector[:, np.argsort(np.real(eigenvalue))[0]] vector_delta = eigenvector[:, np.argsort(np.real(eigenvalue))[0]]
vector_array.append(vector_delta) vector_array.append(vector_delta)
# line_4 # line_4
for i2 in range(n2-1): for i2 in range(n2):
H_delta = hamiltonian(kx, ky+delta-delta/n2*(i2+1)) H_delta = hamiltonian(kx, ky+delta-delta/n2*i2)
eigenvalue, eigenvector = np.linalg.eig(H_delta) eigenvalue, eigenvector = np.linalg.eig(H_delta)
vector_delta = eigenvector[:, np.argsort(np.real(eigenvalue))[0]] vector_delta = eigenvector[:, np.argsort(np.real(eigenvalue))[0]]
vector_array.append(vector_delta) vector_array.append(vector_delta)