Search NASA⌕ Search

Engineering topics

Harding, Lawrence B.

Publications and source records attributed to Harding, Lawrence B..

Competing Radical and Molecular Channels in the Unimolecular Dissociation of Methylformate

The thermal dissociation of methylformate has been characterized by numerous experimental and theoretical studies that all seem to agree that the lowest energetically accessible process is a 3-center H-atom transfer that leads to the molecular products CH3OH 3 OH + CO. However, these literature studies seem to be at odds with regards to the role of other competing molecular eliminations and bond-fission processes and therefore a complete and resolved mechanistic picture of methylformate thermal dissociation still eludes us. In this work, we have performed high-level electronic structure theory calculations to characterize the energetics of other overlooked molecular and radical processes that can originate on the complex CH3OCHO 3 OCHO potential energy surface (PES). The present calculations do indeed confirm that the lowest energy process accessible on this PES is molecular elimination to form CH3OH 3 OH + CO. However, unlike prior theoretical studies, the present calculations reveal that the second lowest energy process is a 5-center concerted elimination process that leads to the direct formation of H2 2 + CH2O 2 O + CO. We also note that HCO2 2 from the lowest-lying bond fission has two energetically comparable electronic states (of 2 B 2 and 2 A 1 symmetry). Furthermore, the barrier for H-atom migration of CH3OCHO 3 OCHO to form the carbene (CH3OCOH) 3 OCOH) is similar to that for the direct 4-center elimination process leading to CH2O 2 O + CH2O 2 O characterized in prior literature studies. Radical and molecular pathways to CH3 3 + HOCO and CH4 4 + CO2 2 can also be facilitated from this carbene. Master equation calculations were performed to characterize the competition between the various molecular and radical processes on this more elaborate CH3OCHO 3 OCHO PES. The results of the present theoretical analyses were used to resolve outstanding questions on the role of secondary radical- initiated reactions in characterizing prior literature experimental studies.

Chemical Kinetics↗

Substitution Reactions in the Pyrolysis of Acetone Revealed through a Modeling, Experiment, Theory Paradigm

The development of high-fidelity mechanisms for chemically reactive systems is a challenging process that requires the compilation of rate descriptions for a large and somewhat ill-defined set of reactions. The present unified combination of modeling, experiment, and theory provides a paradigm for improving such mechanism development efforts. Here we combine broadband rotational spectroscopy with detailed chemical modeling based on rate constants obtained from automated ab initio transition state theory-based master equation calculations and high-level thermochemical parametrizations. Broadband rotational spectroscopy offers quantitative and isomer-specific detection by which branching ratios of polar reaction products may be obtained. Using this technique, we observe and characterize products arising from H atom substitution reactions in the flash pyrolysis of acetone (CH 3 C(O)CH 3 ) at a nominal temperature of 1800 K. The major product observed is ketene (CH 2 CO). Minor products identified include acetaldehyde (CH 3 CHO), propyne (CH 3 CCH), propene (CH 2 CHCH 3 ), and water (HDO). Literature mechanisms for the pyrolysis of acetone do not adequately describe the minor products. The inclusion of a variety of substitution reactions, with rate constants and thermochemistry obtained from automated ab initio kinetics predictions and Active Thermochemical Tables analyses, demonstrates an important role for such processes. Furthermore, the pathway to acetaldehyde is shown to be a direct result of substitution of acetone's methyl group by a free H atom, while propene formation arises from OH substitution in the enol form of acetone by a free H atom.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reaction Profiles and Kinetics for Radical–Radical Hydrogen Abstraction via Multireference Coupled Cluster Theory

Radical-radical abstractions in hydrocarbon oxidation chemistry are disproportionation reactions that are generally exothermic with little or no barrier yet are underappreciated and poorly studied. Such challenging multireference electronic structure problems are tackled here using the recently developed state-specific multireference coupled cluster methods Mk-MRCCSD and Mk-MRCCSD(T), as well as the companion perturbation theory Mk-MRPT2 and the established MRCISD, MRCISD+Q, and CASPT2 approaches. Reaction paths are investigated for five prototypes involving radical-radical hydrogen abstraction: H + BeH → H 2 + Be, H + NH 2 → H 2 + NH, CH 3 + C 2 H 5 → CH 4 + C 2 H 4 , H + C2H 5 → H 2 + C 2 H 4 , and H + HCO → H 2 + CO. Full configuration interaction (FCI) benchmark computations for the H + BeH, H + NH 2 , and H + HCO reactions prove that Mk-MRCCSD(T) provides superior accuracy for the interaction energies in the entrance channel, with mean absolute errors less than 0.3 kcal mol(-1) and percentage deviations less than 10% over the fragment separations of relevance to kinetics. To facilitate combustion studies, energetics for the CH 3 + C 2 H 5 , H + C 2 H 5 , and H + HCO reactions were computed at each level of theory with correlation-consistent basis sets (cc-pVXZ, X = T, Q 5) and extrapolated to the complete basis set (CBS) limit. These CBS energies were coupled with CASPT2 projected vibrational frequencies along a minimum energy path to obtain rate constants for these three reactions. The rigorous Mk-MRCCSD(T)/CBS results demonstrate unequivocally that these three reactions proceed with no barrier in the entrance channel, contrary to some earlier predictions. Mk-MRCCSD(T) also reveals that the economical CASPT2 method performs well for large interfragment separations but may deteriorate substantially at shorter distances.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An improved long range potential for O(1D)+H2

Several dynamics studies have indicated that the product isotopic distributions and vibrational energy distributions in the reaction O(1D) + H2(HD) are very sensitive to the long range interaction of the reactants. In this paper, an improved calculation of the reactant region of the potential energy surface is reported. In agreement with previous work, no barrier is found to edge-on insertion (1 A-prime ground state water surface), but a smaller (less than 0.2 kcal/mol) barrier to collinear addition, 1 Sigma (+) surface, is found than in previous work. The long range potential obtained in the present work most closely resembles the SL3 surface.

Walch, Stephen P.↗