<div dir="ltr"><div>See this note in file Doc/release-notes:<br><br>  * Some constants in the definition of PBE functionals were truncated to <br>    6 significant digits. While not a bug, this could lead to tiny differences<br>    with respect to previous results and other XC implementations (r13592)<br><br></div><div>More explicitly: some constants that previously were truncated to 6 digits, are now defined with full precision. This removes minor numerical discrepancies with the results obtained using libxc.<br><br>Paolo<br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Dec 21, 2017 at 3:11 AM, WF <span dir="ltr"><<a href="mailto:kalamaillist@gmail.com" target="_blank">kalamaillist@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div link="#0563C1" vlink="#954F72" lang="ZH-CN"><div class="m_-2017852202558475452WordSection1"><p class="MsoNormal"><span style="font-size:10.5pt" lang="EN-US">Dear all,<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US">In tests a normal DFT calculations with Si, I noticed that QE 6.0 and 6.1 gives exactly same total energy with 6.0 and 6.1 (-15.75190880 Ry), but with QE 6.2.1, the total energy is -15.75191309Ry. I tested this on two different machines with different compilers, the results are the same across all of them (at least to last digit). I did some tests on larger systems also and found even larger difference. This 1e-5 Ry order of difference is not negligible. Is there anything introduced in 6.2 leads to this difference? Thanks very much.<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt" lang="EN-US">Best regards,<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt" lang="EN-US">Feng<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US">The input file is attached:<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US">====================<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US">&control<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US">    calculation='scf'<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US">    outdir='./tmp'<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US">    prefix='si'<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US">    wf_collect = .true.,<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US">/<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US"> &system<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US">    ibrav=2,<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US">    celldm(1) = 10.263,<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US">    nat=2,<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US">    ntyp= 1,<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US">    ecutwfc = 40,<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US">/<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US"> &electrons<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US">    mixing_mode='plain'<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US">    mixing_beta = 0.7,<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US">    conv_thr =  1.0d-10,<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US"> /<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US">ATOMIC_SPECIES<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US">Si  1.0 Si_ONCV_PBE-1.1.upf<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US">ATOMIC_POSITIONS {alat}<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US">Si 0.00 0.00 0.00<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US">Si 0.25 0.25 0.25<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US">K_POINTS (automatic)<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US">4 4 4 0 0 0<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US">==============<u></u><u></u></span></p><p class="MsoNormal" style="text-indent:21.0pt"><span style="font-size:10.5pt" lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal" style="text-align:justify;text-justify:inter-ideograph"><span style="font-size:10.5pt" lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal" style="text-align:justify;text-justify:inter-ideograph"><span style="font-size:10.5pt" lang="EN-US">------------------------------<wbr>--------------------<u></u><u></u></span></p><p class="MsoNormal" style="text-align:justify;text-justify:inter-ideograph"><span style="font-size:10.5pt" lang="EN-US">Dr. Feng Wu<u></u><u></u></span></p><p class="MsoNormal" style="text-align:justify;text-justify:inter-ideograph"><span style="font-size:10.5pt" lang="EN-US">Postdoctoral researcher<u></u><u></u></span></p><p class="MsoNormal" style="text-align:justify;text-justify:inter-ideograph"><span style="font-size:10.5pt" lang="EN-US">Department of Chemistry and Biochemistry<u></u><u></u></span></p><p class="MsoNormal" style="text-align:justify;text-justify:inter-ideograph"><span style="font-size:10.5pt" lang="EN-US">University of California, Santa Cruz<u></u><u></u></span></p><p class="MsoNormal" style="text-align:justify;text-justify:inter-ideograph"><span style="font-size:10.5pt" lang="EN-US">Phone: <a href="tel:(831)%20459-2874" value="+18314592874" target="_blank">831-459-2874</a><u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt" lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p></div></div><br>______________________________<wbr>_________________<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/<wbr>listinfo/pw_forum</a><br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="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>