On Thu, Jul 16, 2009, Hemant Kumar wrote:
> I want to calculate free energy change during the process of insertion of
> a molecules inside a pore.So I am planning to to umbrella sampling by
> defining distance between COM of pore and molecules as the reaction
> coordinate.
> Now I am unable to figure out ,how to define this distance as my reaction
> coordinate in my input file(DISANG).To be more specific,I have problem in
> understaning iat(n) variable in &rst namelist.
To have a coordinate that is the center of geometry of several atoms, set
iat1<0 and use the igr1 to specify the atoms involved. See section 6.11.5
in the manual for an example of this.
Don't use too many atoms in igr1, or things may become unstable. You may have
to just use some selected atoms of the pore to compute its COM.
....dac
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Jul 16 2009 - 18:08:25 PDT
