update

academic_codes/2021.03.01_Wilson_loop_in_SSH_model/Wilson_loop_in_SSH_model.py

@@ -11,7 +11,7 @@ from math import *
 def hamiltonian(k):  # SSH模型哈密顿量
     gamma = 0.5
     lambda0 = 1
-    h = np.zeros((2, 2))*(1+0j)
+    h = np.zeros((2, 2), dtype=complex)
     h[0,0] = 0
     h[1,1] = 0
     h[0,1] = gamma+lambda0*cmath.exp(-1j*k)

academic_codes/2021.03.01_Wilson_loop_in_SSH_model/Wilson_loop_in_SSH_model_with_better_code.py

@@ -12,7 +12,7 @@ def hamiltonian(k):
     gamma = 0.5
     lambda0 = 1
     delta = 0
-    h = np.zeros((2, 2))*(1+0j)
+    h = np.zeros((2, 2), dtype=complex)
     h[0,0] = delta
     h[1,1] = -delta
     h[0,1] = gamma+lambda0*cmath.exp(-1j*k)

academic_codes/2022.04.21_Berry_curvature/Berry_curvature_distribution_of_graphene_under_broken_inversion_symmery_with_Kubo_formula.py

@@ -10,7 +10,7 @@ import cmath
 
 
 def hamiltonian(k1, k2, t1=2.82, a=1/sqrt(3)):  # 石墨烯哈密顿量（a为原子间距，不赋值的话默认为1/sqrt(3)）
-    h = np.zeros((2, 2))*(1+0j)
+    h = np.zeros((2, 2), dtype=complex)
     h[0, 0] = 0.28/2
     h[1, 1] = -0.28/2
     h[1, 0] = t1*(cmath.exp(1j*k2*a)+cmath.exp(1j*sqrt(3)/2*k1*a-1j/2*k2*a)+cmath.exp(-1j*sqrt(3)/2*k1*a-1j/2*k2*a))

academic_codes/2022.04.21_Berry_curvature/Berry_curvature_distribution_of_graphene_under_broken_inversion_symmery_with_Wilson_loop.py

@@ -10,7 +10,7 @@ import cmath
 
 
 def hamiltonian(k1, k2, t1=2.82, a=1/sqrt(3)):  # 石墨烯哈密顿量（a为原子间距，不赋值的话默认为1/sqrt(3)）
-    h = np.zeros((2, 2))*(1+0j)
+    h = np.zeros((2, 2), dtype=complex)
     h[0, 0] = 0.28/2
     h[1, 1] = -0.28/2
     h[1, 0] = t1*(cmath.exp(1j*k2*a)+cmath.exp(1j*sqrt(3)/2*k1*a-1j/2*k2*a)+cmath.exp(-1j*sqrt(3)/2*k1*a-1j/2*k2*a))

academic_codes/2022.04.21_Berry_curvature/Berry_curvature_distribution_of_graphene_under_broken_inversion_symmery_with_the_efficient_method.py

@@ -10,7 +10,7 @@ import cmath
 
 
 def hamiltonian(k1, k2, t1=2.82, a=1/sqrt(3)):  # 石墨烯哈密顿量（a为原子间距，不赋值的话默认为1/sqrt(3)）
-    h = np.zeros((2, 2))*(1+0j)
+    h = np.zeros((2, 2), dtype=complex)
     h[0, 0] = 0.28/2
     h[1, 1] = -0.28/2
     h[1, 0] = t1*(cmath.exp(1j*k2*a)+cmath.exp(1j*sqrt(3)/2*k1*a-1j/2*k2*a)+cmath.exp(-1j*sqrt(3)/2*k1*a-1j/2*k2*a))

academic_codes/2022.08.13_Berry_curvature_distribution_in_function_form/Berry_curvature_distribution_with_Wilson_loop_(function_form).py

@@ -10,7 +10,7 @@ import math
 
 
 def hamiltonian(k1, k2, t1=2.82, a=1/sqrt(3)):  # 石墨烯哈密顿量（a为原子间距，不赋值的话默认为1/sqrt(3)）
-    h = np.zeros((2, 2))*(1+0j)
+    h = np.zeros((2, 2), dtype=complex)
     h[0, 0] = 0.28/2
     h[1, 1] = -0.28/2
     h[1, 0] = t1*(cmath.exp(1j*k2*a)+cmath.exp(1j*sqrt(3)/2*k1*a-1j/2*k2*a)+cmath.exp(-1j*sqrt(3)/2*k1*a-1j/2*k2*a))

academic_codes/2022.08.13_Berry_curvature_distribution_in_function_form/Berry_curvature_distribution_with_Wilson_loop_for_degenerate_case_(function_form).py

@@ -10,7 +10,7 @@ import math
 
 
 def hamiltonian(k1, k2, t1=2.82, a=1/sqrt(3)):  # 石墨烯哈密顿量（a为原子间距，不赋值的话默认为1/sqrt(3)）
-    h = np.zeros((2, 2))*(1+0j)
+    h = np.zeros((2, 2), dtype=complex)
     h[0, 0] = 0.28/2
     h[1, 1] = -0.28/2
     h[1, 0] = t1*(cmath.exp(1j*k2*a)+cmath.exp(1j*sqrt(3)/2*k1*a-1j/2*k2*a)+cmath.exp(-1j*sqrt(3)/2*k1*a-1j/2*k2*a))

academic_codes/2022.08.13_Berry_curvature_distribution_in_function_form/Berry_curvature_distribution_with_the_efficient_method_(function_form).py

@@ -9,7 +9,7 @@ import cmath
 import math
 
 def hamiltonian(k1, k2, t1=2.82, a=1/sqrt(3)):  # 石墨烯哈密顿量（a为原子间距，不赋值的话默认为1/sqrt(3)）
-    h = np.zeros((2, 2))*(1+0j)
+    h = np.zeros((2, 2), dtype=complex)
     h[0, 0] = 0.28/2
     h[1, 1] = -0.28/2
     h[1, 0] = t1*(cmath.exp(1j*k2*a)+cmath.exp(1j*sqrt(3)/2*k1*a-1j/2*k2*a)+cmath.exp(-1j*sqrt(3)/2*k1*a-1j/2*k2*a))

academic_codes/2022.08.13_Berry_curvature_distribution_in_function_form/Berry_curvature_distribution_with_the_efficient_method_for_degenerate_case_(function_form).py

@@ -10,7 +10,7 @@ import math
 
 
 def hamiltonian(k1, k2, t1=2.82, a=1/sqrt(3)):  # 石墨烯哈密顿量（a为原子间距，不赋值的话默认为1/sqrt(3)）
-    h = np.zeros((2, 2))*(1+0j)
+    h = np.zeros((2, 2), dtype=complex)
     h[0, 0] = 0.28/2
     h[1, 1] = -0.28/2
     h[1, 0] = t1*(cmath.exp(1j*k2*a)+cmath.exp(1j*sqrt(3)/2*k1*a-1j/2*k2*a)+cmath.exp(-1j*sqrt(3)/2*k1*a-1j/2*k2*a))