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-01-04 15:04:15 +08:00
parent 41cc03bb72
commit 7b92011a54
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
academic_codes/2021.12.28_calculation_of_Chern_number_by_Wilson_loop/calculation_of_Chern_number_by_Wilson_loop.py
View File

@ -24,7 +24,7 @@ def hamiltonian(kx, ky): # 量子反常霍尔QAH模型（该参数对应的陈
def main(): def main():
start_time = time.time() start_time = time.time()
n = 200 # 积分密度 n = 100 # 积分密度
delta = 2*pi/n delta = 2*pi/n
chern_number = 0 chern_number = 0
for kx in np.arange(-pi, pi, delta): 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_3 = np.dot(vector_delta_kx_ky.transpose().conj(), vector_delta_ky)
line_4 = np.dot(vector_delta_ky.transpose().conj(), vector) 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 arg = np.log(np.dot(np.dot(np.dot(line_1, line_2), line_3), line_4))/1j
chern_number = chern_number + angle chern_number = chern_number + arg
chern_number = chern_number/(2*pi) chern_number = chern_number/(2*pi)
print('Chern number = ', chern_number) print('Chern number = ', chern_number)
end_time = time.time() end_time = time.time()