Understanding Structural and Dynamic Properties of Well-Defined Rhenium-Based Olefin Metathesis Catalysts, Re(⋮CR)(CHR)(X)(Y), from DFT and QM/MM Calculations

Effect of Lewis Acids on the Catalyst Activity for Alkene Metathesis, Z-/E- Selectivity and Stability of Tungsten Oxo Alkylidenes

Lewis acids increase the catalytic activity of classical heterogeneous catalysts and molecular d0 tungsten oxo alkylidenes in a variety of olefin metathesis processes. The formation of labile adducts between the metal complex and the Lewis acid has been observed experimentally and suggested to be involved in the catalyst activity increase. In this contribution, DFT (M06) calculations have been performed to determine the role of Lewis acids on catalyst activity, Z-/E- selectivity and stability by comparing three W(E)(CHR)(2,5-dimethylpyrrolide)(O-2,6-dimesithylphenoxide) (E = oxo, imido or oxo-Lewis acid adduct) alkylidenes. Results show that the formation of the alkylidene—Lewis acid adducts influences the reactivity of tungsten oxo alkylidenes due to both steric and electronic effects. The addition of the Lewis acid on the E group increases its bulkiness and this decreases catalyst Z-selectivity. Moreover, the interaction between the oxo ligand and the Lewis acid decreases the donating ability of the former toward the metal. This is important when the oxo group has either a ligand in trans or in the same plane that is competing for the same metal d orbitals. Therefore, the weakening of oxo donating ability facilitates the cycloaddition and cycloreversion steps and it stabilizes the productive trigonal bipyramid metallacyclobutane isomer. The two factors increase the catalytic activity of the complex. The electron donating tuneability by the coordination of the Lewis acid also applies to catalyst deactivation and particularly the key β-hydride elimination step. In this process, the transition states show a ligand in pseudo trans to the oxo. Therefore, the presence of the Lewis acid decreases the Gibbs energy barrier significantly. Overall, the optimization of the E group donating ability in each step of the reaction makes tungsten oxo alkylidenes more reactive and this applies both for the catalytic activity and catalyst deactivation.

Olefin metathesis: what have we learned about homogeneous and heterogeneous catalysts from surface organometallic chemistry?

Since its early days, olefin metathesis has been in the focus of scientific discussions and technology development. While heterogeneous olefin metathesis catalysts based on supported group 6 metal oxides have been used for decades in the petrochemical industry, detailed mechanistic studies and the development of molecular organometallic chemistry have led to the development of robust and widely used homogeneous catalysts based on well-defined alkylidenes that have found applications for the synthesis of fine and bulk chemicals and are also used in the polymer industry. The development of the chemistry of high-oxidation group 5-7 alkylidenes and the use of surface organometallic chemistry (SOMC) principles unlocked the preparation of so-called well-defined supported olefin metathesis catalysts. The high activity and stability (often superior to their molecular analogues) and molecular-level characterisation of these systems, that were first reported in 2001, opened the possibility for the first direct structure-activity relationships for supported metathesis catalysts. This review describes first the history of SOMC in the field of olefin metathesis, and then focuses on what has happened since 2007, the date of our last comprehensive reviews in this field.

Embedded Mean-Field Theory for Solution-Phase Transition-Metal Polyolefin Catalysis

Decreasing the wall-clock time of quantum mechanics/molecular mechanics (QM/MM) calculations without sacrificing accuracy is a crucial prerequisite for widespread simulation of solution-phase dynamical processes. In this work, we demonstrate the use of embedded mean-field theory (EMFT) as the QM engine in QM/MM molecular dynamics (MD) simulations to examine polyolefin catalysts in solution. We show that employing EMFT in this mode preserves the accuracy of hybrid-functional DFT in the QM region, while providing up to 20-fold reductions in the cost per SCF cycle, thereby increasing the accessible simulation time-scales. We find that EMFT reproduces DFT-computed binding energies and optimized bond lengths to within chemical accuracy, as well as consistently ranking conformer stability. Furthermore, solution-phase EMFT/MM simulations provide insight into the interaction strength of strongly coordinating and bulky counterions.

Understanding surface site structures and properties by first principles calculations: an experimental point of view!

Computational Chemistry is key for the molecular-level understanding of active sites in heterogeneous catalysis paving the way to the rational design and development.

The Nature of Secondary Interactions at Electrophilic Metal Sites of Molecular and Silica-Supported Organolutetium Complexes from Solid-State NMR Spectroscopy

