On 5/6/13 5:58 PM, jhon michael espinosa duran wrote:
Hi all
I am doing an NVE simulation of a protein immersed in water, and I want
to keep track of the potential energy in the protein. Do you know how can I do
it?.
I tried saving the energy of the protein ( the protein as an energy group) in
the .edr file, but when I checked the file with g_energy, it does not give me
the option of potential energy only for the protein, it gives me the total
potential energy.
The potential energy of a subset of atoms is not stored in the .edr file; only
short-range nonbonded interactions can be decomposed by using energygrps. There
are a variety of reasons for doing so. The easiest way to get the quantity you
want is to use tpbconv create a .tpr file that contains only the protein and
then re-calculate energies using mdrun -rerun on the existing trajectory.
-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 gmx-users@gromacs.org
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 gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists