From the Multiple-Well Master Equation to Phenomenological Rate Coefficients: Reactions on a C3H4 Potential Energy Surface

Chemistry of energetically activated cumulenes - from allene (H2CCCH2) to hexapentaene (H2CCCCCCH2)

During the last decade, experimental and theoretical studies on the unimolecular decomposition of cumulenes (H(2)C(n)H(2)) from propadiene (H(2)CCCH(2)) to hexapentaene (H(2)CCCCCCH(2)) have received considerable attention due to the importance of these carbon-bearing molecules in combustion flames, chemical vapor deposition processes, atmospheric chemistry, and the chemistry of the interstellar medium. Cumulenes and their substituted counterparts also have significant technical potential as elements for molecular machines (nanomechanics), molecular wires (nano-electronics), nonlinear optics, and molecular sensors. In this review, we present a systematic overview of the stability, formation, and unimolecular decomposition of chemically, photo-chemically, and thermally activated small to medium-sized cumulenes in extreme environments. By concentrating on reactions under gas phase thermal conditions (pyrolysis) and on molecular beam experiments conducted under single-collision conditions (crossed beam and photodissociation studies), a comprehensive picture on the unimolecular decomposition dynamics of cumulenes transpires.

Pressure and temperature dependence of the reaction of vinyl radical with ethylene

This work reports measurements of absolute rate coefficients and Rice-Ramsperger-Kassel-Marcus (RRKM) master equation simulations of the C2H3+C2H4 reaction. Direct kinetic studies were performed over a temperature range of 300-700 K and pressures of 20 and 133 mbar. Vinyl radicals (H2C=CH) were generated by laser photolysis of vinyl iodide (C2H3I) at 266 nm, and time-resolved absorption spectroscopy was used to probe vinyl radicals through absorption at 423.2 nm. Measurements at 20 mbar are in good agreement with previous determinations at higher temperature. A weighted three-parameter Arrhenius fit to the experimental rate constant at 133 mbar, with the temperature exponent fixed, gives k=(7+/-1)x10(-14) cm3 molecule(-1) s(-1) (T/298 K)2 exp[-(1430+/-70) K/T]. RRKM master equation simulations, based on G3 calculations of stationary points on the C4H7 potential energy surface, were carried out to predict rate coefficients and product branching fractions. The predicted branching to 1-methylallyl product is relatively small under the conditions of the present experiments but increases as the pressure is lowered. Analysis of end products of 248 nm photolysis of vinyl iodide/ethylene mixtures at total pressures between 27 and 933 mbar provides no direct evidence for participation of 1-methylallyl.

The Reaction of n- and i-C4H5 Radicals with Acetylene

In this article, we discuss the reactions of i-C4H5 and n-C4H5 with acetylene. Both have been proposed as possible cyclization steps, forming benzene or fulvene, in rich flames burning aliphatic fuels. The relevant parts of the potential energy surface were determined from rQCISD(T) calculations extrapolated to the infinite-basis-set limit. Using this information in a Rice-Ramsperger-Kassel-Marcus-based master equation, we have calculated thermal rate coefficients and product distributions for both reactions as a function of temperature and pressure. The results are cast in forms that can be used in modeling, and the implications of the results for flame chemistry are discussed.

Thermochemistry of Key Soot Formation Intermediates: C3H3 Isomers

Accurate standard enthalpies of formation for allene, propyne, and four C3H3 isomers involved in soot formation mechanisms have been determined through systematic focal point extrapolations of ab initio energies. Auxiliary corrections have been applied for anharmonic zero-point vibrational energy, core electron correlation, the diagonal Born-Oppenheimer correction (DBOC), and scalar relativistic effects. Electron correlation has been accounted for via second-order Z-averaged perturbation theory (ZAPT2) and primarily through coupled-cluster theory, including single, double, and triple excitations, as well as a perturbative treatment of connected quadruple excitations [ROCCSD, ROCCSD(T), ROCCSDT, and UCCSDT(Q)]. The correlation-consistent hierarchy of basis sets, cc-pVXZ (X = D, T, Q, 5, 6), was employed. The CCSDT(Q) corrections do not exceed 0.12 kcal mol(-)1 for the relative energies of the systems considered here, indicating a high degree of electron correlation convergence in the present results. Our recommended values for the enthalpies of formation are as follows: Delta(f)H(o)(0)(propargyl) = 84.76, Delta(f)H(o)(0) (1-propynyl) = 126.60, Delta(f)H(o)(0) (cycloprop-1-enyl) = 126.28, Delta(f)H(o)(0)(cycloprop-2-enyl) = 117.36, Delta(f)H(o)(0)(allene) = 47.41, and Delta(f)H(o)(0)(propyne) = 46.33 kcal mol(-1), with estimated errors no larger than 0.3 kcal mol(-1). The corresponding C3H3 isomerization energies are about 1 kcal mol(-1) larger than previous coupled-cluster results and several kcal mol(-1) below those previously obtained using density functional theory.

