Dear Anton and other Kwant users, I am currently using the KPM method to calculate the density of states: spectrum = kwant.kpm.SpectralDensity(fsyst) dos = spectrum(energies) In the tutorial, this method is successfully applied to calculate the dos of a large closed system. But I am wondering that here it can also be applied to the case of an open system (i.e., a system composed of scattering region and leads). Thanks a lot ! Best regards, Kuangyia
Dear Kuangyia,
KPM method can be efficiently applied to any tight-binding system,
including infinite ones. Unfortunately, I didn't encounter an
implementation of KPM suitable for these geometries (e.g. it has to
deal with variable-sized vectors). In addition to that, the Kwant's
implementation of KPM is relatively basic, it does not implement
possible optimizations and it does not implement several useful
generalizations of KPM.
Best,
Anton
On Wed, Oct 4, 2017 at 10:53 AM, kuangyia lee
Dear Anton and other Kwant users,
I am currently using the KPM method to calculate the density of states:
spectrum = kwant.kpm.SpectralDensity(fsyst) dos = spectrum(energies)
In the tutorial, this method is successfully applied to calculate the dos of a large closed system. But I am wondering that here it can also be applied to the case of an open system (i.e., a system composed of scattering region and leads).
Thanks a lot !
Best regards, Kuangyia
Dear Anton,
Thank you for your reply. I would make a test to see what's going on.
Best,
Kuangyia
On Sun, Oct 8, 2017 at 4:07 PM, Anton Akhmerov
Dear Kuangyia,
KPM method can be efficiently applied to any tight-binding system, including infinite ones. Unfortunately, I didn't encounter an implementation of KPM suitable for these geometries (e.g. it has to deal with variable-sized vectors). In addition to that, the Kwant's implementation of KPM is relatively basic, it does not implement possible optimizations and it does not implement several useful generalizations of KPM.
Best, Anton
On Wed, Oct 4, 2017 at 10:53 AM, kuangyia lee
wrote: Dear Anton and other Kwant users,
I am currently using the KPM method to calculate the density of states:
spectrum = kwant.kpm.SpectralDensity(fsyst) dos = spectrum(energies)
In the tutorial, this method is successfully applied to calculate the dos of a large closed system. But I am wondering that here it can also be applied to the case of an open system (i.e., a system composed of scattering region and leads).
Thanks a lot !
Best regards, Kuangyia
participants (2)
-
Anton Akhmerov
-
kuangyia lee