Computational Design of an Iron Catalyst for Olefin Metathesis

1 hour ago - ... considering two possible deformation pathways: the transformation of the carbene ligand to a coordinated alkene molecule and the form...
1 downloads 0 Views 2MB Size
Article Cite This: Organometallics XXXX, XXX, XXX−XXX

pubs.acs.org/Organometallics

Computational Design of an Iron Catalyst for Olefin Metathesis Bo Yang* and Donald G. Truhlar* Department of Chemistry, Inorganometallic Catalyst Design Center, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States

Downloaded via UNIV OF NEW ENGLAND on October 27, 2018 at 01:56:03 (UTC). See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles.

S Supporting Information *

ABSTRACT: Olefin metathesis is a versatile reaction in synthetic chemistry. The use of iron compounds as olefin metathesis catalysts has attracted interest but has not succeeded because most Fe complexes catalyze cyclopropanation reactions preferentially to cycloreversion, which is a key step in the metathesis reaction cycle. In the present study, we performed density functional theory calculations to study the catalytic effect of ancillary pincer ligands in an Fe organometallic catalyst for olefin metathesis. We show here that the partial atomic charge on Fe in the metallocyclobutane intermediate is a descriptor for the tendency of the Fe complex to favor olefin metathesis over cyclopropanation. The results suggest that the pincer ligands can decrease the partial atomic charge on the Fe ion and that decreasing this charge can make the desired cycloreversion reaction more energetically favorable than the cyclopropanation side reaction. In this way, we found that the Fe carbene stabilized by Nheterocyclic dicarbene ligands (Fe(C−N−C)(CHCH3), where C−N−C denotes 2,6-bis(methylimidazol-2-ylidene)pyridine), is a potential candidate for catalyzing the olefin metathesis reaction. We then calculated the whole catalytic cycle for the propylene metathesis reaction catalyzed by the Fe(C−N−C) carbene, and the results indicate that the cyclopropanation reaction is disfavored by ∼4 kcal/mol. We verified the stability of the proposed Fe(C−N−C) catalyst by considering two possible deformation pathways: the transformation of the carbene ligand to a coordinated alkene molecule and the formation of a hydride carbyne intermediate followed by the insertion of a propylene molecule. We found that neither deformation pathway is energetically favorable. Finally we propose a viable route to synthesize the Fe(C−N−C) carbene complex by utilizing the already available Fe(N2)2 complex stabilized by an N-heterocyclic dicarbene ligand.



INTRODUCTION Olefin metathesis is a versatile reaction that has extensive applications ranging from large-scale industrial processes to fine chemical manufacturing.1−5 It has been applied to a wide variety of substrates in the synthesis of natural products6−8 and polymers9,10 and has also found applications in peptide and protein modifications.11−14 New C−C double bonds can be formed efficiently through transition-metal-catalyzed olefin metathesis,15−19 as illustrated in Scheme 1. The generally accepted Chauvin mechanism20 is shown in Scheme 2; it involves the repetition of three steps: the coordination of an alkene molecule to a transition metal carbene complex, [2 + 2] cycloaddition to give a metallocyclobutane as an intermediate, and [2 + 2] cycloreversion to

produce a new alkene molecule and regenerate the transition metal carbene complex. Mo-, W-, Re-, and Ru-based homogeneous metal carbene complexes are well-known as active olefin metathesis catalysts.21−26 The 3d transition metals, though desirable for their abundance, are rarely utilized for metathesis because they often lead to cyclopropanation rather than cycloreversion; these two outcomes are contrasted in Scheme 3.27,28 In some cases, the side reaction is favored due to 3d transition metal complexes having low-energy highspin states that can favor cyclopropanation by a radical mechanism.29−36 One of the biggest challenges for the utilization of 3d transition metals is to improve the selectivity for metathesis in competition with the cyclopropanation side reaction. It has proven difficult to find a successful Fe-based olefin metathesis catalyst,37−41 so we embarked on an attempt to use quantum chemistry for the rational design of such a catalyst. It is known that fine-tuning of the metal coordination environment in homogeneous organometallic complexes can result in high olefin metathesis activities and selectivities,27,28,42,43 and our design of new catalysts starts with the establishment of

Scheme 1. Olefin Metathesis Reactiona

Received: August 14, 2018

a

L1 and L2 denote substituents. © XXXX American Chemical Society

A

DOI: 10.1021/acs.organomet.8b00583 Organometallics XXXX, XXX, XXX−XXX

Article

Organometallics

Scheme 3) utilizes the metallocyclobutane intermediate to produce an alkene as the product and Fe alkylidene as the catalytically active complex. The cyclopropanation reaction (R2 in Scheme 3) proceeds through the metallocyclobutane as an intermediate and generates a cyclopropane as the side product. On the basis of our calculations, the cyclopropanation reaction is always more exergonic than the cycloreversion reaction, and this is the case not only for Fe catalysts but also for those that utilize a Ru atom, such as the Grubbs II catalyst.27 Nevertheless, the cyclopropanation reaction can be prevented when its free energy of activation is higher than that of the cycloreversion reaction. It has been demonstrated that tuning the electronic and steric properties of the ancillary ligand of an Fe metallocyclobutane complex can alter its activity toward cycloreversion and cyclopropanation reactions.27,28 Here we considered three pincer ligands. The structures of the pincer ligands with coordinated Fe atoms are shown in Scheme 4.

Scheme 2. Chauvin Mechanism for the Olefin Metathesis Reactiona

Scheme 4. Examples of Pincer-Ligand-Stabilized Fe Complexes

a

M denotes a transition metal atom; L1 and L2 denote substituents.

Scheme 3. Cycloreversion (R1) and Cyclopropanation (R2)

structure−reactivity relationships27,44 that can help us to finetune ligands. We define the theoretical selectivity in terms of transition state theory activation parameters as S = exp{[ΔG⧧(R2) − ΔG⧧(R1)] /RT}

They all involve a pyridine with two identical “arms” ortho to the N atom, and these arms anchor the Fe active site. The atoms that are bonded to the Fe ion are N−N−N in the left structure of Scheme 4, P−N−P in the middle structure, and C−N−C on the right. These ligands modify the electronic characteristics of the Fe center and consequently change its reactivity toward the cycloreversion and cyclopropanation reactions. To further fine-tune the properties of the Fe active site, various substitutions on the three pincer ligands were also tested. The geometries of all of the tested pincer structures are provided in the Supporting Information (SI). Both singlet and triplet spin states are considered for the reactants and transition states of R1 and R2. The quintet spin state was found to be at least 7 kcal/mol higher than the singlet and triplet states for several tested examples, and it was not considered further. We found that for seven of the 17 tested complexes, namely, complexes S3 to S8 and S16 in the SI, the lowest-energy reactant is a singlet. For the cycloreversion reaction (R1), we found that the lowest-energy transition state is always the singlet. For the cyclopropanation reaction (R2), the lowest-energy transition state of 14 of the 17 tested complexes is a triplet, with the exceptions being complex (3) in Scheme 4 and complexes S1 and S8 in the SI. This spin-state propensity for the transition states agrees with previous results on similar systems.27,28 For each complex, the free energies of activation for R1 and R2 were calculated from the free energies of the most stable reactant structure and the lowest-energy transition structure of each reaction. This includes the situation in which the reactant

