<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none"><!-- p { margin-top: 0px; margin-bottom: 0px; }--></style>
</head>
<body dir="ltr" style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">
<p>Dear QE developers,</p>
<p><br>
</p>
<p>I'm using the PWSCF v.6.2 (svn rev. 13949:13950) interface with LibXC (version 4.0.2) to perform some calculations with meta-GGA functionals. I noticed that even for a very simple case, diamond computed with TPSS, QE interfaced with LibXC yields NaNs and
stops early; however, the native TPSS implementation runs successfully.</p>
<div><br>
</div>
<div>I started to look the source code for clues about this odd behavior and I discovered an odd 1/2 factor in the kinetic energy density when LibXC is used.<br>
</div>
<div>As far as I understand the code, that factor should not be present.<br>
</div>
<div><br>
</div>
<div>I kept digging into the code and found the following:</div>
<div><br>
</div>
<div>1. The NaNs were coming from the correlation part of TPSS. </div>
<div><br>
</div>
<div>2. By skipping the TPSS correlation call, I could manage to converge the SCF. However, the results using only TPSS exchange for both LibXC and the native implementation are different.</div>
<div><br>
</div>
<div>2. Getting rid of the 1/2 factor in tau (and the corresponding 1/2 factor in v3x) yields essentially the same result for LibXC and the native implementation using only TPSS exchange.<br>
</div>
<div><br>
</div>
<div>3. Getting rid of the 1/2 factor allows to use of LibXC TPSS correlation. Furthermore, both the LibXC and the native implementation yield the same result.<br>
</div>
<div><br>
</div>
<div>4. From my work with LibXC with other codes, I know that this interface can be rather tricky. In order to further stabilize the PWSCF and LibXC interface I added the same density threshold as in the native case, namely, skipping the computation of the
potential when the density gets below 1.0E-10 a.u.<br>
</div>
<div><br>
</div>
<div><br>
</div>
<div>Note that the 1/2 factor appear in all calls to LibXC meta-GGA potentials, so it should be corrected for all functionals.<br>
</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div>I hope that my explanation is not too cumbersome. I'm also attaching the input file I'm using (C.in) and the output files for the native (C.out) and the corrected LibXC implementation (C.out.2). <br>
</div>
<div><br>
</div>
<div>If you have further questions, please do not hesitate to ask.<br>
</div>
<p><br>
</p>
<p><br>
</p>
<div id="Signature">
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0" name="divtagdefaultwrapper">
<font face="Times,'Times New Roman',Serif" size="3"><font color="#000000"><b><i>Daniel Mejia Rodriguez</i><br>
<font size="2"><br>
Postdoctoral Associate<br>
Q<font face="Times,'Times New Roman',Serif">TP <font face="Times,'Times New Roman',Serif">
Department of Physics</font></font></font></b></font><font size="2"><br>
University of Florida<br>
<font face="Times,'Times New Roman',Serif">2302 New Physics Building<br>
<font face="Times,'Times New Roman',Serif">Gainesville, FL 32611</font><br>
</font>(352)<font face="Times,'Times New Roman',Serif"> </font>392-9227</font></font><br>
</div>
</div>
</body>
</html>