<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt"><div><span></span></div><div style="font-family: times new roman, new york, times, serif; font-size: 12pt;"><div style="font-family: times new roman, new york, times, serif; font-size: 12pt;"><div id="yiv687496978"><div><div style="color:#000;background-color:#fff;font-family:times new roman, new york, times, serif;font-size:12pt;"><div>Dear wien2k users,</div><div> I am trying to do a LDA+U calculation on a Fe2O3 with 50 atoms per unit cell. My struct file shows that all 50 atoms are non-equivalent.I had run <span class="yiv687496978yshortcuts yiv687496978cs4-visible" id="yiv687496978lw_1332896168_0">SCF</span> cycles separately&nbsp; with only spin polarized option and then take into account&nbsp; LDA+U only for Fe atoms which are 20 atoms (index in struct file&nbsp; 31 -50).</div><div><br></div><div>i prepared my case.inorb

 and case.indmc like this</div><div style="color:rgb(64, 127, 0);"><span style="font-weight:bold;text-decoration:underline;">case.indmc</span><br></div><div style="color:rgb(64, 127, 0);">-9.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Emin cutoff

 energy<br>&nbsp;20&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; number of atoms for which density matrix is calculated<br>&nbsp;31&nbsp; 1&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp; index of 1st atom, number of L's, L1<br>&nbsp;32&nbsp; 1&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp; dtto for 2nd atom, repeat NATOM times&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;33&nbsp; 1&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;34&nbsp; 1&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;35&nbsp; 1&nbsp;

 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;36&nbsp; 1&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;37&nbsp; 1&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;38&nbsp; 1&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;39&nbsp; 1&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;40&nbsp; 1&nbsp;

 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;41&nbsp; 1&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;42&nbsp; 1&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;43&nbsp; 1&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;44&nbsp; 1&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;45&nbsp; 1&nbsp;

 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;46&nbsp; 1&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;47&nbsp; 1&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;48&nbsp; 1&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;49&nbsp; 1&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;50&nbsp; 1&nbsp;

 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;0 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; r-index, (l,s)index&nbsp;&nbsp;</div><div style="color:rgb(64, 127, 0);"><br></div><div style="color:rgb(64, 127, 0);"><span style="font-weight:bold;text-decoration:underline;">case .inorb</span></div><div style="color:rgb(64, 127, 0);"><br></div><div><span style="color:rgb(64, 127, 0);">&nbsp; 1 20&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nmod, natorb, ipr</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">PRATT&nbsp; 1.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; BROYD/PRATT, mixing</span><br style="color:rgb(64, 127,

 0);"><span style="color:rgb(64, 127, 0);">&nbsp; 31 1 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iatom nlorb, lorb</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp; 32 1 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iatom nlorb, lorb</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp; 33 1 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iatom nlorb, lorb</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp; 34 1

 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iatom nlorb, lorb</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp; 35 1 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iatom nlorb, lorb</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp; 36 1 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iatom nlorb, lorb</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp; 37 1 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iatom nlorb, lorb</span><br

 style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp; 38 1 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iatom nlorb, lorb</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp; 39 1 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iatom nlorb, lorb</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp; 40 1 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iatom nlorb, lorb</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp; 41 1

 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iatom nlorb, lorb</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp; 42 1 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iatom nlorb, lorb</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp; 43 1 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iatom nlorb, lorb</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp; 44 1 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iatom nlorb, lorb</span><br

 style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp; 45 1 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iatom nlorb, lorb</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp; 46 1 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iatom nlorb, lorb</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp; 47 1 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iatom nlorb, lorb</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp; 48 1

 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iatom nlorb, lorb</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp; 49 1 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iatom nlorb, lorb</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp; 50 1 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iatom nlorb, lorb</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nsic

 0..AFM, 1..SIC, 2..HFM</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp; 0.29 0.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; U J (Ry)&nbsp;&nbsp; Note: we recommend to use U_eff = U-J and J=0</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp; 0.29 0.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; U J</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp; 0.29 0.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; U J</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp; 0.29 0.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; U J</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp; 0.29 0.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; U J</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp; 0.29

 0.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; U J</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp; 0.29 0.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; U J</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp; 0.29 0.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; U J</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp; 0.29 0.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; U J</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp; 0.29 0.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; U J</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp; 0.29 0.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; U J</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp; 0.29 0.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; U J</span><br style="color:rgb(64,

 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp; 0.29 0.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; U J</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp; 0.29 0.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; U J</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp; 0.29 0.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; U J</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp; 0.29 0.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; U J</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp; 0.29 0.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; U J</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp; 0.29 0.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; U J</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127,

 0);">&nbsp;&nbsp; 0.29 0.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; U J</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp; 0.29 0.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; U J</span></div><div><br></div><div>but when i run scf cycle with LDA+U it stops and when i checked case.outputorbup file it gives this error.</div><div><br></div><div>&nbsp;C<span style="color:rgb(64, 127, 0);">alculation of orbital potential for spin block: up&nbsp; </span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;Type of potential:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; LDA+U&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;Vorb applied to atom&nbsp; 31 orbit.

 numbers&nbsp;&nbsp; 2</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;Vorb applied to atom&nbsp; 32 orbit. numbers&nbsp;&nbsp; 2</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;Vorb applied to atom&nbsp; 33 orbit. numbers&nbsp;&nbsp; 2</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;Vorb applied to atom&nbsp; 34 orbit. numbers&nbsp;&nbsp; 2</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;Vorb applied to atom&nbsp; 35 orbit. numbers&nbsp;&nbsp; 2</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;Vorb applied to atom&nbsp; 36 orbit. numbers&nbsp;&nbsp; 2</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;Vorb applied to atom&nbsp; 37 orbit. numbers&nbsp;&nbsp; 2</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127,

 0);">&nbsp;Vorb applied to atom&nbsp; 38 orbit. numbers&nbsp;&nbsp; 2</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;Vorb applied to atom&nbsp; 39 orbit. numbers&nbsp;&nbsp; 2</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;Vorb applied to atom&nbsp; 40 orbit. numbers&nbsp;&nbsp; 2</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;Vorb applied to atom&nbsp; 41 orbit. numbers&nbsp;&nbsp; 2</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;Vorb applied to atom&nbsp; 42 orbit. numbers&nbsp;&nbsp; 2</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;Vorb applied to atom&nbsp; 43 orbit. numbers&nbsp;&nbsp; 2</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;Vorb applied to atom&nbsp; 44 orbit. numbers&nbsp;&nbsp; 2</span><br style="

