A Mechanistic Study of the Migratory Insertion Reaction: A

Aug 19, 2013 - Introduction to Density Functional Theory: Calculations by Hand on the Helium Atom. Kyle A. Baseden and Jesse W. Tye. Journal of Chemic...
3 downloads 0 Views 497KB Size
Laboratory Experiment pubs.acs.org/jchemeduc

A Mechanistic Study of the Migratory Insertion Reaction: A Computational Chemistry Exercise Craig D. Montgomery* Department of Chemistry, Trinity Western University, 7600 Glover Road, Langley, British Columbia, Canada, V2Y 1Y1 S Supporting Information *

ABSTRACT: An exercise in computational chemistry is presented whereby the student models the mechanism of migratory insertion. In the case of migratory insertion of CO into the Mn−CH3 bond of [Mn(CH3)(CO)5], semiempirical calculations (PM3) are used to evaluate the concerted mechanism of methyl migration, followed by coordination of CO. The three-centered transition state is determined along with the activation energy, calculated to be 131 kJ/mol. The 16-electron intermediate [Mn(C(O)CH3)(CO)4] is suggested to display an agostic Mn−H interaction as has been postulated previously. The exercise is appropriate as a take-home assignment in a course on inorganic or organometallic chemistry.

KEYWORDS: Upper-Division Undergraduate, Inorganic Chemistry, Laboratory Instruction, Computer-Based Learning, Hands-On Learning/Manipulatives, Coordination Compounds, IR Spectroscopy, Mechanisms of Reactions, Molecular Modeling, Organometallics

M

igratory insertion reactions are of immense importance in organometallic chemistry1 and are involved in a variety of catalytic processes, both homogeneous and heterogeneous. For example, hydrogenation by Wilkinson’s catalyst requires the insertion of the olefin into the Rh−H bond, as does hydroformylation. The Monsanto process involves the insertion of CO into a Rh−CH3 bond, and the heterogeneous process of Ziegler−Natta polymerization involves the insertion of the olefin into a Ti−alkyl bond.1 Insight into the mechanism is gained by considering the thermal reaction of [Mn(CH3)(CO)5] and 13CO as in Scheme 1, a case of 1,1-migratory insertion (where, in the product, both

An alternative mechanism would involve migration of the CO (rather than the methyl group) and its insertion into the Mn−CH3 bond in step 1. However, it has been found that the insertion of CO into cis-[Mn(CH3)(13CO)(CO)4] results in a product mixture containing cis and trans products in a 2:1 ratio (as well as [Mn(13C(O)CH3)(CO)5]), as shown in Scheme 2. The presence of trans product eliminates the possibility of intramolecular CO insertion as the operative mechanism, but is consistent with the mechanism of intramolecular methyl migration.2,4 A final experiment supportive of the methyl migration mechanism (versus the intramolecular CO insertion mecha-

Scheme 1. Migratory Insertion of CO in [Mn(CH3)(CO)5]

Scheme 2. Migratory Insertion of CO in cis[Mn(CH3)(13CO)(CO)4]

the Mn and CH3 groups are bound to the same atom of the inserting CO ligand). The observation that the incoming 13CO ligand is not directly inserted into the Mn−CH3 bond but rather occupies a position cis to the newly formed acyl moiety suggests that the mechanism proceeds by intramolecular migration of the CH3 ligand onto a CO ligand in a cis position (step 1), followed by coordination of the labeled 13CO in the position vacated by the methyl group (step 2).2 This methyl migration likely proceeds in a concerted manner and results in an intermediate that features an agostic Mn−H interaction.3 © XXXX American Chemical Society and Division of Chemical Education, Inc.

A

dx.doi.org/10.1021/ed300526k | J. Chem. Educ. XXXX, XXX, XXX−XXX

Journal of Chemical Education

Laboratory Experiment

nism) involves the use of the chiral Fe(II) complex, CpFe(Et)(PPh3)(CO), as the starting material as illustrated in Scheme 3. In this case insertion of CO proceeds with inversion at the Fe(II) center.5 This is again consistent with the mechanism of alkyl migration. Scheme 3. Migratory Insertion of CO in Chiral Metal Complexes

