<div dir="ltr">On Tue, Apr 12, 2016 at 2:54 PM, Uri Argaman <span dir="ltr"><<a href="mailto:argamanu@post.bgu.ac.il" target="_blank">argamanu@post.bgu.ac.il</a>></span> wrote:<br><br><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div style="font-size:12.8px"><div><div><div><div><div><div><div>I wish to separate the "xc contribution" at the scf.out file into "exchange contribution" and "correlation contribution". What are you recommend?<br></div></div>1. to use the subroutine exxenergy2 from exx.f90 and substract the results from the "xc contribution"?<br></div></div></div></div></div></div></div></blockquote><div><br></div>exxenergy2 compute "the Hartree-Fock exchange", not "the exchange contribution" in "the exchange-correlation energy". Is this what you want? <br><br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div style="font-size:12.8px"><div><div><div><div><div></div>2. to do something similar to what is written in v_of_rho.f90 for each of the contributions?<br></div></div></div></div></div></div></blockquote><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div style="font-size:12.8px"><div><div><div><div></div>3. maybe something smarter?<br></div></div></div></div></div></blockquote><div><br></div><div>it is not easy to disentangle the two contributions because they are intertwined and you have to go quite deep into the code compute them separately. The simplest way is I think to set to zero the correlation terms (variables icorr, igcc in Modules/funct.f90) and recompute the energy using v_of_rho<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div style="font-size:12.8px"><div><div><div><br></div>I found similar question in the forum from 2014 with the answer that I have to modify the code. </div></div></div></div></blockquote><div><br></div><div>this is still valid (and will always be, I think)<br><br></div><div>Paolo <br><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div style="font-size:12.8px">thanks in advance,<br></div><div style="font-size:12.8px"><br></div><span style="font-size:12.8px">Uri Argaman</span><br><div><span style="font-size:12.8px">Materials engineering department </span></div><div><span style="font-size:12.8px">Ben-Gurion university </span></div><div><span style="font-size:12.8px">Israel</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" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>
</div></div>