Hi,

You could do that, but this is already done in Kwant, and you can read the velocities off from the modes object [1].

I was able to obtain the velocities from the modes object. I confirmed that modes.velocities is just the inverse of integral |𝜓|2 over the unit cell. If I did the math correctly, I got the units of modes.velocities to be eV (energy units of my Hamiltonian). How can I convert this to lets say m/s?

Thanks,

Harshad

On Mon, Jan 23, 2017 at 10:18 AM, Harshad Sahasrabudhe <hsahasra@purdue.edu> wrote:

Great, thanks a lot!

On Mon, Jan 23, 2017 at 10:15 AM, Anton Akhmerov <anton.akhmerov+kd@gmail.com> wrote:
> Thank you. Then to calculate the velocity, should I just divide the
> probability current by the integral of |𝜓|2 over the unit cell?

You could do that, but this is already done in Kwant, and you can read
the velocities off from the modes object [1].

Best,
Anton

[1]: https://kwant-project.org/doc/1/reference/generated/kwant.physics.PropagatingModes#kwant.physics.PropagatingModes