Computational chemistry is increasingly finding application in chemical education as a means of wedding the theory of lecture material with the practical nature of laboratory experimentation.6−31 The purpose of this exercise is to model the mechanism of migratory insertion involving intramolecular migration of the methyl group, followed by coordination of CO. A number of such studies3,32 have been undertaken using DFT calculations; however, such calculations are not practical for an undergraduate exercise, due to the time required (see discussion accompanying Table 1 below). Herein the mechanism is modeled using semiempirical (PM3) calculations. The exercise is appropriate as a take-home assignment in a course on inorganic or organometallic chemistry. The use of computational chemistry can allow students to better visualize the mechanism and deepen their understanding of the process of migratory insertion. Indeed that has been the experience of students who have done this assignment; specifically, using the calculations to construct an energy profile diagram was cited as a significant aid to understanding the mechanism.

■ ■

2. Students open the [Mn(CH3)(CO)5] file and rename it as the “Transition State”. The “Guess Transition State” command in Spartan Student (or similar commands in other software) is used to draw a reaction arrow moving the electron pair in the Mn−CH3 bond to the carbon atom of a CO ligand in a cis position. The transition state is calculated using a semiempirical calculation (PM3), including a calculation of the IR spectrum. 3. From the IR spectrum, students animate the atomic motion that would bring about the transition state from the reactant molecule [Mn(CH3)(CO)5]. 4. Students note the energy of each molecule (in kJ/mol) and record the data in a table. 5. From the results, students calculate the activation energy and draw a reaction profile. A figure of the transition state should be included and discussed, and the atomic motions necessary to move from the reactant molecule to the transition state should be identified. Likewise a figure of the intermediate [Mn(C(O)Me)(CO)4] should be included and discussed.

HAZARDS There are no hazards associated with this lab. RESULTS AND DISCUSSION It is always important when employing computational chemistry as a means of illustrating chemical principles that one first determine the accuracy of the models. Table 1 displays Table 1. Comparison of Calculated and Experimental Bond Lengths for [Mn(CH3)(CO)5] Bond Distance in [MnCH3(CO)5]/Å



PROCEDURE The software employed here is Spartan Student 5.0.0.33 Although it is assumed that the students have a basic understanding of the software and of computational chemistry, more detailed step-by-step instructions are provided in the Supporting Information. It may be possible to employ other computational software; however, the method and step-by-step instructions may have to be adapted, particularly with regard to the determination of the transition state. The exercise can be utilized as a take-home assignment, where the students are asked to calculate the various models along with their energies, construct an energy profile diagram for the reaction, and determine the activation energy as well as the likely transition state and intermediate in the mechanism. The time taken by students (who are somewhat familiar with Spartan) to build the models and perform the calculations is approximately one hour. The steps are as follows: 1. Students build the following molecules and then optimized, first with molecular mechanics, then by semiempirical methods (PM3): [Mn(CH3)(CO)5], [Mn(C(O)CH3)(CO)4], [Mn(C(O)CH3)(CO)5], [Mn(CO)5]+, CO, and CH3−. The final PM3 optimization should include the calculation of molecular orbitals and atomic charges. For purposes of comparison, the calculations are repeated for [Mn(CH3)(CO)5] using DFT (B3LYP-6-31G*). Students save the models and output of the calculations in each case.

Bond Mn− CH3 Mn− COcis Mn− COtrans C−O

Model Calculation (PM3)

Model Calculation DFT (B3LYP/ 6-31G*)

Lit. Model Calculation DFT (B3B86/split valence)3

Experimental34

2.129

2.191

2.17

2.20

1.860−1.863

1.852−1.855

1.84

1.85

1.826

1.829

1.83

1.83

1.157 (cis) 1.166 (trans)

1.148 (cis) 1.152 (trans)

1.15

1.16 (cis) 1.17 (trans)

a comparison of some selected bond lengths for the reactant complex [Mn(CH3)(CO)5] obtained in this study by both semiempirical (PM3) and DFT (B3LYP with 6-31G* basis set) with previous DFT theoretical studies 3,32 as well as experimental values.34 The differences between the semiempirical calculated values and the experimental values range from 0.004 Å in the case of the Mn−COtrans bond length to 0.071 Å in the case of the Mn−CH3 bond. When DFT (B3LYP with 6-31G*) calculations were employed in this study, the error in these calculated values decreased to 0.001 Å and 0.009 Å for the Mn−COtrans and Mn−CH3 bonds, respectively. Although it is clear that more accurate values are obtained by utilizing DFT calculations, nevertheless, this also results in an unnecessary increase in the time required for the exercise. For example, starting from a model previously optimized by molecular mechanics, the CPU time required to optimize the final product [Mn(C(O)CH 3)(CO)5] by semiempirical B

