MOLECULAR MODELING
Chemists frequently use models to better understand chemical phenomena. In this laboratory models are used to analyze chemical bonding and molecular structure. We will start with the Lewis electron dot (Localized Electron) model1 and use it to predict the electronic structure of some simple molecules. We will then use ball-and-stick models to construct three-dimensional models of those molecules. The three-dimensional models will be refined using Gillespie’s Valence State Electron Pair Repulsion (VSEPR) model2. Finally, we will use Spartan3, a computational chemistry program, to calculate the molecular geometry. The results from Spartan will then be compared with the VSEPR results.
Computational chemistry – In theory, all of chemistry can be understood as solutions of Schrödinger’s wave equation. In practice, solution of the Schrödinger equation is so complicated that chemically accurate calculations for molecules involving more than three electrons and three nuclei, e.g. H3, are beyond the reach of current computer hardware/algorithms. To model the properties of larger (and more interesting) molecules, chemists must use some form of approximate computational method.
The simplest, and thus fastest, computational method is called molecular mechanics. In molecular mechanics the bonds of a molecule are treated as springs and the force constants of the springs are determined from experimental measurements. The overall molecular energy is the sum of the energies of stretching, bending and torsion (the molecular force field) for all of the bonds in the molecule. Although molecular mechanics ignores electronic properties, it gives quite good structural results for large molecules such as proteins. Drug companies performed much of the pioneering work on molecular mechanics. Spartan uses a molecular mechanics program called MMFF (for Merck Molecular Force Field) developed by Merck Pharmaceuticals and one called SYBYL developed by Tripos, Inc., a St. Louis drug discovery company.
Semi-empirical molecular orbital methods are the next level of computational accuracy. Semi-empirical methods take advantage of the fact that some of the quantities that are hardest to calculate contribute very little to molecular energy. They can thus be ignored without a significant loss of accuracy. Other hard-to-calculate quantities can be replaced by experimental values. There is a veritable alphabet soup of semi-empirical methods. Spartan uses MNDO (Modified Neglect of Diatomic Overlap), AM1 (Austin Model 1) and PM3 (Parametric Model 3). All three were developed by M. J. S. Dewar’s research group at the University of Texas, Austin between 1975 and 1990. Each method has its own strengths and weaknesses. A lot of semi-empirical computational chemistry is deciding which model works best for the problem you’re working on. In this laboratory, we will use AM1.
At a higher level of computational accuracy is the Hartree-Fock molecular orbital method. This was developed in the 1930s by Douglas Hartree, a British mathematician and Vladimir Fock, a Russian physicist, but was not widely used until digital computers became available in the 1960s. We will not be using the Hartree-Fock method, but when you get to the Calculations… screen in Spartan’s Setup menu, you will see Hartree-Fock listed as one of the choices.
To put this in perspective, for a molecule like SF6, a molecular mechanics calculation takes about one second, a semi-empirical calculation takes about five seconds and a Hartree-Fock calculation takes several minutes. Molecular mechanics does not give molecular energies but can do a reasonable job of predicting which molecular arrangement is most stable. Semi-empirical methods can calculate molecular energies to an accuracy of a few percent, while Hartree-Fock calculations give molecular energies accurate to a few tenths of a percent. The difficulty is that most chemical problems, such as predicting the path of a reaction, require accuracy of ~10-3 percent. Computational procedures that can work at this accuracy are available, but they are very complex and computationally intensive. When Cray introduced their 52.4 teraflop X1 system last week, http://www.cray.com/news/0211/x1announce.html, they stated that the system was aimed at the computing needs of the chemical and pharmaceutical markets.
Procedure - The attached worksheet lists ten molecules that you will study. For each molecule:
1. Write the Lewis electron dot structure,
2. Using the “ball-and-stick” model kit, construct a 3-dimensional model,
3. Use the VSEPR model to predict the bond angles,
4. Log onto your server space and create a folder in which to save your molecular models,
5. Using Spartan, calculate the theoretical bond angles and compare them to the experimental bond angles listed on the work sheet.
When you are ready to calculate the structure, Select Calculations… from the Setup menu. In the Calculate: block, select Eqilibrium Geometry with Semi-Empirical and AM1.
In the Start from: block, select Initial geometry.
In the Subject to: block, check Symmetry.
For a neutral molecule, the Total Charge should be Neutral and the Multiplicity should be Singlet.
In the Compute block, check Elect Charges.
Check Global Calculations.
Click on Submit. Spartan will prompt you to save your model. After saving the model, a window will open notifying you that Spartan has started. Click on OK to close that window. After 5 – 30 seconds, a window will open notifying you that Spartan has completed. Click on OK to close that window and you will see the calculated equilibrium geometry of the molecule.
6. From your knowledge of electron structure, explain any differences between the bond angles predicted by the VSEPR model and those predicted by Spartan,
7. For molecules 9 and 10 (H2O and H2S), comment on the effect of changing the electronegativity of the central atom. Can this be explained using the idea of electron-pair repulsion?
Your lab report (to be turned in at the end of lab) will consist of the completed worksheet annotated with the answers to 6 and 7 above.
References and Notes
Instructions for
Spartan
Spartan runs on the campus network. Start Spartan by double clicking on the icon on your PC’s desktop. When Spartan opens, maximize the screen so that you can read information at the bottom of the active window. Choose New from the File heading on the button bar or click on the blank page icon. A partial periodic chart will appear to the right of the screen. Actually four tabs will appear there: Ent.; Exp.; Pep.; and Nuc. Choose the Exp. tab, since, after all, we do want to be experts in molecular modeling. A more complete periodic chart with several different choices for bonding arrangements for each atom will comes up. These bonding arrangements correspond to the VSEPR theory geometries. You can now build any molecule atom-by-atom by choosing an atom and a configuration and clicking in the open space to the left. These instructions walk you through building two of the molecules on your list. Additionally, they tell you how to measure bond angles.
H2O
Choose O from the periodic table and choose the tetrahedral bonding arrangement, since oxygen has four electron groups around it in the water molecule. Left-click on the big screen, the active window, to the left, and an O atom will appear with four empty bonding sites. You can rotate the atom by left-clicking anywhere on the blank part of the screen and moving the mouse. Choose a H atom and the configuration with one bond from the periodic chart. Left-click on one of the open bond of the O atom in the active window, and a H atom will appear bonded to O atom. Left-click again on one of the other open bonds, and a second O-H bond will appear. You now need to get rid of the two remaining bonding sites on the O atom, or Spartan will automatically fill them with H atoms. To remove these bonding sites, choose the Delete tool, the red asterisk next to the Add Fragment icon, on the button bar. Left-click on the empty bonding sites on at a time and they will disappear. The complete H2O molecule should now be displayed in the active window. You might try choosing the Space Filling model from the Model option on the button bar. This shows the actual size of the electron clouds associated with each atom. You will probably find it easier to build molecules using the Ball and Spoke model, but the Space Filling model will be instructive for viewing.
This structure serves as a starting point for the calculation to determine the lowest-energy geometry for the molecule. Measure the bond angle by choosing the ? angle icon on the button bar. Left-click on a H atom atom, the O atom, and the other H atom, in that order. You should see “Angle (H3,O1,H1) = 109.47º” in the bottom right-hand corner of the screen. You may instead see “Angle (H1,O1,H3) = 109.47º”, depending on which H atom you chose first. This display tells you that your initial guess for the bond angle is 109.47º, that predicted by VSEPR theory.
To have Spartan determine the lowest-energy geometry, choose Setup from the button bar and choose Calculations. In the dialogue box, choose the options outlined earlier in the lab handout.
SF4
Clear your screen by choosing Edit from the button bar and selecting Clear. Make the SF4 molecule by selecting the Add Fragment icon, a bold plus sign, and choosing S from the periodic table to the right. Choose the trigonal bipyramidal configuration for S and click in the active window. (The Lewis dot structure shows 5 groups of electrons around the S atom, four bonding pairs and one lone pair.) Rotate the atom until you can see all the available bonding sites. Add F atoms (choose the single bond configuration) to 4 of the available bonding sites of the S atom; leave one of the equatorial sites vacant. Select the Delete tool and remove the remaining bonding site, which is actually occupied by a lone pair of electrons around the S atom. Determine the lowest-energy geometry as you did for H2O above. Measure the bond angles. Are they different from what VSEPR theory predicts? Why?