Lu[CH(SiMe3)2]3 reacts with [SiO2-700] to give [(≡SiO)Lu[CH(SiMe3)2]2] and CH2(SiMe3)2. [(≡SiO)Lu[CH(SiMe3)2]2] is characterized by solid-state NMR and EXAFS spectroscopy, which show that secondary Lu···C and Lu···O interactions, involving a γ-CH3 and a siloxane bridge, are present. From X-ray crystallographic analysis, the molecular analogues Lu[CH(SiMe3)2]3-x[O-2,6-tBu-C6H3]x (x = 0-2) also have secondary Lu···C interactions. The (1)H NMR spectrum of Lu[CH(SiMe3)2]3 shows that the -SiMe3 groups are equivalent to -125 °C and inequivalent below that temperature, ΔG(⧧)(Tc = 148 K) = 7.1 kcal mol(-1). Both -SiMe3 groups in Lu[CH(SiMe3)2]3 have (1)JCH = 117 ± 1 Hz at -140 °C. The solid-state (13)C CPMAS NMR spectrum at 20 °C shows three chemically inequivalent resonances in the area ratio of 4:1:1 (12:3:3); the J-resolved spectra for each resonance give (1)JCH = 117 ± 2 Hz. The (29)Si CPMAS NMR spectrum shows two chemically inequivalent resonances with different values of chemical shift anisotropy. Similar observations are obtained for Lu[CH(SiMe3)2]3-x[O-2,6-tBu-C6H3]x (x = 1 and 2). The spectroscopic data point to short Lu···Cγ contacts corresponding to 3c-2e Lu···Cγ-Siβ interactions, which are supported by DFT calculations. Calculated natural bond orbital (NBO) charges show that Cγ carries a negative charge, while Lu, Hγ, and Siβ carry positive charges; as the number of O-based ligands increases so does the positive charge at Lu, which in turns shortens the Lu···Cγ distance. The change in NBO charges and the resulting changes in the spectroscopic and crystallographic properties show how ligands and surface-support sites rearrange to accommodate these changes, consistent with Pauling's electroneutrality concept.

Quantitatively analyzing metathesis catalyst activity and structural features in silica-supported tungsten imido-alkylidene complexes

A broad series of fully characterized, well-defined silica-supported W metathesis catalysts with the general formula [(≡SiO)W(═NAr)(═CHCMe2R)(X)] (Ar = 2,6-iPr2C6H3 (AriPr), 2,6-Cl2C6H3 (ArCl), 2-CF3C6H4 (ArCF3), and C6F5 (ArF5); X = OC(CF3)3 (OtBuF9), OCMe(CF3)2 (OtBuF6), OtBu, OSi(OtBu)3, 2,5-dimethylpyrrolyl (Me2Pyr) and R = Me or Ph) was prepared by grafting bis-X substituted complexes [W(NAr)(═CHCMe2R)(X)2] on silica partially dehydroxylated at 700 °C (SiO2-(700)), and their activity was evaluated with the goal to obtain detailed structure-activity relationships. Quantitative influence of the ligand set on the activity (turnover frequency, TOF) in self-metathesis of cis-4-nonene was investigated using multivariate linear regression analysis tools. The TOF of these catalysts (activity) can be well predicted from simple steric and electronic parameters of the parent protonated ligands; it is described by the mutual contribution of the NBO charge of the nitrogen or the IR intensity of the symmetric N-H stretch of the ArNH2, corresponding to the imido ligand, together with the Sterimol B5 and pKa of HX, representing the X ligand. This quantitative and predictive structure-activity relationship analysis of well-defined heterogeneous catalysts shows that high activity is associated with the combination of X and NAr ligands of opposite electronic character and paves the way toward rational development of metathesis catalysts.

DFT study on the reaction mechanism of the ring closing enyne metathesis (RCEYM) catalyzed by molybdenum alkylidene complexes

DFT calculations show that the RCEYM reaction catalyzed by Mo-based catalysts proceeds preferentially through an yne-then-ene mechanism and that the endo-/exo- selectivity mainly depends on sterics.

Stabilizing reactive intermediates through site isolation

