On 12/26/13, 4:52 PM, lloyd riggs wrote:
Yes I know what you mean, however when your looking over a unit cell and trying to define a single molecule in the solvent (even as represented by the entire solvent makup), does it properly alighn the energy based on the averidged "zero" point. I ask, not knowing, but playing with this (single small molecules), I found the systems standard error swamps the potential energy, and you then get a proper shaped graph, where the axis has to be corrected by an energy term (is what I meant by zeroing). I did this with just h2o playing around out of boardom at some point...as I had run across something showing perfect sinusoidal energy graphs of waters/solvent...but the trajectory used contained 5 ions and a single protein....DMSO is a tiny molecule...
We're not dealing with a single molecule, at least assuming I understand what the OP is discussing. One has a homogeneous box of solvent, such that the potential energy of the system can be divided evenly between equivalent molecules and averaged over time. The literature often reports the "configurational energy" of such a liquid.
I'm not entirely sure what it is that you were playing around with, but such decompositions as these are not suited for heterogeneous systems, and what you've come up with really just sounds like a way of looking at the interconversion of kinetic and potential energy.
-Justin
Stephan *Gesendet:* Donnerstag, 26. Dezember 2013 um 14:04 Uhr *Von:* "Justin Lemkul" <jalem...@vt.edu> *An:* gmx-us...@gromacs.org *Betreff:* Re: [gmx-users] Potential energy calculations On 12/26/13, 5:30 AM, lloyd riggs wrote: > Zero it--> If you plot a potential energy value strait from the simulation it > will be something like point a) -7000 point b) -7050 across a single > simulation. Thus you have to find a beginning point and subtract it across the > run. Additionally, fluctuations may dictate 10 or so runs to determine a good > mean, as the drifet or error might be 20 kcal/mol in a system. As I read your > other post, I would assume the experiment would be non random placement of DMSO > in water (say a wall system), otherwise you would be looking at a single > molecules energy dictated by conformational space over time...ie something like > Justin's post...which you still may have to calculate a mean for, but I assume > this is done through the -nmol option. You still also may have to "zero it", > but I am not sure, direct use of g_energy gives abitrary starting points as > above from my experience. > The values printed by g_energy are not arbitrary. They are the energy of a mole of equivalent systems of a given configuration, hence their large magnitude. -Justin -- ================================================== Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 ================================================== -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.
-- ================================================== Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 ================================================== -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.
