Crossed-beam dynamics studies of the radical-radical combustion reaction O((3)P) + CH3 (methyl)

The dynamics of the radical-radical reaction O((3)P) + CH(3), a prototypical case for the reactions of atomic oxygen with alkyl radicals of great relevance in combustion chemistry, has been investigated by means of the crossed molecular beam technique with mass spectrometric detection at a collision energy of 55.9 kJ mol(-1). The results have been examined in the light of previous kinetic and theoretical work. From product angular and velocity distribution measurements, the dynamics of the predominant H-displacement channel leading to formaldehyde formation has been characterized. This channel has been found to proceed via the formation of an osculating complex; a significant coupling between the product centre-of-mass angular and translational energy distributions has been noted. Experimental attempts to characterize the dynamics of the channel leading to HCO + H(2) have failed and it remains unclear whether HCO is formed by the reaction and/or, if formed, a part of HCO does not dissociate quickly into CO + H.     

Seven-dimensional quantum dynamics study of the O(3P)+CH4 reaction

The initial state selected time-dependent wave packet calculations have been carried out to study the title reaction with seven degrees of freedom included by restricting the nonreacting CH(3) group under C(3V) symmetry and the CH bond length in the group. Total reaction probabilities as well as integral cross sections were calculated for the ground and four vibrationally excited reagent states. Our calculation shows that the reactivity is very small for the reaction for collision energy up to 1.0 eV for all the initial states. Initial vibration excitation of CH(4), in particular, the CH stretch excitation, enhances the reactivity, but only part of the excitation energy deposited can be used to reduce the reaction threshold. The rate constant for the ground initial state agrees rather well with that from a recent quasiclassical trajectory study and is larger than that from the semirigid vibrating rotor target calculations, in particular, in the low temperature region. On the other hand, the thermal rate constant calculated from the integral cross sections for these five vibrational states is about a factor of 20 smaller than that from the multiconfiguration time-dependent Hartree study.     

Collision dynamics of O(3P) + DMMP using a specific reaction parameters potential form

Starting from previous benchmark CBS-QB3 electronic structure calculations (Conforti, P. F.; Braunstein, M.; Dodd, J. A. J. Phys. Chem. A 2009, 113, 13752), we develop two global potential energy surfaces for O((3)P) + DMMP collisions, using the specific reaction parameters approach. Each surface is simultaneously fit along the three major reaction pathways: hydrogen abstraction, hydrogen elimination, and methyl elimination. We then use these surfaces in classical dynamics simulations and compute reactive cross sections from 4 to 10 km s(-1) collision velocity. We examine the energy disposal and angular distributions of the reactive and nonreactive products. We find that for reactive collisions, an unusually large amount of the initial collision energy is transformed into internal energy. We analyze the nonreactive and reactive product internal energy distributions, many of which fit Boltzmann temperatures up to ~2000 K.     

Trajectory surface hopping molecular dynamics on Chemiluminescence of cyclic peroxides

Chemiluminescence (CL) of peroxides is one of the most highly sensitive and most useful analytical techniques. Although the mechanisms of CL were studied experimentally and theoretically in the past decades, the chemiexcitation that a ground-state specie being excited from its electronic ground state to yield excited-state products by a chemical reaction is still not completely understood. Direct dynamics simulation on CL reaction which takes into account the full complexity of the relevant potential energy surface characterizes nonadiabatic processes involved in chemiexcitation and could provide access not only to the available reaction channels but also to statistical quantities such as reaction times and quantum yields. In the last decade, the trajectory surface hopping (TSH) molecular dynamics (MD) which is one of the mixed quantum-classical approaches is hence adopted to simulate the nonadiabatic process in CL of cyclic peroxides. In this article, the basic principle of TSH-MD and the successful applications on the CL reactions were shortly review.     

Quantum chemical and theoretical kinetics study of the O(3P) + CS2 reaction

The triplet potential energy surface of the O((3)P) + CS(2) reaction is investigated by using various quantum chemical methods including CCSD(T), QCISD(T), CCSD, QCISD, G3B3, MPWB1K, BB1K, MP2, and B3LYP. The thermal rate coefficients for the formation of three major products, CS + SO ((3)Σ(-)), OCS + S ((3)P) and CO + S(2) ((3)Σ(-)(g)) were computed by using transition state and RRKM statistical rate theories over the temperature range of 200-2000 K. The computed k(SO + CS) by using high-level quantum chemical methods is in accordance with the available experimental data. The calculated rate coefficients for the formation of OCS + S ((3)P) and CO + S(2) ((3)Σ(-)(g)) are much lower than k(SO + CS); hence, it is predicted that these two product channels do not contribute significantly to the overall rate coefficient.     

