Methanediol from cloud-processed formaldehyde is only a minor source of atmospheric formic acid.

Atmospheric formic acid is severely underpredicted by models. A recent study proposed that this discrepancy can be resolved by abundant formic acid production from the reaction (1) between hydroxyl radical and methanediol derived from in-cloud formaldehyde processing and provided a chamber-experiment-derived rate constant, k1 = 7.5 × 10-12 cm3 s-1. High-level accuracy coupled cluster calculations in combination with E,J-resolved two-dimensional master equation analyses yield k1 = (2.4 ± 0.5) × 10-12 cm3 s-1 for relevant atmospheric conditions (T = 260-310 K and P = 0-1 atm). We attribute this significant discrepancy to HCOOH formation from other molecules in the chamber experiments. More importantly, we show that reversible aqueous processes result indirectly in the equilibration on a 10 min. time scale of the gas-phase reaction [Formula: see text] (2) with a HOCH2OH to HCHO ratio of only ca. 2%. Although HOCH2OH outgassing upon cloud evaporation typically increases this ratio by a factor of 1.5-5, as determined by numerical simulations, its in-cloud reprocessing is shown using a global model to strongly limit the gas-phase sink and the resulting production of formic acid. Based on the combined findings in this work, we derive a range of 1.2-8.5 Tg/y for the global HCOOH production from cloud-derived HOCH2OH reacting with OH. The best estimate, 3.3 Tg/y, is about 30 times less than recently reported. The theoretical equilibrium constant Keq (2) determined in this work also allows us to estimate the Henry's law constant of methanediol (8.1 × 105 M atm-1 at 280 K).

Refining the thermochemical properties of CF, SiF, and their cations by combining photoelectron spectroscopy, quantum chemical calculations, and the Active Thermochemical Tables approach.

Fluorinated species have a pivotal role in semiconductor material chemistry and some of them have been detected beyond the Earth's atmosphere. Achieving good energy accuracy on fluorinated species using quantum chemical calculations has long been a challenge. In addition, obtaining direct experimental thermochemical quantities has also proved difficult. Here, we report the threshold photoelectron and photoion yield spectra of SiF and CF radicals generated with a fluorine reactor. The spectra were analysed with the support of ab initio calculations, resulting in new experimental values for the adiabatic ionisation energies of both CF (9.128 ± 0.006 eV) and SiF (7.379 ± 0.009 eV). Using these values, the underlying thermochemical network of Active Thermochemical Tables was updated, providing further refined enthalpies of formation and dissociation energies of CF, SiF, and their cationic counterparts.

Active Thermochemical Tables: Enthalpies of Formation of Bromo- and Iodo-Methanes, Ethenes and Ethynes.

The thermochemistry of halocarbon species containing iodine and bromine is examined through an extensive interplay between new Feller-Peterson-Dixon (FPD) style composite methods and a detailed analysis of all available experimental and theoretical determinations using the thermochemical network that underlies the Active Thermochemical Tables (ATcT). From the computational viewpoint, a slower convergence of the components of composite thermochemistry methods is observed relative to species that solely contain first row elements, leading to a higher computational expense for achieving comparable levels of accuracy. Potential systematic sources of computational uncertainty are investigated, and, not surprisingly, spin-orbit coupling is found to be a critical component, particularly for iodine containing molecular species. The ATcT analysis of available experimental and theoretical determinations indicates that prior theoretical determinations have significantly larger uncertainties than originally reported, particularly in cases where molecular spin-orbit effects were ignored. Accurate and reliable heats of formation are reported for 38 halogen containing systems, based on combining the current computations with previous experimental and theoretical work via the ATcT approach.

Sub 20 cm-1 computational prediction of the CH bond energy - a case of systematic error in computational thermochemistry.

The bond dissociation energy of methylidyne, D0(CH), is studied using an improved version of the High-Accuracy Extrapolated ab initio Thermochemistry (HEAT) approach as well as the Feller-Peterson-Dixon (FPD) model chemistry. These calculations, which include basis sets up to nonuple (aug-cc-pCV9Z) quality, are expected to be capable of providing results substantially more accurate than the ca. 1 kJ mol-1 level that is characteristic of standard high-accuracy protocols for computational thermochemistry. The calculated 0 K CH bond energy (27 954 ± 15 cm-1 for HEAT and 27 956 ± 15 cm-1 for FPD), along with equivalent treatments of the CH ionization energy and the CH+ dissociation energy (85 829 ± 15 cm-1 and 32 946 ± 15 cm-1, respectively), were compared to the existing benchmarks from Active Thermochemical Tables (ATcT), uncovering an unexpected difference for D0(CH). This has prompted a detailed reexamination of the provenance of the corresponding ATcT benchmark, allowing the discovery and subsequent correction of a systematic error present in several published high-level calculations, ultimately yielding an amended ATcT benchmark for D0(CH). Finally, the current theoretical results were added to the ATcT Thermochemical Network, producing refined ATcT estimates of 27 957.3 ± 6.0 cm-1 for D0(CH), 32 946.7 ± 0.6 cm-1 for D0(CH+), and 85 831.0 ± 6.0 cm-1 for IE(CH).

