<p style="margin:0;padding:0;" id="result_box" dir="ltr">Dear stephane, <br /><br />Once again thank you for your help. <br /><br />I checked my files and actually before lunching runsp_lapw-orbc the case.dmat is empty. I thought the script runsp_lapw execute automatically lapwdm due to the presence of the inorbc&nbsp; option and the appropriate case.indmc file. I will recalculate by adding lapwdm before executing&nbsp; runsp_lapw -orbc.</p><br />For the Ef, I do simply a grep command and&nbsp; I get:<br /><br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.60908<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.60908<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.54087<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.54275<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.53488<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.54824<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.54019<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.54627<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.53275<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.54530<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.49869<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.53494<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.48346<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.68766<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.48440<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.66995<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.48226<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.65435<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.48261<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.66488<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.48292<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.65511<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.49169<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.67817<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.48419<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.67826<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.48373<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.68742<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.48172<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.68709<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.48223<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.68729<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.48286<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.67837<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.48289<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.68727<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.48254<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.68725<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.48251<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.68725<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.48270<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.68720<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.48279<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.68726<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.48245<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.68731<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.48230<br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.68735<br />&nbsp;<br />When I look in the&nbsp; case.scf2up/dn files, I get also: <br /><br /> case.scf2up&nbsp; :<br /><br />&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; GAUSS-SMEARING WITH&nbsp;&nbsp;&nbsp; 0.00100 Ry<br /><br />&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; BANDENERGY CORRECTION:&nbsp; 0.000000<br />&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Bandranges (emin - emax) and occupancy:<br />&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Energy to separate low and high energystates:&nbsp;&nbsp;&nbsp; 0.11380<br /><br /><br />:NOE&nbsp; : NUMBER OF ELECTRONS&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 376.000<br /><br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.48230<br /><br />case.scf2dn&nbsp; :<br /><br />&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; GAUSS-SMEARING WITH&nbsp;&nbsp;&nbsp; 0.00100 Ry<br /><br />&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; BANDENERGY CORRECTION:&nbsp; 0.000000<br />&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Bandranges (emin - emax) and occupancy:<br />&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Energy to separate low and high energystates:&nbsp;&nbsp;&nbsp; 0.12854<br /><br /><br />:NOE&nbsp; : NUMBER OF ELECTRONS&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 376.000<br /><br />:FER&nbsp; : F E R M I - ENERGY(GAUSS-.M.)=&nbsp;&nbsp; 0.68735<br /><br /><br />Laurent<br /><br /><br /><br />&gt; 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 />Yes, that file seems to be fine.<br /><br />&gt;&nbsp;  1&nbsp; 2&nbsp; 0&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  nmod, natorb, ipr<br />&gt; PRATT&nbsp; 1.0&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; BROYD/PRATT, mixing<br />&gt;&nbsp;  1 1 3&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; iatom nlorb, lorb<br />&gt;&nbsp;  2 1 3&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; iatom nlorb, lorb<br />&gt;&nbsp; &nbsp; 1&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; nsic 0..AFM, 1..SIC, 2..HFM<br />&gt;&nbsp; &nbsp; 0.625 0.00&nbsp; &nbsp; &nbsp; &nbsp; U J (Ry)&nbsp;  Note: we recommend to use U_eff = U-J and J=0<br />&gt;&nbsp; &nbsp; 0.625 0.00&nbsp; &nbsp; &nbsp; &nbsp; U J (Ry)&nbsp;  Note: we recommend to use U_eff = U-J and J=0<br /><br />&gt; the calculation are performed as follow :<br />&gt; <br />&gt; runsp_lapw -p -fc 0.1<br />&gt; save_lapw case1<br /><br />&lt;here lapwdm?&gt;<br /><br />&gt; runsp_lapw -orbc -p -fc 0.1<br />&gt; save_lapw case2<br />&gt; runsp_lapw -orb -p -fc 0.1.<br />&gt; <br />&gt; x lapw2 -qtl -c -up&nbsp; -p<br />&gt; x lapw2 -qtl -c -dn&nbsp; -p<br /><br />This is OK, but I'm not sure whether your -orbc step has any effect. With orbc, the orbital potential is calculated from the existing <span class="yshortcuts">density matrix</span> (case.dmat*), without updating these later on. In your description, there are no dmat files present yet (I did not verify how wien2k deals with that). You could run lapwdm (with the same options as lapw2) once at that stage.<br /><br />Can you show by output where you get this information on 'two different fermi energies'? You have to look in case.scf, not elsewhere.<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>