Reaction dynamics of CN + O2 --> NCO + O(3P2)

We have used oxygen Rydberg time-of-flight spectroscopy to carry out a crossed molecular beam study of the CN + O2 reaction at collision energies of 3.1 and 4.1 kcal/mol. The O(3P2) products were tagged by excitation to high-n Rydberg levels and subsequently field ionized at a detector. The translational energy distributions were broad, indicating that the NCO is formed with a wide range of internal excitation, and the angular distribution was forward-backward symmetric, indicating the participation of NCOO intermediates with lifetimes comparable to or longer than their rotational periods. Rice-Ramsperger-Kassel-Marcus modeling of the dissociation of NCOO to NCO + O suggests that Do(NC-OO) > or = 38 kcal/mol, which is consistent with several theoretical calculations. Implications for the competing CO + NO channel are discussed.     

The effects of surface temperature on the gas-liquid interfacial reaction dynamics of O(3P)+squalane

OH/OD product state distributions arising from the reaction of gas-phase O(3P) atoms at the surface of the liquid hydrocarbon squalane C30H62/C30D62 have been measured. The O(3P) atoms were generated by 355 nm laser photolysis of NO2 at a low pressure above the continually refreshed liquid. It has been shown unambiguously that the hydroxyl radicals detected by laser-induced fluorescence originate from the squalane surface. The gas-phase OH/OD rotational populations are found to be partially sensitive to the liquid temperature, but do not adapt to it completely. In addition, rotational temperatures for OH/OD(v'=1) are consistently colder (by 34+/-5 K) than those for OH/OD(v'=0). This is reminiscent of, but less pronounced than, a similar effect in the well-studied homogeneous gas-phase reaction of O(3P) with smaller hydrocarbons. We conclude that the rotational distributions are composed of two different components. One originates from a direct abstraction mechanism with product characteristics similar to those in the gas phase. The other is a trapping-desorption process yielding a thermal, Boltzmann-like distribution close to the surface temperature. This conclusion is consistent with that reached previously from independent measurements of OH product velocity distributions in complementary molecular-beam scattering experiments. It is further supported by the temporal profiles of OH/OD laser-induced fluorescence signals as a function of distance from the surface observed in the current experiments. The vibrational branching ratios for (v'=1)/(v'=0) for OH and OD have been found to be (0.07+/-0.02) and (0.30+/-0.10), respectively. The detection of vibrationally excited hydroxyl radicals suggests that secondary and/or tertiary hydrogen atoms may be accessible to the attacking oxygen atoms.     

The mechanism of O(3P) atom reaction with ethylene and other simple olefins

Reactions of oxygen atoms with ethylene, propene, and 2-butene were studied at room temperature under discharge flow conditions by resonance fluorescence spectroscopy of O and H atoms at pressures of 0.08 to 12 torr. The measured total rate constants of these reactions are K = (7.8 ± 0.6)·10^?13cm³s^?1,K = (4.3 ± 0.4) ± 10^?12 cm³ s^?1, K = (1.4 ± 0.4) · 10^?11 cm³ s^?1. The branching ratios of H atom elimination channels were measured for reactions of O atoms with ethylene and propene. No H-atom elimination was found for the reaction of O-atoms with 2-butene. A redistribution of reaction O + C₂ channels with pressure was found. A mechanism of the O + C₂ reaction was proposed and the possibility of its application to other olefins is discussed. On the basis of mechanism the pressure dependence of the total rate constant for reaction O + C₂ was predicted and experimentally confirmed in the pressure range 0.08–1.46 torr.     

Ab initio investigations of the radical-radical reaction of O(3P) + C3H3

We present ab initio calculations of the reaction of ground-state atomic oxygen [O((3)P)] with a propargyl (C(3)H(3)) radical based on the application of the density-functional method and the complete basis-set model. It has been predicted that the barrierless addition of O((3)P) to C(3)H(3) on the lowest doublet potential-energy surface produces several energy-rich intermediates, which undergo subsequent isomerization and decomposition steps to generate various exothermic reaction products: C(2)H(3)+CO, C(3)H(2)O+H, C(3)H(2)+OH, C(2)H(2)+CHO, C(2)H(2)O+CH, C(2)HO+CH(2), and CH(2)O+C(2)H. The respective reaction pathways are examined extensively with the aid of statistical Rice-Ramsperger-Kassel-Marcus calculations, suggesting that the primary reaction channel is the formation of propynal (CHCCHO)+H. For the minor C(3)H(2)+OH channel, which has been reported in recent gas-phase crossed-beam experiments [H. Lee et al., J. Chem. Phys. 119, 9337 (2003); 120, 2215 (2004)], a comparison on the basis of prior statistical calculations is made with the nascent rotational state distributions of the OH products to elucidate the mechanistic and dynamic characteristics at the molecular level.     