Ring-Opening Dynamics of the Cyclopropyl Radical and Cation: the Transition State Nature of the Cyclopropyl Cation.

We provide compelling experimental and theoretical evidence for the transition state nature of the cyclopropyl cation. Synchrotron photoionization spectroscopy employing coincidence techniques together with a novel simulation based on high-accuracy ab initio calculations reveal that the cation is unstable via its allowed disrotatory ring-opening path. The ring strains of the cation and the radical are similar, but both ring opening paths for the radical are forbidden when the full electronic symmetries are considered. These findings are discussed in light of the early predictions by Longuet-Higgins alongside Woodward and Hoffman; we also propose a simple phase space explanation for the appearance of the cyclopropyl photoionization spectrum. The results of this work allow the refinement of the cyclopropane C-H bond dissociation energy, in addition to the cyclopropyl radical and cation cyclization energies, via the Active Thermochemical Tables approach.

Mechanism, thermochemistry, and kinetics of the reversible reactions: C2H3 + H2 ⇌ C2H4 + H ⇌ C2H5.

High-level coupled cluster theory, in conjunction with Active Thermochemical Tables (ATcT) and E,J-resolved master equation calculations, was used in a study of the title reactions, which play an important role in the combustion of hydrocarbons. In the set of radical/radical reactions leading to soot formation in flames, the addition of H-atoms to alkenes is likely a common reaction, triggering the isomerization of complex hydrocarbons to aromatics. The heats of formation of C2H3, C2H4, and C2H5 are established to be 301.26 ± 0.30 at 0 K (297.22 ± 0.30 at 298 K), 60.89 ± 0.11 (52.38 ± 0.11), and 131.38 ± 0.22 (120.63 ± 0.22) kJ mol-1, respectively. The calculated rate constants from first principles agree well with experiments where they are available. Under conditions typical of high temperature combustion - where experimental work is very challenging with a consequent dearth of accurate data - we provide high-level theoretical results for kinetic modeling.

Elaborated thermochemical treatment of HF, CO, N2, and H2O: Insight into HEAT and its extensions.

Empirical, highly accurate non-relativistic electronic total atomization energies (eTAEs) are established by combining experimental or computationally converged treatments of the nuclear motion and relativistic contributions with the total atomization energies of HF, CO, N2, and H2O obtained from the Active Thermochemical Tables. These eTAEs, which have estimated (2σ) uncertainties of less than 10 cm-1 (0.12 kJ mol-1), form the basis for an analysis of high-level ab initio quantum chemical calculations that aim at reproducing these eTAEs for the title molecules. The results are then employed to analyze the performance of the high-accuracy extrapolated ab initio thermochemistry, or High-Accuracy Extrapolated Ab Initio Thermochemistry (HEAT), family of theoretical methods. The method known as HEAT-345(Q), in particular, is found to benefit from fortuitous error cancellation between its treatment of the zero-point energy, extrapolation errors in the Hartree-Fock and coupled cluster contributions, neglect of post-(T) core-correlation, and the basis-set error involved in higher-level correlation corrections. In addition to shedding light on a longstanding curiosity of the HEAT protocol-where the cheapest HEAT-345(Q) performs comparably to the theoretically more complete HEAT-456QP procedure-this study lays the foundation for extended HEAT variants that offer substantial improvements in accuracy relative to the established approaches.

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 (CH3C(O)CH3) at a nominal temperature of 1800 K. The major product observed is ketene (CH2CO). Minor products identified include acetaldehyde (CH3CHO), propyne (CH3CCH), propene (CH2CHCH3), 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. 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.

An Automated Thermochemistry Protocol Based on Explicitly Correlated Coupled-Cluster Theory: The Methyl and Ethyl Peroxy Families.

