[gmx-users] Shift and Switch are obsolete!?

Alan Wilter Sousa da Silva alan at biof.ufrj.br
Mon Aug 19 17:18:30 CEST 2002


Hi, pardon if I go straightforward to the question.
	I'm trying to set coulomtype to Shift (or Switch).  This is my
file.mdp:
constraints              = h-bonds
constraint-algorithm     = shake
nstenergy                = 0
ns_type                  = simple
coulombtype              = Switch
rlist                    = 1.3
rcoulomb                 = 1.2
rcoulomb_switch          = 1.1
vdwtype                  = Switch
rvdw                     = 1.2
rvdw_switch              = 1.1

GMX manual recommends rlist to be 0.1 to 0.3 nm larger than rcoulomb.  But
I got it from grompp:

ERROR: rcoulomb must be >= rlist
ERROR: rvdw must be >= rlist

Please, one question more.  If I use pure cut-off,  I still do not follow
why rlist must be larger or equal to rcoulomb.  How does GMX couple with
an atom at the exactly distance of cut-off, and atoms beyond?  One cannot
have discontinuity!  So, is a shift function implicit with cut-off?

Best regards,
-----------------------
Alan Wilter S. da Silva
-----------------------
 Laboratório de Física Biológica
  Instituto de Biofísica Carlos Chagas Filho
   Universidade do Brasil/UFRJ
    Rio de Janeiro, Brasil




More information about the gromacs.org_gmx-users mailing list