Search

AlO is relatively abundant around oxygen-rich Asymptotic Giant Branch (AGB) stars, where it can react with major gas-phase species such as H2 and H2O to form AlOH. These Al-containing species are t...

ConspectusHomogeneous metal-mediated organometallic reactions represent a very large and diverse reaction class. Density functional theory calculations are now routinely carried out and reported for analyzing organometallic mechanisms and reaction pathways. While density functional theory calculations are extremely powerful to understand the energy and structure of organometallic reactions, there are several assumptions in their use and interpretation to define reaction mechanisms and to analyze reaction selectivity. Almost always it is assumed that potential energy structures calculated with density functional theory adequately describe mechanisms and selectivity within the framework of statistical theories, for example, transition state theory and RRKM theory. However, these static structures and corresponding energy landscapes do not provide atomic motion information during reactions that could reveal nonstatistical intermediates without complete intramolecular vibrational redistribution and nonintrinsic reaction coordinate (non-IRC) pathways. While nonstatistical intermediates and non-IRC reaction pathways are now relatively well established for organic reactions, these dynamic effects have heretofore been highly underexplored in organometallic reactions. Through a series of quasiclassical density functional theory direct dynamics trajectory studies, my group has recently demonstrated that dynamic effects occur in a variety of fundamental organometallic reactions, especially bond activation reactions. For example, in the C-H activation reaction between methane and [Cp*(PMe3)IrIII(CH3)]+, while the density functional theory energy landscape showed a two-step oxidative cleavage and reductive coupling mechanism, trajectories revealed a mixture of this two-step mechanism and a dynamic one-step mechanism that skipped the [Cp*(PMe3)IrV(H)(CH3)2]+ intermediate. This study also showed that despite a methane σ-complex being located on the density functional theory surface before oxidative cleavage and after reductive coupling, this intermediate is always skipped and should not be considered an intermediate during reactive trajectories. For non-IRC reaction pathways, quasiclassical direct dynamics trajectories showed that for the isomerization of [Tp(NO)(PMe3)W(η2-benzene)] to [Tp(NO)(PMe3)W(H)(Ph)], there are many dynamic reaction pathway connections due to a relatively flat energy landscape and π coordination is not necessary for C-H bond activation through oxidative cleavage. Trajectories also showed that dynamic effects are important in selectivity for ethylene C-H activation versus π coordination in reaction with Cp(PMe3)2Re, and trajectories provide a more quantitative model of selectivity than transition state theory. Quasiclassical trajectories examining Au-catalyzed monoallylic diol cyclizations showed dynamic coupling of several reaction steps that include alkoxylation π bond addition, proton shuttling, and water elimination reaction steps. Overall, these studies highlight the need to use direct dynamics trajectory simulations to consider atomic motion during reactions to understand organometallic reaction mechanisms and selectivity.

The oxidation mechanism of CH3CHCO initiated by OH is systematically studied using quantum theoretical methods. According to thermodynamic research, the dominant channel is the addition of CH3CHCO with OH generating C-IM1 (CH3CHOHCO) and C-IM2 (CH3CHCOOH), and then dissociate to the dominant products P1 (CH3CHO + CO) from C-IM1 with the lowest barrier. The rate constants for CH3CHCO + OH → C-IM1 → P1/CH3CHCO + OH → C-IM2 channels are computed through RRKM theory at 200–2000 K. The obtained overall rate constant at 298 K (7.60 × 10−11 cm3 molecule−1 s−1) is well consistent with the reported experimental value. The atmospheric lifetime of CH3CHCO is estimated to be around 1.83 h.

Decomposition of heptafluoro-iso-butyronitrile (C 3 F 7 CN), a novel eco-friendly dielectric gas for high voltage electrical applications, is investigated theoretically in the presence of water vapor using various computational methods including DFT, CBS-Q, G4, CCSD(T)-F12 and RRKM theory. Hydrolysis of C 3 F 7 CN proceeds dominantly via the stepwise addition/isomerization or elimination mechanisms to form amide or HF. The reaction pathways are bifurcated into cis and trans conformations. Activation barriers are strongly dependent on the number of water molecules in hydrolysis due to catalytic effect. The intermediate imidic acid prefers to produce formamido imine via recombination or secondary reaction with C 3 F 7 CN. Hydrolysis of CF 3 and CF 3 CFCN undergoes via the direct H-abstraction mechanism to form OH radicals enhancing hydrolysis of C 3 F 7 CN. Kinetic calculations predict that imidic acid is the major decomposition product below 500 K while amide is generated exclusively above 1000 K with a small fraction of HF. The present work not only provides a theoretical basis for the increasing risk of dielectric failure of C 3 F 7 CN with the increment of the high humidity levels, but also identifies the potential characteristic decomposition products to develop failure diagnostic techniques.