An automated computational thermochemistry protocol based on explicitly correlated coupled-cluster theory was designed to produce highly accurate enthalpies of formation and atomization energies for small- to medium-sized molecular species (3-12 atoms). Each potential source of error was carefully examined, and the sizes of contributions to the total atomization enthalpies were used to generate uncertainty estimates. The protocol was first used to generate total atomization enthalpies for a family of four molecular species exhibiting a variety of charges, multiplicities, and electronic ground states. The new protocol was shown to be in good agreement with the Active Thermochemical Tables database for the four species: the methyl peroxy radical, methoxyoxoniumylidene (methyl peroxy cation), methyl peroxy anion, and methyl hydroperoxide. Updating the Active Thermochemical Tables to include those results yielded significantly improved accuracy for the formation enthalpies of those species. The derived protocol was then used to predict formation enthalpies for the larger ethyl peroxy family of species.

High-accuracy extrapolated ab initio thermochemistry. IV. A modified recipe for computational efficiency.

A number of economical modifications to the high-accuracy extrapolated ab initio thermochemistry (HEAT) model chemistry are evaluated. The two resulting schemes, designated as mHEAT and mHEAT+, are designed for efficient and pragmatic evaluation of molecular energies in systems somewhat larger than can be practically studied by the unapproximated HEAT scheme. It is found that mHEAT+ produces heats of formation with nearly subchemical (±1 kJ/mol) accuracy at a substantially reduced cost relative to the full scheme. Total atomization energies calculated using the new thermochemical recipes are compared to the results of the HEAT-345(Q) model chemistry, and enthalpies of formation for the three protocols are also compared to Active Thermochemical Tables. Finally, a small selection of transition states is studied using mHEAT and mHEAT+, which illuminates some interesting features of reaction barriers and serves as an initial benchmark of the performance of these model chemistries for chemical kinetics applications.

Active Thermochemical Tables: The Partition Function of Hydroxymethyl (CH2OH) Revisited.

The best currently available set of temperature-dependent nonrigid rotor anharmonic oscillator (NRRAO) thermochemical and thermophysical properties of hydroxymethyl radical is presented. The underlying partition function relies on a critically evaluated complement of accurate experimental and theoretical data and is constructed using a two-pronged strategy that combines contributions from large amplitude motions obtained from direct counts, with contributions from the other internal modes of motion obtained from analytic NRRAO expressions. The contributions from the two strongly coupled large-amplitude motions of CH2OH, OH torsion and CH2 wag, are based on energy levels obtained by solving the appropriate two-dimensional projection of a fully dimensional potential energy surface that was recently obtained at the CCSD(T)/cc-pVTZ level of theory. The contributions of the remaining seven, more rigid, vibrational modes and of the external rotations are captured by NRRAO corrections to the standard rigid rotor harmonic oscillator (RRHO) treatment, which include corrections for vibrational anharmonicities, rotation-vibration interaction, Coriolis effects, and low temperature. The basic spectroscopic constants needed for the construction of the initial RRHO partition function rely on experimental ground-state rotational constants and the best available experimental fundamentals, additionally complemented by fundamentals obtained from the variational solution of the full-dimensional potential energy surface using a recently developed two-component multilayer Lanczos algorithm. The higher-order spectroscopic constants necessary for the NRRAO corrections are extracted from a second-order variational perturbation treatment (VPT2) of the same potential energy surface. The Lanczos solutions of the fully dimensional surface are validated against available experimental data, and the VPT2 results and the solutions of the reduced dimensionality surface are validated both against the Lanczos solutions and available experiments. The NRRAO thermophysical and thermochemical properties, given both in tabular form and as seven- and nine-coefficient NASA polynomials, are compared to previous results. In addition, the latest ATcT values for the enthalpy of formation of CH2OH at 298.15 K (0 K), -16.75 ± 0.27 kJ/mol (-10.45 ± 0.27 kJ/mol), and of other related CH nO m species ( n = 0-4, m = 0,1) are reported, together with a plethora of related bond dissociation enthalpies (BDEs), such as the C-H, O-H, and C-O bond dissociation enthalpies of methanol, 402.16 ± 0.26 kJ/mol (395.61 ± 0.26 kJ/mol), 440.34 ± 0.26 kJ/mol (434.86 ± 0.26 kJ/mol), and 384.85 ± 0.15 kJ/mol (377.14 ± 0.15 kJ/mol), respectively, and analogous BDEs for hydroxymethyl, 343.67 ± 0.37 kJ/mol (339.16 ± 0.37 kJ/mol), 125.54 ± 0.28 kJ/mol (121.11 ± 0.28 kJ/mol), and 445.86 ± 0.29 kJ/mol (438.76 ± 0.29 kJ/mol), respectively. The reasons governing the alternation between strong and weak sequential H atom BDEs of methanol are also discussed.

