<html xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style>
</head>
<body lang="EN-US" link="#0563C1" vlink="#954F72">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt">Hello,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">I am developing a code, Sella (<a href="https://github.com/zadorlab/sella">https://github.com/zadorlab/sella</a> ), that can interface with a variety of electronic structure theory packages to perform saddle
point optimizations. Sella performs finite difference of the gradient vectors to approximate Hessian-vector products for an iterative eigensolver (operating on the Hessian matrix). Sella therefore requires a high degree of precision in the gradients. If the
error in the gradients is at or close to the same order of magnitude as the difference between gradients for two configurations separated by the finite displacement distance (typically 1e-4 Angstrom), then the iterative eigensolver will fail to converge. Increasing
the finite displacement distance also results in a poor approximation of Hessian-vector products due to nonlinearity, which also can cause the eigensolver to fail to converge.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">Some other codes such as GPAW enable use of the gradient error as an SCF convergence criterion. As far as I can tell, Quantum Espresso does not have this ability. In order to get the accuracy necessary for
my code to function properly, I have so far needed to specify an energy convergence threshold of conv_thr=1e-11. However, even with this absurdly tight convergence threshold, my iterative eigensolver still occasionally fails to converge.
<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">I am not particularly concerned about limiting the number of SCF iterations to achieve convergence. For my purposes, it is more important to have gradients with a consistent and tunable level of precision.
Lowering conv_thr increases the precision of the gradients, but not in a consistent and predictable way.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">Is there a better way to improve gradient precision in Quantum Espresso?<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">Thanks,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">Eric Hermes<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">Postdoctoral Researcher<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">Sandia National Labs<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">-- <o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</body>
</html>