Thermochemical Properties, Pathway, and Kinetic Analysis on the Reactions of Benzene with OH: An Elementary Reaction Mechanism

On-The-Fly Kinetics of the Hydrogen Abstraction by Hydroperoxyl Radical: An Application of the Reaction Class Transition State Theory

A Reaction Class Transition State Theory (RC-TST) is applied to calculate thermal rate constants for hydrogen abstraction by OOH radical from alkanes in the temperature range of 300–2500 K. The rate constants for the reference reaction C₂H₆ + ∙OOH → ∙C₂H₅ + H₂O₂, is obtained with the Canonical Variational Transition State Theory (CVT) augmented with the Small Curvature Tunneling (SCT) correction. The necessary parameters were obtained from M06-2X/aug-cc-pVTZ data for a training set of 24 reactions. Depending on the approximation employed, only the reaction energy or no additional parameters are needed to predict the RC-TST rates for other class representatives. Although each of the reactions can in principle be investigated at higher levels of theory, the approach provides a nearly equally reliable rate constant at a fraction of the cost needed for larger and higher level calculations. The systematic error is smaller than 50% in comparison with high level computations. Satisfactory agreement with literature data, augmented by the lack of necessity of tedious and time consuming transition state calculations, facilitated the seamless application of the proposed methodology to the Automated Reaction Mechanism Generators (ARMGs) programs.

Benchmark calculations for bond dissociation energies and enthalpy of formation of chlorinated and brominated polycyclic aromatic hydrocarbons

Thermodynamic properties, i.e., bond dissociation energies and enthalpy of formation, of chlorinated and brominated polycyclic aromatic hydrocarbons play a fundamental role in understanding their formation mechanisms and reactivity. Computational electronic structure calculations routinely used to predict thermodynamic properties of various species are limited for these compounds due to large computational cost to obtain accurate results by employing high-level wave function theory methods. In this work, a number of composite model chemistry methods (CBS-QB3, G3MP2, G3, and G4) are used to compute bond dissociation energies and enthalpies of formation of small to medium-size chlorinated and brominated polycyclic aromatic hydrocarbon compounds. The enthalpy of formation is derived via the atomization method and compared against the recommended values. Statistical analysis indicates that G4 is the best method. For comparison, three commonly used density functional theory (DFT) methods (M06-2X, ωB97X-D and B2PLYP-D3) with various basis sets including 6-311++G(d, p), cc-pVTZ, and cc-pVQZ in the prediction of bond dissociation energies and enthalpies of formation have been tested using the optimized geometries at the same M06-2X/6-311++G(d, p) level of theory. It is found that ωB97X-D/6-311++G(d, p) shows the best performance in computing the bond dissociation energies, while ωB97X-D/cc-pVTZ exhibits the best prediction in enthalpy of formation of the studied reaction systems. The structural effect on the bond dissociation energies and enthalpy of formation of chlorinated and brominated polycyclic aromatic hydrocarbons are then systematically analyzed. Based on comparisons of the various methods, reliable DFT methods are recommended for future theoretical studies on large chlorinated and brominated polycyclic aromatic hydrocarbons considering both accuracy and computational cost. This work, to the authors' knowledge, is the first to systematically benchmark theoretical methods for the accurate prediction of thermodynamic properties for chlorinated and brominated polycyclic aromatic hydrocarbons.

Benchmark calculations using state-of-the-art DFT functionals and composite methods for bond dissociation energy and enthalpy of formation of halogenated polycyclic aromatic hydrocarbons are performed.

Electronic structure-based rate rules for ipso addition-elimination reactions on mono-aromatic hydrocarbons with single and double OH/CH3/OCH3/CHO/C2H5 substituents: a systematic theoretical investigation

