Use of Direct Dynamics Simulations to Determine Unimolecular Reaction Paths and Arrhenius Parameters for Large Molecules.

In a previous study (J. Chem. Phys.2008, 129, 094701) it was shown that for a large molecule, with a total energy much greater than its barrier for decomposition and whose vibrational modes are harmonic oscillators, the expressions for the classical Rice-Ramsperger-Kassel-Marcus (RRKM) (i.e., RRK) and classical transition-state theory (TST) rate constants become equivalent. Using this relationship, a molecule's unimolecular rate constants versus temperature may be determined from chemical dynamics simulations of microcanonical ensembles for the molecule at different total energies. The simulation identifies the molecule's unimolecular pathways and their Arrhenius parameters. In the work presented here, this approach is used to study the thermal decomposition of CH3-NH-CH═CH-CH3, an important constituent in the polymer of cross-linked epoxy resins. Direct dynamics simulations, at the MP2/6-31+G* level of theory, were used to investigate the decomposition of microcanonical ensembles for this molecule. The Arrhenius A and Ea parameters determined from the direct dynamics simulation are in very good agreement with the TST Arrhenius parameters for the MP2/6-31+G* potential energy surface. The simulation method applied here may be particularly useful for large molecules with a multitude of decomposition pathways and whose transition states may be difficult to determine and have structures that are not readily obvious.

[1]  Po-han Lin,et al.  Molecular Simulation of Cross-Linked Epoxy and Epoxy−POSS Nanocomposite , 2009 .

[2]  William L. Hase,et al.  An analytic function describing the H+C2H4?C2H5 potential energy surface , 1978 .

[3]  Don L. Bunker,et al.  Monte Carlo Calculation of Triatomic Dissociation Rates. I. N2O and O3 , 1962 .

[4]  W. Hase.  Theoretical Critical Configuration for Ethane Decomposition and Methyl Radical Recombination , 1972 .

[5]  S. Klippenstein,et al.  Predictive theory for the combination kinetics of two alkyl radicals. , 2006, Physical chemistry chemical physics : PCCP.

[6]  W. Hase.  Variational unimolecular rate theory , 1983 .

[7]  G. Schatz,et al.  Hyperthermal chemistry in the gas phase and on surfaces: Theoretical studies , 2004 .

[8]  Haobin Wang,et al.  A Model Multidimensional Analytic Potential Energy Function for the Cl- + CH3Br .fwdarw. ClCH3 + Br- Reaction , 1994 .

[9]  K. Song,et al.  Direct dynamics study of N-protonated diglycine surface-induced dissociation. Influence of collision energy , 2003, Journal of the American Society for Mass Spectrometry.

[10]  S. Meroueh,et al.  Direct Dynamics Simulations of Collision- and Surface-Induced Dissociation of N-Protonated Glycine. Shattering Fragmentation† , 2002 .

[11]  R. Marcus,et al.  On the Statistical Theory of Unimolecular Processes , 2007 .

[12]  Lawrence B Harding,et al.  Predictive theory for hydrogen atom-hydrocarbon radical association kinetics. , 2005, The journal of physical chemistry. A.

[13]  Chuh‐Yung Chen,et al.  Effect of reactive layered double hydroxides on the thermal and mechanical properties of LDHs/epoxy nanocomposites , 2007 .

[14]  W. Hase,et al.  Born–Oppenheimer Direct Dynamics Classical Trajectory Simulations , 2003 .

[15]  W. Hase,et al.  Role of computational chemistry in the theory of unimolecular reaction rates , 2005 .

[16]  M. Frisch,et al.  Ab initio potential energy curve for CH bond dissociation in methane , 1982 .

[17]  W. Hase.  Properties of variational transition states for association reactions , 1987 .

[18]  Upakarasamy Lourderaj,et al.  Theoretical and computational studies of non-RRKM unimolecular dynamics. , 2009, The journal of physical chemistry. A.

[19]  W. Hase.  The criterion of minimum state density in unimolecular rate theory. An application to ethane dissociation , 1976 .

[20]  Arindam Banerjee,et al.  Submitted for publication , 1981 .

[21]  D. M. Hirst,et al.  Thermal Rate Constants for H + CH3 CH4 Recombination. II. Comparison of Experiment and Canonical Variational Transition State Theory , 1987 .

