<div dir="ltr"><div class="gmail_default" style="font-size:large;color:#073763">Dear all, </div><div class="gmail_default" style="font-size:large;color:#073763">I am trying to do 'SCF' calculations for CrI3 by PAW-PBE pseudopotentials.  (I used full-relativistic pseudopotential from QE site). I need to consider SOC effect in my calculations. When I use starting_magnetization for Cr atoms, the SCF is not converged after 400 iterations!! I will copy my input file and the last lines of my output at the end. Is there any suggestion to overcome this? </div><div class="gmail_default" style="font-size:large;color:#073763">Please noticed that this calculation was converged after 12 iterations without starting magnetization. However, I know that Cr is a magnetic atom and I should do the calculations with starting magnetization.</div><div class="gmail_default" style="font-size:large;color:#073763"><br></div><div class="gmail_default" style="font-size:large;color:#073763">Best regards, </div><div class="gmail_default" style="font-size:large;color:#073763">Sahar</div><div class="gmail_default" style="font-size:large;color:#073763">Postdoc</div><div class="gmail_default" style="font-size:large;color:#073763">IPM</div><div class="gmail_default" style="font-size:large;color:#073763"><br></div><div class="gmail_default" style="font-size:large;color:#073763">!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!</div><div class="gmail_default" style="font-size:large;color:#073763"><b>My input:</b></div><div class="gmail_default" style="font-size:large;color:#073763"><b>!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!</b></div><div class="gmail_default" style="font-size:large;color:#073763">&CONTROL<br>  calculation='scf',<br>  outdir='./',<br>  prefix='CrI3SOC',<br>  pseudo_dir='../',<br>/<br>&SYSTEM<br>  ibrav=4,<br>  celldm(1)=13.2422747152d0, celldm(3)=3.5676010452d0,<br>  nat=8,<br>  ntyp=3,<br>  ecutwfc=50,<br>  ecutrho=500,<br>  occupations='smearing',<br>  smearing='mv',<br>  degauss=0.01d0,<br>  lspinorb=.true.<br>  noncolin=.true.<br>  starting_magnetization(1)=0.3,<br>  starting_magnetization(2)=0.3,<br>  angle1(1)=0,<br>  angle1(2)=0,<br> /<br><br>&ELECTRONS<br>  conv_thr=1d-06,<br>  electron_maxstep=400,<br>  mixing_beta=0.3d0,<br> /<br><br>ATOMIC_SPECIES<br>  Cr1 51.996100d0 Cr.rel-pbe-spn-kjpaw_psl.1.0.0.UPF<br>  Cr2 51.996100d0 Cr.rel-pbe-spn-kjpaw_psl.1.0.0.UPF<br>  I 126.904000d0 I.rel-pbe-n-kjpaw_psl.0.2.2.UPF<br><br>ATOMIC_POSITIONS {crystal}<br>  Cr1   0.7443381782d0   0.2568556861d0   0.5000002867d0<br>  Cr2   0.4101686019d0   0.5897708980d0   0.5000002867d0<br>   I   0.0769393899d0   0.5644218819d0   0.4372216127d0<br>   I   0.0769393899d0   0.2818907021d0   0.5627783787d0<br>   I   0.4358136183d0   0.2821277022d0   0.4370458976d0<br>   I   0.4358136183d0   0.9230591101d0   0.5629541024d0<br>   I   0.7179617978d0   0.5641268817d0   0.5627224713d0<br>   I   0.7179617978d0   0.9232081102d0   0.4372775287d0<br><br>K_POINTS {automatic}<br>  4 4 1 0 0 0<br></div><div class="gmail_default" style="font-size:large;color:#073763"><br></div><div class="gmail_default" style="font-size:large;color:#073763">!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!</div><div class="gmail_default" style="font-size:large;color:#073763"><b>the last lines of my output</b> <br></div><div class="gmail_default" style="font-size:large;color:#073763">!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!</div><div class="gmail_default" style="font-size:large;color:#073763">total cpu time spent up to now is    43505.0 secs<br><br>     total energy              =   -2784.23985010 Ry<br>     Harris-Foulkes estimate   =   -2784.14428568 Ry<br>     estimated scf accuracy    <      45.17285183 Ry<br><br>     total magnetization       =    -0.44     0.19    -1.95 Bohr mag/cell<br>     absolute magnetization    =    56.91 Bohr mag/cell<br><br>     End of self-consistent calculation<br><br>     convergence NOT achieved after 400 iterations: stopping<br><br></div><div class="gmail_default" style="font-size:large;color:#073763">!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!</div></div>