<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">
ok <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">
gli chiediamo gli pseudi, per il resto abbiamo tutto. <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">
sembrerebbe un problema collegato allo scan, alle volte se si usano pseudi PBE i sistemi tendono a esplodere.
<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 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="BodyFragment"><font size="2"><span style="font-size:11pt;">
<div class="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>
</body>
</html>