<br>My student and I have been trying the latest Gromacs version (3.3.2) for the simulation of a CNT. We noticed that the bonded energies calculated with 3.3.2 are quite different compared to those from 3.3.1 (see resulting output below). Has anyone observed similar results? Was anything changed in 
3.3.2 in how the bonded potentials were calculated? I tried another case for a system with harmonic bond and angle potentials, and the bonded energies were identical for both versions, which makes me believe that the Morse and G96Angle potentials are somehow different in 
3.3.2. Could this be a bug?<br><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; --Amadeu<br><br><br>Initial set of energy obtained from the same TPR file. Both output were obtained with the single precision versions. <br><br><span style="font-family: courier new,monospace;">
<u>Output from Gromacs-3.3.2:</u><br>&nbsp;&nbsp; Energies (kJ/mol)</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Morse&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; G96Angle&nbsp;&nbsp;&nbsp; Proper Dih.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; LJ (SR)&nbsp; Disper. corr.
</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">&nbsp;&nbsp;&nbsp; 5.78994e+03&nbsp;&nbsp;&nbsp; 1.44076e+04&nbsp;&nbsp;&nbsp; 5.47403e+02&nbsp;&nbsp; -1.57601e+02&nbsp;&nbsp; -2.46982e+01</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">&nbsp;&nbsp; Coulomb (SR)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Potential&nbsp;&nbsp;&nbsp; Kinetic En.&nbsp;&nbsp; Total Energy&nbsp;&nbsp;&nbsp; Temperature</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">
&nbsp;&nbsp;&nbsp; 0.00000e+00&nbsp;&nbsp;&nbsp; 2.05627e+04&nbsp;&nbsp;&nbsp; 3.32197e+01&nbsp;&nbsp;&nbsp; 2.05959e+04&nbsp;&nbsp;&nbsp; 1.11447e+01</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">&nbsp;Pressure (bar)</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">&nbsp;&nbsp; -3.06280e+02<br><br><br><u>Output from Gromacs-3.3.1:</u><br>&nbsp;&nbsp; Energies (kJ/mol)<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Morse&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; G96Angle&nbsp;&nbsp;&nbsp; Proper Dih.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; LJ (SR)&nbsp; Disper. corr.<br>&nbsp;&nbsp;&nbsp; 
4.38162e+01&nbsp;&nbsp;&nbsp; 4.15264e+01&nbsp;&nbsp;&nbsp; 5.53398e+02&nbsp;&nbsp; -1.57601e+02&nbsp;&nbsp; -2.46982e+01<br>&nbsp;&nbsp; Coulomb (SR)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Potential&nbsp;&nbsp;&nbsp; Kinetic En.&nbsp;&nbsp; Total Energy&nbsp;&nbsp;&nbsp; Temperature<br>&nbsp;&nbsp;&nbsp; 0.00000e+00&nbsp;&nbsp;&nbsp; 4.56441e+02&nbsp;&nbsp;&nbsp; 9.54848e+00&nbsp;&nbsp;&nbsp; 4.65990e+02&nbsp;&nbsp;&nbsp; 3.20338e+00
<br>&nbsp;Pressure (bar)<br>&nbsp;&nbsp; -5.04502e+02<br><br style="font-family: courier new,monospace;"></span><br>