From: Christopher Hartshorn (cmhartshorn_at_wsu.edu)
Date: Mon Feb 23 2009 - 12:53:15 CST
Hi everyone. After running an FEP calculation (see parameters below
message),I am obtaining an FEP graph (dG vs Lambda step) which seems
odd. I write to this forum because I think it is my system (or
parameters) and not an actual physical result. Essentially I am
seeing, as lambda goes from 0 to 1, the dG drop below zero, hit a
minimum (around -120), and then rises again to around lambda step
0.98. This is what I expect. But at around lambda step 0.98, it
reverts going negative again (e.g it rises from -120 to about -44, at
lambda=0.98, then reverts and ends around -50 at lambda=1). Is this
one of those end point anomalies? Does my config. file seem okay (e.g
using Langevin coupling to pressure/temp baths, etc.)? I'm using a
CVS build from Dec 2, 2008 is this a good build for FEP? etc.
Thank you to all for any help you can offer.
Chris Hartshorn
WSU
structure dmpc_mid1.psf
coordinates mdsim2ns_out.coor
velocities mdsim2ns_out.vel
extendedSystem mdsim2ns_out.xsc
set outputname dmpcmid1FEPGA155ksteps04_out
# firsttimestep 1000000
#############################################################
## SIMULATION PARAMETERS ##
#############################################################
# Input parameters
paraTypeCharmm on
parameters par_all27_prot_lipid_sil.inp
parameters par-extraterms.inp
# Force-Field Parameters
exclude scaled1-4
1-4scaling 1.0
cutoff 10.0
switching on
switchdist 8.5
pairlistdist 11.5
# Integrator Parameters
timestep 1.0 ;# 1fs/step
nonbondedFreq 1
fullElectFrequency 2
stepspercycle 20
margin 1
wrapWater on
wrapNearest off
# PME (for full-system periodic electrostatics) only with PBC's
# Must get cell dimensions from PBC basisvectors for accurate PME
# Remember multiples of 2,3, and/or 5
PME on
PMEGridSizeX 120
PMEGridSizeY 120
PMEGridSizeZ 100
# Constant Pressure Control (variable volume) only for PBC's
useGroupPressure yes # smaller fluctuations
useFlexibleCell no # allow dimensions to fluctuate independently
# Constant Temperature Control
langevin on ;# do langevin dynamics
langevinDamping 10 ;# damping coefficient (gamma) of 10/ps
langevinTemp 310
langevinHydrogen off ;# don't couple langevin bath to hydrogens
# Nos´e-Hoover Langevin piston pressure control
langevinPiston on
langevinPistonTarget 1.01325 ;# in bar -> 1 atm
langevinPistonPeriod 200.0
langevinPistonDecay 100.0
langevinPistonTemp 310
rigidbonds all
rigidtolerance 0.000001
rigiditerations 400
fixedAtoms on
fixedAtomsForces on
fixedAtomsFile fixQ011.pdb
fixedAtomsCol B
#############################################################
## EXECUTION SCRIPT ##
#############################################################
# FEP
source fep.tcl
fep on
fepFile dmpc_mid1_FEPGA1.fep
fepCol B
fepOutFile DMPCmid1GA155ksteps04.fepout
fepOutFreq 10
# LOOP OVER LAMBDA-STATES
fepEquilSteps 5000
set nSteps 110000
set dLambda 0.04
set init {0.0 0.00001 0.0001 0.001 0.01 0.02}
set end {0.98 0.99 0.999 0.9999 0.99999 1.0}
runFEPlist $init $nSteps
runFEP 0.02 0.98 $dLambda $nSteps
runFEPlist $end $nSteps
This archive was generated by hypermail 2.1.6 : Wed Feb 29 2012 - 15:52:24 CST