This overview describes the reaction of organometallics with oxide surfaces and the formation of highly reactive species. In the case of silica, the surface can be seen as a large siloxy ligand, which helps to stabilize reactive intermediates through site isolations. This is translated into very highly reactive and stable well-defined alkene metathesis catalysts as well as the formation of hydrides species, which display unusual reactivities toward alkanes (e.g., low-temperature hydrogenolysis and metathesis of alkanes). In the case of alumina, it allows the formation of highly reactive, but stable cationic species or masked carbenic species whose structures are unusual by comparison with molecular chemistry.

Beta-H transfer from the metallacyclobutane: a key step in the deactivation and byproduct formation for the well-defined silica-supported rhenium alkylidene alkene metathesis catalyst

The surface complex [([triple bond]SiO)Re([triple bond]CtBu)(=CHtBu)(CH2tBu)] (1) is a highly efficient propene metathesis catalyst with high initial activities and a good productivity. However, it undergoes a fast deactivation process with time on stream, which is first order in active sites and ethene. Noteworthy, 1-butene and pentenes, unexpected products in the metathesis of propene, are formed as primary products, in large amount relative to Re (>>1 equiv/Re), showing that their formation is not associated with the formation of inactive species. DFT calculations on molecular model systems show that byproduct formation and deactivation start by a beta-H transfer trans to the weak sigma-donor ligand (siloxy) at the metallacyclobutane intermediate having a square-based pyramid geometry. This key step has an energy barrier slightly higher than that calculated for olefin metathesis. After beta-H transfer, the most accessible pathway is the insertion of ethene in the Re-H bond. The resulting pentacoordinated trisperhydrocarbyl complex rearranges via either (1) alpha-H abstraction yielding the unexpected 1-butene byproduct and the regeneration of the catalyst or (2) beta-H abstraction leading to degrafting. These deactivation and byproduct formation pathways are in full agreement with the experimental data.

Dinitrogen Dissociation on an Isolated Surface Tantalum Atom

Both industrial and biochemical ammonia syntheses are thought to rely on the cooperation of multiple metals in breaking the strong triple bond of dinitrogen. Such multimetallic cooperation for dinitrogen cleavage is also the general rule for dinitrogen reductive cleavage with molecular systems and surfaces. We have observed cleavage of dinitrogen at 250 degrees C and atmospheric pressure by dihydrogen on isolated silica surface-supported tantalum(III) and tantalum(V) hydride centers [(identical with Si-O)2Ta(III)-H] and [(identical with Si-O)2Ta(V)H3], leading to the Ta(V) amido imido product [(identical with SiO)2Ta(=NH)(NH2)]: We assigned the product structure based on extensive characterization by infrared and solid-state nuclear magnetic resonance spectroscopy, isotopic labeling studies, and supporting data from x-ray absorption and theoretical simulations. Reaction intermediates revealed by in situ monitoring of the reaction with infrared spectroscopy support a mechanism highly distinct from those previously observed in enzymatic, organometallic, and heterogeneous N2 activating systems.

Understanding d0-Olefin Metathesis Catalysts: Which Metal, Which Ligands?

Density functional theory (DFT, B3PW91) calculations have been carried out on the reactivity of ethene with model systems M(NR)(=CHCH3)(X)(Y) for M = Mo or W, R = methyl or phenyl, X = CH2CH3, OCH3, or OSiH3, and Y = CH2CH3, OCH3, or OSiH3, which are representative of experimental olefin metathesis catalysts, and the results are compared to those previously obtained for Re(CCH3)(=CHCH3)(X)(Y). The general pathway comprises four steps: olefin coordination, [2+2] cycloaddition, cycloreversion, and olefin de-coordination. Two key factors have been found to control the detailed shape of the energy profiles: the energy of distortion of the tetrahedral catalyst and the stability of the metallacycle intermediate, which is controlled by the M-C bond strength. The efficiency has been evaluated by calculating the turnover frequency (TOF) based on the steady-state approximation, and the most striking feature is that the unsymmetrical catalysts (X not equal to Y) are systematically more efficient for all systems (Mo, W, and Re). Overall, the Re complexes have been found to be less efficient than the Mo and W catalysts, except when Re is unsymmetrically substituted: it is then calculated to be as efficient as the best Mo and W catalysts.

Formation and Decomposition of C3 Metallacycles from Ethylene and Methylene on MoAl Alloy Thin Films

