[gmx-users] charged group problem

Adriana Garro adrianagarrosl at gmail.com
Tue Aug 11 21:29:44 CEST 2015


Dear users,

My system contains peptide + TFE + water in an octahedron box.
I am using oplsaa.ff
I am trying to do EM but running this command

grompp -f em.mdp -c  C22_TFE-W.gro -p topol.top -o em.tpr

the next error message come up

WARNING 1 [file em.mdp]:
  The sum of the two largest charge group radii (2.045171) is larger than
  rlist (1.000000)

I am afraid of changing the em.mdp parameters because as I understand these
are the correct for my force field.
I'd really appreciate suggestions!
Thanks in advance


Adriana


here a part of the em.mdp file

nstlist        = 1
ns_type        = grid
rlist          = 1.0
coulombtype    = PME                        ; Use particle-mesh ewald
rcoulomb       = 1.0
rvdw           = 1.0

I read that large charge groups indicate that a fair number of atoms have
been included in the same charge group. Usually only two or three atoms are
in a charge group, rendering them fairly small.

Looking in the TFE.itp, it's true, the whole molecule is in the same
charged group

     1   opls_160      1    TFE   CA3F      1      0.126     12.011   ;
qtot 0.126
     2   opls_161      1    TFE   CF3F      1      0.532     12.011   ;
qtot 0.658
     3   opls_162      1    TFE   OH3F      1     -0.635    15.9994   ;
qtot 0.023
     4   opls_163      1    TFE   HO3F      1      0.429      1.008   ;
qtot 0.452
     5   opls_164      1    TFE   F13F      1     -0.206    18.9984   ;
qtot 0.246
     6   opls_164      1    TFE   F23F      1     -0.206    18.9984   ;
qtot 0.04
     7   opls_164      1    TFE   F33F      1     -0.206    18.9984   ;
qtot -0.166
     8   opls_165      1    TFE   H13F      1      0.083      1.008   ;
qtot -0.083
     9   opls_165      1    TFE   H23F      1      0.083      1.008   ;
qtot 0


***********************************************************
Dra. Adriana D. Garro
Química Medicinal
Facultad de Química, Bioquímica y Farmacia
Universidad Nacional de San Luis
IMASL-CONICET
San Luis, Argentina

Tel..............:+54 266 4424689  int 6153
e-mail...: adgarro at unsl.edu.ar
e-mail...:  adrianagarrosl at gmail.com


**********************************************************


More information about the gromacs.org_gmx-users mailing list