Re: [AMBER] How to calculate the solvation free energy using Amber?

From: <hannes.loeffler.stfc.ac.uk>
Date: Mon, 17 Apr 2017 10:19:43 +0000

You want to set scmask to the empty string for one end state. If you
want to separate electrostatics from vdW transformation you can use
crgmask and run two simulations. You create two topology files: one
with Ni2+ and the second one without the Ni2+ (just remove the ion from
the coordinates). There should be no need for any TI/FEP specific
setup in leap.

Example mdin files:

vdW transformation:
1a)
ifsc = 1
scmask = ''
crgmask = '.Ni' ! set charges to zero

1b)
ifsc = 1
scmask = '.Ni' ! switch on the vdW interactions
crgmask = '.Ni' ! set charges to zero


electrostatic transformation.
2a)
ifsc = 0
crgmask = '.Ni'

2b)
ifsc = 0


Cheers,
Hannes.

_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Apr 17 2017 - 03:30:02 PDT
Custom Search