<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Hello Marko,<div class=""><br class=""></div><div class="">Thank you for your prompt response.<br class=""><div class=""><br class=""></div><div class="">I copy-pasted the wrong input file. The right one is a SCF one (the nscf should be subsequent to scf):</div><div class=""><div class="">&CONTROL</div><div class=""> calculation ='scf',</div><div class=""> nstep = 200,</div><div class=""> etot_conv_thr = 1.d-7,</div><div class="">! forc_conv_thr = 1.d-4,</div><div class=""> wfcdir = './WFC' ,</div><div class="">! tstress = .true.,</div><div class=""> prefix = 'STe_scf_LS',</div><div class=""> pseudo_dir = '/scratch/cnt0022/pmc6881/paboulet/pseudo/',</div><div class=""> verbosity = 'high',</div><div class=""> restart_mode = ‘from_scratch',</div><div class=""> wf_collect = .true.,</div><div class=""> disk_io = 'high',</div></div><div class=""><div class=""><br class="webkit-block-placeholder"></div><div class="">The rest is the same as in my previous email, except for the k-points: 12x12x4 1 1 1, as in the non-SO optimisation.</div><div class=""><br class=""></div><div class="">Pascal</div><div class=""><br class=""></div>
<div><br class=""><blockquote type="cite" class=""><div class="">Le 19 déc. 2019 à 17:42, Mladenovic Marko <<a href="mailto:marko.mladenovic@epfl.ch" class="">marko.mladenovic@epfl.ch</a>> a écrit :</div><br class="Apple-interchange-newline"><div class=""><div id="divtagdefaultwrapper" dir="ltr" style="caret-color: rgb(0, 0, 0); font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-size: 12pt; font-family: Calibri, Helvetica, sans-serif;" class=""><div id="divtagdefaultwrapper" dir="ltr" style="font-size: 12pt; font-family: Calibri, Helvetica, sans-serif, EmojiFont, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols;" class=""><div style="margin-top: 0px; margin-bottom: 0px;" class="">Hello Pascal,</div><div style="margin-top: 0px; margin-bottom: 0px;" class=""><br class=""></div><div style="margin-top: 0px; margin-bottom: 0px;" class="">it seems that you are performing a non-scf calculation (<span class="">calculation ='nscf'</span>) and you are restarting an other calculation (<span class="">restart_mode = 'restart'</span>). Is this what you want?</div><div style="margin-top: 0px; margin-bottom: 0px;" class=""><br class=""></div><div style="margin-top: 0px; margin-bottom: 0px;" class="">Best,</div><div style="margin-top: 0px; margin-bottom: 0px;" class="">Marko<br class=""></div></div><hr tabindex="-1" style="display: inline-block; width: 733.03125px;" class=""><div id="divRplyFwdMsg" dir="ltr" class=""><font face="Calibri, sans-serif" style="font-size: 11pt;" class=""><b class="">From:</b><span class="Apple-converted-space"> </span>users <<a href="mailto:users-bounces@lists.quantum-espresso.org" class="">users-bounces@lists.quantum-espresso.org</a>> on behalf of pboulet <<a href="mailto:pascal.boulet@univ-amu.fr" class="">pascal.boulet@univ-amu.fr</a>><br class=""><b class="">Sent:</b><span class="Apple-converted-space"> </span>Thursday, December 19, 2019 5:03:40 PM<br class=""><b class="">To:</b><span class="Apple-converted-space"> </span>Quantum Espresso users Forum<br class=""><b class="">Subject:</b><span class="Apple-converted-space"> </span>[QE-users] Problem with SCF including spin-orbit coupling</font><div class=""> </div></div><div class="">Dear All,<div class=""><br class=""></div><div class="">I have some troubles with a SCF calculation using SOC (QE 6.1).</div><div class="">I have optimised a structure with non-SOC pseudopotentials and now I want to perform a SCF calculation by including SOC to get the electronic band structure. The problem is that the SCF diverges right away as you can see here:</div><div class="">>>>></div><div class=""><div class="">Self-consistent Calculation</div><div class=""><br class=""></div><div class=""> iteration # 1 ecut= 35.00 Ry beta=0.10</div><div class=""> Davidson diagonalization with overlap</div><div class=""> ethr = 1.00E-05, avg # of iterations = 14.2</div><div class=""><br class=""></div><div class=""> negative rho (up, down): 8.197E+02 0.000E+00</div><div class=""><br class=""></div><div class=""> total cpu time spent up to now is 2293.7 secs</div><div class=""><br class=""></div><div class=""> total energy = -16093.90449193 Ry</div><div class=""> Harris-Foulkes estimate = -18277.15379034 Ry</div><div class=""> estimated scf accuracy < 5077678.52441197 Ry</div><div class=""><br class=""></div><div class=""> iteration # 2 ecut= 35.00 Ry beta=0.10</div><div class=""> Davidson diagonalization with overlap</div><div class=""> ethr = 1.00E-02, avg # of iterations = 5.0</div><div class=""><br class=""></div><div class=""> negative rho (up, down): 2.677E+03 0.000E+00</div><div class=""><br class=""></div><div class=""> total cpu time spent up to now is 3119.6 secs</div><div class=""><br class=""></div><div class=""> total energy = 79024.56634461 Ry</div><div class=""> Harris-Foulkes estimate = -17101.52000727 Ry</div><div class=""> estimated scf accuracy < 4036754.22296032 Ry</div><div class=""><br class=""></div><div class=""> iteration # 3 ecut= 35.00 Ry beta=0.10</div><div class=""> Davidson diagonalization with overlap</div><div class=""> ethr = 1.00E-02, avg # of iterations = 6.0</div><div class=""><br class=""></div><div class=""> negative rho (up, down): 2.815E+03 0.000E+00</div><div class=""><br class=""></div><div class=""> total cpu time spent up to now is 4658.4 secs</div><div class=""><br class=""></div><div class=""> total energy = 132135.46722241 Ry</div></div><div class=""><div class=""> Harris-Foulkes estimate = -142275.32127673 Ry</div><div class=""> estimated scf accuracy < 6116470.38137350 Ry</div></div><div class=""><<<<<</div><div class=""> </div><div class="">I am using PAW pseudopotentials (Ba, Rh, Ge and S) from the PSL library. The cg diagonalisation does not solve the problem. </div><div class=""><br class=""></div><div class="">Here is a piece of the input file:</div><div class=""><div class="">&CONTROL</div><div class=""> calculation ='nscf',</div><div class=""> nstep = 300,</div><div class=""> etot_conv_thr = 1.d-7,</div><div class=""> forc_conv_thr = 1.d-4,</div><div class=""> wfcdir = './WFC' ,</div><div class=""> prefix = 'STe_nscf_LS',</div><div class=""> pseudo_dir = '/scratch/cnt0022/pmc6881/paboulet/pseudo/',</div><div class=""> verbosity = 'high',</div><div class=""> restart_mode = 'restart',</div><div class=""> wf_collect = .true.,</div><div class=""> disk_io = 'high',</div><div class="">/</div><div class="">&SYSTEM</div><div class="">! celldm(1) = 1.0,</div><div class=""> nat = 52,</div><div class=""> ntyp = 5,</div><div class=""> ibrav = 0,</div><div class=""> ecutwfc = 35.d0,</div><div class=""> ecutrho = 350.d0,</div><div class=""> occupations = 'fixed',</div><div class=""> nbnd = 630,</div><div class=""> lspinorb = .true.,</div><div class=""> noncolin = .true.,</div><div class=""> nr1 = 72, nr2 = 72, nr3 = 360,</div><div class="">/</div></div><div class=""><div class="">&ELECTRONS</div><div class=""> electron_maxstep = 200,</div><div class=""> conv_thr = 1.d-10,</div><div class=""> mixing_beta = 0.2d0,</div><div class=""> diagonalization = 'david',</div><div class="">/</div><div class=""><br class=""></div><div class=""><div class="">K_POINTS automatic</div><div class="">16 16 8 0 0 0</div><div class=""><br class=""></div><div class=""><div class="">ATOMIC_SPECIES</div><div class="">Ba 137.3270 Ba.rel-pbe-spn-kjpaw_psl.1.0.0.UPF,</div><div class="">Rh 102.9055 Rh.rel-pbe-spn-kjpaw_psl.1.0.0.UPF,</div><div class="">Ge 72.6400 Ge.rel-pbe-dn-kjpaw_psl.1.0.0.UPF,</div><div class="">S 32.0650 S.rel-pbe-nl-kjpaw_psl.1.0.0.UPF,</div><div class="">Te 127.6000 Te.rel-pbe-dn-kjpaw_psl.0.2.2.UPF</div></div><div class=""><br class=""></div><div class="">etc…</div><div class=""><br class=""></div><div class=""><br class=""></div></div></div><div class="">Do you have a special recipe for this type of calculation?</div><div class=""><br class=""></div><div class="">Thank you,</div><div class="">Best regards</div><div class=""><br class=""></div><div class=""><div dir="auto" class="" style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; word-wrap: break-word; line-break: after-white-space;"><div dir="auto" class="" style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; word-wrap: break-word; line-break: after-white-space;"><div dir="auto" class="" style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; word-wrap: break-word; line-break: after-white-space;"><div dir="auto" class="" style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; word-wrap: break-word; line-break: after-white-space;"><div dir="auto" class="" style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; word-wrap: break-word; line-break: after-white-space;"><div dir="auto" class="" style="text-align: start; text-indent: 0px; word-wrap: break-word; line-break: after-white-space;"><div class="" style="letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none;"><span class="" style="font-size: 16pt; font-family: Mistral;">Pascal Boulet</span></div><div class="" style="letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none;"><font class="" face="Lucida Handwriting"><span class="" style="font-size: 21px;">—</span></font></div><div class="" style="letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none;"><em class="" style="color: rgb(34, 187, 234); font-family: verdana; font-size: 11px; line-height: 22px;">Professor in computational materials - DEPARTMENT OF CHEMISTRY</em></div><div class=""><div class="" style="letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; font-family: Tahoma, Verdana, Arial, Helvetica, sans-serif; font-size: 15px; line-height: 22px;"><span class="" style="font-size: 11px; font-family: verdana;">Aix-Marseille University </span><span class="" style="font-size: 11px; font-family: verdana;">- Avenue Escadrille Normandie Niemen - F-13013 Marseille - FRANCE</span></div><div class="" style="letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; font-family: Tahoma, Verdana, Arial, Helvetica, sans-serif; font-size: 15px; line-height: 22px;"><span class="" style="font-size: 11px; font-family: verdana;">Tél: +33(0)4 13 55 18 10 - Fax : +33(0)4 13 55 18 50</span></div><div class="" style="letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none;"><font class="" face="verdana"><span class="" style="font-size: 12px; line-height: 22px;">Email : </span></font><font class="" face="verdana" color="#22bbea"><span class="" style="font-size: 12px; line-height: 20px;"><a href="mailto:pascal.boulet@univ-amu.fr" class="">pascal.boulet@univ-amu.fr</a></span></font></div><div class="" style="font-size: 14px;"><br class=""></div></div></div></div></div><br class="Apple-interchange-newline"></div><br class="Apple-interchange-newline"></div><br class="Apple-interchange-newline"></div><br class="Apple-interchange-newline"><br class="Apple-interchange-newline"></div><br class=""></div></div><span style="caret-color: rgb(0, 0, 0); font-family: Papyrus; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;" class="">_______________________________________________</span><br style="caret-color: rgb(0, 0, 0); font-family: Papyrus; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class=""><span style="caret-color: rgb(0, 0, 0); font-family: Papyrus; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;" class="">Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" class="">www.max-centre.eu/quantum-espresso</a>)</span><br style="caret-color: rgb(0, 0, 0); font-family: Papyrus; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class=""><span style="caret-color: rgb(0, 0, 0); font-family: Papyrus; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;" class="">users mailing list <a href="mailto:users@lists.quantum-espresso.org" class="">users@lists.quantum-espresso.org</a></span><br style="caret-color: rgb(0, 0, 0); font-family: Papyrus; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class=""><span style="caret-color: rgb(0, 0, 0); font-family: Papyrus; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;" class=""><a href="https://lists.quantum-espresso.org/mailman/listinfo/users" class="">https://lists.quantum-espresso.org/mailman/listinfo/users</a></span></div></blockquote></div><br class=""></div></div></body></html>