Calculating accurate barriers for olefin insertion and related reactions

Predicting the propene stereoselectivity on transition metal catalysts: A daunting task for density functional theory

Thanks to recent developments in hardware and software, quantum chemical methods are increasingly used for interpreting the complex mechanisms underlying polymerization reaction by homogeneous catalysis. Unfortunately, the dimensions of even the smallest realistic models are too large to permit the use of state-of-the-art composite wave function methods. Under these circumstances, density functional theory still offers the best compromise between cost and accuracy. However, comprehensive benchmarks of different functionals are not yet available for this important research field. The main aim of the present paper is to fill this gap by performing an unbiased comparison of several density functionals and continuum solvent models for the stereo-control in the propylene polymerization on prototypical catalysts inducing different reaction mechanisms. While it was not possible to define a unique computational protocol providing the best results in all the situations, the B3PW91 functional in conjunction with D3 empirical dispersions and the solvent model density solvent model performs remarkably well for three out of the four investigated catalysts. Under such circumstances, it is recommended to compare the results delivered by different models when approaching additional classes of catalysts.

Triptycene as a scaffold in metallocene catalyzed olefin polymerization

A set of metallocene olefin polymerization catalysts bearing triptycene moieties in either position 4-5 (complexes Ty1-Ty5) or in position 5-6 (complexes Ty6-Ty8) of the basic dimethylsilyl-bridged bis(indenyl) system has been tested in propene polymerization and in ethene/1-hexene copolymerization. Comparison of the results with QSPR (quantitative structure-property relationship) predictions not parametrized for these exotic ligand variations demonstrates that trends can still be identified by extrapolation. Interestingly, Ty7, upon suitable activation, provides a highly isotactic polypropylene with an exceptional amount of 2,1 regio-errors (8%). The previously developed QSPR type models successfully predicted the low regioselectivity of this catalyst, despite the fact that the catalyst structure differs significantly from the benchmark set.

Sheet Models for Methylaluminoxane (MAO) Activators? A Theoretical Case Study involving rac-Me2Si(η5-C9H6)2Zr (SBIZr) Complexes

