<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>