<meta http-equiv="Content-Type" content="text/html; charset=UTF-8"><div>Dear all,</div><div>I am doing a scf calculation for a FM system using SCAN functional. But, the energy in output file is positive, and the run shown no converged tendency. </div><div><br></div><div>The input and out files are shown below:. </div><div>input:</div><div><div>&system</div><div>   ibrav =4</div><div>   a=2.9063917063813576</div><div>   c=25 </div><div>   nat = 7 </div><div>   ntyp = 3</div><div>   ecutwfc = 100.0</div><div>   nbnd=48</div><div>   nspin=2</div><div>   ! starting_magnetization(1)= 0,</div><div>   ! starting_magnetization(2)=0,</div><div>   ! starting_magnetization(3)=1,</div><div>     tot_magnetization =1</div><div>   ! occupations='smearing',</div><div>   ! smearing='gauss',</div><div>   ! degauss=1.0d-2,</div><div>    input_dft='scan'</div><div>/</div><div>&electrons</div><div>   electron_maxstep = 100</div><div>   conv_thr = 1.0d-8</div><div>   mixing_mode = 'plain'</div><div>   mixing_beta = 0.7</div><div>   mixing_ndim = 8</div><div>   diagonalization = 'david'</div><div>   diago_david_ndim = 4</div><div>   diago_full_acc = .true.</div><div>/</div><div>ATOMIC_SPECIES</div><div>   Al  26.98  Al.upf</div><div>   Si  28.09   Si.upf</div><div>   N   14.01   N.upf</div><div>ATOMIC_POSITIONS (crystal)</div><div>Al 0.6666666666666643  0.3333333333333357  0.5000000000000000   </div><div>Si  0.0000000000000000  0.0000000000000000  0.3860058177698406 </div><div>Si  0.3333333429999996  0.6666666870000029  0.6139941822301592  </div><div>N  0.0000000000000000  0.0000000000000000  0.4563348732468989  </div><div>N  0.3333333429999996  0.6666666870000029  0.5436651267531011  </div><div>N  0.3333333429999996  0.6666666870000029  0.3651769401128088  </div><div>N 0.0000000000000000 -0.0000000000000000  0.6348230598871913  </div><div>K_POINTS automatic</div><div> 11 11 1 0 0 0</div></div><div>.........................................................................................................</div><div>output:</div><div><div>     ...............</div><div>     iteration # 99     ecut=   100.00 Ry     beta= 0.70</div><div>     Davidson diagonalization with overlap</div><div>     c_bands:  3 eigenvalues not converged</div><div>     c_bands:  5 eigenvalues not converged</div><div>     c_bands:  5 eigenvalues not converged</div><div>     c_bands:  5 eigenvalues not converged</div><div>     c_bands:  5 eigenvalues not converged</div><div>     c_bands:  5 eigenvalues not converged</div><div>     ethr =  1.00E-02,  avg # of iterations = 24.8</div><div><br></div><div>     negative rho (up,down):  5.052E+00 2.651E+00</div><div><br></div><div>     total cpu time spent up to now is     1964.0 secs</div><div><br></div><div>     total energy              =    3309.58902363 Ry</div><div>     estimated scf accuracy    <  304153.85545920 Ry</div><div><br></div><div>     total magnetization       =     1.00 Bohr mag/cell</div><div>     absolute magnetization    =    28.84 Bohr mag/cell</div><div><br></div><div>     iteration #100     ecut=   100.00 Ry     beta= 0.70</div><div>     Davidson diagonalization with overlap</div><div>     ethr =  1.00E-02,  avg # of iterations = 12.3</div><div><br></div><div>     negative rho (up,down):  3.516E+00 2.027E+00</div><div><br></div><div>     total cpu time spent up to now is     1974.9 secs</div><div><br></div><div>     total energy              =    2760.95716210 Ry</div><div>     estimated scf accuracy    <  488812.03105752 Ry</div><div><br></div><div>     total magnetization       =     1.00 Bohr mag/cell</div><div>     absolute magnetization    =    30.68 Bohr mag/cell</div><div><br></div><div>     End of self-consistent calculation</div><div><br></div><div>     convergence NOT achieved after 100 iterations: stopping</div><div>     Writing meta-gga kinetic term</div></div><div>...............................................................................................................</div><div>Can you give me some help?</div><div>Thanks in advance.</div><div>Dr. Yi-min Ding</div>