<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;">hello again<br><br>I had a problem with a sim of some molecules in vacuum and a gmx-user told me to do a simulation only in one molecule just to see what gromacs calculates and what does not. My problem concerned the sum of all energies in the system. After i did the sim with one molecule i found out that either gromacs either me calculate the Coulomb SR(short range) energy wrong. so for step 0 of this simulation I had this result:<br><br>Started mdrun on node 0 Fri Jun 11 12:52:33 2010<br><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Step&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Time&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Lambda<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
 0.00000<br><br>Grid: 3 x 3 x 3 cells<br>Long Range LJ corr.: &lt;C6&gt; 3.8522e-03<br>Long Range LJ corr.: Epot -0.00732091, Pres: -0.00379895, Vir: 0.00732091<br>&nbsp;&nbsp; Energies (kJ/mol)<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Angle Ryckaert-Bell.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; LJ-14&nbsp;&nbsp;&nbsp;&nbsp; Coulomb-14&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; LJ (SR)<br>&nbsp;&nbsp;&nbsp; 5.70627e+00&nbsp;&nbsp;&nbsp; 8.37083e+00&nbsp;&nbsp;&nbsp; 6.80805e-03&nbsp;&nbsp;&nbsp; 1.72679e+02&nbsp;&nbsp; -3.25801e+00<br>&nbsp; Disper. corr.&nbsp;&nbsp; Coulomb (SR)&nbsp;&nbsp; Coul. recip.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Potential&nbsp;&nbsp;&nbsp; Kinetic En.<br>&nbsp;&nbsp; -7.32091e-03&nbsp;&nbsp; -2.93587e+01&nbsp;&nbsp; -1.21273e+02&nbsp;&nbsp;&nbsp; 3.28661e+01&nbsp;&nbsp;&nbsp; 3.94838e+01<br>&nbsp;&nbsp; Total Energy&nbsp;&nbsp;&nbsp; Temperature Pressure (bar)&nbsp; Cons. rmsd ()<br>&nbsp;&nbsp;&nbsp;
 7.23499e+01&nbsp;&nbsp;&nbsp; 3.65291e+02&nbsp;&nbsp;&nbsp; 4.25496e+01&nbsp;&nbsp;&nbsp; 3.53180e-07<br><br>so as you can see the value of&nbsp; Coulomb (SR) is<br>-2.93587e+01 and my calculation is <br>-1302.68<br><br>I believe that even i have done wrong calculations the value that gromacs finds is not correct. <br>I found this value from the formula sum(f*qi*qj/rij) for every atom pair like<br>1-2,1-3,(not the 1-4),1-5....etc 2-3,2-4(not 1-2 because it is already calculated) ....etc ...3-.....etc......... inside the rcoulomb (the whole molecule)<br>where f is the electric conversion factor f=138,935kJmol-1*nm*e-2<br>IS COULOMB (SR) calculated that way by hand??? i an using PME btw<br><br>i have not done arithmetic mistake i am sure for that (at least :) )<br>Please consider that:<br>1) I dont use ffXXX.itp libraries from gromacs i have my own<br>2) I have tried every possible way to do this simulation for the past 3 months (its my problem i know
 that !!! :) )<br>3) I have calculated the coulomb 1-4 with success and great accuracy<br>4) I have also caclulated the LJ (SR) with success (irrelevant!!!) :)<br>5) i am using&nbsp; gromacs 4.05<br><br>Please if you have smth to advise me about this simulation except committing suicide&nbsp; of course :) please feel free to do it<br><br>Thanks in advance<br><br>especially&nbsp; <span class="sender">Mark Abraham who helped my at first place<br></span></td></tr></table><br>