CH_3SGN与O_2气相反应机理的理论研究

在G3(MP2)水平上,通过对CH_3S与O_2rcyi2rvylce dm (PES)上关键驻点的能 量计算,共找到4种中间体,9个过渡态,6种产物通道,并对这些气相反应机理进 行了讨论,同时应用TST－RRKM理论对主要反应的速率进行计算。结果表明:CH_3S 与O_2反应在低温下以生成CH_3SOO为主,并与实验结果吻合;在中高温下以消去和 抽提反应为主,分别生成CH_3 + SO_2和CH_2S + HO_2,其它产物较少。     

The effects of water vapor on the CH3O2 self-reaction and reaction with HO2

The gas phase reactions of CH3O2 + CH3O2, HO2 + HO2, and CH3O2 + HO2 in the presence of water vapor have been studied at temperatures between 263 and 303 K using laser flash photolysis coupled with UV time-resolved absorption detection at 220 and 260 nm. Water vapor concentrations were quantified using tunable diode laser spectroscopy operating in the mid-IR. The HO2 self-reaction rate constant is significantly enhanced by water vapor, consistent with what others have reported, whereas the CH3O2 self-reaction and the cross-reaction (CH3O2 + HO2) rate constants are nearly unaffected. The enhancement in the HO2 self-reaction rate coefficient occurs because of the formation of a strongly bound (6.9 kcal mol(-1)) HO2 x H2O complex during the reaction mechanism where the H2O acts as an energy chaperone. The nominal impact of water vapor on the CH3O2 self-reaction rate coefficient is consistent with recent high level ab initio calculations that predict a weakly bound CH3O2 x H2O complex (2.3 kcal mol(-1)). The smaller binding energy of the CH3O2 x H2O complex does not favor its formation and consequent participation in the methyl peroxy self-reaction mechanism.     

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.     

Electronic structure study of the initiation routes of the dimethyl sulfide oxidation by OH

In the present work the potential energy surface (PES) corresponding to the different initiation routes of the oxidation mechanism of DMS by hydroxyl radical in the absence of O(2) has been studied, and connections among the different stationary points have been established. Single-point high level electronic structure calculations at lower level optimized geometries have been shown to be necessary to assure convergence of energy barriers and reaction energies. Our results demonstrate that the oxidation of DMS by OH turns out to be initiated via three channels: a hydrogen abstraction channel that through a saddle point structure finally leads to CH(3)SCH(2) + H(2)O, an addition-elimination channel that firstly leads to an adduct complex (AD) and then via an elimination saddle point structure finally gives CH(3)SOH and CH(3) products, and a third channel that through a concerted pathway leads to CH(3)OH and CH(3)S. The H-abstraction and the addition-elimination channels initiate by a common pathway that goes through the same reactant complex (RC). Our theoretical results agree quite well with the branching ratios experimentally assigned to the formation of the different products. Finally, the calculated equilibrium constants of the formation of the complex AD and the hexadeuterated complex AD from the corresponding reactants, as a function of the temperature, are in good accordance with the experimental values.     

Direct ab initio molecular dynamics study on a microsolvated SN2 reaction of OH-(H2O) with CH3Cl

Reaction dynamics for a microsolvated SN2 reaction OH-(H2O)+CH3Cl have been investigated by means of the direct ab initio molecular dynamics method. The relative center-of-mass collision energies were chosen as 10, 15, and 25 kcal/mol. Three reaction channels were found as products. These are (1) a channel leading to complete dissociation (the products are CH3OH+Cl- +H2O: denoted by channel I), (2) a solvation channel (the products are Cl-(H2O)+CH3OH: channel II), and (3) a complex formation channel (the products are CH3OH...H2O+Cl-: channel III). The branching ratios for the three channels were drastically changed as a function of center-of-mass collision energy. The ratio of complete dissociation channel (channel I) increased with increasing collision energy, whereas that of channel III decreased. The solvation channel (channel II) was minor at all collision energies. The selectivity of the reaction channels and the mechanism are discussed on the basis of the theoretical results.     

Ab initio study of the F + CH3NHNH2 reaction mechanism