Quantum chemical study of the reaction of trichloroethylene with O(3P)

The adiabatic mechanism of the reaction of trichloroethylene with O(³P), exploring the various O-atom addition and H-atom abstraction channels, is theoretically studied at the MP2/6-311++G(2d, 2p), MP2/aug-cc-pVTZ, CCSD/6-31G(d), G3, and CBS-QB3 levels of theory. From a kinetic point of view, the addition to the less substituted carbon atom of the double bond is more favorable than the addition to the more substituted carbon. Such O-atom addition reactions are favored over the one possible hydrogen-abstraction reaction. Calculations of the present study showed that five products are obtained: HCCl + C(O)Cl₂ (P1), Cl + ClC(O)CHCl (P2), H + ClC(O)CCl₂ (P3), Cl + HC(O)CCl₂ (P4), and CH(O)Cl + CCl₂ (P5). The products P2 and P4 are found to be the most favored ones. The kinetic calculations of rate constant in the range of 285–395 K are performed at the CBS-QB3 level of theory and are in conformity with the experimental outcomes.     

A theoretical study of the reaction of O(3P) with isobutene

The reaction of O((3)P) with isobutene ((CH(3))(2)C=CH(2)) is investigated using the unrestricted second-order M?ller-Plesset perturbation (UMP2) and complete basis set CBS-4M level methods. The minimum energy crossing point (MECP) between the singlet and triplet potential energy surfaces is located using the Newton-Lagrange method, and it is shown that the MECP plays a key role. The calculational results indicate that the site selectivity of the addition of O((3)P) to either carbon atom of the double bond of isobutene is weak, and the major product channels are CH(2)C(O)CH(3) + CH(3,) cis-/trans-CH(3)CHCHO + CH(3), (CH(3))(2)CCO + H(2), and CH(3)C(CH(2))(2) + OH, among which (CH(3))(2)CCO + H(2) is predicted to be the energetically most favorable one. The complex multichannel reaction mechanisms are revealed, and the observations in several recent experiments could be rationalized on the basis of the present calculations. The formation mechanisms of butenols are also discussed.     

A kinetic study of the reaction of O(3P) with toluene

The absolute rate constant of the reaction of O(³P) with toluene was measured by the microwave discharge-fast-flow method to obtain k_T = 10^{9.7–2.7/2.303RT} l./mole·sec at 100–375°C. This was in good agreement with the rate constant calculated from the combination of the relative rate constant obtained by Jones and Cvetanovic with the recently determined absolute rate constants of the reaction of O(³P) + olefins. The extrapolation of the above Arrhenius plot to 27°C was also in good agreement with the absolute value of k_T = 4.5 × 10⁷ l./mole·sec determined recently by Atkinson and Pitts at 27°C. The rate constant of the reaction of chlorobenzene with O(³P), obtained at 238°C as 10^8.3 l./mole·sec by a competitive method, was smaller than k_T by a factor of about two at the same temperature.     

Quasiclassical trajectory study of the O(3P) + CH4 --> OH + CH3 reaction with a specific reaction parameters semiempirical Hamiltonian

We present a theoretical study of the O(3P) + CH4 --> OH + CH3 reaction using electronic structure, kinetics, and dynamics calculations. We calculate a grid of ab initio points at the PMP2/AUG-cc-pVDZ level to characterize the potential energy surface in regions of up to 1.3 eV above reagents. This grid of ab initio points is used to derive a set of specific reaction parameters (SRP) for the MSINDO semiempirical Hamiltonian. The resulting SRP-MSINDO Hamiltonian improves the quality of the standard Hamiltonian, particularly in regions of the potential energy surface beyond the minimum-energy reaction path. Quasiclassical-trajectory calculations are used to study the reaction dynamics with the original and the improved MSINDO semiempirical Hamiltonians, and a prior surface. The SRP-MSINDO semiempirical Hamiltonian yields OH rotational distributions in agreement with experimental results, improving over the results of the other surfaces. Thermal rate constants estimated with Variational Transition State Theory using the SRP-MSINDO Hamiltonian are also in agreement with experiments. Our results indicate that reparametrized semiempirical Hamiltonians are a good alternative to generating potential energy surfaces for accurate dynamics studies of polyatomic reactions.     

