The influence of the counterion B(C6F5)3CH3– and solvent effects on the propagation and termination steps of ethylene polymerization catalyzed by Cp2ZrR+ (R=Me,Pr). A density functional study

(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.

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.

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.

Free energy profiles for monomer capture in Grubbs- and SHOP-type olefin polymerization catalysts: a constraint ab initio molecular dynamics study

Density functional theory together with Car-Parrinello ab initio molecular dynamics simulation has been used to investigate the free energy profiles (FEP) of monomer capture in Grubbs- and SHOP-type olefin polymerization catalysts. The FEPs along the reaction coordinates at 300 K were determined directly by a point wise thermodynamic integration technique. Comparison between potential energy profile (PEP) and the FEP has been made. The results show that, for both catalysts, the PEP for the monomer ethylene uptake by the metal center is a typical Morse curve without energy barrier. However, a small barrier (1.8 kcal/mol for Grubbs catalyst and 2.4 kcal/mol for SHOP catalyst) exists on the FEP. The pi complexation energy on the FES at 300 K is higher by 10-12 kcal/mol over that on the PES. The differences between FES and PES are due to entropy contribution. Slow growth simulations on the ethylene capture process show that the ethylene attacks the metal center by an asynchronous mode. This indicates that the forming of the pi-bonding between the metal and ethylene is initiated by electrophilic attack of the metal to one of the ethylene carbons.