[gmx-users] Lincs warning in free energy calculation (FAISAL NABI)
Sadaf Rani
sadafrani6 at gmail.com
Sun May 3 02:50:13 CEST 2020
Dear Gromacs users
I am using amber99sb-ildn.ff force field for protein and I used GAFF for
ligand parameterization.
I tried different options e.g using constraint =all bonds as well as
h-bonds. Reducing the time step to 1 fs but still, I could not solve the
problem and get lincs warning production run as below:-
Step 106202, time 212.404 (ps) LINCS WARNING
relative constraint deviation after LINCS:
rms 0.035343, max 2.112628 (between atoms 3770 and 3773)
bonds that rotated more than 30 degrees:
atom 1 atom 2 angle previous, current, constraint length
3437 3438 90.0 2.2392 0.1552 0.1080
3770 3771 90.0 0.1045 0.1696 0.1090
3770 3772 90.0 2.9369 0.1212 0.1090
3770 3773 90.0 0.1043 0.3393 0.1090
Wrote pdb files with previous and current coordinates
Step 106202, time 212.404 (ps) LINCS WARNING
relative constraint deviation after LINCS:
rms 0.035038, max 2.096736 (between atoms 3770 and 3773)
bonds that rotated more than 30 degrees:
atom 1 atom 2 angle previous, current, constraint length
3437 3438 90.0 2.2392 0.1540 0.1080
3770 3771 90.0 0.1045 0.1689 0.1090
3770 3772 90.0 2.9369 0.1198 0.1090
3770 3773 90.0 0.1043 0.3375 0.1090
I would really appreciate your help in fixing this problem I am
facing since long.
Thanks.
Sadaf
More information about the gromacs.org_gmx-users
mailing list