Re: [Kwant] 1D Leads

Hi Thomas,

Thomas Schmidt thomas.schmidt@unibas.ch writes:

I'm trying out Kwant, and I've encountered a problem when investigating 1D leads.

The problem that you encounter can be boiled down to the following program

import kwant, numpy kwant.physics.modes(numpy.array([[2]]), numpy.array([[-1]]))

You are asking for modes exactly at an energy where a new mode opens. This is something that is ill-defined. For now, to avoid this problem, you can simply add a small epsilon to your energies, like this:

energies=[0.01 * i + 1e-6 for i in xrange(1000)]

I'm not sure why the same problem does not appear for wider systems. (I'm not the one who wrote that part of Kwant.)

Cheers, Christoph