<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><br></div><div>Dear all,&nbsp;</div><div><br></div><div>When I try to simulation with NVE ensemble, the&nbsp;</div><div>total energy keeps decreasing. Geometry optimization&nbsp;</div><div>and solvent equilibration are done&nbsp;before the NVE simulation.&nbsp;</div><div><br></div><div><div>I follow the requirements provided in&nbsp;</div><div><font class="Apple-style-span" color="#167F7F"><a href="http://www.gromacs.org/Documentation/Terminology/NVE">http://www.gromacs.org/Documentation/Terminology/NVE</a></font></div><div><font class="Apple-style-span" color="#167F7F"><br></font></div></div><div>Would any one help me to figure out what was wrong with</div><div>my settings?</div><div><br></div><div>Millions of thanks and best wishes,</div><div>Yun-an</div><div><br></div><div><font class="Apple-style-span" color="#167F7F"><br></font></div><div>BTW: Please check the following for the .mdp file</div><div><br></div><div><div>;Run control: A leap-frog algorithm for integrating Newton's equations.</div><div>integrator &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = md</div><div>;time step in femtoseconds</div><div>dt &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 0.001</div><div>;number of steps</div><div>nsteps &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 10000000</div><div>;No constraints</div><div>constraints &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= none</div><div>constraint_algorithm &nbsp; &nbsp; = shake</div><div>shake_tol &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1e-08</div><div><br></div><div>; group(s) for center of mass motion removal</div><div>comm-grps &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= AA CCl4</div><div>; frequency for center of mass motion removal</div><div>nstcomm &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &nbsp;1000</div><div>;frequency to write coordinates to output trajectory file</div><div>nstxout &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1000</div><div>;frequency to write velocities to output trajectory file</div><div>nstvout &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1000</div><div>;frequency to write energies to log file</div><div>nstlog &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 1000</div><div>;frequency to write energies to energy file</div><div>nstenergy &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1000</div><div>;group(s) to write to energy file</div><div>energygrps &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = System</div><div><br></div><div>;Frequency to update the neighbor list (and the long-range forces,</div><div>;when using twin-range cut-off's).</div><div>nstlist &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 10</div><div>ns_type &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= grid</div><div>rlist &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1.5</div><div><br></div><div>; Use periodic boundary conditions in all directions.</div><div>pbc = xyz</div><div><br></div><div>;cut-off distance for the short-range neighbor list</div><div>;treatment of electrostatic interactions</div><div>coulombtype &nbsp; &nbsp; &nbsp;= PME-switch</div><div>fourierspacing &nbsp; = 0.12</div><div>pme_order &nbsp; &nbsp; &nbsp; &nbsp;= 4</div><div>ewald_rtol &nbsp; &nbsp; &nbsp; = 1e-5</div><div>rcoulomb_switch &nbsp;= 1.0</div><div><div>rcoulomb &nbsp; &nbsp; &nbsp; &nbsp; = 1.2</div><div><br></div><div>;treatment of van der waals interactions</div><div>vdwtype &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= Switch</div><div>rvdw-switch &nbsp; &nbsp; &nbsp;= 1.0</div><div>rvdw &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 1.2</div><div>dispcorr &nbsp; &nbsp; &nbsp; &nbsp; = EnerPres</div><div><br></div><div>;Temperature coupling</div><div>tcoupl &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= no</div><div><br></div><div>;Pressure coupling</div><div>pcoupl &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = no</div><div><br></div><div>;Velocity generation</div><div>gen_vel &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= no</div><div><br></div></div><div><br></div></div></body></html>