Dear Iwan,<div><br></div><div>You can use a large Hubbard_alpha on the Ti atoms to reduce the number of electrons on them and a large negative Hubbard_alpha to increase the number of electrons. Screen the trace of the occupation matrices to check whether you have the right number of electrons on 3d orbitals. Then you can compute U either by the linear response method, or the method presented in PRB 74, 045202 (2006). The second method presented in the paper requires computation of the dielectric tensor, which you can perform with the phonon code.</div>
<div><br></div><div>Best regards,</div><div>Burak<br><br><div class="gmail_quote">On Wed, Dec 19, 2012 at 7:59 PM, Iwan Darmadi <span dir="ltr"><<a href="mailto:iwan_darmadi@rocketmail.com" target="_blank">iwan_darmadi@rocketmail.com</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="font-size:12pt;font-family:garamond,new york,times,serif"><div><span>Dear All,<br></span></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:"Garamond"">
<span><br></span></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:"Garamond""><span>I am trying to calculate Hubbard U for Titanium like Anderson et al's paper (<a href="http://prb.aps.org/abstract/PRB/v74/i4/e045202" target="_blank">http://prb.aps.org/abstract/PRB/v74/i4/e045202</a>).<br>
</span></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:"Garamond""><span><br></span></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:"Garamond"">
<span>To do that, I need to calculate energy of Titanium with one
 electron removed/added from/to 3d orbital.</span></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:"Garamond""><span>My question is, Is manipulating pseudopotential of Titanium  the only way to remove/add electron in 3d orbital ?</span></div>
<div style="font-style:normal;font-size:16px;background-color:transparent;font-family:"Garamond""><span>Or Is there another feature in QE in which I could do so ?</span><span><br></span></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:"Garamond"">
<span><br></span></div><div> </div><div style="font-style:normal;font-size:13px;background-color:transparent"><span><font face="lucida console, sans-serif"><font><b>***<span class="HOEnZb"><font color="#888888"><br></font></span></b></font></font></span></div>
<span class="HOEnZb"><font color="#888888"><div style="font-style:normal;font-size:13px;background-color:transparent"><span><font face="lucida console, sans-serif"><font><b>Iwan Darmadi</b></font></font></span><br><span><font face="lucida console, sans-serif"> Undergrad.Student - Department of Physics</font></span><font face="lucida console, sans-serif"><br>
</font></div><div style="font-style:normal;font-size:13px;background-color:transparent"><span><font face="lucida console, sans-serif"> Universitas Indonesia</font></span><span><br></span></div></font></span></div></div><br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br><br clear="all"><div><br></div>-- <br>Burak Himmetoglu<div>Post-Doctoral Research Associate</div>
<div>University of Minnesota MN, USA</div>
</div>