The aim of this study is to calculate the unimolecular rate coefficients for the unimolecular decomposition reaction of the industrially important molecule, SF6. The energies of stationary-points involved in the title reaction are calculated by the combination W1 method. Two main reaction paths are considered: SF6 → SF5 + F (R1) and SF6 → SF4 + F2 (R2). Having information on energies and molecular properties of reactants and transition-states, RRKM statistical rate theory is used to compute the rate coefficients as a function of temperature and pressure. For the bond dissociation process R1, special version of RRKM theory, i. e., Variable reaction coordinate-transition state theory (VRC-TST) is employed. Although the reaction R1 is the dominant process over a wide range of pressure and temperature, but the reaction R2 could be significant at high temperatures. The following Arrhenius expressions are obtained for high-pressure limiting rate constants of reaction paths R1 and R2: k∞,1 = 5.71 × 1016 s−1 exp (−429.8 kJ mol−1 /RT) k∞,2 = 2.14 × 1016 s−1 exp (−590.6 kJ mol−1 /RT).

A detailed chemical kinetic model of aluminum oxidation in gas phase has been constructed. Quantum chemical calculations are performed for obtaining molecular structures of intermediates on the reaction pathways. DFT level calculation is applied to search intermediates and transition states, and CBS-QB3 method is applied to calculate highly accurate potential energies. The rate coefficients of each reaction paths were also calculated based on VTST or RRKM theory for constructing the detailed chemical kinetic model. 76 elementary reactions were investigated for AlxOy intermediates, and chemical kinetics analysis was performed to find the major route of oxidation to Al8O12 as the final product of this model. The mainstream of oxidation and formation of large size AlxOy molecules were discussed on the model, and it was found that AlO, Al2O, and Al2O2 are important intermediates to form Al2O3 or Al2O4 as precursors for Al4O6.

A single-pulse shock tube study of the four pentene isomers is carried out at 2¿±¿0.16¿bar and 900¿1600¿K. C1 to C6 species profiles were recorded using gas chromatography mass spectrometry analyses. The species are identified using mass spectrometry and quantified by flame ionization detection. High-pressure limiting and pressure-dependent rate constants for 2M1B, 2M2B and 3M1B¿+¿¿ were calculated using RRKM theory with a Master Equation (ME) analysis using the Master Equation System Solver, MESS. A mechanism was formulated based on rate rules and theoretical calculations. Comparisons between experimental results and model simulations are provided for all of the five pentene isomers investigated with satisfactory agreement. Furthermore, an insight is provided into the influence of molecular structure on the reactivity of pyrolysis chemistry. Interestingly, it is found that the HACA mechanism is much less prominent for benzene formation compared to the role of cyclopentadienyl radical recombination with methyl radicals and also the recombination of propargyl radicals. The authors would like to acknowledge Science Foundation Ireland for funding via project numbers 15/IA/3177 and 16/SP/3829. The work at LLNL was performed under the auspices of the U.S. Department of Energy (DOE), Contract DE-AC52-07NA27344, and was supported by the U.S. Department of Energy, Vehicle Technologies Office, program managers, Mike Weismiller and Gurpreet Singh. peer-reviewed

2-Methyl-3-buten-2-ol (MBO232) is a biogenic volatile organic compound (BVOC), and has a large percentage of emission into the atmosphere. The vacuum ultraviolet (VUV) photochemistry of BVOCs is of great importance for atmospheric chemistry. Studies have been carried out on several BVOCs but have not extended to MBO232. In the present report, the photoionization and dissociation processes of MBO232 in the energy range of 8.0-15.0 eV have been studied by tunable VUV synchrotron radiation coupled with a time-of-flight mass spectrometer. By measuring the photoionization spectra, the adiabatic ionization energy (AIE) of MBO232 and the appearance energies (AEs) of the eight identified fragment ions (i.e., C4H7O+, C3H7O+, C5H9+, C3H6O+, CH3CO+, CH3O+, C4H5+, and C3H5+) were determined. High-level quantum chemistry calculations suggest that there are 3 direct channels and 5 indirect channels via transition states and intermediates accountable for these fragments. Among the reaction channels, the direct elimination of CH3 is the most dominant channel and produces the resonance-stabilized radical cation. Most interestingly, our results show that the CH3 selectively migrates towards the cation, which leads to the different indirect channels. The CH3 migration is a rare process in the dissociative photoionization of metal-free organic molecules. We explain the process by molecular orbital calculations and electron localization function analysis and explore the non-conventional dissociation channels via the CH3 roaming mechanism. We further perform kinetics analysis using RRKM theory for the channels of interest. The activation barrier, and rate constants are analyzed for the branching fractions of the products. These results provide important implications for the VUV photochemistry of BVOCs in the atmosphere.

