force cannot converge

guitar · Post by **guitar** » Fri Mar 25, 2011 3:34 am

Dear All:

I am doing a calculation about slab relaxation. But the force cannot converge to target. What the matter with it?
My input file is as follows:

Code: Select all

kptopt 1
prtwf 0
istwfk *1
npkpt 8
npband 1
npfft 1
bandpp 1
paral_kgb 1
ngkpt 10 16 2
chkprim 0
chksymbreak 0
toldfe 1.0d-6
occopt 3
iscf 5
nstep 60
tsmear 0.01
dilatmx 1.10
ecut 30
enunit 2
typat 1 1 1 1 1 1 1 1 1 1 1 1  
znucl 41
natom 12
acell 6.316 6.316 6.316
optcell 0
ecutsm 0.5
ionmov 3
ntime 60
tolmxf 2.5d-4
natfixx 12
natfixy 12
iatfixx 1 2 3 4 5 6 7 8 9 10 11 12
iatfixy 1 2 3 4 5 6 7 8 9 10 11 12
rprim 1.00 -1.00 0.00
      0.50 0.50 0.50
      -2.00 -2.00 4.00
xred  0.00 0.00 0.00
      0.00 0.10 0.50
      0.50 0.00 0.25
      0.50 0.10 0.75
      0.00 0.6666666666 0.1666666666
      0.00 0.7666666666 0.6666666666
      0.50 0.6666666666 0.4166666666
      0.50 0.7666666666 0.9166666666
      0.50 0.3333333333 0.0833333333
      0.50 0.4333333333 0.5833333333
      0.00 0.3333333333 0.3333333333
      0.00 0.4333333333 0.8333333333

I really appreciate if someone can help me.

Ken