<div dir="ltr"><div dir="ltr">Hello,<div><br></div><div>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:</div><div><br></div><div><div>&CONTROL</div><div>    calculation   = "relax"</div><div>    forc_conv_thr =  1.00000e-03</div><div>    max_seconds   =  4.32000e+05</div><div>    nstep         = 300</div><div>    pseudo_dir    = "/home/coralee/projects/def-jkopysci/coralee/.pseudopot"</div><div>/</div><div><br></div><div>&SYSTEM</div><div>    a                         =  8.32716e+00</div><div>    angle1(1)                 =  0.00000e+00</div><div>    angle1(2)                 =  0.00000e+00</div><div>    angle2(1)                 =  0.00000e+00</div><div>    angle2(2)                 =  0.00000e+00</div><div>    b                         =  8.98689e+00</div><div>    c                         =  2.52767e+01</div><div>    cosab                     =  6.12323e-17</div><div>    cosac                     =  6.12323e-17</div><div>    cosbc                     = -1.85547e-01</div><div>    degauss                   =  2.00000e-02</div><div>    ecutrho                   =  4.75221e+02</div><div>    ecutwfc                   =  5.03902e+01</div><div>    ibrav                     = 12</div><div>    nat                       = 54</div><div>    nbnd                      = 510</div><div>    nspin                     = 2</div><div>    ntyp                      = 4</div><div>    occupations               = "smearing"</div><div>    smearing                  = "gaussian"</div><div>    starting_magnetization(1) =  2.00000e-01</div><div>    starting_magnetization(2) =  2.00000e-01</div><div>/</div><div><br></div><div>&ELECTRONS</div><div>    conv_thr          =  1.00000e-06</div><div>    diagonalization   = "david"</div><div>    electron_maxstep  = 250</div><div>    mixing_beta       =  6.00000e-02</div><div>    mixing_mode       = "local-TF"</div><div>    scf_must_converge = .FALSE.</div><div>    startingpot       = "atomic"</div><div>    startingwfc       = "atomic+random"</div><div>/</div><div><br></div><div>&IONS</div><div>    ion_dynamics = "bfgs"</div><div>/</div><div><br></div><div>&CELL</div><div>/</div><div><br></div><div>K_POINTS {automatic}</div><div> 2  2  1  0 0 0</div><div><br></div><div>ATOMIC_SPECIES</div><div>N      14.00674  N.pbe-n-kjpaw_psl.1.0.0.UPF</div><div>W     183.84000  W.pbe-spn-kjpaw_psl.1.0.0.UPF</div><div>C      12.01070  C.pbe-n-rrkjus_psl.1.0.0.UPF</div><div>H       1.00794  H.pbe-rrkjus_psl.1.0.0.UPF</div><div><br></div><div>ATOMIC_POSITIONS {angstrom}</div><div>N       3.208696   2.759008  14.777304</div><div>N       7.372291   2.759056  14.777299</div><div>N       3.208091  -1.734480  14.777313</div><div>N       7.371666  -1.734409  14.777266</div><div>W       3.206835   4.918476  14.576735</div><div>W       7.370342   4.918461  14.576681</div><div>W       3.207409   0.424961  14.576616</div><div>W       7.371008   0.425000  14.576555</div><div>N       1.125214   5.237649  14.647393</div><div>N       5.288843   5.237678  14.647535</div><div>N       1.125689   0.744163  14.647267</div><div>N       5.289301   0.744198  14.647391</div><div>W       1.127155   3.039266  14.427235</div><div>W       5.290708   3.039279  14.427235</div><div>W       1.126340  -1.454299  14.427190</div><div>W       5.289930  -1.454181  14.427145</div><div>N       3.207623   5.190064  12.480026</div><div>N       7.371213   5.190146  12.479947</div><div>N       3.207739   0.696664  12.479837</div><div>N       7.371242   0.696568  12.479880</div><div>W       3.207774   3.077014  12.539874</div><div>W       7.371320   3.077042  12.539864</div><div>W       3.207757  -1.416448  12.539825</div><div>W       7.371342  -1.416485  12.539806</div><div>N       1.125882   3.569158  12.369997</div><div>N       5.289399   3.569247  12.369990</div><div>N       1.125941  -0.924226  12.369961</div><div>N       5.289548  -0.924284  12.369791</div><div>W       1.125927   5.677973  12.167315</div><div>W       5.289472   5.678004  12.167326</div><div>W       1.125908   1.184568  12.167344</div><div>W       5.289510   1.184576  12.167245</div><div>N       3.207882   3.506729  10.422372  0 0 0</div><div>N       7.371462   3.506729  10.422372  0 0 0</div><div>N       3.207882  -0.986715  10.422372  0 0 0</div><div>N       7.371462  -0.986715  10.422372  0 0 0</div><div>W       3.207882   5.716584  10.333236  0 0 0</div><div>W       7.371462   5.716584  10.333236  0 0 0</div><div>W       3.207882   1.223140  10.333236  0 0 0</div><div>W       7.371462   1.223140  10.333236  0 0 0</div><div>N       1.126093   6.199382  10.052766  0 0 0</div><div>N       5.289672   6.199382  10.052766  0 0 0</div><div>N       1.126093   1.705938  10.052766  0 0 0</div><div>N       5.289672   1.705938  10.052766  0 0 0</div><div>W       1.126093   3.871912  10.000000  0 0 0</div><div>W       5.289672   3.871912  10.000000  0 0 0</div><div>W       1.126093  -0.621532  10.000000  0 0 0</div><div>W       5.289672  -0.621532  10.000000  0 0 0</div><div>C       4.649862   5.094916  16.119430</div><div>H       4.905618   4.303189  16.752099</div><div>H       4.836000   5.987155  16.769592</div><div>C       3.816805   2.863195  16.169239</div><div>H       4.072561   1.971468  16.801908</div><div>H       4.002943   3.655434  16.819401</div><div><br></div><div><br></div></div><div><br></div><div><br></div></div></div>