The temperature- and pressure-dependence of rate constants for several radicals and unsaturated hydrocarbons reactions (1,3-C5H8/1,4-C5H8/cyC(5)H(8) + (H)over dot, C2H4 + (C)over dot(3)H5-a, C3H6 + (C)over dot(2)H(3)) are analyzed in this paper. The abstraction reactions of these systems are also calculated and compared with available literature data. (C)over dot(5)H(9) radicals can be produced via (H)over dot atom addition reactions to the pentadiene isomers and cyclopentene, and also by H-atom abstraction reactions from 1- and 2-pentene and cyclopentane. Comprehensive (C)over dot(5)H(9) potential energy surface (PES) analyses and high-pressure limiting rate constants for related reactions have been explored in part I of this work (J. Phys. Chem. A 2019, 123 (22), 9019-9052). In this work, a chemical kinetic model is constructed based on the computed thermochemistry and high-pressure limiting rate constants from part I, to further understand the chemistry of different C5H8 molecules. The most important channels for these addition reactions are discussed in the present work based on reaction pathway analyses. The dominant reaction pathways for these five systems are combined together to generate a simplified (C)over dot(5)H(9) PES including nine reactants, 25 transition states (TSs), and nine products. Spin-restricted single point energies are calculated for radicals and TSs on the simplified PES at the ROCCSD(T)/aug-cc-pVTZ level of theory with basis set corrections from MP2/aug-cc-pVXZ (where X = T and Q). Temperature- and pressure-dependent rate constants are calculated using RRKM theory with a Master Equation analysis, with restricted energies used for minima on the simplified (C)over dot(5)H(9) PES and unrestricted energies for other species, over a temperature range of 300-2000 K and in the pressure range 0.01-100 atm. The rate constants calculated are in good agreement with those in the literature. The chemical kinetic model is updated with pressure-dependent rate constants and is used to simulate the species concentration profiles for H. atom addition to cyclopentane and cyclopentene. Through detailed analyses and comparisons, this model can reproduce the experimental measurements of species qualitatively and quantitatively with reasonably good agreement. This study is supported by Science Foundation Ireland and the China Scholarship Council (CSC). The authors want to acknowledge the financial support of Science Foundation Ireland under Grant No. 15/IA/3177 and 16/SP/3829, and the provision of computational resources from ICHEC under the NUI Galway shared condominium accounts. The Computational resources are provided by the Irish Centre for High-End Computing (ICHEC), under project number ngche063c, ngcom006c and ngche058c. The authors are grateful to Stephen Klippenstein for the help with MESS code study. Chong-Wen Zhou acknowledges the support from National Science and Technology Major Project (2017-III-0004-0028) and Beihang University under the Fundamental Research Funds. peer-reviewed 2021-05-12

The traditional understanding is that the Hinshelwood-Lindemann mechanism for thermal unimolecular reactions, and the resulting unimolecular rate constant versus temperature and collision frequency ω (i.e., pressure), requires the Rice-Ramsperger-Kassel-Marcus (RRKM) rate constant k(E) to represent the unimolecular reaction of the excited molecule versus energy. RRKM theory assumes an exponential N(t)/N(0) population for the excited molecule versus time, with decay given by RRKM microcanonical k(E), and agreement between experimental and Hinshelwood-Lindemann thermal kinetics is then deemed to identify the unimolecular reactant as a RRKM molecule. However, recent calculations of the Hinshelwood-Lindemann rate constant kuni(ω,E) has brought this assumption into question. It was found that a biexponential N(t)/N(0), for intrinsic non-RRKM dynamics, gives a Hinshelwood-Lindemann kuni(ω,E) curve very similar to that of RRKM theory, which assumes exponential dynamics. The RRKM kuni(ω,E) curve was brought into agreement with the biexponential kuni(ω,E) by multiplying ω in the RRKM expression for kuni(ω,E) by an energy transfer efficiency factor βc. Such scaling is often done in fitting Hinshelwood-Lindemann-RRKM thermal kinetics to experiment. This agreement between the RRKM and non-RRKM curves for kuni(ω,E) was only obtained previously by scaling and fitting. In the work presented here, it is shown that if ω in the RRKM kuni(ω,E) is scaled by a βc factor there is analytic agreement with the non-RRKM kuni(ω,E). The expression for the value of βc is derived.

