Things to remember: To run an AMSOL job you will need to have an input file named myfile.dat (where myfile is whatever you want, but .dat is critical). Remember, to create this you need to use the vi editor. You can either create a file directly while in vi, or go into input mode and paste something from the Mac. Once your input file is complete, type amsoli myfile (don't include the .dat suffix) and it will be submitted to run interactively. When the prompt returns, your run is complete (although that doesn't necessarily mean it was successful!) If a job seems to be running a long time, you probably made some sort of mistake (a very bad geometry, for instance), and you should ^C to interrupt the process and look at whatever output file may have been created to see what you did wrong. If everything worked well, you will have a myfile.out and myfile.arc file with the results of the run found therein (the .arc file is very condensed compared to the .out file).

You can look at the output files by using more myfile.out, in which case it will scroll one page at a time, or by using vi, in which case you can use all of the editor commands to move around (be careful not to change the file, though-if you accidentally delete or change something, you can leave the file unaffected by entering :q! when in command mode).

The Problems (do any two of three):

1. Jorgensen et al. reported in JACS 120 (1998) 7637 the synthesis of a tertiary silicenium ion in solution. This assignment is almost certainly wrong (even JACS referees aren't perfect!)-silicenium ions are terrifically unstable compared to carbonium ions. Indeed, the only unambiguously established silicenium ion, tris-mesitylsilicenium (made by Lambert et al.), owes its stability primarily to steric crowding about silicon, not aryl stabilization; triphenylsilicenium cannot be formed in solution. This is in marked contrast to triphenylcarbenium, which is rock stable, and begs the question of: Why doesn't aryl substitution stabilize a silicenium ion like it does a carbenium ion?

To answer this question more fully, let's look at the MOs for phenyl carbonium (PhCH2+) and phenylsilicenium (PhSiH2+). Optimize the structures of these two cations at the AM1 level. To get the MO's you will need to use the keyword VECTORS in your input decks. This causes the program to print out the coefficients for each MO near the end of the file. In the NDDO approximation, recall that the sum of the squares of all coefficients of a MO is 1, so the square of the coefficient tells you the percentage contribution of each AO basis function to a given MO. Sketch the p MO's (how many are there?) noting the percentage contribution of each AO to the MO. To sketch a molecular orbital, you need to adopt a convention for how you shade the orbitals for positive vs. negative coefficients (consistency is all that matters), and you need to make sure that you have the cartesian coordinates for the molecule so you know how the orbitals relate to the atomic positions, and you need to scale your various atomic orbitals according to the MO coefficients, i.e., a big coefficient means lots of that atomic orbital, a small coefficient means just a little bit of that atomic orbital, and finally, of course, you need to know which atom is which in your input deck! [Modern programs will do all of this for you with glitzy graphics, but doing it by hand once or twice is good for the soul and provides more insight into what's going on.] Label your MO's with the calculated MO energy. If any of these p MOs correspond to the HOMO or the LUMO, so label them. Notice anything you think strange about the LUMO energies?

Also report the partial atomic charges on the carbonium carbon and the silicenium silicon. With all your analysis in hand, why is an aryl silicenium ion so much less stable than its analogous carbonium ion?

[Some of you, exercising the powerful skepticism that characterizes a scientist, may be wondering how the above problem can prove that Jorgensen and co-workers are wrong. Of course it doesn't-it just suggests why it's so hard to get a stable silicenium ion. We'll defer more convincing proof until we can compute NMR chemical shifts ab initio (next problem set).]

2. Let's return to one of the molecules found in Problem Set 1. For the case of 2-methoxy-1,3-tetrahydropyran, MMX gave a rotational potential consistent with a large anomeric effect. The anomeric effect is a stereoelectronic effect, i.e., it derives from orbital-orbital interactions. Repeat your calculations for this torsional potential at the PM3 level.

AMSOL has a way to simplify the calculation of this potential. Instead of flagging the dihedral angle with a "1" (meaning optimize) flag it with a "-1". A -1 flag means to hold that degree of freedom fixed at whatever the value is in the Z-matrix, optimize everything else, and then repeat that process for every value found on the line(s) following the first blank line after the Z-matrix. All results are provided in a single output file For instance:

AM1

phosphinous acid, rotation coordinate for P-O bond (gas phase)

O

P 1.5 1

H 1.3 1 115. 1

H 1.3 1 115. 1 120. 1 2 1 3

H 1. 1 110. 1 0. -1 3 1 2

30. 60. 90. 120. 150. 180.

would calculate the rotational potential for phosphinous acid (H2POH) by 30 degree intervals starting from zero.

Compare the PM3 potential to the MMX potential. What, if any, are the differences? If there are important differences, what could you do to establish which level of theory (if either) is the more correct?

3. Design a problem of your own that uses semiempirical MO theory to illustrate some chemically interesting concept. Write down the problem, and then provide the answer. Note that you do not need to pick something that semiempirical theory performs well for! However, if you create a problem where NDDO theory clearly gives the wrong prediction, provide some discussion in your answer of why the level fails to be accurate. To get more of a feel for typical problems, feel free to drop by the website and look at second problem sets from previous years. Grading will be based on quality of the problem and originality.