Dear Sir,<br>
           I want to calculate virial, hence pressure. <br>
PV=NKT-1/3 * &lt; sum of (r * f)&gt;<br>
          where r is
displacement and f is internal force.  Now, how can I proceed?<br>
1. g_traj -f name.trr -af gives me force and each velocity with time? <br>
Now, question, is the force given in data is only due to internal force or total force on each atom? <br>
Am I right if I use this force and position of the atom to calculate virial.<br>
<br>
2. If I am wrong, then please tell me how I can calculate Virial, and hence pressure. <br>
3. I saw that g_energy -f name.edr -vis gives pressure, virial and all
the stuffs. I want to know that how I can calculate it myself.<br>
4. Can you tell me how I can calculate pressure change with z - axis where I am simulating Bilayer Membrane.<br>
Thanking You,<br>
Anirban<br>
   <br>
         <br>
<br>
         <br>