1,3-Metal–Carbon Bonding and Alkyne Metathesis: DFT Investigations

Oct 3, 2012 - Cβ, as the Cα−Cβ bonds possessed largely formal Csp2−Csp2 single-bond character. Electron density analysis showed a. “catastrop...
0 downloads 0 Views 4MB Size
Article pubs.acs.org/Organometallics

1,3-Metal−Carbon Bonding and Alkyne Metathesis: DFT Investigations on Model Complexes of Group 4, 5, and 6 Transition Metals Cherumuttathu H. Suresh*,†,‡ and Gernot Frenking† †

Fachbereich Chemie, Philipps-Universität Marburg, Hans-Meerwein-Strasse, D-35032 Marburg, Germany Chemical Sciences and Technology Division, National Institute for Interdisciplinary, Science and Technology, CSIR, Trivandrum 695 019, India



S Supporting Information *

ABSTRACT: The formation of metallacyclobutadienes (MCBs) from chloro-ligated alkylidyne complexes of group 4, 5, and 6 transition metals (MCln(C3H3)) has been studied at the BP86/ def2-TZVPP level. All the MCBs showed M−Cβ distances (∼2.1 Å) very close to M−Cα distances (1.8−2.0 Å), suggesting a bonding interaction between the metal and the β-carbon (1,3-MC bond). Energy decomposition analysis using C2v symmetric structures revealed that a b2 orbital composed of mainly metal dπ and Cβ pπ overlap and an agostic a1 orbital contributed to the orbital interaction of the 1,3-MC bond. The bond order of the 1,3MC bond is a minimum of 0.26 for M = Cr and a maximum of 0.43 for M = Ta. Further, all the MCBs showed a characteristic δ orbital interaction through an a2 orbital, which contributed to the double-bond character of M−Cα bonds (bond order 1.27− 1.44). Although the formation of b2 and a2 orbitals increased the M−C interactions, they significantly reduced the π interactions within the C3H3 fragment (C−C bond order 1.09−1.18). 1,3-MC bonding suggested a planar tetracoordinate configuration for Cβ, as the Cα−Cβ bonds possessed largely formal Csp2−Csp2 single-bond character. Electron density analysis showed a “catastrophic” character of the 1,3-MC bond. In groups 4 and 5, MCBs were more stable than the isomeric η3-structures (metallatetrahedranes). A mechanistic study on the reaction between acetylene and alkylidyne complex MClnCH showed that a nearly barrierless and exothermic pathway exists for MCB formation (exothermic value 75−102 kcal/mol for groups 4 and 5; 6− 27 kcal/mol for group 6). The rich metathesis chemistry associated with Mo and W is attributed mainly to the moderate activation energy required for the alkyne disproportionation step of metathesis. A mechanistic possibility other than Katz's is also proposed for alkyne metathesis that showed that the 1,3-MC bonded MCB complex can act as a metathesis catalyst by reacting with alkyne to form a bicyclic intermediate and subsequently disproportionating to yield the alkyne and the MCB. For this mechanism to be effective, rearrangement of the bicyclic intermediate to a more stable cyclopentadienyl complex has to be prevented.



theory.22,27−33 All the reported MCB intermediates of alkyne metathesis showed a characteristic short M−Cβ distance typically in the W−C single bond range 2.1−2.3 Å, and this was also very close to the value of the M−Cα distance (1.8−2.0 Å). The short M−Cβ distance suggested significant bonding interaction between the metal and the Cβ. Fenske−Hall molecular orbital calculations on [W(C3H3)]3+ reported by Bursten in 1983 confirmed that the short M−Cβ distance is due to a strong orbital interaction between dπ on the W and pπ on the Cβ.34 Later theoretical studies by Ziegler et al. and Lin and Hall have supported the dπ−pπ-type bonding for the short M− Cβ distance in MCB structures.35,36 Ziegler et al. also argued that the double-bond character of the M−Cα bond as well as partial double-bond character of Cα−Cβ bonds imparts

INTRODUCTION Alkyne metathesis provides an important homogeneous catalytic procedure for the creation of CC triple bonds.1−11 Although the first alkyne metathesis catalyst12 was discovered in 1974, only a limited number of well-defined catalysts have been reported since then for the metathesis of CC triple bonds, and several of them were alkylidyne complexes of molybdenum and tungsten.13−22 Alkyne metathesis often passes through stable metallacyclobutadiene (MCB) intermediates, and in the Katz mechanism,23−25 a [2 + 2] cycloaddition step between alkyne and metal alkylidyne is responsible for the formation of the MCB intermediate. In 1982 Schrock et al. experimentally established this mechanism by isolating and characterizing the X-ray structure of the MCB intermediate Cl3W[C-t-BuCMeCMe] in the reaction between 2-butyne and the alkylidyne complex W(C-t-Bu)(dme)Cl3.26,27 Since then, a few more stable MCB structures are reported by both experiment and © 2012 American Chemical Society

Received: July 26, 2012 Published: October 3, 2012 7171

dx.doi.org/10.1021/om3007097 | Organometallics 2012, 31, 7171−7180

Organometallics

Article

Figure 1. Optimized structures of MCln(C3H3) systems at the BP86/def2-TZVPP level. (a) 1,3-MC bonded metallacyclobutadienes of group 4, 5, and 6 metals. (b) η3-Complexes (metallatetrahedranes) of group 4, 5, and 6 metals. (c) Normal metallacyclobutadienes (without 1,3-MC bonding) of group 6 metals. Relative energy in kcal/mol, and bond lengths in Å.

and also suggested that Cβ exists in the planar tetracoordinate state in the MCB complex. The aim of this paper is to show that 1,3-MC bonding not only is unique to group 6 metals but also can be observed in other groups, particularly the earlier ones. To prove this, the metallacyclobutadienes of the group 4 and 5 transition metals are modeled and compared with complexes of group 6 metals. The 1,3-MC bonding has been examined with EDA and electron density analysis. Further, the importance of this bonding has been discussed in the context of alkyne metathesis by studying the potential energy surface for the [2 + 2] type cycloaddition and cycloreversion processes. Moreover, the utilization of this bonding in designing realistic models of MCB complexes and a proposal to use them as potential catalysts for alkyne metathesis are discussed.

