vivek sharma wrote:
Dear Justin,
Thanks for your reply..
But I am not getting how and where should I give the option of
energy_excl. For freezing a part of molecule ?
In the .mdp file.
what are the other checks to be done for the same ?
I don't know what you're asking.
-Justin
With Thanks,
Vivek
2008/10/14 Justin A. Lemkul <[EMAIL PROTECTED] <mailto:[EMAIL PROTECTED]>>
vivek sharma wrote:
2008/10/14 Justin A. Lemkul <[EMAIL PROTECTED]
<mailto:[EMAIL PROTECTED]> <mailto:[EMAIL PROTECTED]
<mailto:[EMAIL PROTECTED]>>>
vivek sharma wrote:
Hi there,
I tried the mdrun by keeping a group of residue position
restrained, and it is working well (verified by comparing
RMSD
plot for the same), thanks for your suggestions regarding
the same.
But, I have seen that the part of the molecule I kept for
PR is
still having some motion, how can I keep that part rigid ?
Position restraints do not necessarily guarantee that atomic
positions stay absolutely fixed. Instead, there is an energy
penalty for moving them, which could in some cases be overcome.
I am not getting any idea of energy_excl option mentioned
above,
please explain.
The freezegrps option can be used to fix atomic positions, but in
this case large forces can be generated within the frozen
group. As
such, you can apply energygrp_excl to exclude energetic terms
within
the frozen group.
Will it help me in reducing the time taken for the simulation if
I'll keep the part of molecule frozen ?
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.
-Justin
Thanks,
~Vivek
-Justin
With Thanks,
Vivek
2008/10/12 Erik Lindahl <[EMAIL PROTECTED]
<mailto:[EMAIL PROTECTED]>
<mailto:[EMAIL PROTECTED] <mailto:[EMAIL PROTECTED]>>
<mailto:[EMAIL PROTECTED] <mailto:[EMAIL PROTECTED]>
<mailto:[EMAIL PROTECTED] <mailto:[EMAIL PROTECTED]>>>>
Hi,
On Oct 10, 2008, at 12:46 PM, Justin A. Lemkul wrote:
I want to run MD over a part of my molecule ,
for few
residues only (not the whole molecule).
Can I do it using GROMACS ?
I searched for the online documentation and
mailing list,
but unable to get appropriate information.
If somebody has already tried such things
earlier, please
suggest and direct me for appropriate link and
address.
Well, if your goal is to keep certain parts fixed
and allow
others to move, probably the easiest way to do it
is to apply
position restraints to the "fixed" part.
You can also set parts of the system as a freeze
group, in which
case you can exclude all nonbonded interactions inside
the freeze
group with the energygrp_excl option in your mdp file.
The main advantage of this is of course that you will
improve
performance if 99% of your system is frozen (although all
interactions between the frozen and non-frozen parts still
have to
be calculated). On the other hand, completely freezing
part
of the
system is not very realistic, and you're likely to get
strange
behavior in the interface...
Cheers,
Erik
_______________________________________________
gmx-users mailing list gmx-users@gromacs.org
<mailto:gmx-users@gromacs.org>
<mailto:gmx-users@gromacs.org <mailto:gmx-users@gromacs.org>>
<mailto:gmx-users@gromacs.org
<mailto:gmx-users@gromacs.org> <mailto:gmx-users@gromacs.org
<mailto:gmx-users@gromacs.org>>>
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at
http://www.gromacs.org/search before
posting!
Please don't post (un)subscribe requests to the list.
Use thewww
interface or send it to [EMAIL PROTECTED]
<mailto:[EMAIL PROTECTED]>
<mailto:[EMAIL PROTECTED]
<mailto:[EMAIL PROTECTED]>>
<mailto:[EMAIL PROTECTED]
<mailto:[EMAIL PROTECTED]>
<mailto:[EMAIL PROTECTED]
<mailto:[EMAIL PROTECTED]>>>.
Can't post? Read
http://www.gromacs.org/mailing_lists/users.php
------------------------------------------------------------------------
_______________________________________________
gmx-users mailing list gmx-users@gromacs.org
<mailto:gmx-users@gromacs.org>
<mailto:gmx-users@gromacs.org <mailto:gmx-users@gromacs.org>>
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search
before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to [EMAIL PROTECTED]
<mailto:[EMAIL PROTECTED]>
<mailto:[EMAIL PROTECTED]
<mailto:[EMAIL PROTECTED]>>.
Can't post? Read
http://www.gromacs.org/mailing_lists/users.php
-- ========================================
Justin A. Lemkul
Graduate Research Assistant
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu <http://vt.edu> <http://vt.edu> | (540)
231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin
========================================
_______________________________________________
gmx-users mailing list gmx-users@gromacs.org
<mailto:gmx-users@gromacs.org>
<mailto:gmx-users@gromacs.org <mailto:gmx-users@gromacs.org>>
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before
posting!
Please don't post (un)subscribe requests to the list. Use the www
interface or send it to [EMAIL PROTECTED]
<mailto:[EMAIL PROTECTED]>
<mailto:[EMAIL PROTECTED]
<mailto:[EMAIL PROTECTED]>>.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
--
========================================
Justin A. Lemkul
Graduate Research Assistant
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu <http://vt.edu> | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin
========================================
--
========================================
Justin A. Lemkul
Graduate Research Assistant
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin
========================================
_______________________________________________
gmx-users mailing list gmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php