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

Update quantum_transport_of_square_lattice_under_magnetic_fields_in_multi_lead_systems_with_guan.py

Browse Source
This commit is contained in:
guanjihuan 2022-01-10 16:49:09 +08:00
parent 444daabd6e
commit 1e5988662f
1 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
academic_codes/2021.12.27_Chern_numbers_of_Landau_levels/quantum_transport_of_square_lattice_under_magnetic_fields_in_multi_lead_systems_with_guan.py
View File

@ -36,9 +36,9 @@ def get_center_hamiltonian(Nx, Ny, B):
def main(): def main():
start_time = time.time() start_time = time.time()
width = 30 width = 20
length = 80 length = 70
fermi_energy_array = np.arange(-4, 4, .02) fermi_energy_array = np.arange(-4, 4, .05)
# 中心区的哈密顿量 # 中心区的哈密顿量
H_cetner = get_center_hamiltonian(Nx=length, Ny=width, B=1/width) H_cetner = get_center_hamiltonian(Nx=length, Ny=width, B=1/width)
@ -68,7 +68,7 @@ def main():
h_lead4_to_center = np.zeros((width, width*length), dtype=complex) h_lead4_to_center = np.zeros((width, width*length), dtype=complex)
h_lead5_to_center = np.zeros((width, width*length), dtype=complex) h_lead5_to_center = np.zeros((width, width*length), dtype=complex)
h_lead6_to_center = np.zeros((width, width*length), dtype=complex) h_lead6_to_center = np.zeros((width, width*length), dtype=complex)
move = 10 move = 10 # the step of leads 2,3,6,5 moving to center
for i0 in range(width): for i0 in range(width):
h_lead1_to_center[i0, i0] = 1 h_lead1_to_center[i0, i0] = 1
h_lead2_to_center[i0, width*(move+i0)+(width-1)] = 1 h_lead2_to_center[i0, width*(move+i0)+(width-1)] = 1