The recent interest in bio-oils combustion and the key role of mono-aromatic hydrocarbons (MAHs) in existing kinetic frameworks, both in terms of poly-aromatic hydrocarbons growth and surrogate fuels formulation, motivates the current systematic theoretical investigation of one of the relevant reaction classes in MAHs pyrolysis and oxidation: ipso substitution by hydrogen. State-of-the-art theoretical methods and protocols implemented in automatized computational routines allowed to investigate 14 different potential energy surfaces involving MAHs with hydroxy and methyl single (phenol and toluene) and double (o-,m-,p-C6H4(OH)2, o-,m-,p-CH3C6H4OH, and o-,m-,p-C6H4(CH3)2) substituents, providing rate constants for direct implementation in existing kinetic models. The accuracy of the adopted theoretical method was validated by comparison of the computed rate constants with the available literature data. Systematic trends in energy barriers, pre-exponential factors, and temperature dependence of the Arrhenius parameters were found, encouraging the formulation of rate rules for ipso substitutions on MAHs. The rules here proposed allow to extrapolate from a reference system the necessary activation energy and pre-exponential factor corrections for a large number of reactions from a limited set of electronic structure calculations. We were able to estimate rate constants for other 63 ipso addition-elimination reactions on di-substituted MAHs, reporting in total 75 rate constants for ipso substitution reactions o-,m-,p-R'C6H4R + → C6H5R + ', with R,R' = OH/CH3/OCH3/CHO/C2H5, in the 300-2000 K range. Additional calculations performed for validation showed that the proposed rate rules are in excellent agreement with the rate constants calculated using the full computational protocol in the 500-2000 K range, generally with errors below 20%, increasing up to 40% in a few cases. The main results of this work are the successful application of automatized electronic structure calculations for the derivation of accurate rate constants for ipso substitution reactions on MAHs, and an efficient and innovative approach for rate rules formulation for this reaction class.

Kinetics of the Hydrogen Abstraction PAH + •OH → PAH Radical + H2O Reaction Class: An Application of the Reaction Class Transition State Theory (RC-TST) and Structure-Activity Relationship (SAR)

A reaction class transition state theory (RC-TST) augmented with structure-activity relationship (SAR) methodology is applied to predict high-pressure limit thermal rate constants for hydrogen abstraction by ^•OH radical from polycyclic aromatic hydrocarbons (PAHs) reaction class in the temperature range of 300-3000 K. The rate constants for the reference reaction of C₆H₆ + ^•OH → C₆H₅ + H₂O is calculated by the canonical variational transition state theory (CVT) with small curvature tunneling (SCT). Only the reaction energy is needed to predict RC-TST rates for other processes within the family, the parameters needed were obtained from M06-2X/cc-pVTZ data for a training set of 34 reactions. The systematic error of the resulting RC-TST rates is smaller than 50% in comparison with explicit rate calculations, which facilitates application of the proposed methodology to the automated reaction mechanism generators (ARMGs) schemes.

Highly efficient gas-phase reactivity of protonated pyridine radicals with propene

