<br><br><div class="gmail_quote">On 31 March 2011 12:58, Justin A. Lemkul <span dir="ltr"><<a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br>
<br>
Elisabeth wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br>
<br>
in your mail:<div class="im"><br>
<br>
On 30 March 2011 15:30, Justin A. Lemkul <<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a> <mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>>> wrote:<br>
<br>
<br>
<br>
Elisabeth wrote:<br>
<br>
Dear all,<br>
<br>
I intend to obtain vaporization heat per volume for a /pure<br>
alkane system/. Here is the steps I am taking. Please correct me.<br>
<br>
1- Obtain total energy of system (kinetic+potential) and divide<br>
by number of molecules to obtain energy per mol of molecules.<br>
g_energy -f *.edr -nmol XXX<br>
2- Obtain total energy of a single molecule (use pbc).<br>
3- Subtract step 2 from step 1.<br>
4- Divide by simulation box volume.<br>
<br>
My questions is:<br>
<br>
in step 2 : what should be the box size? The same size as in 1<br>
or it does not matter? (step 1 is done for the actual denstiy)<br>
<br>
<br>
More troubling, how does one define the energy of a molecule? If<br>
you use any sort of long-range algorithms (especially PME, but also<br>
dispersion correction), you can't simply decompose the system like this.<br>
<br>
Thanks Justin and David.<br>
<br>
I have been trying to find the article in which this has been presented. If you have time Please see page 5937, right column, equation 11. I think I made a mistake and I dont have to include kinetic energy, Only nonboded energies!?<br>
<br>
<a href="http://pubs.acs.org/doi/pdfplus/10.1021/jp0707539" target="_blank">http://pubs.acs.org/doi/pdfplus/10.1021/jp0707539</a><br>
<br>
</div></blockquote>
<br>
What is cohesive energy and how does it relate to the quantity you're trying to calculate?<div class="im"><br>
It is delta Hvap/volume. It is directly related to Hvap. What is happening is that they are calculating nonbonded energy of some chains, divide by number of chains and substract from nonbonded energy of a single chain in vacuum. These are the steps I wrote in my first post but I think I should not have included kinetic and should just look at LJ-SR and Coulomb-SR.<br>
</div></blockquote><div><br>I am using PME..If I remember correctly LR is included in Coulomb-SR and can not get decomposed? But I dont think this doesnt matter since if I am to take nonbonded energies this should not hurt,,,<br>
<br>Please comment ... </div><div><br> </div><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div class="im">
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br>
In the derivation of recent Gromos96 parameter sets, the heat of<br>
vaporization is quite simple:<br>
<br>
DHvap = <Ugas> - <Uliq> + RT<br>
<br>
<br>
1- So <Uliq> is the total energy or only potential (no kinetic)<br>
</blockquote>
<br></div>
Potential.<div class="im"><br>
<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
2- How can I compute <Ugas>? I have liquid now...<br>
<br>
</blockquote>
<br></div>
Run a simulation in the gas phase.<br></blockquote><div><br>Sorry, but how can I do this? :( I have box of molecules with density of actual liquid..How can I shift to gas phase ..I mean how many molecules I need to keep in the box..<br>
<br>Many thanks... <br></div><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br>
-Justin<br>
<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div class="im">
<br>
Thank you,<br>
Regards,<br>
<br>
<br>
-- ========================================<br>
<br>
Justin A. Lemkul<br>
Ph.D. Candidate<br>
ICTAS Doctoral Scholar<br>
MILES-IGERT Trainee<br>
Department of Biochemistry<br>
Virginia Tech<br>
Blacksburg, VA<br></div>
jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> <<a href="http://vt.edu" target="_blank">http://vt.edu</a>> | (540) 231-9080<div class="im"><br>
<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br>
<br>
========================================<br>
<br>
-- gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br></div>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><div class="im"><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at<br>
<a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
Please don't post (un)subscribe requests to the list. Use the www<br>
interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br></div>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>>.<div class="im"><br>
Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
<br>
<br>
</div></blockquote><div><div></div><div class="h5">
<br>
-- <br>
========================================<br>
<br>
Justin A. Lemkul<br>
Ph.D. Candidate<br>
ICTAS Doctoral Scholar<br>
MILES-IGERT Trainee<br>
Department of Biochemistry<br>
Virginia Tech<br>
Blacksburg, VA<br>
jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> | (540) 231-9080<br>
<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br>
<br>
========================================<br>
-- <br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
Please don't post (un)subscribe requests to the list. Use the www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
</div></div></blockquote></div><br>