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

Update calculation_of_Chern_number_by_Wilson_loop_for_degenerate_case.py

Browse Source
This commit is contained in:
guanjihuan
2022-08-12 22:19:33 +08:00
parent a9c90d96b6
commit 8a61baf98e
1 changed files with 5 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

10
academic_codes/2022.08.12_calculation_of_Chern_number_by_Wilson_loop_for_degenerate_case/calculation_of_Chern_number_by_Wilson_loop_for_degenerate_case.py
View File

@@ -77,13 +77,13 @@ def calculate_chern_number_for_square_lattice_with_Wilson_loop_for_degenerate_ca
vector_delta = eigenvector[:, np.argsort(np.real(eigenvalue))] vector_delta = eigenvector[:, np.argsort(np.real(eigenvalue))]
vector_array.append(vector_delta) vector_array.append(vector_delta)
Wilson_loop = 1 Wilson_loop = 1
dim = len(num_of_bands) dim = len(index_of_bands)
for i0 in range(len(vector_array)-1): for i0 in range(len(vector_array)-1):
dot_matrix = np.zeros((dim , dim), dtype=complex) dot_matrix = np.zeros((dim , dim), dtype=complex)
i01 = 0 i01 = 0
for dim1 in num_of_bands: for dim1 in index_of_bands:
i02 = 0 i02 = 0
for dim2 in num_of_bands: for dim2 in index_of_bands:
dot_matrix[i01, i02] = np.dot(vector_array[i0][:, dim1].transpose().conj(), vector_array[i0+1][:, dim2]) dot_matrix[i01, i02] = np.dot(vector_array[i0][:, dim1].transpose().conj(), vector_array[i0+1][:, dim2])
i02 += 1 i02 += 1
i01 += 1 i01 += 1
@@ -91,9 +91,9 @@ def calculate_chern_number_for_square_lattice_with_Wilson_loop_for_degenerate_ca
Wilson_loop = Wilson_loop*det_value Wilson_loop = Wilson_loop*det_value
dot_matrix_plus = np.zeros((dim , dim), dtype=complex) dot_matrix_plus = np.zeros((dim , dim), dtype=complex)
i01 = 0 i01 = 0
for dim1 in num_of_bands: for dim1 in index_of_bands:
i02 = 0 i02 = 0
for dim2 in num_of_bands: for dim2 in index_of_bands:
dot_matrix_plus[i01, i02] = np.dot(vector_array[len(vector_array)-1][:, dim1].transpose().conj(), vector_array[0][:, dim2]) dot_matrix_plus[i01, i02] = np.dot(vector_array[len(vector_array)-1][:, dim1].transpose().conj(), vector_array[0][:, dim2])
i02 += 1 i02 += 1
i01 += 1 i01 += 1