<div dir="ltr"><div><div>Dear QE Users,</div><div><br></div><div>I am running a AFM case, QE.in is undermentioned, aiming to get the atomic magnetic moment for Cr. In literature, it is reported ~1.5 mue_b for each Cr.</div><div><br></div><div></div><div>I tried with many configurations for starting magnetization and for each configuration I am getting total (MMT)  and well atomic magnetic (MMI) moment zero.</div><div><br></div><div>I am only getting the magnetic moment in the vicinity of reported value if I take ecutwfc 20 and ecutrho 200 which is below and in tune with literature. <br></div><div class="gmail_extra"><br></div><div class="gmail_extra">I took</div><div class="gmail_extra">  nspin=2,<br>  starting_magnetization(3)=1,<br>  starting_magnetization(4)=-1,</div><div class="gmail_extra"><br clear="all"></div><div class="gmail_extra"><div><div class="gmail-m_152708887253533859gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><span style="font-size:12.8px">     Magnetic moment per site:<br>     atom:    1    charge:    2.1886    magn:   -0.0429    constr:    0.0000<br>     atom:    2    charge:   10.6671    magn:    0.0204    constr:    0.0000<br>     atom:    3    charge:   12.3448    magn:   -1.7550    constr:    0.0000<br>     atom:    4    charge:   12.3094    magn:    1.8137    constr:    0.0000<br></span></div><div dir="ltr"><span style="font-size:12.8px"><br></span></div><div dir="ltr"><span style="font-size:12.8px"><br></span></div>for rest, all configurations, the MMT (total) or MMI (atom dependent) is always zero.<br></div><div dir="ltr"><br></div><div dir="ltr"><br></div>other configurations, that I used.<br><div><br></div><div></div><div>  nspin=2,<br>  starting_magnetization(3)=2,<br>  starting_magnetization(4)=-2,</div><div></div><div>or</div><div></div><div></div><div>  nspin=2,<br> starting_magnetization(3)=1,<br></div><div></div><div></div><div>or</div><div></div><div>  nspin=2,<br>  starting_magnetization=05d,  <br></div></div></div></div></div></div></div></div></div><div class="gmail_extra"><br></div><div class="gmail_extra"></div><div class="gmail_extra"><b>for a range of ecutwfc and ecutrho</b></div></div><div><br></div><div>Could you please correct me if I am doing any mistake somewhere? I am using rrkjus PPs.<br></div><div><br></div><div>#================ QE.in ==================<br></div><div><br></div><div>&CONTROL</div>  calculation='scf',<br>  outdir='.',<br>  prefix='pwscf',<br>  pseudo_dir='.',<br>  verbosity='low',<br>/<br><br>&SYSTEM<br>  ibrav=2,<br>  celldm(1)=11.0560959151d0,<br>  nat=4,<br>  ntyp=3,<br>  ecutwfc=$k,<br>  ecutrho=$p,<br>  input_dft='pbe',<br>  occupations='smearing',<br>  smearing='mv',<br>  degauss=0.005d0,<br>/<br><br>&ELECTRONS<br>  conv_thr=1d-08,<br>  mixing_beta=0.7d0,<br>/<br><br>ATOMIC_SPECIES<br>  Cr 51.996100d0 Cr.UPF<br>  Si 28.085500d0 Si.UPF<br>  Zn 65.409000d0 Zn.UPF<br><br>ATOMIC_POSITIONS {crystal}<br>  Si   0.0000000000d0   0.0000000000d0   0.0000000000d0<br>  Zn   0.7500000000d0   0.7500000000d0   0.7500000000d0<br>  Cr   0.2500000000d0   0.2500000000d0   0.2500000000d0<br>  Cr   0.5000000000d0   0.5000000000d0   0.5000000000d0<br><br>K_POINTS {automatic}<br>  12 12 12 0 0 0</div>