From: Hannes Loeffler (Hannes.Loeffler_at_stfc.ac.uk)
Date: Tue Jan 10 2017 - 09:59:35 CST
On Tue, 10 Jan 2017 15:09:52 +0000
"B.W.J. Irwin" <bwji2_at_cam.ac.uk> wrote:
> This .pdf is useful (sorry if you already know this)
> Particularly the heading:
> Why a soft-core potential should always be included
Except that this here may very well be a case where that should not be
necessary. The dual-topology approach in NAMD means that there is
always a part vanishing and one appearing. Softcores however are
really only effective/needed at short distances(*) but Abhishek seems
to keep both ions at the same coordinates. This means there may always
be enough excluded volume available to ensure water molecules not to
come too close to the ion(s). Actually, I would think that a
dual-topology approach in this paticular case is not really a clever
idea because a simple linear transformation of the three non-bonded
parameters in a single-topology approach would not require the
restraint (there is already the constraint of a common coordinate set).
But that's how FEP/TI is currently implemented in NAMD(**).
(*) The code appears to have some WCA functionality included but that
does not seem to be documented. It is functional?
(**) I heard there were plans to work on this?
> If you haven't included one.
> The fepout fileshould show the contribution from each λ-window
> "Delta Delta A" and a measure of uncertainty.
Abishek said TI not FEP.
> If a particular window has a large contribution you could try
> splitting it into two smaller windows.
> On 2017-01-10 12:15, Hannes Loeffler wrote:
> > It would be useful it you kept the discussion in one thread.
> > Your attachment suggests that you plot λ vs <U> but your really need
> > <∂U/∂λ>. Is that what it is? How have you extracted this?
> > Maybe you send all relevant input files (control file, parmfile,
> > coordinates/alchfile) and output files (alchOutFile) so that
> > somebody on
> > the list can take a look.
> > On Tue, 10 Jan 2017 16:39:01 +0530
> > Abhishek Kumar Singh <bo13m1002_at_iith.ac.in> wrote:
> >> Hello,
> >> I have alchemically transformed K+ to Na+ at the center of water
> >> boxes (40, 60 and 80 A edge length respectively) using NAMD
> >> software and CHARMM36 forcefield. A harmonic spring of force
> >> constant = 10 kcal/mol.A^-2 has been applied throughout the
> >> transformation. Temp and Pressure has been kept at 310K and 1 atm
> >> respectively. Free energy of alchemical transformation has been
> >> calculated using thermodynamic integaration. The derivative has
> >> been computed as (∂G/∂λ) = ⟨ U(λ=1)- U)(λ=0) ⟩λ, where λ is the
> >> coupling parameter such that λ=1 and λ=0 are K+ and Na+ vdw
> >> parameters respectively. Free energy derivatives has been computed
> >> for 11 λvdw values ( 1.0, 0.9, 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, 0.2,
> >> 0.1 and 0.0). The area under the plot (derivate vs λ) gives me =
> >> -23.40, -23.46 and -23.4 8
> >> kcal/mol for the box sizes of 40A, 60A and 80A respectively. I
> >> am attaching alongwith the derivate vsλ values for 40A water box.
> >> The standard CHARMM36 prm file actually states that the value
> >> should be -18.6 kcal/mol.
> >> So what can be reason for this difference. Is there any problem in
> >> this procedure ?
> >> Your suggestions are most welcome.
> >> Sincerely
> >> Abhishek
This archive was generated by hypermail 2.1.6 : Sun Dec 31 2017 - 23:20:58 CST