This commit is contained in:
guanjihuan 2022-04-01 16:40:05 +08:00
parent c80a3e2c1e
commit f77d69dfb8
4 changed files with 5 additions and 5 deletions

View File

@ -4,4 +4,4 @@ import numpy as np
k_array = np.linspace(-np.pi, np.pi, 100) k_array = np.linspace(-np.pi, np.pi, 100)
hamiltonian_function = guan.one_dimensional_fourier_transform_with_k(unit_cell=0, hopping=1) # one dimensional chain hamiltonian_function = guan.one_dimensional_fourier_transform_with_k(unit_cell=0, hopping=1) # one dimensional chain
eigenvalue_array = guan.calculate_eigenvalue_with_one_parameter(k_array, hamiltonian_function) eigenvalue_array = guan.calculate_eigenvalue_with_one_parameter(k_array, hamiltonian_function)
guan.plot(k_array, eigenvalue_array, xlabel='k', ylabel='E', type='-k') guan.plot(k_array, eigenvalue_array, xlabel='k', ylabel='E', style='-k')

View File

@ -3,6 +3,6 @@ import numpy as np
k_array = np.linspace(-np.pi, np.pi, 100) k_array = np.linspace(-np.pi, np.pi, 100)
eigenvalue_array = guan.calculate_eigenvalue_with_one_parameter(k_array, guan.hamiltonian_of_square_lattice_in_quasi_one_dimension) eigenvalue_array = guan.calculate_eigenvalue_with_one_parameter(k_array, guan.hamiltonian_of_square_lattice_in_quasi_one_dimension)
guan.plot(k_array, eigenvalue_array, xlabel='k', ylabel='E', type='-k') guan.plot(k_array, eigenvalue_array, xlabel='k', ylabel='E', style='-k')
eigenvalue_array = guan.calculate_eigenvalue_with_one_parameter(k_array, guan.hamiltonian_of_graphene_with_zigzag_in_quasi_one_dimension) eigenvalue_array = guan.calculate_eigenvalue_with_one_parameter(k_array, guan.hamiltonian_of_graphene_with_zigzag_in_quasi_one_dimension)
guan.plot(k_array, eigenvalue_array, xlabel='k', ylabel='E', type='-k') guan.plot(k_array, eigenvalue_array, xlabel='k', ylabel='E', style='-k')

View File

@ -4,7 +4,7 @@ import numpy as np
hamiltonian = guan.hamiltonian_of_finite_size_system_along_two_directions_for_square_lattice(2,2) hamiltonian = guan.hamiltonian_of_finite_size_system_along_two_directions_for_square_lattice(2,2)
fermi_energy_array = np.linspace(-4, 4, 400) fermi_energy_array = np.linspace(-4, 4, 400)
total_dos_array = guan.total_density_of_states_with_fermi_energy_array(fermi_energy_array, hamiltonian, broadening=0.1) total_dos_array = guan.total_density_of_states_with_fermi_energy_array(fermi_energy_array, hamiltonian, broadening=0.1)
guan.plot(fermi_energy_array, total_dos_array, xlabel='E', ylabel='Total DOS', type='-') guan.plot(fermi_energy_array, total_dos_array, xlabel='E', ylabel='Total DOS', style='-')
fermi_energy = 0 fermi_energy = 0
N1 = 3 N1 = 3

View File

@ -5,7 +5,7 @@ fermi_energy_array = np.linspace(-4, 4, 400)
h00 = guan.hamiltonian_of_finite_size_system_along_one_direction(4) h00 = guan.hamiltonian_of_finite_size_system_along_one_direction(4)
h01 = np.identity(4) h01 = np.identity(4)
conductance_array = guan.calculate_conductance_with_fermi_energy_array(fermi_energy_array, h00, h01) conductance_array = guan.calculate_conductance_with_fermi_energy_array(fermi_energy_array, h00, h01)
guan.plot(fermi_energy_array, conductance_array, xlabel='E', ylabel='Conductance', type='-') guan.plot(fermi_energy_array, conductance_array, xlabel='E', ylabel='Conductance', style='-')
fermi_energy = 0 fermi_energy = 0
guan.print_or_write_scattering_matrix(fermi_energy, h00, h01) guan.print_or_write_scattering_matrix(fermi_energy, h00, h01)