dx.doi.org/10.1021/ed300526k | J. Chem. Educ. XXXX, XXX, XXX−XXX

Journal of Chemical Education

Laboratory Experiment

Figure 2. These results suggest the activation energy for the three-centered concerted mechanism to be 131 kJ/mol. This compares to an experimental value of 119.3 kJ/mol for ΔG⧧.32,35 Previous DFT calculations resulted in a calculated value of 106 kJ/mol for ΔG⧧.32 Alternately one could propose a dissociative mechanism whereby the Mn−CH3 bond is broken and the methyl reattaches on the cis carbonyl ligand. The activation energy for such a mechanism can be estimated by comparing the energies of the [Mn(CH3)(CO)5] model with that of the [Mn(CO)5]+ and CH3− models; the activation energy in that case would be 1033 kJ/mol. Comparison of these two values calculated here for the activation energy also supports the three-centered concerted mechanism. Of course, the stereochemistry observed in Scheme 1 also eliminates the possibility of the dissociative mechanism, in that such a mechanism would be expected to result in some trans product. It is also of interest to consider the intermediate [Mn(C(O)CH3)(CO)4] as it is a 16-electron species. A number of structures have been postulated previously for this proposed intermediate, including structures displaying an η1-acyl ligand, an η2-acyl ligand, and an η1-acyl ligand featuring an agostic Mn−H interaction, as in Figure 3.

methods was 0.21 s (actual time of approximately 1−2 s) compared with 0.79 s (actual time of approximately 28 min) when DFT methods are employed. Furthermore the activation energy calculated herein as well as the transition state (as reported below) compare very favorably with previous theoretical studies as well as experimental work. Therefore it seems that the semiempirical calculations employed in this study are sufficient for the purposes of an undergraduate exercise. The calculated transition state is shown in Figure 1. From the animations of atomic motion, two vibrational modes are

Figure 1. The transition state for the migratory insertion reaction [Mn(Me)(CO)5] + CO → [Mn(C(O)Me)(CO)5].

identified that are necessary for the molecule to adopt the transition state structure. Obviously the carbon atoms of the methyl ligand and one of the cis carbonyl ligands must move together via a C−Mn−C bending mode, and the methyl group also needs to rotate. The energies of the various models are included in Table 2 and are then used to construct the reaction profile shown in

Figure 3. Proposed structures for the intermediate [Mn(C(O)CH3)(CO)4].

Previous theoretical studies3,32 had proposed the likelihood of the agostic structure for the intermediate, and indeed, when the [Mn(C(O)CH3)(CO)4] model was optimized in this study using PM3 methods, such an agostic interaction was observed (see Figure 4). The Mn−H distance was calculated to be 2.29 Å

Table 2. Calculated (PM3) Energies for the Molecular Models Compound

Energy /(kJ/mol)

[Mn(CH3)(CO)5] CO [Mn(CH3)(CO)5] transition state [Mn(CO)5]+ CH3− [Mn(C(O)CH3)(CO)4] [Mn(C(O)CH3)(CO)5]

−1513.75 −82.61 −1382.68 −695.77 +215.52 −1427.73 −1718.93

Figure 4. Structure of the intermediate [Mn(C(O)CH3)(CO)4] consistent with an agostic Mn−H interaction.

by DFT calculations previously3 and 1.938 Å herein using PM3 calculations. DFT calculations yielded a Cmethyl−Ccarbonyl−Mn bond angle of 83°3 whereas the PM3 calculations in this study yielded a value of 84.80°.



CONCLUSIONS The computational chemistry exercise described herein uses PM3 calculations to consider the mechanism of migratory insertion. Using the example of [Mn(CH3)(CO)5], the exercise demonstrates that the mechanism is likely a concerted one,

Figure 2. Reaction profile for [MnMe(CO)5] + CO → [Mn(C(O)Me)(CO)5]. C

dx.doi.org/10.1021/ed300526k | J. Chem. Educ. XXXX, XXX, XXX−XXX

Journal of Chemical Education

Laboratory Experiment

