AMBER Archive (2005)

Subject: RE: AMBER: Antechamber/Divcon assigned charges

From: Junmei Wang (jwang_at_encysive.com)
Date: Thu Mar 03 2005 - 23:20:38 CST


In the tutorial for simulating sustiva, the charges of sustiva are am1-bcc charges, not Mulliken or CM1 or CM2. AM1-bcc, which mimics HF/6-31G* RESP charges, is obtained by performing bond charge correction on AM1 Mulliken charges.

Best

Junmei

I have a question regarding Antechamber and Divcon.

If one looks at the tutorial for simulating Sustiva using the General Amber
Force Field, the Divcon program is called to optimize the geometry then
assign Partial Mulliken, Partial CM1, and Partial CM2 charges (at least this
is what I gather from the Divcon.out file). However; when looking at the
charges that end up in the Sustiva.prepin file they are different from the
ones in the Divcon.out file.

For instance the charges on the chlorine from divcon.out

ATOMIC CHARGES:

   ATOM ELEMENTAL PARTIAL PARTIAL PARTIAL
    NO. SYMBOL MULLIKEN CM1 CM2
                     CHARGE CHARGE CHARGE

30 CL 0.00328 -0.02671 -0.02622

and the charge from the sustiva.prepin (-0.07512)

8 Cl1 cl E 7 5 4 1.742 119.676 -179.857 -0.07512

Sorry if this is a naïve question, but how are the charges that end up in
the prepin file being assigned - as a function of the ones calculated in
divcon or another source?

I have searched the reflector, DivconLite manual, and AMBER 8.0 manual to no
avail.

Thanks in advance...

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu