Dear all:

I want to know some details about how SIESTA computes LDOS, because I found
something weired in my results.

Since in the input file, only energy range is required to control the LDOS,
my intuitive thinking is that SIESTA first chooses eigenstates within the
user-specified energy range from the k mesh and then adds the norm square of
all the wavefunctions together to get the LDOS.

However, I saw some discussions on the mailing-list saying broadening
has been made in LDOS calculation. Also, I found when the energy range is
extremely narrow, the LDOS plot doesn't approach the wavefunction plot.

So I'm confused.  Would you please give me some ideas or references about
how LDOS is obtained in SIESTA?

I'll really appreciate that.

Zheng Liu
Center for Advanced Study
Tsinghua University
Beijing, China, 100084

Reply via email to