<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=gb2312">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
Hello <br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
Sorry for the previous message in Italian. <br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
One of the issues could be that the pseudos have been generated with PBE and then used with the SCAN functional, this can in some cases create issues as the one you have reported. Could you provide more information about the pseudo that you have used?
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
thanks and regards <br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
Pietro <br>
</div>
<div id="appendonsend"></div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>Da:</b> users <users-bounces@lists.quantum-espresso.org> per conto di Pietro Davide Delugas <pdelugas@sissa.it><br>
<b>Inviato:</b> mercoled¨¬ 15 marzo 2023 09:38<br>
<b>A:</b> users@lists.quantum-espresso.org <users@lists.quantum-espresso.org><br>
<b>Oggetto:</b> [QE-users] R: unconverged issue in spin-polarized SCAN calculations</font>
<div> </div>
</div>
<style type="text/css" style="display:none">
<!--
p
        {margin-top:0;
        margin-bottom:0}
-->
</style>
<div dir="ltr">
<div class="x_elementToProof" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">
ok <br>
</div>
<div class="x_elementToProof" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">
gli chiediamo gli pseudi, per il resto abbiamo tutto. <br>
</div>
<div class="x_elementToProof" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">
sembrerebbe un problema collegato allo scan, alle volte se si usano pseudi PBE i sistemi tendono a esplodere.
<br>
</div>
<div id="x_appendonsend"></div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>Da:</b> users <users-bounces@lists.quantum-espresso.org> per conto di Giuseppe Mattioli <giuseppe.mattioli@ism.cnr.it><br>
<b>Inviato:</b> marted¨¬ 14 marzo 2023 12:35<br>
<b>A:</b> users@lists.quantum-espresso.org <users@lists.quantum-espresso.org><br>
<b>Oggetto:</b> Re: [QE-users] unconverged issue in spin-polarized SCAN calculations</font>
<div> </div>
</div>
<div class="x_BodyFragment"><font size="2"><span style="font-size:11pt">
<div class="x_PlainText"><br>
Dear Yi-min Ding<br>
Supposing that the pristine silicon nitride cell converges without  <br>
problems, try first to use smearing (uncomment the lines) and to  <br>
reduce mixing_beta to a low or very low value (0.1-0.01). I don't know  <br>
silicon nitride very well, so excuse me for the possibly stupid  <br>
question. Have you tried to use a larger supercell? 33% Al doping  <br>
might be not only unphysical, but also the source of wild instability  <br>
in scf iterations.<br>
HTH<br>
Giuseppe<br>
<br>
Quoting ¶¡ÒÔÃñ via users <users@lists.quantum-espresso.org>:<br>
<br>
> Dear all,<br>
> I am doing a scf calculation for a FM system using SCAN functional.  <br>
> But, the energy in output file is positive, and the run shown no  <br>
> converged tendency.<br>
><br>
><br>
> The input and out files are shown below:.<br>
> input:<br>
> &amp;system<br>
> &nbsp; &nbsp;ibrav =4<br>
> &nbsp; &nbsp;a=2.9063917063813576<br>
> &nbsp; &nbsp;c=25<br>
> &nbsp; &nbsp;nat = 7<br>
> &nbsp; &nbsp;ntyp = 3<br>
> &nbsp; &nbsp;ecutwfc = 100.0<br>
> &nbsp; &nbsp;nbnd=48<br>
> &nbsp; &nbsp;nspin=2<br>
> &nbsp; &nbsp;! starting_magnetization(1)= 0,<br>
> &nbsp; &nbsp;! starting_magnetization(2)=0,<br>
> &nbsp; &nbsp;! starting_magnetization(3)=1,<br>
> &nbsp; &nbsp; &nbsp;tot_magnetization =1<br>
> &nbsp; &nbsp;! occupations='smearing',<br>
> &nbsp; &nbsp;! smearing='gauss',<br>
> &nbsp; &nbsp;! degauss=1.0d-2,<br>
> &nbsp; &nbsp; input_dft='scan'<br>
> /<br>
> &amp;electrons<br>
> &nbsp; &nbsp;electron_maxstep = 100<br>
> &nbsp; &nbsp;conv_thr = 1.0d-8<br>
> &nbsp; &nbsp;mixing_mode = 'plain'<br>
> &nbsp; &nbsp;mixing_beta = 0.7<br>
> &nbsp; &nbsp;mixing_ndim = 8<br>
> &nbsp; &nbsp;diagonalization = 'david'<br>
> &nbsp; &nbsp;diago_david_ndim = 4<br>
> &nbsp; &nbsp;diago_full_acc = .true.<br>
> /<br>
> ATOMIC_SPECIES<br>
> &nbsp; &nbsp;Al&nbsp; 26.98&nbsp; Al.upf<br>
> &nbsp; &nbsp;Si&nbsp; 28.09&nbsp; &nbsp;Si.upf<br>
> &nbsp; &nbsp;N&nbsp; &nbsp;14.01&nbsp; &nbsp;N.upf<br>
> ATOMIC_POSITIONS (crystal)<br>
> Al 0.6666666666666643&nbsp; 0.3333333333333357&nbsp;  <br>
> 0.5000000000000000&nbsp; &nbsp;<br>
> Si&nbsp; 0.0000000000000000&nbsp; 0.0000000000000000&nbsp; 0.3860058177698406<br>
> Si&nbsp; 0.3333333429999996&nbsp; 0.6666666870000029&nbsp;  <br>
> 0.6139941822301592&nbsp;<br>
> N&nbsp; 0.0000000000000000&nbsp; 0.0000000000000000&nbsp;  <br>
> 0.4563348732468989&nbsp;<br>
> N&nbsp; 0.3333333429999996&nbsp; 0.6666666870000029&nbsp;  <br>
> 0.5436651267531011&nbsp;<br>
> N&nbsp; 0.3333333429999996&nbsp; 0.6666666870000029&nbsp;  <br>
> 0.3651769401128088&nbsp;<br>
> N 0.0000000000000000 -0.0000000000000000&nbsp; 0.6348230598871913&nbsp;<br>
> K_POINTS automatic<br>
> &nbsp;11 11 1 0 0 0<br>
><br>
> .........................................................................................................<br>
> output:<br>
> &nbsp; &nbsp; &nbsp;...............<br>
> &nbsp; &nbsp; &nbsp;iteration # 99&nbsp; &nbsp; &nbsp;ecut=&nbsp;  <br>
> &nbsp;100.00 Ry&nbsp; &nbsp; &nbsp;beta= 0.70<br>
> &nbsp; &nbsp; &nbsp;Davidson diagonalization with overlap<br>
> &nbsp; &nbsp; &nbsp;c_bands:&nbsp; 3 eigenvalues not converged<br>
> &nbsp; &nbsp; &nbsp;c_bands:&nbsp; 5 eigenvalues not converged<br>
> &nbsp; &nbsp; &nbsp;c_bands:&nbsp; 5 eigenvalues not converged<br>
> &nbsp; &nbsp; &nbsp;c_bands:&nbsp; 5 eigenvalues not converged<br>
> &nbsp; &nbsp; &nbsp;c_bands:&nbsp; 5 eigenvalues not converged<br>
> &nbsp; &nbsp; &nbsp;c_bands:&nbsp; 5 eigenvalues not converged<br>
> &nbsp; &nbsp; &nbsp;ethr =&nbsp; 1.00E-02,&nbsp; avg # of iterations = 24.8<br>
><br>
><br>
> &nbsp; &nbsp; &nbsp;negative rho (up,down):&nbsp; 5.052E+00 2.651E+00<br>
><br>
><br>
> &nbsp; &nbsp; &nbsp;total cpu time spent up to now is&nbsp; &nbsp;  <br>
> &nbsp;1964.0 secs<br>
><br>
><br>
> &nbsp; &nbsp; &nbsp;total energy&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  <br>
> &nbsp; &nbsp; =&nbsp; &nbsp; 3309.58902363 Ry<br>
> &nbsp; &nbsp; &nbsp;estimated scf accuracy&nbsp; &nbsp; <&nbsp;  <br>
> 304153.85545920 Ry<br>
><br>
><br>
> &nbsp; &nbsp; &nbsp;total magnetization&nbsp; &nbsp; &nbsp;  <br>
> &nbsp;=&nbsp; &nbsp; &nbsp;1.00 Bohr mag/cell<br>
> &nbsp; &nbsp; &nbsp;absolute magnetization&nbsp; &nbsp; =&nbsp;  <br>
> &nbsp; 28.84 Bohr mag/cell<br>
><br>
><br>
> &nbsp; &nbsp; &nbsp;iteration #100&nbsp; &nbsp; &nbsp;ecut=&nbsp;  <br>
> &nbsp;100.00 Ry&nbsp; &nbsp; &nbsp;beta= 0.70<br>
> &nbsp; &nbsp; &nbsp;Davidson diagonalization with overlap<br>
> &nbsp; &nbsp; &nbsp;ethr =&nbsp; 1.00E-02,&nbsp; avg # of iterations = 12.3<br>
><br>
><br>
> &nbsp; &nbsp; &nbsp;negative rho (up,down):&nbsp; 3.516E+00 2.027E+00<br>
><br>
><br>
> &nbsp; &nbsp; &nbsp;total cpu time spent up to now is&nbsp; &nbsp;  <br>
> &nbsp;1974.9 secs<br>
><br>
><br>
> &nbsp; &nbsp; &nbsp;total energy&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  <br>
> &nbsp; &nbsp; =&nbsp; &nbsp; 2760.95716210 Ry<br>
> &nbsp; &nbsp; &nbsp;estimated scf accuracy&nbsp; &nbsp; <&nbsp;  <br>
> 488812.03105752 Ry<br>
><br>
><br>
> &nbsp; &nbsp; &nbsp;total magnetization&nbsp; &nbsp; &nbsp;  <br>
> &nbsp;=&nbsp; &nbsp; &nbsp;1.00 Bohr mag/cell<br>
> &nbsp; &nbsp; &nbsp;absolute magnetization&nbsp; &nbsp; =&nbsp;  <br>
> &nbsp; 30.68 Bohr mag/cell<br>
><br>
><br>
> &nbsp; &nbsp; &nbsp;End of self-consistent calculation<br>
><br>
><br>
> &nbsp; &nbsp; &nbsp;convergence NOT achieved after 100 iterations: stopping<br>
> &nbsp; &nbsp; &nbsp;Writing meta-gga kinetic term<br>
><br>
> ...............................................................................................................<br>
> Can you give me some help?<br>
> Thanks in advance.<br>
> Dr. Yi-min Ding<br>
<br>
<br>
<br>
GIUSEPPE MATTIOLI<br>
CNR - ISTITUTO DI STRUTTURA DELLA MATERIA<br>
Via Salaria Km 29,300 - C.P. 10<br>
I-00015 - Monterotondo Scalo (RM)<br>
Mob (*preferred*) +39 373 7305625<br>
Tel + 39 06 90672342 - Fax +39 06 90672316<br>
E-mail: <giuseppe.mattioli@ism.cnr.it><br>
<br>
_______________________________________________<br>
The Quantum ESPRESSO community stands by the Ukrainian<br>
people and expresses its concerns about the devastating<br>
effects that the Russian military offensive has on their<br>
country and on the free and peaceful scientific, cultural,<br>
and economic cooperation amongst peoples<br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu">www.max-centre.eu</a>)<br>
users mailing list users@lists.quantum-espresso.org<br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users">https://lists.quantum-espresso.org/mailman/listinfo/users</a></div>
</span></font></div>
</div>
</body>
</html>