The reaction between CFCl2CH2O2 radicals and ClO was studied using the B3LYP and CCSD(T) methods associated with the 6-311++G(d,p) and cc-pVTZ basis sets, and subsequently RRKM-TST theory was used to predict the thermal rate constants and product distributions. On the singlet PES, the dominant reaction is the addition of the ClO oxygen atom to the terminal-O of CFCl2CH2O2 to generate adduct IM1 (CFCl2CH2OOOCl), and then dissociation to final products P1 (CFCl2CHO + HO2 + Cl) occurs. RRKM theory is employed to calculate the overall and individual rate constants over a wide range of temperatures and pressures. It is predicted that the collision-stabilized IM1 (CFCl2CH2OOOCl) dominates the reaction at 200–500 K (accounting for about 60–100%) and the dominant products are P1 (CFCl2CHO + HO2 + Cl). The yields of the other products are very low and insignificant for the title reaction. The total rate constants exhibit typical “falloff” behavior. The pathways on the triplet PES are less competitive than that on the singlet PES. The calculated overall rate constants are in good agreement with the experimental data. The atmospheric lifetime of CFCl2CH2O2 in ClO is around 2.04 h. TD-DFT calculations imply that IM1 (CFCl2CH2OOOCl), IM2 (CFCl2CH2OOClO) and IM3 (CFCl2CH2OClO2) will photolyze under sunlight.

The kinetic and electron density flows are studied theoretically for the gas phase pyrolysis of cubane via its cage opening to reach bicyclooctatriene and then thermal rearrangement of bicyclooctatriene to produce [8]annulene which is the experimentally observed major product. The observed kinetic data at the MN15-L/maug-cc-pVTZ level of theory were in good agreement with the experimental results as compared to the CBS-QB3 method. The cage opening and the thermal rearrangement steps at the experimentally employed temperature of 520 K were exergonic and exothermic. The atmospheric rate constants calculated by means of the RRKM theory show that the cage opening is the rate-determining step. The temperature dependence of the rate constant for the cage opening step at the MN15-L level can be expressed as log(k/s−1)1barMN15-L = (15.63) − (48.99 kcal mol−1)/RT ln 10. The molecular mechanism of the reactions has been investigated by means of the bonding evolution theory (BET) at the B3LYP/6-311G (d,p) level of theory. The cage opening course is described topologically by cleaving of C1–C2, C4–C8, and C5–C6 single bonds and electron saturation of the C1–C4, C2–C6, and C5–C8 bonds, while the rearrangement of bicyclooctatriene is described by C3–C7 bond rupture, depopulation of C1–C4 and C5–C8 double bonds, and electron saturation of C1–C5, C3–C4, and C7–C8 bonds. Electron density rearrangement along the two successive steps are asynchronous and the sequence of catastrophes can be represented as: η-1-13-C†C†FFFC†C†FFFC†C†-2-6-[C]2C†[F]2[C†]2C†-0.

Criegee intermediates are reactive species formed in the ozonolysis of alkenes. Their subsequent chemistry is critical to an accounting of OH production, aerosol formation, and the oxidative capaci...

