<div dir="ltr"><div>Dear users,</div><div><br></div><div>I have done vc-relax jobs using pw.x binary v.6.0 (svn rev. 13079). It completed  with this information in the final step<br></div><div><br></div><div>...<br>!    total energy              =   -2893.96255976 Ry<br></div><div>...</div><br><div>          total   stress  (Ry/bohr**3)                   (kbar)     P=    5.50<br></div><div>...<br><br>     Final enthalpy =   -2893.4953806269 Ry<br>Begin final coordinates<br>     new unit-cell volume =  12495.34754 a.u.^3 (  1851.61945 Ang^3 )<br></div><div><br></div><div>...</div><div><br></div><div>So I got the final enthalpy as reported by pw.x (call this H_reported) and I calculate the enthalpy from total energy using <br></div><div><br></div><div>H = E+PV</div><div>(called H_calc)<br></div><div><br></div><div><br></div><div>It turned out that H_reported and H_calc are not equal. <br></div><div><br></div><div>H_reported = -2893.4953806269 Ry = -39367.97664 eV</div><div>H_calc = E+PV=-39368.05688 eV<br></div><div><br></div><div>H_reported - H_calc = -0.08024 eV<br></div><div><br></div><div>So my questions are, <br></div><div><br></div><div>1) how does the pw.x determine final enthalpy in its report? <br></div><div>2) Is it theoretically correct to determine H=E+P*V directly from E, P and V reported in pw.x output. (Or should I determine P(V) by fitting the equation of state to E-V curve then derive H from that?)<br></div><div>3) If both the H_reported and H_calc are theoretically correct, why do the values differ?</div><div><br></div><div>My problem is not in this single case, it is that I did several relax jobs and the discrepancies in H_reported and H_calc are not systematic.<br></div><div><br></div><div>Thank you,</div><div>Vichawan Sakulsupich</div><div><br></div><div>MSc, Chulalongkorn University, THAILAND<br></div><div><br></div><div><br>-- <br><div class="gmail_signature">Vichawan Sakulsupich [Print]</div>
</div></div>