A master equation simulation for the •OH + CH3OH reaction.

A combined (fixed-J) two-dimensional master-equation/semi-classical transition state theory/variational Rice-Ramsperger-Kassel-Marcus approach has been used to compute reaction rate coefficients of •OH with CH3OH over a wide range of temperatures (10-2500 K) and pressures (10-1-104 Torr) based on a potential energy surface that has been constructed using a modification of the high accuracy extrapolated ab initio thermochemistry (HEAT) protocol. The calculated results show that the title reaction is nearly pressure-independent when T > 250 K but depends strongly on pressure at lower temperatures. In addition, the preferred mechanism and rate constants are found to be very sensitive to temperature. The reaction pathway CH3OH + •OH → CH3O• + H2O proceeds exclusively through tunneling at exceedingly low temperatures (T ≤ 50 K), typical of those established in interstellar environments. In this regime, the rate constant is found to increase with decreasing temperature, which agrees with low-temperature experimental results. The thermodynamically favored reaction pathway CH3OH + •OH → •CH2OH + H2O becomes dominant at higher temperatures (T ≥ 200 K), such as those found in Earth's atmosphere as well as combustion environments. By adjusting the ab initio barrier heights slightly, experimental rate constants from 200 to 1250 K can be satisfactorily reproduced.

Enthalpy of Formation of C2H2O4 (Oxalic Acid) from High-Level Calculations and the Active Thermochemical Tables Approach.

High-level coupled cluster calculations obtained with the Feller-Peterson-Dixon (FPD) approach and new data from the most recent version of the Active Thermochemical Tables (ATcT) are used to reassess the enthalpy of formation of gas-phase C2H2O4 (oxalic acid). The theoretical value was further calibrated by comparing FPD and ATcT gas-phase enthalpies of formation for H2CO (formaldehyde) and the two low-lying conformations of C2H4O2 ( syn and anti acetic acid). The FPD approach produces a theoretical enthalpy of formation of gas-phase oxalic acid of -732.2 ± 4.0 kJ/mol at 298.15 K (-721.8 ± 4.0 kJ/mol at 0 K). An independently obtained ATcT value, based on reassessing the existent experimental determinations and expanding the resulting thermochemical network with select mid-level composite theoretical results, disagrees with several earlier recommendations that were based solely on experimental determinations but is in excellent accord with the current FPD value. The inclusion of the latter in the most recent ATcT thermochemical network produces a further refined value for the gas-phase enthalpy of formation, -731.6 ± 1.2 kJ/mol at 298.15 K (-721.0 ± 1.2 kJ/mol at 0 K). The condensed-phase ATcT enthalpy of formation of oxalic acid is -829.7 ± 0.5 kJ/mol, and the resulting sublimation enthalpy is 98.1 ± 1.3 kJ/mol, both at 298.15 K.

Unimolecular Reaction of Methyl Isocyanide to Acetonitrile: A High-Level Theoretical Study.

A combination of high-level coupled-cluster calculations and two-dimensional master equation approaches based on semiclassical transition state theory is used to reinvestigate the classic prototype unimolecular isomerization of methyl isocyanide (CH3NC) to acetonitrile (CH3CN). The activation energy, reaction enthalpy, and fundamental vibrational frequencies calculated from first-principles agree well with experimental results. In addition, the calculated thermal rate constants adequately reproduce those of experiment over a large range of temperature and pressure in the falloff region, where experimental results are available, and are generally consistent with statistical chemical kinetics theory (such as Rice-Ramsperger-Kassel-Marcus (RRKM) and transition state theory (TST)).

Time-Resolved Kinetic Chirped-Pulse Rotational Spectroscopy in a Room-Temperature Flow Reactor.

