<br><br>----- Original Message -----<br>From: Thanasis Koukoulas &lt;koukoulas_th@yahoo.gr&gt;<br>Date: Tuesday, June 15, 2010 4:55<br>Subject: [gmx-users] Positive Total energy<br>To: gmx-users@gromacs.org<br><br><table border="0" cellpadding="0" cellspacing="0"><tbody><tr><td style="font: inherit;" valign="top"><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>Hello&nbsp;<div><br></div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>I have sent several emails and i am sorry. It is my last attempt to finish this simulation.</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>The molecule i have is C-C-C-C-C-C-O-C-C-O-C-C-O-H or C10H1O3. (ciej molecule)</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font><br></div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>I dont use gromacs libraries</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>the files i use are attached to the email in which i have descriped step by step the forcefield i use.&nbsp;</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>as you can see in [pairs] i have an extra non Lj 1-4 pair. i did it because the model i use has a potential between atom 10 and 14 ether oxygen and alchohol hydrogen interaction. It is described in the paper as a repulsive energy between these two atoms. and it has a formula:</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font><br></div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>Erepulsive =&nbsp;3.325789e-7/r^12</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font><br></div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>the only way to import in gromacs i think is this:</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font><br></div><div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>[ pairs ]</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>; ai aj  funct&nbsp;</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>10 14 &nbsp; &nbsp;2 0 0 0 0 3.325789e-7</div></div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font><br></div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>it is an intra-molecular only interaction. so i couldnt think another way to do it except this. as you can see i have several zeros for the charges, C6=o and C12 =&nbsp;3.325789e-7 fidgeQQ=0 for this case</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font><br></div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>the problem is i get positive total energy. i did this sim for only one molecule just to see what is wrong and i think that coulomb (SR) are wrong for sure.</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font><br></div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>so i want&nbsp;</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font><br></div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>coulomb 1-4 /2 &nbsp; !thats why i put fudgeQQ = 0.5</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font><br></div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>the values i get at the step 0 are</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font><br></div><div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>Energies (kJ/mol)</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>&nbsp;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;Angle &nbsp; &nbsp;Proper Dih. &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;LJ-14 &nbsp; &nbsp; Coulomb-14 &nbsp; &nbsp; &nbsp; &nbsp;LJ (SR)</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>&nbsp;&nbsp; &nbsp;2.30340e+03 &nbsp; &nbsp;1.08417e+00 &nbsp;  &nbsp;3.16570e+00 &nbsp; &nbsp;8.79204e+04 &nbsp; -2.54164e+03</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>&nbsp;&nbsp;Disper. corr. &nbsp; Coulomb (SR) &nbsp; Coul. recip. &nbsp; &nbsp; &nbsp;Potential &nbsp; &nbsp;Kinetic En.</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>&nbsp;&nbsp; -2.51250e+00 &nbsp; -4.81042e+04 &nbsp; -2.84931e+04 &nbsp; &nbsp;1.10867e+04 &nbsp; &nbsp;1.84852e+04</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>&nbsp;&nbsp; Total Energy &nbsp; &nbsp;Temperature Pressure (bar) &nbsp;Cons. rmsd ()</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>&nbsp;&nbsp; &nbsp;2.95719e+04 &nbsp; &nbsp;2.99528e+02 &nbsp; -3.03583e+02 &nbsp; &nbsp;3.13477e-06<br><br></div></div></td></tr></tbody></table>OK, now that you've simplified, it seems reasonably clear that Coulomb-14 is badly broken. My leading theory is that grompp is somehow reading (part of) the text of your C12 value into a charge field.<br><br>I would test with this [pairs] line commented out, with zero in place of your intended C12 value, and with your C12 value expressed as 0.00000325789. If that doesn't shed light on the nature of the problem, please report the Step 0 energies (as above) in each case.<br><br>Then, use<br><br>gmxdump -s topol.tpr | less<br><br>and search for 3.257 and maybe attach 50 lines either side of the line that looks like it matches that value in a parameter entry. Don't attach the whole file, it'll be huge.<br><br>As a side point, the contents of your #included files seems to achieve a correct result, although your organization of the content with respect to the filenames is badly at odds with the correct approach, and will give you problems later. Compare with a simple tutorial example.<br><br>Mark