[gmx-users] VDW Cut-off vs Switch with CHARMM force field

James Krieger jkrieger at mrc-lmb.cam.ac.uk
Sun Aug 25 01:31:33 CEST 2013


Here is a relevant previous discussion 
http://lists.gromacs.org/pipermail/gmx-users/2013-May/081425.html

Best wishes
James

On 25 Aug 2013, at 00:28, James Krieger <jkrieger at mrc-lmb.cam.ac.uk> wrote:

> CHARMM is generally used with vdw switching between 1.0 and 1.2 nm - I presume it must be more accurate that way. There are a number of discussions on this already. Justin recommended using rlistlong 1.4 as well I think. You might want to use a short range cutoff (rcoulomb) of 1.2 for PME as well. That is what was done in the Bjelkmar et al paper where CHARMM27 was implemented in GROMACS (http://pubs.acs.org/doi/full/10.1021/ct900549r).
> 
> Best wishes
> James
> 
> On 24 Aug 2013, at 23:08, Gianluca Interlandi <gianluca at u.washington.edu> wrote:
> 
>> Hi!
>> 
>> I am running explicit solvent simulations with the CHARMM force field implemented in gromacs 4.6.3. For VDW I have been using:
>> 
>> rlist    = 1.0
>> vdw-type = Cut-off
>> rvdw     = 1.0
>> 
>> in conjunction with PME and rcoulomb=1 for electrostatics.
>> 
>> I wonder whether it would be more accurate to use instead a switching function for VDW between 0.8 and 1.0 and a larger rlist:
>> 
>> rlist       = 1.2
>> vdw-type    = Switch
>> rvdw_switch = 0.8
>> rvdw        = 1.0
>> 
>> My question is: does switching between 0.8 and 1.0 generate less artifacts than abrupt truncation at 1.0 nm for VDW?
>> 
>> Thanks!
>> 
>>    Gianluca
>> 
>> -----------------------------------------------------
>> Gianluca Interlandi, PhD gianluca at u.washington.edu
>>                   http://artemide.bioeng.washington.edu/
>> 
>> Research Scientist at the Department of Bioengineering
>> at the University of Washington, Seattle WA U.S.A.
>> -----------------------------------------------------
>> -- 
>> gmx-users mailing list    gmx-users at gromacs.org
>> http://lists.gromacs.org/mailman/listinfo/gmx-users
>> * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
>> * Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-request at gromacs.org.
>> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
> -- 
> gmx-users mailing list    gmx-users at gromacs.org
> http://lists.gromacs.org/mailman/listinfo/gmx-users
> * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
> * Please don't post (un)subscribe requests to the list. Use the 
> www interface or send it to gmx-users-request at gromacs.org.
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists



More information about the gromacs.org_gmx-users mailing list