Ynt: Re: Datta’s ballistic transport formalism vs KWANT
Thanks Joe! That was indeed the case! It became much closer with small a. Best, Ran � Hi Ran, > I tried to compare KWANT’s results for transmission with Datta’s ballistic transport formalism where total transmission is written as > > Ttot=T(E)M(E) > > Here Datta takes T(E)=1 for ballistic transport (please see: J. Appl. Phys. 105, 034506, 2009) and M(E) is the number of modes in transverse direction. When I compared KWANT's results with Datta’s expression, for the system given in “quantum_wire_revisited.py”, I found different results (please see the attached figure where I tried to put every relevant thing in the calculation). Since the reflectance is zero for that system and so transmission is 1 for each mode, shouldn’t it give the same results with Datta’s transmission expression? > Nice question! Looking at your results it seems that the energies at which new modes open is shifted with respect to Datta's result. I believe that this is simply due to the fact that your discretization is not fine enough. Datta's result is valid in the continuum limit, whereas the Kwant simulation (in the case presented) uses a finitedifference discretization to render the problem discrete. If you decrease the 'a' parameter, you should see the discrepancy between the two result decrease. Happy Kwanting, Joe
Dear Mutcran, I would like just to make an additional remark. Indeed, the formula in Datta's book is for a continuous model where the maximum number of modes is given by M=2 W/ \lambda_f (\lambda_f is the Fermi wavelength). In the case of a discrete model, like in Kwant, you can deduce M from the band structure. In fact, we have: E=4t 2t cos(kx)2t cos(m pi/(N+1)), where N is the number of sites for one transversal cell. The maximum number of transmitted modes is obtained when the longitudinal kinetic energy is 0 (kx~0). so we get M= (N+1)/pi arcos((E2t)/(2t)) (you need to take the integer part of this expression). This result is exact for a discrete square lattice model. A simple Taylor expansion around E~0 will give you the Datta's formula (where W=(N+1)a). As a conclusion, the number of transmitted modes for a given discrete model is obtained from the form of the Band structure. This allows you to calculate it for any complex lattice (as long as you can get analytically the band structure). I hope this helps, Adel On Tue, Jan 22, 2019 at 11:22 PM mutcran <mutcran@mynet.com> wrote:
Thanks Joe! That was indeed the case! It became much closer with small a.
Best,
Ran
Hi Ran,
I tried to compare KWANT’s results for transmission with Datta’s ballistic transport formalism where total transmission is written as
Ttot=T(E)M(E)
Here Datta takes T(E)=1 for ballistic transport (please see: J. Appl. Phys. 105, 034506, 2009) and M(E) is the number of modes in transverse direction. When I compared KWANT's results with Datta’s expression, for the system given in “quantum_wire_revisited.py”, I found different results (please see the attached figure where I tried to put every relevant thing in the calculation). Since the reflectance is zero for that system and so transmission is 1 for each mode, shouldn’t it give the same results with Datta’s transmission expression?
Nice question!
Looking at your results it seems that the energies at which new modes open is shifted with respect to Datta's result.
I believe that this is simply due to the fact that your discretization is not fine enough. Datta's result is valid in the continuum limit, whereas the Kwant simulation (in the case presented) uses a finitedifference discretization to render the problem discrete. If you decrease the 'a' parameter, you should see the discrepancy between the two result decrease.
Happy Kwanting,
Joe
 Abbout Adel
participants (2)

Abbout Adel

mutcran