Reaction of Hydrogen Atoms with Propyne at High Temperatures: An Experimental and Theoretical Study

The kinetics of the reaction of hydrogen atoms with propyne (pC3H4) was experimentally studied in a shock tube at temperatures ranging from 1200 to 1400 K and pressures between 1.3 and 4.0 bar with Ar as the bath gas. The hydrogen atoms (initial mole fraction 0.5-2.0 ppm) were produced by pyrolysis of C2H5I and monitored by atomic resonance absorption spectrometry under pseudo-first-order conditions with respect to propyne (initial mole fraction 5-20 ppm). From the hydrogen atom time profiles, overall rate coefficients k(ov) identical with -([pC3H4][H])(-1) x d[H]/dt for the reaction H + pC3H4 --> products ( not equal H) were deduced; the following temperature dependence was obtained: kov = 1.2 x 10(-10) exp(-2270 K/T) cm(3) s(-1) with an estimated uncertainty of +/-20%. A pressure dependence was not observed. The results are analyzed in terms of statistical rate theory with molecular and transition state data from quantum chemical calculations. Geometries were optimized using density functional theory at the B3LYP/6-31G(d) level, and single-point energies were computed at the QCISD(T)/cc-pVTZ level of theory. It is confirmed that the reaction proceeds via an addition-elimination mechanism to yield C2H2 + CH3 and via a parallel direct abstraction to give C3H3 + H2. Furthermore, it is shown that a hydrogen atom catalyzed isomerization channel to allene (aC3H4), H + pC3H4 --> aC3H4 + H, is also important. Kinetic parameters to describe the channel branching of these reactions are deduced.

Kinetics of the Reaction of Methyl Radical with Hydroxyl Radical and Methanol Decomposition

The CH3 + OH bimolecular reaction and the dissociation of methanol are studied theoretically at conditions relevant to combustion chemistry. Kinetics for the CH3 + OH barrierless association reaction and for the H + CH2OH and H + CH3O product channels are determined in the high-pressure limit using variable reaction coordinate transition state theory and multireference electronic structure calculations to evaluate the fragment interaction energies. The CH3 + OH --> 3CH2 + H2O abstraction reaction and the H2 + HCOH and H2 + H2CO product channels feature localized dynamical bottlenecks and are treated using variational transition state theory and QCISD(T) energies extrapolated to the complete basis set limit. The 1CH2 + H2O product channel has two dynamical regimes, featuring both an inner saddle point and an outer barrierless region, and it is shown that a microcanonical two-state model is necessary to properly describe the association rate for this reaction over a broad temperature range. Experimental channel energies for the methanol system are reevaluated using the Active Thermochemical Tables (ATcT) approach. Pressure dependent, phenomenological rate coefficients for the CH3 + OH bimolecular reaction and for methanol decomposition are determined via master equation simulations. The predicted results agree well with experimental results, including those from a companion high-temperature shock tube determination for the decomposition of methanol.

Direct Measurement and Theoretical Calculation of the Rate Coefficient for Cl + CH3 in the Range from T = 202−298 K

