[gmx-users] Enthalpy of Vaporization
jalemkul at vt.edu
Mon Oct 5 20:04:03 CEST 2015
On 10/5/15 2:01 PM, Rasoul Nasiri wrote:
> On Mon, Oct 5, 2015 at 6:52 PM, Justin Lemkul <jalemkul at vt.edu> wrote:
>> On 10/5/15 1:50 PM, Pallavi Banerjee wrote:
>>> Correct me if I am wrong, please. Ugas would be obtained from a simulation
>>> of a single molecule in vacuum. And U liq would be the total potential
>>> energy of the system (which is given by g_energy) minus the total bonded
>>> potential energy (bonds+angles+dihedrals). Am I thinking right?
>> No. It's the actual potential energy. You need the whole potential
>> energy to account for flexibility in the molecule.
> Hopefully FFs can do this because of the anharmonicity effects caused by
> multi-structural effects (e.g., many coupling of torsions). This
> scenario will be serious when T goes up.
Indeed it depends to some extent on the size, flexibility, and chemical
composition of the molecule. Additive FFs will inherently overestimate dipole
moments to account for a lack of polarization response, so generally very polar
molecules have poor DHvap values with additive models. Simple hydrophobics are
easier and often agree better with experimental data.
Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow
Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 629
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201
jalemkul at outerbanks.umaryland.edu | (410) 706-7441
More information about the gromacs.org_gmx-users