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

Update calculation_of_Berry_curvature_and_Chern_number_by_another_method.py

Browse Source
This commit is contained in:
guanjihuan 2021-07-29 20:13:44 +08:00
parent 7b784f3888
commit 46dda84c9e
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.07.26_calculation_of_Berry_curvature_and_Chern_number_by_another_method/calculation_of_Berry_curvature_and_Chern_number_by_another_method.py
View File

@ -8,7 +8,7 @@ from math import *
import time import time
def hamiltonian(kx, ky): # 量子反常霍尔QAH模型（该参数对应的陈数为2） def hamiltonian(kx, ky): # one QAH model with Chern number = 2
t1 = 1.0 t1 = 1.0
t2 = 1.0 t2 = 1.0
t3 = 0.5 t3 = 0.5
@ -23,11 +23,11 @@ def hamiltonian(kx, ky): # 量子反常霍尔QAH模型（该参数对应的陈
def main(): def main():
start_time = time.time() start_time = time.time()
n = 200 n = 200 # integration
delta = 2*pi/n delta = 1e-6 # derivation
chern_number = 0 chern_number = 0
for kx in np.arange(-pi, pi, delta): for kx in np.arange(-pi, pi, 2*pi/n):
for ky in np.arange(-pi, pi,delta): for ky in np.arange(-pi, pi,2*pi/n):
H = hamiltonian(kx, ky) H = hamiltonian(kx, ky)
eigenvalue, eigenvector = np.linalg.eig(H) eigenvalue, eigenvector = np.linalg.eig(H)
vector_0 = eigenvector[:, np.argsort(np.real(eigenvalue))[0]] vector_0 = eigenvector[:, np.argsort(np.real(eigenvalue))[0]]