[22]  Xiche Hu,et al.  Vectorization of the general Monte Carlo classical trajectory program VENUS , 1991 .

[23]  J. Simons,et al.  An experimental chemist's guide to ab initio quantum chemistry , 1991 .

[24]  W. Hase,et al.  Kinetics of Chemically Activated Isobutane and Neopentane from the 4358‐ and 3660‐Å Photolyses of Diazomethane with Propane and Isobutane , 1970 .

[25]  K. Fukui.  Formulation of the reaction coordinate , 1970 .

[26]  K. Song,et al.  Importance of shattering fragmentation in the surface-induced dissociation of protonated octaglycine , 2009, Journal of the American Society for Mass Spectrometry.

[27]  Hase,et al.  Comparison of molecular dynamics and variational transition-state-theory calculations of the rate constant for H-atom association with the diamond {111} surface. , 1995, Physical review. B, Condensed matter.

[28]  W. Hase,et al.  Unimolecular reaction dynamics : theory and experiments , 1996 .

[29]  H. W. Schranz,et al.  An efficient microcanonical sampling procedure for molecular systems , 1991 .

[30]  W. Hase,et al.  Thermal rate constant for H+CH3 → CH4 recombination. Comparison of quasiclassical trajectory and variational transition state theory , 1985 .

[31]  W. Hase,et al.  Potential energy surface and unimolecular dynamics of stretched n-butane. , 2008, The Journal of chemical physics.

[32]  Timothy K. Minton,et al.  DYNAMICS OF ATOMIC-OXYGEN-INDUCED POLYMER DEGRADATION IN LOW EARTH ORBIT , 2001 .

[33]  S. Klippenstein.  Variational optimizations in the Rice–Ramsperger–Kassel–Marcus theory calculations for unimolecular dissociations with no reverse barrier , 1992 .

[34]  W. Hase,et al.  Sensitivity of the H + CH3 → CH4 recombination rate constant to the shape of the CH stretching potential , 1984 .

[35]  Gilles H. Peslherbe,et al.  Monte Carlo Sampling for Classical Trajectory Simulations , 2007 .

[36]  W. Hase,et al.  Complete multidimensional analytic potential energy surface for chloride + chloroform SN2 nucleophilic substitution , 1990 .

[37]  Xiche Hu,et al.  Modification of the Duchovic–Hase–Schlegel potential energy function for H+CH3↔CH4. Comparison of canonical variational transition state theory, trajectory, and experimental association rate constants , 1991 .

[38]  J. Futrell,et al.  Surface-induced dissociation of peptide ions: Kinetics and dynamics , 2003, Journal of the American Society for Mass Spectrometry.

[39]  William L. Hase,et al.  DIRECT DYNAMICS SIMULATIONS OF REACTIVE SYSTEMS , 1998 .

[40]  W. Hase,et al.  Monte carlo sampling of a microcanonical ensemble of classical harmonic oscillators , 1980 .

[41]  Mark S. Gordon,et al.  General atomic and molecular electronic structure system , 1993, J. Comput. Chem..

[42]  M. Persico,et al.  Algorithms for sampling a quantum microcanonical ensemble of harmonic oscillators at potential minima and conical intersections. , 2011, The journal of physical chemistry. A.

[43]  Y. Zhuang,et al.  Higher-accuracy schemes for approximating the Hessian from electronic structure calculations in chemical dynamics simulations. , 2010, The Journal of chemical physics.

[44]  Upakarasamy Lourderaj,et al.  Direct dynamics simulations using Hessian-based predictor-corrector integration algorithms. , 2007, The Journal of chemical physics.

[45]  Jimmie C. Oxley,et al.  Evidence for homolytic decomposition of ammonium nitrate at high temperature , 1989 .

[46]  S. Levchik,et al.  Thermal decomposition, combustion and flame‐retardancy of epoxy resins—a review of the recent literature , 2004 .

[47]  William L. Hase,et al.  On non‐RRKM unimolecular kinetics: Molecules in general, and CH3NC in particular , 1973 .

[48]  W. Hase,et al.  A Dynamical Study of the H + CH3 CH4 Recombination Reaction , 1985 .

[49]  T. Yan,et al.  A PM3-SRP + Analytic Function Potential Energy Surface Model for O(3P) Reactions with Alkanes. Application to O(3P) + Ethane† , 2004 .