[gmx-developers] Energy not conserved in GROMACS replicate

#CHEN WENJIE# CHEN0819 at e.ntu.edu.sg
Tue May 10 09:50:47 CEST 2016


Hi Berk,

Sorry for misusing this mailing list.
I really appreciate if you or any other member is willing to help me, please contact me at chen0819 at e.ntu.edu.sg<mailto:chen0819 at e.ntu.edu.sg>, because it's an urgent problem for my research project.
Thank you so much!

Best regards,
Wenjie


On 10 May 2016, at 1:42 PM, Berk Hess <hess at kth.se<mailto:hess at kth.se>> wrote:

Hi,

The total energy will never be conserved perfectly. How good it is conserved depends on the time step.
But it seems, although it's not clear from your mail, that you are talking about your own program. This is a mailing list for problems with Gromacs relatead problems, not for help with your own code. All I can add is that most bugs in an MD code will cause to total energy to increase.

Cheers,

Berk

On 05/10/2016 05:57 AM, #CHEN WENJIE# wrote:
Hi GROMACS team,

I am an undergraduate student from Nanyang Technological University. I am working on my MD physics engine project that replicates some of the GROMACS features. The programming language is Java. The project is to simulate some small molecules in atomic level. The topology file is generated by GROMACS with pdb2gmx program. The force field I choose is AMBER03, and  force field parameters are also from GROMACS. I use Velocity Verlet integrator.
The problem I encounter is that when I run the simulation, the total energy doesn't conserve. It increases a little bit every frame and eventually becomes huge and the molecule deforms.

Attached is my GRO file and TOP file, and some log output files for the simulation of water molecule. One only involves bond stretching energy. One only involves angle rotation energy. One involves both bond and angle energy. For each frame, I log the kinetic energy, potential energy, total energy and position of each atom. Although when Ep increases, Ek will decrease, and vice versa, the total energy increases a little bit each frame.

I really want to know if there is any part I miss out or didn't do it correctly. I really appreciate it if you can give me some advice. Thank you in advice.

Best regards,
Wenjie




--
Gromacs Developers mailing list

* Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-developers_List before posting!

* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-developers or send a mail to gmx-developers-request at gromacs.org<mailto:gmx-developers-request at gromacs.org>.

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-developers/attachments/20160510/c2298609/attachment-0001.html>


More information about the gromacs.org_gmx-developers mailing list