On Sat, Nov 23, 2013, X Sun wrote:
>
> I am trying to fix an atom of a polymer in my MD simulation to make a
> circumstance that the polymer is immoblized on some subtrate. At the
> beginning, I added a force constant 500kcal/mol to the atom, in tutorial
> B1, it said this was a very big number. But in my simulation, I found
> the atom moved around slowly. Should I increase the force constant? or
> is there any other method to fix one atom of a polymer in amber?
500 is indeed a big number. If you "fixed" atom actually moves (more than
a tiny amount, like 0.1 Ang) then there is probably something wrong with
how you are trying to restrain it. We'd need to know details of what you did
to be of any help.
>
> Additionally, is it possible to calculate free energy of some
> comformations in Amber?
Have you consulted Chap 4 ("Sampling and free energies") in the manual?
...dac
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Sat Nov 23 2013 - 13:30:02 PST
