[gmx-users] Lincs warning in free energy calculation

Sadaf Rani sadafrani6 at gmail.com
Tue Mar 24 12:30:21 CET 2020


Dear Gromacs users
I ran an MD simulation for 3ns and from the last coordinates started a free
energy calculation for 3ns. My system is well equilibrated but after 710725
steps I am getting Lincs warning every time as below:-

WARNING: There are no atom pairs for dispersion correction
starting mdrun 'GROtesk MACabre and Sinister in water'
1500000 steps,   3000.0 ps.

Step 710726, time 1421.45 (ps)  LINCS WARNING
relative constraint deviation after LINCS:
rms 0.000003, max 0.000023 (between atoms 5470 and 5472)
bonds that rotated more than 30 degrees:
 atom 1 atom 2  angle  previous, current, constraint length
   3437   3438   43.8    0.1080   0.1080      0.1080

Step 710726, time 1421.45 (ps)  LINCS WARNING
relative constraint deviation after LINCS:
rms 0.000003, max 0.000021 (between atoms 5470 and 5472)
bonds that rotated more than 30 degrees:
 atom 1 atom 2  angle  previous, current, constraint length
   3437   3438   43.1    0.1080   0.1080      0.1080

I have added restraints between atoms of ligand and protein and getting
Lincs warning between 1 atom included in restraints as below:-

[ bonds ]
;    i     j  type     r0A     r1A     r2A    fcA    r0B     r1B     r2B
 fcB
  *3437*  7908    10     0.418   0.418   10.0   0.0    0.418   0.418   10.0
  41840.00

[ angle_restraints ]
;   ai    aj    ak    al  type    thA      fcA    multA  thB      fcB
 multB
  3437  7908  7906  7908       1    100.29   0.0    1     100.29   418.40
1
  7908  3437  3439  *3437*       1    146.50   0.0    1     146.50   418.40
  1

[ dihedral_restraints ]
;   ai    aj    ak    al   type    phiA     dphiA  fcA    phiB      dphiB
 fcB
  7906  7908  *3437*  3439     1    156.16   0.0    0.0    156.16    0.0
 418.40
  7907  7906  7908  *3437*     1    -56.93   0.0    0.0    -56.93    0.0
 418.40
  7909  *3437*  3430  3435     1    131.32   0.0    0.0    131.32     0.0
418.40

Can anyone please suggest to me how Should I fix this error? Any help will
be really appreciated.

Thanks.

Sadaf


More information about the gromacs.org_gmx-users mailing list