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 <hsaha...@purdue.edu> wrote: > Great, thanks a lot! > > On Mon, Jan 23, 2017 at 10:15 AM, Anton Akhmerov < > anton.akhmerov...@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.ph >> ysics.PropagatingModes#kwant.physics.PropagatingModes >> > >