On 2/7/13 2:13 PM, Yun Shi wrote:
Hi all,
I want to rerun a trajectory to calculate interaction energies between
each residue of protein A and protein B. In other words, I want to
calculate interaction energies for pairs A1 - B, A2-B, A3-B, A4-B...,
A399-B.
So instead of making an index file with 399 groups of each residue in
A and typing in rerun.mdp file 400 group names as
"energygrps = A1 A2 A3 A4 ... A399 B", and issuing g_energy command 399 times,
Writing such lines could probably be scripted in a loop of your favorite
programming language, but be aware that your .edr file is going to become
enormous, probably several GB, depending on the length of the trajectory. You
won't have 400 energy terms, you will have 400*400 energy terms.
is there an easier way to do this in Gromacs?
Not in Gromacs, but using scripts to call Gromacs functions and/or write lines
to files, yes.
-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
* 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
