<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Dear developers,</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
I apologize if this is incorrect or discussed somewhere else, but I couldn't find anything.<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
I have a a few comments/questions about the exx routines in QE.</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<div>
<div>
<ol>
<li>On routines PAW_newdxx and PAW_xx_energy, the code is only executed on ionode which I believe is a bug when k-point and/or band parallelization is used.
<br>
</li><li>From what I see, tqr=.true. is not implemented for exx routines when k-points are used. The phases of the pair densities need be added, probably inside addusxx_r for this to be correct. In a handful of cases I've tested, I get significantly different results
with/without tqr with hybrid functionals and kpoints. Gamma point or single kpoint runs, of course, is fine. It would probably be better to perform a runtime check and abort if necessary.<br>
</li><li>Many of the executables in the PP/src folder rely on read_file to initialize/setup the calculation. The routine init_us_1() is called inside. When uspp/paw are used in combination with finite kpoints in qvan_init(...) (e.g. as in exx routines), qnorm needs
to be set correctly to get the interpolation grid of qrad set properly. Unfortunately, I can't calculate the necessary value of qnorm without knowing the kpoints, which are also read in read_file. Right now, I am forced to use an upper bound to qnorm (using
3.5), which I set before calling read_file. This might be fine, since the overhead is negligible and it probably doesn't matter in the bigger scheme of things, but it would be nice to be able to control this (the value of qnorm or the range of interpolation)
more transparently/easily. Just a wish! <br>
</li></ol>
</div>
<div><br>
</div>
<div>best regards and apologies if I'm wrong in any of my comments above,<br>
</div>
</div>
<div><br>
</div>
</div>
<div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div id="Signature">
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:#000000; font-family:Calibri,Helvetica,sans-serif">
<div style="font-family:Tahoma; font-size:13px">*************************************<br>
Miguel A. Morales<br>
Quantum Simulations Group, Physics Division<br>
Lawrence Livermore National Laboratory<br>
phone: 925-423-4956<br>
*************************************<br>
</div>
</div>
</div>
</div>
</body>
</html>