>
> Yes, that is correct. The equation that you quoted is indeed the probability
> current
Thank you. Then to calculate the velocity, should I just divide the
probability current by the integral of |휓|2 over the unit cell?
On Mon, Jan 23, 2017 at 10:00 AM, Anton Akhmerov <
rov...@gmail.com
> wrote:
> Dear Hardshad,
>
> I might be missing some factors here, but the unit of current is
> elementary charge / unit of time. The unit of time is hbar / unit of energy
> that you used in defining your tight-binding model.
>
> Best,
> Anton
>
> On M
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
Hi,
I was wondering what units would the wavefunction obtained from Kwant
have? I was thinking they would have the units 1/sqrt(nm.eV) (since my
energies are in eV and lengths are in nm) as the modes are normalized
according to
[image: Inline image 1]
How would I calculate the velocity of the
ath 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
Can I just do 1/2 m*v^2 = modes.velocities, assuming parabolic dispersion?
On Mon, Jan 23, 2017 at 3:13 PM, Harshad Sahasrabudhe <hsaha...@purdue.edu>
wrote:
> Hi,
>
> You could do that, but this is already done in Kwant, and you can read the
>> velocities off from the mode
Hi All,
I am trying to calculate current density from the wavefunctions using the
following code:
for i in range(Np):
for j in range(Nc):
lat_idx_i = i-floor(Np/2)
lat_idx_j = j-floor(Nc/2)
site_i = lat(lat_idx_i, lat_idx_j)
idx_i = sys.sites.index(site_i)
color is for the opisit directions
> def Bond_current(site1,site2,phi=phi):
> i,j = sites.index(site1),sites.index(site2)
> return 6*abs(I[i,j])
>
> def Current_color(site1,site2,phi=phi):
> i,j = sites.index(site1),sites.index(site2)
> if (site1.pos[0]>site2.