Density of states for a lead

Dear all, is it implemented a way to calculate the density or the density of states of a lead (defined as kwant.TranslationalSymmetry ...)? Thank you. Regards, Jacopo

Dear Jacopo, There's nothing out of the box, however the calculation is rather straightforward. If you want to compute the LDOS of a lead at a given energy, compute its modes (via lead.modes), then take the propagating ones and square their sum. Since the lead modes are normalized to carry unit current, this procedure is sufficient. If you only wan the overall DOS, sum the inverse velocities, as returned by the modes. Come to think of it, it would make a lot of sense to implement density of states as a convenience method of propagating modes. I've opened https://gitlab.kwant-project.org/kwant/kwant/-/issues/376 ; anyone interested? Best, Anton On Thu, 16 Apr 2020 at 22:07, Jacopo Settino <jac.settino@gmail.com> wrote:
Dear all, is it implemented a way to calculate the density or the density of states of a lead (defined as kwant.TranslationalSymmetry ...)? Thank you. Regards, Jacopo

Thank you very much. I didn't think about the possibility of having velocities already calculated. Regards, Jacopo Il giorno sab 18 apr 2020 alle ore 12:01 Anton Akhmerov < anton.akhmerov+kd@gmail.com> ha scritto:
Dear Jacopo,
There's nothing out of the box, however the calculation is rather straightforward. If you want to compute the LDOS of a lead at a given energy, compute its modes (via lead.modes), then take the propagating ones and square their sum. Since the lead modes are normalized to carry unit current, this procedure is sufficient. If you only wan the overall DOS, sum the inverse velocities, as returned by the modes. Come to think of it, it would make a lot of sense to implement density of states as a convenience method of propagating modes. I've opened https://gitlab.kwant-project.org/kwant/kwant/-/issues/376 ; anyone interested?
Best, Anton
On Thu, 16 Apr 2020 at 22:07, Jacopo Settino <jac.settino@gmail.com> wrote:
Dear all, is it implemented a way to calculate the density or the density of
states of a lead (defined as kwant.TranslationalSymmetry ...)?
Thank you. Regards, Jacopo
participants (2)
-
Anton Akhmerov
-
Jacopo Settino