Dear Stefan,<br /><br />Thank you for the explanations.<br /><br />I actually meant de-localized states instead of localized. I had already tried an Ueff of 8.5 ev, the result is that the split of the f&nbsp; states was more important and in the same manner as you had described.<br /><br /><br />Regarding the two different values of&nbsp; Ef, I checked&nbsp; case.inorb and it seems ok. U is applied only on the&nbsp; f states of Nd:<br /><br />&nbsp; 1&nbsp; 2&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<br />PRATT&nbsp; 1.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; BROYD/PRATT, mixing<br />&nbsp; 1 1 3&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<br />&nbsp; 2 1 3&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<br />&nbsp;&nbsp; 1&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<br />&nbsp;&nbsp; 0.625 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<br />&nbsp;&nbsp; 0.625 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<br />&nbsp;<br />the calculation are performed as follow :<br /><br />runsp_lapw -p -fc 0.1 <br />save_lapw case1<br />runsp_lapw -orbc -p -fc 0.1<br />save_lapw case2<br />runsp_lapw -orb -p -fc 0.1. <br /><br />x lapw2 -qtl -c -up&nbsp; -p<br />x lapw2 -qtl -c -dn&nbsp; -p<br /><br />Is it correct ?<br /><br /><br />Thank you for your help,<br />Laurent Baudery<br /><br />&gt;&gt; I performed spin-polarized calculations on NdTiO3 compound and as excepted the <span class="yshortcuts" style="border-bottom: 1px dashed #0066cc; cursor: pointer">fermi level</span>&nbsp; lies at the 4f states (due to the strong localisation).<br /><br />&gt;(Probably a word is missing in your last sentence, which gives it the opposite meaning? ==&gt; de-localization) The reason why LDA puts the 4f states at the fermi level is that it cannot sufficiently localize these states. They are hence (incorrectly) described as itinerant states.<br /><br />&gt;&gt; Thereafter, I performed spin-polarized LDA+U (Ueff=6.5 ev) calculations <br /><br />&gt;This value of U might be somewhat at the lower side for lanthanides. Try 8.5 eV (see <a href="http://link.aps.org/doi/10.1103/PhysRevB.74.014409" target="_blank"><span class="yshortcuts">http://link.aps.org/doi/10.1103/PhysRevB.74.014409</span></a> and <a href="http://link.aps.org/doi/10.1103/PhysRevB.77.155101" target="_blank"><span class="yshortcuts">http://link.aps.org/doi/10.1103/PhysRevB.77.155101</span></a> , and references therein).<br /><br />&gt;&gt; which split the 4f (up) states : the occupied states down by 3 ev below the Ef and empty states rises by about 5 ev above the Ef. I got two different values of Ef(up) and Ef(dn).<br />&gt; &gt;<br />&gt; &gt;By applying Ueff=6.5 ev I expected a little split that will let the 4f occupied states at the top of the <span class="yshortcuts" style="border-bottom: 1px dashed #0066cc; cursor: pointer">valence band</span> (around the Ef).&nbsp; &nbsp; Can someone give me a comment about this ? Are also normal to have two different values of Ef (up and dn)?<br /><br />&gt;Such two different values of Ef I've never seen. Are you sure you aren't by mistake making an applied field calculation (which is another input number in the same case.inorb file you use for LDA+U).<br /><br />&gt;What should happen for Nd is that all dn-states (which are unoccupied) move up and find themselves all above Ef, while the up-states split in an occupied part several eV below Ef and an unoccupied part several eV above Ef.<br /><br />Stefaan<BR><BR><table><tbody><tr><td valign="middle"><img style="margin: 0pt;" src="http://webmail.laposte.net/webmail/fr_FR/panels/images/logo10ans.jpg" alt="Logo10ans" width="122" height="150"></td><td valign="middle"><p style="border-top: 1px dashed rgb(215, 0, 103); border-bottom: 1px dashed rgb(215, 0, 103); padding: 20px 0pt;">Gratuite, garantie à vie et déjà utilisée par des millions d'internautes... <br>vous aussi, pour votre adresse e-mail, choisissez <a style="color: rgb(215, 0, 103); text-decoration: underline;" target="_blank" href="https://compte.laposte.net/inscription/index.do?jeux=10ANSFOOTER">laposte.net</a>. <br><br><span style="font-family: Time New Roman,serif; font-size: 16px; color: rgb(215, 0, 103);"><b>Laposte.net</b>, bien + qu'une messagerie</span></p></td></tr></tbody></table></body>