Counteranion Effects on the Zirconocene Polymerization Catalyst Olefin Complex from QM/MM Molecular Dynamics Simulations

Cocatalyst effects in Hf-catalysed olefin polymerization: taking well-defined Al-alkyl borate salts into account

Hafnium catalysts for olefin polymerization are often very sensitive to the nature of cocatalysts, especially if they contain "free" aluminium trialkyls. Herein, cocatalyst effects in Hf-catalysed propene polymerization are examined for four Hf catalysts belonging to the family of CS-symmetric (Hf-CS-Met) and C2-symmetric (Hf-C2-Met) metallocenes, as well as of octahedral (Hf-OOOO) and pentacoordinated (Hf-PyAm) "post-metallocenes". The performance of the recently developed {[iBu2(PhNMe2)Al]2(μ-H)}+[B(C6F5)4]- (AlHAl) cocatalyst is compared with that of established systems like methylalumoxane, phenol-modified methylalumoxane and trityl borate/tri-iso-butylaluminium. The worst catalytic performance is observed with MAO. Conversely, the best cocatalyst varies depending on the Hf catalyst used and the performance indicator of interest, highlighting the complexity and importance of selecting the right precatalyst/cocatalyst combination. AlHAl proved to be a suitable system for all catalysts tested and, in some cases, it provides the best performance in terms of productivity (e.g. with hafnocenes). Furthermore, it generally leads to high molecular weight polymers, also with catalysts enabling easy chain transfer to Al like Hf-PyAm. This suggests that AlHAl has a low tendency to form heterodinuclear adducts with the cationic active species, therefore preventing the formation of dormant sites and/or termination events by chain transfer to Al.

(Pyridylamido)Hf(IV)-Catalyzed 1-Octene Polymerization Reaction Interwoven with the Structural Dynamics of the Ion-Pair-Active Species: Bridging from Microscopic Simulation to Chemical Kinetics with the Red Moon Method

We performed the atomistic simulation of 1-octene polymerization reaction catalyzed by the ionic pair (IP) consisting of the cationic active species of (pyridylamido)Hf(IV) catalyst, HfCat^Pn+, and different counteranions (CAs), B(C₆F₅)₄^- and MeB(C₆F₅)₃^-, at different monomer concentrations. Using a hybrid Monte Carlo/molecular dynamics method, that is, the Red Moon (RM) method, the reaction progress measured by the "RM cycle" was transformed into effective real time using the time transformation theory. Then, the degree of polymerization was found to be consistent with that in the chemical kinetics, a macroscopic theory, and experimental ones. Remarkably, the current simulation has revealed the different dynamical features in the polymerization behavior originating from the CA. Namely, the HfCat^Pn+-B(C₆F₅)₄^- IP mainly forms an outer-sphere IP (OSIP) throughout the polymerization. The HfCat^Pn+-MeB(C₆F₅)₃^- IP, on the other hand, forms an inner-sphere IP (ISIP) in the initial stage of polymerization, and the ratio of ISIP steeply drops after the first monomer insertion because the IP interaction is reduced by the steric hindrance between the inserted monomers and the CA. In conclusion, we have shown that the microscopic IP dynamics interwoven with the polymerization reaction can be computationally observed in the real-time domain by using the RM method. Therefore, our current work demonstrates the promising potential of the RM method in studying catalytic olefin polymerization and complex chemical reaction systems.

Solution-Phase Conformational/Vibrational Anharmonicity in Comonomer Incorporation Polyolefin Catalysis

The prediction of comonomer incorporation statistics in polyolefin catalysis necessitates an accurate calculation of free energies corresponding to monomer binding and insertion, often requiring sub-kcal/mol resolution to resolve experimental free energies. Batch reactor experiments are used to probe incorporation statistics of ethene and larger α-olefins for three constrained geometry complexes which are employed as model systems. Herein, over 6 ns of quantum mechanics/molecular mechanics (QM/MM) molecular dynamics is performed in combination with the zero-temperature string method to characterize the solution-phase insertion barrier and to analyze the contributions from conformational and vibrational anharmonicity arising both in vacuum and in solution. Conformational sampling in the solution-phase results in 0-2 kcal/mol corrections to the insertion barrier which are on the same scale necessary to resolve experimental free energies. Anharmonic contributions from conformational sampling in the solution phase are crucial energy contributions missing from static density functional theory calculations and implicit solvation models, and the accurate calculation of these contributions is a key step toward the quantitative prediction of comonomer incorporation statistics.

Polymerization via Insertion of Ethylene into Al-C bond under Mild Conditions: Mechanistic Studies on the Promotion Exerted by a Catalytic Amount of Cationic Gadolinium Metallocene

The cationic gadolinium metallocene [(C₅ Me₅ )₂ Gd][B(C₆ F₅ )₄ ], when combined with an excess amount of Al(ⁱ Bu)₃ , efficiently produces polyethylene at 80 °C under 0.8 MPa pressure of ethylene. After quenching, the resulting polyethylene has ethyl group at one end and isobutyl group at the other terminal. Because no Gd-alkyl species appears to be involved, a mechanism with conventional coordinative chain transfer polymerization (CCTP) is not feasible. Density functional theory (DFT) analyses indicate a novel mechanism in which the cationic Gd plays a crucial role by coordinating ethylene and assists the insertion of the coordinated ethylene into Al-C bond.

Atomistic Simulation of the Polymerization Reaction by a (Pyridylamido)hafnium(IV) Catalyst: Counteranion Influence on the Reaction Rate and the Living Character of the Catalytic System