DFT investigation of the mechanism of CH2CO + O(3P) reaction

A detailed theoretical survey of the potential energy surface (PES) for the CH₂CO + O(³P) reaction is carried out at the QCISD(T)/6‐311+G(3df,2p)//B3LYP/6‐311+G(d,p) level. The geometries, vibrational frequencies, and energies of all stationary points involved in the reaction are calculated at the B3LYP/6‐311+G(d,p) level. More accurate energy information is provided by single‐point calculations at the QCISD(T)/6‐311+G(3df,2p) level. Relationships of the reactants, transition states, intermediates, and products are confirmed by the intrinsic reaction coordinate (IRC) calculations. The results suggest that P1(CH₂+CO₂) is the most important product. This study presents highlights of the mechanism of the title reaction. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Direct dynamics simulations of O(3P) + HCl at hyperthermal collision energies

The dynamics of the O(3P) + HCl reaction at hyperthermal collision energies were investigated using the quasiclassical trajectory method. Stationary points on the OClH 3A" and 3A' potential energy surfaces (PESs) were also examined. The lowest transition state leading to OCl + H on the 3A" surface is 2.26 eV above the reagents at the CCSD(T)/cc-pVTZ level of theory. This saddle point is bent and product-like. Direct dynamics calculations at the MP2/cc-pVTZ level of theory were used to investigate the excitation functions for OH + Cl, OCl + H, and O + H + Cl formation. OCl is formed mainly from small-impact-parameter collisions, and the OCl + H excitation function peaks around 5 eV, where it is similar in magnitude to the OH + Cl excitation function. The shape of the OCl + H excitation function is discussed, and features are identified that should be general to hyperthermal collision dynamics.     

Direct ab initio dynamics study on the rate constants and kinetics isotope effects of CH(3)O+H-->CH(2)O+H(2) reaction

We present a direct ab initio dynamics study of thermal rate constants of the hydrogen abstraction reaction of CH(3)O+H-->CH(2)O+H(2). The unrestricted Becke's half-and-half hybrid functional using the Lee-Yang-Parr correlation functional with Dunning's correlation consistent polarized valence double-zeta basis set, the unrestricted quadratic configuration interaction calculation including single and double substitutions with Dunning's correlation consistent polarized valence double-zeta basis set, and the unrestricted quadratic configuration interaction calculation including single and double substitutions with a triples contribution with Dunning's correlation consistent polarized valence triple-zeta basis set methods were employed to optimize the structures and to calculate frequencies for all stationary points. Minimum energy paths were obtained by the unrestricted Becke's half-and-half hybrid functional using the Lee-Yang-Parr correlation functional and the unrestricted quadratic configuration interaction calculation including single and double substitutions with the same Dunning's correlation consistent polarized valence double-zeta basis set levels of theory. No barrier is found at the unrestricted Becke's half-and-half hybrid functional using the Lee-Yang-Parr correlation functional with Dunning's correlation consistent polarized valence double-zeta basis set level of theory in contrast to a small barrier of 1.43 kcal mol(-1) at the unrestricted quadratic configuration interaction calculation including single and double substitutions with Dunning's correlation consistent polarized valence double-zeta basis set level of theory. In particular, the barrier vanishes as the energies along the minimum energy path MEP are refined at the unrestricted quadratic configuration interaction calculation including single and double substitutions with a triples contribution with Dunning's correlation consistent polarized valence triple-zeta basis set level of theory. Smaller barriers of 0.47 and 0.17 kcal mol(-1) were obtained at the unrestricted quadratic configuration interaction calculation including single and double substitutions with a triples contribution with Dunning's correlation consistent polarized valence triple-zeta basis set and the unrestricted quadratic configuration interaction calculation including single and double substitutions with a triples contribution with Dunning's correlation consistent polarized valence triple-zeta basis set based on the geometries at the unrestricted quadratic configuration interaction calculation including single and double substitutions with Dunning's correlation consistent polarized valence triple-zeta basis set levels of theory, respectively. The forward rate constants are evaluated with the canonical variational transition state theory in the temperature range of 300-2500 K. The calculated forward rate constants at the unrestricted quadratic configuration interaction calculation including single and double substitutions with a triples contribution with Dunning's correlation consistent polarized valence triple-zeta basis set based on the geometries at the unrestricted quadratic configuration interaction calculation including single and double substitutions with Dunning's correlation consistent polarized valence double-zeta basis set level of theory are in good agreement with the available experimental data. The kinetic isotope effects are estimated.     

