<div dir="ltr">Dear QE users and developers<br><br>I do Fermi surface calculation and I got an error message: <br><div>"Error in routine cdiaghg (126):<br>  S matrix not positive definite"<br></div><div>I note that in some crystal geometries I got this message and in other geometries I didn't got this error and the Fermi surface looks fine and reasonable.  <br><br></div><div>This is the input file:<br>&control<br>    calculation='bands'<br>    restart_mode='from_scratch',<br>    prefix='ti',<br>    pseudo_dir='/home/uriargaman/QE/espresso-5.1/pseudo/',<br>    outdir='/home/uriargaman/tmp31971/TMP_FS/',<br>    tstress = .true. ,<br>    tprnfor = .true. ,<br>    verbosity='high',<br> /<br> &system<br>    ibrav=  4, celldm(1) =5.48090015121, celldm(3) =1.6120416161 , nat= 2, ntyp= 1,<br>    ecutwfc =50,ecutrho=500, <br>    occupations='smearing', smearing="methfessel-paxton", degauss=0.01,<br>     nbnd=16 <br> /<br> &electrons<br>    conv_thr =  1.0d-9,<br>    electron_maxstep=1000<br>    mixing_beta = 0.7d0,<br> /<br><br>ATOMIC_SPECIES<br>   Ti   47.9  Ti.bp-sp-van_ak.UPF<br>ATOMIC_POSITIONS (crystal)<br>   Ti      0.666666666    0.333333333    0.250000000    0  0  0 <br>   Ti      0.333333333    0.666666666    0.750000000    <br>K_POINTS <br>  4913<br>....<br><br></div><div>What can be the problem and how can I solve it?<br><br></div><div>Thank you in advance<br></div><div>Uri Argaman<br></div><div>Ben-Gurion University<br></div><div>Israel<br></div><div><br></div></div>