Ab initio study on the mechanism of the HCO+O2→HO2+CO reaction

Low temperature reaction kinetics inside an extended Laval nozzle: REMPI characterization and detection by broadband rotational spectroscopy

Chirped-Pulse Fourier-Transform millimeter wave (CP-FTmmW) spectroscopy is a powerful method that enables detection of quantum state specific reactants and products in mixtures. We have successfully coupled this technique with a pulsed uniform Laval flow system to study photodissociation and reactions at low temperature, which we refer to as CPUF ("Chirped-Pulse/Uniform flow"). Detection by CPUF requires monitoring the free induction decay (FID) of the rotational coherence. However, the high collision frequency in high-density uniform supersonic flows can interfere with the FID and attenuate the signal. One way to overcome this is to sample the flow, but this can cause interference from shocks in the sampling region. This led us to develop an extended Laval nozzle that creates a uniform flow within the nozzle itself, after which the gas undergoes a shock-free secondary expansion to cold, low pressure conditions ideal for CP-FTmmW detection. Impact pressure measurements, commonly used to characterize Laval flows, cannot be used to monitor the flow within the nozzle. Therefore, we implemented a REMPI (resonance-enhanced multiphoton ionization) detection scheme that allows the interrogation of the conditions of the flow directly inside the extended nozzle, confirming the fluid dynamics simulations of the flow environment. We describe the development of the new 20 K extended flow, along with its characterization using REMPI and computational fluid dynamics. Finally, we demonstrate its application to the first low temperature measurement of the reaction kinetics of HCO with O2 and obtain a rate coefficient at 20 K of 6.66 ± 0.47 × 10-11 cm3 molec-1 s-1.

Dynamics of the reaction CH2I + O2 probed via infrared emission of CO, CO2, OH and H2CO

The reaction CH2I + O2 has been widely employed recently for the production of the simplest Criegee intermediate CH2OO in laboratories, but the detailed dynamics of this reaction have been little explored. Infrared emission of several products of this reaction, initiated on irradiation of CH2I2 and O2 (∼8 Torr) in a flowing mixture at 308 or 248 nm, was recorded with a step-scan Fourier-transform spectrometer; possible routes of formation were identified according to the observed vibrational distribution of products and published theoretical potential-energy schemes. Upon irradiation at 308 nm, Boltzmann distributions of CO (v ≤ 5, J ≤ 19) with an average vibrational energy of 32 ± 3 kJ mol-1 and OH (v ≤ 3, J ≤ 5.5) with an average vibrational energy of 29 ± 4 kJ mol-1 were observed and assigned to the decomposition of HCOOH* to form CO + H2O and OH + HCO, respectively. The broadband emission of CO2 was simulated with two vibrational distributions of average energies (91 ± 4) and (147 ± 8) kJ mol-1 and assigned to be produced from the decomposition of HCOOH* and methylene bis(oxy), respectively. Upon irradiation of samples at 248 nm, the emission of OH and CO2 showed similar distributions with slightly greater energies, but the distribution of CO (v ≤ 11, J ≤ 19) became bimodal with average vibrational energies of (23 ± 4) and (107 ± 29) kJ mol-1, and branching (56 ± 5) : (44 ± 5). The additional large-v component is assigned to be produced from a secondary reaction HCO + O2 to form CO + HO2; HCO is a coproduct of OH. The branching between CO and OH is (50 ± 5) : (50 ± 5) at 308 nm and (64 ± 5) : (36 ± 4) at 248 nm, consistent with the mechanism according to which an additional channel to produce CO opens at 248 nm. Highly internally excited H2CO was also observed. With O2 at 16 Torr, the extrapolated nascent internal distributions are similar to those with O2 at 8 Torr except for a slight quenching effect.

Experimental Evidence of Dioxole Unimolecular Decay Pathway for Isoprene-Derived Criegee Intermediates

