From 794aad4a8e9d4424d3367ca949051c40d6fdd917 Mon Sep 17 00:00:00 2001 From: guanjihuan Date: Wed, 10 Aug 2022 12:58:55 +0800 Subject: [PATCH] Update calculation_of_Chern_number_by_Wilson_loop.py --- .../calculation_of_Chern_number_by_Wilson_loop.py | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/academic_codes/2021.12.28_calculation_of_Chern_number_by_Wilson_loop/calculation_of_Chern_number_by_Wilson_loop.py b/academic_codes/2021.12.28_calculation_of_Chern_number_by_Wilson_loop/calculation_of_Chern_number_by_Wilson_loop.py index 137bdc9..ffd4bff 100644 --- a/academic_codes/2021.12.28_calculation_of_Chern_number_by_Wilson_loop/calculation_of_Chern_number_by_Wilson_loop.py +++ b/academic_codes/2021.12.28_calculation_of_Chern_number_by_Wilson_loop/calculation_of_Chern_number_by_Wilson_loop.py @@ -32,7 +32,7 @@ def main(): for ky in np.arange(-pi, pi, delta): vector_array = [] # line_1 - for i2 in range(n2+1): + for i2 in range(n2): H_delta = hamiltonian(kx+delta/n2*i2, ky) eigenvalue, eigenvector = np.linalg.eig(H_delta) vector_delta = eigenvector[:, np.argsort(np.real(eigenvalue))[0]] @@ -40,19 +40,19 @@ def main(): vector_array.append(vector_delta) # line_2 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) vector_delta = eigenvector[:, np.argsort(np.real(eigenvalue))[0]] vector_array.append(vector_delta) # line_3 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) vector_delta = eigenvector[:, np.argsort(np.real(eigenvalue))[0]] vector_array.append(vector_delta) # line_4 - for i2 in range(n2-1): - H_delta = hamiltonian(kx, ky+delta-delta/n2*(i2+1)) + for i2 in range(n2): + H_delta = hamiltonian(kx, ky+delta-delta/n2*i2) eigenvalue, eigenvector = np.linalg.eig(H_delta) vector_delta = eigenvector[:, np.argsort(np.real(eigenvalue))[0]] vector_array.append(vector_delta)