The rate coefficient has been measured under pseudo-first-order conditions for the Cl+CH3 association reaction at T=202, 250, and 298 K and P=0.3-2.0 Torr helium using the technique of discharge-flow mass spectrometry with low-energy (12-eV) electron-impact ionization and collision-free sampling. Cl and CH3 were generated rapidly and simultaneously by reaction of F with HCl and CH4, respectively. Fluorine atoms were produced by microwave discharge in an approximately 1% mixture of F2 in He. The decay of CH3 was monitored under pseudo-first-order conditions with the Cl-atom concentration in large excess over the CH3 concentration ([Cl]0/[CH3]0=9-67). Small corrections were made for both axial and radial diffusion and minor secondary chemistry. The rate coefficient was found to be in the falloff regime over the range of pressures studied. For example, at T=202 K, the rate coefficient increases from 8.4x10(-12) at P=0.30 Torr He to 1.8x10(-11) at P=2.00 Torr He, both in units of cm3 molecule-1 s-1. A combination of ab initio quantum chemistry, variational transition-state theory, and master-equation simulations was employed in developing a theoretical model for the temperature and pressure dependence of the rate coefficient. Reasonable empirical representations of energy transfer and of the effect of spin-orbit interactions yield a temperature- and pressure-dependent rate coefficient that is in excellent agreement with the present experimental results. The high-pressure limiting rate coefficient from the RRKM calculations is k2=6.0x10(-11) cm3 molecule-1 s-1, independent of temperature in the range from 200 to 300 K.

Direct dynamics study on the reaction of acetaldehyde with ozone

The hydrogen abstraction reaction of ozone with acetaldehyde has been studied theoretically over the temperature range 250-2500 K. Two different reactive sites of acetaldehyde molecule, CH(3) and CHO groups have been investigated, and results confirm that the CHO group is a highly reactive site. In this study, the geometries and harmonic vibrational frequencies of all stationary points are calculated at the MPW1K, BHandHLYP, and MPWB1K levels of theory. The minimum energy paths (MEPs) were obtained at the MPW1K/6-31+G(d,p) level of theory. To refine the energies along the MEPs of each channel, single-point energy calculations were performed by a higher-level energy calculation method (denoted as HL). The rate constants were evaluated based on the MEPs from the HL method in the temperature range 250-2500 K by using the conventional transition state theory (TST), the canonical variational transition state theory (CVT), the microcanonical variational transition state theory (muVT), the CVT coupled with small-curvature tunneling (SCT) correction (CVT/SCT), and the muVT coupled with Eckart tunneling correction (muVT/Eckart). The fitted three-parameter Arrhenius expressions of the calculated CVT/SCT and muVT/Eckart rate constants of the H abstraction from CHO group are k CVT/SCT(T) = 4.92 x 10(-27).T 3.77.e(-7867.0/T) and k muVT/Eckart(T) = 2.10 x 10(-27).T(3.90).e(-7706.2/T), respectively. The fitted three-parameter Arrhenius expressions of the calculated CVT/SCT and muVT/Eckart rate constants of the H abstraction from CH3 group are k(CVT/SCT)(T) = 1.27 x 10(-27).T(3.94).e(-14554.1/T) and k muVT/Eckart(T) = 1.62 x 10(-26).T(3.66).e(-15459.8/T), respectively.

Theory, measurements, and modeling of OH and HO2 formation in the reaction of cyclohexyl radicals with O2

The production of OH and HO(2) in Cl-initiated oxidation of cyclohexane has been measured using pulsed-laser photolytic initiation and continuous-laser absorption detection. The experimental data are modeled by master equation calculations that employ new G2(MP2)-like ab initio characterizations of important stationary points on the cyclo-C(6)H(11)O(2) surface. These ab initio calculations are a substantial expansion on previously published characterizations, including explicit consideration of conformational changes (chair-boat, axial-equatorial) and torsional potentials. The rate constants for the decomposition and ring-opening of cyclohexyl radical are also computed with ab initio based transition state theory calculations. Comparison of kinetic simulations based on the master equation results with the present experimental data and with literature determinations of branching fractions suggests adjustment of several transition state energies below their ab initio values. Simulations with the adjusted values agree well with the body of experimental data. The results once again emphasize the importance of both direct and indirect components of the kinetics for the production of both HO(2) and OH in radical + O(2) reactions.

Master equation methods for multiple well systems: application to the 1-,2-pentyl system

