Transition Metals in Gaussian Summary

From: "Benjamin J. Moritz" <benjamin.moritz { *at * } asu.edu>
Subject: Transition Metals in Gaussian Summary
Date: Thu, 02 Nov 2000 09:37:35 -0700

I thank all who helped me in overcoming my problem with Gaussian 98 and transition metals basis sets. I was able to overcome the problem with the GEN keyword and the LANL2DZ and am trying other basis as well. Below is a summary of the responses I received.

****************************************************************************
If you are performing some calculations on heavy atoms you have to use
some pseudo-potential. You cannot simply use a basis set such as e.g.
6-311G**.
You can use the following statements (for example):
#B3LYP/GEN Opt Pseudo(SDD)
or
#B3LYP/GEN Opt Pseudo(Read)
In the first case G98 will use the Stuttgart pseudo-potential. The GEN
keyword specifies that you have to put, at the end of your input file, the
basis set for all atoms. I suggest you to have a look to the web site:
http://www.emsl.pnl.gov:2080/forms/basisform.html
In the second case you have to specify both the basis set (GEN) and the
pseudo-potentials.
I guess you can also have such a line command:
#B3LYP/6-31G Opt Pseudo(SDD)
In this case, the 6-31G will be used for light atoms. But you still have
to put the pseudo-potentials for heavy ones.
You can retrieve the well-formatted basis sets and pseudo potentials in
the web site I have just mentioned.
***************************************************************
For Hg (and most transition metals) you will need to use a basis set
modified with an electrostatic core potential (ECP) which basically removes
the core electrons from the atom and replaces them with an empirical
potential. G98 comes with the LANL2DZ basis and ECP built-in. If you want
to use this basis then set your input like the following:
#P RHF/GEN Pseudo=read
comments...
charge multiplicity
molecular specifications
C N 0
6-31+G(d')
****
Hg 0
LANL2DZ
****
Hg 0
LANL2DZ
The important parts are the Gen and Pseudo=read - this tells gaussian to
look at the end of the input file for the basis set (which can be different
for each atom) and the ECP. In the example above I used carbon and
nitrogen, you would put the element symbol for whatever atoms are in your
molecule there. Let me know if this works for you!
ATTACHMENT AT END OF FILE!

******************************************************************
LanL2DZ basis set should be enough as it goes up to atomic number 83 I
have used it for a number of transition metals already and had no
problems.
******************************************************************
According to the Basis Set Order Form, mercury is only available with the
following basis sets:
WTBS
SBKJC VDZ ECP
CRENBL ECP
CRENBS ECP
Stuttgart RLC ECP
Stuttgart RSC ECP
Ahlrichs Coulomb Fitting
http://www.emsl.pnl.gov:2080/forms/basisform.html
They can be used in Gaussian with the GEN keyword. I do not have any
experience with any of these basis sets.

**************************************************************
You've discovered one of the annoying things about calculations involving
elements past chlorine: not all basis sets have been "parameterized" for all
elements. Somewhere near the front of the G98 manual is a table of which
built-in basis sets are defined for which elements.
For your calculations on Hg compounds, your choices are to use one of the
built-in basis sets that is defined for that element, or to use gen and
pseudo=cards to read in a basis set and pseudopotential that is not built
into Gaussian.

****************************************************************
Benjamin J. Moritz                                       home (602) 453-9411
Dept. of Chemistry & Biochemistry                 work (480) 965-8509
Arizona State University                                Mailto:bmoritz { *at * } asu.edu
P. O. Box 871604
Tempe, AZ 85287-1604
****************************************************************

"Karma is justice without the satisfaction."

%chk=PtCO2CN2-LANL2DZ-RHF.chk
 #P RHF/Gen Pseudo=Read OPT(Calcfc) pop=full gfinput
 iop(6/7=3) Density=Current SCF=Tight Nosymm
 Optimize the Ground State of Pt(CO)2(CN)2 at RHF/(LANL2DZ,STO-3G)
 0 1
  C
  N,1,R2
  Pt,1,R3,2,A3
  C,3,R4,1,A4,2,D4,0
  C,3,R5,1,A5,4,D5,0
  C,3,R6,4,A6,1,D6,0
  N,4,R7,3,A7,1,D7,0
  O,5,R8,3,A8,1,D8,0
  O,6,R9,3,A9,1,D9,0
       Variables:
  R2=1.15108055
  R3=1.98976461
  R4=1.98976461
  R5=2.04407697
  R6=2.04407697
  R7=1.15108055
  R8=1.11956234
  R9=1.11956234
  A3=178.61942868
  A4=89.30938056
  A5=86.63735515
  A6=86.63735515
  A7=178.61942868
  A8=177.66388611
  A9=177.66388611
  D4=179.99793404
  D5=180.
  D6=180.
  D7=180.
  D8=0.
  D9=-0.00125552
 O C N 0
 STO-3G
 ****
  PT  0
 LANL2DZ
  ****
 PT  0
 LANL2DZ