Hi Harshad, > 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?

Sorry that I forgot to clarify that the unit of length in this case is the side of the lead unit cell. Multiplying by that size should allow to easily translate your velocity into m/s. Anton