<div dir="ltr">Meaning not too large. Large forces would lead to bad relaxation to unrealistic structure that cause convergence problems<div><br></div><div>Duy Le</div><div>(UCF)</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Jul 10, 2020 at 1:13 PM Coralie Khabbaz <<a href="mailto:khabbaz.coralie@gmail.com">khabbaz.coralie@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hello Duy Le,<div><br></div><div>Thank you so much for your answer. What do you mean by reasonable forces? I have performed several scf calculations before and I get energy values that are consistent. </div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, 10 Jul 2020 at 13:08, Duy Le <<a href="mailto:ttduyle@gmail.com" target="_blank">ttduyle@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">It could be because of non-realistic structure as a result of relaxation. Please check the force in previously converged scf to see if they are reasonable. If not you may need to remake the starting structure so it has reasonable forces.<div><br></div><div>Duy Le</div><div>(UCF)</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Jul 10, 2020 at 12:53 PM Coralie Khabbaz <<a href="mailto:khabbaz.coralie@gmail.com" target="_blank">khabbaz.coralie@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Also, when I perform scf calculation (and not relaxation), my energy value converges!!</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, 10 Jul 2020 at 12:30, Coralie Khabbaz <<a href="mailto:khabbaz.coralie@gmail.com" target="_blank">khabbaz.coralie@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><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>
</blockquote></div>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div>
</blockquote></div>