[gmx-users] rlist, rcolumb, and rvdw

[Justin Lemkul]

On 1/11/17 2:57 AM, Mark Abraham wrote:
> Hi,
>
> Those ancient comments pertain only to the deprecated "group" cutoff
> scheme. You should look at the extensive documentation of both schemes in
> the current reference manuals. Particularly vdw cutoffs are baked into the
> forcefield and should not be varied without extensive testing. See
> http://www.gromacs.org/Documentation/Terminology/Force_Fields/CHARMM
>

To emphasize, this is especially true for CHARMM lipid parameters, as is the 
case here.

Repeat after me: "I will not mess around with cutoffs for lipids." :)

-Justin

> Mark
>
> On Wed, 11 Jan 2017 04:51 Mohsen Ramezanpour <ramezanpour.mohsen at gmail.com>
> wrote:
>
>> Dear gromacs users,
>>
>> Reading through mailing list I found a nice discussion on the relation
>> between rlist, rcolumb, and rvdw:
>>
>> https://www.mail-archive.com/gmx-users@gromacs.org/msg08387.html
>>
>>
>> https://mailman-1.sys.kth.se/pipermail/gromacs.org_gmx-users/2011-February/058507.html
>>
>> If I understood correctly,
>>
>> rcoulomb < rlist and rvdw < rlist is the most accurate way,
>>
>> and
>> rlist=rcoulomb=rvdw is the commonly used way for these parameters.
>>
>> I was wondering what would be the case if I want to change the cut-offs?
>>
>> For instance, to use charmm36 in Groamcs, it is recommended to use 1.0-1.2
>>
>> as the cut-off for LJ, and rlist=rcoulomb=rvdw
>>
>> However, this setting might not result in good behaviour for some
>> lipid bilayers.
>>
>> Therefore, I want to check if other cut-offs works better in gromacs.
>>
>> If I change "1.2 nm" to "A nm" for rvdw, do I need to change rcolumb and
>> rlist to A as well? i.e. rlist=rcolumb=rvdw=A?
>>
>> I am using following parameters already:
>>
>> cutoff-scheme           = Verlet
>> nstlist                       = 20
>> rlist                          = *1.2*
>> coulombtype             = *pme*
>> rcoulomb                  = *1.2*
>> vdwtype                   = Cut-off
>> vdw-modifier            = *Force-switch*
>> rvdw_switch             = *1.0*
>> rvdw                        = *1.2*
>>
>> Thanks in advance for your comments
>> Cheers
>> Mohsen
>> --
>> *Rewards work better than punishment ...*
>> --
>> Gromacs Users mailing list
>>
>> * Please search the archive at
>> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
>> posting!
>>
>> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>>
>> * For (un)subscribe requests visit
>> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
>> send a mail to gmx-users-request at gromacs.org.
>>

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

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

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