<div dir="ltr"><div>Hi all ,</div><div>I haven't encountered this kind of a problem with QE before.I was performing an optimization on degauss value. Upon a relaxation ,I obtained the following data :</div><div><table class="gmail-ajC" cellpadding="0"><tbody><tr class="gmail-ajv"><td colspan="2" class="gmail-gG"><span class="gmail-gI"></span><br></td><td colspan="2" class="gmail-gL"><br></td></tr></tbody></table>degauss = 0.02  Final energy = -4970.7298625102 Ry (kbar) P= 3.80<br>degauss = 0.04  Final energy = -4970.7331133155 Ry (kbar) P= 0.02<br>degauss = 0.05  Final energy = -4970.7357925053 Ry (kbar) P= -0.66<br>degauss = 0.06  Final energy = -4970.7395506085 Ry (kbar) P= -4.54<br>degauss = 0.07  Final energy = -4970.7451525972 Ry (kbar) P= -8.43.</div><div><br></div><div>I have already optimized the lattice parameters , ecut , mixing beta , k-pts etc. following a vc-relax. The threshold parameters e.g., etot_conv_thr , forc_conv_thr  & conv_thr in &electrons were set  accordingly to satisfy scf corrections to forces. As you can see that the final energy and the pressure are changing in an inverse way as the degauss decreases. <br></div><div>Is there any reason behind this kind of unstable behavior of P and total energy?<br></div><div>Further which set should I choose to proceed ? Is it the one with P=0.02 with somewhat larger energy  or the one with the minimum energy but with a larger negative pressure ? <br></div><div>Is there any reason behind this kind of unstable behavior of P and total energy?<br></div><div><br></div><div>QE version - 5.4.0</div><div><br></div><div>Thank you.<br></div><div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr">Sudip Kumar Mondal<div>DST-INSPIRE Junior Research Fellow</div><div>High Pressure Temperature Laboratory<br></div><div>Dept. Of Physics/Geological Sciences<br></div><div>Jadavpur University</div></div></div></div></div></div></div></div></div>