The reaction between adsorbed ethylene and methylene species has been investigated on a molybdenum-aluminum alloy grown from Mo(CO)(6) on a planar alumina film formed on a Mo(100) single crystal in ultrahigh vacuum. Di-sigma-bonded ethylene reacts with carbene species, formed on the surface from methylene iodide, to form a C(3) metallacycle. This predominantly decomposes to yield propylene, while a smaller portion yields cross-metathesis products since (12)C(13)CH(4) is formed from reaction between (13)C(2)H(4) and (12)CH(2). This work demonstrates for the first time that the reaction proceeds in heterogeneous phase via a C(3) metallacycle as proposed in the Hérisson-Chauvin mechanism.

d0 Re-based olefin metathesis catalysts, Re([triple bond]CR)(=CHR)(X)(Y): the key role of X and Y ligands for efficient active sites

DFT(B3PW91) calculations show that the reaction pathways for ethylene metathesis with Re([triple bond]CMe)(=CHMe)(X)(Y) (X/Y = CH2CH3/CH2CH3; CH2CH3/OSiH3; OSiH3/CH2CH3; OCH3/OCH3, CH2CH3/OCH3, and OCF3/OCF3) occur in two steps: first, the pseudo-tetrahedral d0 Re complexes distort to a trigonal pyramid to open a coordination site for ethylene, which remains far from Re (early transition state for C-C bond formation). The energy barrier, determined by the energy required to distort the catalyst, is the lowest for unsymmetrical ligands (X not equal Y) when the apical site of the TBP is occupied by a good sigma-donor ligand (X) and the basal site by a poor sigma-donor (Y). Second, the formation of metallacyclobutanes (late transition state for C-C bond formation) has a low energy barrier for any type of ligands, decreasing for poor sigma-donor X and Y ligands, because they polarize the Re-C alkylidene bond as Re(+delta)=C(-delta), which favors the reaction with ethylene, itself polarized by the metal center in the reverse way. The metallacyclobutane is also a TBP, with apical alkylidyne and Y ligands, and it is stabilized by poor sigma-donor X and Y. The best catalyst will have the more shallow potential energy surface, and will thus be obtained for the unsymmetrical set of ligands with X = a good sigma-donor (alkyl) and Y = a poor sigma-donor (O-based ligand). This rationalizes the high efficiency of well-defined Re alkylidene supported on silica, compared to its homogeneous equivalent, Re([triple bond]CMe)(=CHMe)(OR)2.

DFT calculations of d0M(NR)(CHtBu)(X)(Y) (M = Mo, W; R = CPh3, 2,6-iPr–C6H3; X and Y = CH2tBu, OtBu, OSi(OtBu)3) olefin metathesis catalysts: structural, spectroscopic and electronic properties

DFT(B3PW91) calculations have been carried out to rationalise the structural, electronic and spectroscopic properties of Mo and W imido M(NR1)(CHR2)(X)(Y) olefin metathesis catalysts by using either simplified or actual ligands of the experimental complexes. The calculated structures, energetics (preference for the syn isomer and alkylidene rotational barrier for the syn/anti interconversion), and spectroscopic properties (NMR J(C-H) coupling constants) are in good agreement with available experimental data. Additionally, the alkylidene nu(C-H) stretching frequencies, not available experimentally, have been calculated. These quasi-tetrahedral complexes have a linear imido group and a C-H alkylidene agostic interaction, which stabilizes the syn isomer. Whether looking at M(NR1)(CHR2)(X)(Y), M = Mo, W, or the isolobal Re complexes, Re(CR1)(CHR2)(X)(Y), a linear correlation is obtained between both the alkylidene nu(C-H) stretching frequencies and J(C-H) coupling constants with the calculated alkylidene C-H bond lengths. These correlations show that the strength of the alpha-C-H agostic interaction increases from alkylidyne Re to imido group 6 complexes and from Mo to W. The NBO and AIM Bader analyses show firstly that the imido and alkylidyne groups are both triply bonded to the metal, but that the triply bonded imido ligand is a weaker electron donor than the alkylidyne, hence the stronger alpha-C-H agostic interaction for group 6 imido complexes. Secondly, one of the pi bonds of the triply bonded ligand is weakened at the transition state of the alkylidene rotation: while no lone pair is formed, the metal-ligand triple bond is polarized. This is more favourable for an imido than for an alkylidyne ligand, hence the lower alkylidene rotational barrier for the former complexes. Conversely, the aryl imido is even less of an electron donor than the alkyl imido group, which in turn strengthens the alpha-C-H agostic interaction and lowers the alkylidene rotational barrier even more.