Hi, <br><br>I have a 10 ns job. But when using the g_energy to calculate the component of the energy, I got a positive total-energy!! <br><br>Who can give some suggestions?&nbsp; Thanks a lot!<br><br><br>The NTP ensemble is used. Before the final run, it has been run for 2ns to get the equilibrium state. The system has about 10000 atoms.
<br>As shown in the following, the rmsd of the total energy and density are not quite big.<br><br>Energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Average&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; RMSD&nbsp;&nbsp;&nbsp;&nbsp; Fluct.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Drift&nbsp; Tot-Drift<br>-------------------------------------------------------------------------------
<br>Bond&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 8121.16&nbsp;&nbsp;&nbsp; 129.568&nbsp;&nbsp;&nbsp; 129.564 0.000320592&nbsp;&nbsp;&nbsp; 3.20592<br>Angle&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 16284.9&nbsp;&nbsp;&nbsp; 174.598&nbsp;&nbsp;&nbsp; 174.578 -0.00090818&nbsp;&nbsp;&nbsp; -9.0818<br>Ryckaert-Bell.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 4489.18&nbsp;&nbsp;&nbsp; 147.842&nbsp;&nbsp;&nbsp; 
144.019 -0.0115717&nbsp;&nbsp; -115.717<br>LJ-14&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1169.97&nbsp;&nbsp;&nbsp; 34.0818&nbsp;&nbsp;&nbsp; 33.9909 -0.000861911&nbsp;&nbsp; -8.61912<br>Coulomb-14&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2666.03&nbsp;&nbsp;&nbsp; 191.396&nbsp;&nbsp;&nbsp; 186.964 -0.0141849&nbsp;&nbsp; -141.849<br>LJ-(SR)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -
32381.8&nbsp;&nbsp;&nbsp; 146.511&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 146.4 -0.0019765&nbsp;&nbsp;&nbsp; -19.765<br>Disper.-corr.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -2462.41&nbsp;&nbsp;&nbsp; 5.29279&nbsp;&nbsp;&nbsp; 5.28651&nbsp; 8.934e-05&nbsp;&nbsp;&nbsp;&nbsp; 0.8934<br>Coulomb-(SR)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -11046.1&nbsp;&nbsp;&nbsp; 213.609&nbsp;&nbsp;&nbsp; 210.484&nbsp; 0.0126112&nbsp;&nbsp;&nbsp; 126.112<br>
Coul.-recip.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -4664.24&nbsp;&nbsp;&nbsp;&nbsp; 141.26&nbsp;&nbsp;&nbsp; 136.108&nbsp; 0.0130948&nbsp;&nbsp;&nbsp; 130.948<br>Potential&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -17823.3&nbsp;&nbsp;&nbsp; 250.663&nbsp;&nbsp;&nbsp; 250.472 -0.00338716&nbsp;&nbsp; -33.8716<br>Kinetic-En.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 32678.4&nbsp;&nbsp;&nbsp; 188.545&nbsp;&nbsp;&nbsp; 188.545
 -6.01562e-06 -0.0601562<br>Total-Energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 14855.1&nbsp;&nbsp;&nbsp; 163.416&nbsp;&nbsp;&nbsp; 163.122 -0.00339318&nbsp;&nbsp; -33.9318<br>Temperature&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 299.965&nbsp;&nbsp;&nbsp; 1.73071&nbsp;&nbsp;&nbsp; 1.73071&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0<br>Pressure-(bar)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
1.02259&nbsp;&nbsp;&nbsp; 336.285&nbsp;&nbsp;&nbsp; 336.285 -5.31626e-06 -0.0531626<br>Density-(SI)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1583.31&nbsp;&nbsp;&nbsp; 3.40323&nbsp;&nbsp;&nbsp; 3.39918 -5.7445e-05&nbsp;&nbsp; -0.57445<br><br>