<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>
</head>
<body ocsi="0" fpstyle="1">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Dear all,<br>
My two cents on this infamous negative dr2 issue:<br>
<br>
I am not a fan of commenting it out, it often is a sign of some convergence/pseudo related issue,<br>
unless it is at the very beginning of an scf run..<br>
<br>
Now I am looking at Jack's files:<br>
<br>
1.The error is indeed at the beginning of a run, after a bfgs movement of atoms<br>
2.his ecutrho = ecutwfc*4, the default value.. which is surely not guaranteed to be a converged value for a PAW pseudo..<br>
<br>
By commenting out that error, we prevent him to look further into his systems convergence and check his pseudopotentials.<br>
<br>
My suggestion would be any sensible combination of the following:<br>
<br>
1.we put a check and only for the first steps, ignore if this error occurs, but report dr2 anyways<br>
2.write a more explanatory error message or put this in FAQ with  more explanation<br>
3.encourage users to report these pseudopotentials because i think it is most of the time due to a very hard pseudo..<br>
4.sometimes it is possible that the assumption that we make while writing paw contribution to dr2, that exchange could be neglected, is not accurate enough, and no matter how well behaving your pseudo is, high convergence, you might stil have a small negative
 dr2. only after one makes sure that that is the case, there could be an input option to ignore this error, so that people dont have to modify the code that way.
<br>
<br>
I would love to hear what you think.<br>
ciao<br>
emine kucukbenli, postdoc, theos, epfl, switzerland<br>
<br>
Quoting Paolo Giannozzi <giannozz@democritos.it>:<br>
<br>
> I sent the link to the related svn commit:<br>
><br>
>> http://qe-forge.org/gf/project/q-e/scmsvn/?action=browse&path=%2Ftrunk%2Fespresso%2FPW%2Fsrc%2Fscf_mod.f90&sortby=file&r1=9744&r2=9743&pathrev=9744<br>
><br>
> P.<br>
> --<br>
> Paolo Giannozzi, IOM-Democritos and University of Udine, Italy<br>
<br>
<br>
</div>
</body>
</html>