Mixed Basis Sets with Gaussian

Obtaining a Basis Set at EMSL Basis Set Exchange

Many basis sets for atoms can be obtained from the EMSL Basis Set Exchange website. Open this page now to use for the remainder of this tutorial.

Imagine a fictitious molecule that contains C, H, O, and Pt. You decide to use a 6-31+G* basis set. On the EMSL page, find this basis set in the list located on the left and select it. You can type the name in to narrow the search results. You should see the following:

The elements on the periodic table that this basis set exists for have a gold triangle in the bottom left. Notice how C, H, and O have basis functions defined by the 6-31+G* basis set, however, none of the transition metals including Pt do. Select the C, H, and O elements (the selected boxes turn green) and in the “Format:” drop-down menu, select “Gaussian94”. The “Optimized General Contractions” can be checked. Finally, click the “Get Basis Set” button. A new window will appear.

Before moving on, deselect all atoms and select the Pt atom only. The list on the left automatically omits any basis set that does not exist for Pt. For this tutorial, find the “modified LANL2DZ” basis set for Pt and select it in the list. Finally, click “Get Basis Set” after ensuring that the Format is set to “Gaussian94”.

Dissecting the Basis Set File

Building a Gaussian Input File with a Mixed Basis Set

Click here to see an example Gaussian input file implementing our mixed basis set. The contents of this file is discussed below.

Lines 1-13 contain the normal Gaussian input syntax. Note the route line (line 4) contains the keyword GENECP instead of the name of a basis set. GENECP was chosen because 1.) we have decided to implement a mixed basis set and 2.) an ECP is being included (for Pt). If no ECP is being used, GEN is the appropriate keyword.

Below the atom specification block is a blank line (line 13). Following this is the section that defines the basis functions for all the atoms. The order does not matter. Here I’ve chosen to list the basis functions for H first, followed by C, O, and Pt. Each set of basis functions are separated by four asterisks. I simply pasted in the appropriate lines from the EMSL Basis Set files obtained earlier.

The basis function section must always be terminated by a blank line (line 89).

Since the modified LANL2DZ basis set contains an ECP for Pt, we must now include this in the ECP section that follows the basis function section (lines 90-126). Again, I have pasted in the bottom part of the basis set output from EMSL. Terminate the section with a single blank line (line 127).

You now have implemented a mixed basis set where there are basis functions for each atom given in the input file.