<div class="gmail_quote">On Sep 11, 2015 11:44 PM, "Peter Cybertron" <<a href="mailto:petercybertron@gmail.com">petercybertron@gmail.com</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
<br>
I am working in University of Florida as a postdoc.<br>
I am exploring the spin-orbital calculations using pwscf, and I<br>
encounter a very slow convergence like this:<br>
 > grep estimated scf.out<br>
      estimated scf accuracy    <       0.00000004 Ry<br>
      ... (100 similar lines)<br>
      estimated scf accuracy    <       0.00000003 Ry<br>
 > grep "total energy" scf.out<br>
     ...<br>
      total energy              =   -3667.91077719 Ry<br>
      total energy              =   -3667.91078176 Ry<br>
      total energy              =   -3667.91077642 Ry<br>
      total energy              =   -3667.91078300 Ry<br>
<br>
Although the convergence of total energy seems OK (down to 0.00001 Ry ~<br>
0.0001 eV = 0.1 meV), I still feel uncomfortable since I never encounter<br>
such a problem without spin-orbit interactions.<br>
<br>
I downloaded the ultrasoft and PAW pseudopotentials from PSLib and I am<br>
using the PBE flavor of XC functional. The system I am studying is a two<br>
dimensional magnetic insulator.<br>
<br>
The parameters I am using are:<br>
  constrained_magnetization = 'total direction'<br>
  fixed_magnetization(3) = 0.0<br>
  lambda = 0.2<br>
  ecutwfc = 30.0 ,<br>
  ecutrho = 300.0 ,<br>
  occupations = 'smearing' ,<br>
  degauss = 0.005,<br>
  smearing = 'mv' ,<br>
  nosym=.true.<br>
  mixing_mode = 'local-TF' ,<br>
  mixing_beta = 0.9 ,<br>
  diagonalization = 'david' ,<br>
<br>
Which parameter on earth keep my calculations non-converged?<br>
<br>
Best,<br>
Yun-Peng<br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
</blockquote></div>