Please note the following regarding Gaussian98 and
partial charge generation for Amber uses:
The FINAL set of molecular properties in a G98 output
is the one to use for charge generation.
If you run a single point, that is the only one.
If you run a geomertry opt, Gaussian does a property
analysis on the INITIAL geometry and then again on
the FINAL geometry.
The reason I'm bringing this up is that the script
distributed with Amber7 and on the Web site for a
standalone RESP setup assumes a SINGLE point G98 run.
On the other hand, Antechamber is programed to look for
the last property set and so it doesn't matter if you
feed it a single point or optmization run.
To clarify the Antechamber G98 input in the sample run:
#HF/6-31G* SCF=tight Test Pop=MK iop(6/33=2) iop(6/42=6) opt
HF/6-31G* calculation type and basis set
SCF=tight normally G98 does optimizations at a
reduced convergence level (a single
point run is "tight" by default)
Test extra output
Pop=MK Do a Merz-Kollman ESP fit to get partial charges
iop(6/33=2) Output the ESP data to the standard output file
for Gaussian.
iop(6/42=6) calculate 6 points per unit area for the ESP, the
default is 1 (I have no idea what the unit area
actually is)
opt do a geometry vs. energy optimization of the input
coordinates.
--
----------------------------------------------------------------------------
James W. Caldwell (voice) 415-476-8603
Department of Pharmaceutical Chemistry (fax) 415-502-1411
Mail Stop 0446 (email) caldwell_at_heimdal.ucsf.edu
513 Parnassus Avenue
University of California
San Francisco, CA 94143-0446
----------------------------------------------------------------------------
Received on Mon Sep 30 2002 - 15:27:37 PDT