Direct gas-liquid interfacial dynamics: the reaction between O(3P) and a liquid hydrocarbon

We report the first measurements of internal energy distributions of the OH produced via a direct mechanism, isolated from other components on the basis of time-of-flight, in the interfacial reaction between gas-phase O((3)P) atoms and the liquid hydrocarbon squalane, C(30)H(62). O((3)P) atoms were generated by laser photolysis of NO(2) above the liquid. Resulting hydroxyl radicals that escape from the surface were detected by laser-induced fluorescence. Time-of-flight profiles demonstrate that the kinetic energy of the fastest OH (nu' = 1) is lower than that of (nu' = 0). Rotational distributions were measured at the rising edge of their appearance for both OH (nu' = 0) and (nu' = 1). They were found to differ substantially more than at the peak of their profiles. They were also less dependent on the bulk liquid temperature. We conclude that the new data confirm strongly that at least two mechanisms contribute to the production of OH. The higher-velocity component has translational and rotational energy distributions, observed cleanly for the first time, consistent with a direct mechanism. The close correspondence of these rotational distributions to those from the corresponding homogeneous gas-phase reaction of O((3)P) with smaller hydrocarbons suggests a very similar, near collinear direct abstraction. This is accompanied by a slower component with kinetic energy and rotational (but not vibrational) distributions reflecting the temperature of the liquid, consistent with a distinct trapping-desorption mechanism.     

Br+CH3S(O)CH3反应机理的直接动力学研究

采用直接动力学的方法，对多通道反应体系Br＋CH3S（O）CH3进行了理论研究．在BH＆H-LYP／6-311G（2d，2p）水平下获得了优化几何构型、频率及最小能量路径（MEP），能量信息的进一步确认在MC-QCISD（单点）水平下完成．利用正则变分过渡态理论，结合小曲率隧道效应校正（CVT／SCT）方法计算了该反应的两个可行的反应通道在200K～2000K温度范围内的速率常数．在整个反应区间内，生成HBr的反应通道与生成CHa的反应通道存在着竞争，前者是主反应通道，后者是次反应通道．变分效应和小曲率隧道效应对反应速率常数的计算影响都很小．理论计算得到的两个反应通道的反应速率常数与实验值符合得很好．     

The reaction of O(3P) with OClO

The laser photolysis-long path transient absorption technique was used to study the mechanism and products of the reaction O(³P) + OClO (→^M) Products (3) at 260 K in 100 to 400 torr of He, N₂, and Ar. ClO was not detected as a reaction product, <5% yield, from this reaction at 400 torr and 260 K. A broad UV absorption feature associated with a product of this reaction, with a peak located at 260 nm, was observed. The peak absorption cross section of this species was measured to be (1.72 ± 0.12) × 10⁻¹⁷ cm² molecule⁻¹. The rate coefficient for the appearance of this species was measured to be (1.69 ± 0.46) × 10⁻¹³ cm³ molecule⁻¹ s⁻¹ and was independent of pressure. The rate coefficient for the appearance of the species is ca. 10 times lower than that for the disappearance of O(³P), indicating that the observed species is not a direct product of the reaction of O(³P) with OClO. Mechanistic considerations and the possible identity of the absorber are discussed. © 1997 John Wiley & Sons, Inc.     

A computational investigation on the potential energy surface of thiosulfeno with O(3P) reaction

The reaction paths of thiosulfeno radical (HS₂) with O(³P) have been investigated at the UB3LYP/aug-cc-pV(T + d)Z and UCCSD(T)/aug-cc-pV(T + d)Z//B3LYP levels. Two stable collision intermediates, HSSO and SS(H)O, have been considered for the HS₂ + O(³P) reaction. Four products of S + HSO, H + SSO, HS + SO, and S₂ + OH are obtained by starting from HSSO and SS(H)O. The calculated results show that the most feasible paths for the formation of S + HSO, H + SSO, and HS + SO products include no transition states in reaction path, while that of S₂ + OH product includes relatively high energy barriers of 23.0 kcal/mol. Therefore, S + HSO, H + SSO, and HS + SO are main products (with the stability other of HS + SO > H + SSO > S + HSO) and S₂ + OH is the second product in HS₂ + O(³P) reaction. Because, all intermediates, transition states, and products involved in the reaction paths lie below the initial reactants, the HS₂ + O(³P) reaction is expected to be rapid even at low temperatures.