Small nitrogen containing heteroaromatics are fundamental building blocks for many biological molecules, including the DNA nucleotides. Pyridine, as a prototypical N-heteroaromatic, has been implicated in the chemical evolution of many extraterrestrial environments, including the atmosphere of Titan. This paper reports on the gas-phase ion-molecule reactions of the three dehydro-N-pyridinium radical cation isomers with propene. Photodissociation ion-trap mass spectrometry experiments are used to measure product branching ratios and reaction kinetics. Reaction efficiencies for 2-dehydro-N-pyridinium, 3-dehydro-N-pyridinium and 4-dehydro-N-pyridinium with propene are 70%, 47% and 41%, respectively. The m/z 106 channel is the major product channel across all cases and assigned 2-, 3-, and 4-vinylpyridinium for each reaction. The m/z 93 channel is also significant and assigned the 2-, 3-, and 4-N-protonated-picolyl radical cation for each case. H-Abstraction from propene is not competitive under experimental conditions. Potential energy schemes, at the M06-2X/6-31(2df,p) level of theory and basis set, are described to assist in rationalising observed product branching ratios and elucidating possible reaction mechanisms. Reaction barriers to the production of vinylpyridinium (m/z 106) + CH₃ are the lowest identified for the 3- and 4-dehydro-N-pyridinium reactions, in support of the observed dominance of the m/z 106 ion signal. Ethylene loss via ring-mediated H-transfer along the propyl group is found to be the lowest energy pathway for the 2-dehydro-N-pyridinium reaction, suggesting a preference toward m/z 93 (N-protonated-picolyl radical cation) over the experimentally observed products. Entropic bottle-necks along the m/z 93 pathway however, associated with ring-mediated H-atom transfer, are responsible for the dominance of m/z 106 in the 2-dehydro-N-pyridinium + propene reaction. For all three isomers, computed barriers for all observed reaction channels were below the entrance channel, suggesting these reactions can contribute to molecular weight growth in extraterrestrial environments with accelerated reaction rates in low temperature regions of space.

Atmospheric oxidation of halogenated aromatics: comparative analysis of reaction mechanisms and reaction kinetics

Atmospheric transport is the major route for global distribution of semi-volatile compounds such as halogenated aromatics as well as their major exposure route for humans. Their major atmospheric removal process is oxidation by hydroxyl radicals. There is very little information on the reaction mechanism or reaction-path dynamics of atmospheric degradation of halogenated benzenes. Furthermore, the measured reaction rate constants are missing for the range of environmentally relevant temperatures, i.e. 230-330 K. A series of recent theoretical studies have provided those valuable missing information for fluorobenzene, chlorobenzene, hexafluorobenzene and hexachlorobenzene. Their comparative analysis has provided additional and more general insight into the mechanism of those important tropospheric degradation processes as well as into the mobility, transport and atmospheric fate of halogenated aromatic systems. It was demonstrated for the first time that the addition of hydroxyl radicals to monohalogenated as well as to perhalogenated benzenes proceeds indirectly, via a prereaction complex and its formation and dynamics have been characterized including the respective transition-state. However, in fluorobenzene and chlorobenzene reactions hydroxyl radical hydrogen is pointing approximately to the center of the aromatic ring while in the case of hexafluorobenzene and hexachlorobenzene, unexpectedly, the oxygen is directed towards the center of the aromatic ring. The reliable rate constants are now available for all environmentally relevant temperatures for the tropospheric oxidation of fluorobenzene, chlorobenzene, hexafluorobenzene and hexachlorobenzene while pentachlorophenol, a well-known organic micropollutant, seems to be a major stable product of tropospheric oxidation of hexachlorobenzene. Their calculated tropospheric lifetimes show that fluorobenzene and chlorobenzene are easily removed from the atmosphere and do not have long-range transport potential while hexafluorobenzene seems to be a potential POP chemical and hexachlorobenzene is clearly a typical persistent organic pollutant.

Hydrogen-atom attack on phenol and toluene is ortho-directed

