Dear Pablo,

 in the example here:

https://kwant-project.org/doc/dev/tutorial/kpm

you write:

..

fsyst_staggered  =  make_syst_staggered().finalized()
# find 'A' and 'B' sites in the unit cell at the center of the disk
center_tag  =  np.array([0,  0])
where  =  lambda  s:  s.tag  ==  center_tag
# make local vectors
vector_factory  =  kwant.kpm.LocalVectors(fsyst_staggered,  where)

and plot their respective local density of states.

# 'num_vectors' can be unspecified when using 'LocalVectors'
local_dos  =  kwant.kpm.SpectralDensity(fsyst_staggered,  num_vectors=None,
                                      vector_factory=vector_factory,
                                      mean=False,
                                      rng=0)
energies,  densities  =  local_dos()

.

My question is:  how do I know which sublattice LDOS comes out first in the result?

I plan to compute LDOS in a lot of point (say, 50 points) at a given set of energies.  So, what will be the ordering of vectors given by vector_factory.   Or, can I create those vectors manually without the vector factory?

Thank you,

Sergey


Reply via email to