[AMBER] Actual force constant in SMD

From: Dmitri Nilov <nilovdm.gmail.com>
Date: Mon, 16 Mar 2009 18:07:09 +0000

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
Coud you help please?!

Dmitri Nilov,
Lomonosov Moscow State University
AMBER mailing list
Received on Wed Mar 18 2009 - 01:10:57 PDT
Custom Search