Chirped-pulse Fourier transform millimeter-wave spectroscopy is a potentially powerful tool for studying chemical reaction dynamics and kinetics. Branching ratios of multiple reaction products and intermediates can be measured with unprecedented chemical specificity; molecular isomers, conformers, and vibrational states have distinct rotational spectra. Here we demonstrate chirped-pulse spectroscopy of vinyl cyanide photoproducts in a flow tube reactor at ambient temperature of 295 K and pressures of 1-10 µbar. This in situ and time-resolved experiment illustrates the utility of this novel approach to investigating chemical reaction dynamics and kinetics. Following 193 nm photodissociation of CH2CHCN, we observe rotational relaxation of energized HCN, HNC, and HCCCN photoproducts with 10 µs time resolution and sample the vibrational population distribution of HCCCN. The experimental branching ratio HCN/HCCCN is compared with a model based on RRKM theory using high-level ab initio calculations, which were in turn validated by comparisons to Active Thermochemical Tables enthalpies.

Post-transition state dynamics and product energy partitioning following thermal excitation of the F⋯HCH2CN transition state: Disagreement with experiment.

Born-Oppenheimer direct dynamics simulations were performed to study atomistic details of the F + CH3CN → HF + CH2CN H-atom abstraction reaction. The simulation trajectories were calculated with a combined M06-2X/MP2 algorithm utilizing the 6-311++G** basis set. The experiments were performed at 300 K, and assuming the accuracy of transition state theory (TST), the trajectories were initiated at the F⋯HCH2CN abstraction TS with a 300 K Boltzmann distribution of energy and directed towards products. Recrossing of the TS was negligible, confirming the accuracy of TST. HF formation was rapid, occurring within 0.014 ps of the trajectory initiation. The intrinsic reaction coordinate (IRC) for reaction involves rotation of HF about CH2CN and then trapping in the CH2CN⋯HF post-reaction potential energy well of ∼10 kcal/mol with respect to the HF + CH2CN products. In contrast to this IRC, five different trajectory types were observed: the majority proceeded by direct H-atom transfer and only 11% approximately following the IRC. The HF vibrational and rotational quantum numbers, n and J, were calculated when HF was initially formed and they increase as potential energy is released in forming the HF + CH2CN products. The population of the HF product vibrational states is only in qualitative agreement with experiment, with the simulations showing depressed and enhanced populations of the n = 1 and 2 states as compared to experiment. Simulations with an anharmonic zero-point energy constraint gave product distributions for relative translation, HF rotation, HF vibration, CH2CN rotation, and CH2CN vibration as 5%, 11%, 60%, 7%, and 16%, respectively. In contrast, the experimental energy partitioning percentages to HF rotation and vibration are 6% and 41%. Comparisons are made between the current simulation and those for other F + H-atom abstraction reactions. The simulation product energy partitioning and HF vibrational population for F + CH3CN → HF + CH2CN resemble those for other reactions. A detailed discussion is given of possible origins of the difference between the simulation and experimental energy partitioning dynamics for F + CH3CN → HF + CH2CN. The F + CH3CN reaction also forms the CH3C(F)N intermediate, in which the F-atom adds to the C≡N bond. However, this intermediate and F⋯CH3CN and CH3CN⋯F van der Waals complexes are not expected to affect the F + CH3CN → HF + CH2CN product energy partitioning.

Active Thermochemical Tables: The Adiabatic Ionization Energy of Hydrogen Peroxide.

The adiabatic ionization energy of hydrogen peroxide (HOOH) is investigated, both by means of theoretical calculations and theoretically assisted reanalysis of previous experimental data. Values obtained by three different approaches: 10.638 ± 0.012 eV (purely theoretical determination), 10.649 ± 0.005 eV (reanalysis of photoelectron spectrum), and 10.645 ± 0.010 eV (reanalysis of photoionization spectrum) are in excellent mutual agreement. Further refinement of the latter two values to account for asymmetry of the rotational profile of the photoionization origin band leads to a reduction of 0.007 ± 0.006 eV, which tends to bring them into even closer alignment with the purely theoretical value. Detailed analysis of this fundamental quantity by the Active Thermochemical Tables approach, using the present results and extant literature, gives a final estimate of 10.641 ± 0.006 eV.

Ab Initio Computations and Active Thermochemical Tables Hand in Hand: Heats of Formation of Core Combustion Species.

