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

Update Wilson_loop_in_SSH_model.py

Browse Source
This commit is contained in:
guanjihuan 2022-02-22 18:44:07 +08:00
parent 8d732bf85e
commit 74707df3fc
1 changed files with 3 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
academic_codes/2021.03.01_Wilson_loop_in_SSH_model/Wilson_loop_in_SSH_model.py
View File

@ -26,7 +26,7 @@ def main():
for k in k_array: for k in k_array:
vector = get_occupied_bands_vectors(k, hamiltonian) vector = get_occupied_bands_vectors(k, hamiltonian)
vector_array.append(vector) vector_array.append(vector)
# vector_array.append(vector*cmath.exp(1j*np.random.uniform(0, pi))) # vector_array.append(vector*cmath.exp(1j*np.random.uniform(0, pi))) # 随机相位测试
# 波函数固定一个规范 # 波函数固定一个规范
vector_sum = 0 vector_sum = 0
@ -38,12 +38,7 @@ def main():
# # 波函数固定一个规范 # # 波函数固定一个规范
# import guan # import guan
# vector_sum = 0 # vector_array = guan.find_vector_array_with_fixed_gauge_by_making_one_component_real(vector_array, precision=0.005)
# for i0 in range(Num_k):
# vector_sum += np.abs(vector_array[i0])
# index = np.argmax(np.abs(vector_sum))
# for i0 in range(Num_k):
# vector_array[i0] = guan.find_vector_with_fixed_gauge_by_making_one_component_real(vector_array[i0], index=index)
# 计算Wilson loop # 计算Wilson loop
W_k = 1 W_k = 1
@ -79,4 +74,4 @@ def find_vector_with_fixed_gauge_by_making_one_component_real(vector, precision=
if __name__ == '__main__': if __name__ == '__main__':
main() main()