화학공학소재연구정보센터
Journal of Physical Chemistry A, Vol.102, No.2, 395-411, 1998
Statistical modeling of gas-phase organometallic reactions based on density functional theory: Ni++C3H8
Density functional theory in its B3LYP variant has been used to explore quantitative details of the adiabatic potential energy surface leading from Ni+ + C3H8 reactants through a deep N+(C3H8) well to NiC2H4+ + CH4 and NiC3H6+ + H-2 elimination products. The lowest energy path to CH4 elimination involves facile CC bond insertion followed by a high multicenter transition state (MCTS) leading directly to the exit-channel complex Ni+(C2H4)(CH4). The lowest energy path to H-2 elimination involves comparably facile secondary CH insertion followed by a comparably high MCTS leading directly to the Ni+(C3H6)(H-2) complex. Primary CH insertion leads to significantly higher barriers to both CH4 and H-2 elimination; in particular, beta-methyl migration is energetically very costly. These results support a mechanism significantly different from the stepwise mechanisms invoked earlier but the same as that found in recent calculations on the Fe+ + C3H8 reaction by Holthausen and Koch. The geometries suggest that agostic interactions are important in stabilizing the key MCTSs. We use the B3LYP geometry (moments of inertia) and harmonic vibrational frequencies at each stationary point to construct a detailed rate model of the reaction, applying RRKM theory to each reaction step on the adiabatic ground-state surface. A steady-state approximation holds well and leads to a simple parallel decay model for the long-lived Ni+(C3H8) complexes. By adjusting the energies of the key MCTSs downward by 5-7 kcal/mol from the values from B3LYP theory, we can explain the range of experimental time scales, the product branching fractions, total cross section vs kinetic energy, and deuterium isotope effects, Differential centrifugal effects arising from the substantial variation of the mass distribution along the reaction coordinates lead to a strong J-dependence of the Ni+(C3H8) decay rate and of product branching fractions as well. The resulting mechanistic picture indicates that at low energy only low-J complexes (formed at small impact parameter) can overcome the centrifugal barriers atop the MCTSs and produce elimination products. High-J complexes live as Ni+(C3H8) for nanoseconds-microseconds, repeatedly insert in CC and CH bonds, but eventually revert to Ni+ + C3H8 reactants. We suggest possible reasons why the new model cannot explain the bimodal kinetic energy release distribution observed by Bowers and co-workers in the NiC3H6+ + H-2 channel.