[Pw_forum] convergence not achieved .related problem

m rostami rostami7.m at gmail.com
Fri Apr 27 16:39:14 CEST 2012


Dear users
I have been trying to relax a graphene structure with a vacancy but
many erros stoped optimization and convergence not achieved.like

1. Efermy cannot braket :
2. to0 many bands not coveraged

for solving the errors I have chenged these parametrs

 increasing ecutoff to 60 and ecuthro to 720 .increasing edeguss to
0.07 and changing diagnolization from cg to david  . I have tested
smearing for m-v and mp .increasing nbands from   70 to 120  too .but
optimization didnt achived .I appreciate if any one helpe me .

my input file is :
&control
    calculation='relax',
    prefix='graphene'
    pseudo_dir = '/media/espresso-4.2/pseudo/',
    outdir='./tmp/',
    etot_conv_thr=1.0D-3 ,
    forc_conv_thr=1.0D-2 ,

 /
 &system
    ibrav=  4, celldm(1) =13.9506, celldm(3)=1.667,nat=35, ntyp= 1,nbnd=120
   ecutwfc =80.0,ecutrho = 900,
    occupations='smearing',smearing='m-v',degauss=0.07,
 /
 &electrons
  diagonalization='cg',
  conv_thr=1.D-8 ,
  mixing_beta=0.1,
 /
 &ions
 upscale=20,
 /
ATOMIC_SPECIES
C  1.00  C.blyp-mt.UPF
ATOMIC_POSITIONS {angstrom}
    C    0.0071683270    0.0041386380    1.6740000000
 .........
......
K_POINTS  AUTOMATIC
8 8 1 0 0 0


.............................
I am not sure about these questions to0 :.

1.in example 3 the atomic weight for different atoms is set to 1 in
relax process.is this correct? as I did for carbon (ATOMIC_SPECIES
C  1.00  C.blyp-mt.UPF)
I am not sure about k-point .is 8 8 1 0 0 0  enough for this structure
with 10 A` vaccum .

how many nbands I should for coverage ?

wich diagnbolization is better david or cg ?

mixing_beta should be around 0.1-0.3?


thanks alot

Melika Rostami
Master Student of Sharif University


More information about the Pw_forum mailing list