Ozonolysis of isoprene, one of the most abundant volatile organic compounds emitted into the Earth's atmosphere, generates two four-carbon unsaturated Criegee intermediates, methyl vinyl ketone oxide (MVK-oxide) and methacrolein oxide (MACR-oxide). The extended conjugation between the vinyl substituent and carbonyl oxide groups of these Criegee intermediates facilitates rapid electrocyclic ring closures that form five-membered cyclic peroxides, known as dioxoles. This study reports the first experimental evidence of this novel decay pathway, which is predicted to be the dominant atmospheric sink for specific conformational forms of MVK-oxide (anti) and MACR-oxide (syn) with the vinyl substituent adjacent to the terminal O atom. The resulting dioxoles are predicted to undergo rapid unimolecular decay to oxygenated hydrocarbon radical products, including acetyl, vinoxy, formyl, and 2-methylvinoxy radicals. In the presence of O₂, these radicals rapidly react to form peroxy radicals (ROO), which quickly decay via carbon-centered radical intermediates (QOOH) to stable carbonyl products that were identified in this work. The carbonyl products were detected under thermal conditions (298 K, 10 Torr He) using multiplexed photoionization mass spectrometry (MPIMS). The main products (and associated relative abundances) originating from unimolecular decay of anti-MVK-oxide and subsequent reaction with O₂ are formaldehyde (88 ± 5%), ketene (9 ± 1%), and glyoxal (3 ± 1%). Those identified from the unimolecular decay of syn-MACR-oxide and subsequent reaction with O₂ are acetaldehyde (37 ± 7%), vinyl alcohol (9 ± 1%), methylketene (2 ± 1%), and acrolein (52 ± 5%). In addition to the stable carbonyl products, the secondary peroxy chemistry also generates OH or HO₂ radical coproducts.

Ab initio kinetics of the HOSO2 + 3O2 → SO3 + HO2 reaction

The detailed kinetic mechanism of the HOSO₂ + ³O₂ reaction, which plays a pivotal role in the atmospheric oxidation of SO₂, was investigated using accurate electronic structure calculations and novel statistical thermodynamic/kinetic models. Explored using the accurate composite method W1U, the detailed potential energy surface (PES) revealed that the addition of O₂ to a HOSO₂ radical to form the adduct (HOSO₄) proceeds via a transition state with a slightly positive barrier (i.e., 0.7 kcal mol^-1 at 0 K). Such a finding compromises a long-term hypothesis about this channel of being a barrierless process. Moreover, the overall reaction was found to be slightly exothermic by 1.7 kcal mol^-1 at 0 K, which is in good agreement with recent studies. On the newly-constructed PES, the temperature- and pressure-dependent behaviors of the title reaction were characterized in a wide range of conditions (T = 200-1000 K & P = 10-760 Torr) using the integrated deterministic and stochastic master equation/Rice-Ramsperger-Kassel-Marcus (ME/RRKM) rate model in which corrections for hindered internal rotation (HIR) and tunneling treatments were included. The calculated numbers were found to be in excellent agreement with literature data. The sensitivity analyses on the derived rate coefficients with respect to the ab initio input parameters (i.e., barrier height and energy transfer) were also performed to further understand the kinetic behaviors of the title reaction. The detailed kinetic mechanism, consisting of thermodynamic and kinetic data (in NASA polynomial and modified Arrhenius formats, respectively), was also provided at different T & P for further use in the modeling/simulation of any related systems.

Silane-initiated nucleation in chemically active plasmas: validation of density functionals, mechanisms, and pressure-dependent variational transition state calculations

Pressure-dependent rate constants for nucleation in nanodusty plasmas are calculated by variational transition state theory with system-specific quantum RRK theory.

Multi-path variational transition state theory for chiral molecules: the site-dependent kinetics for abstraction of hydrogen from 2-butanol by hydroperoxyl radical, analysis of hydrogen bonding in the transition state, and dramatic temperature dependence of the activation energy

The goal of the present work is modeling the kinetics of a key reaction involved in the combustion of the biofuel 2-butanol. To accomplish this we extended multi-path variational transition state theory (MP-VTST) with the small curvature tunneling (SCT) approximation to include multistructural anharmonicity factors for molecules with chiral carbons. We use the resulting theory to predict the site-dependent rate constants of the hydrogen abstraction from 2-butanol by hydroperoxyl radical. The generalized transmission coefficients were averaged over the four lowest-energy reaction paths. The computed forward reaction rate constants indicate that hydrogen abstraction from the C-2 site has the largest contribution to the overall reaction from 200 K to 2400 K, with a contribution ranging from 99.9988% at 200 K to 88.9% at 800 K to 21.2% at 3000 K, while hydrogen abstraction from the oxygen site makes the lowest contribution at all temperatures, ranging from 2.5 × 10^-9% at 200 K to 0.65% at 800 K to 18% at 3000 K. This work highlights the importance of including the multiple-structure and torsional potential anharmonicity in the computation of the thermal rate constants. We also analyzed the role played by the hydrogen bond at the transition state, and we illustrated the risks of (a) considering only the lowest-energy conformations in the calculations of the rate constants or (b) ignoring the nonlinear temperature dependence of the activation energies. A hydrogen bond at the transition state can lower the enthalpy of activation, but raise the free energy of activation. We find an energy of activation that increases from 11 kcal mol^-1 at 200 K to more than 36 kcal mol^-1 at high temperature for this radical reaction with a biofuel molecule.

