<div dir="ltr"><div dir="ltr"><div>Dear Wien community,</div><div><br></div><div>I am performing LDA+U
 type calculation for SmCo5, and interested in look in the role of Co in
 generating large magnetic anisotropy for this materials.</div><div> </div><div>(i)  First i ran a LDA+U calculation with U of 5 eV on Sm, with the following case.inorb file. Everything worked fine.</div><div><br></div><div>1  1  0                         nmod, natorb, ipr<br>PRATT  1.0                    BROYD/PRATT, mixing<br>1 1 3                          iatom nlorb, lorb<br>1                              nsic 0..AFM, 1..SIC, 2..HFM<br>0.367  0.00        U J (Ry)   Note: we recommend to use U_eff = U-J and J=0</div><div><br></div><div>(ii) then I included U of 3 eV at Co sites, in the following case.inorb  input file. Everything worked fine.</div><div>1  6  0                         nmod, natorb, ipr<br>PRATT  1.0                    BROYD/PRATT, mixing<br>1 1 3                          iatom nlorb, lorb</div><div>2 1 2                          iatom nlorb, lorb</div><div>3 1 2                          iatom nlorb, lorb</div><div>4 1 2                          iatom nlorb, lorb</div><div>5 1 2                          iatom nlorb, lorb</div><div>6 1 2                          iatom nlorb, lorb</div><div>1                              nsic 0..AFM, 1..SIC, 2..HFM<br>0.367  0.00        U J (Ry)   Note: we recommend to use U_eff = U-J and J=0</div><div>0.220  0.00        U J (Ry)   Note: we recommend to use U_eff = U-J and J=0</div><div>0.220  0.00        U J (Ry)   Note: we recommend to use U_eff = U-J and J=0</div><div>0.220  0.00        U J (Ry)   Note: we recommend to use U_eff = U-J and J=0</div><div>0.220  0.00        U J (Ry)   Note: we recommend to use U_eff = U-J and J=0</div><div>0.220  0.00        U J (Ry)   Note: we recommend to use U_eff = U-J and J=0</div><div><br></div><div>(iii) At last I wanted to do a calculation with U just at Co site. I use the following case.inorb  input file:</div><div><br></div><div><div>1  5  0                         nmod, natorb, ipr<br>PRATT  1.0                    BROYD/PRATT, mixing<br></div><div>2 1 2                          iatom nlorb, lorb</div><div>3 1 2                          iatom nlorb, lorb</div><div>4 1 2                          iatom nlorb, lorb</div><div>5 1 2                          iatom nlorb, lorb</div><div>6 1 2                          iatom nlorb, lorb</div><div>1                              nsic 0..AFM, 1..SIC, 2..HFM<br><div>0.220  0.00        U J (Ry)   Note: we recommend to use U_eff = U-J and J=0</div><div>0.220  0.00        U J (Ry)   Note: we recommend to use U_eff = U-J and J=0</div><div>0.220  0.00        U J (Ry)   Note: we recommend to use U_eff = U-J and J=0</div><div>0.220  0.00        U J (Ry)   Note: we recommend to use U_eff = U-J and J=0</div><div>0.220  0.00        U J (Ry)   Note: we recommend to use U_eff = U-J and J=0</div><div><br></div><div>This run HOWEVER gives an error in case.outputorbup/dn file</div><div><br></div><div><br></div><div><br></div><div> Conflict in atom indexes: iatom           2 ne jatom           1</div></div><div><br></div><div><br></div><div><br></div><div>This will works if I put U = 0eV at Sm in step (2). But before proceeding I wanted to clarify the following</div><div><br></div><div>(I) Is this expected behavior ?</div><div><br></div><div>(II) If I put a U of 0 eV at Sm, will that mimic a calculation with U only at Co site?</div><div><br></div><div><br></div><div>I
 can make this work by swapping the atom (putting the Co first) in the 
structure files, but then I will have to restart the entire run.</div><div><br></div><div>Your suggestions and thoughts are much appreciated.</div><div><br></div><div>Thanks in advance.</div><div><br></div><div><br></div><div>Regards,</div><div>David Yang</div></div></div></div>