Re: Setting up softcore parameters for FEP

From: Brian Radak (bradak_at_anl.gov)
Date: Thu Aug 03 2017 - 08:20:36 CDT

Hi John,

It is perfectly possible to overlap the regions or have them be
completely concurrent. For example, the following will give you
concurrent /linear/ coupling:

alchVdwLambdaEnd 1.0
alchVdwShiftCoeff 0.0
alchElecLambdaStart 0.0

I prefer to think of the coupling coefficients as functions dependent on
alchLambda:

vdw_couple = min[1, alchLambda / alchVdwLambdaEnd]

elec_couple = max[0, (alchLambda - alchElecLambdaStart) / (1 -
alchElecLambdaStart)]

The group 1 interactions are scaled by xxx_couple(alchLambda) and the
group -1 interactions are scaled by xxx_couple(1 - alchLambda).

New in NAMD2.12 (or maybe it was 2.11) there is now a coupling for
bonded terms that is completely analogous to vdW coupling. This is off
by default because alchBondLambdaEnd is set to 0.0. You can also turn
off vdW scaling with the same trick or turn off elec scaling by setting
alchElecLambdaStart 1.0.

HTH,
Brian

On 08/02/2017 04:18 PM, John Green wrote:
> Hi,
>
> I am trying to set up soft-core parameters for an FEP simulation. This
> is how this part of input file looks like:
>
> alchVdwLambdaEnd1.0
>
> alchElecLambdaStart 0.5
>
> alchVdWShiftCoeff5.0
>
> alchDecoupleno
>
>
> Based on tutorial, my understanding is that coupling and decoupling of
> electrostatic interaction for disappearing and appearing particles
> _always_ occur in two separate periods.
>
> [ 0 (1-alchElecLambdaStart) ] & [ (1-alchElecLambdaStart) 1 ]
>
>
> While for VDW you can have coupling and decoupling at the same time.
>
> [ (1-alchVdwLambdaEnd) 1 ] & [0 (1-alchVdwLambdaEnd) ]
>
> Is there any option to do that for electrostatic interactions as well?
> The manual says the default for alchElecLambdaStart is 0.5, so even if
> I don't use it in my script it wouldn't be helpful.
>
>
> Best,
>
> John
>
>
>

-- 
Brian Radak
Postdoctoral Appointee
Leadership Computing Facility
Argonne National Laboratory
9700 South Cass Avenue, Bldg. 240
Argonne, IL 60439-4854
(630) 252-8643
brian.radak_at_anl.gov

This archive was generated by hypermail 2.1.6 : Sun Dec 31 2017 - 23:21:31 CST