<div dir="ltr"><div>Dear Lorenzo,</div><div><br></div><div>Thank you so much.</div><div><br></div><div>Best<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Sep 16, 2020 at 1:19 PM Lorenzo Paulatto <<a href="mailto:paulatz@gmail.com">paulatz@gmail.com</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">>> By changing the threshold values you mentioned, nscf calculation finished (no "cholesky" error)...<br>
This is a common problem and maybe would require an upstream fixing:<br>
<br>
when doing an SCF calculation, the diagonalisation threshold is set to a <br>
relatively loose value at the beginning, then tightened.<br>
<br>
      IF ( iter == 2 ) ethr = 1.D-2<br>
      ethr = MIN( ethr, 0.1D0*dr2 / MAX( 1.D0, nelec ) )<br>
      ! ... do not allow convergence threshold to become too small:<br>
      ! ... iterative diagonalization may become unstable<br>
      ethr = MAX( ethr, 1.D-13 )<br>
(from electrons.f90)<br>
<br>
On the other hand, in an NSCF calculation, the threshold for <br>
diagonalisation is proportional to the input convergence threshold.<br>
<br>
IF ( ethr == 0.D0 ) ethr = 0.1D0 * MIN( 1.D-2, tr2 / nelec )<br>
(from setup.f90)<br>
<br>
The latter can turn out to be much tighter than the former, and even <br>
impossibly tight sometimes which makes the algorithm fail. Maybe adding <br>
a minimu of 1.d-13 also for NSCF would fix this kind of problems, but <br>
I'm not 100% sure this is exactly how it is done, as there are several <br>
places where ethr is set.<br>
<br>
cheers<br>
<br>
<br>
<br>
-- <br>
Lorenzo Paulatto - Paris<br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
</blockquote></div>