Activation of SBIZrMe2 or SBIZrMeCl and a sheet model for an active component of hydrolytic MAO, (MeAlO)16(Me3Al)6, (16,6) has been studied by DFT. Contact ion-pair formation occurs through the intermediacy of SBIZrMe(Cl) or SBIZrMe2 reacting with sheet 16,6 to furnish SBIZrMe-μ-X(MeAlO)16(Me3Al)6 (2, X=Me, Cl). Contact ion-pairs 2 would be in equilibrium with heterodinuclear catalyst precursors [SBIZrMe2AlMe2][(MeAlO)16(Me3Al)6X] (3 (X=Me, Cl) through reversible binding of Me3Al at higher Al : Zr ratios. Calculations show that formation of ion-pairs 3 from contact ion-pairs 2 is more favourable for the SBIZr compared with the parent Cp2Zr complexes. TD-DFT calculations were conducted on relevant SBIZr complexes to relate the results to earlier spectroscopic studies of catalyst activation using UV-Vis spectroscopy. Finally, propene insertion into ion-pairs 2, SBIZrMe-μ-MeB(C6F5)3 (6) and [SBIZrMe][B(C6F5)4] (7) was studied at M06-2X/TZVP level of theory. These studies suggest that contact ion-pairs 2 are significantly less reactive towards insertion than 6 or 7, in disagreement with experiment.

Ethylene/Polar Monomer Copolymerization by [N, P] Ti Complexes: Polar Copolymers with Ultrahigh-Molecular Weight

A series of novel titanium complexes (2a-2e) bearing [N, P] aniline-chlorodiphenylphosphine ligands (1a-1e) featuring CH3 and F substituents have been synthesized and characterized. Surprisingly, in the presence of polar additive, the complexes (2a-2e) all displayed high catalytic activities (up to 1.04 × 106 gPolymer (mol·Ti)-1·h-1 and produced copolymer with the ultrahigh molecular weight up to 1.37 × 106 g/mol. The catalytic activities are significantly enhanced by introducing electron-withdrawing group (F) into the aniline aromatic ring. Especially, the increase in activity based on different complexes followed the order of 2e > 2d > 2c > 2b > 2a. Simultaneously, density functional theory (DFT) calculations have been performed to probe the polymerization mechanism as well as the electronic and steric effects of various substituents on the catalyst backbone. DFT computation revealed that the polymerization behaviors could be adjusted by the electronic effect of ligand substituents; however, it has little to do with the steric hindrance of the substituents. Furthermore, theoretical calculation results keep well in accordance with experimental measurement results. The article provided an appealing design method that the employment of fluorine atom as electron-withdrawing to be studied is the promotive effect of transition-metal coordination polymerization.

Octahedral Zirconium Salan Catalysts for Olefin Polymerization: Substituent and Solvent Effects on Structure and Dynamics

Group 4 metal-Salan olefin polymerization catalysts typically have relatively low activity, being slowed down by a pre-equilibrium favoring a non-polymerization active resting state identified as a mer-mer isomer (MM); formation of the polymerization active fac-fac species (FF) requires isomerization. We now show that the chemistry is more subtle than previously realized. Salan variations bearing large, flat substituents can achieve very high activity, and we ascribe this to the stabilization of the FF isomer, which becomes lower in energy than MM. Detailed in situ NMR studies of a fast (o-anthracenyl) and a slow (o-tBu) Salan precursors, suitably activated, indicate that preferred isomers in solution are different: the fast catalyst prefers FF while the slow catalyst prefers a highly distorted MM geometry. Crystal structures of the activated o-anthracenyl substituted complex with a moderately (chlorobenzene) and, more importantly, a weakly coordinating solvent (toluene) in the first coordination sphere emphasize that the active FF isomer is preferred, at least for the benzyl species. Site epimerization (SE) barriers for the fast catalyst (ΔS > 0, dissociative) and the slow catalyst (ΔS < 0, associative) in toluene corroborate the solvent role. Diagnostic NMe 13C chemical shift differences allow unambiguous detection of FF or MM geometries for seven activated catalysts in different solvents, highlighting the role of solvent coordination strength and bulkiness of the ortho-substituent on the isomer equilibrium. For the first time, active polymeryl species of Zr-Salan catalysts were speciated. The slow catalyst is effectively trapped in the inactive MM state, as previously suggested. Direct observation of fast catalysts is hampered by their high reactivity, but the product of the first 1-hexene insertion maintains its FF geometry.

Arene Insertion with Pincer-Supported Molybdenum-Hydrides: Determination of Site Selectivity, Relative Rates, and Arene Complex Formation

The synthesis of phosphino(oxazoline)pyridine-supported molybdenum(0) cycloocta-1,5-diene complexes is described. Exposure of these complexes to dihydrogen in the presence of an arene resulted in insertion of the substrate into the molybdenum hydride bond and afforded the corresponding molybdenum cyclohexadienyl hydrides. For mono- and disubstituted arenes, the site selectivity for insertion of the most substituted bond increases with increasing size of the substituent from methyl to ethyl, iso-propyl, and tert-butyl. In contrast, 1,3,5-trisubstituted arenes underwent insertion with exclusive site selectivity. Relative rates of insertion were determined by competition experiments and established faster insertions for electron-rich arenes. Introduction of electron-withdrawing trifluoromethyl groups on the arene resulted in decreased relative rates of insertion and an increased rate for H2 reductive elimination, favoring formation of the corresponding molybdenum η6-arene complex. Studies on the reductive elimination of the cyclohexadienyl ligand with the hydride enabled the synthesis of an enantioenriched cyclohexa-1,3-diene. This study provides new insights into the ligand requirements for catalytic arene hydrogenation and a new strategy for selective arene reduction.

Density functional theory and CCSD(T) evaluation of ionization potentials, redox potentials, and bond energies related to zirconocene polymerization catalysts

Quantum-mechanical-based computational design of molecular catalysts requires accurate and fast electronic structure calculations to determine and predict properties of transition-metal complexes. For Zr-based molecular complexes related to polyethylene catalysis, previous evaluation of density functional theory (DFT) and wavefunction methods only examined oxides and halides or select reaction barrier heights. In this work, we evaluate the performance of DFT against experimental redox potentials and bond dissociation enthalpies (BDEs) for zirconocene complexes directly relevant to ethylene polymerization catalysis. We also examined the ability of DFT to compute the fourth atomic ionization potential of zirconium and the effect the basis set selection has on the ionization potential computed with CCSD(T). Generally, the atomic ionization potential and redox potentials are very well reproduced by DFT, but we discovered relatively large deviations of DFT-calculated BDEs compared to experiment. However, evaluation of BDEs with CCSD(T) suggests that experimental values should be revisited, and our CCSD(T) values should be taken as most accurate.

Hafnium vs. Zirconium, the Perpetual Battle for Supremacy in Catalytic Olefin Polymerization: A Simple Matter of Electrophilicity?

The performance of C2-symmetric ansa-hafnocene catalysts for isotactic polypropylene typically deteriorates at increasing temperature much faster than that of their zirconium analogues. Herein, we analyze in detail a set of five Hf/Zr metallocene pairs-including some of the latest generation catalysts-at medium- to high-polymerization temperature. Quantitative structure-activity relationship (QSAR) models for stereoselectivity, the ratio allyl/vinyl chain ends, and 2,1/3,1 misinsertions in the polymer indicate a strong dependence of polymerization performance on electrophilicity of the catalyst, which is a function of the ligand framework and the metal center. Based on this insight, the stronger performance decline of hafnocenes is ascribed to electrophilicity-dependent stabilization effects.

Probing β-alkyl elimination and selectivity in polyolefin hydrogenolysis through DFT

A long chain substrate with [(SiO)3ZrH] has been investigated to elucidate selectivity rules in β-alkyl elimination. DFT studies indicate that polypropylene preferentially undergoes β-Me elimination.

Comparison of the Reactivity and Structures for the Neutral and Cationic Bis(imino)pyridyl Iron and Cobalt Species by DFT Calculations

Density Functional Theory (DFT) method was adopted to investigate and compare the reaction mechanisms of ethylene polymerization catalyzed by neutral, cationic bis(imino)pyridyl (PDI) iron and cobalt derivatives. The electronic structure and the oxidation states of the metal center and the PDI ligand were analyzed by taking spin states, natural bond orbital (NBO) charge distribution, etc. into consideration, revealing that the reactivity is closely related to the valence electron numbers instead of the charge numbers. The neutral Co(0) had the lowest reactivity as it possessed the most electrons. During the formation of the cationic Co(+)/Fe(+), one electron was mainly lost from PDI ligand rather than the metal center while the metal center maintained +II valence state through the process. Moreover, a special unsymmetrically bidentate N^N coordination manner was found to provide the deficient metal surroundings with 14e, which may initiate the reactivity of some unsymmetrical species with rich electrons. Finally, an anion [AlMe4]− participating process was proposed to explain the presence of the experimentally observed LCo(+)B(C2H4). A special intermediate, Co(+)B(C2H4) [AlMe4]− with Co in +I and absence of Co–C σ bond, was obtained. These calculation results may provide fundamental information for further understanding and designing the ethylene polymerization catalysts.

An Integrated High Throughput Experimentation/Predictive QSAR Modeling Approach to ansa-Zirconocene Catalysts for Isotactic Polypropylene

Compared to heterogenous Ziegler-Natta systems (ZNS), ansa-metallocene catalysts for the industrial production of isotactic polypropylene feature a higher cost-to-performance balance. In particular, the C2-symmetric bis(indenyl) ansa-zirconocenes disclosed in the 1990s are complex to prepare, less stereo- and/or regioselective than ZNS, and lose performance at practical application temperatures. The golden era of these complexes, though, was before High Throughput Experimentation (HTE) could contribute significantly to their evolution. Herein, we illustrate a Quantitative Structure - Activity Relationship (QSAR) model trained on a robust and highly accurate HTE database. The clear-box QSAR model utilizes, in particular, a limited number of chemically intuitive 3D geometric descriptors that screen various regions of space in and around the catalytic pocket in a modular way thus enabling to quantify individual substituent contributions. The main focus of the paper is on the methodology, which should be of rather broad applicability in molecular organometallic catalysis. Then again, it is worth emphasizing that the specific application reported here led us to identify in a comparatively short time novel zirconocene catalysts rivaling or even outperforming all previous homologues which strongly indicates that the metallocene story is not over yet.

Reactivity Trends of Lewis Acidic Sites in Methylaluminoxane and Some of Its Modifications

The established model cluster (AlOMe)16(AlMe3)6 for methylaluminoxane (MAO) cocatalyst has been studied by density functional theory, aiming to rationalize the different behaviors of unmodified MAO and TMA-depleted MAO/BHT (TMA = trimethylaluminum; BHT = 2,6-di-tert-butyl-4-methylphenol), highlighted in previous experimental studies. The tendency of the three model Lewis acidic sites A-C to release neutral Al fragments (i.e., AlMe2R; R = Me or bht) or transient aluminum cations (i.e., [AlMeR]+) has been investigated both in the absence and in the presence of neutral N-donors. Sites C are most likely responsible for the activation capabilities of TMA-rich MAO, but TMA depletion destabilizes them, possibly inducing structural rearrangements. The remaining sites A and B, albeit of lower Lewis acidity, should be still able to release cationic Al fragments when TMA-depleted modified MAOs are treated with N-donors (e.g. [AlMe(bht)]+ from MAO/BHT). These findings provide tentative interpretations for earlier observations of donor-dependent ionization tendencies of MAO and MAO/BHT and how TMA depleted MAOs can still be potent activators.

C1-Symmetric Si-bridged (2-indenyl)(1-indenyl) ansa-metallocenes as efficient ethene/1-hexene copolymerization catalysts

Synthesis, ethene/1-hexene copolymerization screening, and QSAR study of 28 silicon bridged (2-indenyl)(1-indenyl) zirconocenes: new promising catalysts for LLDPE were identified.

On the limits of tuning comonomer affinity of ‘Spaleck-type’ ansa-zirconocenes in ethene/1-hexene copolymerization: a high-throughput experimentation/QSAR approach

A change in rate-limiting step imparts a natural limit for comonomer affinity of C₂-symmetric zirconocenes.

Extent of structural change during the reaction and its relationship to isoselectivity in polypropylene polymerization with ansa-zirconocene/borate catalyst: A computational study

The mechanism of isotactic polypropylene (iPP) polymerization with an (R,R)-ansa-zirconocene/borate catalyst system was analyzed using quantum chemistry (QC) calculations by focusing on the extent of structural change during monomer insertion. The activation energy for migratory insertion, E_a , was compared for four possible reaction paths with regard to monomer coordination, that is, 1,2-re, 1,2-si, 2,1-si, and 2,1-re, until the seventh monomer insertion step, explicitly including a borate anion cocatalyst. This indicated that the 1,2-re path was most favorable, except for the first step, which favored 1,2-si. As far as the first step, the product of 1,2-si is a conformational isomer to that of the 1,2-re path, and the exceptional favorability of 1,2-si does not affect the isoselectivity. These results support previous studies, except that our results address the unexplored seventh insertion step with a borate anion cocatalyst by QC calculations. The isoselectivity correlated with the extent of structural change in the whole system during the reaction. It was proved from our detail analysis that the advantage of 1,2-re with a small E_a is attributed to its smaller structural changes due to low steric repulsion in the system compared with other paths. Conversely, larger repulsion in the systems involved in other paths results in larger structural changes to minimize the structural strain. However, the relaxation appears insufficient due to structural restriction of the enforced four-membered ring transition state structure. A borate anion cocatalyst broke the C₂ symmetry of the electronic structures of zirconocene, resulting in an odd-even E_a frequency for the monomer insertion. Molecular orbital analysis demonstrated that the d-π orbital overlaps can explain the approach direction of the olefin coordination and the bent structure of zirconocene, providing a different viewpoint from previous studies. The potential for catalyst control was discussed based on our results. © 2019 Wiley Periodicals, Inc.

Selective Copper Complex-Catalyzed Hydrodefluorination of Fluoroalkenes and Allyl Fluorides: A Tale of Two Mechanisms

The transition to more economically friendly small-chain fluorinated groups is leading to a resurgence in the synthesis and reactivity of fluoroalkenes. One versatile method to obtain a variety of commercially relevant hydrofluoroalkenes involves the catalytic hydrodefluorination (HDF) of fluoroalkenes using silanes. In this work it is shown that copper hydride complexes of tertiary phosphorus ligands (L) can be tuned to achieve selective multiple HDF of fluoroalkenes. In one example, HDF of the hexafluoropropene dimer affords a single isomer of heptafluoro-2-methylpentene in which five fluorines have been selectively replaced with hydrogens. DFT computational studies suggest a distinct HDF mechanisms for L₂CuH (bidentate or bulky monodentate phosphines) and L₃CuH (small cone angle monodentate phosphines) catalysts, allowing for stereocontrol of the HDF of trifluoroethylene.

Separating Electronic from Steric Effects in Ethene/α-Olefin Copolymerization: A Case Study on Octahedral [ONNO] Zr-Catalysts

Four Cl/Me substituted [ONNO] Zr-catalysts have been tested in ethene/α-olefin polymerization. Replacing electron-donating methyl with isosteric but electron-withdrawing chlorine substituents results in a significant increase of comonomer incorporation. Exploration of steric and electronic properties of the ancillary ligand by DFT confirm that relative reactivity ratios are mainly determined by the electrophilicity of the metal center. Furthermore, quantitative DFT modeling of propagation barriers that determine polymerization kinetics reveals that electronic effects observed in these catalysts affect relative barriers for insertion and a capture-like transition state (TS).

A Benchmark of Density Functional Approximations For Thermochemistry and Kinetics of Hydride Reductions of Cyclohexanones

The performance of density functionals and wavefunction methods for describing the thermodynamics and kinetics of hydride reductions of 2-substituted cyclohexanones has been evaluated for the first time. A variety of exchange correlation functionals ranging from generalized gradient approximations to double hybrids have been tested and their performance to describe the facial selectivity of hydride reductions of cyclohexanones has been carefully assessed relative to the CCSD(T) method. Among the tested methods, an approach in which single-point energy calculations using the double hybrid B2PLYP-D3 functional on ωB97X-D optimized geometries provides the most accurate transition state energies for these kinetically-controlled reactions. Moreover, the role of torsional strain, temperature, solvation, noncovalent interactions on the stereoselectivity of these reductions was elucidated. Our results indicate a prominent role of the substituent on the cis/trans ratios driven by the delicate interplay between torsional strain and dispersion interactions.

Adsorption of Titanium Tetrachloride on Magnesium Dichloride Clusters

Magnesium dichloride and titanium tetrachloride are key components in heterogeneous Ziegler-Natta olefin polymerization catalysis. We have determined the preferred binding modes of titanium tetrachloride on magnesium dichloride clusters at the M06-2X level of theory, thus accounting for dispersion. A systematic study was carried out to locate the lowest energy isomers of (MgCl₂) _n (TiCl₄) _m complexes. Altogether ca. 1000 complexes with n ranging from 1 to 19 and m from 1 to 13 were studied. In line with the previous literature, the results consistently show that adsorption of TiCl₄ onto MgCl₂ preferably leads to octahedral six-coordination for both Mg and Ti atoms. This can be achieved by mononuclear binding of TiCl₄ and binuclear binding of Ti₂Cl₈ on (110) and (104) surfaces of MgCl₂, respectively. The preferred octahedral six-coordination is also achieved by binding of TiCl₄ on defect sites, of which the most striking example is trinuclear binding mode at corners of the clusters. Overall, the results highlight the relevance of multinuclear binding of TiCl₄ on MgCl₂.

Role(s) of TMA in polymerization

Quenched-flow data for propene polymerization with rac-Me2Si(2-Me-4-Ph-1-indenyl)2ZrCl2/MAO support a picture where removal of MAO qualitatively changes the kinetic profile from a mainly enthalpic to a mainly entropic barrier. DFT studies suggest that a not previously recognized singly-bridged end-on coordination mode of Me6Al2 to catalytically active centers may be kinetically relevant as a resting state. In contrast, the more traditional doubly-bridged complex of Me3Al is proposed to be more relevant to chain transfer to cocatalyst.

Olefin polymerisation catalysts: when perfection is not enough

Despite decades of thorough mechanistic investigations, it is still hard to predict the activity of a novel olefin polymerisation catalyst, even when the precursor is a well-defined molecular entity. In the present study, we highlight the crucial importance of activation entropy on the polymerisation rate and how weak interactions of the catalytic species with electron donating species in the reaction pool can ultimately lower the activation free energy.