(12) Cook, A. G.; Kreeger, P. K. Reaction of morpholine with t-butyl acetoacetate: a study in kinetic vs thermodynamic control, product identification, and molecular modeling. J. Chem. Educ. 2000, 77, 90− 92. (13) Feller, S. E.; Dallinger, R. F.; McKinney, P. C. A program of computational chemistry exercises for the first-semester general chemistry course. J. Chem. Educ. 2004, 81, 283−287. (14) Fuhrer, T. J.; Franciose, N. E. Toward a course in molecular modeling for upper division chemistry majors. Chem. Educ. 2010, 15, 39−41. (15) Graham, K. J.; Skoglund, K.; Schaller, C. P.; Muldoon, W. P.; Klassen, J. B. Molecular modeling to predict regioselectivity of hydration reactions. J. Chem. Educ. 2000, 77, 396−397. (16) Habata, Y.; Akabori, S. Teaching 1H NMR spectrometry using computer modeling. J. Chem. Educ. 2001, 78, 121−123. (17) Haddy, A. Molecular modeling exercises and experiments: using a molecular modeling program to calculate electron paramagnetic resonance hyperfine couplings in semiquinone anion radicals. J. Chem. Educ. 2001, 78, 1206−1208. (18) Helms, E.; Arpaia, N.; Widener, M. Molecular modeling exercises and experiments assigning the NMR spectrum of glycidol: An advanced organic chemistry exercise. J. Chem. Educ. 2007, 84, 1328− 1330. (19) Keller, J. W. Lewis acid-catalyzed Diels-Alder reaction of carvone with isoprene. Using two-dimensional NMR and molecular modeling to solve a stereo- and regiochemical puzzle. Chem. Educ. 2006, 11, 262−266. (20) Linenberger, K. J.; Cole, R. S.; Sarkar, S. Looking Beyond Lewis Structures: A General Chemistry Molecular Modeling Experiment Focusing on Physical Properties and Geometry. J. Chem. Educ. 2011, 88, 962−965. (21) Montgomery, C. D. π-backbonding in carbonyl complexes and carbon-oxygen stretching frequencies: A molecular modeling exercise. J. Chem. Educ. 2006, 84, 102−105. (22) Montgomery, C. D. Mechanisms of pentacoordinate pseudorotation. A molecular modeling study of PF5. J. Chem. Educ. 2001, 78, 844−846. (23) Montgomery, C. D. Integrating molecular modeling into the inorganic chemistry laboratory. J. Chem. Educ. 2001, 78, 840−844. (24) Wang, L. Using Molecular Modeling in Teaching Group Theory Analysis of the Infrared Spectra of Organometallic Compounds. J. Chem. Educ. 2012, 89, 360−364. (25) Nataro, C.; Ferguson, M. A.; Bocage, K. M.; Hess, B. J.; Ross, V. J.; Swarr, D. T. Lewis acid-base, molecular modeling, and isotopic labeling in a sophomore inorganic chemistry laboratory. J. Chem. Educ. 2004, 81, 722−724. (26) Spessard, G. O. Use of molecular modeling in teaching organic and organometallic chemistry. Abstracts of Papers, 221st ACS National Meeting, San Diego, CA, United States, April 1−5, 2001; American Chemical Society: Washington, DC, 2001; CHED-129. (27) Simpson, J. M.; Rivera, O. Gas chromatography and molecular modeling: a correlation experiment for the undergraduate laboratory. J. Chem. Educ. 2001, 78, 942−943. (28) Voss, E. J.; Parlatore, I. A.; Baker, D. S. Synthesis, IR spectroscopy, and vibrational analysis of benzenetricarbonylchromium(0): An inorganic experiment incorporating molecular modeling. Abstracts of Papers, 221st ACS National Meeting, San Diego, CA, United States, April 1−5, 2001; American Chemical Society: Washington, DC, 2001; COMP-216. (29) Shattuck, T. W.; Dunham, S. U.; King, D. W.; Dunham, S. U.; Thamattoor, D. M.; Millard, J. T.; Mundy, B. P.; Conry, R. R.; Hennessy, M. H. Molecular modeling in the undergraduate curriculum. Abstracts of Papers, 221st ACS National Meeting, San Diego, CA, United States, April 1−5, 2001; American Chemical Society: Washington, DC, 2001; CHED-1050. (30) Palmer, D. R. J. Integration of computational and preparative techniques to demonstrate physical organic concepts in synthetic organic chemistry: An example using Diels-Alder Reactions. J. Chem. Educ. 2004, 81, 1633−1635.

proceeding by migration of the CH3 ligand to a cis-CO via a three-centered transition state, rather than a dissociative one. The activation energy was calculated to be 131 kJ/mol. The structure of the [Mn(C(O)CH3)(CO)4] intermediate features an agostic Mn−H interaction, as has been suggested previously. Students are readily able to perform the calculations and the exercise and find it helpful in visualizing the mechanism of migratory insertion and confirming the viability of such a mechanism.



