<!doctype html>

<html>

 <head> 

  <meta charset="UTF-8"> 

 </head>

 <body>

  <div class="default-style" style="font-size: 12pt; font-family: helvetica,arial,sans-serif; color: #333333;">

   <span style="font-family: helvetica; font-size: 12pt;">Hello,</span>

  </div> 

  <div class="default-style" style="font-size: 12pt; font-family: helvetica,arial,sans-serif; color: #333333;">

   <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 class="default-style" style="font-size: 12pt; font-family: helvetica,arial,sans-serif; color: #333333;">

  </div> 

  <div class="default-style" style="font-size: 12pt; font-family: helvetica,arial,sans-serif; color: #333333;">

   Error in routine cdiaghg (797):

   <br>S matrix not positive definite

  </div> 

  <div class="default-style" style="font-size: 12pt; font-family: helvetica,arial,sans-serif; color: #333333;">

   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 class="default-style" style="font-size: 12pt; font-family: helvetica,arial,sans-serif; color: #333333;">

   Thank you for any potential help.

  </div> 

 </body>

</html>