[gmx-users] Re: Enthalpy Confusion

[Jeffery Perkins]

>Your calculation seems correct. Which gmx version did you use? 
>The correlation between the numbers is almost 100% so there must be a 
>simple explanation.

gmx version is 4.5.4, and yeah the correlation is odd, in the code you
listed:

>>              pv = vol*md->ref_p/PRESFAC; 
>> 
>>              add_ebin(md->ebin, md->ipv, 1, &pv, bSum); 
>>              enthalpy = pv + enerd->term[F_ETOT]

so pv = vol * md->ref_p/PRESFAC; is the vol * ref_p/conversion factor,
right?
the only thing there that seems like it could be off is if PRESFAC isn't
correct for some reason

(correct me if I'm wrong of course!)

Jeffery



--
View this message in context: http://gromacs.5086.x6.nabble.com/Enthalpy-Confusion-tp5009053p5009168.html
Sent from the GROMACS Users Forum mailing list archive at Nabble.com.