The master equation (ME) provides a powerful technique for modeling reactions that involve at least one potential energy well. It can be widely applied to reactions with several connected energy wells and multiple product channels. The application of the technique is reviewed by reference to the H + SO(2) reaction, where phenomenological rate constants for use, for example, in a combustion model can be extracted through an analysis of the eigenvalues and eigenvectors of the collision matrix, M, that describes formation of the adducts HSO(2) and HOSO from the source H + SO(2), collisional energy transfer in the adduct wells and reaction via the product channel (sink) OH + SO. The approach is extended to systems with more than one sink and it is demonstrated that macroscopic (phenomenological) rate coefficients derived from a ME obey detailed balance if the original ME is appropriately constructed. The method has been applied to the 1-, 2-pentyl radical system, that includes isomerisation and dissociation via two channels to form C(3)H(6) + C(2)H(5) and C(2)H(4) + C(3)H(7). The calculations clearly demonstrate the importance of indirect dissociation channels, in which an isomer can dissociate to form the product set to which it is not directly connected, e.g. formation of C(3)H(6) + C(2)H(5) from 1-pentyl, via the energized states of 2-pentyl. As in previous studies of pentyl dissociation, there is a convergence of the chemically significant eigenvalues and the internal energy relaxation eigenvalues above approximately 1000 K; the consequences of this convergence are discussed.

Shock tube study of dissociation and relaxation in 1,1-difluoroethane and vinyl fluoride

This paper reports measurements of the thermal dissociation of 1,1-difluoroethane in the shock tube. The experiments employ laser-schlieren measurements of rate for the dominant HF elimination using 10% 1,1-difluoroethane in Kr over 1500-2000 K and 43 < P < 424 torr. The vinyl fluoride product of this process then dissociates affecting the late observations. We thus include a laser schlieren study (1717-2332 K, 75 < P < 482 torr in 10 and 4% vinyl fluoride in Kr) of this dissociation. This latter work also includes a set of experiments using shock-tube time-of-flight mass spectrometry (4% vinyl fluoride in neon, 1500-1980 K, 500 < P < 1300 torr). These time-of-flight experiments confirm the theoretical expectation that the only reaction in vinyl fluoride is HF elimination. The dissociation experiments are augmented by laser schlieren measurements of vibrational relaxation (1-20% C(2)H(3)F in Kr, 415-1975 K, 5 < P < 50 torr, and 2 and 5% C(2)H(4)F(2) in Kr, 700-1350 K, 6 < P < 22 torr). These experiments exhibit very rapid relaxation, and incubation delays should be negligible in dissociation. An RRKM model of dissociation in 1,1-difluoroethane based on a G3B3 calculation of barrier and other properties fits the experiments but requires a very large DeltaE(down) of 1600 cm(-1), similar to that found in a previous examination of 1,1,1-trifluoroethane. Dissociation of vinyl fluoride is complicated by the presence of two parallel HF eliminations, both three-center and four-center. Structure calculations find nearly equal barriers for these, and TST calculations show almost identical k(infinity). An RRKM fit to the observed falloff again requires an unusually large DeltaE(down) and the experiments actually support a slightly reduced barrier. These large energy-transfer parameters now seem routine in these large fluorinated species. It is perhaps a surprising result for which there is as yet no explanation.

A combined experimental and theoretical study of the reaction between methylglyoxal and OH/OD radical: OH regeneration

Experimental studies have been conducted to determine the rate coefficient and mechanism of the reaction between methylglyoxal (CH(3)COCHO, MGLY) and the OH radical over a wide range of temperatures (233-500 K) and pressures (5-300 Torr). The rate coefficient is pressure independent with the following temperature dependence: k(3)(T) = (1.83 +/- 0.48) x 10(-12) exp((560 +/- 70)/T) cm(3) molecule(-1) s(-1) (95% uncertainties). Addition of O(2) to the system leads to recycling of OH. The mechanism was investigated by varying the experimental conditions ([O(2)], [MGLY], temperature and pressure), and by modelling based on a G3X potential energy surface, rovibrational prior distribution calculations and master equation RRKM calculations. The mechanism can be described as follows: Addition of oxygen to the system shows that process (4) is fast and that CH(3)COCO completely dissociates. The acetyl radical formed from reaction (4) reacts with oxygen to regenerate OH radicals (5a). However, a significant fraction of acetyl radical formed by reaction (R4) is sufficiently energised to dissociate further to CH(3) + CO (R4b). Little or no pressure quenching of reaction (R4b) was observed. The rate coefficient for OD + MGLY was measured as k(9)(T) = (9.4 +/- 2.4) x 10(-13) exp((780 +/- 70)/T) cm(3) molecule(-1) s(-1) over the temperature range 233-500 K. The reaction shows a noticeable inverse (k(H)/k(D) < 1) kinetic isotope effect below room temperature and a slight normal kinetic isotope effect (k(H)/k(D) > 1) at high temperature. The potential atmospheric implications of this work are discussed.

