Abstract

Ethylene polymerization mediated by methyl-bis(cyclopentadienyl)-zirconium or zirconocene catalyst, [ZrCp2CH3](+), is one of the most popular catalytic reaction for polyethylene production. Rationalizing the major effects that control the polymer growth result in a challenge for computational studies. Through quantum chemical calculations, we characterized the zirconocene ethylene polymerization reaction mechanism: chain initiation (I; first ethylene insertion) [ZrCp2CH2CH2CH3](+), chain propagation (P; from second (P-1) to ninth (P-9) ethylene insertion) [ZrCp2(CH2)(20)CH3](+), and chain termination processes (T; beta-hydrogen elimination from P-5 or P-9) [ZrHCp2(H2C = CH(CH2)(18)CH3](+) are analyzed through the potential energy surface (PES) and reaction force analysis (RFA). The RFA approach involves pulling out the portion of an activation barrier that corresponds to distorting reactants into the geometries they adopt in a transition state structure until it reaches the structural relaxation toward the equilibrium geometry of the product. Because the interactions between the zirconocene and the ethylene molecule are influenced by a combination of several kinds of steric and electronic effects, it is indispensable to understanding these interactions in order to rationalize and predict in a quantitative manner the reaction barrier heights and the concomitant polymer growth. In the present work, we employ a simple procedure within the framework of the RFA and the density functional steric energy decomposition analysis (EDA) approach to quantitatively separate the different types of interactions; steric (Delta E-s), electrostatic (Delta E-e), and quantum (Delta E-q) effects in order to predict the impact of each factor on the course of the polymerization process as well as for the polymer control and design.