AMBER: units

From: Gustavo Pierdominici Sottile <gsottile.unq.edu.ar>
Date: Wed, 11 Aug 2004 12:03:38 -0300

Hi, I am trying to obtain parameters for a non-standar aminoacid. For bond
stretching, I have already done an energy scan with Gaussian using different
values for the distance. With this ab initio energy plot I can calculate k.
The potential energy for bonds is : k (x-x)**2.
The question is: are the units of this k Kcal/A**2?
are those for bending Kcal/rad**2?
This question is because when I try to obtain k for a bendind C1-N-C2 , the
result is 157 Kcal/rad**2 . This number is higher respect to other
similar bendings and this make me think I am doing something wrong.

pd : C1 is sp2, N is sp2, and C2 is sp3
Thanking in advance
Gustavo


-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Wed Aug 11 2004 - 16:53:00 PDT
Custom Search