# Huge difference in results depending on Gromacs version

**URL:** <https://gromacs.bioexcel.eu/t/huge-difference-in-results-depending-on-gromacs-version/1392>\
**Category:** User discussions\
**Created:** [January 6, 2021, 9:33am UTC](https://gromacs.bioexcel.eu/t/huge-difference-in-results-depending-on-gromacs-version/1392 "2021-01-06T09:33:38Z")\
**Posts on this page:** 7\
**Page:** 1

<div class="post-metadata">

**Author:** ![Sasha](https://avatars.discourse-cdn.com/v4/letter/s/54ee81/32.png) [@Sasha](https://gromacs.bioexcel.eu/u/Sasha)\
**Post date:** [January 6, 2021, 9:33am UTC](https://gromacs.bioexcel.eu/t/huge-difference-in-results-depending-on-gromacs-version/1392/1 "2021-01-06T09:33:39Z")

</div>

GROMACS version: 2018.1, 2020.4  
GROMACS modification: No

Hi all,

We have a porous solid membrane positioned in the XY-plane, which is automatically solvated with ions added, using gmx solvate and genion, respectively. We then proceed with EM, dynamic relaxation (semiisotropic barostat, box size constant in XY) and a production simulation under external field along Z. The final result is the ionic flux through the membrane, as calculated using a custom perl script. The whole process is scripted, the forcefield is custom (broadly based on OPLS-AA) and is inside the local directory every time the batch is started.

Up until we installed Gromacs v2020.4, all the results have always been consistent across versions (2018x, 2019x) on various architectures. With 2020.4, the batch proceeds happily and produces a massively different result. I know that some degree of variation is a must, but this looks very problematic. Here is an [example](https://www.dropbox.com/s/u4nugnpdlqq93d0/fluxes.jpg?dl=0).

Do we have a problem? I can share the entire set of inputs, if needed (preferably privately).

Thanks!

---

<div class="post-metadata">

**Author:** ![pbauer](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/pbauer/32/13_2.png) [@pbauer](https://gromacs.bioexcel.eu/u/pbauer)\
**Post date:** [January 6, 2021, 9:55am UTC](https://gromacs.bioexcel.eu/t/huge-difference-in-results-depending-on-gromacs-version/1392/2 "2021-01-06T09:55:55Z")

</div>

Hello,

I think this is related to a bug we just got fixed for the 2020.5 point  
release that should go out today or tomorrow.

Please look here for the description and how we fixed it:

> **[Bug with electric field simulation using gromacs2020 (#3800) · Issues ·...](https://gitlab.com/gromacs/gromacs/-/issues/3800)**
>
> Dear all, I’m trying to run EF simulation of membrane channels using -1V with gromacs2020.1 but as soon as the simulation starts a pore is formed...

And here for the fix:

> **[Fix indexing bug in electric field module (!920) · Merge Requests · GROMACS /...](https://gitlab.com/gromacs/gromacs/-/merge_requests/920)**
>
> Improper refactoring in c8eca5d1 lead to the module using the wrong maximum atom...

Cheers

Paul

---

<div class="post-metadata">

**Author:** ![Sasha](https://avatars.discourse-cdn.com/v4/letter/s/54ee81/32.png) [@Sasha](https://gromacs.bioexcel.eu/u/Sasha)\
**Post date:** [January 6, 2021, 9:51pm UTC](https://gromacs.bioexcel.eu/t/huge-difference-in-results-depending-on-gromacs-version/1392/3 "2021-01-06T21:51:02Z")

</div>

Thanks Paul, I will wait for the new release. Did this bug affect anything not involving external fields? Basically, I have a set of simulation results entirely unrelated to what I described, and those used constant-force COM pulling, but no external fields. Do you think I should worry about those results? Thanks.

---

<div class="post-metadata">

**Author:** ![pbauer](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/pbauer/32/13_2.png) [@pbauer](https://gromacs.bioexcel.eu/u/pbauer)\
**Post date:** [January 7, 2021, 6:41am UTC](https://gromacs.bioexcel.eu/t/huge-difference-in-results-depending-on-gromacs-version/1392/4 "2021-01-07T06:41:22Z")

</div>

The new release is out now, so I would ask you to try it to see if the  
issue has also been fixed for you.  
Please check the release notes to make sure, but I don’t think any other  
simulations had been affected.

Cheers

Paul

---

<div class="post-metadata">

**Author:** ![Sasha](https://avatars.discourse-cdn.com/v4/letter/s/54ee81/32.png) [@Sasha](https://gromacs.bioexcel.eu/u/Sasha)\
**Post date:** [January 7, 2021, 9:41am UTC](https://gromacs.bioexcel.eu/t/huge-difference-in-results-depending-on-gromacs-version/1392/5 "2021-01-07T09:41:38Z")

</div>

I was hoping to. I just built it, but unfortunately, just like it was with 2020.4, sudo make install produces an inaccessible gromacs folder in /usr/local. I could change the permissions myself, but will defer to our sysadmin.

We never had this issue with earlier versions. Is there anything we should be aware of?

Thank you!

---

<div class="post-metadata">

**Author:** ![pbauer](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/pbauer/32/13_2.png) [@pbauer](https://gromacs.bioexcel.eu/u/pbauer)\
**Post date:** [January 7, 2021, 10:20am UTC](https://gromacs.bioexcel.eu/t/huge-difference-in-results-depending-on-gromacs-version/1392/6 "2021-01-07T10:20:20Z")

</div>

Hello,

can you open a bug for this in our issue tracker ([Issues · GROMACS / GROMACS · GitLab](https://gitlab.com/gromacs/gromacs/-/issues))?  
I haven’t heard of this issue yet, can you check with your sysadmin on what is happening there and to help us with resolving this?

Cheers

Paul

---

<div class="post-metadata">

**Author:** ![Sasha](https://avatars.discourse-cdn.com/v4/letter/s/54ee81/32.png) [@Sasha](https://gromacs.bioexcel.eu/u/Sasha)\
**Post date:** [January 7, 2021, 8:11pm UTC](https://gromacs.bioexcel.eu/t/huge-difference-in-results-depending-on-gromacs-version/1392/7 "2021-01-07T20:11:09Z")

</div>

Hi Paul,

Our sysadmin’s solution is simple: sudo find /usr/local/gromacs -type d -exec chmod 755 {} ;  
If you still want me to post on gitlab, I will.  
In the meantime, 2020.5 does _not_ have the issue I described, so it does indeed seem like you fixed it – thank you.

Sasha
