ReaxFF reactive force field for molecular dynamics simulations of hydrocarbon oxidation

To investigate the initial chemical events associated with high-temperature gas-phase oxidation of hydrocarbons, we have expanded the ReaxFF reactive force field training set to include additional transition states and chemical reactivity of systems relevant to these reactions and optimized the force field parameters against a quantum mechanics (QM)-based training set. To validate the ReaxFF potential obtained after parameter optimization, we performed a range of NVT-MD simulations on various hydrocarbon/O2 systems. From simulations on methane/O2, o-xylene/O2, propene/O2, and benzene/O2 mixtures, we found that ReaxFF obtains the correct reactivity trend (propene > o-xylene > methane > benzene), following the trend in the C-H bond strength in these hydrocarbons. We also tracked in detail the reactions during a complete oxidation of isolated methane, propene, and o-xylene to a CO/CO2/H2O mixture and found that the pathways predicted by ReaxFF are in agreement with chemical intuition and our QM results. We observed that the predominant initiation reaction for oxidation of methane, propene, and o-xylene under fuel lean conditions involved hydrogen abstraction of the methyl hydrogen by molecular oxygen forming hydroperoxyl and hydrocarbon radical species. While under fuel rich conditions with a mixture of these hydrocarbons, we observed different chemistry compared with the oxidation of isolated hydrocarbons including a change in the type of initiation reactions, which involved both decomposition of the hydrocarbon or attack by other radicals in the system. Since ReaxFF is capable of simulating complicated reaction pathways without any preconditioning, we believe that atomistic modeling with ReaxFF provides a useful method for determining the initial events of oxidation of hydrocarbons under extreme conditions and can enhance existing combustion models.     

Development of the ReaxFF reactive force field for describing transition metal catalyzed reactions, with application to the initial stages of the catalytic formation of carbon nanotubes

With the aim of developing a computationally inexpensive method for modeling the high-temperature reaction dynamics of transition metal catalyzed reactions we have developed a ReaxFF reactive force field in which the parameters are fitted to a substantial quantum mechanics (QM) training set, containing full reaction pathways for relevant reactions. In this paper we apply this approach to reactions involving carbon materials plus Co, Ni, and Cu atoms. We find that ReaxFF reproduces the QM reaction data with good accuracy while also reproducing the binding characteristics of Co, Ni, and Cu atoms to hydrocarbon fragments. To demonstrate the applicability of ReaxFF we performed high-temperature (1500 K) molecular dynamics simulations on a nonbranched all-carbon feedstock in the presence and absence of Co, Ni, and Cu atoms. We find that the presence of Co and Ni leads to substantial amounts of branched carbon atoms, leading eventually to the formation of carbon-nanotube-like species. In contrast, we find that under the same simulation conditions Cu leads to very little branching and leads to products with no nanotube character. In the absence of metals no branching is observed at all. These results suggest that Ni and Co catalyze the production of nanotube-like species whereas Cu does not. This is in excellent agreement with experimental observations, demonstrating that ReaxFF can provide a useful and computational tractable tool for studying the dynamics of transition metal catalytic chemistry.     

A detailed reaction mechanism for hexamethyldisiloxane combustion via experiments and ReaxFF molecular dynamics simulations

Hexamethyldisiloxane (HMDSO) is one of the main impurities in the syngas produced from sewage and landfill plants. In order to utilize this syngas or control the characteristics of the generated silica particles, it is crucial to understand the chemical kinetics of HMDSO combustion. This study investigated the process of HMDSO combustion using synchrotron radiation mass spectrometry (SRMS), gas chromatography (GC), and ReaxFF molecular dynamics simulations. First, the force field used for ReaxFF simulation was validated by comparing the energies of different bond lengths, bond angles, and dihedral angles with the ones from DFT calculations. Good agreements were found. Then, ReaxFF simulations of HMDSO combustion with this force field were conducted under various conditions, which include different equivalence ratios (0.67, 1.0, and 1.5) and temperatures ranging from 2000 to 3500 K. The oxidation characteristics of HMDSO were analyzed, including the evolution of gas products and particle formation. Finally, based on the results from experiments and ReaxFF simulations, the reaction pathways, reaction lists, and reaction kinetics data during HMDSO combustion were obtained. A detailed reaction mechanism was proposed and validated by applying it in modeling the H₂/HMDSO/O₂ combustion systems. The temperature and part of the gas products such as CO and CO₂ as well as SiO could be well predicted.     

