<br><br>----- Original Message -----<br>From: sapna sarupria &lt;sapna.sarupria@gmail.com&gt;<br>Date: Thursday, September 16, 2010 1:35<br>Subject: [gmx-users] Wrtiing pair interaction energy on the fly<br>To: Discussion list for GROMACS users &lt;gmx-users@gromacs.org&gt;<br><br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>Dear all<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font><br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>I understand that pair interactions in gromacs are calculated through assembly loops. Can some one guide me to the resources using which I can write out the interaction energy/force for each pair directly from gromacs during the simulation (without using the rerun option). The rerun option is especially difficult when you want pair interaction between all atoms of one's system since one has to define that many index groups. I most cases rerun is good, but in my case I have to test interaction between several pairs. <br><br>You are limited to 256 energy groups with or without using -rerun. If you need more, then perhaps you should consider some other simulation software that is better suited to this task. With GROMACS, you could use an environment variable to force the use of the generic C kernels, and inside that kernel use printf to write out the energy for each interaction. This will be horribly slow, and will be vastly simpler if not attempted in parallel.<br><br>Mark