[Pw_forum] Too many bands are not converged from nscf calculation
Zhiting Tian
zhiting at MIT.EDU
Sun Sep 25 18:22:44 CEST 2011
Dear all,
I first did scf calculation and it finished. Then I run nscf, but always got error message: "Too many bands are not converged" and then it stopped.
I had read a few posts before and tried to reduce mixing-beta from 0.7 to 0.4 and even to 0.05, but same error. Also, I have tried to change the diagolization='david' to 'cg', same error.
My input files are following:
scf_in:
&control
calculation = 'scf'
restart_mode='from_scratch'
prefix='PbSeAl'
tstress = .true.
tprnfor = .true.
pseudo_dir = "..."
outdir = "..."
/
&system
ibrav= 1, celldm(1)= 22.52, nat= 64, ntyp= 3 ,
ecutwfc = 60,
noncolin = .TRUE.
lspinorb = .TRUE.
occupations = 'smearing'
smearing = 'mv'
degauss = 0.01
/
&electrons
mixing_beta = 0.7D0
conv_thr = 1.0d-10
/
ATOMIC_SPECIES
Al 26.98 Al.rel-pz-nc.UPF
Pb 207.21 Pb.rel-pz-nc.UPF
Se 78.96 Se.rel-pz-nc.UPF
ATOMIC_POSITIONS
Al 0.000000 0.000000 0.000000
Pb 0.500000 0.000000 0.000000
Pb 0.250000 0.250000 0.000000
Pb 0.750000 0.250000 0.000000
Pb 0.000000 0.500000 0.000000
Pb 0.500000 0.500000 0.000000
Pb 0.250000 0.750000 0.000000
Pb 0.750000 0.750000 0.000000
Pb 0.250000 0.000000 0.250000
Pb 0.750000 0.000000 0.250000
Pb 0.000000 0.250000 0.250000
Pb 0.500000 0.250000 0.250000
Pb 0.250000 0.500000 0.250000
Pb 0.750000 0.500000 0.250000
Pb 0.000000 0.750000 0.250000
Pb 0.500000 0.750000 0.250000
Pb 0.000000 0.000000 0.500000
Pb 0.500000 0.000000 0.500000
Pb 0.250000 0.250000 0.500000
Pb 0.750000 0.250000 0.500000
Pb 0.000000 0.500000 0.500000
Pb 0.500000 0.500000 0.500000
Pb 0.250000 0.750000 0.500000
Pb 0.750000 0.750000 0.500000
Pb 0.250000 0.000000 0.750000
Pb 0.750000 0.000000 0.750000
Pb 0.000000 0.250000 0.750000
Pb 0.500000 0.250000 0.750000
Pb 0.250000 0.500000 0.750000
Pb 0.750000 0.500000 0.750000
Pb 0.000000 0.750000 0.750000
Pb 0.500000 0.750000 0.750000
Se 0.250000 0.000000 0.000000
Se 0.750000 0.000000 0.000000
Se 0.000000 0.250000 0.000000
Se 0.500000 0.250000 0.000000
Se 0.250000 0.500000 0.000000
Se 0.750000 0.500000 0.000000
Se 0.000000 0.750000 0.000000
Se 0.500000 0.750000 0.000000
Se 0.000000 0.000000 0.250000
Se 0.500000 0.000000 0.250000
Se 0.250000 0.250000 0.250000
Se 0.750000 0.250000 0.250000
Se 0.000000 0.500000 0.250000
Se 0.500000 0.500000 0.250000
Se 0.250000 0.750000 0.250000
Se 0.750000 0.750000 0.250000
Se 0.250000 0.000000 0.500000
Se 0.750000 0.000000 0.500000
Se 0.000000 0.250000 0.500000
Se 0.500000 0.250000 0.500000
Se 0.250000 0.500000 0.500000
Se 0.750000 0.500000 0.500000
Se 0.000000 0.750000 0.500000
Se 0.500000 0.750000 0.500000
Se 0.000000 0.000000 0.750000
Se 0.500000 0.000000 0.750000
Se 0.250000 0.250000 0.750000
Se 0.750000 0.250000 0.750000
Se 0.000000 0.500000 0.750000
Se 0.500000 0.500000 0.750000
Se 0.250000 0.750000 0.750000
Se 0.750000 0.750000 0.750000
K_POINTS (automatic)
4 4 4 0 0 0
nscf_in:
&control
calculation = 'nscf'
prefix='PbSeAl'
tstress = .true.
tprnfor = .true.
pseudo_dir = "..."
outdir = "..."
/
&system
ibrav= 1, celldm(1)= 22.52, nat= 64, ntyp= 3, nbnd=18
ecutwfc = 60,
noncolin = .TRUE.
lspinorb = .TRUE.
occupations = 'tetrahedra'
/
&electrons
mixing_beta = 0.05D0
conv_thr = 1.0d-10
diagonalization = 'david'
/
ATOMIC_SPECIES
Al 26.98 Al.rel-pz-nc.UPF
Pb 207.21 Pb.rel-pz-nc.UPF
Se 78.96 Se.rel-pz-nc.UPF
ATOMIC_POSITIONS
Al 0.000000 0.000000 0.000000
Pb 0.500000 0.000000 0.000000
Pb 0.250000 0.250000 0.000000
Pb 0.750000 0.250000 0.000000
Pb 0.000000 0.500000 0.000000
Pb 0.500000 0.500000 0.000000
Pb 0.250000 0.750000 0.000000
Pb 0.750000 0.750000 0.000000
Pb 0.250000 0.000000 0.250000
Pb 0.750000 0.000000 0.250000
Pb 0.000000 0.250000 0.250000
Pb 0.500000 0.250000 0.250000
Pb 0.250000 0.500000 0.250000
Pb 0.750000 0.500000 0.250000
Pb 0.000000 0.750000 0.250000
Pb 0.500000 0.750000 0.250000
Pb 0.000000 0.000000 0.500000
Pb 0.500000 0.000000 0.500000
Pb 0.250000 0.250000 0.500000
Pb 0.750000 0.250000 0.500000
Pb 0.000000 0.500000 0.500000
Pb 0.500000 0.500000 0.500000
Pb 0.250000 0.750000 0.500000
Pb 0.750000 0.750000 0.500000
Pb 0.250000 0.000000 0.750000
Pb 0.750000 0.000000 0.750000
Pb 0.000000 0.250000 0.750000
Pb 0.500000 0.250000 0.750000
Pb 0.250000 0.500000 0.750000
Pb 0.750000 0.500000 0.750000
Pb 0.000000 0.750000 0.750000
Pb 0.500000 0.750000 0.750000
Se 0.250000 0.000000 0.000000
Se 0.750000 0.000000 0.000000
Se 0.000000 0.250000 0.000000
Se 0.500000 0.250000 0.000000
Se 0.250000 0.500000 0.000000
Se 0.750000 0.500000 0.000000
Se 0.000000 0.750000 0.000000
Se 0.500000 0.750000 0.000000
Se 0.000000 0.000000 0.250000
Se 0.500000 0.000000 0.250000
Se 0.250000 0.250000 0.250000
Se 0.750000 0.250000 0.250000
Se 0.000000 0.500000 0.250000
Se 0.500000 0.500000 0.250000
Se 0.250000 0.750000 0.250000
Se 0.750000 0.750000 0.250000
Se 0.250000 0.000000 0.500000
Se 0.750000 0.000000 0.500000
Se 0.000000 0.250000 0.500000
Se 0.500000 0.250000 0.500000
Se 0.250000 0.500000 0.500000
Se 0.750000 0.500000 0.500000
Se 0.000000 0.750000 0.500000
Se 0.500000 0.750000 0.500000
Se 0.000000 0.000000 0.750000
Se 0.500000 0.000000 0.750000
Se 0.250000 0.250000 0.750000
Se 0.750000 0.250000 0.750000
Se 0.000000 0.500000 0.750000
Se 0.500000 0.500000 0.750000
Se 0.250000 0.750000 0.750000
Se 0.750000 0.750000 0.750000
K_POINTS (automatic)
8 8 8 0 0 0
Can anyone help me with this?
Thanks a lot,
Zhiting Tian
PhD student
Mechanical Engineering Department
Massachusetts Institute of Technology
Email: zhiting at mit.edu
More information about the users
mailing list