Mechanism and Thermodynamics of Multichannel 1:1 Ammonia and Ozone Tropospheric Oxidation Reaction

The reaction of NH3 + O3 on the singlet potential energy surface has been investigated theoretically using CCSD(T) and G3B3 methods with the 6-311++G(3df,3pd) basis set based on the geometric parameters optimised at the B3LYP/6-311++G(3df,3pd) level of theory. The calculated results revealed that the reactants are firstly associated with the adduct NH3–O3, through a barrier-less process and two stable collision complexes, C1 and C2 have been considered. Six kinds of products P1 (HOO + H2NO), P2 (OH + H2NO2), P3 (O2 + H3NO), P4 (HNO + H2O2), P5 (HNO2 + H2O) and P6 (NO + HO2 + H2) are obtained through a variety of transformations of C1 and C2. Thermodynamic parameters of all possible products show five production channels which are spontaneous. After taking into account the reaction barrier and enthalpy, the most probable reaction pathways are P1(2) and P4. Also, P4 is the most stable product in comparison with the others and P1 yields the other important products from the kinetic viewpoint; this is relatively stable, with a negative Gibbs free energy. Branching ratio results show P1 is the major product at low temperature. Finally, the rates constant have been computed only for the P1(2) and P4 pathways by related theories in the temperature range of 200–2500 K.

RADICAL-MOLECULE REACTIONS HCO/HOC + C2H4: A MECHANISTIC STUDY

A detailed computational study is performed on the radical-molecule reactions between HCO/HOC and ethylene ( C 2 H 4) at the Gaussian-3//B3LYP/6-31G(d) level. For the HCO + C 2 H 4 reaction, the most favorable pathway is the direct C -addition forming the intermediate H 2 CCH 2 CHO , followed by a 1,2- H -shift leading to H 3 CCHCHO . Subsequently, there are two highly competitive dissociation pathways for H 3 CCHCHO : one is the formation of the direct H -extrusion product H 2 CCHCHO + H , and the other is the formation of C 2 H 5 + CO via the intermediate H 3 CCH 2 CO . The overall reaction barrier is 14.1 and 14.6 kcal/mol respectively, at the G3B3 level. The quasi-direct H -donation process to produce C 2 H 5 + CO with the barrier 16.5 kcal/mol is less competitive. Thus, only at higher temperatures, the HCO + C 2 H 4 reaction could play a role. In contrast, the HOC + C 2 H 4 reaction just need to overcome a small barrier 2.0 kcal/mol to generate C 2 H 5 + CO via the quasi-direct H -donation mechanism. This is suggestive of the potential importance of the HOC + C 2 H 4 reaction in combustion processes. However, the direct C -addition channel is much less competitive. The present kinetic data and orbital analysis show that the HCO radical has much higher reactivity than HOC , although the latter is more energetic. Till now, no kinetic study on the HOC radical has been reported, the present study can provide useful information on understanding the reactivity and depletion mechanism of the energetic HOC radical.

Bond-forming reactions of dications with molecules: a computational and experimental study of the mechanisms for the formation of HCF2+ from CF3(2+) and H2