ASSOCIATED CONTENT

* Supporting Information S

A student handout is provided that includes detailed step-bystep instructions. This material is available via the Internet at http://pubs.acs.org.



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS C.D.M. acknowledges the Institute of Chemistry at the Hebrew University of Jerusalem for a Visiting Professorship and is grateful to Trinity Western University and the Natural Sciences and Engineering Research Council of Canada for financial support.



REFERENCES

(1) Meissler, G. L.; Tarr, D. A. In Inorganic Chemistry; Pearson Prentice Hall: Upper Saddle River, NJ, 2004; pp 528−532. (2) Noack, K.; Calderazzo, F. Carbon monoxide insertion reactions. V. The carbonylation of methylmanganese pentacarbonyl with carbon13 monoxide. J. Organomet. Chem. 1967, 10, 101−104. (3) Derecskei-Kovacs, A.; Marynick, D. S. A New Look at an Old Reaction: The Potential Energy Surface for the Thermal Carbonylation of Mn(CO)5CH3. The Role of Two Energetically Competitive Intermediates on the Reaction Surface, and Comments on the Photodecarbonylation of Mn(CO)5(COCH3). J. Am. Chem. Soc. 2000, 122, 2078−2086. (4) Flood, T. C.; Jensen, J. E.; Statler, J. A. Stereochemistry at manganese of the carbon monoxide insertion in pentacarbonylmethylmanganese(I). The geometry of the intermediate. J. Am. Chem. Soc. 1981, 103, 4410−4414. (5) Brunner, H.; Hammer, B.; Bernal, I.; Draux, M. Optically active transition-metal complexes. 85. The BF3-promoted carbonylation of CpFe(CO)(L)Me compounds - a reaction stereospecific at low temperatures. Organometallics 1983, 2, 1595−1603. (6) Amarne, H. Y.; Bain, A. D.; Neumann, K.; Zelisko, P. M. Extensions of a basic laboratory experiment: [4 + 2] and [2 + 2] cycloadditions. J. Chem. Educ. 2008, 85, 104−106. (7) Barrows, S. E.; Eberlein, T. H. Cis and trans isomerization in cyclic alkenes: A topic for discovery using the results of molecular modeling. J. Chem. Educ. 2004, 81, 1529−1532. (8) Box, V. G. S. Using Molecular Modeling To Understand Some of the More Subtle Aspects of Aromaticity and Antiaromaticity. J. Chem. Educ. 2011, 88, 898−906. (9) Clausen, T. P. Combining a Standard Fischer Esterification Experiment with Stereochemical and Molecular-Modeling Concepts. J. Chem. Educ. 2011, 88, 1007−1009. (10) Coleman, W. F. Molecular Models of Alkyl Carboxylic Acids and Amines. J. Chem. Educ. 2010, 87, 457−458. (11) Coleman, W. F. Molecular Models of Phthalocyanine and Porphyrin Complexes. J. Chem. Educ. 2010, 87, 346. D

dx.doi.org/10.1021/ed300526k | J. Chem. Educ. XXXX, XXX, XXX−XXX

Journal of Chemical Education

Laboratory Experiment

(31) Parnis, J. M.; Thompson, M. G. K. Modeling stretching modes of common organic molecules with the quantum mechanical harmonic oscillator: an undergraduate vibrational spectroscopy laboratory exercise. J. Chem. Educ. 2004, 81, 1196−1198. (32) Wang, X.; Weitz, E. A density functional theory study of alkyl group migration in RMn(CO)5 complexes. J. Organomet. Chem. 2004, 689, 2354−2360. (33) Spartan Student 5.0.0; Wavefunction, Inc.: Irvine, CA, 2011. (34) Seip, H. M.; Seip, R. Structure of methylmanganese pentacarbonyl, CH3Mn(CO)5. Acta Chem. Scand. 1970, 24, 3431− 3433. (35) Andersen, J. M.; Moss, J. R. Synthesis of an Extensive Series of Manganese Pentacarbonyl Alkyl and Acyl Compounds: Carbonylation and Decarbonylation Studies on [Mn(R)(CO)5] and [Mn(COR)(CO)5]. Organometallics 1994, 13, 5013−5020.

E

dx.doi.org/10.1021/ed300526k | J. Chem. Educ. XXXX, XXX, XXX−XXX