# Stability for long NVE run

**URL:** <https://gromacs.bioexcel.eu/t/stability-for-long-nve-run/10948>\
**Category:** User discussions\
**Tags:** mdp-parameters, simulation-setup\
**Created:** [January 6, 2025, 9:47am UTC](https://gromacs.bioexcel.eu/t/stability-for-long-nve-run/10948 "2025-01-06T09:47:16Z")\
**Posts on this page:** 3\
**Page:** 1

<div class="post-metadata">

**Author:** ![nitinmalapally](https://avatars.discourse-cdn.com/v4/letter/n/a9adbd/32.png) [@nitinmalapally](https://gromacs.bioexcel.eu/u/nitinmalapally)\
**Post date:** [January 6, 2025, 9:47am UTC](https://gromacs.bioexcel.eu/t/stability-for-long-nve-run/10948/1 "2025-01-06T09:47:16Z")

</div>

GROMACS version: 2024-dev-20241213-e782eff84f-dirty-unknown  
GROMACS modification: Yes/No

I’ve implemented the Path MD algorithm (MD in trajectory space) in a GROMACS copy based on the above version. It’s branching off from the commit that integrates the multi-simulation feature.

Now, I’m testing the implementation for correctness i.e., proper equilibration in NVT and total energy conservation in NVE, using the system acetone in explicit water (TIP3P) with 50k atoms.

I’ve run energy minimization to a state with energy -11013.0 and maximum force 0.7455. To equilibrate, I’m using the `sd` integrator running 500,000 steps with 5 step annealing control and final temperature 150. For the NVE run, I’m using the `md` integrator running 5,000,000 steps. For both runs, I’m _not_ using DD.

Now, the NVE run seems to be unstable, crashing at step 2,615,753 with LINCS warnings. The graph looks like this:

 ![image](https://europe1.discourse-cdn.com/flex017/uploads/bioexcel1/original/2X/c/cc1d8081a29fb1cbf836b18748789d356c425073.png)

I would be very grateful for any tips/wisdoms on how to prepare for and successfully run long NVE simulations such that the energy is conserved.

Note: I’ll gladly provide more (specific) information on request.

Thanks in advance,  
Nitin Malapally

---

<div class="post-metadata">

**Author:** ![obZehn](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/obzehn/32/2629_2.png) [@obZehn](https://gromacs.bioexcel.eu/u/obZehn)\
**Post date:** [January 6, 2025, 4:25pm UTC](https://gromacs.bioexcel.eu/t/stability-for-long-nve-run/10948/2 "2025-01-06T16:25:03Z")

</div>

Dear @nitinmalapally

Is the same behaviour present in a clean GROMACS installation (i.e. without your modification)?

---

<div class="post-metadata">

**Author:** ![nitinmalapally](https://avatars.discourse-cdn.com/v4/letter/n/a9adbd/32.png) [@nitinmalapally](https://gromacs.bioexcel.eu/u/nitinmalapally)\
**Post date:** [January 6, 2025, 4:49pm UTC](https://gromacs.bioexcel.eu/t/stability-for-long-nve-run/10948/3 "2025-01-06T16:49:37Z")

</div>

@obZehn I’m testing the behaviour of the new feature itself. In theory, it must be similar to standard MD, so I guess I’m asking for tips for preparing long NVE simulations.