<div dir="ltr"><br><div class="gmail_extra"><div class="gmail_quote">On Sat, Jun 6, 2015 at 4:40 PM, Jiqiang Li <span dir="ltr"><<a href="mailto:jqli14@fudan.edu.cn" target="_blank">jqli14@fudan.edu.cn</a>></span> wrote:<div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div><font face="Arial"><span style="font-size:10.5pt;line-height:1.5;background-color:window"></span><span style="font-size:10.5pt;line-height:1.5;background-color:window"></span></font><span style="font-size:10.5pt;line-height:1.5;background-color:window"><font face="Arial">a) I am wondering whether the integral can be simply handled as \sum_i {V(i) * rho(i)} * volume / N, where volume of unit cell and N the size of array for V and rho. If no, any other formula?</font></span></div></div></blockquote><div><br><div><div>yes: all integrals in QE are computed that way, on the Discrete Fourier Transform grid<br></div></div><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><div><font face="Arial"><span style="font-size:10.5pt;line-height:1.5;background-color:window">b) Is the integral above equal to the sum of '</span><span style="font-size:10.5pt;line-height:1.5;background-color:window">hartree contribution" and "</span><span style="font-size:10.5pt;line-height:1.5;background-color:window">xc contribution"  in the stdout file?</span></font></div></div></blockquote><div><br></div><div>no: "xc contribution" =  variable etxc = <span style="font-size:10.5pt;line-height:1.5;background-color:window"><font face="Arial"> \int \epsilon_{xc}(r) * rho(r) d^3r <br>while </font></span><span style="font-size:10.5pt;line-height:1.5;background-color:window"><font face="Arial">\int V_{xc}(r) * rho(r) d^3r is contained in variable vtxc (not printed in output).<br></font></span> <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><div><font face="Arial"><span style="font-size:10.5pt;line-height:1.5;background-color:window"> Or where can I find some information in the stdout file as benchmark data?</span></font></div><div><span style="font-size:10.5pt;line-height:1.5;background-color:window"><font face="Arial">  </font></span></div><div><span style="font-size:10.5pt;line-height:1.5;background-color:window"><font face="Arial">     Any consideration will be gratefully appreciated.</font></span></div><div><span style="font-size:10.5pt;line-height:1.5;background-color:window"><font face="Arial"><br></font></span></div><div><span style="font-size:10.5pt;line-height:1.5;background-color:window"><font face="Arial"><br></font></span></div><div><span style="font-size:10.5pt;line-height:1.5;background-color:window"><font face="Arial">Best regards!</font></span></div><span class=""><font color="#888888"><div><span style="font-size:10.5pt;line-height:1.5;background-color:window"><font face="Arial"><br></font></span></div><div><span style="font-size:10.5pt;line-height:1.5;background-color:window"><font face="Arial">Jiqiang Li</font></span></div><div><span style="font-size:10.5pt;line-height:1.5;background-color:window"><font face="Arial">Fudan University, China</font></span></div><hr style="width:210px;min-height:1px" color="#b5c4df" align="left" size="1">
<div><span></span></div>
</font></span></div><br>_______________________________________________<br>
Q-e-developers mailing list<br>
<a href="mailto:Q-e-developers@qe-forge.org">Q-e-developers@qe-forge.org</a><br>
<a href="http://qe-forge.org/mailman/listinfo/q-e-developers" target="_blank">http://qe-forge.org/mailman/listinfo/q-e-developers</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr"><span><span><font color="#888888">Paolo Giannozzi, Dept. Chemistry&Physics&Environment,<br>
Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
hone <a href="tel:%2B39-0432-558216" value="+390432558216" target="_blank">+39-0432-558216</a>, fax <a href="tel:%2B39-0432-558222" value="+390432558222" target="_blank">+39-0432-558222</a></font></span></span></div></div>
</div></div>