<div dir="ltr"><div dir="ltr"><div dir="ltr">Hello,<div><br></div><div>I am performing a relaxation calculation on my Tungsten nitride (WN) slab with a methane molecule adsorbed on its surface. My first scf calculations always converges, but all the other scf calculations keep oscillating a little bit, like this:</div><div><br></div><div><div>iteration #107 Â  Â  ecut= Â  Â 50.39 Ry Â  Â  beta=0.10</div><div>  Â  Â Davidson diagonalization with overlap</div><div>  Â  Â ethr = Â 1.00E-02, Â avg # of iterations = Â 3.2</div><div><br></div><div>  Â  Â negative rho (up, down): Â 1.877E+02 1.832E+02</div><div><br></div><div>  Â  Â total cpu time spent up to now is Â  Â  6820.6 secs</div><div><br></div><div>  Â  Â total energy Â  Â  Â  Â  Â  Â  Â = Â <span style="background-color:rgb(255,255,0)">-22310.34633946 Ry</span></div><div>  Â  Â Harris-Foulkes estimate Â  = Â -22286.51427887 Ry</div><div>  Â  Â estimated scf accuracy Â  Â < Â  36291.44935321 Ry</div><div><br></div><div>  Â  Â total magnetization Â  Â  Â  = Â  Â  0.88 Bohr mag/cell</div><div>  Â  Â absolute magnetization Â  Â = Â  Â 13.09 Bohr mag/cell</div><div><br></div><div>  Â  Â iteration #108 Â  Â  ecut= Â  Â 50.39 Ry Â  Â  beta=0.10</div><div>  Â  Â Davidson diagonalization with overlap</div><div>  Â  Â ethr = Â 1.00E-02, Â avg # of iterations = Â 2.0</div><div><br></div><div>  Â  Â negative rho (up, down): Â 1.828E+02 1.844E+02</div><div><br></div><div>  Â  Â total cpu time spent up to now is Â  Â  6849.8 secs</div><div><br></div><div>  Â  Â total energy Â  Â  Â  Â  Â  Â  Â = Â <span style="background-color:rgb(255,255,0)">-22482.46538384</span> Ry</div><div>  Â  Â Harris-Foulkes estimate Â  = Â -22312.40216721 Ry</div><div>  Â  Â estimated scf accuracy Â  Â < Â  36350.61037433 Ry</div><div><br></div><div>  Â  Â total magnetization Â  Â  Â  = Â  Â  0.84 Bohr mag/cell</div><div>  Â  Â absolute magnetization Â  Â = Â  Â 11.93 Bohr mag/cell</div><div><br></div><div>  Â  Â iteration #109 Â  Â  ecut= Â  Â 50.39 Ry Â  Â  beta=0.10</div><div>  Â  Â Davidson diagonalization with overlap</div><div>  Â  Â ethr = Â 1.00E-02, Â avg # of iterations = Â 3.0</div><div><br></div><div>  Â  Â negative rho (up, down): Â 1.848E+02 1.794E+02</div><div><br></div><div>  Â  Â total cpu time spent up to now is Â  Â  6888.4 secs</div><div><br></div><div>  Â  Â total energy Â  Â  Â  Â  Â  Â  Â = Â -22531.89474337 Ry</div><div>  Â  Â Harris-Foulkes estimate Â  = Â -22546.66368909 Ry</div><div>  Â  Â estimated scf accuracy Â  Â < Â  34615.41965371 Ry</div><div><br></div><div>  Â  Â total magnetization Â  Â  Â  = Â  Â  0.49 Bohr mag/cell</div><div>  Â  Â absolute magnetization Â  Â = Â  Â  8.33 Bohr mag/cell</div><div><br></div><div>  Â  Â iteration #110 Â  Â  ecut= Â  Â 50.39 Ry Â  Â  beta=0.10</div><div>  Â  Â Davidson diagonalization with overlap</div><div>  Â  Â ethr = Â 1.00E-02, Â avg # of iterations = Â 3.5</div><div><br></div><div>  Â  Â negative rho (up, down): Â 1.882E+02 1.776E+02</div><div><br></div><div>  Â  Â total cpu time spent up to now is Â  Â  6947.2 secs</div><div><br></div><div>  Â  Â total energy Â  Â  Â  Â  Â  Â  Â = Â -22697.64012878 Ry</div><div>  Â  Â Harris-Foulkes estimate Â  = Â -22561.26019900 Ry</div><div>  Â  Â estimated scf accuracy Â  Â < Â  35110.57986884 Ry</div><div><br></div><div>  Â  Â total magnetization Â  Â  Â  = Â  Â  1.02 Bohr mag/cell</div><div>  Â  Â absolute magnetization Â  Â = Â  Â  8.80 Bohr mag/cell</div><div><br></div><div>  Â  Â iteration #111 Â  Â  ecut= Â  Â 50.39 Ry Â  Â  beta=0.10</div><div>  Â  Â Davidson diagonalization with overlap</div><div>  Â  Â ethr = Â 1.00E-02, Â avg # of iterations = Â 3.0</div><div><br></div><div>  Â  Â negative rho (up, down): Â 1.953E+02 1.846E+02</div><div><br></div><div>  Â  Â total cpu time spent up to now is Â  Â  6984.5 secs</div><div><br></div><div>  Â  Â total energy Â  Â  Â  Â  Â  Â  Â = Â -22595.85493501 Ry</div><div>  Â  Â Harris-Foulkes estimate Â  = Â -22716.54469069 Ry</div><div>  Â  Â estimated scf accuracy Â  Â < Â  35636.87841917 Ry</div><div><br></div><div>  Â  Â total magnetization Â  Â  Â  = Â  Â  1.33 Bohr mag/cell</div><div>  Â  Â absolute magnetization Â  Â = Â  Â 11.84 Bohr mag/cell</div><div><br></div><div>  Â  Â iteration #112 Â  Â  ecut= Â  Â 50.39 Ry Â  Â  beta=0.10</div><div>  Â  Â Davidson diagonalization with overlap</div><div>  Â  Â ethr = Â 1.00E-02, Â avg # of iterations = Â 3.0</div><div><br></div><div>  Â  Â negative rho (up, down): Â 2.248E+02 1.996E+02</div><div><br></div><div>  Â  Â total cpu time spent up to now is Â  Â  7023.3 secs</div><div><br></div><div>  Â  Â total energy Â  Â  Â  Â  Â  Â  Â = Â -22958.24642440 Ry</div><div>  Â  Â Harris-Foulkes estimate Â  = Â -22626.00971729 Ry</div><div>  Â  Â estimated scf accuracy Â  Â < Â  38324.66640261 Ry</div><div><br></div><div>  Â  Â total magnetization Â  Â  Â  = Â  Â  1.31 Bohr mag/cell</div><div>  Â  Â absolute magnetization Â  Â = Â  Â 10.12 Bohr mag/cell</div><div><br></div><div>  Â  Â iteration #113 Â  Â  ecut= Â  Â 50.39 Ry Â  Â  beta=0.10</div><div>  Â  Â Davidson diagonalization with overlap</div><div>  Â  Â ethr = Â 1.00E-02, Â avg # of iterations = Â 3.0</div><div><br></div><div>  Â  Â negative rho (up, down): Â 2.614E+02 2.124E+02</div><div><br></div><div>  Â  Â total cpu time spent up to now is Â  Â  7070.7 secs</div><div><br></div><div>  Â  Â total energy Â  Â  Â  Â  Â  Â  Â = Â -23388.04862201 Ry</div><div>  Â  Â Harris-Foulkes estimate Â  = Â -23409.72136428 Ry</div><div>  Â  Â estimated scf accuracy Â  Â < Â  40154.38739523 Ry</div><div><br></div><div>  Â  Â total magnetization Â  Â  Â  = Â  Â  1.90 Bohr mag/cell</div><div>  Â  Â absolute magnetization Â  Â = Â  Â 15.87 Bohr mag/cell</div><div><br></div><div>  Â  Â iteration #114 Â  Â  ecut= Â  Â 50.39 Ry Â  Â  beta=0.10</div><div>  Â  Â Davidson diagonalization with overlap</div><div>  Â  Â ethr = Â 1.00E-02, Â avg # of iterations = Â 3.5</div><div><br></div><div>  Â  Â negative rho (up, down): Â 2.824E+02 2.176E+02</div><div><br></div><div>  Â  Â total cpu time spent up to now is Â  Â  7120.5 secs</div><div><br></div><div>  Â  Â total energy Â  Â  Â  Â  Â  Â  Â = Â -23737.17074269 Ry</div><div>  Â  Â Harris-Foulkes estimate Â  = Â -23665.30705552 Ry</div><div>  Â  Â estimated scf accuracy Â  Â < Â  35266.44456712 Ry</div><div><br></div><div>  Â  Â total magnetization Â  Â  Â  = Â  Â  2.56 Bohr mag/cell</div><div>  Â  Â absolute magnetization Â  Â = Â  Â 16.20 Bohr mag/cell</div></div><div><br></div><div><b>How can I make my scf calculation converge?</b></div><div><br></div><div>This is my <b>input</b> 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 Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  = 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>/</div><div><br></div><div>&ELECTRONS</div><div>  Â  conv_thr Â  Â  Â  Â  Â = Â 1.00000e-05</div><div>  Â  diagonalization Â  = "david"</div><div>  Â  electron_maxstep Â = 500</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.203704 Â  2.747410 Â 14.771098</div><div>N Â  Â  Â  7.372735 Â  2.759212 Â 14.776453</div><div>N Â  Â  Â  3.208642 Â -1.731177 Â 14.774215</div><div>N Â  Â  Â  7.371486 Â -1.734287 Â 14.777353</div><div>W Â  Â  Â  3.214859 Â  4.922015 Â 14.312292</div><div>W Â  Â  Â  7.369397 Â  4.914334 Â 14.580589</div><div>W Â  Â  Â  3.207301 Â  0.425737 Â 14.576255</div><div>W Â  Â  Â  7.371014 Â  0.424829 Â 14.576439</div><div>N Â  Â  Â  1.113188 Â  5.239676 Â 14.637084</div><div>N Â  Â  Â  5.300454 Â  5.246400 Â 14.637081</div><div>N Â  Â  Â  1.125960 Â  0.744616 Â 14.647342</div><div>N Â  Â  Â  5.289215 Â  0.745329 Â 14.647446</div><div>W Â  Â  Â  1.122438 Â  3.033042 Â 14.425674</div><div>W Â  Â  Â  5.295367 Â  3.033769 Â 14.426077</div><div>W Â  Â  Â  1.124035 Â -1.452510 Â 14.426388</div><div>W Â  Â  Â  5.289229 Â -1.455341 Â 14.427007</div><div>N Â  Â  Â  3.210511 Â  5.196004 Â 12.452975</div><div>N Â  Â  Â  7.370831 Â  5.189844 Â 12.478804</div><div>N Â  Â  Â  3.207641 Â  0.696907 Â 12.479359</div><div>N Â  Â  Â  7.371300 Â  0.696480 Â 12.479853</div><div>W Â  Â  Â  3.206497 Â  3.070844 Â 12.536437</div><div>W Â  Â  Â  7.370775 Â  3.078954 Â 12.541331</div><div>W Â  Â  Â  3.207829 Â -1.415777 Â 12.539146</div><div>W Â  Â  Â  7.371280 Â -1.416331 Â 12.539239</div><div>N Â  Â  Â  1.125925 Â  3.569153 Â 12.369437</div><div>N Â  Â  Â  5.289946 Â  3.569319 Â 12.367599</div><div>N Â  Â  Â  1.125855 Â -0.924614 Â 12.368587</div><div>N Â  Â  Â  5.289457 Â -0.924540 Â 12.368612</div><div>W Â  Â  Â  1.125854 Â  5.677831 Â 12.169785</div><div>W Â  Â  Â  5.290755 Â  5.676885 Â 12.165625</div><div>W Â  Â  Â  1.126115 Â  1.184835 Â 12.168057</div><div>W Â  Â  Â  5.289247 Â  1.185090 Â 12.168067</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 Â  Â  Â  3.206835 Â  4.918476 Â 16.276735 Â 0 0 0</div><div>H Â  Â  Â  2.573835 Â  4.285476 Â 16.909735 Â 0 0 0</div><div>H Â  Â  Â  3.839835 Â  4.285476 Â 15.643735 Â 0 0 0</div><div>H Â  Â  Â  2.573835 Â  5.551476 Â 15.643735 Â 0 0 0</div><div>H Â  Â  Â  3.839835 Â  5.551476 Â 16.909735 Â 0 0 0</div><div><br></div><div><br></div></div><div><br></div></div></div></div>