The reaction of H + phenol and H/D + toluene has been studied in a supersonic expansion after electric discharge. The (1 + 1') resonance-enhanced multiphoton ionization (REMPI) spectra of the reaction products, at m/z = parent + 1, or parent + 2 amu, were measured by scanning the first (resonance) laser. The resulting spectra are highly structured. Ionization energies were measured by scanning the second (ionization) laser, while the first laser was tuned to a specific transition. Theoretical calculations, benchmarked to the well-studied H + benzene → cyclohexadienyl radical reaction, were performed. The spectrum arising from the reaction of H + phenol is attributed solely to the ortho-hydroxy-cyclohexadienyl radical, which was found in two conformers (syn and anti). Similarly, the reaction of H/D + toluene formed solely the ortho isomer. The preference for the ortho isomer at 100-200 K in the molecular beam is attributed to kinetic, not thermodynamic effects, caused by an entrance channel barrier that is ∼5 kJ mol(-1) lower for ortho than for other isomers. Based on these results, we predict that the reaction of H + phenol and H + toluene should still favour the ortho isomer under elevated temperature conditions in the early stages of combustion (200-400 °C).

Kinetic study of the OH radical reaction with phenylacetylene

The reaction of the OH radical with phenylacetylene is studied over the 298-423 K temperature range and 1-7.5 Torr pressure range in a quasi-static reaction cell. The OH radical is generated by 266 nm photolysis of hydrogen peroxide (H2O2) or 355 nm photolysis of nitrous acid (HONO), and its concentration monitored using laser-induced fluorescence. The measured reaction rates are found to strongly depend on laser fluence at 266 nm. The 266 nm absorption cross-section of phenylacetylene is measured to be 1.29 (±0.71) × 10(-17) cm(2), prohibiting any accurate kinetic measurements at this wavelength. The rates are independent of laser fluence at 355 nm with an average value of 8.75 (±0.73) × 10(-11) cm(3) s(-1). The reaction exhibits no pressure or temperature dependence over the studied experimental conditions and is much faster than the estimated values presently used in combustion models. These results are consistent with the formation of a short lifetime intermediate that stabilizes by collisional quenching with the buffer gas. The structures of the most likely formed products are discussed based on both the computed energies for the OH-addition intermediates and previous theoretical investigations on similar chemical systems.

Mechanisms and reaction-path dynamics of hydroxyl radical reactions with aromatic hydrocarbons: the case of chlorobenzene

All geometries and energies significant for the first step of tropospheric degradation of chlorobenzene were characterized using the MP2/6-31+G(d,p) and G3 methods. A pre-reaction complex for the addition of OH radical to chlorobenzene was found and the associated transition state was determined for the first time. The reaction path for the association of OH radical and chlorobenzene into the pre-reaction complex was extrapolated from the selected low frequency normal mode of pre-reaction complex. The reaction rate constant for addition of OH radical to chlorobenzene was determined for the temperature range 230-330K, using RRKM theory and G3 energies. The calculated rate constants are in agreement with the experimental results. Regioselectivity was also determined for the title reaction from the ratio of respective reaction rates and our results are in very good agreement with the experimental results, which show the dominance of the ortho and para channels as well as a negligible contribution by the ipso channel.

A mechanistic and kinetic study on the formation of PBDD/Fs from PBDEs

This study presents a detailed mechanistic and kinetic investigation that explains the experimentally observed high yields of formation of polybrominated dibenzo-p-dioxins and dibenzofurans (PBDD/Fs) from the polybrominated diphenyl ethers (PBDEs), commonly deployed in brominated flame retardants (BFRs). Theoretical calculations involved the accurate meta hybrid functional of M05-2X. The previously suggested pathways of debromination and generation of bromophenols/bromophenoxys/bromobenzenes were found to be unimportant corridors for the formation of PBDD/Fs. A loss of an ortho Br or H atom from PBDEs, followed by a ring-closure reaction, is the most accessible pathway for the production of PBDFs via modest reaction barriers. The initially formed peroxy-type adduct (RO₂) is found to evolve in a complex, nevertheless very exoergic, mechanism to produce PBDDs. Results indicate that, degree and pattern of bromination, in the vicinity of the ether oxygen bridge, has a minor influence on governing mechanisms and that even fully brominated isomers of BFRs are capable of forming PBDD/Fs. We thoroughly discuss bimolecular reactions of PBDEs with Br and H, as well as the Br-displacement reaction by triplet oxygen. The rate of the Br-displacement reaction significantly exceeds that of the unimolecular inititiation reactions due to loss of ortho Br or H. Results presented herein address conclusively the intriguing question of how PBDEs form PBDD/Fs, a matter that has been in the center of much debate among environmental chemists.

Kinetics of the gas-phase reaction of OH with chlorobenzene

The kinetics of the reaction of hydroxyl radicals with chlorobenzene was studied experimentally using a pulsed laser photolysis/pulsed laser induced fluorescence technique over a wide range of temperatures, 298-670 K, and at pressures between 13.33 and 39.92 kPa. The bimolecular rate constants demonstrate different behavior at low and high temperatures. At room temperature, T = 298.8 +/- 1.5 K, the rate constant is equal to (6.02 +/- 0.34) x 10(-13) cm3 molecule(-1) s(-1); at high temperatures (474-670 K), the rate constant values are significantly lower and have a positive temperature dependence that can be described by an Arrhenius expression k1(T) = (1.01 +/- 0.35) x 10(-11) exp[(-2490 +/- 170 K)/T] cm3 molecule(-1) s(-1). This behavior is consistent with the low-temperature reaction being dominated by reversible addition and the high-temperature reaction representing abstraction and addition-elimination channels. The potential energy surface of the reaction was studied using quantum chemical methods, and a transition state theory model was developed for all reaction channels. The temperature dependences of the high-temperature rate constants obtained in calculations using the method of isodesmic reactions for transition states (IRTS) and the CBS-QB3 method are in very good agreement with experiment, with deviations smaller than the estimated experimental uncertainties. The G3//B3LYP-based calculated rate constants are in disagreement with the experimental values. The IRTS-based model was used to provide modified Arrhenius expressions for the temperature dependences of the rate constant for the abstraction and addition-elimination (Cl replacement) channels of the reaction.

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.

Thermochemical and Kinetic Analysis on the Reactions of O2 with Products from OH Addition to Isobutene, 2-Hydroxy-1,1-dimethylethyl, and 2-Hydroxy-2-methylpropyl Radicals: HO2 Formation from Oxidation of Neopentane, Part II

Unimolecular dissociation of a neopentyl radical to isobutene and methyl radical is competitive with the neopentyl association with O2 ((3)Sigma(g)-) in thermal oxidative systems. Furthermore, both isobutene and the OH radical are important primary products from the reactions of neopentyl with O2. Consequently, the reactions of O2 with the 2-hydroxy-1,1-dimethylethyl and 2-hydroxy-2-methylpropyl radicals resulting from the OH addition to isobutene are important to understanding the oxidation of neopentane and other branched hydrocarbons. Reactions that correspond to the association of radical adducts with O2((3)Sigma(g)-) involve chemically activated peroxy intermediates, which can isomerize and react to form one of several products before stabilization. The above reaction systems were analyzed with ab initio and density functional calculations to evaluate the thermochemistry, reaction paths, and kinetics that are important in neopentyl radical oxidation. The stationary points of potential energy surfaces were analyzed based on the enthalpies calculated at the CBS-Q level. The entropies, S(degrees)298, and heat capacities, C(p)(T), (0 <or= T/K <or= 1500), from vibration, translation, and external rotation contributions were calculated using statistical mechanics based on the vibrational frequencies and structures obtained from the density functional study. The hindered internal rotor contributions to S(degrees)298 and C(p)(T) were calculated by solving the Schrödinger equation with free rotor wave functions, and the partition coefficients were treated by direct integration over energy levels of the internal rotation potentials. Enthalpies of formation (DeltaH(f)(degrees)298) were determined using isodesmic reaction analysis. The DeltaH(f)(degrees)298 values of (CH3)2C*CH(2)OH, (CH3)2C(OO*)CH(2)OH, (CH3)2C(OH)C*H2, and (CH3)2C(OH)CH(2)OO* radicals were determined to be -23.3, -62.2, -24.2, and -61.8 kcal mol(-1), respectively. Elementary rate constants were calculated from canonical transition state theory, and pressure-dependent rate constants for multichannel reaction systems were calculated as functions of pressure and temperature using multifrequency quantum Rice-Ramsperger-Kassel (QRRK) analysis for k(E) and a master equation for pressure falloff. Kinetic parameters for intermediate and product formation channels of the above reaction systems are presented as functions of temperature and pressure.

Reaction of the Hydroxyl Radical with Phenol in Water Up to Supercritical Conditions

The rate constants for the reactions of phenol with the hydroxyl radical (OH*) in water have been measured from room temperature to 380 degrees C using electron pulse radiolysis and transient absorption spectroscopy. The reaction scheme designed to fit the data shows the importance of an equilibrium, giving back reactants (OH* radical and phenol) from the dihydroxycyclohexadienyl radical formed by their reaction, and the non-negligible contribution of the hydroxycyclohexadienyl radical absorption from H* atom addition. The accuracy of the reaction scheme and the reaction rate constants determined from it have been determined by the analysis of two different experiments, one under pure N2O atmosphere and the second under a mixture a N2O and O2. We report reaction rates for the H* and OH* radical addition to phenol, the formation of phenoxyl, the second-order recombination, the reaction of dihydroxycyclohexadienyl with O2, and the decay of the peroxyl adduct. Nearly all of the reaction rates deviate strongly from Arrhenius behavior.

High-Temperature Reactions of OH Radicals with Benzene and Toluene

The rate constants for the reactions of OH radicals with benzene and toluene have been measured directly by a shock tube/pulsed laser-induced fluorescence imaging method at high temperatures. The OH radicals were generated by the thermal decomposition of nitric acid or tert-butyl hydroperoxide. The derived Arrhenius expressions for the rate constants were k(OH + benzene) = 8.0 x 10(-11) exp(-26.6 kJ mol(-1)/RT) [908-1736 K] and k(OH + toluene) = 8.9 x 10(-11) exp(-19.7 kJ mol(-1)/RT) [919-1481 K] in the units of cubic centimeters per molecule per second. Transition-state theory (TST) calculations based on quantum chemically predicted energetics confirmed the dominance of the H-atom abstraction channel for OH + benzene and the methyl-H abstraction channel for OH + toluene in the experimental temperature range. The TST calculation indicated that the anharmonicity of the C-H-O bending vibrations of the transition states is essential to reproduce the observed rate constants. Possible implications to the other analogous H-transfer reactions were discussed.

Formation of phenol and carbonyls from the atmospheric reaction of OH radicals with benzene

The gas-phase reaction of OH radicals with benzene has been studied in a flow tube operated at 295 +/- 2 K and 950 mbar of synthetic air or O2. Ozonolysis of tetramethylethylene (dark reaction) with a measured OH radical yield of 0.92 +/- 0.08 or photolysis of methyl nitrite in the presence of NO served as the OH sources. For investigations in the presence of NOx, the conditions were chosen so that more than 95% of the OH/benzene adduct reacted with O2 even for the highest NO2 concentration occurring in the experiment. In the absence of NOx, a phenol yield from the reaction of OH radicals with benzene of 0.61 +/- 0.07 was measured by means of long-path FT-IR and UV spectroscopy over a wide range of experimental conditions. This yield was confirmed by measurements performed in the presence of NOx. Detected carbonyls were glyoxal, cis-butenedial and trans-butenedial with formation yields of 0.29 +/- 0.10, 0.08 +/- 0.03 and 0.023 +/- 0.007, respectively, measured in synthetic air and in the presence of NOx. There was no significant difference in the product yields applying both experimental approaches for OH generation (dark reaction or photolysis). Nitrobenzene and o-nitrophenol were detected in traces. The yield of nitrobenzene increased with increasing NOx resulting in a maximum formation yield of 0.007. The detected products in the presence of NOx account for approximately 78% of the reacted carbon. Butenedial yields from benzene degradation are reported for the first time. In the absence of NOx, glyoxal, cis-butenedial and trans-butenedial were also detected, but with distinctly lower yields compared to the experiments with NOx.