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

Update example_rotation_of_degenerate_vectors.py

Browse Source
This commit is contained in:
guanjihuan 2022-07-14 14:05:15 +08:00
parent 88aba53d92
commit e89a3908c7
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/2022.07.12_gauge_fixing/example_rotation_of_degenerate_vectors.py
View File

@ -34,7 +34,7 @@ def rotation_of_degenerate_vectors(vector1, vector2, index1=None, index2=None, p
break break
return vector1, vector2 return vector1, vector2
def hamiltonian_of_BBH_model(kx, ky, gamma_x=0.5, gamma_y=0.5, lambda_x=1, lambda_y=1): def hamiltonian_of_bbh_model(kx, ky, gamma_x=0.5, gamma_y=0.5, lambda_x=1, lambda_y=1):
# label of atoms in a unit cell # label of atoms in a unit cell
# (2) —— (0) # (2) —— (0)
# | | # | |
@ -52,14 +52,14 @@ def hamiltonian_of_BBH_model(kx, ky, gamma_x=0.5, gamma_y=0.5, lambda_x=1, lambd
# For kx=0 # For kx=0
print('\nFor kx=0:\n') print('\nFor kx=0:\n')
eigenvalue, eigenvector = np.linalg.eigh(hamiltonian_of_BBH_model(kx=0, ky=0)) eigenvalue, eigenvector = np.linalg.eigh(hamiltonian_of_bbh_model(kx=0, ky=0))
print(eigenvalue, '\n') print(eigenvalue, '\n')
print(eigenvector[:, 0]) print(eigenvector[:, 0])
print(eigenvector[:, 1], '\n') print(eigenvector[:, 1], '\n')
# For kx=0.005 # For kx=0.005
print('\nFor kx=0.005:\n') print('\nFor kx=0.005:\n')
eigenvalue, eigenvector = np.linalg.eigh(hamiltonian_of_BBH_model(kx=0.005, ky=0)) eigenvalue, eigenvector = np.linalg.eigh(hamiltonian_of_bbh_model(kx=0.005, ky=0))
print(eigenvalue, '\n') print(eigenvalue, '\n')
print(eigenvector[:, 0]) print(eigenvector[:, 0])
print(eigenvector[:, 1], '\n\n') print(eigenvector[:, 1], '\n\n')
@ -75,4 +75,4 @@ print(vector2, '\n')
# # 安装命令：pip install --upgrade guan。 # # 安装命令：pip install --upgrade guan。
# import guan # import guan
# vector1, vector2 = guan.rotation_of_degenerate_vectors(vector1, vector2, index1=None, index2=None, precision=0.01, criterion=0.01, show_theta=0) # vector1, vector2 = guan.rotation_of_degenerate_vectors(vector1, vector2, index1=None, index2=None, precision=0.01, criterion=0.01, show_theta=0)
# hamiltonian = guan.hamiltonian_of_BBH_model(kx, ky, gamma_x=0.5, gamma_y=0.5, lambda_x=1, lambda_y=1) # hamiltonian = guan.hamiltonian_of_bbh_model(kx, ky, gamma_x=0.5, gamma_y=0.5, lambda_x=1, lambda_y=1)