RP-3高温氧化初始阶段反应机理的ReaxFF MD模拟

本文采用ReaxFF MD方法对一种较新的RP-3四组分替代燃料模型的高温氧化过程进行了研究。利用作者所在课题组研发的独特分析工具VARxMD,对燃烧过程中主要物种(燃料分子、O₂、C₂H₄、·CH₃)随时间和温度的演变规律及其化学反应进行了系统分析。ReaxFF MD模拟得到的燃料和氧气消耗量、乙烯和甲基自由基的生成量与相同温度和初始压力条件下CHEMKIN的计算结果处于同一量级,同时获得了详细的物质结构信息和反应列表。进一步对模拟得到的反应机理形式进行观察后发现,模拟获得的机理形式与文献中的描述一致。对燃料分子第一步反应数量的统计发现,其类型主要为攫氢反应和分子内断裂反应,且后者占主导;燃料分子第一步反应数量的统计也定性展现了不同燃烧条件下各类反应发生的可能性。对氧元素相关的反应分析发现,氧分子和C₁-C₃小分子发生的反应所占比例较大,能在一定程度上为机理简化提供有益线索。在对反应机理分析的基础上获得了RP-3四组分替代燃料体系高温氧化过程的化学反应网络。我们认为,ReaxFF MD反应分子动力学模拟、结合VARxMD对模拟结果深入分析的方法是有潜力系统认识燃料氧化反应机理的新方法,对构建燃料的燃烧反应机理库有一定的帮助。     

Atomistic-scale simulations of the initial chemical events in the thermal initiation of triacetonetriperoxide

To study the initial chemical events related to the detonation of triacetonetriperoxide (TATP), we have performed a series of molecular dynamics (MD) simulations. In these simulations we used the ReaxFF reactive force field, which we have extended to reproduce the quantum mechanics (QM)-derived relative energies of the reactants, products, intermediates, and transition states related to the TATP unimolecular decomposition. We find excellent agreement between the QM-predicted reaction products and those observed from 100 independent ReaxFF unimolecular MD cookoff simulations. Furthermore, the primary reaction products and average initiation temperature observed in these 100 independent unimolecular cookoff simulations match closely with those observed from a TATP condensed-phase cookoff simulation, indicating that unimolecular decomposition dominates the thermal initiation of the TATP condensed phase. Our simulations demonstrate that thermal initiation of condensed-phase TATP is entropy-driven (rather than enthalpy-driven), since the initial reaction (which mainly leads to the formation of acetone, O(2), and several unstable C(3)H(6)O(2) isomers) is almost energy-neutral. The O(2) generated in the initiation steps is subsequently utilized in exothermic secondary reactions, leading finally to formation of water and a wide range of small hydrocarbons, acids, aldehydes, ketones, ethers, and alcohols.     

Dynamics of the dissociation of hydrogen on stepped platinum surfaces using the ReaxFF reactive force field

