<meta http-equiv="Content-Type" content="text/html; charset=GB18030"><div>Dear QE Developer,</div><div><br></div><div>I apologize for bothering you with a question that has likely been asked many times on this forum. I am running into the following error message with the input_dft='scan' interface using the libxc library: "S matrix not positive definite or / the charge is wrong". I have searched the forum for solutions to this issue and have found that the problem does not arise when using the PBE functional.</div><div><br></div><div>Unfortunately, my goal is to accurately describe the structure and energy of my system using the SCAN functional. I have attempted to modify the encut energy and diagonalization methods with no success. I am using ONCV norm conserving pseudopotentials, and have included my input file below for reference.</div><div><br></div><div>&CONTROL</div><div><br></div><div>    calculation='vc-relax', disk_io='low', prefix='CrB',</div><div><br></div><div>    pseudo_dir='./', outdir='./tmp', verbosity='high'</div><div><br></div><div>    tprnfor=.true., tstress=.true., forc_conv_thr=1.0d-5</div><div><br></div><div>/</div><div><br></div><div>&SYSTEM</div><div><br></div><div>    ibrav= 0,</div><div><br></div><div>    celldm(1) = 1.8897261328856432, ! a.u. to Angst</div><div><br></div><div>    nat= 3, ntyp= 2,</div><div><br></div><div>    input_dft='scan0',</div><div><br></div><div>    occupations = 'smearing', smearing = 'gauss', degauss = 0.03</div><div><br></div><div>    ecutwfc= 150, ecutrho = 600,</div><div><br></div><div>/</div><div><br></div><div>&ELECTRONS</div><div><br></div><div>    conv_thr = 1.0d-9</div><div><br></div><div>    mixing_beta = 0.8</div><div><br></div><div>    diagonalization = 'ppcg'</div><div><br></div><div>/</div><div><br></div><div> &IONS</div><div><br></div><div>    ion_dynamics = 'bfgs'</div><div><br></div><div>    upscale = 20</div><div><br></div><div>/</div><div><br></div><div> &CELL</div><div><br></div><div>    cell_dynamics = 'bfgs'</div><div><br></div><div>    press = 0</div><div><br></div><div>/</div><div><br></div><div>ATOMIC_SPECIES</div><div><br></div><div>  Cr 51.966 Cr.UPF</div><div><br></div><div>  B  10.81 B.UPF</div><div><br></div><div>CELL_PARAMETERS (alat=  1.88972613)</div><div><br></div><div>     2.904617312   0.000000000   0.000000000</div><div><br></div><div>  -1.452308656   2.515472381   0.000000000</div><div><br></div><div>  -0.000000000  -0.000000000   2.825090285</div><div><br></div><div>ATOMIC_POSITIONS (crystal)</div><div><br></div><div>Cr 0  0 0</div><div><br></div><div>B 0.3333333429999996  0.6666666870000029  0.5000000000000000</div><div><br></div><div>B 0.6666666269999979  0.3333333129999971  0.5000000000000000</div><div><br></div><div>K_POINTS {automatic}</div><div><br></div><div>  20 20 12 0 0 0</div><div><br></div><div><br></div><div><br></div><div>I would greatly appreciate any advice or guidance you can offer on how to address this issue and successfully use the SCAN functional to accurately describe my system's structure and energy.</div><div><br></div><div>Thank you for your time and expertise.</div><div><br></div><div>Best regards,</div><div><br></div><div>zhouchao</div><div><br></div><div><br></div><div><br></div>