Home
Basis Sets Lab Activity
Basis Set Main Page
Students
Introduction
Objectives
Background Reading
Procedure
Questions
Further Work
References/Support Materials
Teachers
Additional Background
Materials
Standards
First year chemistry curriculum concepts
Second year chemistry curriculum concepts
Readings
Overview
Atomic Orbitals
Lab Activities
Zmatrices
Basis Sets
Geometry Optimizations
Ionization Energies
Support Materials
Interactive Tools
Glossary of Terms
Quick Guide to DISCO Output File
Related Links
ChemViz
Computational Chemistry
SUCCEED's Computational Chemistry
Developers' Tools
What's New?
Discussion Board
Team Members
Email the Group
Contact Webmaster

Procedure:

For this lab, you will need access to a spreadsheet.
Part I:
Below is the equation for a Gaussian Type Orbital where alpha () is 0.4166 and r is the radius from the nucleus.
Data Entry:

At the top of the spreadsheet, label the first cell "alpha." In the next cell
define the constant.

Leaving a few rows blank, identify one cell "Radius." This is where we will
give the equation various values of r.

Beneath Radius, you are going to create a list of values from 3 to +3. Divide
the values into 0.01 increments. In order to do this, you are going to use a
"Fill" comand. First, type a 3 in the first cell beneath the "Radius"
label. Next, highlight the 3 and all the cells down to at least the 600th
line. In the menu bar, look under Edit to find the "Fill" command. You want to
choose the "Series" option. This will then prompt you for the increment by
which you want your column to increase. The default is set to 1, but you will
type 0.01. Click on OK. If your values do not reach +3, then repeat these
steps until you reach this value.

Next to the Radius column, title the next cell "Gaussian" or "GTO."

Beneath this label you will need to enter you GTO function as given above.
Remeber to identify alpha and the radius by their cell names and that pi=PI(
), e^{x}=EXP(x).
Hints:

When referring to a cell, make sure you use, $, before the column or row if it
remains constant.

The first value should be 0.00869.

Use a "Fill" command, (this time use "Down" instead of "Series") in order
to get values for all the radii.
Graphing:

Highlight the lists "Radius" and "GTO" including labels.

In the toolbar, click on the button with a wand and bar graph.

Select the area on the spreadsheet where you want the graph to appear.

Follow the steps to create a line graph.
Hints:

In Excel, make sure you choose the XY scatter plot and the second type
of line graph.
Part II:
Now we will compare this approximation to the STO. Below is the STO
equation. The "squiggle" is a Greek character standing for an orbital
coefficient, similar to alpha in gaussian curves. For hydrogen, the "squiggle"
(better known as "xi", is equal to 1.24.

Under alpha, type "xi." In the next cell, type the value.

Now, label the cell next to your GTO heading "Slater" or "STO." In the
cell beneath the heading, enter the Slater equation. Fill down your column.

With the "Radius," "GTO," and "STO" columns highlighted, create a new graph
to compare the GTO and STO calculations.

Notice the major discrepency for small radii values. Also, notice how the
approximations become similar for larger radii.
Part III:
Now we would like to use several gaussians to model the hydrogen orbital.
To get all of the numerical constants required, we will have to use
a basis set as a reference. These are stored in a large database, from
which you can "order" a specific basis set.

Go to the Gaussian Basis Set Order Format http://www.emsl.pnl.gov:2080/forms/basisform.html.
Your browser should look like this:

Order your constants:

Select the STO3G basis set. By selecting the STO3G we are requesting
to receive constants for 3 gaussians.

Next type an h for hydrogen in the "Elements" blank.

Don't put in an email!!

Click Submit.
The result should be two columns of three numbers each. This is a basis set
for hydrogen. The left list is alpha values, and the right list is
coefficients which will be used to sum up the three gaussians.

Copy these constants so you can enter them into your spreadsheet.

At the top of your spreadsheet, next to the other constants, label each
constant and define the value in the next cell.

Next to the STO column, put three new columns, labelled "g1," "g2," and
"g3".

Under each new label, enter the gaussian equation using the values of alpha
from the basis set and the radius as needed.

Fill down the column, and repeat the procedure for the other two columns.

Change the first number under Radius to be 0. When you do this, the whole list
should change to go from 0 to 6 instead of 3 to 3. Also, the values in your
gaussian lists should change, since we are computing them for different radii
now. Check your values against those of your classmates to make sure you are
all getting correct answers.
***Since we looked at the equations on both sides of the nucleus before,
you should have observed that it is exactly the same on both sides. This is
important! We are only interested in the positive side of the graph now, since
we know that both sides will be exactly the same (and since negative and
positive do not have any true physical meaning here!).***

Compare the answers from each new gaussian equation to those of the slater
equation.
Notice that no single one of these three gaussians is an approximation on its
own; they must be combined. The way of combining them is called "Linear
Combination of Atomic Orbitals" (LCAO). The other list of numbers in the basis
set are coefficients for adding up the three gaussians. They are used in the
following equation.

Create another column next to the three gaussians, labelled "LCAO".

Put the equation, referring to the appropriate constants, in the first
cell.

Fill the column.

Create a graph of all three new gaussians and the LCAO.

To make this contrast clearer, we are going to examine several gaussian
approximations and compare them to the STO.

Label three new columns for the two equations needed for a STO2G model,
and the STO2G model itself.

Do the same for a STO6G

Now go back to the order form for basis sets and find the basis sets, both
alpha and d values, for STO2G and STO6G.

Put these values into your spreadsheet as before.

Set up each column with the neccessary equations.

In two new columns, use the LCAO equation for both of the 2G and 6G
approximations.

Graph the 1 gaussian (the first one created), 2 G, 3 G , 6G, and STO all on a
single graph. **DO NOT put all of the gaussians for each approximation on
there, ONLY THE SUMMATIONS.**