color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;Vorb applied to atom&nbsp; 45 orbit. numbers&nbsp;&nbsp; 2</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;Vorb applied to atom&nbsp; 46 orbit. numbers&nbsp;&nbsp; 2</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;Vorb applied to atom&nbsp; 47 orbit. numbers&nbsp;&nbsp; 2</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;Vorb applied to atom&nbsp; 48 orbit. numbers&nbsp;&nbsp; 2</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;Vorb applied to atom&nbsp; 49 orbit. numbers&nbsp;&nbsp; 2</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;Vorb applied to atom&nbsp; 50 orbit. numbers&nbsp;&nbsp; 2</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp; Around the mean field

 method</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Atom 31 L=&nbsp; 2 U=&nbsp; 0.290 J=&nbsp; 0.000 Ry</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Atom 32 L=&nbsp; 2 U=&nbsp; 0.290 J=&nbsp; 0.000 Ry</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Atom 33 L=&nbsp; 2 U=&nbsp; 0.290 J=&nbsp; 0.000 Ry</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Atom 34 L=&nbsp; 2 U=&nbsp; 0.290 J=&nbsp; 0.000 Ry</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Atom 35 L=&nbsp; 2 U=&nbsp; 0.290 J=&nbsp; 0.000 Ry</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127,

 0);">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Atom 36 L=&nbsp; 2 U=&nbsp; 0.290 J=&nbsp; 0.000 Ry</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Atom 37 L=&nbsp; 2 U=&nbsp; 0.290 J=&nbsp; 0.000 Ry</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Atom 38 L=&nbsp; 2 U=&nbsp; 0.290 J=&nbsp; 0.000 Ry</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Atom 39 L=&nbsp; 2 U=&nbsp; 0.290 J=&nbsp; 0.000 Ry</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Atom 40 L=&nbsp; 2 U=&nbsp; 0.290 J=&nbsp; 0.000 Ry</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Atom 41 L=&nbsp; 2 U=&nbsp; 0.290 J=&nbsp;

 0.000 Ry</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Atom 42 L=&nbsp; 2 U=&nbsp; 0.290 J=&nbsp; 0.000 Ry</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Atom 43 L=&nbsp; 2 U=&nbsp; 0.290 J=&nbsp; 0.000 Ry</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Atom 44 L=&nbsp; 2 U=&nbsp; 0.290 J=&nbsp; 0.000 Ry</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Atom 45 L=&nbsp; 2 U=&nbsp; 0.290 J=&nbsp; 0.000 Ry</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Atom 46 L=&nbsp; 2 U=&nbsp; 0.290 J=&nbsp; 0.000 Ry</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127,

 0);">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Atom 47 L=&nbsp; 2 U=&nbsp; 0.290 J=&nbsp; 0.000 Ry</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Atom 48 L=&nbsp; 2 U=&nbsp; 0.290 J=&nbsp; 0.000 Ry</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Atom 49 L=&nbsp; 2 U=&nbsp; 0.290 J=&nbsp; 0.000 Ry</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Atom 50 L=&nbsp; 2 U=&nbsp; 0.290 J=&nbsp; 0.000 Ry</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp; end of OP input</span><br style="color:rgb(64, 127, 0);"><span style="color:rgb(64, 127, 0);">&nbsp;<span style="font-weight: bold;">STRUCT file read</span></span><br style="color: rgb(64, 127, 0); font-weight: bold;"><span style="color: rgb(64, 127, 0);

 font-weight: bold;">&nbsp; VSP read</span><br style="color: rgb(64, 127, 0); font-weight: bold;"><span style="color: rgb(64, 127, 0); font-weight: bold;">&nbsp; Conflict in atom indexes: iatom&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 31 ne jatom&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1</span><br></div><div><br></div><div><br></div><div>&nbsp;i don't know how to solve this,&nbsp; please guide me about this. i am really unable to solve this.</div><div><br></div><div>i will be very grateful to you.</div><div><br></div><div>Best Regards<br><br>Arqum Hashmi<br><br><br><br></div></div></div></div> </div> </div>  </div></body></html>