structural restrictions within the four-membered ring of MCB to yield a short M−Cβ distance.35 Also a study by Goddard et al. based on generalized valence bond theory did not support the M−Cβ bonding in the molybdenum MCB complex.37 Very recently, we have made a thorough analysis of the nature of the M−C bonding in the MCB complexes (Cl 3 WC 3 H 3 ), (Cl3MoC3H3), and (Cl3CrC3H3) using energy decomposition analysis (EDA) and atoms-in-molecule analysis38 and showed that σ, π, and δ M−C orbital interactions contribute to the overall M−C bonding, and among them the π-type interaction from a b2 orbital was mainly responsible for the M−Cβ bonding. Further, the Mayer bond order values 0.37, 0.31, and 0.26 calculated respectively for (Cl 3 WC 3 H 3 ), (Cl3MoC3H3), and (Cl3CrC3H3) suggested the formation of an M−Cβ bond. We have named this bond as a 1,3-MC bond 7172

dx.doi.org/10.1021/om3007097 | Organometallics 2012, 31, 7171−7180

Organometallics



Article

METHODOLOGY

for this study (ClTi(C3H3), ClZr(C3H3), ClHf(C3H3), Cl2V(C3H3), Cl2Nb(C3H3), Cl2Ta(C3H3), Cl3Cr(C3H3), Cl3Mo(C3H3), and Cl3W(C3H3)) are presented in Figure 1a. Among the group 4 MCB complexes, TiCl(C3H3) is planar and shows C2v symmetry, while ZrCl(C3H3) and HfCl(C3H3) are nonplanar and possess Cs symmetry. All the MCl2(C3H3) systems of group 5 showed C2v symmetry. The structure of CrCl3(C3H3) showed no symmetry, while both MoCl3(C3H3) and WCl3(C3H3) have C2v symmetry (Figure 1a). Optimization of ZrCl(C3H3), HfCl(C3H3), and CrCl3(C3H3) under C2v symmetry restriction gives structures that are only slightly higher in energy (0.5−0.9 kcal/mol) than the respective structures with lower symmetry. This indicates the effect of slight Jahn−Teller distortion in C2v structures. All the MCB complexes depicted in Figure 1a possess short M−Cα bonds (1.80−2.03 Å), indicating strong double-bond character of the bonds. The most remarkable structural feature of these complexes is that they all show short M−Cβ distances (2.01 to 2.25 Å), which is unusual for atoms placed in 1 and 3 positions.33 A bonding interaction between the metal and the β-carbon (1,3-MC bond) can be suspected because the M−Cβ distance is shorter than the sum of the covalent radii of carbon and the metal (2.12−2.48 Å). The Cα−Cβ bond distance is in the range 1.42−1.45 Å, which is close to a formal Csp2−Csp2 single-bond distance. The internal Cα−Cβ−Cα angle is in the range 119−125°, while the external M−Cα−H angles are in the range 154−164°. These structural features are in good agreement with the X-ray structural features reported by Schrock et al. for the MCB complex W(C3Et3)[O-2,6-C6H3(iPr2)]3.33 Since the structural features suggest double-bond character for the M−Cα bond, single-bond character between M and Cβ, and formal Csp2−Csp2 single-bond character for both CC bonds, the drawing given in Scheme 1 is proposed for the

39

Geometry optimizations were carried out with Gaussian03 at the gradient-corrected BP86/def2-TZVPP level of theory40−42 as well as the BP86/TZ2P+ level of theory using the ADF program version 2010.02.43 Essentially, both these levels of theory give almost identical results, and the later theory is chosen for doing the EDA analysis using ADF. The basis sets are triple-ζ quality, which are augmented by polarization functions p and d for hydrogen and d and f for other atoms. Scalar relativistic effects were considered by using the zeroorder regular approximation (ZORA) in ADF calculations.44,45 EDA analysis implemented in ADF employs Morokuma’s energy decomposition scheme46 to Kohn−Sham molecular orbitals as well as uses the extended transition state partitioning scheme of the orbital interactions as developed by Ziegler and Rauk.47−50 Vibrational frequencies were determined by calculating the Hessian matrix analytically. In all cases, the nature of the stationary points and transition states was verified. In the EDA scheme, the complex is divided into two fragments, viz., the inorganic fragment consisting of the metal and the chloro ligands (fragment A) and the organic fragment consisting of the hydrocarbon part (fragment B). The bond dissociation energy (De) required for the fragmentation process can be partitioned into four contributions, viz., electrostatic interaction (ΔEele), Pauli repulsion (ΔEPauli), orbital interaction (ΔEorb), and preparation energy (ΔEprep). ΔEele is calculated with a frozen density distribution of the fragments in the geometry of the complex. ΔEPauli accounts for the energy change associated with the transformation from the superposition of the unperturbed electron densities of fragments ρA + ρB to the wave function Ψ0 = NÂ (ΨA·ΨB), which properly obeys the Pauli principle through explicit antisymmetrization (Â ) and renormalization (N) of the product wave function. ΔEorb gives the stabilizing orbital interaction when the orbitals of the fragments relax to their final form. ΔEprep is the energy that is necessary to promote the fragments A and B from their equilibrium geometry and electronic ground state to the geometry and electronic state that they have in the molecule. The sum of ΔEele, ΔEPauli, and ΔEorb gives the total interaction energy (ΔEint). The sum of ΔEint and ΔEprep gives the bond dissociation energy (De, by definition with opposite sign) given in eq 1.

− De = ΔE int + ΔEprep

Scheme 1. General Structural Drawing for the Metallacycobutadinea Showing 1,3-MC Bonding between the Metal and the β-Carbon

(1)

EDA provides one of the most robust methods for doing quantitative and systematic analysis on the nature of the chemical bonds.51,52 More accounts of the EDA method can be found in the literature.43,51,53 We have also used the ETS-NOCV scheme, a recent extension of the EDA scheme that combines the extended transition state (ETS)47,48,54 approach with the natural orbitals for chemical valence (NOCV)49,50 method. This scheme decomposes the deformation density that is associated with the bond formation, Δρ, into different components of the chemical bond and provides pairwise energy contributions for each pair of interacting orbitals to the total bond energy. The orbital interaction term, ΔEorb, in the ETS-NOCV scheme is given in eq 2. N /2

ΔEorb =

k=1

k=1

M = group 4, 5, and 6 metals.

MCB complexes. This structural drawing also suggests some amount of planar tetracoordinate character for the β-carbon. In our previous paper, the atoms-in-molecule theory56 indicated a catastrophic nature for the 1,3-MC bonding, as a slight contraction of the M−Cβ distance showed the formation of a bond critical point and bond path between M and Cβ as well as two ring critical points, each corresponding to a triangle defined by M, Cα, and Cβ.38 The C3H3 moiety can interact with the metal through η3mode complexation to yield a metallatetrahedrane complex. In the η3-complex (Figure 1b), the M−Cα bond distance is slightly larger than the M−Cα bond distance in the corresponding MCB (Figure 1a). The C−C distance of the η3-complex shows considerable lengthening while going from Ti (1.524 Å) to Zr (1.558 Å) to Hf (1.595 Å). Although, the C−C distance may indicate the formation of a C−C bond between two tetravalent

