# Too many LINCS warnings with Comm-mode Angular

**URL:** <https://gromacs.bioexcel.eu/t/too-many-lincs-warnings-with-comm-mode-angular/620>\
**Category:** User discussions\
**Created:** [August 4, 2020, 1:25am UTC](https://gromacs.bioexcel.eu/t/too-many-lincs-warnings-with-comm-mode-angular/620 "2020-08-04T01:25:13Z")\
**Posts on this page:** 11\
**Page:** 1

<div class="post-metadata">

**Author:** ![Neena](https://avatars.discourse-cdn.com/v4/letter/n/bc79bd/32.png) [@Neena](https://gromacs.bioexcel.eu/u/Neena)\
**Post date:** [August 4, 2020, 1:25am UTC](https://gromacs.bioexcel.eu/t/too-many-lincs-warnings-with-comm-mode-angular/620/1 "2020-08-04T01:25:13Z")

</div>

GROMACS version: 2019.1  
GROMACS modification: No  
Here post your question

I am conducting simulations of a peptide in gas phase using opls forcefield. Box size ~ 999 nm^3 and cutoffs 333.3 nm. Although I understand that angular is the correct comm-mode to be used here, I get too many LINCS warnings when I set comm-mode to angular but none when linear. Any thoughts?

---

<div class="post-metadata">

**Author:** ![cblau](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/cblau/32/5_2.png) [@cblau](https://gromacs.bioexcel.eu/u/cblau)\
**Post date:** [August 4, 2020, 12:19pm UTC](https://gromacs.bioexcel.eu/t/too-many-lincs-warnings-with-comm-mode-angular/620/2 "2020-08-04T12:19:00Z")

</div>

Hi Neena,

I would guess here that large movements are introduced when using the angular comm-mode due to rotating a large box. If you must use comm-mode angular, a first shot would be to equilibrate the system with comm-mode linear for some steps, then start the simulation with comm-mode angular.

---

<div class="post-metadata">

**Author:** ![hess](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/hess/32/416_2.png) [@hess](https://gromacs.bioexcel.eu/u/hess)\
**Post date:** [August 5, 2020, 8:36am UTC](https://gromacs.bioexcel.eu/t/too-many-lincs-warnings-with-comm-mode-angular/620/3 "2020-08-05T08:36:21Z")

</div>

This should work when the molecule is not crossing box edges. Where did your place your molecule in the box?

---

<div class="post-metadata">

**Author:** ![Neena](https://avatars.discourse-cdn.com/v4/letter/n/bc79bd/32.png) [@Neena](https://gromacs.bioexcel.eu/u/Neena)\
**Post date:** [August 5, 2020, 11:38am UTC](https://gromacs.bioexcel.eu/t/too-many-lincs-warnings-with-comm-mode-angular/620/4 "2020-08-05T11:38:31Z")

</div>

@cblau That did not work, but thank you

---

<div class="post-metadata">

**Author:** ![Neena](https://avatars.discourse-cdn.com/v4/letter/n/bc79bd/32.png) [@Neena](https://gromacs.bioexcel.eu/u/Neena)\
**Post date:** [August 5, 2020, 11:39am UTC](https://gromacs.bioexcel.eu/t/too-many-lincs-warnings-with-comm-mode-angular/620/5 "2020-08-05T11:39:16Z")

</div>

@hess using gmx editconf command, placed in the center of box -c -d 496.5 -bt cubic

---

<div class="post-metadata">

**Author:** ![hess](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/hess/32/416_2.png) [@hess](https://gromacs.bioexcel.eu/u/hess)\
**Post date:** [August 5, 2020, 11:52am UTC](https://gromacs.bioexcel.eu/t/too-many-lincs-warnings-with-comm-mode-angular/620/6 "2020-08-05T11:52:22Z")

</div>

Could you attach your tpr and md.log files so I can have look?

---

<div class="post-metadata">

**Author:** ![jalemkul](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/jalemkul/32/18_2.png) [@jalemkul](https://gromacs.bioexcel.eu/u/jalemkul)\
**Post date:** [August 5, 2020, 4:01pm UTC](https://gromacs.bioexcel.eu/t/too-many-lincs-warnings-with-comm-mode-angular/620/7 "2020-08-05T16:01:30Z")

</div>

@hess is there any plan to re-implement proper non-periodic vacuum calculations? The workarounds of large boxes seem to have lots of unexpected pitfalls. I know there is an open issue on this ([https://gitlab.com/gromacs/gromacs/-/issues/3526](https://gitlab.com/gromacs/gromacs/-/issues/3526)) but wanted to see if there is any momentum as the request has sat dormant for a few months.

---

<div class="post-metadata">

**Author:** ![Neena](https://avatars.discourse-cdn.com/v4/letter/n/bc79bd/32.png) [@Neena](https://gromacs.bioexcel.eu/u/Neena)\
**Post date:** [August 6, 2020, 2:53am UTC](https://gromacs.bioexcel.eu/t/too-many-lincs-warnings-with-comm-mode-angular/620/8 "2020-08-06T02:53:08Z")

</div>

@hess The forum does not allow to upload .tpr or .log files, instead I am attaching here my most recent md.mdp file where I just changed from angular to linear comm-mode to circumvent lincs warnings with the former mode. [md.mdp](https://gromacs.bioexcel.eu/uploads/short-url/hrrR1pOx2O3a5chp7nlhlLFTGdw.mdp) (2.0 KB)

---

<div class="post-metadata">

**Author:** ![hess](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/hess/32/416_2.png) [@hess](https://gromacs.bioexcel.eu/u/hess)\
**Post date:** [August 7, 2020, 8:31am UTC](https://gromacs.bioexcel.eu/t/too-many-lincs-warnings-with-comm-mode-angular/620/9 "2020-08-07T08:31:44Z")

</div>

I realized there could be an issue if your peptide undergoes large conformational transitions, as then the orientational fit might suddently change the orientation of the molecule. Do you have large conformational transitions?

---

<div class="post-metadata">

**Author:** ![hess](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/hess/32/416_2.png) [@hess](https://gromacs.bioexcel.eu/u/hess)\
**Post date:** [August 7, 2020, 8:33am UTC](https://gromacs.bioexcel.eu/t/too-many-lincs-warnings-with-comm-mode-angular/620/10 "2020-08-07T08:33:30Z")

</div>

I plan to re-implement non-periodic boundary conditions since a long time, but this has still not happened because other things have had higher priority. The main working is making sure the domain decomposition works properly and efficiently. But I just realized that maybe supporting this without domain decomposition would cover most use cases.

---

<div class="post-metadata">

**Author:** ![Neena](https://avatars.discourse-cdn.com/v4/letter/n/bc79bd/32.png) [@Neena](https://gromacs.bioexcel.eu/u/Neena)\
**Post date:** [August 7, 2020, 9:14pm UTC](https://gromacs.bioexcel.eu/t/too-many-lincs-warnings-with-comm-mode-angular/620/11 "2020-08-07T21:14:26Z")

</div>

I did not try with angular mode as I did not want to override lincs warnings. But yes I see significant conformational transitions with linear mode. Please see this post where I show those transitions [Rmsd shots after recentering with trjconv](https://gromacs.bioexcel.eu/t/rmsd-shots-after-recentering-with-trjconv/621)
