Pseudo

DESCRIPTION

This keyword requests that a model potential be substituted for the core electrons. The Cards option is by far its most-used mode. Gaussian supports a new effective core potential (ECP) input format (similar to that used by ExtraBasis) which is described below. When reading-in pseudopotentials, do not give them the same names as any internally-stored pseudopotentials: CEP, CHF, LANL1, LANL2, LP-31, SDD and SHC.

If used with ONIOM, the Pseudo keyword applies to all layer of the ONIOM. If you want to read in ECPs only for one ONIOM layer, then use the GenECP keyword instead.

Pseudo without any options defaults to Pseudo=Read.

OPTIONS

Read
Read pseudopotential data from the input stream. Input is described in the next subsection below. Cards is a synonym for Read.

Old
Read pseudopotential data using the old format (used by Gaussian 92 and earlier versions).

CHF
Requests the Coreless Hartree-Fock potentials. This option is normally used with the LP-31G basis sets.

SHC
Requests the SHC potentials.

LANL1
Requests the LANL1 potentials.

LANL2
Requests the LANL2 potentials.

FULL ECP INPUT FORMAT

Effective Core Potential operators are sums of products of polynomial radial functions, Gaussian radial functions and angular momentum projection operators. ECP input therefore specifies which potential to use on each atomic center, and then includes a collection of triplets of:

(coefficient, power of R, exponent)

for each potential for each term in each angular momentum of the ECP. Since only the first few angular momentum components have different terms, the potential is expressed as (1) terms for the general case, typically d or f and higher projection, and (2) the extra terms for each special angular momentum. Thus for an LP-31G potential, which includes special s and p projected terms, the input includes the general (d and higher) term, the s-d term (i.e., what to add to the general term to make the s component) and the p-d term.

All ECP input is free-format. Each block is introduced by a line containing the center numbers (from the molecule specification) and/or atomic symbols, specifying the atoms and/or atoms types to which it applies (just as for general basis set input-see the discussion of the Gen keyword). The list ends with a value of 0.

The pseudopotential for those centers/atoms follows:

Name,Max,ICore
Name of the potential, maximum angular momentum of the potential (i.e., 2 if there are special s and p projections, 3 if there are s, p, and d projections), and number of core electrons replaced by the potential. If Name matches the name of a previous potential, that potential is reused and no further input other than the terminator line (see below) is required.

For each component (I=1 to Max) of the current potential, a group of terms is read, containing the following information:

Title
A description of the block, not otherwise used.

NTerm
Number of terms in the block.

NPower,Expon,Coef
Power of R, exponent, and coefficient for each of the NTerm terms. NPower includes the R² Jacobian factor.

An example of an input file which includes a nonstandard ECP with its associated basis set is given below.

SIMPLIFIED ECP INPUT FORMAT

Gaussian adds flexibility to ECP input by allowing it to include pre-defined basis sets names. An ECP definition may be replaced by a line containing the standard keyword for a pre-defined basis set. In this case, the ECPs within the specified basis set corresponding to the specified atom type(s) will be used for that atom (see the examples).

KEYWORDS FOR STUTTGART/DRESDEN ECP INPUT

In Pseudo input, keywords for these ECPs are of the form XYn where n is the number of core electrons which are replaced by the pseudopotential and X denotes the reference system used for generating the pseudopotential (S for a single-valence-electron ion or M for a neutral atom).

Y specifies the theoretical level of the reference data: HF for Hartree-Fock, WB for Wood-Boring quasi-relativistic and DF for Dirac-Fock relativistic. For one- or two-valence electron atoms SDF is a good choice; otherwise MWB or MDF is recommended (although for small atoms or for the consideration of relativistic effects, the corresponding SHF and MHF pseudopotentials may be useful).

AVAILABILITY

The Stuttgart/Dresden ECPs are not uniformly available across the periodic table. The following table shows the availability of the various XY combinations, along with valid values for n. The Defaults columns list the equivalencies for the SDD keyword (which selects an all electron basis set through Cl and ECPs thereafter) and when IOp(3/6) is set to 6 (which selects ECPs for all elements).

				Valid values of n for given values of X and Y
