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