[Pw_forum] Having problem with Total energy in Band structure calculation

Ari P Seitsonen Ari.P.Seitsonen at iki.fi
Sun Jan 12 07:58:19 CET 2014


Dear Alhassan Shuaibu,

   As already stated in the previous answer, there is a problem in your 
structure: I guess that the coordinates are _not_ given Ångström like 
stated in your input. And the k points, ibrav = 5 would suggest that you 
need a more homogeneous k point sampling in all directions, doesn't it?

   Indeed the cut-off energz is far too low; I don't know these 
elements/pseudo potentials too well, but they might also require an 
ecutwfc of > 100 Ry. Please test carefully (eg forces as a function of the 
cut-off energy etc).

     Greetings,

        apsi

-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-
   Ari P Seitsonen / Ari.P.Seitsonen at iki.fi / http://www.iki.fi/~apsi/
   Institut für Chemie der Universität Zürich, CH-8057 Zürich
   Tel: +41 44 63 54 497  /  Mobile: +41 79 71 90 935


On Sun, 12 Jan 2014, Alhassan Shuaibu wrote:

> Hello
> I am trying to do band structure calculations for Ternary Bi2Te2S.But I am getting a positive TOTAL energy . I have checked my crystal structure with Xcrysden, looks
> okay to me. please is this related to my generated coordinate?
> my INPUT FILE
> &CONTROL
>     calculation='scf'
>     restart_mode='from_scratch',
>     prefix='Bi2Te2SRH_exc1',
>     tstress = .true.
>     tprnfor = .true.
>     pseudo_dir ='/home/alhas/Alhas/espresso-5.0.2/pseudo/' ,
>     outdir='/home/alhas/Alhas/espresso-5.0.2/bin/' ,
>  /
>  &SYSTEM
>                     ibrav =  5,
>                      a = 10.931,
>                      cosab = 0.922901387,
>                      nat = 5,
>                      ntyp = 3,
>                      ecutwfc = 20 ,
>                     
> /
>  &ELECTRONS
>                conv_thr = 1.0d-8
>                 startingwfc = 'atomic' ,
>                  mixing_mode = 'plain' ,
>                  mixing_beta = 0.7 ,
>                 
> /
> 
> ATOMIC_SPECIES
> Bi 208.9804 Bi.pbe-hgh.UPF
> Te 127.6    Te.pbe-hgh.UPF
> S  32.065   S.pbe-hgh.UPF
> ATOMIC_POSITIONS (angstrom)
>   Bi      0.38788600000000   0.38788600000000   0.38788600000000  
>   Bi      0.61211400000000   0.61211400000000   0.61211400000000 
>   Te      0.77733500000000   0.77733500000000   0.77733500000000 
>   Te      0.22266500000000   0.22266500000000   0.22266500000000 
>   S       0.00000000000000   0.00000000000000   0.00000000000000 
> K_POINTS (automatic)
>  7 7 1 0 0 0   
> and the OUT PUT FILE LOOK LIKE
>


More information about the users mailing list