Tandem mass spectrometry was used to explore the trends in the unimolecular fragmentation of the ionized hydroxy-substituted polycyclic aromatic hydrocarbons 1-naphthol, 9-hydroxyphenanthrene, and 1-hydroxypyrene. The main dissociation reactions across all ring systems were CO- and HCO-losses, with ionized 1-naphthol also losing H atoms. Both ionized 1-naphthol and 9-hydroxyphenanthrene displayed the sequential loss of C2H2 and C4H2 from the [M-HCO]+ ions, reminiscent of unsubstituted PAH ions. CO-loss is slightly favored for 1-naphthol and 9-hydroxyphenanthrene, at low collision energy, but less so for 1-hydroxypyrene. Reaction mechanisms for HCO- and CO-losses from 1-hydroxypyrene were derived from CCSD/6-31G(d)//B3-LYP/6-31G(d) calculations. The CO-loss mechanism is dominated by two transition states: TS-A governing a 1,3-H shift in the molecular ion and TS-C which governs a ring-closing step to form a five-member ring in the product ion. HCO-loss occurs over a much flatter potential energy surface with the intermediate being the product ion bound to the carbon atom of HCO. Imaging photoelectron photoion coincidence spectroscopy of 1-hydroxypyrene yielded threshold photon-energy resolved breakdown curves and time-of-flight distributions that were modeled with RRKM theory to give 0 K reaction energies for HCO- and CO-losses of 3.92 ± 0.05 and 2.91 ± 0.05 eV, respectively. The entropies of activation for the two channels were very different, 14 and 95 JK-1 mol-1, respectively, a result consistent with the calculated mechanisms. The threshold photoelectron spectrum yielded an IE value of 7.14 ± 0.01 eV for 1-hydroxypyrene.

Oxygenated polycyclic aromatic hydrocarbons (OPAH) have received increasing attention due to their toxic effect on human health. This study comprehensively investigates the evolution of OPAH chemistry at flame temperatures. Jet-stirred reactor (JSR) experiments with benzene/phenol/C2H2/N2 and benzene/C2H2/O2/N2 revealed that OPAH with oxygenated heterocycle can be formed by the addition of C2H2 at 1400 K. To further clarify the evolution of OPAH chemistry in soot systems, OPAH formation and decomposition reaction pathways and kinetic parameters have been theoretically investigated. The potential energy surfaces of 1-naphtholate and 2-naphtholate growth, and thermal decomposition reactions, were calculated by combining the density functional theory B3LYP/6–311+G(d,p) and CCSD(T)/cc-pvdz methods. The reaction rate coefficients in the temperature range of 800–2500 K and pressure range of 0.1–100 atm were calculated using RRKM theory by solving the master equations. The potential energy surface of C2H2+1-naphtholate and C2H2+2-naphtholate growth reactions showed that the O atom could be locked in a naphthofuran molecule with the formation of a C O C oxygenated heterocycle; and the reaction rates were determined by adding the C2H2 elementary step with the energy barrier of 26.0 and 19.9 kcal/mol, respectively. Thermal decomposition reactions of 1-naphtholate and 2-naphtholate yielded an indenyl radical and CO. The thermal decomposition reaction rates were significantly sensitive to the zig-zag site structure next to the C O bond. The decomposition rate of 1-naphtholate at 1500 K, with a zig-zag site near the C O bond, was 14.8 times lower than that of 2-naphtholate with no zig-zag site near the C O bond. Rate comparison results indicate that the C O functional group rapidly converts to a C O C functional group with the addition of C2H2. The formation, growth and thermal decomposition reactions of 1-naphtholate and 2-naphtholate were added to a detailed PAH mechanism to check the effect of OPAH reactions on PAH formation chemistry. The concentration profile of naphthalene predicted by the updated PAH mechanism was lower than current PAH mechanism predictions by 29%, indicating that the OPAH reactions had a significant effect on PAH formation chemistry, and should be included in the PAH mechanism. However, due to the relatively low concentration of OPAH compared to PAH, it is possible to ignore the correlation between OPAH and soot nucleation at flame temperatures; therefore an OPAH evolution pathway (PAH → incipient soot → OPAH formation on soot particle → selective thermal decomposition of OPAH), is proposed to explain the high content of OPAH molecules (e.g., 9,10-anthraquinone, benz(a)anthracene-7,12-dione, and benzanthrone) adsorbed on the soot particle.

A comprehensive experimental and kinetic study of 2,4,4-trimethyl-1-pentene oxidation was conducted in this work. Oxidation in JSR was investigated at the temperature range of 675–1200 K and equivalence ratios of 0.5, 1.0, 2.0 under atmospheric pressure. The high-level quantum calculation was conducted on DLPNO-CCSD(T)/cc-pVTZ//M06-2X/6-311G(d,p) level to get some rate coefficients of high accuracy and choose important pathway. Temperature and pressure-dependent rate coefficients were obtained by solving master equation and RRKM theory. A detailed chemical kinetic model for 2,4,4-trimethyl-1-pentene was developed based on the high-level quantum calculation and validated against the mole fraction data measured in JSR in this work. Available experimental data including ignition delay times and laminar flame speeds in our previous work were also used for the validation of this model. Reaction pathway analysis and sensitivity analysis were performed using the Modified model. The results demonstrate the significance of H abstraction reactions in a wide range of temperature, while fuel unimolecular decomposition reaction are only dominant at high temperature. In addition, at low temperature, O 2 addition to octyl radical reactions are more important than radical decomposition and isomerization reactions. Furthermore, sensitivity analysis indicates that H abstraction reactions of fuel and isobutene have a strong inhibiting effect on low temperature reactivity.