The F + CH(3)NHNH(2) reaction mechanism is studied based on ab initio quantum chemistry methods as follows: the minimum energy paths (MEPs) are computed at the UMP2/6-311++G(d,p) level; the geometries, harmonic vibrational frequencies, and energies of all stationary points are predicted at the same level of theory; further, the energies of stationary points and the points along the MEPs are refined by UCCSD(T)/6-311++g(3df,2p). The ab initio study shows that, when the F atom approaches CH(3)NHNH(2), the heavy atoms, namely N and C atoms, are the favorable combining points. For the two N atoms, two prereaction complexes with C(s) symmetry are generated and there exists seven possible subsequent reaction routes, of which routes 1, 2, 5, and 7 are the main channels. Routes 1, 2, and 5 are associated with HF elimination, with H from the amino group or imido group, and route 7 involves the N-N bond break. Routes 3 and 6 with relation to HF elimination with H from methyl, and route 4 involved the C-N bond break, are all energetically disfavored. For the C atom, the attack of F results in the break of the C-N bond and the products are CH(3)F + NHNH(2). This route is very competitive.     

Photodissociation dynamics of the ethoxy radical investigated by photofragment coincidence imaging

The photodissociation dynamics of the ethoxy radical (CH3CH2O) have been studied at energies from 5.17 to 5.96 eV using photofragment coincidence imaging. The upper state of the electronic transition excited at these energies is assigned to the C2A'state on the basis of electronic structure calculations. Fragment mass distributions show two photodissociation channels, OH + C2H4 and CH3 + CH2O. The presence of an additional photodissociation channel, identified as D + C2D4O, is revealed in time-of-flight distributions from the photodissociation of CD3CD2O. The product branching ratios and fragment translational energy distributions for all of the observed mass channels are nonstatistical. Moreover, the significant yield of OH + C2H4 product suggests that the mechanism for this channel involves isomerization on the excited-state surface. Photodissociation at a much lower yield is seen following excitation at 3.91 eV, corresponding to a vibronic band of the B2A' <-- X2A' transition.     

CH3CH2+O(3P)反应机理的理论研究

The reaction for CH3CH2+O(3P) was studied by ab initio method. The geometries of the reactants, intermediates, transition states and products were optimized at MP2/6-311+G(d,p) level. The corresponding vibration frequencies were calculated at the same level. The single-point calculations for all the stationary points were carried out at the QCISD(T)/6-311+G(d,p) level using the MP2/6-311+G(d,p) optimized geometries. The results of the theoretical study indicate that the major products are the CH2O＋CH3, CH3CHO+H and CH2CH2+OH in the reaction. For the products CH2O＋CH3 and CH3CHO+H, the major production channels are A1: (R)→IM1→TS3→(A) and B1: (R)→IM1→TS4→(B), respectively. The majority of the products CH2CH2+OH are formed via the direct abstraction channels C1 and C2: (R)→TS1(TS2)→(C). In addition, the results suggest that the barrier heights to form the CO reaction channels are very high, so the CO is not a major product in the reaction.     

Overtone-induced dissociation and isomerization dynamics of the hydroxymethyl radical (CH2OH and CD2OH). I. A theoretical study

The dissociation of the hydroxymethyl radical, CH(2)OH, and its isotopolog, CD(2)OH, following the excitation of high OH stretch overtones is studied by quasi-classical molecular dynamics calculations using a global potential energy surface (PES) fitted to ab initio calculations. The PES includes CH(2)OH and CH(3)O minima, dissociation products, and all relevant barriers. Its analysis shows that the transition states for OH bond fission and isomerization are both very close in energy to the excited vibrational levels reached in recent experiments and involve significant geometry changes relative to the CH(2)OH equilibrium structure. The energies of key stationary points are refined using high-level electronic structure calculations. Vibrational energies and wavefunctions are computed by coupled anharmonic vibrational calculations. They show that high OH-stretch overtones are mixed with other modes. Consequently, trajectory calculations carried out at energies about ~3000 cm(-1) above the barriers reveal that despite initial excitation of the OH stretch, the direct OH bond fission is relatively slow (10 ps) and a considerable fraction of the radicals undergoes isomerization to the methoxy radical. The computed dissociation energies are: D(0)(CH(2)OH → CH(2)O + H) = 10,188 cm(-1), D(0)(CD(2)OH → CD(2)O + H) = 10,167 cm(-1), D(0)(CD(2)OH → CHDO + D) = 10,787 cm(-1). All are in excellent agreement with the experimental results. For CH(2)OH, the barriers for the direct OH bond fission and isomerization are: 14,205 and 13,839 cm(-1), respectively.     

