<div dir="ltr">Dear QE users,<div><br></div><div>I just worked example 1 in espresso-5.3.0/Environ/examples/ and obtained the following result for a water molecule in water solvent:</div><div><br></div>!    total energy              =     -34.06234148 Ry<br>     Harris-Foulkes estimate   =     -33.95165940 Ry<br>     estimated scf accuracy    <          6.8E-12 Ry<br><br>     The total energy is the sum of the following terms:<br><br>     one-electron contribution =     -81.79999129 Ry<br>     hartree contribution      =      42.43346037 Ry<br>     xc contribution           =      -8.39845654 Ry<br>     ewald contribution        =      13.81332806 Ry<br>     solvation energy          =      -0.03439636 Ry<br>     cavitation energy         =       0.01360379 Ry<br>     PV energy                 =      -0.00691218 Ry<div><br></div><div>When I sum the energy components, I get -33.97936415 Ry, which is not equal to the total energy above. I read</div><div><span style="color:rgb(0,0,0);white-space:pre-wrap"><br></span></div><div><span style="color:rgb(0,0,0);white-space:pre-wrap">O. Andreussi, I. Dabo and N. Marzari, J. Chem. Phys. 136 064102 (2012)</span></div><div><font color="#000000"><span style="white-space:pre-wrap"><br></span></font></div><div><font color="#000000"><span style="white-space:pre-wrap">and the authors state that </span></font><span style="font-size:12.8px">dGsol = dGel + Gcav + <i>a</i>*S + <i>b</i>*V + PdV where <i>a</i> and <i>b</i> are fitting parameters to reproduce the experimental solvation energies for a set of reference molecules. Is the "solvation energy" above the same as dGsol? Is the total energy the sum of these terms or does the "solvation energy" contain the cavitation and PV energies?</span></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">Thanks in advance for your assistance!</span></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">Best,</span></div><div><span style="font-size:12.8px">Rob</span></div><div><div><div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><span style="font-size:12.8px">---------------------------------------------</span><br></div><div>Robert B. Wexler                          </div><div>PhD Candidate in Chemistry</div><div>University of Pennsylvania, 2018</div><div><a href="mailto:rwexler@sas.upenn.edu" target="_blank">rwexler@sas.upenn.edu</a></div><div>(215) 801-8741</div><div>---------------------------------------------<br></div></div></div></div></div></div>
</div></div></div>