The growth of polycyclic aromatic hydrocarbons (PAH) can proceed via multiple chemical mechanisms. The mechanism of naphthyl radical and vinylacetylene (C4H4) addition reaction has been systematically investigated in this computational study. A combination of DFT/B3LYP/6-311+G(d,p), CCSD/6-311+G(d,p) and CBS-QB3 methods were performed to calculate the potential energy surfaces. It revealed that the products, including phenanthrene, anthracene, a PAH with a five-membered ring structure, and PAH with a C4H3 radical substitution, can be formed in A2-1 (1-naphthyl)+C4H4 and A2-2 (2-naphthyl) +C4H4 reaction networks. The reaction rate constants at 0.1-100 atm were evaluated by RRKM theory by solving the master equation in the temperature range of 800–2500 K, which showed that the rate constants of reactions A2-1 (A2-2)+C4H4→product+H are highly temperature-dependent but nearly pressure-independent. The distribution of products was investigated in a 0-D batch reactor, wherein the initial reactant concentrations were taken from experimental measurements. The results showed that adduct intermediates were the main products at low temperature (T 1000 K). It was observed that a significant amount of phenanthrene is formed from PAH with a C4H3 radical substitution with the assistance of H atom. Reaction pathway sensitivity analysis for the PAH radical+C4H4 reaction system was performed and showed that the new benzene rings are more likely to be generated near the zig-zag edge surface site instead of the free edge. For the development of a PAH mechanism, the analogous treatment of rate constants for larger PAH radical + C4H4 reaction system are discussed. The formation rate of naphthalene from the reaction of phenyl+C4H4 was found to be very close to that of phenanthrene from the reaction of naphthyl+C4H4, suggesting that the analogous treatment of the rates is reasonable in PAH mechanisms.

The Journal of Physical Chemistry A, 125 (9), ISSN:1089-5639, ISSN:1520-5215

We have studied the decomposition and isomerization of furfural in the gas phase using quantum chemical and statistical reaction rate theory techniques. This work uncovers a variety of new reaction channels in furfural pyrolysis that lead to formation of the experimentally observed products, including CO2, which was previously unexplained. In addition to the known mechanism for furan + CO production, furfural is shown to isomerize directly to 2-pyrone, with a barrier height of 69 kcal mol-1, from where it can decompose to vinylketene + CO (highest barrier of 65 kcal mol-1) or to CO2 + 1,3-cyclobutadiene (highest barrier of 66 kcal mol-1). Alternative pathways to vinylketene + CO and 4-pyrone are also described. An RRKM theory/master equation model is developed to describe reactions on the C5O2H4 surface and used to simulate the decomposition kinetics of furfural and 2-pyrone. For both molecules, decomposition at 1400-2100 K is dominated by the formation of furan + CO, which represents around 75% of the total products, compared to around 19% and 6% for vinylketene + CO and total CO2, respectively. The model also predicts significant formation of stabilized 2-pyrone under these conditions. Rate coefficient expressions are reported as a function of both temperature and pressure for the main decomposition and isomerization channels identified in the pyrolysis of furfural and 2-pyrone, to facilitate detailed chemical kinetic modelling of these important oxygenated hydrocarbons.

In this work we report a joint experimental and computational study on the 2,5-dimethylfuran oxidation reaction in the gas phase initiated by atomic oxygen O(3P). The experiments have been performed by using vacuum-ultraviolet synchrotron radiation at the Advanced Light Source (ALS) of the Lawrence Berkeley National Laboratory (LBNL), at a temperature of 550 K and a pressure of 8 Torr. The experimental data were supported by quantum-chemical calculations along with a kinetic model, also taking into account the possible involvement of different magnetic states, performed in the framework of the RRKM theory. Propyne, acetaldehyde, methylglyoxal, dimethylglyoxal, 3-penten-2-one, 2,5-dimethylfuran-3(2H)-one, and 1,2-diacetyl ethylene have been identified as the main primary products arising under the conditions of the experiment. Our computational model suggests that these species can be formed at the concentration and branching ratio experimentally observed only in the presence of a non-negligible fraction of non-thermalized intermediates.