The fidelity of combustion simulations is strongly dependent on the accuracy of the underlying thermochemical properties for the core combustion species that arise as intermediates and products in the chemical conversion of most fuels. High level theoretical evaluations are coupled with a wide-ranging implementation of the Active Thermochemical Tables (ATcT) approach to obtain well-validated high fidelity predictions for the 0 K heat of formation for a large set of core combustion species. In particular, high level ab initio electronic structure based predictions are obtained for a set of 348 C, N, O, and H containing species, which corresponds to essentially all core combustion species with 34 or fewer electrons. The theoretical analyses incorporate various high level corrections to base CCSD(T)/cc-pVnZ analyses (n = T or Q) using H2, CH4, H2O, and NH3 as references. Corrections for the complete-basis-set limit, higher-order excitations, anharmonic zero-point energy, core-valence, relativistic, and diagonal Born-Oppenheimer effects are ordered in decreasing importance. Independent ATcT values are presented for a subset of 150 species. The accuracy of the theoretical predictions is explored through (i) examination of the magnitude of the various corrections, (ii) comparisons with other high level calculations, and (iii) through comparison with the ATcT values. The estimated 2σ uncertainties of the three methods devised here, ANL0, ANL0-F12, and ANL1, are in the range of ±1.0-1.5 kJ/mol for single-reference and moderately multireference species, for which the calculated higher order excitations are 5 kJ/mol or less. In addition to providing valuable references for combustion simulations, the subsequent inclusion of the current theoretical results into the ATcT thermochemical network is expected to significantly improve the thermochemical knowledge base for less-well studied species.

Enthalpy of Formation of N2H4 (Hydrazine) Revisited.

In order to address the accuracy of the long-standing experimental enthalpy of formation of gas-phase hydrazine, fully confirmed in earlier versions of Active Thermochemical Tables (ATcT), the provenance of that value is re-examined in light of new high-end calculations of the Feller-Peterson-Dixon (FPD) variety. An overly optimistic determination of the vaporization enthalpy of hydrazine, which created an unrealistically strong connection between the gas phase thermochemistry and the calorimetric results defining the thermochemistry of liquid hydrazine, was identified as the probable culprit. The new enthalpy of formation of gas-phase hydrazine, based on balancing all available knowledge, was determined to be 111.57 ± 0.47 kJ/mol at 0 K (97.42 ± 0.47 kJ/mol at 298.15 K). Close agreement was found between the ATcT (even excluding the latest theoretical result) and the FPD enthalpy.

Thermal Decomposition of Potential Ester Biofuels. Part I: Methyl Acetate and Methyl Butanoate.

Two methyl esters were examined as models for the pyrolysis of biofuels. Dilute samples (0.06-0.13%) of methyl acetate (CH3COOCH3) and methyl butanoate (CH3CH2CH2COOCH3) were entrained in (He, Ar) carrier gas and decomposed in a set of flash-pyrolysis microreactors. The pyrolysis products resulting from the methyl esters were detected and identified by vacuum ultraviolet photoionization mass spectrometry. Complementary product identification was provided by matrix infrared absorption spectroscopy. Pyrolysis pressures in the pulsed microreactor were about 20 Torr and residence times through the reactors were roughly 25-150 µs. Reactor temperatures of 300-1600 K were explored. Decomposition of CH3COOCH3 commences at 1000 K, and the initial products are (CH2═C═O and CH3OH). As the microreactor is heated to 1300 K, a mixture of CH2═C═O and CH3OH, CH3, CH2═O, H, CO, and CO2 appears. The thermal cracking of CH3CH2CH2COOCH3 begins at 800 K with the formation of CH3CH2CH═C═O and CH3OH. By 1300 K, the pyrolysis of methyl butanoate yields a complex mixture of CH3CH2CH═C═O, CH3OH, CH3, CH2═O, CO, CO2, CH3CH═CH2, CH2CHCH2, CH2═C═CH2, HCCCH2, CH2═C═C═O, CH2═CH2, HC≡CH, and CH2═C═O. On the basis of the results from the thermal cracking of methyl acetate and methyl butanoate, we predict several important decomposition channels for the pyrolysis of fatty acid methyl esters, R-CH2-COOCH3. The lowest-energy fragmentation will be a 4-center elimination of methanol to form the ketene RCH═C═O. At higher temperatures, concerted fragmentation to radicals will ensue to produce a mixture of species: (RCH2 + CO2 + CH3) and (RCH2 + CO + CH2═O + H). Thermal cracking of the ß C-C bond of the methyl ester will generate the radicals (R and H) as well as CH2═C═O + CH2═O. The thermochemistry of methyl acetate and its fragmentation products were obtained via the Active Thermochemical Tables (ATcT) approach, resulting in ΔfH298(CH3COOCH3) = -98.7 ± 0.2 kcal mol-1, ΔfH298(CH3CO2) = -45.7 ± 0.3 kcal mol-1, and ΔfH298(COOCH3) = -38.3 ± 0.4 kcal mol-1.