CHM111L - Computational Chemistry Exercise
Using SPARTAN to Create Potential Curves for Li2, Li22-, O22+, and O22-.
Today we will carry out Molecular Orbital calculations in order to understand the bonding (or lack thereof) in several homonuclear diatomic molecules and ions. We will construct a potential curve showing how the energy changes with internuclear distance and thereby be able to estimate the equilibrium internuclear distance as well as the bond strength. We will also determine the energy-ordering of the molecular orbitals, and consequently will be able to predict which MO’s are occupied (the electron configuration). Finally, a comparison of orbital symmetries will enable us to deduce how the occupation of certain orbitals impacts the length and strength of the bond.
______
Introduction:
What goes on in the chemical world is dictated by the behavior of electrons and nuclei; and what dictates the behavior of electrons and nuclei is their physics, i.e. their potential and kinetic energies and their distribution in space. One of the main goals of Computational Chemistry is to calculate how the electrons and nuclei might be arranged in atoms and molecules and to estimate the energies associated with those arrangements. Specifically with molecules, one can carry out calculations looking for the lowest energy arrangements of the nuclei and electron waves, then predict that those are the ones likely to be seen in the real world. The good news is that we can use sophisticated computer programs to do the energy calculations for us; all we have to do is tell the program what atoms are involved and what sort of calculation we want do.
Review:
You have already studied the behavior of electrons in atoms and atomic ions. You learned four key aspects of electrons in atoms:
1) they are described by waves (orbitals) that have distinct shapes and sizes,
2) they have specific energies,
3) the ground state involves filling the lowest energy orbitals, and
4) no more than 2 electrons can be “in” an orbital at one time (Pauli Principle).
These same four principles are pertinent for electrons in molecules. The only real difference is that atomic orbitals (AOs) are localized on one atom while molecular orbitals (MOs) can be spread over the entire molecule. Consequently, MOs aren’t characterized by the same 3 quantum numbers used for AOs, and they rarely have any meaningful quantum numbers. The good news is that there are some simple ways to describe the MOs. MOs in homonuclear diatomic molecules are among the simplest MOs, and you will gain some experience with them today.
1) With homonuclear diatomic molecules the MOs are spread evenly between the two atoms, and they resemble combinations of the same AO on each atom. These combinations of AOs can be characterized by their symmetry along the internuclear axis and their symmetry around the internuclear axis. MOs that lie on the internuclear axis are referred to as sigma orbitals (σ), and those that lie on either side of the axis are referred to as pi orbitals (π). MOs that are symmetric along the internuclear axis are referred to as bonding and those that are asymmetric (i.e., change sign around the center) are referred to as antibonding and are identified with an aterix *.
These various distinctions lead to the following 6 homonuclear diatomic MOs:
[red lobes represent one sign (say, +) and blue the other (say, -).]
σ2s σ2s*
σ2p σ2p*
π2p π2p*
[The π MOs actually represent 2 MOs, one in the plane of the paper and one coming out at you.]
2) The relative energies of these MOs follows the general notion of “more nodes means higher energy”: σ2s < σ2s* < σ2p < σ2p* and π2p < π2p*. Just where the π MO energies overlap with the σs differes for different molecules – they are usually near the σ2p but less than the σ2p*.
3) & 4) Once the energy ordering is known, one can write a valence electron configuration by listing the occupied MOs and their occupations.
For O2, in which there are 12 valence electrons, the valence electron configuration would be …
σ2s2 σ2s*2 σ2p2 π2p4 π2p*2
Calculations:
One of a molecule’s most important properties is its geometry or shape, and it can be determined and understood through the use of these energy calculations. The program …
…starts with some initial configuration of the nuclei (a guess made by you or the program),
…then it evaluates the energy the molecule would have in that configuration,
…then examines the forces on the nuclei in that configuration, and allows them to move slightly in response to those forces to make a new configuration,
…then repeats the energy calculation to see if the overall energy is any lower.
This process is repeated until a lowest energy geometry is identified, and that minimum-energy configuration is usually a pretty good approximation to the true configuration found in the real world.
One could say that the energies of all the possible arrangements of the nuclei describe a potential energy function, and the actual arrangement of the nuclei in real life tends to find its way to the lowest energy part of that function.
In today’s exercise we will “brute-force” our way through the energy optimization process outlined above.. Specifically, we will attempt to predict the equilibrium internuclear distance and bond strengths for O22+, O22-, Li2, and Li22 . To do this, we will calculate the energies of the nuclei and electrons if the nuclei are held at a wide range of internuclear distances, from very far apart to very close together. The internuclear distance corresponding to the minimum energy will be our prediction of the equilibrium internuclear distance, and the energy difference between when the nuclei are very far apart and when they are at the equilibrium internuclear distance will be our approximation of the bond strength.
An example of what we will be doing is shown in the figure below for HCl. The points on the curve represent the calculated energies of an arbitrarily chosen set of 10 internuclear distances. The line is just a smooth line drawn through the points to show our expectation for the energies of other internuclear distances if we bothered to calculate them. This is called a potential energy diagram for HCl . The lowest energy separation is expected to be the equilibrium internuclear distance, i.e. the true bond distance you’d see for HCl molecules flying around in the gas phase.
Furthermore, the difference in energy between the energy at the equilibrium internuclear distance and the energy of the fully separated nuclei represents the energy needed to pull the two atoms totally apart, i.e. the bond strength or dissociation energy of the H-Cl bond.
The calculations used in computational chemistry cover a large number of different approaches, some faster than others, and some more reliable than others. The complexity of the molecules involved prohibits any calculation involving more than two particles (one nuclei and one electron) from being calculated exactly, so approximations are used.
- MOLECULAR MECHANICS calculations use empirical observation as well as high quality calculations on small molecules to approximate “force fields” around atoms. The program is told how the atoms are bound to each other, and then it uses the force fields to find how they will arrange themselves to reach the most stable configuration. These are among the fastest but least precise calculations. The VSEPR model discussed in lecture is a qualitative approach to a molecular mechanics technique; it uses a force field which simply indicates that groups of electrons (bonding and non-bonding) repel one another and get as far apart as possible without breaking bonds.
- MOLECULAR ORBITAL calculations are probably the most commonly used calculation in Computational Chemistry. In these calculations one assumes that the electrons are paired in molecular orbitals (MO’s) spread over a molecule, in much the same way we think about them being paired in atomic orbitals (AO’s) in atoms. The energy of the system can be calculated using this assumption of the arrangement of the electrons. In addition to the total energy, another interesting feature of a MO calculation is seeing what the MO’s look like. Most MO calculations attempt to approximate the MO’s as combinations of AO’s. The set of AO’s used to describe the MO’s is called the “basis set”. The more AO’s in the basis set, the more accurate the calculation is likely to be. In addition to the system energy, an MO calculation determines (1) which AO’s are used to describe the various MO’s, and (2) the energy of an electron in each MO; both of which are important in developing a qualitative feel for the chemistry surrounding the molecule. In our exercise we will want to see how some orbital energies go up in energy as the atoms come together, and how some go down.
- SEMI-EMPIRCAL MO Calculations are a subset of MO calculations in
which some of the more rigorous applications of quantum physics are
only approximated – usually using some sort of empirical observations.
The program calculates MO shapes and energies, just like other MO
calculations, but these calculations tend to be quite fast and not always
reliable.
- Ab Initio MO calculations only assume that doubly occupied MO’s
exist; then it uses rigorous quantum physics to calculate MO shapes
and energies. These are the sorts of MO calculations we will use
today, and you will see them referred to as “Hartree Fock calculations”.
- ELECTRON CORRELATION calculations recognize the fact that electrons in molecules don’t actually occupy MO’s, but that their motions are correlated with one another (i.e., when one electron moves, others will adjust their distributions accordingly). These calculations describe the motion of the electrons as a combination of occupations of various MO’s while other electrons are in other combinations of other MO’s. These calculations are quite involved computationally, and have the potential to be exactly correct, but they, consequently, take quite long to complete. We won’t do any of these in this experiment.
Today you will be working with ab initio MO calculations on homonuclear diatomic molecules and molecular ions.
Before you begin, you should be aware of the existence of the four most common types of MO’s: σ bonding MO’s, σ antibonding MO’s, π bonding MO’s, and π antibonding MO’s, and how their occupation affects bond strength. You should also be familiar with the concepts of total energy as opposed to orbital energies. They can be reviewed in Chapter 9 (section 9.7) of Gilbert, etal, Chemistry, 3e.
Also, before you come to lab, you need to install the program on your computer. The program name is SPARTAN, and it can be installed by CAREFULLY following the instructions found at the Chemistry department website.
http://college.wfu.edu/chemistry/course-materials/course-software/wavefunction-spartan-student
You will need to get a code from your instructor before you can download and install the program.
General Procedure:
Your goal is to produce potential energy diagrams for one of the pairs of molecules and/or ions, O22+ and O22-, or Li2 and Li22-. You will also examine the differences in the electron configurations of these species, specifically showing how the orbital energies change with internuclear distance for the underlined molecules. To do this, you will fill-in the following table and plot the data accordingly.
Starting with your underlined molecule, you will set an internuclear distance, do the calculation, then record the orbital symmetries (names) and their energies and the total energy. You’ll set a new internuclear distance and repeat the calculation… and repeat… and repeat…
Next, you’ll set up the calculation for your other molecule or ion and repeat the calculations, though you don’t have to pay attention to the orbital names or energies any more, just the total energies.
You will begin by doing calculations at a set of prescribed internuclear distances, but you may supplement the table with additional internuclear distances if you would like to clarify the MO or total energy behavior in certain regions.
Use Excel to plot the energy of each orbital symmetry for the underlined molecule as a function of internuclear distance. Since two pairs of the eight orbitals are degenerate, the π2px or 2py, and π*2px or 2py, that plot will have six lines on it. You will also make a plot of the total energies of the two species as a function of the internuclear distance. From these plots you should be able to make general conclusions about how bonding and antibonding MO energies change as the internuclear distance decreases. You should also be able to estimate the equilibrium internuclear distance and the bond strength. Furthermore, by relating the differences in your two potential curves to the differences in electron configuration of your two species, you should be able to conclude how occupation of different types of orbitals affects bond strength.
Be sure your lab report includes…
… an abstract containing an objective, the results, and a brief mention of the technique used.
… one plot showing MO energies for your underlined molecule as a function of
internuclear separation. [You do not need to include the completed table.]
… a general statement of the way the bonding and antibonding MO energies change with internuclear distance.
… one plot showing the total energy vs internuclear separation for both of your molecules.
… the electron configuration for each molecule at the equilibrium internuclear distance (or some other appropriate distance if there is no equilibrium).
… identification of the bond lengths and strengths, and an explanation of how the electron configurations affect them.
17
Results for _Li2 / Li22-___:
MO # / 1.5ǺOrbE Sym / 2.0Ǻ
OrbE Sym / 2.3Ǻ
OrbE Sym / 2.6Ǻ
OrbE Sym / 3.0 Ǻ
OrbE Sym / 4.0Ǻ
OrbE Sym / 6.00Ǻ
OrbE Sym / 8.00 Ǻ
OrbE Sym
3
4
5
6
7
8
9
10
ETot Li2
ETot Li22-
17
Results for _____O22+_/ O22-____:
MO # / 1.00ǺOrbE Sym / 1.10Ǻ
OrbE Sym / 1.25Ǻ
OrbE Sym / 1.5Ǻ
OrbE Sym / 1.75Ǻ
OrbE Sym / 2.00Ǻ
OrbE Sym / 3.00Ǻ
OrbE Sym / 6.00 Ǻ
OrbE Sym
3
4
5
6
7
8
9
10
ETot O22-
ETot O22+
17