The dissociation of hydrogen on eight platinum surfaces, Pt(111), Pt(100), Pt(110), Pt(211), Pt(311), Pt(331), Pt(332), and Pt(533), has been studied using molecular dynamics and the reactive force field, ReaxFF. The force field, which includes the degrees of freedom of the atoms in the platinum substrate, was used unmodified with potential parameters determined from previous calculations performed on a training set exclusive of the surfaces considered in this work. The energetics of the eight surfaces in the absence of hydrogen at 0 K were first compared to previous density functional theory (DFT) calculations and found to underestimate excess surface energy. However, taking Pt(111) as a reference state, we found that the trend between surfaces was adequately predicted to justify a relative comparison between the various stepped surfaces. To assess the strengths and weaknesses of the force field, we performed detailed simulations on two stepped surfaces, Pt(533) and Pt(211), and compared our findings to published experimental and theoretical results. In general, the absolute magnitude of reaction rate predictions was low, a result of the force field's tendency to underpredict surface energy. However, when normalized, the simulations show the correct linear scaling with incident energy and angular dependence at collision energies where a direct dissociation mechanism is observed. Because ReaxFF includes all degrees of freedom in the substrate, we carried out simulations aimed at understanding surface-temperature effects on Pt(533). On the basis of the results on Pt(533)/Pt(211), we studied the reaction of hydrogen at normal incidence on all eight surfaces in a range of energies where we anticipated the force field to give reasonable qualitative trends. These results were subsequently fit to a simple linear model that predicts the enhanced reactivity of surfaces containing 111-type atomic steps versus 100-type atomic steps. This model provides a simple framework for predicting high-energy/high-temperature kinetics of complex surfaces not vicinal to Pt(111).     

Challenges and perspectives of combustion chemistry research

In the past decades, combustion chemistry research grew rapidly due to the development of combustion diagnostic methods,quantum chemistry methods, kinetic theory, and computational techniques. A lot of kinetic models have been developed for fuels from hydrogen to transportation fuel surrogates. Besides, multi-scale research method has been widely adopted to develop comprehensive models, which are expected to cover combustion conditions in real combustion devices. However, critical gaps still remain between the laboratory research and real engine application due to the insufficient research work on high pressure and low temperature combustion chemistry. Besides, there is also a great need of predictive pollutant formation model. Further development of combustion chemistry research depends on a closer interaction of combustion diagnostics, theoretical calculation and kinetic model development. This paper summarizes the recent progress in combustion chemistry research briefly and outlines the challenges and perspectives.     

Effects of fuel additives on the thermal cracking of n-decane from reactive molecular dynamics

Thermal cracking of n-decane and n-decane in the presence of several fuel additives are studied in order to improve the rate of thermal cracking by using reactive molecular dynamics (MD) simulations employing the ReaxFF reactive force field. From MD simulations, we find the initiation mechanisms of pyrolysis of n-decane are mainly through two pathways: (1) the cleavage of a C-C bond to form smaller hydrocarbon radicals, and (2) the dehydrogenation reaction to form an H radical and the corresponding decyl radical. Another pathway is the H-abstraction reactions by small radicals including H, CH(3), and C(2)H(5). The basic reaction mechanisms are in good agreement with existing chemical kinetic models of thermal decomposition of n-decane. Quantum mechanical calculations of reaction enthalpies demonstrate that the H-abstraction channel is easier compared with the direct C-C or C-H bond-breaking in n-decane. The thermal cracking of n-decane with several additives is further investigated. ReaxFF MD simulations lead to reasonable Arrhenius parameters compared with experimental results based on first-order kinetic analysis. The different chemical structures of the fuel additives greatly affect the apparent activation energy and pre-exponential factors. The presence of diethyl ether (DEE), methyl tert-butyl ether (MTBE), 1-nitropropane (NP), 3,6,9-triethyl-3,6,9-trimethyl-1,2,4,5,7,8-hexaoxonane (TEMPO), triethylamine (TEA), and diacetonediperodixe (DADP) exhibit remarkable promoting effect on the thermal cracking rates, compared with that of pure n-decane, in the following order: NP > TEMPO > DADP > DEE (～MTBE) > TEA, which coincides with experimental results. These results demonstrate that reactive MD simulations can be used to screen for fuel additives and provide useful information for more comprehensive chemical kinetic model studies at the molecular level.     

On the theory of the CO+OH reaction, including H and C kinetic isotope effects

