<div dir="ltr"><div>I think I did the vdW-df porting to CP. I am quite confident that the fixed-cell case works properly (if not, I must expunge a paper from my list of publications!). The variable-cell case may have received no seriously testing, though. I'll have a look sooner or later. If later, I'll disable that case.<br><br></div>Paolo<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Sep 13, 2016 at 8:52 AM, Filippo SPIGA <span dir="ltr"><<a href="mailto:filippo.spiga@quantum-espresso.org" target="_blank">filippo.spiga@quantum-espresso.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello Pier Luigi ,<br>
<br>
I forward your email to the q-e-developers mailing-list. I am not expert of the physics, I am more on the "computing/HPC" side of things. Hopefully someone can clarify this point of you. If it is a bug (it looks like a bug and it smells like a bug.... so probably we need to investigate properly) we will work to fix it by the official release My uninformed guess is that the quantities that varies across different machines is calculated heavily using some libraries or some OpenMP. This is something easy to check. I cannot comment about cp.x versus pw.x correctness.<br>
<br>
Anyone has a guess about these discrepancies or a practical suggestion how to verify this?<br>
<br>
<br>
On Sep 13, 2016, at 7:21 AM, Pier Luigi Silvestrelli <<a href="mailto:pierluigi.silvestrelli@pd.infn.it">pierluigi.silvestrelli@pd.<wbr>infn.it</a>> wrote:<br>
<br>
I have found that in QE (in the new 6.0-beta version<br>
and in previous versions as well) the pressure calculation seems<br>
to be incorrect using 'cp.x' in combination with non-local vdW-corrected<br>
functionals, such as vdW-DF, vdW-DF2, vdW-DF-cx, rVV10,...,<br>
as verified by comparing results obtained using both 'pw.x' and<br>
'cp.x' programs.<br>
For instance, considering a single liquid-water configuration<br>
and the vdW-DF functional, 'pw.x' (both using the 'scf' and 'md'<br>
options) gives P=-3.55 kbar, while 'cp.x' gives instead P=200.99 kbar.<br>
Moreover, P computed by 'cp.x' using the same version of the code<br>
on different machines gives different results, while all the other<br>
quantities (total energy,...) are identical.<br>
<br>
I would like to let you know this problem, which is of some importance<br>
for efficient vdW-corrected ab initio MD simulations of liquids.<br>
<br>
Thanks in advance. Ciao.<br>
<br>
Pier Luigi Silvestrelli<br>
<br>
Dipartimento di Fisica e Astronomia G. Galilei<br>
Via Marzolo 8<br>
I-35131 Padova, Italy<br>
tel. +39-049-8277171<br>
fax  +39-049-8277102<br>
e-mail : <a href="mailto:psil@pd.infn.it">psil@pd.infn.it</a><br>
<br>
--<br>
Filippo SPIGA ~ Quantum ESPRESSO Foundation ~ <a href="http://www.quantum-espresso.org" rel="noreferrer" target="_blank">http://www.quantum-espresso.<wbr>org</a><br>
<br>
<br>
______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">http://qe-forge.org/mailman/<wbr>listinfo/q-e-developers</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, <span>Dip. Scienze Matematiche Informatiche e Fisiche</span>,<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>