Hi Emine, <div><br></div><div>Thanks a ton for your kind reply.  <div><br><div>For lead (Pb), if one gives Hubbard_5d = y eV, is this enough ?</div><div>Or, one also has to have a Pb pseudopotential where 10 electrons </div>
<div>of 5d-state must be treated as valence states ? </div><div><br></div><div>I mean, is this acceptable,  if one has a Pb pseudo with  6s2 6p2 </div><div>(only four electrons in valence)  and then she/he wants to apply </div>
<div>Hubbard_5d = y eV for the 5d-orbital which is fully frozen inside the core ?</div><div><br></div><div>Furthermore, is this at all permitted to give Hubbard U_f  and Hubbard U_d</div><div>together for the same element ? In some cases in the literature, the Hubbard U_p </div>
<div>and Hubbard U_d together have been used for the same Titanium </div><div>(Ti) element.  Is this physically meaningful ? </div><div>  <br><div>Regards, </div><div>Suza</div><div><br><div class="gmail_quote">On Sat, Apr 27, 2013 at 11:16 PM, Kucukbenli Emine <span dir="ltr"><<a href="mailto:emine.kucukbenli@epfl.ch" target="_blank">emine.kucukbenli@epfl.ch</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">




<div>
<div style="direction:ltr;font-size:10pt;font-family:Tahoma">
<div>Hi Suza,</div>
<div><br>
</div>
<div>there are two places in the code which you would need to modify to be able to use a different angular momentum than the default.  Dont be scared though, they are fairly easy to understand and modify:</div>
<div><br>
</div>
<div>espresso/flib/set_Hubbard_l.f90</div>
<div>tells which angular momentum to use for each atom type (very transparent: Hubbard_l 2 is for d and 1 is for p)</div>
<div><br>
</div>
<div>espresso/PW/src/tabd.f90</div>
<div>tells how to distribute initial occupation on the Hubbard projector states </div>
<div><br>
</div>
<div><span style="font-size:10pt">Change these files to your wish and re-make pw and you will be set.</span></div>
<div><br>
</div>
<div>best,</div>
<div>emine kucukbenli, postdoc at theos, epfl, switzerland</div>
<div><br>
</div>
</div>
</div>

<br><br></blockquote></div><br></div></div></div></div>