<div dir="ltr">Dear QE community,<br><br>I have struggled for over a year to obtain reliable SCF elastic constants with thermo_pw when ionic relaxation is enabled. I know this query is more relevant to the thermo_pw community, unfortunately i haven't got any response from the themo_pw community so I am reaching out to this forum in anticipation:<br><br>The themo_control file is:<br>&INPUT_THERMO<br>  what = 'scf_elastic_constants',<br>  frozen_ions = .FALSE.<br>/<br>With `frozen_ions=.TRUE.` everything converges, but setting `.FALSE.` aborts with<br><br>```<br><b>task # 0<br>from bfgs : error # 1<br>dE0s is positive which should never happen</b><br>```<br>Lowering thresholds  (etot_conv_thr)  yields elastic constants, yet the values swing wildly (tens of GPa's) with minute changes in energy tolerances, given that the system I am working on is of low-symmetry (triclinic).<br><br><div>Additionally,following adjustments didn't yield any success:<br><br>- `trust_radius_max` reduced to 0.05 a.u.<br>    <br>- tighter `conv_thr`, `etot_conv_thr`, `forc_conv_thr`<br>    <br>- `ion_dynamics='fire'` and `'damp'` (BFGS still invoked in thermo_pw)<br><br>Has anyone resolved this coupled QE and thermo_pw issue  or identified safer optimization settings (e.g., alternative algorithms, pre-relax procedures, or recommended `trust_radius_max`) that stabilize the calculation?<br><br>Any guidance or example inputs would be greatly appreciated, I will be happy to share my input files if needed.<br><br>Thank you for your time.<br><br>Best regards,  <br>Vinay Maithani<br>IIT Kharagpur  <br><br></div></div><br><div class="gmail_quote gmail_quote_container"><div dir="ltr" class="gmail_attr">On Tue, Mar 4, 2025 at 12:08 PM Vinay Maithani <<a href="mailto:vinaymaithani21@kgpian.iitkgp.ac.in">vinaymaithani21@kgpian.iitkgp.ac.in</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><p>Dear QE community,</p><p>I am currently utilizing thermo_pw to compute the self-consistent field (SCF) elastic constants for various materials. While setting <code>frozen_ions=.TRUE.</code> in the <code>&INPUT_THERMO</code> namelist yields successful results, attempts with <code>frozen_ions=.FALSE.</code> lead to the following error:</p><pre><div><div dir="ltr"><code><b>task #         0
     from bfgs : error #         1
     dE0s is positive which should never happen</b>
</code></div></div></pre><p>To address this issue, I have implemented several strategies without success:</p><ol><li style="margin-left:15px">Reducing the <code>trust_radius_max</code> parameter to as low as 0.05.</li><li style="margin-left:15px">Adjusting convergence thresholds within the input file.</li><li style="margin-left:15px">Modifying <code>ion_dynamics</code> to 'fire' and 'damp'; however, the optimization persists with the BFGS algorithm.</li></ol><p>I seek your guidance on resolving this error to facilitate accurate elastic constant calculations with ionic relaxation.</p><p>Thank you for your assistance.</p></div><div><br></div><span class="gmail_signature_prefix">-- </span><br><div dir="ltr" class="gmail_signature"><div dir="ltr"><span style="color:rgb(34,34,34)">Regards,</span><br style="color:rgb(34,34,34)"><span style="color:rgb(34,34,34)">Vinay Maithani</span><br><div><span style="color:rgb(34,34,34)">Department of Metallurgical and Materials Engineering</span><br style="color:rgb(34,34,34)"><span style="color:rgb(34,34,34)">Indian Institute of Technology, Kharagpur</span><br></div></div></div></div>
</blockquote></div>