[gmx-users] About Physical parameters

vidhya sankar scvsankar_agr at yahoo.com
Tue Jun 4 14:22:23 CEST 2013



Dear Justin Thank you for your Previuos reply

                                                                                I am using gromos53a6 ff When i changed the parameters for cut-off  (r list ) value to 1.2 

I have got Error as follows  What is the Meaning of Note 2 & 3


NOTE 2 [file cntcycpepfull2.mdp]:
  The switch/shift interaction settings are just for compatibility; you
  will get betterperformance from applying potential modifiers to your
  interactions!

NOTE 3 [file cntcycpepfull2.mdp]:
  For energy conservation with switch/shift potentials, rlist should be 0.1
  to 0.3 nm larger than rcoulomb.
NOTE 4 [file cntcycpepfull2.mdp]:
  For energy conservation with switch/shift potentials, rlist should be 0.1
  to 0.3 nm larger than rvdw.

WARNING 1 [file cntcycpepfull2.mdp]:
  The sum of the two largest charge group radii (0.518710) is larger than
  rlist (1.200000) - rvdw (1.200000)

WARNING 2 [file cntcycpepfull2.mdp]:
  The sum of the two largest charge group radii (0.518710) is larger than
  rlist (1.200000) - rcoulomb (1.200000)

nstlist        = 5 

rlist        = 1.2
rcoulomb    = 1.2 

vdwtype         = Shift
rvdw        = 1.2
coulombtype    =  Reaction-Field-zero 

pme_order    = 4
fourierspacing    = 0.16
pcoupl        = Parrinello-Rahman 

pcoupltype    = isotropic   

thus .grompp terminated due to warnings

To Avoid his I have used the rlist=1.8 (in Above valuue) Because the difference betwen rlist-rcoloum should Be greater than 0.518710 

But you maild me the cut-off looks Bizarre .it should be based on parent force field(mine is gromos53a6)

I hope  your  kind Suugestion are useful to get succesfull production mdrun
Thanks In ADVANCE


More information about the gromacs.org_gmx-users mailing list