<div dir="ltr"><div>If the structure is correct and pseudopotentials reasonably good, there isn't much that can be done. You may want to use "ParO" or "PPCG" diagonalization instead of "cg" that is very slow.</div><div><br></div><div>Paolo<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Jun 7, 2021 at 6:44 PM John Stoppelman <<a href="mailto:jstoppelman@comcast.net">jstoppelman@comcast.net</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"><u></u>

  
   
 
 <div>
  <div style="font-size:12pt;font-family:helvetica,arial,sans-serif;color:rgb(51,51,51)">
   <span style="font-family:helvetica;font-size:12pt">Hello,</span>
  </div> 
  <div style="font-size:12pt;font-family:helvetica,arial,sans-serif;color:rgb(51,51,51)">
   <span style="font-family:helvetica;font-size:12pt">I am running a vc-relax calculation for a crystal structure, and I'm trying to compare the results between PBE and PBEsol. The SCF converges with PBE using the default Davidson diagonalization. When I switch to PBEsol pseudopotentials from the PSlibary, I receive the error: </span>
  </div> 
  <div style="font-size:12pt;font-family:helvetica,arial,sans-serif;color:rgb(51,51,51)">
    
  </div> 
  <div style="font-size:12pt;font-family:helvetica,arial,sans-serif;color:rgb(51,51,51)">
   Error in routine cdiaghg (797):
   <br>S matrix not positive definite
  </div> 
  <div style="font-size:12pt;font-family:helvetica,arial,sans-serif;color:rgb(51,51,51)">
   When I switch the diagonalization method to 'cg', the calculation converges as normal, however. Are there any settings I can try in order to use PBEsol with the Davidson diagonalization for the SCF? I have attached my input file.
  </div> 
  <div style="font-size:12pt;font-family:helvetica,arial,sans-serif;color:rgb(51,51,51)">
   Thank you for any potential help.
  </div> 
 </div>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer" target="_blank">www.max-centre.eu</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></blockquote></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 206, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>