Re: [Kwant] Distribution of spin in graphene lattice
Dear Anton,
Thanks for your reply. If I don't include the spin-orbital coupling term
and just want to see how are the spin up /down components of the
wavefunction distributed on the graphene surface, is it possible to do it by
1. change the n.n. hopping from -t to -t*sigma0
2. just plot out the even/odd index of the wavefunction (just like in
https://mailman-mail5.webfaction.com/pipermail/kwant-discuss/2013-September/...
)?
-Originally, I have 86 wavefunctions in total (since there are 86
eigenvalues).
-By using the above approach, I get 172 wavefunctions (now there are 172
eigenvalues, ).
-But when I sum the ground state wavefunction together, I cannot recover
the original ground state wavefunction.
Part of my code is as follows.
####### Case 1: hopping = -t * sigma_0 ########
iplot = 1
evecs = la.eigh(ham)[1]
wf_up = abs(evecs[::2, 2*iplot])**2
wf_down = abs(evecs[1::2,2*iplot])**2
Note: Now the number of wavefunction obtained is doubled (say 2n in total).
But the 1st and the 2nd one seems to be the same, so I just pick the
(2*i)th wavefunction
wf = wf_up + wf_down
def site_size_sum(i):
return 1* wf[i] / wf.max()
def site_size_up(i):
return 1* wf_up[i] / wf.max()
def site_size_down(i):
return 1* wf_down[i] / wf.max()
kwant.plot(sys, site_size=site_size_up, site_color=(0, 0, 1,
0.3),hop_lw=0.1)
kwant.plot(sys, site_size=site_size_down, site_color=(0, 0, 1,
0.3),hop_lw=0.1)
kwant.plot(sys, site_size=site_size_sum, site_color=(0, 0, 1,
0.3),hop_lw=0.1)
####### Case 2: hopping = -t ##########
wf = abs(evecs[:, iplot])**2
def site_size(i):
return 1* wf[i] / wf.max()
kwant.plot(sys, site_size=site_size, site_color=(0, 0, 1, 0.3),hop_lw=0.1)
I noticed that the plot of site_size_sum in Case 1 is not the same as Case
2.
What should be correct way to do it? Any hint for the mistakes is
appreciated. Thanks.
Best,
Johnny
On 26 November 2015 at 19:01, Anton Akhmerov
Dear Johnny,
Sorry for a delayed reply. Please send the email to the mailing list and not only to me personally (that way someone else may reply if I'm slow, and others can benefit from the discussion).
Best, Anton
On Tue, Nov 24, 2015 at 10:54 AM, T.C. Wu
wrote: Dear Anton,
Thanks for your reply. If I don't include the spin-orbital coupling term and just want to see how are the spin up /down components of the wavefunction distributed on the graphene surface, is it possible to do it by 1. change the n.n. hopping from -t to -t*sigma0 2. just plot out the even/odd index of the wavefunction (just like in
https://mailman-mail5.webfaction.com/pipermail/kwant-discuss/2013-September/... )?
-Originally, I have 86 wavefunctions in total (since there are 86 eigenvalues). -By using the above approach, I get 172 wavefunctions (now there are 172 eigenvalues, ). -But when I sum the ground state wavefunction together, I cannot recover
the
original ground state wavefunction.
What should be correct way to do it? Any hint for the mistakes is appreciated. Thanks.
Best, Johnny
On 22 November 2015 at 20:18, Anton Akhmerov
wrote: Dear Johnny,
Yes, you can superimpose more plots on top of each other, by providing axes returned by one plot as an argument to the next one. If you want to do something with spin you probably should include spin-orbit coupling, yes.
Best, Anton
On Sun, Nov 22, 2015 at 10:22 AM, T.C. Wu
wrote: Dear Kwant developers and users,
For the plot of the graphene wavefunction as illustrated in http://kwant-project.org/doc/1.0/tutorial/tutorial6, is it possible
to
show the distribution of the spin in the same plot? Should I include the spin orbital coupling? Thanks.
Best, Johnny
participants (2)
-
Anton Akhmerov
-
T.C. Wu