AMBER Archive (2005)Subject: Re: AMBER: Dummy atoms: Changing the source code
From: Ilyas Yildirim (yildirim_at_pas.rochester.edu)
Date: Fri Mar 04 2005  20:56:34 CST
Hi,
The functional form of the parameter function I want to use is
(1+tanh(k*(2x1)))/2. Here k is 3,4,5... If I use k=3 or k=4, the
function is pretty smooth. Initially it is zero and gradually increases
till it gradually becomes 1 at x=1. I dont see any reason why it should
not work. If we want to use a function which will let us use the dummy
atoms in the molecules, I cannot think of any other function than this one.
Best,
On Fri, 4 Mar 2005, Julien Michel wrote:
> Hi,
>
> > So, we were thinking that if we use a parameter function,
> > tanh[k(2lambda1)], such that the Eq. 5 becomes,
> >
> > V(lambda)= tanh[k(2lambda1)]*V0 + {1  tanh[k(2lambda1)]}*V1
> >
> > the problem might be solved. When we plot the function tanh[k(2lambda1)],
> > we see that the initial and the final transformations look gradual. So, we
> > were wondering how we can change the AMBER code such that we can use this
> > new function rather the old one. Do we have to change just the
> > thermo_int.f file, or something else too? Thanks in advance.
> >
> > Best regards,
> > Ilyas Yildirim
>
> Which values of k would you advocate then ? In principle, one of the
> potential should be zero at the end state, which you don't get if you don't
> set k to a high value. K would have to be large (more than 10) to extinguish
> the LJ of V1 at a lambda value close to 1. Furthermore, if you use this
> functional form, you will simply see that that your Hamiltonian is more or
> less constant until the region of lambda 0.5 where it is rapidly changing,
> thus your perturbation will just be harder to accomplish at lambda 0.5 instead
> of 0.0 or 1.0.
> You might want to read some papers on softcore if you are interested in
> creating/destroying atoms in a free energy simulation (see J.Chem.Phys
> ,100,12,9025 for example).
>
> Regards,
> 
> Julien Michel
> University of Southampton
> Department of Chemistry
> Dr J. W. Essex group
> 
> The AMBER Mail Reflector
> To post, send mail to amber_at_scripps.edu
> To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
>
>

Ilyas Yildirim

 Department of Chemisty  
 University of Rochester  
 Hutchison Hall, # B10  
 Rochester, NY 146270216  Ph.:(585) 275 67 66 (Office) 
 http://www.pas.rochester.edu/~yildirim/ 


The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
