AMBER: About the Amber constraint constant

From: Yuhui Cheng <>
Date: Sun, 15 Aug 2004 01:12:50 -0700 (PDT)

     When Sander set the harmonic constraint using NMR metnod, Dr. Alan Grossfield mentioned in his WHAM program that the constraint energy should be "E=k*(x-x0)^2" instead of "E=(1/2)*k*(x-x0)^2". I checked the fortran source code "_disnrg.f" in the $AMBERHOME/src/sander/ and found that the comments line use the formula "E=(1/2)*k*(x-x0)^2" but the program line neglected the "1/2". Is my understanding correct?

Do you Yahoo!?
New and Improved Yahoo! Mail - 100MB free storage!
The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" to
Received on Mon Aug 16 2004 - 09:53:01 PDT
Custom Search