Questions About Magnetic Field Units in Landau_level.py for Landau Level Calculations

Dear Kwant Community, I am writing to seek assistance with two issues I have encountered while using Kwant for solving transport problems. Your guidance would be greatly appreciated. 1. I am trying to calculate the Landau levels of monolayer graphene using the following Hamiltonian (python): ``` hamiltonian = """ + hbar_vF * (k_x * kron(sigma_0, sigma_x) + k_y * kron(sigma_0, sigma_y)) """ params = dict(hbar_vF=6.582) ``` In the documentation for the `discretize_landau` function, it is stated that "The units of magnetic field are :math:`ϕ₀ / 2 π a²` with :math:`ϕ₀ = h/e` the magnetic flux quantum and :math:`a` the unit length." Is a^2 here the area (A) of a primitive cell of graphene (A = 5.24 Ų)? How can I convert this unit into Tesla? 2. Can I use this code to calculate the Landau levels of multilayer stacking graphene with a sequences of hoppings ? Are there specific considerations or modifications I should be aware of when extending the calculation to multilayer graphene? Thank you in advance for your help. Best regards, ac76888

Attatchment: test code for calculating lls of monolayer graphene hamiltonian = """ + hbar_vF * (k_x * kron(sigma_0, sigma_x) + k_y * kron(sigma_0, sigma_y)) """ params = dict(hbar_vF=6.582) syst = kwant.continuum.discretize_landau(hamiltonian, N=20) syst = syst.finalized() b_values = np.linspace(0, 0.00002, 200) fig = kwant.plotter.spectrum(syst, ('B', b_values), params=params, show=False,file = 'example.png') fig.savefig('example.png')
participants (1)
-
ac76888@sjtu.edu.cn