The effect of pressure, temperature, HD isotopes, and C isotopes on the kinetics of the OH+CO reaction are investigated using Rice-Ramsperger-Kassel-Marcus theory. Pressure effects are treated with a step-ladder plus steady-state model and tunneling effects are included. New features include a treatment of the C isotope effect and a proposed nonstatistical effect in the reaction. The latter was prompted by existing kinetic results and molecular-beam data of Simons and co-workers on incomplete intramolecular energy transfer to the highest vibrational frequency mode in HOCO(*). In treating the many kinetic properties two small customary vertical adjustments of the barriers of the two transition states were made. The resulting calculations show reasonable agreement with the experimental data on (1) the pressure and temperature dependence of the HD effect, (2) the pressure-dependent (12)C(13)C isotope effect, (3) the strong non-Arrhenius behavior observed at low temperatures, (4) the high-temperature data, and (5) the pressure dependence of rate constants in various bath gases. The kinetic carbon isotopic effect is usually less than 10 per mil. A striking consequence of the nonstatistical assumption is the removal of a major discrepancy in a plot of the k(OH+CO)k(OD+CO) ratio versus pressure. A prediction is made for the temperature dependence of the OD+CO reaction in the low-pressure limit at low temperatures.     

工业NiW/Al2O3催化剂上二苯并噻吩的加氢脱硫动力学