Direct Dynamics Study on the Reaction of N2H4 with F Atom: A Hydrogen Abstraction Reaction?

We present a systematic direct ab initio dynamics investigation of the reaction between N2H4 and F atom, which is predicted to have three possible reaction channels. The structures and frequencies at the stationary points and the points along the minimum energy paths (MEPs) of all reaction channels were calculated at the UB3LYP/6-31+G(d,p) level of theory. Energetic information of stationary points and the points along the MEPs was further refined by means of the CCSD(T)/aug-cc-pVTZ method. The calculated results revealed that the first two primary channels (N2H4+F-->N2H3+HF) are equivalent and occur synchronously via the formation of a pre-reaction complex with Cs symmetry rather than via the direct H abstraction. The pre-reaction complex then evolves into a hydrogen-bonding intermediate through a transition state with nearly no barrier and a high exothermicity, which finally makes the intermediate further decompose into N2H3 and HF. Another reaction channel of minor role (N2H4+F-->NH2F+NH2) was also found during the calculations, which has the same Cs pre-reaction complex but forms NH2F and NH2 via another transition state with high-energy barrier and low exothermicity. The rate constants of these channels were calculated using the improved canonical variational transition state theory with the small-curvature tunneling correction (ICVT/SCT) method. The three-parameter ICVT/SCT rate constant expressions of k(ICVT/SCT) at the CCSD(T)/aug-cc-pVTZ//UB3LYP/6-31+G(d,p) level of theory within 220-3000 K were fitted as (7.64x10(-9))T (-0.87) exp(1180/T) cm3 mole-1 s-1 for N2H4+F-->N2H3+HF and 1.45x10(-12)(T/298)(2.17) exp(-1710/T) cm3 mole-1 s-1 for N2H4+F-->NH2F+NH2.

Reaction-Path Dynamics and Theoretical Rate Constants for the CHnF4-n + O3 → HOOO + CHn-1F4-n (n = 2,3) Reactions

We present a theoretical study of the hydrogen abstraction reactions from CH(3)F and CH(2)F(2) by an ozone molecule. The geometries and harmonic vibrational frequencies of all stationary points are calculated at the MPW1K, BHandHLYP, and MPWB1K levels of theory. The energies of all of the stationary points were refined by using both higher-level (denoted as HL) energy calculations and QCISD(T)/6-311++G(2df,2pd) calculations based on the optimized geometries at the MPW1K/6-31+G(d,p) level of theory. The minimum energy paths (MEPs) were obtained by the MPW1K/6-31+G(d,p) level of theory. Energetic information of the points along the MEPs is further refined by the HL method. The rate constants were evaluated on the basis of the MEPs from the HL level of theory in the temperature range 200-2500 K by using the conventional transition-state theory (TST), the canonical variational transition-state theory (CVT), the microcanonical variational transition-state theory (microVT), the CVT coupled with the small-curvature tunneling (SCT) correction (CVT/SCT), and the microVT coupled with the Eckart tunneling correction (microVT/Eckart) based on the ab initio calculations. A general agreement was found among the TST, CVT, and microVT theories. The fitted three-parameter Arrhenius expressions of the calculated forward CVT/SCT and microVT/Eckart rate constants of the ozonolysis of fluoromethane are k(CVT/SCT)(T) = 2.76 x 10(-34)T(5.81)e((-13975/)(T)) and k(microVT/Eckart)(T) = 1.15 x 10(-34)T(5.97)e((-14530.7/)(T)), respectively. The fitted three-parameter Arrhenius expressions of the calculated forward CVT/SCT and microVT/Eckart rate constants of the ozonolysis of difluoromethane are k(CVT/SCT)(T) = 2.29 x 10(-36)T(6.42)e((-15451.6/)(T)) and k(microVT/Eckart)(T) = 1.31 x 10(-36)T(6.45)e((-15465.8/)(T)), respectively.

