From b39bf04dd4b5c74a384877d48c5518c66918c36a Mon Sep 17 00:00:00 2001 From: guanjihuan Date: Tue, 21 Dec 2021 18:04:54 +0800 Subject: [PATCH] update --- .../2D_graphene.py | 6 +++--- .../flat_bands_of_kagome lattice.py | 6 +++--- 2 files changed, 6 insertions(+), 6 deletions(-) diff --git a/academic_codes/2019.10.23_Hamiltonian_and_bands_of_graphene/2D_graphene.py b/academic_codes/2019.10.23_Hamiltonian_and_bands_of_graphene/2D_graphene.py index fa37b5c..a68736a 100755 --- a/academic_codes/2019.10.23_Hamiltonian_and_bands_of_graphene/2D_graphene.py +++ b/academic_codes/2019.10.23_Hamiltonian_and_bands_of_graphene/2D_graphene.py @@ -33,8 +33,8 @@ def hamiltonian(k1, k2, M, t1, a=1/sqrt(3)): # graphene哈密顿量（a为原 def main(): hamiltonian0 = functools.partial(hamiltonian, M=0, t1=1, a=1/sqrt(3)) # 使用偏函数，固定一些参数 - k1 = np.linspace(-2*pi, 2*pi, 800) - k2 = np.linspace(-2*pi, 2*pi, 800) + k1 = np.linspace(-2*pi, 2*pi, 500) + k2 = np.linspace(-2*pi, 2*pi, 500) plot_bands_two_dimension(k1, k2, hamiltonian0) @@ -57,7 +57,7 @@ def plot_bands_two_dimension(k1, k2, hamiltonian): ax = fig.gca(projection='3d') k1, k2 = np.meshgrid(k1, k2) for dim0 in range(dim): - ax.plot_surface(k1, k2, eigenvalue_k[:, :, dim0], cmap=cm.coolwarm, linewidth=0, antialiased=False) + ax.plot_surface(k1, k2, eigenvalue_k[:, :, dim0], rcount=200, ccount=200, cmap=cm.coolwarm, linewidth=0, antialiased=False) plt.show() diff --git a/academic_codes/2021.08.09_flat_bands_of_kagome lattice/flat_bands_of_kagome lattice.py b/academic_codes/2021.08.09_flat_bands_of_kagome lattice/flat_bands_of_kagome lattice.py index 722e8e5..22795be 100644 --- a/academic_codes/2021.08.09_flat_bands_of_kagome lattice/flat_bands_of_kagome lattice.py +++ b/academic_codes/2021.08.09_flat_bands_of_kagome lattice/flat_bands_of_kagome lattice.py @@ -20,7 +20,7 @@ def hamiltonian(kx, ky): # kagome lattice h = -t*h return h -kx_array = np.linspace(-pi ,pi, 100) -ky_array = np.linspace(-pi ,pi, 100) +kx_array = np.linspace(-pi ,pi, 500) +ky_array = np.linspace(-pi ,pi, 500) eigenvalue_array = guan.calculate_eigenvalue_with_two_parameters(kx_array, ky_array, hamiltonian) -guan.plot_3d_surface(kx_array, ky_array, eigenvalue_array, xlabel='kx', ylabel='ky', zlabel='E') \ No newline at end of file +guan.plot_3d_surface(kx_array, ky_array, eigenvalue_array, xlabel='kx', ylabel='ky', zlabel='E', rcount=200, ccount=200) \ No newline at end of file