it's somewhat more complex than that, so I encourage you to look at the
code. The force constant is indeed used for all atoms in the mask, but it's
applied based on the projection of forces along the path tangent vector, so
not as simple as saying it is applied to all atoms. I encourage you to read
the Bergonzo et al IJQC paper on our implementation and the articles that we
cite.
On Mon, Nov 8, 2010 at 5:41 PM, Buddhadev Maiti <chebm.langate.gsu.edu>wrote:
> Dear Amber users,
>
> Basically I would like to know whether in NEB calculations, as far as the
> force constant is concerned,
> the selected part of the different replicas can be regarded as a single
> effective "atom" or the force constant refers to
> individual atoms inside each selected region. Also, I would like to know
> what is the appropriate value of the force
> constant for a very large system such as mine (about 4000 atoms).
>
> Thanks,
> Maiti
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Nov 08 2010 - 17:00:03 PST