N /2

∑ ΔEkorb = ∑ υk[−F −TSk ,−k + FkTS,k]

FTS −k, −k

a

(2)

FTS k,k

and are diagonal transition-state Kohn−Sham matrix where elements corresponding to NOCVs with the eigenvalues −υk and υk, respectively. More details of the ETS-NOCV method can be found in the recent literature.49,50,55 Topological properties associated with the bond critical point and Laplacian distributions56 were obtained by using the AIMAll (version 12.06.03) program.57



RESULTS AND DISCUSSION Geometric Features of Metallacyclobutadienes. Optimized geometries of metallacyclobutadiene complexes selected 7173

dx.doi.org/10.1021/om3007097 | Organometallics 2012, 31, 7171−7180

Organometallics

Article

Table 1. EDA Results (BP86/TZ2P+) of MCln(C3H3) Complexes (M = Ti, Zr, Hf, V, Nb, Ta, Cr, Mo, W)a (all energies in kcal/ mol)

a

energy

Ti

Zr

Hf

V

Nb

Ta

Cr

Mo

W

ΔEint ΔEPauli ΔEele ΔEorb ΔEa1 ΔEa2 ΔEb1 ΔEb2 ΔEprep −De

−183.3 325.2 −206.4 −302.0 −101.5 −52.6 −133.2 −14.8 49.2 −134.1

−201.4 353.1 −236.2 −318.3 −100.3 −58.3 −142.5 −17.2 54.3 −147.1

−216.5 401.0 −288.2 −329.3 −98.2 −63.4 −149.9 −17.8 71.8 −144.7

−169.0 428.6 −260.2 −337.5 −139.7 −53.8 −127.6 −16.4 61.6 −107.4

−205.1 460.0 −301.2 −363.8 −138.2 −65.1 −140.9 −19.7 64.0 −141.0

−227.5 514.1 −355.3 −386.3 −146.3 −70.6 −148.1 −21.3 67.8 −159.7

−117.2 358.1 −191.8 −283.6 −98.0 −34.9 −124.5 −26.2 85.2 −32.0

−164.6 404.7 −239.3 −330 −108.4 −50.3 −142.5 −28.8 82.2 −82.4

−193.2 473.4 −305.3 −361.3 −122.7 −57.6 −150.6 −30.3 84.4 −108.8

Interacting fragments MCln and C3H3 are in the quartet ground state.

carbon atoms in the sp3-hybridized state, the angle defined by H, C, and M (163° for Ti, 168° for Zr, and 172° for Hf) supports mostly the sp-hybridized state of the carbon atoms. One of the sp orbitals will be directed toward the bonding direction of the metal and the other toward the bonding direction of the hydrogen, indicating in-plane overlap of p orbitals for C−C σ bonds. Unlike the group 4 η3-complexes, one of the C−C bonds of group 4 is significantly longer than a typical C−C single bond, whereas the other two C−C bonds show equal bond lengths in the range 1.42 to 1.44 Å. An X-ray structure showing very similar structural features was reported for the η3-vanadium complex.58 The C−C bond length of η3complexes of group 6 is in the range 1.43 to 1.47 Å, which is typically observed in reported X-ray structures. Further, the angle defined by H, C, and M is 137° for Cr, 142° for Mo, and 146° for W. The η3-TiCl(C3H3), η3-ZrCl(C3H3), η3-HfCl(C3H3), η3VCl2(C3H3), η3-NbCl2(C3H3), and η3-TaCl2(C3H3) complexes are 24.0, 28.5, 32.6, 28.2, 37.0, and 41.7 kcal/mol less stable than the respective MCB complexes given in Figure 1a, indicating that the conventional η3-coordination of C3H3 with a metal center is much weaker than the unconventional 1,3-MC bonding. However, the η3-CrCl3(C3H3), η3-MoCl3(C3H3), and η3-WCl3(C3H3) complexes are 27.9, 11.1, and 0.8 kcal/mol more stable than the respective MCBs with a 1,3-MC interaction. Our result qualitatively agreed with a previous study by Goddard et al., who used GVB calculations to show that the metallatetrahedrane complex of Mo is 4.8 kcal/mol more stable than the metallacyclobutadiene complex.37 Our results are also consistent with the results of Ziegler et al., who using nonlocal DFT showed that the metallatetrahedrane complex of Mo is 12.4 kcal/mol more stable than the metallacyclobutadiene complex.35 Previously, Lin and Hall reported an MP2 single-point energy evaluation on HF level optimized geometries of metallatetrahedrane and metallacyclobutadiene complexes of WCl3(C3H3) and found that the latter system is more stable than the former by 16.9 kcal/ mol, which is not in good agreement with the present results.36 In a previous study, we have shown that with alkyl substitution on C, the 1,3-MC-bonded complex of tungsten can attain more stability than the η3 complex.38 The X-ray structure of tungstenacyclobutadiene W[C-t-BuCMeCMe]Cl3 is known,27 and calculations at BP86/def2-TZVPP levels showed that this MCB complex is more stable than the corresponding η3-W[C-tBuCMeCMe]Cl3 complex by 2.5 kcal/mol. We expect that an MCB complex has to be more stable than the η3-complex to show good alkyne metathesis catalytic activity.