Atom		Defaults		MWB	SDF	SHF	MDF	MHF
		IOp(3/6=6)	SDD keyword
1	H	D95	D95
2	He	D95	D95
3	Li	SDF2	D95
4	Be	SDF2	D95		2
5	B	MWB2	D95	2	2
6	C	MWB2	D95	2	2
7	N	MWB2	D95	2	2
8	O	MWB2	D95	2	2
9	F	MWB2	D95	2	2
10	Ne	MWB2	D95	2				2
11	Na	SDF10	6-31G		10
12	Mg	SDF10	6-31G		10
13	Al	MWB10	D95	10	10
14	Si	MWB10	D95	10	10
15	P	MWB10	D95	10	10
16	S	MWB10	D95	10	10
17	Cl	MWB10	D95	10	10
18	Ar	MWB10	6-31G	10				10
19	K	MWB10	MWB10	10	18	18
20	Ca	MWB10	MWB10	10	18	18
21	Sc	MDF10	MDF10				10	10
22	Ti	MDF10	MDF10				10	10
23	V	MDF10	MDF10				10	10
24	Cr	MDF10	MDF10				10	10
25	Mn	MDF10	MDF10				10	10
26	Fe	MDF10	MDF10				10	10
27	Co	MDF10	MDF10				10	10
28	Ni	MDF10	MDF10				10	10
29	Cu	MDF10	MDF10			28	10	10
30	Zn	MDF10	MDF10	28	28		10	10
31	Ga	MWB28	MWB28	28	28
32	Ge	MWB28	MWB28	28	28	28
33	As	MWB28	MWB28	28	28
34	Se	MWB28	MWB28	28	28
35	Br	MWB28	MWB28	28	28
36	Kr	MWB28	MWB28	28				28
37	Rb	MWB28	MWB28	28	36	36
38	Sr	MWB28	MWB28	28	36	36
39	Y	MWB28	MWB28	28				28
40	Zr	MWB28	MWB28	28				28
41	Nb	MWB28	MWB28	28				28
42	Mo	MWB28	MWB28	28				28
43	Tc	MWB28	MWB28	28				28
44	Ru	MWB28	MWB28	28				28
45	Rh	MWB28	MWB28	28				28
46	Pd	MWB28	MWB28	28				28
47	Ag	MWB28	MWB28	28		46		28
48	Cd	MWB28	MWB28	28				28
49	In	MWB46	MWB46	46	46
50	Sn	MWB46	MWB46	46	46
51	Sb	MWB46	MWB46	46	46
52	Te	MWB46	MWB46	46	46
53	I	MWB46	MWB46	46	46		46
54	Xe	MWB46	MWB46	46				46
55	Cs	MWB46	MWB46	46	54	54
56	Ba	MWB46	MWB46	46	54
57	La	MWB28	MWB28	28, 46, 47				46, 47
58	Ce	MWB28	MWB28	28, 47, 48				47, 48
59	Pr	MWB28	MWB28	28, 48, 49				48, 49
60	Nd	MWB28	MWB28	28, 49, 50				49, 50
61	Pm	MWB28	MWB28	28, 50, 51				50, 51
62	Sm	MWB28	MWB28	28, 51, 52				51, 52
63	Eu	MWB28	MWB28	28, 52, 53				52, 53
64	Gd	MWB28	MWB28	28, 53, 54				53, 54
65	Tb	MWB28	MWB28	28, 54, 55				54, 55
66	Dy	MWB28	MWB28	28, 55, 56				55, 56
67	Ho	MWB28	MWB28	28, 56, 57				56, 57
68	Er	MWB28	MWB28	28, 57, 58				57, 58
69	Tm	MWB28	MWB28	28, 58, 59				58, 59
70	Yb	MWB28	MWB28	28, 59				59
71	Lu	MWB60	MWB60	28, 60				60
72	Hf	MWB60	MWB60	60				60
73	Ta	MWB60	MWB60	60				60
74	W	MWB60	MWB60	60				60
75	Re	MWB60	MWB60	60				60
76	Os	MWB60	MWB60	60				60
77	Ir	MWB60	MWB60	60				60
78	Pt	MWB60	MWB60	60				60
79	Au	MWB60	MWB60	60		78	60	60
80	Hg	MWB60	MWB60	60, 78			60	60, 78
81	Tl	MWB78	MWB78	78				78
82	Pb	MWB78	MWB78	78				78
83	Bi	MWB78	MWB78	78				78
84	Po	MWB78	MWB78	78				78
85	At	MWB78	MWB78	78				78
86	Rn	MWB78	MWB78	78				78
89	Ac	MWB60	MWB60	60				60
90	Th	MWB60	MWB60	60				60
91	Pa	MWB60	MWB60	60				60
92	U	MWB60	MWB60	60				60
93	Np	MWB60	MWB60	60				60
94	Pu	MWB60	MWB60	60				60
95	Am	MWB60	MWB60	60				60
96	Cm	MWB60	MWB60	60				60
97	Bk	MWB60	MWB60	60				60
98	Cf	MWB60	MWB60	60				60
99	Es	MWB60	MWB60	60				60
100	Em	MWB60	MWB60	60				60
101	Md	MWB60	MWB60	60				60
102	No	MWB60	MWB60	60				60
103	Lr	MWB60	MWB60	60				60
104	Rf						92

