From 7b92011a54277cc3414c8e661f0bca3125bc7da1 Mon Sep 17 00:00:00 2001 From: guanjihuan Date: Tue, 4 Jan 2022 15:04:15 +0800 Subject: [PATCH] Update calculation_of_Chern_number_by_Wilson_loop.py --- .../calculation_of_Chern_number_by_Wilson_loop.py | 6 +++--- 1 file changed, 3 insertions(+), 3 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 9db23a8..4b3094f 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 @@ -24,7 +24,7 @@ def hamiltonian(kx, ky): # 量子反常霍尔QAH模型（该参数对应的陈 def main(): start_time = time.time() - n = 200 # 积分密度 + n = 100 # 积分密度 delta = 2*pi/n chern_number = 0 for kx in np.arange(-pi, pi, delta): @@ -51,8 +51,8 @@ def main(): line_3 = np.dot(vector_delta_kx_ky.transpose().conj(), vector_delta_ky) line_4 = np.dot(vector_delta_ky.transpose().conj(), vector) - angle = np.log(np.dot(np.dot(np.dot(line_1, line_2), line_3), line_4))/1j - chern_number = chern_number + angle + arg = np.log(np.dot(np.dot(np.dot(line_1, line_2), line_3), line_4))/1j + chern_number = chern_number + arg chern_number = chern_number/(2*pi) print('Chern number = ', chern_number) end_time = time.time()