On 4/18/16 1:19 PM, Sheng Bi wrote:
Dear GMX Users I have some questions when using -rerun in mdrun to get some specific energy for my groups. Let me describe my question in this way. I have a system containing some groups which include A and B and others. My goal is to calculate the energy between A and B. So Here is my steps: first, I use gmx trjconv to get trajectory contains only A and B called A_B.xtc. Second, I grompp a new tpr which is also only contain group A and B called A_B.tpr. Third, I use mdrun -rerun A_B.xtc -s A_B.tpr -v -deffnm energy_A_B Last, I use gmx energy -f energy_A_B.edr and choose "Total Energy" to get total energy of A_B system. My question is, by this way, I can get the total energy (name it E_total), but in my opinion this energy is composed of the energy between A and A (name it E_A_A), energy between B and B (name it E_B_B), and energy between A and B (name it E_A_B). I only care about E_A_B. By now, I have to repeat above four steps to get E_A_A, and E_B_B, then use E_total minus E_A_A and E_B_B to get E_A_B. This is a very tedious work. I am not quite familiar with -rerun, Is there any ingenious method to get E_A_B ? I am not quite familiar with -rerun?
Just set: energygrps = A B and create a new .tpr, which is used for the rerun. -Justin -- ================================================== Justin A. Lemkul, Ph.D. Ruth L. Kirschstein NRSA Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 629 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 http://mackerell.umaryland.edu/~jalemkul ================================================== -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.