<div dir="ltr"><br>Dear Justin,<br>Thanks for your reply..<br>But I am not getting how and where should I give the option of energy_excl. For freezing a part of molecule ?<br>what are the other checks to be done for the same ?<br>
<br>With Thanks,<br>Vivek<br><div class="gmail_quote">2008/10/14 Justin A. Lemkul <span dir="ltr"><<a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>></span><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
<br>
vivek sharma wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
<br>
2008/10/14 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>>><div class="Ih2E3d"><br>
<br>
<br>
<br>
vivek sharma wrote:<br>
<br>
Hi there,<br>
<br>
I tried the mdrun by keeping a group of residue position<br>
restrained, and it is working well (verified by comparing RMSD<br>
plot for the same), thanks for your suggestions regarding the same.<br>
But, I have seen that the part of the molecule I kept for PR is<br>
still having some motion, how can I keep that part rigid ?<br>
<br>
<br>
Position restraints do not necessarily guarantee that atomic<br>
positions stay absolutely fixed. Instead, there is an energy<br>
penalty for moving them, which could in some cases be overcome.<br>
<br>
<br>
I am not getting any idea of energy_excl option mentioned above,<br>
please explain.<br>
<br>
<br>
The freezegrps option can be used to fix atomic positions, but in<br>
this case large forces can be generated within the frozen group. As<br>
such, you can apply energygrp_excl to exclude energetic terms within<br>
the frozen group.<br>
<br>
<br>
Will it help me in reducing the time taken for the simulation if I'll keep the part of molecule frozen ?<br>
<br>
</div></blockquote>
<br>
I've never used freezegrps for any large set of atoms, so maybe. If you've got surrounding solvent, though, the nonbonded interactions are still calculated between the frozen group and solvent, unless you turn that off too. But then it starts getting really unphysical, in my view.<br>
<br>
-Justin<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div class="Ih2E3d">
Thanks,<br>
~Vivek<br>
<br>
<br>
<br>
-Justin<br>
<br>
With Thanks,<br>
Vivek<br>
<br>
2008/10/12 Erik Lindahl <<a href="mailto:lindahl@cbr.su.se" target="_blank">lindahl@cbr.su.se</a><br></div>
<mailto:<a href="mailto:lindahl@cbr.su.se" target="_blank">lindahl@cbr.su.se</a>> <mailto:<a href="mailto:lindahl@cbr.su.se" target="_blank">lindahl@cbr.su.se</a><div><div></div><div class="Wj3C7c"><br>
<mailto:<a href="mailto:lindahl@cbr.su.se" target="_blank">lindahl@cbr.su.se</a>>>><br>
<br>
<br>
Hi,<br>
<br>
<br>
On Oct 10, 2008, at 12:46 PM, Justin A. Lemkul wrote:<br>
<br>
<br>
I want to run MD over a part of my molecule , for few<br>
residues only (not the whole molecule).<br>
Can I do it using GROMACS ?<br>
I searched for the online documentation and mailing list,<br>
but unable to get appropriate information.<br>
If somebody has already tried such things earlier, please<br>
suggest and direct me for appropriate link and address.<br>
<br>
<br>
Well, if your goal is to keep certain parts fixed and allow<br>
others to move, probably the easiest way to do it is to apply<br>
position restraints to the "fixed" part.<br>
<br>
<br>
You can also set parts of the system as a freeze group, in which<br>
case you can exclude all nonbonded interactions inside the freeze<br>
group with the energygrp_excl option in your mdp file.<br>
<br>
The main advantage of this is of course that you will improve<br>
performance if 99% of your system is frozen (although all<br>
interactions between the frozen and non-frozen parts still<br>
have to<br>
be calculated). On the other hand, completely freezing part<br>
of the<br>
system is not very realistic, and you're likely to get strange<br>
behavior in the interface...<br>
<br>
Cheers,<br>
<br>
Erik<br>
<br>
_______________________________________________<br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><br></div></div>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a> <mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>>><div class="Ih2E3d">
<br>
<br>
<a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before<br>
posting!<br>
Please don't post (un)subscribe requests to the list. Use thewww<br>
interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>><br>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>>>.<br>
<br>
<br>
Can't post? Read <a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
<br>
<br>
<br>
------------------------------------------------------------------------<br>
<br>
<br>
_______________________________________________<br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><br>
<a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a><br>
before posting!<br></div><div class="Ih2E3d">
Please don't post (un)subscribe requests to the list. Use the<br>
www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br></div><div class="Ih2E3d">
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>>.<br>
<br>
Can't post? Read <a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
<br>
<br></div><div class="Ih2E3d">
-- ========================================<br>
<br>
Justin A. Lemkul<br>
Graduate Research Assistant<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="Ih2E3d"><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><div class="Ih2E3d">
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><br>
<a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before<br>
posting!<br></div><div class="Ih2E3d">
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><div class="Ih2E3d">
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>>.<br>
<br>
Can't post? Read <a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
<br>
<br>
</div></blockquote>
<br>
-- <br><div><div></div><div class="Wj3C7c">
========================================<br>
<br>
Justin A. Lemkul<br>
Graduate Research Assistant<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>
</div></div></blockquote></div><br></div>