<div dir="ltr"><div dir="ltr">Hello Matteo,<div><br></div><div>Thank you so much for your reply. I didn't understand what do you mean by initial trust radius? Also, at first my convergence threshold was at 1^-6 but the first scf step didn't converge in 1000 steps. Then, I fixed it at 1^-4 and the first step converged. So do you think I should tighten it? I didn't try starting a new scf from the new iconic positions because I don't know how to. </div><div><br></div><div>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                     = 14</div><div>    nat                       = 48</div><div>    nbnd                      = 490</div><div>    nspin                     = 2</div><div>    ntyp                      = 2</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-04</div><div>    diagonalization  = "david"</div><div>    electron_maxstep = 1000</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><br></div><div>ATOMIC_POSITIONS {angstrom}</div><div>N       3.207882   2.672960  14.837784</div><div>N       7.371462   2.672960  14.837784</div><div>N       3.207882  -1.820484  14.837784</div><div>N       7.371462  -1.820484  14.837784</div><div>W       3.207882   4.882815  14.748648</div><div>W       7.371462   4.882815  14.748648</div><div>W       3.207882   0.389371  14.748648</div><div>W       7.371462   0.389371  14.748648</div><div>N       1.126093   5.365612  14.468179</div><div>N       5.289672   5.365612  14.468179</div><div>N       1.126093   0.872168  14.468179</div><div>N       5.289672   0.872168  14.468179</div><div>W       1.126093   3.038143  14.415413</div><div>W       5.289672   3.038143  14.415413</div><div>W       1.126093  -1.455301  14.415413</div><div>W       5.289672  -1.455301  14.415413</div><div>N       3.207882   5.336567  12.630078</div><div>N       7.371462   5.336567  12.630078</div><div>N       3.207882   0.843123  12.630078</div><div>N       7.371462   0.843123  12.630078</div><div>W       3.207882   3.052977  12.540942</div><div>W       7.371462   3.052977  12.540942</div><div>W       3.207882  -1.440467  12.540942</div><div>W       7.371462  -1.440467  12.540942</div><div>N       1.126093   3.535775  12.260473</div><div>N       5.289672   3.535775  12.260473</div><div>N       1.126093  -0.957669  12.260473</div><div>N       5.289672  -0.957669  12.260473</div><div>W       1.126093   5.701750  12.207706</div><div>W       5.289672   5.701750  12.207706</div><div>W       1.126093   1.208305  12.207706</div><div>W       5.289672   1.208305  12.207706</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><br></div><div><br></div></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, 1 Jul 2020 at 10:48, Matteo Cococcioni <<a href="mailto:matteo.cococcioni@unipv.it">matteo.cococcioni@unipv.it</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hello,</div><div><br></div><div>it's very difficult to guess what is wrong in your input if you don't show it.</div><div>In any case, if you have converged forces and stresses on cut-off(s) and number of k-points/smearing, it might be a good idea to tighten your convergence threshold and reduce the initial trust radius. Maybe your forces and stresses are so loosely converged and your ionic/cell step so big that the system moves to a configuration that is difficult to converge. Or maybe your starting lattice parameter/ionic configuration are very far from equilibrium and the system tries to make a big jump at the beginning. Also, have you tried to start a new scf from the new ionic position?</div><div><br></div><div>Please add your affiliation</div><div><br></div><div>Regards,</div><div><br></div><div>Matteo<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno mer 1 lug 2020 alle ore 13:45 Coralie Khabbaz <<a href="mailto:khabbaz.coralie@gmail.com" target="_blank">khabbaz.coralie@gmail.com</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hello,<div><br></div><div>I am doing a relax calculation on my WN (48 atoms) slab. The first step of scf calculations converged but then the second step didn't converge in 1000 steps (which was my limit). So, the calculations stopped, and I got many files (data-file.xml, charge-density.dat, paw.txt, spin-polarization.dat...) but no output file. If my first scf step converged doesn't that mean that all of my other scf steps should converge but need more steps? Or is there a problem in my input file?</div></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 clear="all"><br>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>Matteo Cococcioni<br></div><div>Department of Physics<br></div>University of Pavia<br>Via Bassi 6, I-27100 Pavia, Italy<br>tel +39-0382-987485<br><div>e-mail <a href="mailto:lucio.andreani@unipv.it" target="_blank">matteo.cococcioni@unipv.it</a></div></div></div></div></div></div></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>