Note: These ECPs are not available for elements 87 (Fr), 88 (Ra), and 105 and higher.

RELATED KEYWORDS

ChkBasis, ExtraBasis, Gen, GenECP

EXAMPLES

Specifying an ECP. This input file runs an RHF/LP-31G calculation on hydrogen peroxide, with the basis set and ECP data read from the input file:

# HF/Gen Pseudo=Read Test
 
Hydrogen peroxide

0,1
O
H,1,R2
O,1,R3,2,A3
H,3,R2,1,A3,2,180.,0
 
R2=0.96
R3=1.48
A3=109.47
 
General basis set input
****
 
O 0                        ECPs for the oxygen atoms.
OLP 2 2                    ECP name=OLP, applies to d & higher, replaces 2 electrons.
D component                Description for the general terms.
3                          Number of terms to follow.
1 80.0000000  -1.60000000
1 30.0000000  -0.40000000
2  1.0953760  -0.06623814
S-D projection             Corrections for projected terms (lowest angular momentum).
3
0  0.9212952   0.39552179
0 28.6481971   2.51654843
2  9.3033500  17.04478500
P-D                        Corrections for projected terms (highest angular momentum).
2
2 52.3427019  27.97790770
2 30.7220233 -16.49630500 
                           Blank line indicates end of the ECP block.

The basis set data follows the molecule specification section. The first line of the ECP data requests that a potential be read in (type 7) for atoms number 1 and 3 (the oxygen atoms). No potential is to be used for atoms 2 and 4 (the hydrogen atoms).

The second line of ECP data begins the input for the centers requiring a read-in potential: in this case, oxygen atoms. The potential on these centers is named OLP, it is a general term and applies to angular momentum 2 (d) and higher, and the potential replaces two electrons. Next comes a title for the general term (D component), and the number of components of that term (3); the individual components follow on the next 3 lines. Next come the corrections for the projected terms in two sections, lowest angular momentum first. Each section again consists of a title line, the number of terms to follow, and then the terms themselves.

Using Standard Basis Set Keywords to Specify ECPs. The following input file illustrates the use of the simplified ECP input format:

# Becke3LYP/Gen Pseudo=Read Opt Test
 
HF/6-31G(d) Opt of Cr(CO)6
 
0 1
Cr 0.0  0.0  0.0
molecule specification continues …
 
C O 0
6-31G(d)
****
Cr 0
LANL2DZ
****
 
Cr 0                                  ECP for chromium atom.
LANL2DZ                               Use the ECP in this basis set.

Last updated on: 9 Jun 2009