> None of the current Amber force fields have a special H-bond term. As you
> noticed, they treat this as electrostatic and van der Waals. Aside from
> writing new code, there is nothing that you can switch on.
You could develop a new ff using explicit hbonding terms, or as Dave
suggests (I believe) add code to collect the electrostatic and van
der Waals for designated interactions and report this separately.
The latter can also be done using ANAL, creating a group for donors
and a group for acceptors.
Generally what people want to know about hbonding is to what extent
a given hbond exists, which carnal and probably ptraj can determine
from the mdcrd output.
Bill Ross
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Wed Jan 14 2004 - 15:53:08 PST
