Machine Learning Models for Predicting Zirconocene Properties and Barriers.

Zr metallocenes have significant potential to be highly tunable polyethylene catalysts through modification of the aromatic ligand framework. Here we report the development of multiple machine learning models using a large library (>700 systems) of DFT-calculated zirconocene properties and barriers for ethylene polymerization. We show that very accurate machine learning models are possible for HOMO-LUMO gaps of precatalysts but the performance significantly depends on the machine learning algorithm and type of featurization, such as fingerprints, Coulomb matrices, smooth overlap of atomic positions, or persistence images. Surprisingly, the description of the bonding hapticity, the number of direct connections between Zr and the ligand aromatic carbons, only has a moderate influence on the performance of most models. Despite robust models for HOMO-LUMO gaps, these types of machine learning models based on structure connectivity type features perform poorly in predicting ethylene migratory insertion barrier heights. Therefore, we developed several relatively robust and accurate machine learning models for barrier heights that are based on quantum-chemical descriptors (QCDs). The quantitative accuracy of these models depends on which potential energy surface structure QCDs were harvested from. This revealed a Hammett-type principle to naturally emerge showing that QCDs from the π-coordination complexes provide much better descriptions of the transition states than other potential-energy structures. Feature importance analysis of the QCDs provides several fundamental principles that influence zirconocene catalyst reactivity.

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.

Active Sites in a Heterogeneous Organometallic Catalyst for the Polymerization of Ethylene.

Heterogeneous derivatives of catalysts discovered by Ziegler and Natta are important for the industrial production of polyolefin plastics. However, the interaction between precatalysts, alkylaluminum activators, and oxide supports to form catalytically active materials is poorly understood. This is in contrast to homogeneous or model heterogeneous catalysts that contain resolved molecular structures that relate to activity and selectivity in polymerization reactions. This study describes the reactivity of triisobutylaluminum with high surface area aluminum oxide and a zirconocene precatalyst. Triisobutylaluminum reacts with the zirconocene precatalyst to form hydrides and passivates -OH sites on the alumina surface. The combination of passivated alumina and zirconium hydrides formed in this mixture generates ion pairs that polymerize ethylene.

Reversal of polarization in amidophosphines: neutral- and anionic-kappaP coordination vs. anionic-kappaP,N coordination and the formation of nickelaazaphosphiranes.

Nickel(ii) chloride reacts with the bis(tert-butylamino)diazadiphosphetidine {Bu(t)(H)NP(micro-NBu(t))(2)PN(H)Bu(t)} to form trans-[{Bu(t)(H)NP(micro-NBu(t))(2)PN(H)Bu(t)}(2)NiCl(2)]. In solution and the solid-state each heterocyclic ligand coordinates nickel through one phosphorus atom only. For comparison the solid-state structure of the known trans-[NiCl(2)(PEt(3))(2)] was also determined and it was found that the two complexes have almost identical bond parameters about nickel. The nickel-amidophosphine complexes [{Bu(t)OP(micro-NBu(t))(2)PNBu(t)}NiCl(PBu(n)(3))], [(PBu(n)(3))ClNi{Bu(t)NP(micro-NBu(t))(2)PNBu(t)}NiCl(PBu(n)(3))], and [{Me(2)Si(micro-NBu(t))(2)PNBu(t)}NiCl(PBu(n)(3))] were synthesized and X-ray structurally characterized. In these mono- and di-nuclear nickel complexes the nickel ions are coordinated in pseudo square-planar fashions, by one trialkylphosphine ligand, one chloride ligand and one kappaP,N-coordinated amidophosphine moiety from tert-butylamido-substituted heterocycles. Attempts to create nickel complexes chelated in a kappa(2)P fashion by the o-phenylenediamine-tethered mono- and di-anionic 1-{Me(2)Si(micro-NBu(t))(2)PN} 2-{Me(2)Si(micro-NBu(t))(2)PNH}C(6)H(4) and 1,2-{Me(2)Si(micro-NBu(t))(2)PN}C(6)H(4), respectively, afforded instead [1,2-{Me(2)Si(micro-NBu(t))(2)PN}{Me(2)Si(micro-NBu(t))(2)PN}C(6)H(4)NiCl] and [1,2-{Me(2)Si(micro-NBu(t))(2)PN}{Me(2)Si(micro-NBu(t))(2)PN}C(6)H(4)Ni{PEt(3)}], each complex having kappaP,N and kappaP coordinated amidophosphine ligands.

Inorganic heterocyclic bis(phosphines): syntheses and structures of a 1,2-bis(diazasilaphosphetidino)ethane and its nickel, molybdenum, and rhodium complexes.

Synthesis and characterization of a new, highly electron-rich, chelating bis(phosphine), based on the ethanediyl-linked inorganic heterocycle [Me(2)Si(mu-N(t)Bu)(2)P], are reported. Treatment of nickel chloride with this bis(phosphine) afforded square-planar cis-[[Me(2)Si(mu-N(t)Bu)(2)PCH(2)](2)NiCl(2)], which features isometric nickel-chloride (2.2220(8) A) and nickel-phosphorus (2.1572(8) A) bonds. The ligand reacted with cis-[(piperidine)(2)Mo(CO)(4)] to form colorless cis-[[Me(2)Si(mu-N(t)Bu)(2)PCH(2)](2)Mo(CO)(4)], which has distorted octahedral geometry and long Mo-P bonds (2.5461(18) A). Because of its potential applications in hydrogenation catalysis cis-[[Me(2)Si(mu-N(t)()Bu)(2)PCH(2)](2)Rh(COD)]BF(4) was synthesized. This square-planar, cationic rhodium(I) complex, having symmetrical Rh-P (2.250(2) A) and Rh-C (2.305(6) A) bonds, is structurally related to bis(phospholano)- and bis(phosphetano)rhodium species.