Hi Elodie,
running a TI change in which some molecule completely disappears is no
problem. See here for an example:
Steinbrecher, T. and Mobley, D.L. and Case, D.A., J. Chem. Phys., 127,
214108, 2007
however if your ion is in anisotropic surroundings, i.e. bound to
something, you may get convergence problems towards the end (high lambdas)
of your simulation and need to introduce additional restraints. Look at
e.g. Deng, Roux, J. Chem. Theor. Comput. 2006, 2, 1255-73 on one way to do
that.
You will have to run charge removal and vdw removal as two completely
independent simulations (we're working on a way to allow both in one step
for amber11). Each of these will run with two groups, so that your
starting state is on one processor group and the end state on the other.
See Amber Tutorial A9 for practical consideration on how to run such a
simulation:
http://ambermd.org/tutorials/
Kind Regards,
Thomas
On Wed, May 13, 2009 5:13 am, Elodie Laine wrote:
> Dear all,
>
> I want to perform a FEP using the Thermal Integration functionality of
> sander, with softcore. Instead of mutating my ligand (a calcium ion
> actually), I simply want to make it disappear. Do I need to run the TI
> steps (removing charges + removing vdW) on 2 processes V0 and V1, using
> the "groupfile" option, in this case ? Or can I simply simulate one
> process ?
>
> Thanks in advance.
>
> Elodie Laine
>
> Pasteur Institute
> Structural Bioinformatics Unit
> 25-28, rue du Dr Roux
> 75724 Paris Cedex 15
> FRANCE
> +33 (0)1 44 38 94 82
>
>
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>
Dr. Thomas Steinbrecher
BioMaps Institute
Rutgers University
610 Taylor Rd.
Piscataway, NJ 08854
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed May 20 2009 - 16:21:27 PDT