One of the referees of this paper raised the following important question regarding the structure of the metallacyclobutadiene. How far can a Cβ atom in a four-membered ring get away from the metal? We may speculate that the structural restrictions within the unsaturated four-membered ring will always restrict the position of Cβ within the coordination sphere of the metal, and hence the observed short M−Cβ distance of MCBs is primarily structural and not electronic in nature. However, this argument is not fully correct because metallacyclobutadiene structures with a significant difference in the M−Cα and M−Cβ distances are possible for group 6 MCl3(C3H3) complexes (Figure 1c). In these minimum energy structures, the M−Cβ distance (2.519 to 2.611 Å) is significantly larger (0.549−0.696 Å) than the M−Cα distance. Further, the M−Cβ distance is larger than the sum of the covalent radii of carbon and metal (2.12−2.48 Å). The C− C bond alternation suggests that the shorter C−C bond is double-bond localized and the longer one is a formal Csp2−Csp2 single bond. Hence, the shorter M−Cα bond has to be assigned with more double-bond character than the other. The localized M−Cα double bond, the Cα−Cβ double bond, and the long M−Cβ distance suggest that the p orbital of the Cβ has only a negligible interaction with the metal center for the complexes given in Figure 1c, and therefore a significant 1,3-MC interaction may not exist in such complexes. An intriguing observation is that M−Cα bonds of MCB structures given in Figure 1a are shorter than the M−Cα bonds designated as double bonds for the structures given in Figure 1c. Structures showing bond alternation at the Cα−Cβ bonds or structures showing a larger M−Cβ distance than those reported in Figure 1a are not observed for group 4 and group 5 MCBs. The 1,3-MC-bonded MCB structure of CrCl3(C3H3) is 5.1 kcal/mol less stable than the structure showing C−C bond alternation, whereas 1,3-MC complexes of MoCl3(C3H3) and WCl3(C3H3) are 3.9 and 9.1 kcal/mol more stable than the corresponding structures having C−C bond alternation. Nature of Metal−Carbon Bonding. For all the MCln(C3H3) metallacyclobutadiene complexes (Figure 1a), EDA is carried out with A and B fragments in singlet, doublet, triplet, and quartet electronic configurations. Always C2v symmetry is used for the description of the MCln(C3H3) complexes and the fragments. In the doublet and quartet states A is MCln and B is C3H3, while in the singlet and triplet cases, A and B are respectively MCln+ and C3H3−. The optimized structure and relative energy of C3H3 (doublet and quartet) and C3H3− (singlet and triplet) are presented in Figure S1 of the Supporting Information. The ground state of C3H3 is doublet 7174

dx.doi.org/10.1021/om3007097 | Organometallics 2012, 31, 7171−7180

Organometallics

Article

Figure 2. Plots of the highest occupied molecular orbitals responsible for metal−carbon bonding in VCl2(C3H3). Eigenvalues of the orbitals are given in eV. Isosurface value is 0.08 for MO30 and 0.06 for all other MOs.

Figure 2 for a representative case VCl2(C3H3). The nature of the metal−ligand bonding is very similar for other metals. a1 and a2 are σ and δ orbitals, respectively. b1 orbitals correspond to π orbitals that lie in the plane of the metallacycle, and the b2 orbital corresponds to a π orbital orthogonal to the plane of the metallacycle. The δ orbital a2, the highest lying σ orbital a1, and the inplane π orbital b1 contribute to the strength of M−Cα bonds. The π orbital b2 is mainly responsible for the M−Cβ bond, which is the same as the dπ−pπ orbital reported by Bursten using Fenske−Hall molecular orbital calculations on the [W(C3H3)]3+ fragment.34 Ziegler and co-workers also reported qualitatively similar results on the δ- and π-bonding features of the metal−carbon bonds.35 Although the low-lying b1 and a1 orbitals are largely localized within the C3H3 fragment, they show some minor interactions with the d orbitals on the metal. Since the low-lying b1 and a1 orbitals largely define the C−C σ bonds, their interactions with the metal are agostic in nature.59−61 The EDA data reveal that ΔEorb contributions from a1, a2, b1, and b2 orbitals are respectively in the range 30− 41%, 12−19%, 38−46%, and 5−9%. The largest contributions come from the in-plane a1 and b1 orbitals of the M−Cα bonds, while the contribution of b2 suggests the formation of a weak M−Cβ bond. The stabilizing δ orbital interaction is the highest for WCl 3 (C 3 H 3) (57.6 kcal/mol) and the lowest for CrCl3(C3H3) (34.9 kcal/mol). Group 4 complexes show bond dissociation energies (De) of 134.1, 147.1, and 144.7 kcal/mol for M = Ti, Zr, and Hf, respectively. The vanadium complex in group 5 has a 26.7 kcal/mol lower De value than the titanium complex, whereas in the case of tantalum, De is increased by 15.0 kcal/mol compared to hafnium. Both zirconium and niobium complexes show comparable De values. In group 6, the chromium complex has the lowest De (32.0

and the quartet state is 38.8 kcal/mol higher in energy, whereas the ground state of C3H3− is singlet and its triplet state is 6.1 kcal/mol higher in energy. The structure of the C3H3 fragment in the complex resembles more the optimized quartet structure. The MCln fragment is more stable in the quartet state than the doublet (M = Hf is an exception), whereas the MCln+ fragment is more stable in the triplet state than the singlet ((M = Hf is an exception; Table S1). In the previous study, we have shown that in the case of group 6 complexes ΔEorb was the lowest for the quartet configuration, and hence this configuration was selected for the fragments in the complex. In the case of group 4 and 5 complexes, the neutral fragments in the quartet scheme are more suited for EDA analysis than those in the doublet scheme, whereas with charged fragments, the triplet scheme is slightly more preferred (5−8% on the orbital interaction) than the quartet scheme. In order to give a consistent description of the metal−carbon bonding using neutral fragments C3H3 and MCln, hereafter only the quartet EDA scheme will be discussed unless otherwise stated (EDA data for all other schemes are given in Table S2). Table 1 reports the EDA data for the quartet states for all complexes. The ΔEprep values suggest that the fragments of Ti, Zr, V, and Nb complexes require relatively lower ΔEprep values (49−64 kcal/mol) than the rest. Among the three main components of ΔEint, the destabilizing ΔEPauli is very high for all the cases (325−514 kcal/mol), which can be attributed to the compact ring structure of the complexes. In all the cases, the metal−ligand attraction is slightly dominated by ΔEorb (52− 60%) than ΔEele. Four types of molecular orbitals, which belong to irreducible representations a1, a2, b1, and b2 of the C2v point group, contribute to the ΔEorb, which are designated as ΔEa1, ΔEa2, ΔEb1, and ΔEb2, respectively. These orbitals are given in 7175

dx.doi.org/10.1021/om3007097 | Organometallics 2012, 31, 7171−7180

Organometallics

Article

Å3) for M = Mo (see also Table S3 in the Supporting Information). Figure 3a shows the contour line diagram of the Laplacian of electron density (∇2ρ(r)) in the plane of the metallacycle along with bond paths (dark green lines) and ring paths (pink lines) for the 1,3-MC-bonded vanadacyclobutadiene VCl2(C3H3). This plot is a representative of other complexes, as they all show very similar Laplacian distributions. Areas of relative electronic charge concentration (dotted green lines) and areas of relative charge depletion (solid blue lines) are seen in the Laplacian plot. In the case of α-carbon atoms of the metallacycle, the area of charge concentration is clearly aligned in the direction of the metal center, while in the case of the β-carbon, an appreciable amount of charge concentration in the direction of the metal is visible. All the ring paths show a flat curvature near the rcp, suggesting that the rcp point is a “catastrophe” point.38 The “catastrophic” character of the rcp emerges very clearly when the Laplacian distribution is plotted for a geometry of VCl2(C3H3) with slight contraction (9%) of the V−Cβ distance from the equilibrium geometry (Figure 3b). The V−Cβ distance-restricted VCl2(C3H3) is 8.5 kcal/mol higher in energy than the fully optimized geometry, and it showed a bcp near the ring center flanked by two very close lying rcps along with the bond path between the metal and Cβ. The charge concentration in the direction of the metal from Cβ is more in the restricted geometry than the equilibrium geometry. Also depicted in Figure 3c is the Laplacian plot for the molybdenacyclobutadiene complex MoCl3(C3H3), a representative example for a metallacyclobutadiene complex without a 1,3-MC bond (Mo−Cβ distance is 2.611 Å; Figure 1c). This plot is very helpful to further assess the 1,3-MC bond. In this plot, each ring path between the bcp and rcp appears as a straight line, and unlike MCB with a 1,3-MC bond, a flat curvature is not seen near the rcp. Most significantly, the charge concentration is aligned in the direction of the metal only from Cα atoms and not from Cβ, which indicates that Cβ is not bonded with the metal. Alkyne Metathesis Reaction Models. The MCB systems used in this study can be considered as formed from a [2 + 2] cycloaddition between acetylene and the alkylidyne complex MCln(CH). Hence, transition state modeling for the cycloaddition reaction can yield valuable insight into alkyne metathesis processes of group 4, 5, and 6 metals.35 The energy profiles of these reactions are presented in Figure 4. All of them are exothermic and give rise to the formation of 1,3-MCbonded MCB complexes. All attempts to optimize a reactant van der Waals complex failed because the optimization always gave the MCB complex. Failure to detect van der Waals complexes may be due to a poor description of dispersion effects by the BP86 functional. Since van der Waals complexes are typically stabilized by 1−3 kcal/mol, the reactions can be described as nearly barrierless (Figure 4). The exothermic values of the reactions are very high for group 4 and 5 metals; the lowest is 75.2 kcal/mol for the Zr complex, and the highest is 101.7 kcal/mol for the Hf complex. The exothermic values of Cr, Mo, and W complexes are 6.4, 17.7, and 27.3 kcal/mol, respectively. In 1993, Ziegler and co-workers reported the energy profile for the formation of a molybdenacyclobutadiene complex from MoCl3(CH) and acetylene using nonlocal DFT.35 They showed that the reaction was 16.7 kcal/mol exothermic, and the value obtained in the present study using more sophisticated DFT is only 1 kcal/mol higher. The high exothermic character of the reactions indicates that the loss of alkyne from the MCB (the cycloreversion reaction) will be very

kcal/mol), while the De of molybdenum and tungsten complexes are 82.4 and 108.8 kcal/mol, respectively. We have also used the ETS-NOCV analysis to quantify the orbital interaction energy of the low-lying agostic a1 and b1 orbitals given in Figure 2. This analysis showed that interaction due to b1 is negligible (−2 to −4 kcal/mol), while that due to a1 is significant. The contribution of a1 in the quartet scheme of ETS-NOCV is −10.4, −12.4, −13.9, −11.5, −17.2, −20.2, −16.7, −19.1, and −24.2 kcal/mol for MCB complexes of Ti, Zr, Hf, V, Nb, Ta, Cr, Mo, and W, respectively. In Figure S2 in the Supporting Information, NOCV deformation density maps of all the MOs for the metal−carbon bonding are given for all the MCln(C3H3) complexes. They show that the 1,3-MC bond is largely defined by the sum of the ΔEorb of the agostic a1 orbital and the π-type b2 orbital, and this is estimated to be −25.2, −29.6, −31.7, −27.9, −36.9, −41.5, −42.9, −47.9, and −54.5 kcal/mol for MCB complexes of Ti, Zr, Hf, V, Nb, Ta, Cr, Mo, and W, respectively. It may be noted that these numbers cannot be exclusively used for M−Cβ interactions alone because the corresponding orbitals are delocalized to the α-carbons too. We have also analyzed the MCB complexes in Figure 1a on the basis of their Mayer bond order62 values given in Table 2. Table 2. Mayer Bond Order Values Calculated at the BP86/ TZ2P+ Level for Metallacyclobutadiene Complexes Mayer bond order complex

