On 7/26/12 6:23 AM, akn wrote:
Dear Gromacs users,
Is it easily possible to calculate the number of H-bonds per molecule across
the box by using gromacs? If it is, how can I do this?
One would have to use g_select on a particular region of the box (simple
coordinate boundaries) to work with atoms in that region. The result is an
index file for each frame in the trajectory. You must then loop over all these
index groups (which may be of different size in each frame) and conduct analysis
using a for-loop in a shell script (or something similar) since there is no
automatic way to use the index files produced by g_select in this case.
See the discussion on this topic from last week in the archive.
-Justin
--
========================================
Justin A. Lemkul, Ph.D.
Research Scientist
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 [email protected]
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Only plain text messages are allowed!
* Please search the archive at
http://www.gromacs.org/Support/Mailing_Lists/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/Support/Mailing_Lists