Atomistic simulation of the 1-octene polymerization reaction by a (pyridylamido)Hf(IV) catalyst was conducted on the basis of Red Moon (RM) methodology, focusing on the effect of the counteranions (CAs), MeB(C₆F₅)₃^-, and B(C₆F₅)₄^-, on the catalyst activity and chain termination reaction. We show that RM simulation reasonably reproduces the faster reaction rate with B(C₆F₅)₄^- than with MeB(C₆F₅)₃^-. Notably, the initiation of the polymerization reaction with MeB(C₆F₅)₃^- is comparatively slow due to the difficulty of the first insertion. Then, we investigated the free energy map of the ion pair (IP) structures consisting of each CA and the cationic (pyridylamido)Hf(IV) catalyst with the growing polymer chain (HfCat^Pn+), which determines the polymerization reaction rates, and found that HfCat^Pn+-MeB(C₆F₅)₃^- can keep forming "inner-sphere" IPs even after the polymer chain becomes sufficiently bulky, while HfCat^Pn+-B(C₆F₅)₄^- forms mostly "outer-sphere" IPs. Finally, we further tried to elucidate the origin of the broader molecular weight distribution (MWD) of the polymer experimentally produced with B(C₆F₅)₄^- than that with MeB(C₆F₅)₃^-. Then, through the trajectory analysis of the RM simulations, it was revealed that the chain termination reaction would be more sensitive to the IP structures than the monomer insertion reaction because the former involves a more constrained structure than the latter, which is likely to be a possible origin of the MWDs dependent on the CAs.

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.

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.

Atomistic chemical computation of Olefin polymerization reaction catalyzed by (pyridylamido)hafnium(IV) complex: Application of Red Moon simulation

We have realized the microscopic simulation of olefin polymerization, that is, the simulation of the catalytic polymerization (CP) reaction system composed of (pyridylamido)hafnium(IV) complex as the catalyst. For this purpose, we adopted Red Moon (RM) method, a novel molecular simulation method to simulate the complex reaction system. First, according to the previous research, with the help of the QM calculation, we proposed a model system and elementary processes and explained the theoretical treatment of the simulation by the RM method (the RM simulation). In addition, we also proposed a macroscopic simulation based on chemical kinetics simulation. Then, we performed two simulations and compared them in terms of the effective time evolution of the three macroscopic physical quantities, the number-average molecular weight M_n , the mass-average molecular weight M_w , and the molar-mass dispersity Đ_M . The comparison showed that the two simulations are in quantitative or partially qualitative agreement with each other. Therefore, it is concluded that the RM simulation could not only simulate the CP reaction process microscopically, but also it is connected essentially to reproduce the time evolution of the macroscopic physical quantities on the basis of its microscopic simulation data. © 2018 Wiley Periodicals, Inc.

First-Principles Molecular Dynamics Studies of Organometallic Complexes and Homogeneous Catalytic Processes

Computational chemistry is a valuable aid to complement experimental studies of organometallic systems and their reactivity. It allows probing mechanistic hypotheses and investigating molecular structures, shedding light on the behavior and properties of molecular assemblies at the atomic scale. When approaching a chemical problem, the computational chemist has to decide on the theoretical approach needed to describe electron/nuclear interactions and the composition of the model used to approximate the actual system. Both factors determine the reliability of the modeling study. The community dedicated much effort to developing and improving the performance and accuracy of theoretical approaches for electronic structure calculations, on which the description of (inter)atomic interactions rely. Here, the importance of the model system used in computational studies is highlighted through examples from our recent research focused on organometallic systems and homogeneous catalytic processes. We show how the inclusion of explicit solvent allows the characterization of molecular events that would otherwise not be accessible in reduced model systems (clusters). These include the stabilization of nascent charged fragments via microscopic solvation (notably, hydrogen bonding), transfer of charge (protons) between distant fragments mediated by solvent molecules, and solvent coordination to unsaturated metal centers. Furthermore, when weak interactions are involved, we show how conformational and solvation properties of organometallic complexes are also affected by the explicit inclusion of solvent molecules. Such extended model systems may be treated under periodic boundary conditions, thus removing the cluster/continuum (or vacuum) boundary, and require a statistical mechanics simulation technique to sample the accessible configurational space. First-principles molecular dynamics, in which atomic forces are computed from electronic structure calculations (namely, density functional theory), is certainly the technique of choice to investigate chemical events in solution. This methodology is well established and thanks to advances in both algorithms and computational resources simulation times required for the modeling of chemical events are nowadays accessible, though the computational requirements use to be high. Specific applications reviewed here include mechanistic studies of the Shilov and Wacker processes, speciation in Pd chemistry, hydrogen bonding to metal centers, and the dynamics of agostic interactions.

The CHARMM-TURBOMOLE interface for efficient and accurate QM/MM molecular dynamics, free energies, and excited state properties

The quantum mechanical (QM)/molecular mechanical (MM) interface between Chemistry at HARvard Molecular Mechanics (CHARMM) and TURBOMOLE is described. CHARMM provides an extensive set of simulation algorithms, like molecular dynamics (MD) and free energy perturbation, and support for mature nonpolarizable and Drude polarizable force fields. TURBOMOLE provides fast QM calculations using density functional theory or wave function methods and excited state properties. CHARMM-TURBOMOLE is well-suited for extended QM/MM MD simulations using first principles methods with large (triple-ζ) basis sets. We demonstrate these capabilities with a QM/MM simulation of Mg(2+) (aq), where the MM outer sphere water molecules are represented using the SWM4-NDP Drude polarizable force field and the ion and inner coordination sphere are represented using QM PBE, PBE0, and MP2 methods. The relative solvation free energies of Mg(2+) and Zn(2+) were calculated using thermodynamic integration. We also demonstrate the features for excited state properties. We calculate the time-averaged solution absorption spectrum of indole, the emission spectrum of the indole 1La excited state, and the electronic circular dichroism spectrum of an oxacepham.