AMBER Archive (2005)Subject: AMBER: Dummy atoms: Changing the source code
From: Ilyas Yildirim (yildirim_at_pas.rochester.edu)
Date: Thu Mar 03 2005  18:24:03 CST
Dear David Case,
In my prev. emails, I was asking why AMBER does not let us use dummy
atoms in the initial state when k is not equal to 1 (pp.134, note #2). And
as far as I understand, the reason lies at the parameter function chosen
in Eq. # 5 (pp.133). When I plot the function, (1lambda)^k, it seems
that when we have dummy atoms at the initial state, the transformation
will be sharp. It is not as gradual as it is when we have dummy atoms at
the final state.
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

 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
