[gmx-users] mdp options for charmm27 force field in gromacs version 4.6.7

[Justin Lemkul]

On 5/18/15 2:21 PM, Christopher Neale wrote:
> Dear Users:
>
> I would like to use the charmm27 force field in gromacs version 4.6.7 and I would like to know if it is possible to get the proper treatment of vdw interactions.
>
> Information for gromacs 5 is here:
> http://www.gromacs.org/Documentation/Terminology/Force_Fields/CHARMM
>
> However, Gromacs 4.6.7 does not contain the option for vdw-modifier = force-switch
>
>    Invalid enum 'force-switch' for variable vdw-modifier, using
>    'Potential-shift-Verlet'
>    Next time use one of: 'Potential-shift-Verlet' 'Potential-shift' 'None'
>    'Potential-switch' 'Exact-cutoff'
>
> Is a potential-shift the same as a force-switch and therefore the following options are correct for charmm27 in gromacs 4.6.7?
>
> ns_type = grid
> cutoff-scheme = verlet
> vdwtype = cut-off
> rlist = 1.2
> rvdw = 1.2
> rvdw-switch = 1.0
> rcoulomb = 1.2
> vdw-modifier = Potential-shift-Verlet
> coulombtype = PME
>
> I know this has been discussed before on the list, but I can't find anything that seems authoritative for gromacs versions prior to v5.
>

I cannot personally attest to anything aside from what's on the wiki.  We have 
verified the forces are equivalent between CHARMM's vfswitch and GROMACS' 
force-switch modifier.  This is how we do all our simulations, so for posterity, 
what I wrote on the wiki is what I must emphasize is known to be correct.

 From the manual, it seems that the group scheme + vdwtype = shift is equivalent 
to the modern force-switch modifier.  The language of potential-shift does not 
sound as if it is the same as force-switch, and I don't know why such redundancy 
would be present anyway; perhaps legacy naming reasons.

Any reason why you can't use 5.0 with the known correct settings?  The check I 
would do is a single-point energy with the settings you think are equivalent and 
want to test and print out the forces with gmx traj.

-Justin

-- 
==================================================

Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow

Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 629
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201

jalemkul at outerbanks.umaryland.edu | (410) 706-7441
http://mackerell.umaryland.edu/~jalemkul

==================================================