(1)

which involves the difference in the free energies of activation of reactions R2 and R1, the gas constant R, and the temperature T. Motivated by this expression, we define the free energy selectivity as δ ≡ ΔΔG⧧ ≡ ΔG⧧(R2) − ΔG⧧(R1)

(2)

Tridentate pincer ligands and the resulting Fe(0) carbene complexes were shown to stabilize singlet metallocyclobutane intermediates, and this was hypothesized to hinder the cyclopropanation side reaction.27,45,46 Therefore, we chose to study tridentate pincer ligands, and the plan of this article is as follows. We begin by computing δ for 17 such ligands. We then use these calculations to find a suitable descriptor for the desired selectivity, and we use the descriptor to guide the choice of ligand. We then confirm the predicted selectivity by calculating the whole catalytic cycle for the propylene metathesis reaction over the Fe(C−N−C) carbene catalyst. Finally we examine the stability and possible synthesis of the proposed new catalyst.



RESULTS AND DISCUSSION Cycloreversion versus Cyclopropanation: Effect of the Ancillary Ligand. As can be seen in the Chauvin mechanism (Scheme 2), the cycloreversion reaction (R1 in B

DOI: 10.1021/acs.organomet.8b00583 Organometallics XXXX, XXX, XXX−XXX

Article

Organometallics

where μ is the electronic chemical potential and η is the chemical hardness. For the triplet metallacyclobutanes, we evaluate ϵHOMO as the energy of the highest singly occupied molecular orbital, and for the singlets, it is the energy of the highest doubly occupied orbital. It should be noted that even though the HOMO and LUMO of the metallocyclobutane complex are mainly Fe d orbitals, ω is a property of the whole complex and is not associated with any individual atom. The top panel of Chart 2 shows that the electrophilicity of the metallacyclobutane is correlated with the partial atomic

and transition state possess different spin states, and we assume that the spin coupling is strong enough to allow a transition between the two spin states. The quantitative effect of the spin−orbit coupling on the barrier to reaction could be characterized by the TSSMM program,47 which mixes the two spin states through spin−orbit coupling. However, the spin− orbit coupling is small enough for Fe that refinement of the calculation in that way would not make a significant enough change in the kinetics to change any of our conclusions, so this refinement was not undertaken. The charge distributions were characterized by using charge model 5 (CM548). Chart 1 shows the free energy selectivity as

Chart 2. (top) Electrophilicity Index, ω, as a Function of the CM5 Partial Atomic Charge of the Fe Ion in the Metallocyclobutane Complex, Q(Fe); (bottom) Free Energy Selectivity, δ, Computed as the Difference in the Free Energies of Activation of Reactions R2 and R1, as a Function of the Electrophilicity Index, ωa

Chart 1. Free Energy Selectivity Defined by Equation 2, δ, as a Function of the CM5 Partial Charge of the Fe Ion of the Metallocyclobutane Complex, Q(Fe)a

a

Each data point represents an Fe complex with a certain pincer ligand; three of the 17 ligands (the ones identified by labels with arrows) are shown in Scheme 4, and the others are shown in section S2 in the SI.

a function of the CM5 partial charge of the Fe atom in the metallocyclobutane complex. Since the cyclopropanation reaction formally reduces the Fe center, a less positively charged Fe atom would be expected to make the cyclopropanation reaction less favorable, and the chart shows a clear trend in which the cycloreversion becomes energetically more favorable than the cyclopropanation reaction when the charge of the Fe ion decreases. When the Fe partial charge is lower than 0.4, reaction R2 has a higher free energy of activation than R1 and is accordingly predicted to have a lower reaction rate. The correlation of the free energy selectivity with the partial atomic charge on the Fe ion can be rationalized by noting that the electrophilicity is also correlated with this partial charge. We can make this connection by using the electrophilicity index ω of Parr et al.,49 which has been used to rationalize chemical reactivity in a variety of contexts in the past.50−52 Parr et al. defined ω in terms of the energetic stabilization upon acquiring an electron from an idealized electron sea. They showed that this may be related to the HOMO and LUMO energies (ϵHOMO and ϵLUMO) of a density functional theory calculation as ω = μ2 /(2η)

(3)

μ ≈ (ϵHOMO + ϵLUMO)/2

(4)

η ≈ ϵLUMO − ϵHOMO

(5)

a

In both panels, each data point is for a different ligand environment of the Fe ion.

charge on the Fe ion. Combining this result with that in Chart 1 implies that the selectivity is correlated with the electrophilicity, and this correlation is confirmed in the bottom panel of Chart 2. Our results show that the electrophilicity of the complex increases as the partial atomic charge of the Fe ion increases, and when the Fe complex is less electrophilic, it tends to favor the cycloreversion reaction instead of the cyclopropanation side reaction. Our screened Fe complexes give calculated ω values ranging from 0.8 to 18.0 eV, indicating medium to strong electrophilic character of the complex. Propylene Metathesis over Fe(C−N−C) Carbene Complexes. We then tested the Fe complex with the C− C

DOI: 10.1021/acs.organomet.8b00583 Organometallics XXXX, XXX, XXX−XXX

Article

Organometallics N−C pincer ligand (3), where C−N−C denotes 2,6bis(methylimidazol-2-ylidene)pyridine, for the propylene metathesis reaction. Complex (3) was chosen for its synthetic feasibility and because it is predicted to disfavor the cyclopropanation side reaction according to the results shown in Chart 1. The results are discussed in the following sections, focusing on the activity, selectivity, and stability of the corresponding Fe carbene complex. The optimized structure of the complex is provided in Scheme 5. A complete cycle for the propylene metathesis

Chart 3. Free Energy Profiles for Propylene Metathesis, Including Cycloaddition and Cycloreversion Steps (black) and Cyclopropanation Reactions (red), over the Fe(C−N− C) Complexa

Scheme 5. Structure of the Bare Fe(C−N−C) Complexa

a

C atoms are shown in gray, H in white, N in blue, and Fe in purple.

Scheme 6. Propylene Metathesis Reaction To Give Ethylene and 2-Butene as Products

(reaction R3, shown in Scheme 6) was constructed on the basis of the Chauvin mechanism using the Fe(C−N−C) carbene complexes as catalytically active species (Scheme 7);

a

The catalytic cycle as shown in Scheme 7 is divided into two diagrams; each includes one cycloaddition reaction and one cycloreversion reaction. Results for the singlet (solid lines) and triplet (dashed lines) spin states are presented in the diagrams. Quintet states have higher energies and for simplicity are not presented in the diagram. “Ln−Fe” represents an Fe atom attached to the C−N−C pincer ligand as shown in Scheme 5.

Scheme 7. Catalytic Cycle for Propylene Metathesis (R3) Based on the Chauvin Mechanisma

tanes, and the bare complex (Ln−Fe as shown in the diagrams) favor a triplet spin state. According to our calculations (as explained in Computational Details, except where explicitly stated otherwise, these were carried out with the OPBE53−56 density functional), the alkene metathesis reactions, including cycloadditions and cycloreversions, have free energies of activation between 7.8 and 11.6 kcal/mol. The transition states for the cycloreversions are ∼4 kcal/mol more stable than the transition states for the cyclopropanation reactions throughout the reaction cycle, which predicts that the cyclopropanation side reaction is not favored over the olefin metathesis reactions. For the proposed Fe(C−N−C) metathesis system, we observed that the transition states for the cycloaddition and cycloreversion reactions have structures in which the dihedral angle α, which is illustrated in Chart 4, is about 120° on average. Similar transition structures are also obtained for the Grubbs II catalyst using the OPBE functional. (Detailed results for the Grubbs II catalyst will be discussed below.) We also noticed that all of the metallocyclobutane intermediates have distorted tetragonal pyramid (4P) structures, as shown in Chart 4, and that the Fe−C−C−C four-membered rings are nonplanar in all cases with a dihedral angle β of about 160°. The observed 4P structure is different from the common

“Ln−Fe” represents an Fe atom attached to the C−N−C pincer ligand as shown in Scheme 5.

a

these complexes are denoted as Ln−Fe in the upcoming diagrams. Free energy profiles based on the proposed reaction cycle were computed, and they are presented in Chart 3. Detailed energies for singlet, triplet, and quintet states are provided in the SI. For the Fe(C−N−C) carbene catalytic system, the transition states for cycloaddition, cycloreversion, and cyclopropanation are in their singlet spin states. On the other hand, the stable intermediates, including Fe methylidene, Fe ethylidene, mono- and dimethyl-substituted metallocyclobuD

DOI: 10.1021/acs.organomet.8b00583 Organometallics XXXX, XXX, XXX−XXX

Article

Organometallics Chart 4. (top left) Transition Structure for the Cycloaddition Reaction in Which a Propylene Reacts with the Fe(C−N−C) Methylidene Intermediate and Gives Metallocyclobutane (top right) as the Product; (bottom) Illustration of Tetragonal Pyramid (4P) and Trigonal Bipyramid (3P) Structuresa

Chart 5. Comparison of Free Energies of Activation for Cycloreversion and Cyclopropanation Reactions over the Fe(C−N−C) Catalyst Calculated Using Different Density Functionalsa

a

The cycloreversion and cyclopropanation reactions refer to the those in the upper panel of Chart 3.

for the cyclopropanation reaction. For all of the tested functionals, the cycloreversion reaction favors the singlet transition state over the triplet. Five of the six functionals (the exception is the OPBE functional) predict the triplet to be the most stable spin state for the transition state of the cyclopropanation reaction. Three of the six functionals, namely, OPBE, M06-L, and revM06-L, suggest a clear disfavoring of the cyclopropanation reaction as opposed to the cycloreversion reaction. However, the MN15 and MN15-L functionals give similar reaction barriers for the two reactions, and the M06 functional predicts the cyclopropanation reaction to be energetically more favorable than the cycloreversion reaction. For the M06 functional, the calculated free energies of activation are 23.8 and 15.0 kcal/mol for the cycloaddition and the cyclopropanation reactions, respectively. We conclude that our predictions have uncertainty due to the difficulty of predicting relative reaction rates for complicated transition metal systems. This highlights the fact that at the present time, theory can identify trends and suggest possible candidates for improved catalysts, but the predictions are quantitatively uncertain. We have also compared the proposed Fe(C−N−C) catalyst with the well-established Grubbs II metathesis catalyst, which involves a saturated N-heterocyclic carbene and has Ru as the active center, as shown in Chart 6. Chart 7 shows the OPBEbased free energy profiles for the cycloaddition and cycloreversion reactions over the Fe(C−N−C) and Grubbs II catalysts; the proposed Fe(C−N−C) catalyst exhibits higher activity than the Grubbs II catalyst. The Ru transition states for the cycloaddition and cycloreversion reactions are predicted to be higher in energy than those for the corresponding Fe transition states by 6.6 and 10.4 kcal/mol, respectively. Considering the large variations observed for the computed barriers with different density functionals, we also calculated the free energies of activation for the cycloreversion reaction using three other functionals, namely, M06, M06-L, and MN15-L. The M06 functional was previously found to give the most accurate energetics for a model system of the Grubbs II catalyst in a benchmark study considering 39 density functionals.64 The M06-L functional was also recommended for its similar accuracy and high efficiency for calculations on large systems.64 The M06-L functional has also been applied to

a

Red dotted lines do not represent bonding interactions but have been added to define the dihedral angles and illustrate 3P and 4P structures.

observation that the singlet metallocyclobutane complex usually exhibits a trigonal bipyramid (3P) structure.27 A 3P structure is illustrated in Chart 4; it is unstable and reverts to the 4P structure during geometry optimization. The molecular orbitals of the optimized singlet metallocyclobutane are given in section S1. To test the sensitivity of the computed activation barrier to the choice of density functional, the free energies of activation for the first cycloreversion and cyclopropanation reactions (as shown in Chart 3 for the OPBE exchange−correlation functional) were calculated using five other density functionals: M06-L,57 revM06-L,58 MN15-L,59 M06,60 and MN15.61 The M06-L and MN15-L functionals show high accuracy for transition metal chemistry,59,62 and the revM06-L functional improves upon the M06-L functional and achieves higher accuracy in transition metal coordination reactions.58 The M06-L and M06 functionals also show good performance in predicting transition metal compound structures for olefin metathesis catalysts.63 We considered both singlet and triplet spin states for all of the intermediates and transition states. Free energies of activation were calculated on the basis of intermediates and transition states optimized using each functional and are shown in Chart 5. In general, the OPBE functional predicts lower barriers than the other five functionals; the computed free energies of activation for the cycloreversion reaction are 11.6, 13.8, 20.7, 19.3, 23.8, and 19.7 kcal/mol using the OPBE, M06-L, revM06-L, MN15-L, M06, and MN15 functionals, respectively. The same trend also holds E

DOI: 10.1021/acs.organomet.8b00583 Organometallics XXXX, XXX, XXX−XXX

Article

Organometallics

functional, however, gave a more complicated energy profile for the cycloaddition reaction; it yielded two transition states with a new local minimum between them. The first transition structure found with M06 has α = 118°, and the second has α = 170°; the respective free energies of activation are 0.8 and 3.3 kcal/mol. The local minimum that connects the two saddle points corresponds to a Ru ethylidene complex with an adsorbed propene molecule. For the cycloreversion reaction (the second reaction step shown in Chart 7), the OPBE functional gives a transition state with α = 120°, and the Minnesota functionals give transition states with α ≈ 170°. (For the cycloreversion reaction, only one transition state was found using the M06 functional, as opposed to the two transition states for the cycloaddition reaction). The free energy profile for the Grubbs II catalyst obtained using M06 (Chart S1) and detailed geometries and energies of the optimized intermediates and transition states for the Grubbs II catalyst are provided in the SI. For our tested Fe complexes, the transition structures for the cycloaddition and cycloreversion always have α ≈ 130°. Stability of the Fe(C−N−C)(CHCH3) Complex: Hydride-Elimination-Induced Side Reactions. The Fe(C−N−C)(CHCH3) complex can be regarded as a Fischer carbene67 with the Fe atom in the formal oxidation state of zero. Oxidative addition reactions, such as hydride elimination, can oxidize the Fe(0) center to give Fe(II) complexes. Because of the increase in the charge of the Fe center, such oxidation processes could render the catalyst inactive for olefin metathesis. In addition, hydride elimination also gives Fe hydride as the product, which can sequentially react with alkenes to produce unwanted alkyl groups. Here we evaluate two possible side reaction pathways induced by the hydride elimination reactions. It should be noted that both singlet and triplet spin states were considered for all of the intermediates involved in the side reactions, and the singlet was the more stable spin state in all cases. Scheme 8 illustrates the process in which Fe(C−N−C)( CHCH3) undergoes α-hydride elimination to produce the

Chart 6. Illustration of the Alkylidene Intermediate Reacting with an Alkene Molecule over the Grubbs II Catalysta

The dihedral angle α is the angle between the Ru−C(1)−C(2) plane involving the C(1)−C(2) double bond of the olefin molecule and the Ru−C(2)−C(3) plane involving the Ru−C(3) double bond of the carbene complex.

a

Chart 7. OPBE-Functional-Based Free Energy Profiles for the Cycloaddition and Cycloreversion Reactions over the Fe(C−N−C) Complex (black) and the Grubbs II Catalyst (blue)a

a

Both singlet and triplet spin states were considered for all of the intermediates, and the most stable states are shown in the diagram.

Scheme 8. α-H Elimination of the Fe(C−N−C)(CHCH3) Complex Followed by Insertion of a Propylene Molecule

describe the carbene rotation and non-covalent interactions in the Grubbs II catalyst activation process.65,66 Calculations based on the M06, M06-L, and MN15-L functionals all suggest that the proposed Fe(C−N−C) catalyst would have lower activity than the Grubbs II catalyst. The overall activation barriers over the Grubbs II catalyst considering both the cycloaddition and cycloreversion reactions are predicted to be 14.3, 2.5, and 8.6 kcal/mol lower in energy than those for the corresponding Fe transition states using the M06, M06-L, and MN15-L functionals, respectively. Here we will further discuss the first reaction step shown in Chart 7, namely, the cycloaddition reaction, in which the ethylidene complex reacts with a propylene molecule to give a metallocyclobutane as the product. We used the OPBE, M06, M06-L, and MN15-L density functionals and found that the transition structure obtained for the cycloaddition reaction depends strongly on which density functional is used. The structures can be distinguished by the dihedral angle α defined in Chart 6. The optimized transition structure obtained with the OPBE functional has α = 129°, while the transition structures optimized using the M06-L and MN15-L functionals have α = 170° and 169°, respectively. In these three cases, we verified that the transition structures connect to the same reactant and product in the cycloaddition reaction. The M06

hydride carbyne, Fe(C−N−C)(H)(CCH3), which in turn reacts with a propylene molecule to give the propyl carbyne, Fe(C−N−C)(C3H7)(CCH3), as the final side product. The corresponding free energy profile is provided in the top panel of Chart 8. Since this process competes with the cycloaddition reaction of propylene metathesis, the energy profile for cycloaddition is also included in Chart 8. On the basis of calculations using the OPBE functional, the α-hydride elimination has a free energy of activation of 2.5 kcal/mol, which is much lower than the 11.5 kcal/mol barrier of the cycloaddition reaction, and the subsequent insertion reaction has a barrier comparable to that of the cycloaddition. On the F

DOI: 10.1021/acs.organomet.8b00583 Organometallics XXXX, XXX, XXX−XXX

Article

Organometallics Chart 8. (red) Computed Free Energy Profiles for (top) Carbene α-H Elimination Followed by Insertion of an Alkene Molecule and (bottom) Carbene β-H Elimination To Produce a Hydride Vinyl Intermediate; (black) Computed Free Energy Profile for the Cycloaddition Reaction

Scheme 9. Experimentally Synthesized Fe(N2)2 Complex Stabilized by an N-Heterocyclic Dicarbene Ligand

Scheme 10. Reaction Scheme for Producing and Activating the Fe(C−N−C) Carbene Catalyst: (i) Carbene Formation; (ii) Alkene Association; (iii) Alkene Metathesisa

a

The C−N−C pincer ligand is not shown in the scheme.

used as a precursor for the synthesis of Fe(N2)L-type complexes with L = C2H4 or PMe3.68 To synthesize the desired Fe(C−N−C) carbene complex, it is necessary to replace the coordinated C2H4 molecule from the Fe(N2)(C2H4) complex with a terminal alkyne, which will produce Fe(N2)(terminal alkyne) as the key intermediate. From the coordinated terminal alkyne group, a carbene group can be formed through tautomerization of the coordinated terminal carbyne, which is a common reaction that has also been used for the preparation of Ru vinylidene complexes.70,71 The final product, Fe(N2)(carbene), can be used as the precatalyst for the metathesis reaction. To activate the Fe(N2)(carbene) precatalyst, an alkene molecule needs to be introduced to the Fe active site to replace the coordinated N2 molecule. The alkene molecule is preferred to take the liquid form under reaction conditions to enable easy separation of the dissociated gas-phase N2 molecule. The coordinated alkene can then react with the carbene group through alkene metathesis to produce the desired product. If the intended alkene reactant is a gas under the reaction conditions, it can be introduced to the catalytic system after the N2 molecule has been removed.

other hand, the alkyl product, Fe(C−N−C)(C3H7)(CCH3), is less stable than the cycloaddition product as well as all of the other considered stable intermediates within the metathesis reaction cycle. Accordingly, the cycloaddition reaction is thermodynamically more favored. β-Hydride elimination of the Fe(C−N−C)(CHCH3) complex could result in the formation of a coordinated alkene molecule. This process includes the formation of the vinyl complex, Fe(C−N−C)(H)(−CHCH2) followed by transfer of the hydride to the CH group. Compared with the competing cycloaddition reaction, the β-hydride elimination exhibits a free energy of activation that is 3.9 kcal/mol higher, as shown in the bottom panel of Chart 8. This indicates that the β-hydride elimination reaction is disfavored over the cycloaddition reaction. Accordingly, the subsequent insertion reaction is unlikely to take place because of the lack of starting intermediate. Proposed Route for Synthesizing the Fe(C−N−C) Carbene Complexes. A reaction route for the synthesis of the Fe(C−N−C) carbene molecule is proposed that utilizes the experimentally available (Fe(N2)2) complex stabilized by an N-heterocyclic dicarbene ligand.68 The structure of the N2 complex is provided in Scheme 9. The proposed C−N−C ligand, which has also been synthesized experimentally,69 is similar to the N-heterocyclic dicarbene ligand in Scheme 9, with the two aryl groups (−2,6-iPr2C6H3) replaced by two methyl groups (−CH3). The reaction scheme for Fe(C−N−C) carbene catalyst generation and activation is presented in Scheme 10. Experiments have shown that the Fe(N2)2 complex can be



CONCLUDING REMARKS Iron olefin metathesis has generated interest for a long time. To the best of our knowledge, despite heavy research efforts, there is no successful example of the utilization of an Fe complex for an olefin metathesis reaction. The difficulty lies in the fact that most Fe complexes favor the cyclopropanation reaction rather than olefin metathesis, which includes cycloaddition and cycloreversion. For the Fe catalyst to work for metathesis, the cyclopropanation reaction has to be prevented. G

DOI: 10.1021/acs.organomet.8b00583 Organometallics XXXX, XXX, XXX−XXX

Article

Organometallics

Table 1. Computed Electronic Energies (in kcal/mol) of Three Fe Compounds in Various Spin States; The Experimentally Determined Most Stable Spin State for Each Compound Is Used as the Reference State for Relative Energies

OPBE M06 MN15 MN15-L M06-L revM06-L SOGGA11-X experiment

doublet

quartet

sextet

singlet

triplet

quintet

singlet

triplet

quintet

4.0 25.0 5.7 33.1 14.0 12.5 25.0 >0

6.3 15.7 −4.0 20.2 12.7 10.8 13.7 >0

0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0

0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0

12.2 10.9 8.0 13.1 17.9 18.8 7.3 >0

3.2 −10.4 0.8 −14.2 7.1 8.2 −13.8 >0

0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0

18.9 26.1 3.7 6.6 23.0 23.9 3.6 >0

2.3 −18.5 −7.4 −27.9 −3.9 −0.8 −20.9 >0

metathesis products/intermediates. Accordingly, the combination of the α-H elimination and insertion reactions should have only a minor effect compared with the metathesis reaction. Finally, we have proposed a route to synthesize the Fe(C− N−C) carbene complex utilizing the Fe(N2)2 complex stabilized by an N-heterocyclic dicarbene ligand. The synthesis of the Fe(C−N−C) carbene complex would involve tautomerization of a coordinated terminal alkyne complex followed by the dissociation of the N2 molecule and association of the alkene molecule to produce the active species for the metathesis reaction.

In this study, we performed density functional theory calculations to study the catalytic effect of coordinated pincer ligands in the Fe organometallic catalyst for olefin metathesis. We used the OPBE functional along with the def2-TZVP basis set for most of our calculations. The OPBE density functional was chosen because of its ability to predict the correct spin state for various Fe organometallic complexes. The results suggest that the olefin metathesis reaction will be energetically more favorable than the cyclopropanation side reaction when the Fe active center is less positively charged. On the basis of this result, we propose the Fe(C−N−C) carbene species as a potential olefin metathesis catalyst, where C−N−C = 2,6bis(methylimidazol-2-ylidene)pyridine. Energy profiles were computed using the OPBE functional for propylene metathesis. The calculations showed that the overall free energy of activation for propylene metathesis is 11.6 kcal/mol. Both the cycloaddition and cycloreversion reactions have lower barriers than the cyclopropanation reaction, suggesting that the later is disfavored. Free energies of activation were also calculated using other density functionals, in particular M06-L, revM06-L, MN15-L, M06, and MN15. The results show that the OPBE functional predicts barriers that are about 8 kcal/mol lower than other functionals. Among the six applied functionals, only the M06 functional favors the cyclopropanation side reaction over propylene metathesis. We verified the stability of the proposed Fe carbene catalyst by considering two possible deformation pathways. The first reaction pathway includes the transformation of the carbene group into a coordinated alkene molecule, which requires βhydride elimination of the carbene group. Calculations showed that the free energy of activation for β-hydride elimination is ∼5 kcal/mol higher than that of the metathesis reaction, indicating that the formation of a coordinated alkene molecule is unlikely. The second reaction pathway included the formation of a hydride carbyne intermediate through αhydride elimination followed by insertion of a propylene molecule into the newly formed hydride. The α-hydride elimination is energetically more favorable than the metathesis reactions, but the subsequent insertion would produce a side product that is thermodynamically less stable than any of the



COMPUTATIONAL DETAILS

All of the calculations were carried out using unrestricted Kohn− Sham density functional theory as implemented in the Gaussian software. Gaussian 1672 was employed for all of the calculations except those using the revM06-L functional, which were performed using version 6.8 of MN-GFM (Minnesota Gaussian Functional Module), which is a locally modified version of Gaussian 0973 for incorporating new density functionals into the Gaussian code. All species were studied in vacuum with the def2-TZVP basis set.74,75 Except where stated otherwise, all of the results in this study are based on calculations using the OPBE functional. The OPBE functional combines Handy’s optimized exchange functional (OPTX)53,54 with the PBE correlation functional,55,56 and it was chosen because it correctly predicts the spin states of various Fe complexes in some previous tests.76−78 To evaluate the performance of some more recent density functionals, we calculated the spin-state energetics of three Fe compounds using seven different functionals, namely, OPBE, M06,60 M06-L,57 revM06-L,58 MN15,61 MN15-L,59 and SOGGA11-X.79 The compounds and their experimentally measured ground spin states are given in Table 1. Geometry optimizations were performed for all of the compounds for different spin states using different functionals. The computed energy for each optimized structure is summarized in Table 1. For this test, OPBE was the only functional that consistently predicted the correct ground spin state for Fe complexes. For catalyst screening and energy profile calculations, the geometries of reagents, transition structures, and reactive intermediates were optimized using the OPBE functional. The singlet and triplet spin states were considered in all of the calculations, and we also included quintet spin states when calculating energy profiles for H

DOI: 10.1021/acs.organomet.8b00583 Organometallics XXXX, XXX, XXX−XXX

Article

Organometallics propylene metathesis over the Fe(C−N−C) carbene catalyst. For each case, various starting geometries were considered. Spin contamination was small (the computed expectation value of S2 had less than a 5% difference from the correct value of S(S + 1), where S is half the number of unpaired electrons) in all geometries with the most stable spin state. Transition structures were optimized using the eigenmode-following method (by using the Gaussian keyword TS). Frequency analysis was performed on all of the optimized geometries; for each equilibrium geometry, we verified that all of the frequencies of stable species were real, and for each transition structure, we verified that there was only one imaginary frequency. For each optimized geometry, we computed the Gibbs free energy (G) at 298.15 K at a standard pressure of 1 atm. We used the FREQ program80 to generate scaling factors to compute the frequencies to obtain G. Scaling factors of 0.998, 0.981, 0.977, 0.972, 0.975, and 0.980 were used for the OPBE, M06, M06-L, revM06-L, MN15, and MN15-L functional-based calculations, respectively. Real frequencies below 100 cm−1 were raised to 100 cm−1 to simulate low-frequency anharmonic effects.81 The calculated free energies and optimized coordinates are given in the SI. In the article proper, we discuss both free energies of reaction (ΔG) and free energies of activation (ΔG⧧), but Born−Oppenheimer energies are given in the SI.



(7) Sulake, R. S.; Chen, C. P. Total Synthesis of (+)-Antroquinonol and (+)-Antroquinonol D. Org. Lett. 2015, 17, 1138−1141. (8) Seetharamsingh, B.; Rajamohanan, P. R.; Reddy, D. S. Total Synthesis and Structural Revision of Mycalol, an Anticancer Natural Product from the Marine Source. Org. Lett. 2015, 17, 1652−1655. (9) Elacqua, E.; Weck, M. Fabrication of Supramolecular Semiconductor Block Copolymers by Ring-Opening Metathesis Polymerization. Chem. - Eur. J. 2015, 21, 7151−7158. (10) Martinez, H.; Ren, N.; Matta, M. E.; Hillmyer, M. A. RingOpening Metathesis Polymerization of 8-Membered Cyclic Olefins. Polym. Chem. 2014, 5, 3507−3532. (11) Kim, Y. W.; Grossmann, T. N.; Verdine, G. L. Synthesis of AllHydrocarbon Stapled α-Helical Peptides by Ring-Closing Olefin Metathesis. Nat. Protoc. 2011, 6, 761−771. (12) Walensky, L. D.; Bird, G. H. Hydrocarbon-Stapled Peptides: Principles, Practice, and Progress. J. Med. Chem. 2014, 57, 6275− 6288. (13) Chhabra, S.; Belgi, A.; Bartels, P.; van Lierop, B. J.; Robinson, S. D.; Kompella, S. N.; Hung, A.; Callaghan, B. P.; Adams, D. J.; Robinson, A. J.; Norton, R. S. Dicarba Analogues of α-Conotoxin RgIA. Structure, Stability, and Activity at Potential Pain Targets. J. Med. Chem. 2014, 57, 9933−9944. (14) Lin, Y. Y. A.; Chalker, J. M.; Davis, B. G. Olefin Metathesis for Site-Selective Protein Modification. ChemBioChem 2009, 10, 959− 969. (15) Chauvin, Y. Olefin Metathesis: The Early Days (Nobel Lecture). Angew. Chem., Int. Ed. 2006, 45, 3740−3747. (16) Schrock, R. R. Recent Advances in High Oxidation State Mo and W Imido Alkylidene Chemistry. Chem. Rev. 2009, 109, 3211− 3226. (17) Vougioukalakis, G. C.; Grubbs, R. H. Ruthenium-Based Heterocyclic Carbene-Coordinated Olefin Metathesis Catalysts. Chem. Rev. 2010, 110, 1746−1787. (18) Schrock, R. R. Multiple Metal−Carbon Bonds for Catalytic Metathesis Reactions (Nobel Lecture). Angew. Chem., Int. Ed. 2006, 45, 3748−3759. (19) Grubbs, R. H. Olefin-Metathesis Catalysts for the Preparation of Molecules and Materials (Nobel Lecture). Angew. Chem., Int. Ed. 2006, 45, 3760−3765. (20) Hérisson, P. J.-L.; Chauvin, Y. Catalyse de Transformation des Oléfines par les Complexes Du Tungstène. II. Télomérisation des Oléfines Cycliques en Présence d’Oléfines Acycliques. Makromol. Chem. 1971, 141, 161−176. (21) Murdzek, J. S.; Schrock, R. R. Well-Characterized Olefin Metathesis Catalysts that Contain Molybdenum. Organometallics 1987, 6, 1373−1374. (22) Schrock, R. R.; Murdzek, J. S.; Bazan, G. C.; Robbins, J.; Dimare, M.; O’Regan, M. Synthesis of Molybdenum Imido Alkylidene Complexes and Some Reactions Involving Acyclic Olefins. J. Am. Chem. Soc. 1990, 112, 3875−3886. (23) Schrock, R. R.; Hoveyda, A. H. Molybdenum and Tungsten Imido Alkylidene Complexes as Efficient Olefin-Metathesis Catalysts. Angew. Chem., Int. Ed. 2003, 42, 4592−4633. (24) Toreki, R.; Schrock, R. R. A Well-Defined Rhenium(VII) Olefin Metathesis Catalyst. J. Am. Chem. Soc. 1990, 112, 2448−2449. (25) Nguyen, S. T.; Johnson, L. K.; Grubbs, R. H.; Ziller, J. W. RingOpening Metathesis Polymerization (ROMP) of Norbornene by a Group VIII Carbene Complex in Protic Media. J. Am. Chem. Soc. 1992, 114, 3974−3975. (26) Scholl, M.; Ding, S.; Lee, C. W.; Grubbs, R. H. Synthesis and Activity of a New Generation of Ruthenium-Based Olefin Metathesis Catalysts Coordinated with 1,3-Dimesityl-4,5-dihydroimidazol-2ylidene Ligands. Org. Lett. 1999, 1, 953−956. (27) de Brito Sá, É .; Rodríguez-Santiago, L.; Sodupe, M.; SolansMonfort, X. Toward Olefin Metathesis with Iron Carbene Complexes: Benefits of Tridentate σ-Donating Ligands. Organometallics 2016, 35, 3914−3923. (28) Luque-Urrutia, J. A.; Gimferrer, M.; Casals-Cruañas, È .; Poater, A. In Silico Switch from Second- to First-Row Transition Metals in

ASSOCIATED CONTENT

* Supporting Information S

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.organomet.8b00583. Molecular orbitals of singlet Fe(C−N−C) metallocyclobutane and Cartesian coordinates and corresponding electronic energies (E) and Gibbs free energies (G) for all of the optimized Fe and Ru complexes (PDF) Cartesian coordinate files (ZIP)



AUTHOR INFORMATION

Corresponding Authors

*E-mail: [email protected]. *E-mail: [email protected]. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS The authors are grateful to Laura Gagliardi, Christopher J. Cramer, Andreas A. Danopoulos, and John E. Ellis for helpful discussions. This work was supported as part of the Inorganometallic Catalysis Design Center, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Basic Energy Sciences, under Award DESC0012702.



REFERENCES

(1) Mol, J. C. Industrial Applications of Olefin Metathesis. J. Mol. Catal. A: Chem. 2004, 213, 39−45. (2) Hoveyda, A. H.; Zhugralin, A. R. The Remarkable MetalCatalysed Olefin Metathesis Reaction. Nature 2007, 450, 243−251. (3) Banks, R. L.; Bailey, G. C. Olefin Disproportionation. A New Catalytic Process. Ind. Eng. Chem. Prod. Res. Dev. 1964, 3, 170−173. (4) Chauvin, Y.; Commereuc, D. Chemical Counting and Characterization of the Active Sites in the Rhenium Oxide/Alumina Metathesis Catalyst. J. Chem. Soc., Chem. Commun. 1992, 462−464. (5) Lwin, S.; Wachs, I. E. Olefin Metathesis by Supported Metal Oxide Catalysts. ACS Catal. 2014, 4, 2505−2520. (6) Hog, D. T.; Huber, F. M. E.; Mayer, P.; Trauner, D. The Total Synthesis of (−)-Nitidasin. Angew. Chem., Int. Ed. 2014, 53, 8513− 8517. I

DOI: 10.1021/acs.organomet.8b00583 Organometallics XXXX, XXX, XXX−XXX

Article

Organometallics Olefin Metathesis: From Ru to Fe and from Rh to Co. Catalysts 2017, 7, 389−389. (29) Brookhart, M.; Tucker, J. R.; Husk, G. R. Synthesis, Spectral Characterization, and Alkylidene Transfer Reactions of Electrophilic Iron Carbene Complexes Cp(CO)(L)FeCHR+, L = CO, P(C6H5)3, R = CH3CH2CH3CH(CH3)2. J. Am. Chem. Soc. 1983, 105, 258−264. (30) Hamaker, C. G.; Mirafzal, G. A.; Woo, L. K. Catalytic Cyclopropanation with Iron(II) Complexes. Organometallics 2001, 20, 5171−5176. (31) Li, Y.; Huang, J.-S.; Zhou, Z.-Y.; Che, C.-M.; You, X.-Z. Remarkably Stable Iron Porphyrins Bearing Nonheteroatom-Stabilized Carbene or (Alkoxycarbonyl)carbenes: Isolation, X-ray Crystal Structures, and Carbon Atom Transfer Reactions with Hydrocarbons. J. Am. Chem. Soc. 2002, 124, 13185−13193. (32) Du, G.; Andrioletti, B.; Rose, E.; Woo, L. K. Asymmetric Cyclopropanation of Styrene Catalyzed by Chiral Macrocyclic Iron(II) Complexes. Organometallics 2002, 21, 4490−4495. (33) Edulji, S. K.; Nguyen, S. T. Catalytic Olefin Cyclopropanation Using μ-Oxo−bis[(salen)iron(III)] Complexes. Organometallics 2003, 22, 3374−3381. (34) Chen, Y.; Zhang, X. P. Asymmetric Cyclopropanation of Styrenes Catalyzed by Metal Complexes of D2-Symmetrical Chiral Porphyrin: Superiority of Cobalt over Iron. J. Org. Chem. 2007, 72, 5931−5934. (35) Russell, S. K.; Hoyt, J. M.; Bart, S. C.; Milsmann, C.; Stieber, S. C. E.; Semproni, S. P.; DeBeer, S.; Chirik, P. J. Synthesis, Electronic Structure and Reactivity of Bis(imino)pyridine Iron Carbene Complexes: Evidence for a Carbene Radical. Chem. Sci. 2014, 5, 1168−1174. (36) Dzik, W. I.; Zhang, X. P.; De Bruin, B. Redox Noninnocence of Carbene Ligands: Carbene Radicals in (Catalytic) C−C Bond Formation. Inorg. Chem. 2011, 50, 9896−9903. (37) Vasiliu, M.; Arduengo, A. J.; Dixon, D. A. Role of Electronegative Substituents on the Bond Energies in the Grubbs Metathesis Catalysts for M = Fe, Ru, Os. J. Phys. Chem. C 2014, 118, 13563−13577. (38) Eisenstein, O.; Hoffmann, R.; Rossi, A. R. Some Geometrical and Electronic Features of the Intermediate Stages of Olefin Metathesis. J. Am. Chem. Soc. 1981, 103, 5582−5584. (39) Poater, A.; Vummaleti, S. V. C.; Pump, E.; Cavallo, L. Comparing Ru and Fe-catalyzed olefin metathesis. Dalton Trans. 2014, 43, 11216−11220. (40) Poater, A.; Pump, E.; Vummaleti, S. V. C.; Cavallo, L. The Activation Mechanism of Fe-based Olefin Metathesis Catalysts. Chem. Phys. Lett. 2014, 610−611, 29−32. (41) Poater, A. Environmental Friendly Fe Substitutive of Ru in Water Oxidation Catalysis. Catal. Commun. 2014, 44, 2−5. (42) Chabanas, M.; Baudouin, A.; Coperet, C.; Basset, J. M. A Highly Active Well-Defined Rhenium Heterogeneous Catalyst for Olefin Metathesis Prepared via Surface Organometallic Chemistry. J. Am. Chem. Soc. 2001, 123, 2062−2063. (43) Solans-Monfort, X.; Clot, E.; Coperet, C.; Eisenstein, O. d0 ReBased Olefin Metathesis Catalysts, Re(CR)(CHR)(X)(Y): The Key Role of X and Y Ligands for Efficient Active Sites. J. Am. Chem. Soc. 2005, 127, 14015−14025. (44) Nørskov, J. K.; Bligaard, T.; Rossmeisl, J.; Christensen, C. H. Towards the Computational Design of Solid Catalysts. Nat. Chem. 2009, 1, 37−46. (45) Mauksch, M.; Tsogoeva, B. Iron-Catalyzed Olefin Metathesis with Low-Valen Iron Alkylidenes. Chem. - Eur. J. 2017, 23, 10264− 10269. (46) de Brito Sá, É .; Rodríguez-Santiago, L.; Sodupe, M.; SolansMonfort, X. Influence of Ligands and Oxidation State on the Reactivity of Pentacoordinated Iron Carbenes with Olefins: Metathesis versus Cyclopropanation. Organometallics 2018, 37, 1229− 1241. (47) Yang, B.; Gagliardi, L.; Truhlar, D. G. Transition States of SpinForbidden Reactions. Phys. Chem. Chem. Phys. 2018, 20, 4129−4136.

(48) Marenich, A. V.; Jerome, C.; Cramer, C. J.; Truhlar, D. G. Charge Model 5: An Extension of Hirshfeld Population Analysis for the Accurate Description of Molecular Interactions in Gaseous and Condensed Phases. J. Chem. Theory Comput. 2012, 8, 527−541. (49) Parr, R. G.; Szentpály, L. V.; Liu, S. Electrophilicity Index. J. Am. Chem. Soc. 1999, 121, 1922−1924. (50) Domingo, L. R.; Aurell, M. J.; Pérez, P.; Contreras, R. Quantitative Characterization of the Global Electrophilicity Power of Common Diene/Dienophile Pairs in Diels−Alder Reactions. Tetrahedron 2002, 58, 4417. (51) Domingo, L. R.; Chamorro, E.; Pérez, P. An Understanding of the Electrophilic/Nucleophilic Behavior of Electro-Deficient 2,3Disubstituted 1,3-Butadienes in Polar Diels−Alder Reactions. A Density Functional Theory Study. J. Phys. Chem. A 2008, 112, 4046− 4053. (52) Chattaraj, P. K.; Roy, D. R. Update 1 of: Electrophilicity Index. Chem. Rev. 2007, 107, PR46−PR74. (53) Handy, N. C.; Cohen, A. J. Left−Right Correlation Energy. Mol. Phys. 2001, 99, 403−412. (54) Hoe, W.-M.; Cohen, A.; Handy, N. C. Assessment of a New Local Exchange Functional OPTX. Chem. Phys. Lett. 2001, 341, 319− 328. (55) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865−3868. (56) Perdew, J. P.; Burke, K.; Ernzerhof, M. Errata: Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1997, 78, 1396−1396. (57) Zhao, Y.; Truhlar, D. G. A New Local Density Functional for Main-Group Thermochemistry, Transition Metal Bonding, Thermochemical Kinetics, and Noncovalent Interactions. J. Chem. Phys. 2006, 125, 194101−194101. (58) Wang, Y.; Jin, X.; Yu, H. S.; Truhlar, D. G.; He, X. Revised M06-L Functional for Improved Accuracy on Chemical Reaction Barrier Heights, Noncovalent Interactions, and Solid-State Physics. Proc. Natl. Acad. Sci. U. S. A. 2017, 114, 8487−8492. (59) Yu, H. S.; He, X.; Truhlar, D. G. MN15-L: A New Local Exchange-Correlation Functional for Kohn−Sham Density Functional Theory with Broad Accuracy for Atoms, Molecules, and Solids. J. Chem. Theory Comput. 2016, 12, 1280−1293. (60) Zhao, Y.; Truhlar, D. G. The M06 Suite of Density Functionals for Main Group Thermochemistry, Thermochemical Kinetics, Noncovalent Interactions, Excited States, and Transition Elements: Two New Functionals and Systematic Testing of Four M06-Class Functionals and 12 Other Functionals. Theor. Chem. Acc. 2008, 120, 215−241. (61) Yu, H. S.; He, X.; Li, S. L.; Truhlar, D. G. MN15: A Kohn− Sham Global-Hybrid Exchange-Correlation Density Functional with Broad Accuracy for Multi-Reference and Single-Reference Systems and Noncovalent Interactions. Chem. Sci. 2016, 7, 5032−5051. (62) Cramer, C. J.; Truhlar, D. G. Density Functional Theory for Transition Metals and Transition Metal Chemistry. Phys. Chem. Chem. Phys. 2009, 11, 10757−10816. (63) Minenkov, Y.; Singstad, Å.; Occhipinti, G.; Jensen, V. R. The Accuracy of DFT-Optimized Geometries of Functional Transition Metal Compounds: A Validation Study of Catalysts for Olefin Metathesis and Other Reactions in the Homogeneous Phase. Dalton Trans. 2012, 41, 5526−5541. (64) Zhao, Y.; Truhlar, D. G. Benchmark Energetic Data in a Model System for Grubbs II Metathesis Catalysis and Their Use for the Development, Assessment, and Validation of Electronic Structure Methods. J. Chem. Theory Comput. 2009, 5, 324−333. (65) Zhao, Y.; Truhlar, D. G. Attractive Noncovalent Interactions in the Mechanism of Grubbs Second-Generation Ru Catalysts for Olefin Metathesis. Org. Lett. 2007, 9, 1967−1970. (66) Yang, H.-C.; Huang, Y.-C.; Lan, Y.-K.; Luh, T.-Y.; Zhao, Y.; Truhlar, D. G. Carbene Rotamer Switching Explains the Reverse Trans Effect in Forming the Grubbs Second-Generation Olefin Metathesis Catalyst. Organometallics 2011, 30, 4196−4200. J

DOI: 10.1021/acs.organomet.8b00583 Organometallics XXXX, XXX, XXX−XXX

Article

Organometallics

Structure Model Chemistries. Comput. Phys. Commun. 2017, 210, 132−138. (81) Pratt, L. M.; Truhlar, D. G.; Cramer, C. J.; Kass, S. R.; Thompson, J. D.; Xidos, J. D. Aggregation of Alkyllithiums in Tetrahydrofuran. J. Org. Chem. 2007, 72, 2962−2966.

(67) de Frémont, P.; Marion, N.; Nolan, S. P. Carbenes: Synthesis, Properties, and Organometallic Chemistry. Coord. Chem. Rev. 2009, 253, 862−892. (68) Danopoulos, A. A.; Wright, J. A.; Motherwell, W. B. Molecular N2 Complexes of Iron Stabilised by N-Heterocyclic ‘Pincer’ Dicarbene Ligands. Chem. Commun. 2005, 784−786. (69) Peris, E.; Loch, J. A.; Mata, J.; Crabtree, R. H. A Pd Complex of a Tridentate Pincer CNC Bis-Carbene Ligand as a Robust Homogenous Heck Catalyst. Chem. Commun. 2001, 201−202. (70) Dragutan, V.; Dragutan, I. Ruthenium Vinylidene Complexes. Platinum Met. Rev. 2004, 48, 148−153. (71) Stegmann, R.; Frenking, G. Mechanism of the AcetyleneVinylidene Rearrangement in the Coordination Sphere of a Transition Metal. Organometallics 1998, 17, 2089−2095. (72) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Petersson, G. A.; Nakatsuji, H.; Li, X.; Caricato, M.; Marenich, A. V.; Bloino, J.; Janesko, B. G.; Gomperts, R.; Mennucci, B.; Hratchian, H. P.; Ortiz, J. V.; Izmaylov, A. F.; Sonnenberg, J. L.; Williams-Young, D.; Ding, F.; Lipparini, F.; Egidi, F.; Goings, J.; Peng, B.; Petrone, A.; Henderson, T.; Ranasinghe, D.; Zakrzewski, V. G.; Gao, J.; Rega, N.; Zheng, G.; Liang, W.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Throssell, K.; Montgomery, J. A., Jr.; Peralta, J. E.; Ogliaro, F.; Bearpark, M. J.; Heyd, J. J.; Brothers, E. N.; Kudin, K. N.; Staroverov, V. N.; Keith, T. A.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A. P.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Millam, J. M.; Klene, M.; Adamo, C.; Cammi, R.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Farkas, O.; Foresman, J. B.; Fox, D. J. Gaussian 16, revision A.03; Gaussian, Inc.: Wallingford, CT, 2016. (73) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; Nakatsuji, H.; Caricato, M.; Li, X.; Hratchian, H. P.; Izmaylov, A. F.; Bloino, J.; Zheng, G.; Sonnenberg, J. L.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Montgomery, J. A., Jr.; Peralta, J. E.; Ogliaro, F.; Bearpark, M.; Heyd, J. J.; Brothers, E.; Kudin, K. N.; Staroverov, V. N.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Rega, N.; Millam, J. M.; Klene, M.; Knox, J. E.; Cross, J. B.; Bakken, V.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Zakrzewski, V. G.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Dapprich, S.; Daniels, A. D.; Farkas, Ö .; Foresman, J. B.; Ortiz, J. V.; Cioslowski, J.; Fox, D. J. Gaussian 09, revision C.01; Gaussian, Inc.: Wallingford, CT, 2009. (74) Weigend, F. Accurate Coulomb-Fitting Basis Sets for H to Rn. Phys. Chem. Chem. Phys. 2006, 8, 1057−1065. (75) Weigend, F.; Ahlrichs, R. Balanced Basis Sets of Split Valence, Triple Zeta Valence and Quadruple Zeta Valence Quality for H to Rn: Design and Assessment of Accuracy. Phys. Chem. Chem. Phys. 2005, 7, 3297−3305. (76) Swart, M.; Groenhof, A. R.; Ehlers, A. W.; Lammertsma, K. Validation of Exchange-Correlation Functionals for Spin States of Iron Complexes. J. Phys. Chem. A 2004, 108, 5479−5483. (77) Swart, M. Accurate Spin-State Energies for Iron Complexes. J. Chem. Theory Comput. 2008, 4, 2057−2066. (78) Rousseau, L.; Brémond, E.; Lefèvre, G. Assessment of the Ground Spin State of Iron(I) Complexes: Insights from DFT PredictiveModels. New J. Chem. 2018, 42, 7612−7616. (79) Peverati, R.; Truhlar, D. G. Communication: A gGlobal Hybrid Generalized Gradient Approximation to the Exchange-Correlation Functional that Satisfies the Second-Order Density-Gradient Constraint and has Broad Applicability in Chemistry. J. Chem. Phys. 2011, 135, 191102−191102. (80) Yu, H. S.; Fiedler, L. J.; Alecu, I. M.; Truhlar, D. G. Computational Thermochemistry: Automated Generation of Scale Factors for Vibrational Frequencies Calculated by Electronic K

DOI: 10.1021/acs.organomet.8b00583 Organometallics XXXX, XXX, XXX−XXX