<div dir="ltr"><div dir="ltr">Hello,<div><br></div><div>My structure is a WN slab with molecules adsorbing on its surface. I always have problems with convergence: everytime I try to relax the structure, multiple scf steps converge, but then the energy value starts oscillating and it doesn't converge anymore. My structure gets stuck every time, even though I visualized it after every scf step and the structure seems fine. This is the input file I am using. I also tried decreasing the degauss to 1e-2, increasing ecutrho and ecutwfc by a little bit, decreasing mixing_beta to 6e-2 and also tried changing the mixing_mode to TF, but nothing worked! Does anyone have advice? Thank you!</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         = 200</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>    angle2(1)                 =  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                       = 53</div><div>    nbnd                      = 480</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>    starting_magnetization(3) =  0.00000e+00</div><div>    starting_magnetization(4) =  0.00000e+00</div><div>/</div><div><br></div><div>&ELECTRONS</div><div>    conv_thr         =  1.00000e-05</div><div>    diagonalization  = "david"</div><div>    electron_maxstep = 200</div><div>    mixing_beta      =  1.00000e-01</div><div>    mixing_mode      = "local-TF"</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-kjpaw_psl.1.0.0.UPF</div><div>H       1.00794  H.pbe-kjpaw_psl.1.0.0.UPF</div><div><br></div><div>ATOMIC_POSITIONS {angstrom}</div><div>N       3.210864   2.727236  14.728253</div><div>N       7.377521   2.749437  14.765113</div><div>N       3.206885  -1.748226  14.750672</div><div>N       7.368425  -1.739825  14.763088</div><div>W       3.182078   4.921673  14.584931</div><div>W       7.371236   4.911599  14.572299</div><div>W       3.210468   0.427880  14.572328</div><div>W       7.375435   0.429218  14.573322</div><div>N       1.114319   5.237998  14.590842</div><div>N       5.297509   5.253239  14.643339</div><div>N       1.127456   0.739579  14.633097</div><div>N       5.291670   0.741857  14.637218</div><div>W       1.130917   3.035176  14.411937</div><div>W       5.301873   3.012738  14.402683</div><div>W       1.122976  -1.443293  14.413944</div><div>W       5.288433  -1.443790  14.404037</div><div>N       3.219452   5.199606  12.448287</div><div>N       7.366466   5.188650  12.437657</div><div>N       3.207970   0.704887  12.465361</div><div>N       7.371281   0.697722  12.461527</div><div>W       3.205392   3.068813  12.535650</div><div>W       7.371059   3.075963  12.535225</div><div>W       3.208342  -1.417371  12.526978</div><div>W       7.371423  -1.415702  12.526259</div><div>N       1.126718   3.566876  12.355498</div><div>N       5.287387   3.570864  12.349260</div><div>N       1.127573  -0.926625  12.347835</div><div>N       5.290289  -0.925589  12.354591</div><div>W       1.124934   5.676097  12.168934</div><div>W       5.300233   5.687505  12.124192</div><div>W       1.125212   1.183514  12.166731</div><div>W       5.291832   1.185520  12.166744</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.904436   5.088998  16.341705</div><div>H       4.817236   4.340738  17.152672</div><div>H       4.810563   5.831388  17.155825</div><div>H       3.209512   2.314472  17.745318</div><div>H       3.207503   3.080782  17.746006</div></div><div><br></div></div></div>