以二苯并噻吩(DBT)为含硫模型化合物, 在高压滴流床反应装置中,考察了工业NiW/Al2O3催化剂（RN-10）的加氢脱硫(HDS)动力学规律,研究了氢分压（1.5 MPa～4.5 MPa）、氢油体积比（150~700）、液体质量空速(15 h-1~60 h-1）、反应温度（280 ℃~380 ℃）等对DBT的HDS反应结果的影响。结果表明,当氢分压和氢油体积比较大时,两者变化对DBT的转化率基本无影响;温度对DBT的转化率影响较大,提高温度可有效提高DBT的转化率,随着温度的升高,DBT转化率的增加逐渐变缓。采用2级平推流反应动力学模型对不同温度实验数据进行了拟合,求得了不同温度的表观反应速率常数,模型的相关系数>0.989。活化能计算结果表明,RN-10催化剂在高反应温度区(>330 ℃)的DBT的HDS活化能明显低于较低温度时的活化能,分别为13.4 kJ/mol和121.4 kJ/mol。对于RN-10催化剂,不可单纯地通过提高反应温度来大幅度提高HDS转化率。     

High-pressure rate rules for alkyl + O2 reactions. 2. The isomerization, cyclic ether formation, and β-scission reactions of hydroperoxy alkyl radicals

The unimolecular reactions of hydroperoxy alkyl radicals (QOOH) play a central role in the low-temperature oxidation of hydrocarbons as they compete with the addition of a second O(2) molecule, which is known to provide chain-branching. In this work we present high-pressure rate estimation rules for the most important unimolecular reactions of the β-, γ-, and δ-QOOH radicals: isomerization to RO(2), cyclic ether formation, and selected β-scission reactions. These rate rules are derived from high-pressure rate constants for a series of reactions of a given reaction class. The individual rate expressions are determined from CBS-QB3 electronic structure calculations combined with canonical transition state theory calculations. Next we use the rate rules, along with previously published rate estimation rules for the reactions of alkyl peroxy radicals (RO(2)), to investigate the potential impact of falloff effects in combustion/ignition kinetic modeling. Pressure effects are examined for the reaction of n-butyl radical with O(2) by comparison of concentration versus time profiles that were obtained using two mechanisms at 10 atm: one that contains pressure-dependent rate constants that are obtained from a QRRK/MSC analysis and another that only contains high-pressure rate expressions. These simulations reveal that under most conditions relevant to combustion/ignition problems, the high-pressure rate rules can be used directly to describe the reactions of RO(2) and QOOH. For the same conditions, we also address whether the various isomers equilibrate during reaction. These results indicate that equilibrium is established between the alkyl, RO(2), and γ- and δ-QOOH radicals.     

Evidence of a Kinetic Isotope Effect in Nanoaluminum and Water Combustion

The normally innocuous combination of aluminum and water becomes violently reactive on the nanoscale. Research in the field of the combustion of nanoparticulate aluminum has important implications in the design of molecular aluminum clusters, hydrogen storage systems, as well as energetic formulations which could use extraterrestrial water for space propulsion. However, the mechanism that controls the reaction speed is poorly understood. While current models for micron‐sized aluminum water combustion reactions place heavy emphasis on diffusional limitations, as reaction scales become commensurate with diffusion lengths (approaching the nanoscale) reaction rates have long been suspected to depend on chemical kinetics, but have never been definitely measured. The combustion analysis of nanoparticulate aluminum with H₂O or D₂O is presented. Different reaction rates resulting from the kinetic isotope effect are observed. The current study presents the first‐ever observed kinetic isotope effect in a metal combustion reaction and verifies that chemical reaction kinetics play a major role in determining the global burning rate.     

Reduced combustion time model for methane in gas turbine flow fields

Computational fluid dynamics (CFD) modeling of the complex processes that occur within the burner of a gas turbine engine has become a critical step in the design process. However, due to computer limitations, it is very difficult to completely couple the fluid mechanics solver with the full combustion chemistry. Therefore, simplified chemistry models are required, and the topic of this research was to provide reduced chemistry models for CH4/O2 gas turbine flow fields to be integrated into CFD codes for the simulation of flow fields of natural gas-fueled burners. The reduction procedure for the CH4/O2 model utilized a response modeling technique wherein the full mechanism was solved over a range of temperatures, pressures, and mixture ratios to establish the response of a particular variable, namely the chemical reaction time. The conditions covered were between 1000 and 2500 K for temperature, 0.1 and 2 for equivalence ratio in air, and 0.1 and 50 atm for pressure. The kinetic time models in the form of ignition time correlations are given in Arrhenius-type formulas as functions of equivaience ratio, temperature, and pressure; or fuel-to-air ratio, temperature, and pressure. A single ignition time model was obtained for the entire range of conditions, and separate models for the low-temperature and high-temperature regions as well as for fuel-lean and rich cases were also derived. Predictions using the reduced model were verified using results from the full mechanism and empirical correlations from experiments. The models are intended for (but not limited to) use in CFD codes for flow field simulations of gas turbine combustors in which initial conditions and degree of mixedness of the fuel and air are key factors in achieving stable and robust combustion processes and acceptable emission levels. The chemical time model was utilized successfully in CFD simulations of a generic gas turbine combustor with four different cases with various levels of fuel-air premixing.     

Dependence of transition state structure on substrate: the intrinsic C-13 kinetic isotope effect is different for physiological and slow substrates of the ornithine decarboxylase reaction because of different hydrogen bonding structures

Ornithine decarboxylase is the first and the rate-controlling enzyme in polyamine biosynthesis; it decarboxylates l-ornithine to form the diamine putrescine. We present calculations performed using a combined quantum mechanical and molecular mechanical (QM/MM) method with the AM1 semiempirical Hamiltonian for the wild-type ornithine decarboxylase reaction with ornithine (the physiological substrate) and lysine (a "slow" substrate) and for mutant E274A with ornithine substrate. The dynamical method is variational transition state theory with quantized vibrations. We employ a single reaction coordinate equal to the carbon-carbon distance of the dissociating bond, and we find a large difference between the intrinsic kinetic isotope effect for the physiological substrate, which equals 1.04, and that for the slow substrate, which equals 1.06. This shows that, contrary to a commonly accepted assumption, kinetic isotope effects on slow substrates are not always good models of intrinsic kinetic isotope effects on physiological substrates. Furthermore, analysis of free-energy-based samples of transition state structures shows that the differences in kinetic isotope effects may be traced to different numbers of hydrogen bonds at the different transition states of the different reactions.     

High-pressure rate rules for alkyl + O2 reactions. 1. The dissociation, concerted elimination, and isomerization channels of the alkyl peroxy radical

The reactions of alkyl peroxy radicals (RO(2)) play a central role in the low-temperature oxidation of hydrocarbons. In this work, we present high-pressure rate estimation rules for the dissociation, concerted elimination, and isomerization reactions of RO(2). These rate rules are derived from a systematic investigation of sets of reactions within a given reaction class using electronic structure calculations performed at the CBS-QB3 level of theory. The rate constants for the dissociation reactions are obtained from calculated equilibrium constants and a literature review of experimental rate constants for the reverse association reactions. For the concerted elimination and isomerization channels, rate constants are calculated using canonical transition state theory. To determine if the high-pressure rate expressions from this work can directly be used in ignition models, we use the QRRK/MSC method to calculate apparent pressure and temperature dependent rate constants for representative reactions of small, medium, and large alkyl radicals with O(2). A comparison of concentration versus time profiles obtained using either the pressure dependent rate constants or the corresponding high-pressure values reveals that under most conditions relevant to combustion/ignition problems, the high-pressure rate rules can be used directly to describe the reactions of RO(2).     

Kinetic Evidence for a Non‐Langmuir‐Hinshelwood Surface Reaction: H/D Exchange over Pd Nanoparticles and Pd(111)

The mechanism of hydrogen recombination on a Pd(111) single crystal and well‐defined Pd nanoparticles is studied using pulsed multi‐molecular beam techniques and the H₂/D₂ isotope exchange reaction. The focus of this study is to obtain a microscopic understanding of the role of subsurface hydrogen in enhancing the associative desorption of molecular hydrogen. HD production from H₂ and D₂ over Pd is investigated using pulsed molecular beams, and the temperature dependence and reaction orders are obtained for the rate of HD production under various reaction conditions designed to modulate the amount of subsurface hydrogen present. The experimental data are compared to the results of kinetic modeling based on different mechanisms for hydrogen recombination. We found that under conditions where virtually no subsurface hydrogen species are present, the HD formation rate can be described exceptionally well by a classic Langmuir–Hinshelwood model. However, this model completely fails to reproduce the experimentally observed high HD formation rates and the reaction orders under reaction conditions where subsurface hydrogen is present. To analyze this phenomenon, we develop two kinetic models that account for the role of subsurface hydrogen. First, we investigate the possibility of a change in the reaction mechanism, where recombination of one subsurface and one surface hydrogen species (known as a breakthrough mechanism) becomes dominant when subsurface hydrogen is present. Second, we investigate the possibility of the modified Langmuir–Hinshelwood mechanism with subsurface hydrogen lowering the activation energy for recombination of two hydrogen species adsorbed on the surface. We show that the experimental reaction kinetics can be well described by both kinetic models based on non‐Langmuir–Hinshelwood‐type mechanisms.     

基于ReaxFF力场的对硝基苯酚臭氧氧化分子动力学模拟

本文采用基于ReaxFF反应力场的分子动力学方法(ReaxFF MD),利用自主研发的国际首个基于GPU加速的ReaxFF MD程序系统GMD-Reax和独特的化学反应分析工具VARx MD,探索臭氧氧化对硝基苯酚的反应机理。通过模拟考察了300 K恒温条件下臭氧氧化水中对硝基苯酚的行为,获得了酚结构开环、CO_2生成、主要自由基(·OH、·O_2、·O)及团簇型自由基的数量演变趋势,并可定性描述六元环开环和CO_2生成均遵循伪一级反应动力学规律。反应机理分析表明酚类分子在水溶液中被臭氧氧化的路径主要经过攫氢、六元环开环、碳链的氧化分解三个阶段,也揭示了自由基和团簇型自由基在臭氧降解对硝基苯酚时所发挥的重要作用。本工作是应用ReaxFF MD分子模拟方法对常温水环境下臭氧降解酚类污染物反应机理研究的一个尝试,可为深入认识该机理及相关的实验、理论研究提供一定的参考。     

Toluene combustion: reaction paths, thermochemical properties, and kinetic analysis for the methylphenyl radical + O2 reaction

Aromatic compounds such as toluene and xylene are major components of many fuels. Accurate kinetic mechanisms for the combustion of toluene are, however, incomplete, as they do not accurately model experimental results such as strain rates and ignition times and consistently underpredict conversion. Current kinetic mechanisms for toluene combustion neglect the reactions of the methylphenyl radicals, and we believe that this is responsible, in part, for the shortcomings of these models. We also demonstrate how methylphenyl radical formation is important in the combustion and pyrolysis of other alkyl-substituted aromatic compounds such as xylene and trimethylbenzene. We have studied the oxidation reactions of the methylphenyl radicals with O2 using computational ab initio and density functional theory methods. A detailed reaction submechanism is presented for the 2-methylphenyl radical + O2 system, with 16 intermediates and products. For each species, enthalpies of formation are calculated using the computational methods G3 and G3B3, with isodesmic work reactions used to minimize computational errors. Transition states are calculated at the G3B3 level, yielding high-pressure limit elementary rate constants as a function of temperature. For the barrierless methylphenyl + O2 and methylphenoxy + O association reactions, rate constants are determined from variational transition state theory. Multichannel, multifrequency quantum Rice-Ramsperger-Kassel (qRRK) theory, with master equation analysis for falloff, provides rate constants as a function of temperature and pressure from 800 to 2400 K and 1 x 10(-4) to 1 x 10(3) atm. Analysis of our results shows that the dominant pathways for reaction of the three isomeric methylphenyl radicals is formation of methyloxepinoxy radicals and subsequent ring opening to methyl-dioxo-hexadienyl radicals. The next most important reaction pathway involves formation of methylphenoxy radicals + O in a chain branching process. At lower temperatures, the formation of stabilized methylphenylperoxy radicals becomes significant. A further important reaction channel is available only to the 2-methylphenyl isomer, where 6-methylene-2,4-cyclohexadiene-1-one (ortho-quinone methide, o-QM) is produced via an intramolecular hydrogen transfer from the methyl group to the peroxy radical in 2-methylphenylperoxy, with subsequent loss of OH. The decomposition of o-QM to benzene + CO reveals a potentially important new pathway for the conversion of toluene to benzene during combustion. A number of the important products of toluene combustion proposed in this study are known to be precursors of polyaromatic hydrocarbons that are involved in soot formation. Reactions leading to the important unsaturated oxygenated intermediates identified in this study, and the further reactions of these intermediates, are not included in current aromatic oxidation mechanisms.     

Application of global kinetic models to HMX beta-delta transition and cookoff processes

The reduction of the number of reactions in kinetic models for both the HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine) beta-delta phase transition and thermal cookoff provides an attractive alternative to traditional multi-stage kinetic models due to reduced calibration effort requirements. In this study, we use the LLNL code ALE3D to provide calibrated kinetic parameters for a two-reaction bidirectional beta-delta HMX phase transition model based on Sandia instrumented thermal ignition (SITI) and scaled thermal explosion (STEX) temperature history curves, and a Prout-Tompkins cookoff model based on one-dimensional time to explosion (ODTX) data. Results show that the two-reaction bidirectional beta-delta transition model presented here agrees as well with STEX and SITI temperature history curves as a reversible four-reaction Arrhenius model yet requires an order of magnitude less computational effort. In addition, a single-reaction Prout-Tompkins model calibrated to ODTX data provides better agreement with ODTX data than a traditional multistep Arrhenius model and can contain up to 90% fewer chemistry-limited time steps for low-temperature ODTX simulations. Manual calibration methods for the Prout-Tompkins kinetics provide much better agreement with ODTX experimental data than parameters derived from differential scanning calorimetry (DSC) measurements at atmospheric pressure. The predicted surface temperature at explosion for STEX cookoff simulations is a weak function of the cookoff model used, and a reduction of up to 15% of chemistry-limited time steps can be achieved by neglecting the beta-delta transition for this type of simulation. Finally, the inclusion of the beta-delta transition model in the overall kinetics model can affect the predicted time to explosion by 1% for the traditional multistep Arrhenius approach, and up to 11% using a Prout-Tompkins cookoff model.