M−Cα

M−Cβ

Cα−Cβ

TiCl(C3H3) ZrCl(C3H3) HfCl(C3H3) VCl2(C3H3) NbCl2(C3H3) TaCl2(C3H3) CrCl3(C3H3) MoCl3(C3H3) WCl3(C3H3)

1.33 1.40 1.42 1.30 1.40 1.44 1.27 1.34 1.42

0.31 0.33 0.39 0.31 0.35 0.43 0.26 0.32 0.37

1.14 1.13 1.11 1.16 1.12 1.09 1.18 1.15 1.11

In all the cases, M−Cα bonds show double-bond character (bond order 1.27−1.44), which can be attributed to the δ interaction. The bond order for the 1,3-MC bond is a minimum of 0.26 for M = Cr and a maximum of 0.43 for M = Ta. The bond order of each C−C bond is in the range 1.09−1.18, which is rather small considering that three π-electrons are available for conjugation within the C3H3 fragment. This diminished double-bond character of the C−C bonds can be attributed to the involvement of the carbon p orbitals in δ and 1,3-MCbonded interactions. The bond order analysis strongly supports that the diagram given in Scheme 1 is appropriate for representing MCB complexes. The π orbital interaction of the type b2, which represents an interaction between atoms at the 1 and 3 positions, is rare in chemistry. The nature of M−C bonding in MCB structures is also analyzed in terms of the topology of the electron density (ρ(r)). All the structures showed the expected (3, −1) bond critical points (bcp) for the M−Cα and Cα−Cβ bonds. The ρ(r) at the bcp of M−Cα bonds is in the range 0.126−0.176 e/Å3, while that of Cα−Cβ bonds is in the range 0.272−0.289 e/Å3. The MCBs have not shown a (3, −1) bcp between M and Cβ, whereas they showed a (3, +1) ring critical point (rcp) close to the midpoint region between M and Cβ. The ρ(r) at the rcp is the lowest (0.071 e/Å3) for M = Hf and the highest (0.090 e/ 7176

dx.doi.org/10.1021/om3007097 | Organometallics 2012, 31, 7171−7180

Organometallics

Article

Figure 3. Contour plots of the Laplacian ∇2ρ(r) of the electron density and molecular graphs. (a) C2v symmetric vanadacyclobutadiene VCl2(C3H3). (b) C2v symmetric VCl2(C3H3) with contracted V−Cβ distance (1.843 Å). (c) Molybdenacyclobutadiene MoCl3(C3H3) with C−C bond alternation and large Mo−Cβ distance (2.611 Å). Bond paths are given in dark green lines and ring paths in pink lines. ρ(r) values at the critical points are given in e/Å3.

difficult in the case of group 4 and 5 metals, while the reaction is feasible at room temperature for Cr and Mo. For W, cycloreversion could be achieved with slight tuning of the ligand environment. The low activation barrier observed for the cycloreversion step of the chromium MCB complex suggests the possibility of developing Cr-based systems for alkyne metathesis. Our results also firmly suggest a strong thermodynamic driving force for the formation of MCB complexes of group 4 and 5 metals, indicating that their alkylidyne complexes can be targeted for the development of alkyne metathesis catalysts. In fact, several good attempts have been made in the recent past to produce alkylidyne complexes capable of metathesis using group 4 and 5 metals.63−70 Proposal to Use Metallacyclobutadienes As Alkyne Metathesis Catalysts. A bottleneck in the development of alkyne metathesis catalysts is in synthesizing robust systems with M−C triple bonds. The inherent high reactivity of M−C triple bonds cannot be easily concealed by the modification of the stereoelectronic properties of the ligand environment. In this scenario, synthesis of catalysts without M−C triple bonds

Figure 4. Energetics obtained for the formation of metallacyclobutadiene complexes, through the alkyne metathesis pathway. Models for a representative case; the molybdenum complex is given. A van der Waals complex connecting to the TS could not be located, and the dotted lines indicate this. Bond lengths in Å, and energy in kcal/mol.

Scheme 2. General Mechanistic Scheme for Alkyne Metathesis Using Metallacyclobutadiene I

7177

dx.doi.org/10.1021/om3007097 | Organometallics 2012, 31, 7171−7180

Organometallics

Article

could be considered for alkyne metathesis.3,71−73 The high amount of energy released in the formation of MCB systems in the case of groups 4 and 5 strongly suggests the possibility of the development of such stable complexes. We have considered the use of 1,3-MC-bonded MCBs as potential alkyne metathesis catalysts and propose a general mechanism given in Scheme 2, which is different from Katz's mechanism for alkyne metathesis. At first the alkyne can be oxidatively added to the M−Cα double bond of the MCB to produce a bicyclic intermediate (III). Reductive elimination of the new alkyne can be expected in the next step, which would liberate the catalyst. Using the 1,3-MC-bonded dichloro derivatives of V, Nb, and Ta complexes, the mechanism given in Scheme 2 is explored using the BP86/def2-TZVPP method, and the results are summarized in Figure 5. In the case of vanadium, coordination

metathesis was described on the basis of the loss of alkyne from the MCB to give an alkylidyne complex and regeneration of the MCB by further reaction of the alkylidyne and the alkyne. Essentially Katz's mechanism was followed in that reaction. The theoretical study by Ziegler et al.35 was also on Katz's mechanism, and through transition state modeling they reported an exothermic energy profile (16.7 kcal/mol) for the formation of molybdenacyclobutadiene. Schrock et al. also showed that the reactions of alkyne with alkoxy- and chloroligated tungstenacyclobutadienes were possible, but they yielded only Cp complexes and not alkyne metathesis.24,25 So far, alkyne metathesis is described via Katz's mechanism, and the alternate possibility presented in Scheme 2 is far from reality. The strategy of using the MCB complex as a robust alkyne metathesis catalyst may become successful if the ring closing of the organic fragment leading to the Cp complex is prevented. This may be achieved by steric tuning of the substituents on the organic fragment as well as by steric and electronic tuning of other ligands on the metal center.



CONCLUSIONS MCBs are intermediates of alkyne metathesis in the Katz mechanism, and they possess a characteristic M−Cβ bond distance, ∼2.1 Å for group 4, 5, and 6 metals. The β-carbon of the MCB exists in a unique planar tetracoordinate geometry due to significant bonding interactions (1,3-MC bond) between the metal and the β-carbon. The dπ−pπ orbital interaction between the metal and the β-carbon (b2 orbital) is mainly responsible for the 1,3-MC bond, and such a bond formed between atoms in the 1 and 3 positions is rare in chemistry. Very recently, Roy et al. have shown that 1,3-MC interactions can exist in metallacycloallene systems.74 The MCBs are also characterized by strong δ-bonding interactions (a2 orbital) between a metal d orbital and pπ orbitals on the α-carbon atoms. The orbital interactions due to b2 and a2 orbitals significantly reduce the π-bonding interaction within the C3H3 fragment, and as a result, both the Cα−Cβ bonds exist as formal Csp2−Csp2 single bonds. Hence it can be assessed that metal− carbon δ- and π-bonding interactions have a direct influence on the activation of C−C bonds to induce cycloreversion leading to loss of alkyne from the MCB. Formation of MCBs from metal alkylidyne complexes and alkynes through [2 + 2] cycloaddition reveals a barrierless and exothermic process. High exothermic values obtained for group 4 and group 5 MCBs suggest the possibility of isolating these complexes in an appropriate ligand environment. The disproportionation reaction or cycloreversion of the alkyne part from the MCB is very difficult to achieve for group 4 and group 5 complexes, and therefore the Katz mechanism of metathesis may not occur. Moderate exothermic values obtained for Cr, Mo, and W MCB complexes suggest the feasibility of metathesis through the Katz mechanism at ambient conditions, which also supports the rich metathesis chemistry associated with Mo and W alkylidyne complexes. Chromium is also an attractive choice for alkyne metathesis owing to the low activation energy required for both cycloaddition and cycloreversion steps. The use of MCBs as C−C triple-bond metathesis catalysts is investigated by using a mechanistic scheme very different from that of Katz's, and it reveals that a possible pathway exists for an MCB complex to carry out C−C triple-bond metathesis. However, this pathway becomes effective only if another dominant pathway leading to the formation of a cyclo-

Figure 5. Energetics of alkyne metathesis using vanada-, nioba-, and tantalacyclobutadienes. Blue, red, and green lines and values represent V, Nb, and Ta systems, respectively. Structures of vanadium complexes are depicted. Bond lengths in Å, and energy in kcal/mol.

of alkyne needs 24.1 kcal/mol of activation energy (TS1a), while the subsequent oxidative addition of the alkyne leading to the formation of a bicyclic product requires 17.6 kcal/mol of activation energy (TS2a). On the other hand, both Nb and Ta showed much lower values of activation energy (11.0−12.3 kcal/mol) for both acetylene coordination and oxidative addition. In all three cases, the reaction leads to the formation of a highly exothermic bicyclic product. The exothermic values are 22.2, 28.7, and 27.9 kcal/mol respectively for M = V, Nb, and Ta. Therefore, the reverse reaction, the reductive elimination of alkyne via TS2b, will be very difficult due to the high activation energy barrier (39.8−40.6 kcal/mol). Although a pathway for alkyne metathesis through a 1,3-MCbonded metallacyclobutadiene as active catalyst exists, to achieve realistic models of MCB-based metathesis, appropriate ligation as well as good tuning of stereoelectronic properties of the ligands has to be considered. Further, the bicyclic product has to be prevented from undergoing a rearrangement to a more stable cyclopentadienyl complex (MCln(C5H5)). The Cp complex is more stable than the bicyclic complex by 56.4, 36.8, and 31.2 kcal/mol for M = V, Nb, and Ta, respectively (details of this mechanism are given in Figure S3 in the Supporting Information for a representative case, the M = Nb system). In 1984 Schrock et al. published a paper in Organometallics with the title “Metathesis of Acetylenes by Triphenoxytungstenacyclobutadiene Complexes and the Crystal Structure of W(C3Et3)[O-2,6-C3H3(i-Pr)2]3”.33 In this paper, CC triple-bond 7178

dx.doi.org/10.1021/om3007097 | Organometallics 2012, 31, 7171−7180

Organometallics

Article

(30) Beer, S.; Hrib, C. G.; Jones, P. G.; Brandhorst, K.; Grunenberg, J.; Tamm, M. Angew. Chem., Int. Ed. 2007, 46, 8890. (31) Beer, S.; Brandhorst, K.; Hrib, C. G.; Wu, X.; Haberlag, B.; Grunenberg, J.; Jones, P. G.; Tamm, M. Organometallics 2009, 28, 1534. (32) Geyer, A. M.; Wiedner, E. S.; Gary, J. B.; Gdula, R. L.; Kuhlmann, N. C.; Johnson, M. J. A.; Dunietz, B. D.; Kampf, J. W. J. Am. Chem. Soc. 2007, 129, 3800. (33) Churchill, M. R.; Ziller, J. W.; Freudenberger, J. H.; Schrock, R. R. Organometallics 1984, 3, 1554. (34) Bursten, B. E. J. Am. Chem. Soc. 1983, 105, 121. (35) Woo, T.; Folga, E.; Ziegler, T. Organometallics 1993, 12, 1289. (36) Lin, Z. Y.; Hall, M. B. Organometallics 1994, 13, 2878. (37) Anslyn, E. V.; Brusich, M. J.; Goddard, W. A. Organometallics 1988, 7, 98. (38) Suresh, C. H.; Frenking, G. Organometallics 2010, 29, 4766. (39) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Montgomery, J. A., Jr.; Vreven, T.; Kudin, K. N.; Burant, J. C.; Millam, J. M.; Iyengar, S. S.; Tomasi, J.; Barone, V.; Mennucci, B.; Cossi, M.; Scalmani, G.; Rega, N.; Petersson, G. A.; Nakatsuji, H.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Klene, M.; Li, X.; Knox, J. E.; Hratchian, H. P.; Cross, J. B.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Ayala, P. Y.; Morokuma, K.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Zakrzewski, V. G.; Dapprich, S.; Daniels, A. D.; Strain, M. C.; Farkas, O.; Malick, D. K.; Rabuck, A. D.; Raghavachari, K.; Foresman, J. B.; Ortiz, J. V.; Cui, Q.; Baboul, A. G.; Clifford, S.; Cioslowski, J.; Stefanov, B. B.; Liu, G.; Liashenko, A.; Piskorz, P.; Komaromi, I.; Martin, R. L.; Fox, D. J.; Keith, T.; Al-Laham, M. A.; Peng, C. Y.; Nanayakkara, A.; Challacombe, M.; Gill, P. M. W.; Johnson, B.; Chen, W.; Wong, M. W.; Gonzalez, C.; Pople, J. A. Gaussian 03, revision E.01; Gaussian, Inc.: Wallingford, CT, 2004. (40) Becke, A. D. Phys. Rev. A 1988, 38, 3098. (41) Perdew, J. P. Phys. Rev. B 1986, 33, 8822. (42) Schäfer, A.; Huber, C.; Ahlrichs, R. J. Chem. Phys. 1994, 100, 5829. (43) Bickelhaupt, F. M.; Baerends, E. J. Rev. Comput. Chem. 2000, 15, 1. (44) Van Lenthe, E.; Baerends, E. J.; Snijders, J. G. J. Chem. Phys. 1993, 99, 4597. (45) Van Lenthe, E.; Baerends, E. J.; Snijders, J. G. J. Chem. Phys. 1994, 101, 9783. (46) Morokuma, K. J. Chem. Phys. 1971, 55, 1236. (47) Ziegler, T.; Rauk, A. Theor. Chem. Acc. 1977, 46, 1. (48) Ziegler, T.; Rauk, A. Inorg. Chem. 1979, 18, 1755. (49) Michalak, A.; Mitoraj, M.; Ziegler, T. J. Phys. Chem. A 2008, 112, 1933. (50) Mitoraj, M. P.; Michalak, A.; Ziegler, T. J. Chem. Theor. Comput. 2009, 5, 962. (51) Frenking, G.; Frohlich, N. Chem. Rev. 2000, 100, 717. (52) Frenking, G.; Wichmann, K.; Fröhlich, N.; Loschen, C.; Lein, M.; Frunzke, J.; Rayón, V. M. Coord. Chem. Rev. 2003, 55, 238. (53) Hopffgarten, M. v.; Frenking, G. WIREs Comput. Mol. Sci. 2012, 2, 43. (54) Frenking, G.; Frohlich, N. Chem. Rev. 2000, 100, 717. (55) Mitoraj, M. P. J. Phys. Chem. A 2011, 115, 14708. (56) Bader, R. F. W. Atoms in Molecules: A Quantum Theory; Oxford University Press: Oxford, 1990. (57) Todd, A.; Keith, T. K. AIMAll (version 12.06.03); Gristmill Software: Overland Park, KS, USA, 2012. (58) Blunden, R. B.; Cloke, F. G. N.; Hitchcock, P. B.; Scott, P. Organometallics 1994, 13, 2917. (59) Suresh, C. H. J. Organomet. Chem. 2006, 691, 5366. (60) Suresh, C. H.; Baik, M. H. Dalton Trans. 2005, 2982. (61) Suresh, C. H.; Koga, N. Organometallics 2004, 23, 76. (62) Mayer, I. Chem. Phys. Lett. 1983, 97, 270. (63) Herndon, J. W. Coord. Chem. Rev. 2009, 253, 86.

pentadienyl complex is prevented. We hope that this can be achieved by proper tuning of the steric and electronic environment of the ligands.



ASSOCIATED CONTENT

S Supporting Information *

Optimized geometries and results of EDA and electron density analyses. This material is available free of charge via the Internet at http://pubs.acs.org.

■ ■

AUTHOR INFORMATION

Notes

The authors declare no competing financial interest.

ACKNOWLEDGMENTS C.H.S. is grateful to the Alexander von Humboldt Foundation for the award of the fellowship for experienced researchers. Excellent service by the computer center of the Philipps Universität Marburg is gratefully acknowledged.



REFERENCES

(1) Schrock, R. R. Chem. Rev. 2002, 102, 145. (2) Fürstner, A.; Davies, P. W. Chem. Commun. 2005, 2307. (3) Mortreux, A.; Coutelier, O. J. Mol. Catal. A 2006, 254, 96. (4) Mindiola, D. J. Acc. Chem. Res. 2006, 39, 813. (5) Schrock, R. R.; Czekelius, C. Adv. Synth. Catal. 2007, 349, 55. (6) Zhang, W.; Moore, J. S. Adv. Synth. Catal. 2007, 349, 93. (7) Wu, X. A.; Tamm, M. Beilstein J. Org. Chem. 2011, 7, 82. (8) Xue, Z.-L.; Morton, L. A. J. Organomet. Chem. 2011, 696, 3924. (9) Fürstner, A. Chem. Commun. 2011, 47, 6505. (10) Tonzetich, Z. J.; Lam, Y. C.; Müller, P.; Schrock, R. R. Organometallics 2007, 26, 475. (11) Schrock, R. R. In Handbook of Metathesis; Grubbs, R. H., Ed.; Wiley VCH: Weinheim, 2003; Vol. 1. (12) Mortreux, A.; Blanchard, M. J. Chem. Soc., Chem. Commun. 1974, 787. (13) Sancho, J.; Schrock, R. R. J. Mol. Catal. 1982, 15, 75. (14) Schrock, R. R.; Clark, D. N.; Sancho, J.; Wengrovius, J. H.; Rocklage, S. M.; Pederson, S. F. Organometallics 1982, 1, 1645. (15) Heppekausen, J.; Stade, R.; Goddard, R.; Fuerstner, A. J. Am. Chem. Soc. 2010, 132, 11045. (16) Wiedner, E. S.; Gallagher, K. J.; Johnson, M. J. A.; Kampf, J. W. Inorg. Chem. 2011, 50, 5936. (17) Heppekausen, J.; Fuerstner, A. Angew. Chem., Int. Ed. 2011, 50, 7829. (18) Wu, X.; Daniliuc, C. G.; Hrib, C. G.; Tamm, M. J. Organomet. Chem. 2011, 696, 4147. (19) Lysenko, S.; Haberlag, B.; Daniliuc, C. G.; Jones, P. G.; Tamm, M. Chemcatchem 2011, 3, 115. (20) Lichtscheidl, A. G.; Ng, V. W. L.; Mueller, P.; Takase, M. K.; Schrock, R. R.; Malcolmson, S. J.; Meek, S. J.; Li, B.; Kiesewetter, E. T.; Hoveyda, A. H. Organometallics 2012, 31, 4558. (21) Sarkar, S.; McGowan, K. P.; Kuppuswamy, S.; Ghiviriga, I.; Abboud, K. A.; Veige, A. S. J. Am. Chem. Soc. 2012, 134, 4509. (22) Freudenberger, J. H.; Schrock, R. R.; Churchill, M. R.; Rheingold, A. L.; Ziller, J. W. Organometallics 1984, 3, 1563. (23) Katz, T. J.; McGinnis, J. J. Am. Chem. Soc. 1975, 97, 1592. (24) Schrock, R. R.; Pedersen, S. F.; Churchill, M. R.; Ziller, J. W. Organometallics 1984, 3, 1574. (25) Schrock, R. R. Polyhedron 1995, 14, 3177. (26) Wengrovius, J. H.; Sancho, J.; Schrock, R. R. J. Am. Chem. Soc. 1981, 103, 3932. (27) Pedersen, S. F.; Schrock, R. R.; Churchill, M. R.; Wasserman, H. J. J. Am. Chem. Soc. 1982, 104, 6808. (28) Schrock, R. R.; Weinstock, I. A.; Horton, A. D.; Liu, A. H.; Schofield, M. H. J. Am. Chem. Soc. 1988, 110, 2686. (29) Zhu, J.; Jia, G. C.; Lin, Z. Y. Organometallics 2006, 25, 1812. 7179

dx.doi.org/10.1021/om3007097 | Organometallics 2012, 31, 7171−7180

Organometallics

Article

(64) Mindiola, D. J.; Bailey, B. C.; Basuli, F. Eur. J. Inorg. Chem. 2006, 3135. (65) Basuli, F.; Bailey, B. C.; Brown, D.; Tomaszewski, J.; Huffman, J. C.; Baik, M. H.; Mindiola, D. J. J. Am. Chem. Soc. 2004, 126, 10506. (66) Bailey, B. C.; Fan, H.; Huffman, J. C.; Baik, M.-H.; Mindiola, D. J. J. Am. Chem. Soc. 2006, 128, 6798. (67) Flores, J. A.; Cavaliere, V. N.; Buck, D.; Pinter, B.; Chen, G.; Crestani, M. G.; Baik, M. H.; Mindiola, D. J. Chem. Sci 2011, 2, 1457. (68) Fout, A. R.; Bailey, B. C.; Buck, D. M.; Tan, H. J.; Huffman, J. C.; Baik, M. H.; Mindiola, D. J. Organometallics 2010, 29, 5409. (69) Scott, J.; Mindiola, D. J. Dalton Trans. 2009, 8463. (70) Carnahan, E. M.; Lippard, S. J. J. Am. Chem. Soc. 1992, 114, 4166. (71) Bogoslavsky, B.; Levy, O.; Kotlyar, A.; Salem, M.; Gelman, F.; Bino, A. Angew. Chem., Int. Ed. 2012, 51, 90. (72) Bino, A.; Ardon, M; Shirman, E. Science 2005, 308, 234. (73) Bunz, U. H. F. Science 2005, 308, 216. (74) Roy, S.; Jemmis, E. D.; Schulz, A.; Beweries, T.; Rosenthal, U. Angew. Chem., Int. Ed. 2012, 51, 5347.

7180

dx.doi.org/10.1021/om3007097 | Organometallics 2012, 31, 7171−7180