Direct dynamics classical trajectory simulations of the O+ + CH4 reaction at hyperthermal energies

A Born-Oppenheimer direct dynamics simulation of the O(+) + CH(4) reaction dynamics at hyperthermal energies has been carried out with the PM3 (ground quartet state) Hamiltonian. Calculations were performed at various collision energies ranging from 0.5 to 10 eV with emphasis on high energy collisions where this reaction is relevant to materials erosion studies in low Earth orbit and geosynchronous Earth orbit. Charge transfer to give CH(4)(+) is the dominant channel arising from O(+) + CH(4) collisions in this energy range, but most of the emphasis in our study is on collisions that lead to reaction. All energetically accessible reaction channels were found, including products containing carbon-oxygen bonds, which is in agreement with the results of recent experiments. After correcting for compensating errors in competing reaction channels, our excitation functions show quantitative agreement with experiment (for which absolute magnitudes of cross sections are available) at high collision energies (several eV). More detailed properties, such as translational and angular distributions, show qualitative agreement. The opacity function reveals a high selectivity for producing OH(+) at high impact parameters, CH(3)(+)/CH(2)(+)/H(2)O(+) at intermediate impact parameters, and H(2)CO(+)/HCO(+)/CO(+) at small impact parameters. Angular distributions for CH(3)(+)/CH(2)(+)/OH(+) are forward scattered at high collision energies which implies the importance of direct reaction mechanisms, while reaction complexes play an important role at lower energies, especially for the H(2)O(+) product. Finally, we find that the nominally spin-forbidden product CH(3)(+) + OH can be produced by a spin-allowed pathway that involves the formation of the triplet excited product CH(3)(+)(?(3)E). This explains why CH(3)(+) can have a high cross section, even at very low collision energies. The results of this work suggest that the PM3 method may be applied directly to the study of O(+) reactions with small alkane molecules and polymer surfaces.     

Photodissociation of methane: exploring potential energy surfaces

The potential energy surface for the first excited singlet state (S(1)) of methane is explored using multireference singles and doubles configuration interaction calculations, employing a valence triple zeta basis set. A larger valence quadruple zeta basis is used to calculate the vertical excitation energy and dissociation energies. All stationary points found on the S(1) surface are saddle points and have imaginary frequencies for symmetry-breaking vibrations. By studying several two-dimensional cuts through the potential energy surfaces, it is argued that CH(4) in the S(1) state will distort to planar structures. Several conical intersection seams between the ground state surface S(0) and the S(1) surface have been identified at planar geometries. The conical intersections provide electronically nonadiabatic pathways towards products CH(3)((approximately)X (2)A"(2))+H, CH(2)((approximately)a (1)A(1))+H(2), or CH(2)((approximately)X (3)B(1))+H+H. The present results thereby make it plausible that the CH(3)((approximately)X (2)A"(2))+H and CH(2)((approximately)a (1)A(1))+H(2) channels are major dissociation channels, as has been observed experimentally.     

CH3SS与OH自由基反应机理的理论研究

应用密度泛函理论(DFT)对CH3SS与OH自由基单重态反应机理进行了研究.在B3PW91/6-311+G(d,p)水平上优化了反应通道上各驻点(反应物、中间体、过渡态和产物)的几何构型,用内禀反应坐标(IRC)计算和频率分析方法对过渡态进行了验证.在QCISD(T)/6-311++G(d,p)水平上计算了各物种的单点能,并对总能量进行了零点能校正.研究结果表明,CH3SS与OH反应为多通道反应,有5条可能的反应通道.反应物首先通过不同的S—O键相互作用形成具有竞争反应机理的中间体IM1和IM2.再经过氢迁移、脱氢和裂解等机理得到主要产物P1(CH2SS+H2O),次要产物P2(CH2S+HSOH),P3(CH3SH+1SO)和P4(CH2SSO+H2),其中最低反应通道的势垒为174.6kJ.mol-1.     

Theoretical studies of the reactions of O(~3p) with halogenated methyl (Ⅰ)——Reaction mechanism of the O(~3p) CH_2Cl reaction