The QCISD and QCISD(T) quantum chemical methods have been used to characterize the energetics of various possible mechanisms for the formation of HCF2+ from the bond-forming reaction of CF3(2+) with H2. The stationary points on four different pathways leading to the product combinations HCF2+ + H+ + F and HCF2+ + HF+ have been calculated. All four pathways begin with the formation of a collision complex [H2-CF3]2+, followed by an internal hydrogen atom migration to give HC(FH)F2(2+). In two of the mechanisms, immediate charge separation of HC(FH)F2(2+) via loss of either HF+ or a proton, followed by loss of an F atom, yields the experimentally observed bond-forming product HCF2+. For the other two mechanisms, internal hydrogen rearrangement of HC(FH)F2(2+) to give C(FH)2F(2+), followed by charge separation, yields the product CF2H+. This product can then overcome a 2.04 eV barrier to rearrange to the HCF2+ isomer, which is 1.80 eV more stable. All four calculated mechanisms are in agreement with the isotope effects and collision energy dependencies of the product ion cross sections that have been previously observed experimentally following collisions between CF3(2+) and H2/D2. We find that in this open-shell system, CCSD(T) and QCISD(T) T1-diagnostic values of up to 0.04 are acceptable. A series of angularly resolved crossed-beam scattering experiments on collisions of CF3(2+) with D2 have also been performed. These experiments show two distinct channels leading to the formation of DCF2+. One channel appears to correspond to the pathway leading to the ground state 1DCF2+ + D+ + F product asymptote and the other to the 3DCF2+ + D+ + F product asymptote, which is 5.76 eV higher in energy. The experimental kinetic energy releases for these channels, 7.55 and 1.55 eV respectively, have been determined from the velocities of the DCF2+ product ion and are in agreement with the reaction mechanisms calculated quantum chemically. We suggest that both of these observed experimental channels are governed by the reaction mechanism we calculate in which charge separation occurs first by loss of a proton, without further hydrogen atom rearrangement, followed by loss of an F atom to give the final products 1DCF2+ + D+ + F or 3DCF2+ + D+ + F.

Radical−Molecule Reactions HCO/HOC + C2H2: Mechanistic Study

A detailed computational study is performed on the unknown radical-molecule reactions between HCO/HOC and acetylene (C2H2) at the CCSD(T)/6-311G(2d,p)//B3LYP/6-311G(d,p)+ZPVE, Gaussian-3//B3LYP/6-31G(d), and Gaussian-3//MP2(full)/6-31G(d) levels. For the HCO + C2H2 reaction, the most favorable pathway is direct C-addition forming the intermediate HC=CHCH=O followed by a 1,3-H-shift leading to H2C=CHC=O, which finally dissociates to the product C2H3 + CO. The overall reaction barrier is 13.8, 10.5, and 11.3 kcal/mol, respectively, at the three levels. The quasi-direct H-donation process to produce C2H3 + CO with barriers of 14.0, 14.1, and 14.1 kcal/mol is less competitive. Thus only at higher temperatures could the HCO + C2H2 reaction play a role. In contrast, the HOC + C2H2 reaction can barrierlessly generate C2H3 + CO via the quasi-direct H-donation mechanism proceeding via a prereactive complex with OH...C2 hydrogen bonding. This is suggestive of the potential importance of the HOC + C2H2 reaction in both combustion and interstellar processes. However, the direct C-addition channel is much less competitive. For both reactions, the possible formation of the intriguing interstellar molecules propadiene and propynal is also discussed. The present theoretical study represents the first attempt to probe the reaction mechanism between HOC and pi-systems. Future laboratory investigations on both reactions (particularly HOC + C2H2) are recommended.

Mechanism of HCS + O2 reaction: Hydrogen- or oxygen-transfer?

In spite of the potential importance of the HCS radical in both combustion and interstellar processes, its chemical reactivity has not been tackled previously. In the present paper, the oxidation reaction of the HCS radical is theoretically investigated for the first time at the CCSD(T)/6-311++G(3df,2p)//BH&HLYP/6-311++G(d,p)+ZPVE and Gaussian-3//B3LYP/6-31G(d) levels. It is shown that the most feasible pathway is the O2 addition to the HCS radical forming the intermediate SC(H)OO which can undergo a subsequent O-extrusion leading to SC(H)O + 3O. This features an indirect O-transfer mechanism with the overall barrier of 4.4 and 3.5 kcal mol(-1), respectively, at the two levels. However, formation of the H-transfer product CS + HO2 is kinetically much less feasible, i.e., the direct mechanism has barriers of 14.3 and 8.7 kcal mol(-1), whereas the indirect mechanism has barriers of 12.6 and 10.7 kcal mol(-1), respectively. This result is in sharp contrast to the analogous HCO + O2 reaction, where the direct (with a barrier of 2.98 kcal mol(-1)) and indirect (2.26 kcal mol(-1)) H-transfer processes are highly competitive over the indirect O-transfer process (the least endothermicity is 19.9 kcal mol(-1)). The possible explanations and implications of the present results are provided.