Re: AMBER: which atoms are H's for SHAKE?

From: Bill Ross <>
Date: Tue, 6 Apr 2004 13:11:41 -0700 (PDT)

> can somebody tell me how the SHAKE algorithm in sander determines which
> atoms actually are hydrogens?

The bonds are grouped in the prmtop according to whether an H is
involved (i.e. 2 groups).

> Does the assignment of the element "H" to a new atom type in Leap with
> addAtomTypes {{"X" "H" "sp3"}}
> mean that SHAKE will treat the new type as hydrogen?

Here is the code in unitio.h for selecting bonds w/ H:

        if (iAtomElement(aA) == HYDROGEN || iAtomElement(aD) == HYDROGEN) {

Since the 'Element' is referred to, I expect assigning the element
of "H" with addAtomTypes should do the trick.

The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" to
Received on Tue Apr 06 2004 - 21:53:00 PDT
Custom Search