Hello, I wander what actual force constant is used in SMD method in case of
distance restraint (Amber 10).
What rk2 value should I define if I want to use constant k=10 kcal/mol*A2
for equation V=1/2k[(x-x0)**2]? 5 or 10?
This question is related to my observation of force value in SMD output file
that is twice than expected from formula f=rk2*(x-x0),
and to information found in mailing list:
"The force constant will be twice the values you give for rk2=rk3 because
Amber historically does not include a factor of 1/2 in constraint
energies".
Coud you help please?!
Dmitri Nilov,
Lomonosov Moscow State University
**
**
**
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Mar 18 2009 - 01:10:57 PDT
