AMBER: antechamber and gaff for a homodimer

From: rebeca <>
Date: Thu, 27 Nov 2008 12:12:23 +0100

I am trying to get the parameters for an organic molecule forming a dimer
with itself (by hydrogen bonds). I am using ANTECHAMBER and GAFF. When I
use one monomer only, there is no problem, everything is OK, but when I try
to get the parameters of the dimer as a whole molecule I have problems.

What I do is a HF Gaussian calculation and then

antechamber -i gaussian.log -fi gout -o dimer.prepin -fo prepi -c resp -s 2

This gives me a dimer.prepin file, only with errors of the type "Info: the
atom number exceeds the MAXATOM, reallocate memory automatically", which
are not relevant, as I have seen in the mailing list.

When I look at the NEWPDB.PDB file I get, I see that only consist in one
monomer. The second monomer appears filled by ceros in the coordinates and
the residue number is also 1, instead of two. This is, it is not

Is it not possible to use antechamber and gaff for a molecule consisting in
two units not-covalently joined?

Thank you very much for your help.

Rebeca García Fandiño

The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
Received on Fri Dec 05 2008 - 17:58:17 PST
Custom Search