The reaction of O(~3P) with CH_2Cl radical has been studied using ab initio molecular orbital theory. G2 (MP2) method is used to calculate the geometrical parameters, vibrational frequencies and energies of various stationary points on the potential energy surface. The reaction mechanism is revealed. The addition of O(~3P) with CH_2Cl leads to the formation of an energy rich intermediate OCH_2Cl which can subsequently undergo decomposition or isomerization to the final products. The calculated heat of reaction for each channel is in agreement with the experimental value. The production of H CHClO and Cl CH_2O are predicted to be the major channels. The overall rate constants are calculated using transition state theory on the basis of ab initio data. The rate constant is pressure independent and exhibits negative temperature dependence at lower temperatures, in accordance with the experimental results.     

Theoretical study on the kinetics of OH radical reactions with CH3OOH and CH3CH2OOH

The mechanisms and kinetics studies of the OH radical with alkyl hydroperoxides CH(3)OOH and CH(3)CH(2)OOH reactions have been carried out theoretically. The geometries and frequencies of all the stationary points are calculated at the UBHandHLYP/6-311G(d,p) level, and the energy profiles are further refined by interpolated single-point energies method at the MC-QCISD level of theory. For two reactions, five H-abstraction channels are found and five products (CH(3)OO, CH(2)OOH, CH(3)CH(2)OO, CH(2)CH(2)OOH, and CH(3)CHOOH) are produced during the above processes. The rate constants for the CH(3)OOH/CH(3)CH(2)OOH + OH reactions are corrected by canonical variational transition state theory within 250-1500 K, and the small-curvature tunneling is included. The total rate constants are evaluated from the sum of the individual rate constants and the branching ratios are in good agreement with the experimental data. The Arrhenius expressions for the reactions are obtained.     

Experimental and modeling studies of the pressure and temperature dependences of the kinetics and the OH yields in the acetyl + O2 reaction

The acetyl + O(2) reaction has been studied by observing the time dependence of OH by laser-induced fluorescence (LIF) and by electronic structure/master equation analysis. The experimental OH time profiles were analyzed to obtain the kinetics of the acetyl + O(2) reaction and the relative OH yields over the temperature range of 213-500 K in helium at pressures in the range of 5-600 Torr. More limited measurements were made in N(2) and for CD(3)CO + O(2). The relative OH yields were converted into absolute yields by assuming that the OH yield at zero pressure is unity. Electronic structure calculations of the stationary points of the potential energy surface were used with a master equation analysis to fit the experimental data in He using the high-pressure limiting rate coefficient for the reaction, k(∞)(T), and the energy transfer parameter, (ΔE(d)), as variable parameters. The best-fit parameters obtained are k(∞) = 6.2 × 10(-12) cm(-3) molecule(-1) s(-1), independent of temperature over the experimental range, and (ΔE(d))(He) = 160(T/298?K) cm(-1). The fits in N(2), using the same k(∞)(T), gave (ΔE(d))(N(2)) = 270(T/298?K) cm(-1). The rate coefficients for formation of OH and CH(3)C(O)O(2) are provided in parametrized form, based on modified Troe expressions, from the best-fit master equation calculations, over the pressure and temperature ranges of 1 ≤ p/Torr ≤ 1.5 × 10(5) and 200 ≤ T/K ≤ 1000 for He and N(2) as the bath gas. The minor channels, leading to HO(2) + CH(2)CO and CH(2)C(O)OOH, generally have yields <1% over this range.     

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.     

On the Effect of Tether Composition on cis/trans Selectivity in Intramolecular Diels–Alder Reactions

