Problem in the structure optimization

[lxbyf]

Dear all,

I am trying to optimize the structure, but the calculation will automatic stop with in a step of self-consistent-field convergence and the abinit process will vanish. Moreover, the calculation will stop in different SCF step, if I use different pseudopotentials files.

There are some warnings in the log files, in such a way as:

scfcge: WARNING -
Potential-based CG line minimization has trouble to converge.
The algorithm is restarted with more secure parameters.
scfcge: actual 45-20-1 2.8830E-03 -1.454028524146E+02 1.1153E-03
scfcge: predict -1.1231E-01 restart the algorithm
scfcge:
scfcge: start 45-21-0 0.0000E+00 -1.454028524146E+02 1.1153E-03

and at last, it became:

scfcge : ERROR -
Potential-based CG line minimization not converged after 13 restarts.
Action : read the eventual warnings about lack of convergence.
Some might be relevant. Otherwise, raise nband.
-P-0000
-P-0000 leave_new : decision taken to exit ...
-P-0000 leave_new : synchronization done...
-P-0000 leave_new : exiting...

My structure is 2*2*2 AlGaN with 32 atoms, and the input files is as follows. I guess if I should set the variable "nband" and how to do it?

ndtset 2
chkprim 0

ionmov1 3
optcell1 0
ntime1 100
tolmxf1 5.0d-5

dilatmx2 1.5
getxred2 -1
getwfk2 -1
ionmov2 3
ntime2 100
optcell2 2
tolmxf2 5.0d-5

acell 3*3.189 angstrom
rprim -0.0000000000000000 -0.9898826185140575 -0.0000000000000000
0.8572634936377291 0.4949413092570287 -0.0000000000000000
0.0000000000000000 0.0000000000000000 3.2500000000000000
ntypat 3
znucl 13 31 7
natom 32
typat 4*1 12*2 16*3
ecut 16.0
ecutsm 0.5
kptopt 1
ngkpt 4 4 4
nshiftk 1
shiftk 0.0 0.0 0.5
diemac 9.0
iscf 5
nstep 50
toldff 5.0d-6
xred 0.666666667 0.333333333 0
0.333333333 0.666666667 0.25
0.666666667 0.333333333 0.5
0.333333333 0.666666667 0.75
0.166666667 0.333333333 0
0.166666667 0.833333333 0
0.666666667 0.833333333 0
0.333333333 0.166666667 0.25
0.833333333 0.166666667 0.25
0.833333333 0.666666667 0.25
0.166666667 0.333333333 0.5
0.166666667 0.833333333 0.5
0.666666667 0.833333333 0.5
0.333333333 0.166666667 0.75
0.833333333 0.166666667 0.75
0.833333333 0.666666667 0.75
0.166666667 0.333333333 0.1875
0.666666667 0.333333333 0.1875
0.166666667 0.833333333 0.1875
0.666666667 0.833333333 0.1875
0.333333333 0.166666667 0.4375
0.833333333 0.166666667 0.4375
0.333333333 0.666666667 0.4375
0.833333333 0.666666667 0.4375
0.166666667 0.333333333 0.6875
0.666666667 0.333333333 0.6875
0.166666667 0.833333333 0.6875
0.666666667 0.833333333 0.6875
0.333333333 0.166666667 0.9375
0.833333333 0.166666667 0.9375
0.333333333 0.666666667 0.9375
0.833333333 0.666666667 0.9375

[ilukacevic]

Dear lxbyf,

try with iscf 7. If it still doesn't work, increase (keeping the iscf 7) npulayit, nline and/or nnsclo.

Cheers!

Igor Lukacevic