[Pw_forum] Different band gap value from the 'relax' and 'nscf'calculation.

Giovanni Cantele Giovanni.Cantele at na.infn.it
Thu Nov 3 17:43:37 CET 2005


Dear TAO JUNGUANG,
try to run one scf run at the optimized positions and immediatly after
the nscf.
I tried to do that with espresso-2.1.4 and got


  -47.5120 -47.4940 -24.1405 -24.1405 -23.9669 -23.9669 -23.9396 -23.9015
    -9.1803  -7.9537  -7.1715  -7.1715   2.9016   3.3425   3.7211   4.8903
     4.8903   6.1633   6.1633   6.7573   8.0547   8.3102   8.3102   8.7543
    10.5831  10.8792  11.1367  11.2196  11.2196  13.2185  13.3614  15.4703


So, 10.5831-8.7543=1.83 eV

Perhaps the pot and/or tho files were for some reason the wrong ones.

With previous releases I found some problem, due to the fact that
the potential after the nscf calculation seemed to change because 
overwritten.
This should have been corrected  in 2.1.4 or 2.1.5. In the Changelog of 
the cvs
version I find:

2005-05-27 14:52  giannozz

        * PW/: atomic_rho.f90, electrons.f90, potinit.f90, punch.f90:

        - do not set to zero the negative part of the charge density
          when starting from superpositions of atomic charges: it is useless
        - print HOMO/top of the VB and LUMO/bottom of the CB when available
        - always start from charge density and recalculate the potential,
          also in a non-scf calculation
        - do not overwrite the charge density after a non-scf calculation

So your problem could be due to overwritten charge density after
the nscf calculation.

Giovanni

-- 




Dr. Giovanni Cantele
Coherentia CNR-INFM and Dipartimento di Scienze Fisiche
Universita' di Napoli "Federico II"
Complesso Universitario di Monte S. Angelo - Ed. G
Via Cintia, I-80126, Napoli, Italy
Phone: +39 081 676910
Fax:   +39 081 676346
E-mail: Giovanni.Cantele at na.infn.it
Web: http://people.na.infn.it/~cantele

****************
RICERCA E FUTURO CAMMINANO INSIEME:
  DA 41 MESI I LAVORATORI DELLA
  RICERCA HANNO IL CONTRATTO SCADUTO,
  IL FUTURO NON DEVE SCADERE!
  per maggiori info: http://www.ge.cnr.it/rsu/rsu.htm
****************




More information about the Pw_forum mailing list