<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;">Dear gerrit sir,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; thank you for your reply Here I have got compilation&nbsp; error as follows<br><br>/root/gromacs-4.5.3/src/mdlib/.libs/libmd_d.a(qm_mopac.o): In function `call_mopac_SH':<br>qm_mopac.c:(.text+0x2a7): undefined reference to `domop_'<br>/<span style="font-weight: bold;">root/gromacs-4.5.3/src/mdlib/.libs/libmd_d.a(qm_mopac.o): In function `call_mopac':</span><br style="font-weight: bold;"><span style="font-weight: bold;">qm_mopac.c:(.text+0x70b): undefined reference to `domop_'</span><br style="font-weight: bold;"><span style="font-weight: bold;">/root/gromacs-4.5.3/src/mdlib/.libs/libmd_d.a(qm_mopac.o): In function `init_mopac':</span><br style="font-weight: bold;"><span style="font-weight:
 bold;">qm_mopac.c:(.text+0x964): undefined reference to `domldt_'</span><br style="font-weight: bold;"><span style="font-weight: bold;">collect2: ld returned 1 exit status</span><br style="font-weight: bold;"><span style="font-weight: bold;">make[3]: *** [grompp] Error 1</span><br style="font-weight: bold;"><span style="font-weight: bold;">make[3]: Leaving directory `/root/gromacs-4.5.3/src/kernel'</span><br style="font-weight: bold;"><span style="font-weight: bold;">make[2]: *** [all-recursive] Error 1</span><br style="font-weight: bold;"><span style="font-weight: bold;">make[2]: Leaving directory `/root/gromacs-4.5.3/src'</span><br style="font-weight: bold;"><span style="font-weight: bold;">make[1]: *** [all] Error 2</span><br style="font-weight: bold;"><span style="font-weight: bold;">make[1]: Leaving directory `/root/gromacs-4.5.3/src'</span><br style="font-weight: bold;"><span style="font-weight: bold;">make: *** [all-recursive] Error 1</span><br>i
 am using the following command&nbsp; <br>./configure --disable-float --with-qmmm-mopac CPPFLAGS=-DUSE_MOPAC LIBS=-lmopac LDFLAGS=-L/usr/local/lib<br><br><span style="font-weight: bold;">i am expecting&nbsp; your useful reply Thanks in Advance </span><br style="font-weight: bold;"><br><br><br><br><br><br><br><br><br><br><br>--- On <b>Tue, 7/12/10, Gerrit Groenhof <i>&lt;ggroenh@gwdg.de&gt;</i></b> wrote:<br><blockquote style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;"><br>From: Gerrit Groenhof &lt;ggroenh@gwdg.de&gt;<br>Subject: [gmx-users] Re: gmx-users Digest, Vol 80, Issue 41<br>To: gmx-users@gromacs.org<br>Date: Tuesday, 7 December, 2010, 1:24 AM<br><br><div class="plainMail">Without-qmmm-mopac will not lead to compilation of the mopac interface source code. <br><br>What is the compilation error?<br><br>Gerrit<br><br><br><br>&gt; <br>&gt;&nbsp;&nbsp;&nbsp;1. MOPAC gromacs mdreun error (vidhya sankar) (vidhya
 sankar)<br>&gt;&nbsp; &nbsp; <br>&gt; <br>&gt; Message: 1<br>&gt; Date: Mon, 6 Dec 2010 18:34:00 +0530 (IST)<br>&gt; From: vidhya sankar &lt;<a ymailto="mailto:scvsankar_agr@yahoo.com" href="/mc/compose?to=scvsankar_agr@yahoo.com">scvsankar_agr@yahoo.com</a>&gt;<br>&gt; Subject: [gmx-users] MOPAC gromacs mdreun error (vidhya sankar)<br>&gt; To: <a ymailto="mailto:gmx-users@gromacs.org" href="/mc/compose?to=gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>&gt; Message-ID: &lt;<a ymailto="mailto:665887.73042.qm@web95501.mail.in.yahoo.com" href="/mc/compose?to=665887.73042.qm@web95501.mail.in.yahoo.com">665887.73042.qm@web95501.mail.in.yahoo.com</a>&gt;<br>&gt; Content-Type: text/plain; charset="iso-8859-1"<br>&gt; <br>&gt; Dear Gerrit sir,<br>&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;&nbsp;&nbsp;Thank you for your atonce reply <br>&gt; my mopac stand alone binary works well&nbsp; i have tested using ./run_mopac7 tests/force<br>&gt; but mdrun
 is not properly linked aganist libmopac.a <br>&gt;&nbsp; i am using the following command to&nbsp; configure <br>&gt; ./configure --disable-float --without-qmmm-mopac CPPFLAGS=-DUSE_MOPAC LIBS=-lmopac LDFLAGS=-L/usr/local/lib<br>&gt; if i use --with-qmmm-mopac i end in&nbsp; make error (compilation&nbsp; error) <br>&gt;&nbsp; also how to check how to check the link of mdrun aganist libmopac.a&nbsp; is there is any command for this?<br>&gt; i am expecting your precious reply thanks in advance<br>&gt; <br>&gt; <br>&gt; <br>&gt; <br>&gt; <br>&gt; <br>&gt; <br>&gt; <br>&gt; <br>&gt; <br>&gt; <br>&gt; <br>&gt; <br>&gt; <br>&gt; <br>&gt; --- On Mon, 6/12/10, Gerrit Groenhof &lt;<a ymailto="mailto:ggroenh@gwdg.de" href="/mc/compose?to=ggroenh@gwdg.de">ggroenh@gwdg.de</a>&gt; wrote:<br>&gt; <br>&gt; From: Gerrit Groenhof &lt;<a ymailto="mailto:ggroenh@gwdg.de" href="/mc/compose?to=ggroenh@gwdg.de">ggroenh@gwdg.de</a>&gt;<br>&gt; Subject: Subject: [gmx-users] :
 MOPAC gromacs mdreun error (vidhya sankar)<br>&gt; To: <a ymailto="mailto:gmx-users@gromacs.org" href="/mc/compose?to=gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>&gt; Date: Monday, 6 December, 2010, 12:57 PM<br>&gt; <br>&gt; Please check if your stand-alone mopac binary works properly. Also check if the mdrun binary is linked against the libmopac.a<br>&gt; <br>&gt; Gerrit<br>&gt;&gt; Message: 2<br>&gt;&gt; Date: Mon, 6 Dec 2010 11:55:58 +0530 (IST)<br>&gt;&gt; From: vidhya sankar&lt;<a ymailto="mailto:scvsankar_agr@yahoo.com" href="/mc/compose?to=scvsankar_agr@yahoo.com">scvsankar_agr@yahoo.com</a>&gt;<br>&gt;&gt; Subject: [gmx-users] : MOPAC gromacs mdreun error (vidhya sankar)<br>&gt;&gt; To: <a ymailto="mailto:gmx-users@gromacs.org" href="/mc/compose?to=gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>&gt;&gt; Message-ID:&lt;<a ymailto="mailto:157886.66733.qm@web95506.mail.in.yahoo.com"
 href="/mc/compose?to=157886.66733.qm@web95506.mail.in.yahoo.com">157886.66733.qm@web95506.mail.in.yahoo.com</a>&gt;<br>&gt;&gt; Content-Type: text/plain; charset="iso-8859-1"<br>&gt;&gt; <br>&gt;&gt; Dear Gerrit&nbsp; sir,<br>&gt;&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Thank you for your previous reply i tried as u said with&nbsp; -nt 1 option in mdrun but still i have got hte following error when i run mdrun_d<br>&gt;&gt; QM/MM calculation requested.<br>&gt;&gt; there we go!<br>&gt;&gt; Layer 0<br>&gt;&gt; nr of QM atoms 2<br>&gt;&gt; QMlevel: PM3/STO-3G<br>&gt;&gt; Program mdrun_d, VERSION 4.5.3<br>&gt;&gt; Source code file: qmmm.c, line: 697<br>&gt;&gt; Fatal error:<br>&gt;&gt; Semi-empirical QM only supported with Mopac.<br>&gt;&gt;&nbsp; &nbsp; but i configured gromacs properly with mopac. how to solve the problem<br>&gt;&gt;&nbsp; &nbsp; is ther is any link that i have not made available ? i am expecting your worthfull
 reply thanks isn advance<br>&gt;&gt; <br>&gt;&gt; <br>&gt; <br>&gt; -- gmx-users mailing list&nbsp; &nbsp; <a ymailto="mailto:gmx-users@gromacs.org" href="/mc/compose?to=gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>&gt; <a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>&gt; 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>&gt; Please don't post (un)subscribe requests to the list. Use the www interface or send it to <a ymailto="mailto:gmx-users-request@gromacs.org" href="/mc/compose?to=gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>&gt; Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>&gt; <br>&gt; <br>&gt;
 -------------- next part --------------<br>&gt; An HTML attachment was scrubbed...<br>&gt; URL: <a href="http://lists.gromacs.org/pipermail/gmx-users/attachments/20101206/a0dfd1a9/attachment-0001.html" target="_blank">http://lists.gromacs.org/pipermail/gmx-users/attachments/20101206/a0dfd1a9/attachment-0001.html</a><br>&gt; <br>&gt; ------------------------------<br>&gt; <br>&gt; Message: 2<br>&gt; Date: Mon, 06 Dec 2010 18:40:36 +0100<br>&gt; From: Sebastian Fritsch &lt;<a ymailto="mailto:fritsch@mpip-mainz.mpg.de" href="/mc/compose?to=fritsch@mpip-mainz.mpg.de">fritsch@mpip-mainz.mpg.de</a>&gt;<br>&gt; Subject: Re: [gmx-users] Constraining two virtual sites<br>&gt; To: <a ymailto="mailto:gmx-users@gromacs.org" href="/mc/compose?to=gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>&gt; Message-ID: &lt;<a ymailto="mailto:4CFD2014.1080802@mpip-mainz.mpg.de"
 href="/mc/compose?to=4CFD2014.1080802@mpip-mainz.mpg.de">4CFD2014.1080802@mpip-mainz.mpg.de</a>&gt;<br>&gt; Content-Type: text/plain; charset=ISO-8859-1; format=flowed<br>&gt; <br>&gt; Hi again,<br>&gt; <br>&gt; <br>&gt; unfortunately I would need the constraint for a large number of molecules, say ~1000 i.e. 2000 vsites representing a few atoms each.<br>&gt; The constraint length should however be the same for each molecule. <br>&gt; <br>&gt; Assigning groups for each constraint group would not be possible due to the group limit I guess.<br>&gt; The do_constraint function in pull.c does not look too complicated, do you think it would be reasonable to adapt it a bit and call it<br>&gt; from some function like calc_bonds which comes after vsite construct and before vsite spread? The constraint forces in my case should be<br>&gt; be completely determined by the vsites, or am I wrong? Is there an easier solution?<br>&gt; <br>&gt; <br>&gt; Thank
 you,<br>&gt; Sebastian<br>&gt; <br>&gt;&gt; Hi,<br>&gt;&gt; <br>&gt;&gt; Unfortunately constraints can not be applied to virtual sites and grompp apparently<br>&gt;&gt; does not not check for this. I will add a check.<br>&gt;&gt; Constraints between virtual sites can lead to very complex constraint equations<br>&gt;&gt; between the masses involved. Thus the general case if difficult to implement.<br>&gt;&gt; Your cases seems quite simple and could be done with the pull code,<br>&gt;&gt; if it would support more than one reference group, which I was thinking to implement.<br>&gt;&gt; <br>&gt;&gt; How many molecules do you have?<br>&gt;&gt; <br>&gt;&gt; Berk<br>&gt;&gt; <br>&gt; <br>&gt; <br>&gt; <br>&gt; ------------------------------<br>&gt; <br>&gt; Message: 3<br>&gt; Date: Mon, 06 Dec 2010 13:56:04 -0500<br>&gt; From: <a ymailto="mailto:nishap.patel@utoronto.ca" href="/mc/compose?to=nishap.patel@utoronto.ca">nishap.patel@utoronto.ca</a><br>&gt;
 Subject: [gmx-users] Energy-groups?<br>&gt; To: <a ymailto="mailto:gmx-users@gromacs.org" href="/mc/compose?to=gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>&gt; Message-ID: &lt;<a ymailto="mailto:20101206135604.dgzv1f90bk08ogw4@webmail.utoronto.ca" href="/mc/compose?to=20101206135604.dgzv1f90bk08ogw4@webmail.utoronto.ca">20101206135604.dgzv1f90bk08ogw4@webmail.utoronto.ca</a>&gt;<br>&gt; Content-Type: text/plain;&nbsp;&nbsp;&nbsp; charset=ISO-8859-1;&nbsp;&nbsp;&nbsp; DelSp="Yes";<br>&gt; &nbsp;&nbsp;&nbsp; format="flowed"<br>&gt; <br>&gt; Hello,<br>&gt; <br>&gt;&nbsp; &nbsp; I want to plot the interaction potential energy between my solute&nbsp; <br>&gt; and solvent. In my .mdp file I did not mention anything under&nbsp; <br>&gt; energygrps,so I am thinking it calculates the energies for the whole&nbsp; <br>&gt; system. But is there a way I can extract say for example LJ-14 term&nbsp; <br>&gt; between my solute and solvent using the same .edr
 file? Or would I&nbsp; <br>&gt; have to specify my energygrps and run the simulation again.<br>&gt; <br>&gt; Thanks.<br>&gt; <br>&gt; Nisha<br>&gt; <br>&gt; <br>&gt; <br>&gt; ------------------------------<br>&gt; <br>&gt; Message: 4<br>&gt; Date: Mon, 06 Dec 2010 13:59:23 -0500<br>&gt; From: "Justin A. Lemkul" &lt;<a ymailto="mailto:jalemkul@vt.edu" href="/mc/compose?to=jalemkul@vt.edu">jalemkul@vt.edu</a>&gt;<br>&gt; Subject: Re: [gmx-users] Energy-groups?<br>&gt; To: Discussion list for GROMACS users &lt;<a ymailto="mailto:gmx-users@gromacs.org" href="/mc/compose?to=gmx-users@gromacs.org">gmx-users@gromacs.org</a>&gt;<br>&gt; Message-ID: &lt;<a ymailto="mailto:4CFD328B.6040308@vt.edu" href="/mc/compose?to=4CFD328B.6040308@vt.edu">4CFD328B.6040308@vt.edu</a>&gt;<br>&gt; Content-Type: text/plain; charset=ISO-8859-1; format=flowed<br>&gt; <br>&gt; <br>&gt; <br>&gt; <a ymailto="mailto:nishap.patel@utoronto.ca"
 href="/mc/compose?to=nishap.patel@utoronto.ca">nishap.patel@utoronto.ca</a> wrote:<br>&gt;&gt; Hello,<br>&gt;&gt; <br>&gt;&gt;&nbsp;&nbsp;&nbsp;I want to plot the interaction potential energy between my solute and <br>&gt;&gt; solvent. In my .mdp file I did not mention anything under energygrps,so <br>&gt;&gt; I am thinking it calculates the energies for the whole system. But is <br>&gt;&gt; there a way I can extract say for example LJ-14 term between my solute <br>&gt;&gt; and solvent using the same .edr file? Or would I have to specify my <br>&gt;&gt; energygrps and run the simulation again.<br>&gt;&gt; <br>&gt; <br>&gt; 1-4 interactions are intramolecular, so there should be no solute-solvent 1-4 <br>&gt; term.&nbsp; If you want nonbonded potentials decomposed, yes, you have to rerun your <br>&gt; trajectory (mdrun -rerun, not from scratch) with a .tpr file specifying the <br>&gt; desired groups.<br>&gt; <br>&gt; -Justin<br>&gt; <br>&gt;&gt;
 Thanks.<br>&gt;&gt; <br>&gt;&gt; Nisha<br>&gt;&gt; <br>&gt; <br>&gt; -- <br>&gt; ========================================<br>&gt; <br>&gt; Justin A. Lemkul<br>&gt; Ph.D. Candidate<br>&gt; ICTAS Doctoral Scholar<br>&gt; MILES-IGERT Trainee<br>&gt; Department of Biochemistry<br>&gt; Virginia Tech<br>&gt; Blacksburg, VA<br>&gt; jalemkul[at]vt.edu | (540) 231-9080<br>&gt; <a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br>&gt; <br>&gt; ========================================<br>&gt; <br>&gt; <br>&gt; ------------------------------<br>&gt; <br>&gt; Message: 5<br>&gt; Date: Mon, 06 Dec 2010 14:31:38 -0500<br>&gt; From: <a ymailto="mailto:nishap.patel@utoronto.ca" href="/mc/compose?to=nishap.patel@utoronto.ca">nishap.patel@utoronto.ca</a><br>&gt; Subject: Re: [gmx-users] Energy-groups?<br>&gt; To: <a ymailto="mailto:gmx-users@gromacs.org"
 href="/mc/compose?to=gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>&gt; Message-ID: &lt;<a ymailto="mailto:20101206143138.irw4yddc6c4gkscs@webmail.utoronto.ca" href="/mc/compose?to=20101206143138.irw4yddc6c4gkscs@webmail.utoronto.ca">20101206143138.irw4yddc6c4gkscs@webmail.utoronto.ca</a>&gt;<br>&gt; Content-Type: text/plain;&nbsp;&nbsp;&nbsp; charset=ISO-8859-1;&nbsp;&nbsp;&nbsp; DelSp="Yes";<br>&gt; &nbsp;&nbsp;&nbsp; format="flowed"<br>&gt; <br>&gt; Quoting "Justin A. Lemkul" &lt;<a ymailto="mailto:jalemkul@vt.edu" href="/mc/compose?to=jalemkul@vt.edu">jalemkul@vt.edu</a>&gt;:<br>&gt; <br>&gt; Thanks Justin! I do want the non-bonded potential between my solute&nbsp; <br>&gt; and solvent. So in my .mdp file I put my solute and solvent as&nbsp; <br>&gt; energygrps and ran mdrun using this command:<br>&gt;&nbsp;&nbsp;&nbsp;mdrun -s md1.tpr(including energygrps) -rerun md.xtc (my trajectory)<br>&gt; <br>&gt; Is that correct? I don't understand how
 it would be faster or any&nbsp; <br>&gt; different from running the simulation from scratch, because in my .mdp&nbsp; <br>&gt; file it still has time of 100ns. I just modified the energygrps.<br>&gt; <br>&gt; <br>&gt;&gt; <br>&gt;&gt; <br>&gt;&gt; <a ymailto="mailto:nishap.patel@utoronto.ca" href="/mc/compose?to=nishap.patel@utoronto.ca">nishap.patel@utoronto.ca</a> wrote:<br>&gt;&gt;&gt; Hello,<br>&gt;&gt;&gt; <br>&gt;&gt;&gt;&nbsp; I want to plot the interaction potential energy between my solute&nbsp; <br>&gt;&gt;&gt; and solvent. In my .mdp file I did not mention anything under&nbsp;&nbsp;&nbsp;<br>&gt;&gt;&gt; energygrps,so I am thinking it calculates the energies for the&nbsp;&nbsp;&nbsp;<br>&gt;&gt;&gt; whole system. But is there a way I can extract say for example&nbsp;&nbsp;&nbsp;<br>&gt;&gt;&gt; LJ-14 term between my solute and solvent using the same .edr file?&nbsp;&nbsp;&nbsp;<br>&gt;&gt;&gt; Or would I have to specify my energygrps and run
 the simulation&nbsp;&nbsp;&nbsp;<br>&gt;&gt;&gt; again.<br>&gt;&gt;&gt; <br>&gt;&gt; <br>&gt;&gt; 1-4 interactions are intramolecular, so there should be no<br>&gt;&gt; solute-solvent 1-4 term.&nbsp; If you want nonbonded potentials decomposed,<br>&gt;&gt; yes, you have to rerun your trajectory (mdrun -rerun, not from scratch)<br>&gt;&gt; with a .tpr file specifying the desired groups.<br>&gt;&gt; <br>&gt;&gt; -Justin<br>&gt;&gt; <br>&gt;&gt;&gt; Thanks.<br>&gt;&gt;&gt; <br>&gt;&gt;&gt; Nisha<br>&gt;&gt;&gt; <br>&gt;&gt; <br>&gt;&gt; -- <br>&gt;&gt; ========================================<br>&gt;&gt; <br>&gt;&gt; Justin A. Lemkul<br>&gt;&gt; Ph.D. Candidate<br>&gt;&gt; ICTAS Doctoral Scholar<br>&gt;&gt; MILES-IGERT Trainee<br>&gt;&gt; Department of Biochemistry<br>&gt;&gt; Virginia Tech<br>&gt;&gt; Blacksburg, VA<br>&gt;&gt; jalemkul[at]vt.edu | (540) 231-9080<br>&gt;&gt; <a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin"
 target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br>&gt;&gt; <br>&gt;&gt; ========================================<br>&gt;&gt; -- <br>&gt;&gt; gmx-users mailing list&nbsp; &nbsp; <a ymailto="mailto:gmx-users@gromacs.org" href="/mc/compose?to=gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>&gt;&gt; <a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>&gt;&gt; Please search the archive at<br>&gt;&gt; <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>&gt;&gt; Please don't post (un)subscribe requests to the list. Use the www<br>&gt;&gt; interface or send it to <a ymailto="mailto:gmx-users-request@gromacs.org" href="/mc/compose?to=gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>&gt;&gt; Can't post? Read <a
 href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>&gt; <br>&gt; <br>&gt; <br>&gt; <br>&gt; <br>&gt; ------------------------------<br>&gt; <br>&gt; Message: 6<br>&gt; Date: Mon, 06 Dec 2010 14:52:05 -0500<br>&gt; From: "Justin A. Lemkul" &lt;<a ymailto="mailto:jalemkul@vt.edu" href="/mc/compose?to=jalemkul@vt.edu">jalemkul@vt.edu</a>&gt;<br>&gt; Subject: Re: [gmx-users] Energy-groups?<br>&gt; To: Discussion list for GROMACS users &lt;<a ymailto="mailto:gmx-users@gromacs.org" href="/mc/compose?to=gmx-users@gromacs.org">gmx-users@gromacs.org</a>&gt;<br>&gt; Message-ID: &lt;<a ymailto="mailto:4CFD3EE5.6020001@vt.edu" href="/mc/compose?to=4CFD3EE5.6020001@vt.edu">4CFD3EE5.6020001@vt.edu</a>&gt;<br>&gt; Content-Type: text/plain; charset=ISO-8859-1; format=flowed<br>&gt; <br>&gt; <br>&gt; <br>&gt; <a ymailto="mailto:nishap.patel@utoronto.ca"
 href="/mc/compose?to=nishap.patel@utoronto.ca">nishap.patel@utoronto.ca</a> wrote:<br>&gt;&gt; Quoting "Justin A. Lemkul" &lt;<a ymailto="mailto:jalemkul@vt.edu" href="/mc/compose?to=jalemkul@vt.edu">jalemkul@vt.edu</a>&gt;:<br>&gt;&gt; <br>&gt;&gt; Thanks Justin! I do want the non-bonded potential between my solute and <br>&gt;&gt; solvent. So in my .mdp file I put my solute and solvent as energygrps <br>&gt;&gt; and ran mdrun using this command:<br>&gt;&gt;&nbsp; mdrun -s md1.tpr(including energygrps) -rerun md.xtc (my trajectory)<br>&gt;&gt; <br>&gt;&gt; Is that correct? I don't understand how it would be faster or any <br>&gt;&gt; different from running the simulation from scratch, because in my .mdp <br>&gt;&gt; file it still has time of 100ns. I just modified the energygrps.<br>&gt;&gt; <br>&gt; <br>&gt; It should be significantly faster.&nbsp; In this case, mdrun is not doing the <br>&gt; integration, it's simply using the known positions to
 re-calculate energies.<br>&gt; <br>&gt; -Justin<br>&gt; <br>&gt;&gt; <br>&gt;&gt;&gt; <br>&gt;&gt;&gt; <br>&gt;&gt;&gt; <a ymailto="mailto:nishap.patel@utoronto.ca" href="/mc/compose?to=nishap.patel@utoronto.ca">nishap.patel@utoronto.ca</a> wrote:<br>&gt;&gt;&gt;&gt; Hello,<br>&gt;&gt;&gt;&gt; <br>&gt;&gt;&gt;&gt;&nbsp; I want to plot the interaction potential energy between my solute <br>&gt;&gt;&gt;&gt; and solvent. In my .mdp file I did not mention anything under&nbsp; <br>&gt;&gt;&gt;&gt; energygrps,so I am thinking it calculates the energies for the&nbsp; whole <br>&gt;&gt;&gt;&gt; system. But is there a way I can extract say for example&nbsp; LJ-14 term <br>&gt;&gt;&gt;&gt; between my solute and solvent using the same .edr file?&nbsp; Or would I <br>&gt;&gt;&gt;&gt; have to specify my energygrps and run the simulation&nbsp; again.<br>&gt;&gt;&gt;&gt; <br>&gt;&gt;&gt; <br>&gt;&gt;&gt; 1-4 interactions are intramolecular, so there should be
 no<br>&gt;&gt;&gt; solute-solvent 1-4 term.&nbsp; If you want nonbonded potentials decomposed,<br>&gt;&gt;&gt; yes, you have to rerun your trajectory (mdrun -rerun, not from scratch)<br>&gt;&gt;&gt; with a .tpr file specifying the desired groups.<br>&gt;&gt;&gt; <br>&gt;&gt;&gt; -Justin<br>&gt;&gt;&gt; <br>&gt;&gt;&gt;&gt; Thanks.<br>&gt;&gt;&gt;&gt; <br>&gt;&gt;&gt;&gt; Nisha<br>&gt;&gt;&gt;&gt; <br>&gt;&gt;&gt; <br>&gt;&gt;&gt; -- <br>&gt;&gt;&gt; ========================================<br>&gt;&gt;&gt; <br>&gt;&gt;&gt; Justin A. Lemkul<br>&gt;&gt;&gt; Ph.D. Candidate<br>&gt;&gt;&gt; ICTAS Doctoral Scholar<br>&gt;&gt;&gt; MILES-IGERT Trainee<br>&gt;&gt;&gt; Department of Biochemistry<br>&gt;&gt;&gt; Virginia Tech<br>&gt;&gt;&gt; Blacksburg, VA<br>&gt;&gt;&gt; jalemkul[at]vt.edu | (540) 231-9080<br>&gt;&gt;&gt; <a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin"
 target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br>&gt;&gt;&gt; <br>&gt;&gt;&gt; ========================================<br>&gt;&gt;&gt; -- <br>&gt;&gt;&gt; gmx-users mailing list&nbsp; &nbsp; <a ymailto="mailto:gmx-users@gromacs.org" href="/mc/compose?to=gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>&gt;&gt;&gt; <a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>&gt;&gt;&gt; Please search the archive at<br>&gt;&gt;&gt; <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>&gt;&gt;&gt; Please don't post (un)subscribe requests to the list. Use the www<br>&gt;&gt;&gt; interface or send it to <a ymailto="mailto:gmx-users-request@gromacs.org"
 href="/mc/compose?to=gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>&gt;&gt;&gt; Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>&gt;&gt; <br>&gt;&gt; <br>&gt;&gt; <br>&gt; <br>&gt; -- <br>&gt; ========================================<br>&gt; <br>&gt; Justin A. Lemkul<br>&gt; Ph.D. Candidate<br>&gt; ICTAS Doctoral Scholar<br>&gt; MILES-IGERT Trainee<br>&gt; Department of Biochemistry<br>&gt; Virginia Tech<br>&gt; Blacksburg, VA<br>&gt; jalemkul[at]vt.edu | (540) 231-9080<br>&gt; <a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br>&gt; <br>&gt; ========================================<br>&gt; <br>&gt; <br>&gt; ------------------------------<br>&gt; <br>&gt; -- <br>&gt; gmx-users mailing list<br>&gt; <a ymailto="mailto:gmx-users@gromacs.org"
 href="/mc/compose?to=gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>&gt; <a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>&gt; 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>&gt; <br>&gt; End of gmx-users Digest, Vol 80, Issue 41<br>&gt; *****************************************<br><br>--<br>Gerrit Groenhof<br>MPI biophysical chemistry<br>Goettingen<br>Germany<br><a href="http://wwwuser.gwdg.de/%7Eggroenh/" target="_blank">http://wwwuser.gwdg.de/~ggroenh/</a><br><br>--<br>gmx-users mailing list&nbsp; &nbsp; <a ymailto="mailto:gmx-users@gromacs.org" href="/mc/compose?to=gmx-users@gromacs.org">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<br>www interface or send it to <a ymailto="mailto:gmx-users-request@gromacs.org" href="/mc/compose?to=gmx-users-request@gromacs.org">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></blockquote></td></tr></table><br>