Uncovering the Fundamental Chemistry of Alkyl + O2Reactions via Measurements of Product Formation

The reactions of alkyl radicals (R) with molecular oxygen (O(2)) are critical components in chemical models of tropospheric chemistry, hydrocarbon flames, and autoignition phenomena. The fundamental kinetics of the R + O(2) reactions is governed by a rich interplay of elementary physical chemistry processes. At low temperatures and moderate pressures, the reactions form stabilized alkylperoxy radicals (RO(2)), which are key chain carriers in the atmospheric oxidation of hydrocarbons. At higher temperatures, thermal dissociation of the alkylperoxy radicals becomes more rapid and the formation of hydroperoxyl radicals (HO(2)) and the conjugate alkenes begins to dominate the reaction. Internal isomerization of the RO(2) radicals to produce hydroperoxyalkyl radicals, often denoted by QOOH, leads to the production of OH and cyclic ether products. More crucially for combustion chemistry, reactions of the ephemeral QOOH species are also thought to be the key to chain branching in autoignition chemistry. Over the past decade, the understanding of these important reactions has changed greatly. A recognition, arising from classical kinetics experiments but firmly established by recent high-level theoretical studies, that HO(2) elimination occurs directly from an alkylperoxy radical without intervening isomerization has helped resolve tenacious controversies regarding HO(2) formation in these reactions. Second, the importance of including formally direct chemical activation pathways, especially for the formation of products but also for the formation of the QOOH species, in kinetic modeling of R + O(2) chemistry has been demonstrated. In addition, it appears that the crucial rate coefficient for the isomerization of RO(2) radicals to QOOH may be significantly larger than previously thought. These reinterpretations of this class of reactions have been supported by comparison of detailed theoretical calculations to new experimental results that monitor the formation of products of hydrocarbon radical oxidation following a pulsed-photolytic initiation. In this article, these recent experiments are discussed and their contributions to improving general models of alkyl + O(2) reactions are highlighted. Finally, several prospects are discussed for extending the experimental investigations to the pivotal questions of QOOH radical chemistry.

Determination of the temperature and pressure dependence of the reaction OH + C2H4from 200–400 K using experimental and master equation analyses

