<div dir="ltr">Dear Joshua<br><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Mar 2, 2015 at 4:27 PM, Joshua Davis <span dir="ltr"><<a href="mailto:davis101@chemistry.msu.edu" target="_blank">davis101@chemistry.msu.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div>Dear QE Users<br><br>I want to know how to calculate U based on <a href="http://journals.aps.org/prb/abstract/10.1103/PhysRevB.71.035105" target="_blank">linear response</a>, but I  am having trouble understanding what occupation values to use from a calculation for a response matrix.  <br><br></div>I have fallowed the tutorials<br><br><a href="http://media.quantum-espresso.org/santa_barbara_2009_07/slides-exercices/Handson_ldau.pdf" target="_blank">http://media.quantum-espresso.org/santa_barbara_2009_07/slides-exercices/Handson_ldau.pdf</a><br><br></div>and <br><br><a href="http://hjklol.mit.edu/content/calculating-hubbard-u" target="_blank">http://hjklol.mit.edu/content/calculating-hubbard-u</a><br><br><br></div></div></div></blockquote><div><br></div><div>you have to use the trace of the occupation matrix.<br></div><div><br> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div></div>Do I use the Trace of the d-electron occupation matrix (like in the first link) or am I to use sum of all the elements of the matrix (like in the second link)?<br><br></div><div>Secondly when calculating U, if the atoms that require a U are not very close to each other, can you get away with not calculating the off-diagonal response matrix elements and just invert the diagonal response matrix elements directly? <br></div></div></blockquote><div><br><br></div><div>I guess it depend on the system you are considering. just try in one case and see how results compare.<br><br></div><div>Matteo<br></div><div><br> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div></div><div><div><br></div><div>I am using QE 5.1.</div><div>Thanks for any help<br clear="all"><div><div><div><div><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><font face="arial, helvetica, sans-serif">----------------------------------------------------------------------------------------------------------------</font></div><div dir="ltr"><font face="arial, helvetica, sans-serif">Joshua D. Davis</font></div><div dir="ltr"><font face="arial, helvetica, sans-serif"><br>Graduate Assistant</font></div><div dir="ltr"><div><span style="background-color:transparent"><font face="arial, helvetica, sans-serif">Department of Chemistry</font></span></div><div><span style="background-color:transparent"><font face="arial, helvetica, sans-serif">Michigan State University</font></span></div><div><span style="background-color:transparent"><p style="margin:0in 0in 0pt"><font face="arial, helvetica, sans-serif"><span lang="EN-US"></span></font></p></span></div><div><font face="arial, helvetica, sans-serif"><span style="background-color:transparent">East Lansing, MI 48824</span><br></font></div><div><span style="background-color:transparent"><font face="arial, helvetica, sans-serif">-----------------------------------------------------------------------------------------------------------------</font></span></div></div></div></div></div></div>
</div></div></div></div></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></div></div>