[QE-users] relaxation calculation not converging
Coralie Khabbaz
khabbaz.coralie at gmail.com
Fri Jul 24 23:04:27 CEST 2020
Hello,
I am performing a relaxation calculation on a tungsten nitride catalyst
with molecules that are adsorbing on its surface. My system is difficult to
converge, and since I added the adsorbed molecules it has not been
converging (scf calculations converge, it's just the relaxation
calculation). I have tried varying the degauss, the ecutrho and the
ecutwfc, the nbnd, decreasing the conv_thr, changing mixing_beta, but
nothing worked. Can someone help me make it converge? This is my input file:
&CONTROL
calculation = "relax"
forc_conv_thr = 1.00000e-03
max_seconds = 4.32000e+05
nstep = 300
pseudo_dir = "/home/coralee/projects/def-jkopysci/coralee/.pseudopot"
/
&SYSTEM
a = 8.32716e+00
angle1(1) = 0.00000e+00
angle1(2) = 0.00000e+00
angle2(1) = 0.00000e+00
angle2(2) = 0.00000e+00
b = 8.98689e+00
c = 2.52767e+01
cosab = 6.12323e-17
cosac = 6.12323e-17
cosbc = -1.85547e-01
degauss = 2.00000e-02
ecutrho = 4.75221e+02
ecutwfc = 5.03902e+01
ibrav = 12
nat = 54
nbnd = 510
nspin = 2
ntyp = 4
occupations = "smearing"
smearing = "gaussian"
starting_magnetization(1) = 2.00000e-01
starting_magnetization(2) = 2.00000e-01
/
&ELECTRONS
conv_thr = 1.00000e-06
diagonalization = "david"
electron_maxstep = 250
mixing_beta = 6.00000e-02
mixing_mode = "local-TF"
scf_must_converge = .FALSE.
startingpot = "atomic"
startingwfc = "atomic+random"
/
&IONS
ion_dynamics = "bfgs"
/
&CELL
/
K_POINTS {automatic}
2 2 1 0 0 0
ATOMIC_SPECIES
N 14.00674 N.pbe-n-kjpaw_psl.1.0.0.UPF
W 183.84000 W.pbe-spn-kjpaw_psl.1.0.0.UPF
C 12.01070 C.pbe-n-rrkjus_psl.1.0.0.UPF
H 1.00794 H.pbe-rrkjus_psl.1.0.0.UPF
ATOMIC_POSITIONS {angstrom}
N 3.208696 2.759008 14.777304
N 7.372291 2.759056 14.777299
N 3.208091 -1.734480 14.777313
N 7.371666 -1.734409 14.777266
W 3.206835 4.918476 14.576735
W 7.370342 4.918461 14.576681
W 3.207409 0.424961 14.576616
W 7.371008 0.425000 14.576555
N 1.125214 5.237649 14.647393
N 5.288843 5.237678 14.647535
N 1.125689 0.744163 14.647267
N 5.289301 0.744198 14.647391
W 1.127155 3.039266 14.427235
W 5.290708 3.039279 14.427235
W 1.126340 -1.454299 14.427190
W 5.289930 -1.454181 14.427145
N 3.207623 5.190064 12.480026
N 7.371213 5.190146 12.479947
N 3.207739 0.696664 12.479837
N 7.371242 0.696568 12.479880
W 3.207774 3.077014 12.539874
W 7.371320 3.077042 12.539864
W 3.207757 -1.416448 12.539825
W 7.371342 -1.416485 12.539806
N 1.125882 3.569158 12.369997
N 5.289399 3.569247 12.369990
N 1.125941 -0.924226 12.369961
N 5.289548 -0.924284 12.369791
W 1.125927 5.677973 12.167315
W 5.289472 5.678004 12.167326
W 1.125908 1.184568 12.167344
W 5.289510 1.184576 12.167245
N 3.207882 3.506729 10.422372 0 0 0
N 7.371462 3.506729 10.422372 0 0 0
N 3.207882 -0.986715 10.422372 0 0 0
N 7.371462 -0.986715 10.422372 0 0 0
W 3.207882 5.716584 10.333236 0 0 0
W 7.371462 5.716584 10.333236 0 0 0
W 3.207882 1.223140 10.333236 0 0 0
W 7.371462 1.223140 10.333236 0 0 0
N 1.126093 6.199382 10.052766 0 0 0
N 5.289672 6.199382 10.052766 0 0 0
N 1.126093 1.705938 10.052766 0 0 0
N 5.289672 1.705938 10.052766 0 0 0
W 1.126093 3.871912 10.000000 0 0 0
W 5.289672 3.871912 10.000000 0 0 0
W 1.126093 -0.621532 10.000000 0 0 0
W 5.289672 -0.621532 10.000000 0 0 0
C 4.649862 5.094916 16.119430
H 4.905618 4.303189 16.752099
H 4.836000 5.987155 16.769592
C 3.816805 2.863195 16.169239
H 4.072561 1.971468 16.801908
H 4.002943 3.655434 16.819401
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20200724/34f0bb10/attachment.html>
More information about the users
mailing list