The pressure and temperature dependence for the reaction of OH + C(2)H(4) was studied over the range of conditions: 200-400 K and 5-600 Torr by laser flash photolysis, laser-induced fluorescence (FP-LIF). Additional experiments were conducted at room temperature by laser flash photolysis, cavity ring-down spectroscopy to facilitate determination of the high pressure limit. One-dimensional master equation calculations were conducted to test the temperature and pressure dependence of the reaction in He and in N(2). The energetics of the reaction and geometries of intermediate species were calculated by ab initio calculations (DFT-BH&HLYP/6-311+G(3df,2p) and CBS-APNO level along DFT-IRC, respectively. An investigation into the importance of a pre-reaction van der Waals complex on the kinetics over the pressure range of the troposphere is discussed. The high pressure rate coefficient was extracted by fitting the master equation calculations to the data and yields k(infinity) = 5.01 x 10(-12) exp(148/T) cm(3) molecule(-1) s(-1). The master equation calculations were then optimized for the pressure fall-off in He and N(2) by varying the average downward energy transfer parameter (DeltaE(down)) for the different collision partners and finally fitted to a Troe expression to determine k(o) and F(cent) for use in atmospheric modeling.

Combustion Chemistry of Enols: Possible Ethenol Precursors in Flames

Before the recent discovery that enols are intermediates in many flames, they appeared in no combustion models. Furthermore, little is known about enols' flame chemistry. Enol formation in low-pressure flames takes place in the preheat zone, and its precursors are most likely fuel species or the early products of fuel decomposition. The OH + ethene reaction has been shown to dominate ethenol production in ethene flames although this reaction has appeared insufficient to describe ethenol formation in all hydrocarbon oxidation systems. In this work, the mole fraction profiles of ethenol in several representative low-pressure flames are correlated with those of possible precursor species as a means for judging likely formation pathways in flames. These correlations and modeling suggest that the reaction of OH with ethene is in fact the dominant source of ethenol in many hydrocarbon flames, and that addition-elimination reactions of OH with other alkenes are also likely to be responsible for enol formation in flames. On this basis, enols are predicted to be minor intermediates in most flames and should be most prevalent in olefinic flames where reactions of the fuel with OH can produce enols directly.

Direct Dynamics Study on the Hydrogen Abstraction Reaction CH2O + HO2 → CHO + H2O2

We present a direct ab initio dynamics study on the hydrogen abstraction reaction CH2O + HO2 --> CHO + H2O2, which is predicted to have four possible reaction channels caused by different attacking orientations of HO2 radical to CH2O. The structures and frequencies at the stationary points and the points along the minimum energy paths (MEPs) of the four reaction channels are calculated at the B3LYP/cc-pVTZ level of theory. Energetic information of stationary points and the points along the MEPs is further refined by means of some single-point multilevel energy calculations (HL). The rate constants of these channels are calculated using the improved canonical variational transition-state theory with the small-curvature tunneling correction (ICVT/SCT) method. The calculated results show that, in the whole temperature range, the more favorable reaction channels are Channels 1 and 3. The total ICVT/SCT rate constants of the four channels at the HL//B3LYP/cc-pVTZ level of theory are in good agreement with the available experiment data over the measured temperature ranges, and the corresponding three-parameter expression is k(ICVT/SCT) = 3.13 x 10(-20) T(2.70) exp(-11.52/RT) cm3 mole(-1) s(-1) in the temperature range of 250-3000 K. Additionally, the flexibility of the dihedral angle of H2O2 is also discussed to explain the different experimental values.

The reaction of the aminoboranylidene-iminoborane isomerization: a CASSCF direct dynamics study

The barrier and the potential-energy surface of the isomerization from aminoboranylidene (BNH2) to iminoborane (HBNH) have been studied using complete active space self-consistent field (CASSCF) with the 6-31 + G(d, p) basis set and higher-level energy methods. The rate constants of the isomerization reaction are reported by employing the direct ab initio dynamics method. The geometries of all the stationary points were optimized using the B3LYP and CCSD methods with the cc-pVTZ and cc-pVQZ basis sets. The information along the intrinsic reaction coordinate (IRC) was also calculated at the CASSCF/6-31 + G (d,p) level of theory. The energies were refined at the G3, G3MP2, G3MP2B3, CBS-Q, CBS-QB3, and two high-level (HL) methods based on the geometries optimized using CASSCF/6-31 + G(d,p). The rate constants were evaluated using conventional transition-state theory (TST), canonical variational transition-state theory (CVT), and canonical variational transition-state theory with small curvature tunneling correction (CVT/SCT) and conventional transition-state theory with Eckart tunneling correction (TST/Eckart). According to the calculated results, we conclude that the tunneling effect is very important to this isomerization reaction.

Isomeric Product Distributions from the Self-Reaction of Propargyl Radicals

We have investigated the isomeric C6H6 product distributions of the self-reaction of propargyl (C3H3) radicals at two nominal pressures of 25 and 50 bar over the temperature range 720-1350 K. Experiments were performed using propargyl iodide as the radical precursor in a high-pressure single-pulse shock tube with a residence time of 1.6-2.0 ms. The relative yields of the C6H6 products are strongly temperature dependent, and the main products are 1,5-hexadiyne (15HD), 1,2-hexadiene-5-yne (12HD5Y), 3,4-dimethylenecyclobutene (34DMCB), 2-ethynyl-1,3-butadiene (2E13BD), fulvene, and benzene, with the minor products being cis- and trans-1,3-hexadiene-5-yne (13HD5Y). 1,2,4,5-Hexatetraene (1245HT) was observed below 750 K but the concentrations were too low to be quantified. The experimentally determined entry branching ratios are: 44% 15HD, 38% 12HD5Y, and 18% 1245HT, which is efficiently converted to 34DMCB. Following the initial recombination step, various C6H6 isomers are formed by thermal rearrangement. The experimentally observed concentrations for the C6H6 species are in good agreement with earlier experiments on 15HD thermal rearrangement.

Global Mapping of Equilibrium and Transition Structures on Potential Energy Surfaces by the Scaled Hypersphere Search Method: Applications to ab Initio Surfaces of Formaldehyde and Propyne Molecules

Technical details of a new global mapping technique for finding equilibrium (EQ) and transition structures (TS) on potential energy surfaces (PES), the scaled hypersphere search (SHS) method (Ohno, K.; Maeda, S. Chem. Phys. Lett. 2004, 384, 277), are presented. On the basis of a simple principle that reaction pathways are found as anharmonic downward distortions of PES around an EQ point, the reaction pathways can be obtained as energy minima on the scaled hypersphere surface, which would have a constant energy when the potentials are harmonic. Connections of SHS paths between each EQ are very similar to corresponding intrinsic reaction coordinate (IRC) connections. The energy maximum along the SHS path reaches a region in close proximity to the TS of the reaction pathway, and the subsequent geometry optimization from the SHS maximum structure easily converges to the TS. The SHS method, using the one-after-another algorithm connecting EQ and TS, considerably reduces the multidimensional space to be searched to certain limited regions around the pathways connecting each EQ with the neighboring TS. Applications of the SHS method have been made to ab initio surfaces of formaldehyde and propyne molecules to obtain systematically five EQ and nine TS for formaldehyde and seven EQ and 32 TS for propyne.

Shock Wave Study on the Thermal Unimolecular Decomposition of Allyl Radicals

We present the first direct study on the thermal unimolecular decomposition of allyl radicals. Experiments have been performed behind shock waves, and the experimental conditions covered temperatures ranging from 1125 K up to 1570 K and pressures between 0.25 and 4.5 bar. Allyl radicals have been generated by thermal decomposition of allyl iodide, and H-atom resonance absorption spectroscopy has been used to monitor the reaction progress. A marked pressure dependence of the rate constant has been observed which is in agreement with the results from a master equation analysis. However, our experimental results as well as our Rice-Ramsperger-Kassel-Marcus calculations seem to contradict the results of Deyerl et al. (J. Chem. Phys. 1999, 110, 1450) who investigated the unimolecular decomposition of allyl radicals upon photoexcitation and tried to deduce specific rate constants for the unimolecular dissociation in the electronic ground state. At pressures around 1 bar we extracted the following rate equation: k(T) = 5.3 x 10(79)(T/K)(-19.29) exp[(-398.9 kJ/mol)/RT] s(-1). The uncertainty of the rate constant calculated from this equation is estimated to be 30%.

Thermochemistry of disputed soot formation intermediates C4H3 and C4H5

Accurate isomeric energy differences and standard enthalpies of formation for disputed intermediates in soot formation, C(4)H(3) and C(4)H(5), have been determined through systematic extrapolations of ab initio energies. Electron correlation has been included through second-order Z-averaged perturbation theory (ZAPT2), and spin-restricted, open-shell coupled-cluster methods through triple excitations [ROCCSD, ROCCSD(T), and ROCCSDT] utilizing the correlation-consistent hierarchy of basis sets, cc-pVXZ (X = D, T, Q, 5, and 6), followed by extrapolations to the complete basis set limit via the focal point method of Allen and co-workers. Reference geometries were fully optimized at the ROCCSD(T) level with a TZ(2d1f,2p1d) basis set. Our analysis finds that the resonance-stabilized i-C(4)H(3) and i-C(4)H(5) isomers lie 11.8 and 10.7 kcal mol(-1) below E-n-C(4)H(3) and E-n-C(4)H(5), respectively, several kcal mol(-1) (more, less) than reported in recent (diffusion Monte Carlo, B3LYP density-functional) studies. Moreover, in these systems Gaussian-3 (G3) theory suffers from large spin contamination in electronic wave functions, poor reference geometries, and anomalous vibrational frequencies, but fortuitous cancellation of these sizable errors leads to isomerization energies apparently accurate to 1 kcal mol(-1). Using focal-point extrapolations for isodesmic reactions, we determine the enthalpies of formation (delta(f)H(0) (composite function)) for i-C(4)H(3), Z-n-C(4)H(3), E-n-C(4)H(3), i-C(4)H(5), Z-n-C(4)H(5), and E-n-C(4)H(5) to be 119.0, 130.8, 130.8, 78.4, 89.7, and 89.1 kcal mol(-1), respectively. These definitive values remove any remaining uncertainty surrounding the thermochemistry of these isomers in combustion models, allowing for better assessment of whether even-carbon pathways contribute to soot formation.