Intramolecular Diels–Alder (IMDA) transition structures (TSs) and energies have been computed at the B3LYP/6‐31+G(d) and CBS‐QB3 levels of theory for a series of 1,3,8‐nonatrienes, H₂C?CH? CH?CH? CH₂? X? Z? CH?CH₂ [? X? Z? =? CH₂? CH₂? ( 1 ); ? O? C(?O)? ( 2 ); ? CH₂? C(?O)? ( 3 ); ? O? CH₂? ( 4 ); ? NH? C(?O)? ( 5 ); ? S? C(?O)? ( 6 ); ? O? C(?S)? ( 7 ); ? NH? C(?S)? ( 8 ); ? S? C(?S)? ( 9 )]. For each system studied ( 1 – 9 ), cis‐ and trans‐TS isomers, corresponding, respectively, to endo‐ and exo‐positioning of the ? C? X? Z? tether with respect to the diene, have been located and their relative energies (E_rel^TS) employed to predict the cis/trans IMDA product ratio. Although the E_rel^TS values are modest (typically <3 kJ mol^?1), they follow a clear and systematic trend. Specifically, as the electronegativity of the tether group X is reduced (X?O→NH or S), the IMDA cis stereoselectivity diminishes. The predicted stereochemical reaction preferences are explained in terms of two opposing effects operating in the cis‐TS, namely (1) unfavorable torsional (eclipsing) strain about the C4? C5 bond, that is caused by the ? C? X? C(?Y)? group’s strong tendency to maintain local planarity; and (2) attractive electrostatic and secondary orbital interactions between the endo‐(thio)carbonyl group, C?Y, and the diene. The former interaction predominates when X is weakly electronegative (X?N, S), while the latter is dominant when X is more strongly electronegative (X?O), or a methylene group (X?CH₂) which increases tether flexibility. These predictions hold up to experimental scrutiny, with synthetic IMDA reactions of 1 , 2 , 3 , and 4 (published work) and 5 , 6 , and 8 (this work) delivering ratios close to those calculated. The reactions of thiolacrylate 5 and thioamide 8 represent the first examples of IMDA reactions with tethers of these types. Our results point to strategies for designing tethers, which lead to improved cis/trans‐selectivities in IMDAs that are normally only weakly selective. Experimental verification of the validity of this claim comes in the form of fumaramide 14 , which undergoes a more trans‐selective IMDA reaction than the corresponding ester tethered precursor 13 .     

Pathways and rate coefficients for the decomposition of vinoxy and acetyl radicals

The potential in the vicinity of the stationary points on the surface for the decomposition of ground-state vinoxy and acetyl radicals has been calculated using the RQCISD(T) method extrapolated to the infinite-basis set limit. Rate coefficients for the decomposition pathways of these two radicals were computed using the master equation and variational transition state theory. Agreement between our calculated rate coefficients for H + CH(2)CO <--> CH(3) + CO and experimental data is very good, without the need for empirical adjustments to the ab initio energy barriers. Multireference configuration-interaction calculations indicate two competitive channels for vinoxy decomposition, with the channel leading to H + CH(2)CO being preferred at photodissociation energies. However, at typical combustion conditions, vinoxy decomposes primarily to CO and methyl. In contrast, decomposition of acetyl shows only one decomposition channel, leading to CO and methyl. The implications of a low-lying exit channel for the calculation of theoretical rate coefficients are discussed briefly.     

Approach to the atmospheric chemistry of methyl nitrate and methylperoxy nitrite. Chemical mechanisms of their formation and decomposition reactions in the gas phase

Potential energy surfaces, minimum energy reaction paths, minima, transition states, reaction barriers, and conical intersections for the most important atmospheric reactions of methyl nitrate (CH(3)ONO(2)) and methylperoxy nitrite (C(3)HOONO) on the electronic ground state have been studied (i) with the second-order multiconfigurational perturbation theory (CASPT2) by computation of numerical energy gradients for stationary points and (ii) with the density functional theory (DFT). The proposed mechanism explains the conversion of unreactive alkyl peroxy radicals into alkoxy radicals: CH(3)O(2) + NO <=> CH(3)OONO <=> CH(3)O + NO(2) left arrow over right arrow CH(3)ONO(2). Additionally, several discrepancies found in the comparison of the results obtained from the two employed approaches are analyzed. CASPT2 predicts that all dissociation reactions into radicals occur without an extra exit energy barrier. In contrast, DFT finds transition states for the dissociations of cis- and trans-methylperoxy nitrite into CH(3)O + NO(2). Furthermore, multiconfigurational methods [CASPT2 and complete active space SCF (CAS-SCF)] predict the isomerization of CH3ONO2 to CH3OONO to occur in a two-step mechanism: (i) CH(3)ONO(2) --> CH(3)O + NO(2); and (ii) CH(3)O + NO(2) --> CH(3)OONO. The reason for this has to do with the coupling of the ground electronic state with the first excited state. Therefore, it is demonstrated that DFT methods based on single determinantal wave functions give an incorrect picture of the aforementioned reaction mechanisms.     

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.