<div dir="ltr"><div><div><p>Dear 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 style="text-wrap-mode: wrap;"><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. Also, this issue is persistent with many noncubic materials.</p><p>Thank you for your assistance.</p></div><font color="#888888"><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></div></div></div></font></div><div><br></div></div>