Kinetics and mechanism of the C6H5 + CH3CHO reaction: experimental measurement and theoretical prediction of the reactivity toward four molecular sites.

The kinetics and mechanism of the reaction of C6H5 with CH3CHO have been investigated experimentally and theoretically. The total rate constant for the reaction has been measured by means of the cavity ring-down spectrometry (CRDS) in the temperature range 299-501 K at pressures covering 20-75 Torr. The overall bimolecular rate constant can be represented by the expression k = (2.8 +/- 0.2) x 10(11) exp[-(700 +/- 30)/T] cm3 mol-1 s-1, which is slightly faster than for the analogous C6H5 + CH2O reaction determined with the same method in the same temperature range. The reaction mechanism for the C6H5 + CH3CHO reaction was also explored with quantum-chemical calculations at various hybrid density functional theories (DFTs) and using ab initio high-level composite methods. The theories predict that the reaction may occur by two hydrogen-abstraction and two addition channels with the aldehydic hydrogen-abstraction reaction being dominant. The rate constant calculated by the transition state theory for the aldehydic hydrogen-abstraction reaction is in good agreement with the experimental result after a very small adjustment of the predicted reaction barrier (+0.3 kcal mol-1). Contributions from other product channels are negligible under our experimental conditions. For combustion applications, we have calculated the rate constants for key product channels in the temperature range of 298-2500 K under atmospheric-pressure conditions; they can be represented by the following expressions in units of cm 3mol-1 s-1: k1,cho = 8.8 x 10(3)T2.6 exp(-90/T), k2,ch3 = 6.0 x 10(1)T3.3 exp(-950/T), k3a(C6H5COCH3 + H) = 4.2 x 10(5)T0.6 exp(-410/T) and k3b(C6H5CHO + CH3) = 6.6 x 10(9)T-0.5 exp(-310/T).     

Mechanism of cysteine oxidation by a hydroxyl radical: a theoretical study.

Cysteine oxidation by HO(.) was studied at a high level of ab initio theory in both gas phase and aqueous solution. Potential energy surface scans in the gas phase performed for the model system methanethiol+HO(.) indicate that the reactants can form two intermediate states: a sulfur-oxygen adduct and a hydrogen bound reactant complex. However these states appear to play a minor role in the reaction mechanism as long as they are fast dissociating states. Thus the main reaction channel predicted at the QCISD(T)/6-311+G(2df,2pd) level of theory is the direct hydrogen atom abstraction. The reaction mechanism is not perturbed by solvation which was found to induce only small variations in the Gibbs free energy of different reactant configurations. The larger size reactant system cysteine+HO* was treated by the integrated molecular orbital+molecular orbital (IMOMO) hybrid method mixing the QCISD(T)/6-311+G(2df,2pd) and the UMP2/6-311+G(d,p) levels of theory. The calculated potential energy, enthalpy, and Gibbs free energy barriers are slightly different from those of methanethiol. The method gave a rate constant for cysteine oxidation in aqueous solution, k=2.4 x 10(9) mol(-1) dm(3) s(-1), which is in good agreement with the experimental rate constant. Further analysis showed that the reaction is not very sensitive to hydrogen bonding and electrical polarity of the molecular environment.     

Ab initio studies of CIO, reactions: prediction of the rate constants of CIO + NO for the forward and reverse processes.

The mechanisms for ClO+NO and its reverse reactions were investigated by means of ab initio molecular orbital and statistical theory calculations. The species involved were optimized at the B3LYP/6-311 +G(3df) level, and their energies were refined at the CCSD(T)/6-311+ G(3df)//B3LYP/6-311 + G(3df) level. Five isomers and the transition states among them were located. The relative stability of these isomers is ClNO2 > cis-ClONO > trans-ClONO > cis-OClNO>trans-OClNO. The heats of formation of the three most-stable isomers were predicted using isodesmic reactions by different methods. The predicted bimolecular reaction rate constant shows that, below 100 atm, the formation of Cl+NO2 is dominant and pressure-independent. The total rate constant can be expressed as: k(ClO+NO)= 1.43 x 10(-9)T(-083)exp(92/ T) cm3 molecule(-1)s(-1) in the temperature range of 200-1000 K, in close agreement with experimental data. For the reverse reaction, Cl+NO2-->ClNO2 and ClONO (cis and trans isomers), the sum of the predicted rate constants for the formation of the three isomers and their relative yields also reproduce the experimental data well. The predicted total third-order rate constants in the temperature range of 200-1000 K can be represented by: k0(He) = 4.89 x 10(-6)T(-5.85) exp(-796/T) cm6 molecule(-1)s(-1) and k0(N2) =5.72 x 10(-15)T(-5.80) exp(-814/T) cm6 molecule(-1)s(-1). The predicted high- and low-pressure limit decomposition rates of CINO2 in Ar in the temperature range 400-1500 K can be expressed, respectively, by: k-(ClNO2) = 7.25 x 10(19)T(-1.89) exp(-16875/T) s(-1) and kd(ClNO2) = 2.51 x 10(38)T(-6.8) exp(-18409/T) cm3 molecule(-1) s(-1). The value of k0(ClNO2) is also in reasonable agreement with available experimental data.     

A theoretical investigation of the gas-phase oxidation reaction of the saturated tert-butyl radical.

The radical-radical reaction mechanisms and dynamics of ground-state atomic oxygen [O(3P)] with the saturated tert-butyl radical (t-C4H9) are investigated using the density functional method and the complete basis set model. Two distinctive reaction pathways are predicted to be in competition: addition and abstraction. The barrierless addition of O(3P) to t-C4H9 leads to the formation of an energy-rich intermediate (OC4H9) on the lowest doublet potential energy surface, which undergoes subsequent direct elimination or isomerization-elimination leading to various products: C3H6O + CH3, iso-C4H8O + H, C3H7O + CH2, and iso-C4H8 + OH. The respective microscopic reaction processes examined with the aid of statistical calculations, predict that the major addition pathway is the formation of acetone (C3H6O) + CH3 through a low-barrier, single-step cleavage. For the direct, barrierless H-atom abstraction mechanism producing iso-C4H8 (isobutene) + OH, which was recently reported in gas-phase crossed-beam investigations, the reaction is described in terms of both an abstraction process (major) and a short-lived addition dynamic complex (minor).     

Mechanism and rate constants of the CH2 + CH2CO reactions in triplet and singlet states: A theoretical study

Ab initio and density functional CCSD(T)-F12/cc-pVQZ-f12//B2PLYPD3/6-311G** calculations have been performed to unravel the reaction mechanism of triplet and singlet methylene CH₂ with ketene CH₂CO. The computed potential energy diagrams and molecular properties have been then utilized in Rice–Ramsperger–Kassel–Marcus-Master Equation (RRKM-ME) calculations of the reaction rate constants and product branching ratios combined with the use of nonadiabatic transition state theory for spin-forbidden triplet-singlet isomerization. The results indicate that the most important channels of the reaction of ketene with triplet methylene lead to the formation of the HCCO + CH₃ and C₂H₄ + CO products, where the former channel is preferable at higher temperatures from 1000 K and above. In the C₂H₄ + CO product pair, the ethylene molecule can be formed either adiabatically in the triplet electronic state or via triplet-singlet intersystem crossing in the singlet electronic state occurring in the vicinity of the CH₂COCH₂ intermediate or along the pathway of CO elimination from the initial CH₂CH₂CO complex. The predominant products of the reaction of ketene with singlet methylene have been shown to be C₂H₄ + CO. The formation of these products mostly proceeds via a well-skipping mechanism but at high pressures may to some extent involve collisional stabilization of the CH₃CHCO and cyclic CH₂COCH₂ intermediates followed by their thermal unimolecular decomposition. The calculated rate constants at different pressures from 0.01 to 100 atm have been fitted by the modified Arrhenius expressions in the temperature range of 300–3000 K, which are proposed for kinetic modeling of ketene reactions in combustion. © 2018 Wiley Periodicals, Inc.     

Ab initio studies of ClO(x) reactions: prediction of the rate constants of ClO+NO2 for the forward and reverse processes.

The potential-energy surface for the reaction of ClO with NO2 has been constructed at the CCSD(T)/6-311+G(3df)//B3LYP/6-311+G(3df) level of theory. Six ClNO3 isomers are located; these are ClONO2, pc-ClOONO, pt-ClOONO, OClNO2, pt-OClONO, pc-OClONO, with predicted energies relative to the reactants of -25.6, -0.5, 1.0, 1.9, 12.2 and 13.6 kcal mol-1, and heats of formation at 0 K of 7.8, 32.9, 34.4, 35.5, 45.6 and 47.0 kcal mol-1, respectively. Isomerizations among them are also discussed. The rate constants for the low-energy pathways have been computed by statistical theory calculations. For the association reaction producing exclusively ClONO2, the predicted low- and high-pressure-limit rate constants in N2 for the temperature range of 200-600 K can be represented by: (N2)=3.19 x 10-17 T-5.54 exp(-384 K/T) cm6 molecule-2 s-1 and =3.33 x 10-7 T-1.48 exp(-18 K/T) cm3 molecule-1 s-1. The predicted low- and high-pressure-limit rate constants for the decomposition of ClONO2 in N2 at 200-600 K can be expressed, respectively, by =6.08 x 1013 T-6.54 exp(-13813 K/T) cm3 molecule-1 s-1 and =4.59 x 1023 T-2.43 exp(-13437 K/T) s-1. The predicted values compare satisfactorily with available experimental data. The reverse Cl+NO3 reaction was found to be independent of the pressure, giving exclusively ClO+NO2; the predicted rate constant can be expressed as k(Cl+NO3)=1.19 x 10-9 T-0.60 exp(58 K/T) cm3 molecule-1 s-1..     

Theoretical dynamic studies on the reactions of CH3C(O)CH3-nCl(n) (n = 0-3) with the chlorine atom

The theoretical investigations were performed on the reaction mechanisms for the title reactions CH(3)C(O)CH(3) + Cl --> products (R1), CH(3)C(O)CH(2)Cl + Cl --> products (R2), CH(3)C(O)CHCl(2) + Cl --> products (R3), and CH(3)C(O)CCl(3) + Cl --> products (R4) by ab initio direct dynamics approach. Two different reaction channels have been found: abstract of the H atom from methyl (--CH(3)) group or chloromethyl (--CH(3-n)Cl(n)) group of chloroacetone and addition of a Cl atom to the carbon atom of the carbonyl group of chloroacetone followed by methyl or chloromethyl eliminations. Because of the higher potential energy barrier, the contribution of addition-elimination reaction pathway to the total rate constants is very small and thus this pathway is insignificant in atmospheric conditions. The rate constants for the H-abstraction reaction channels are evaluated by using canonical variational transition state theory incorporating with the small-curvature tunneling correction. Theoretical overall rate constants are in good agreement with the available experimental values and decrease in the order of k(1) > k(2) > k(3) > k(4). The results indicate that for halogenated acetones the substitution of halogen atom (F or Cl) leads to the decrease in the C--H bond reactivity and more decrease of reactivity is caused by F-substitution.     

Experimental and Theoretical Study of Helium Broadening and Shift of HCO+ Rotational Lines

An experimental and theoretical study of pressure broadening and pressure shift of HCO⁺ rotational lines perturbed by collisions with He is presented. Results are reported from measurements at 88 K for the lines j=4←3, 5←4 and 6←5 with frequencies ranging from 0.35 to 0.54 THz. Using a new CCSD(T)/aug‐cc‐pVQZ potential energy surface for the He–HCO⁺ interaction, the collisional line shape parameters are studied from fully quantum and semiclassical calculations. Results from the quantum treatment are in satisfactory agreement with experiments whereas the semiclassical approach can lead to appreciable differences. A study of the dependence of line width Γ and shift s as a function of the translational energy shows the presence of quantum oscillations. Calculations on a previous Hartree–Fock‐based potential energy surface lead to quite similar results for the collisional line shape parameters. Using a simplified version of the potential morphing method it is found that the line width Γ is particularly sensitive to the long‐range part of the potential energy surface. This also explains the success of the first line‐broadening calculations which date back to the 1950s.     

A theoretical ab initio study on the H(2)NO + O(3) reaction

The deviation of the NH(2) pseudo-first-order decay Arrhenius plots of the NH(2) + O(3) reaction at high ozone pressures measured by experimentalists, has been attributed to the regeneration of NH(2) radicals due to the subsequent reactions of the products of this reaction with ozone. Although these products have not yet been characterized experimentally, the radical H(2)NO has been postulated, because it can regenerate NH(2) radicals through the reactions: H(2)NO + O(3) --> NH(2) + O(2) and H(2)NO + O(3) --> HNO + OH + O(2). With the purpose of providing a reasonable explanation from a theoretical point of view to the kinetic observed behaviour of the NH(2) + O(3) system, we have carried ab initio electronic structure calculations on both H(2)NO + O(3) possible reactions. The results obtained in this article, however, predict that of both reactions proposed, only the H(2)NO + O(3) --> NH(2) + O(2) reaction would regenerate indeed NH(2) radicals, explaining thus the deviation of the NH(2) pseudo-first-order decay observed experimentally.     

A theoretical study of the reaction of HCO+ with C2H2

A theoretical study was performed for the reaction of formyl cation and acetylene to give C₃H+O in flames and C₂H (nonclassical)+CO, both in flames and in interstellar clouds. The corresponding Potential Energy Surface (PES) was studied at the B3LYP/cc‐pVTZ level of theory, and single‐point calculations on the B3LYP geometries were carried out at the CCSD(T)/cc‐pVTZ level. Our results display a route to propynal evolving energetically under C₂H (nonclassical)+CO and, consequently, accessible in interstellar clouds conditions. This route connects the most stable C₃H₃O⁺ isomer (C2‐protonated propadienone) with a species from which propynal may be produced in a dissociative electron recombination reaction. The reaction channel to produce the C₃H+O evolves basically through two TSs and presents an endothermicity of 63.9 kcal/mol at 2000 K. According to our Gibbs energy profiles, the C2‐protonated propadienone is the most stable species at low–moderate temperatures and, consequently, could play a certain role in interstellar chemistry. On the contrary, in combustion chemistry conditions (2000 K) the C₂H (nonclassical)+CO products are the most thermodynamically favored species. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 35–42, 2000     

Ab initio study of the mechanism of the atmospheric reaction: NO2 + O3 --> NO3 + O2

The atmospheric reaction NO2 + O3 --> NO3 + O2 (1) has been investigated theoretically by using the MP2, G2, G2Q, QCISD, QCISD(T), CCSD(T), CASSCF, and CASPT2 methods with various basis sets. The results show that the reaction pathway can be divided in two different parts at the MP2 level of theory. At this level, the mechanism proceeds along two transition states (TS1 and TS2) separated by an intermediate, designated as A. However, when the single-reference higher correlated QCISD methodology has been employed, the minimum A and the transition state TS2 are not found on the hypersurface of potential energy, which confirms a direct reaction mechanism. Single-reference high correlated and multiconfigurational methods consistently predict the barrier height of reaction (1) to be within the range 2.5-6.1 kcal mol(-1), in reasonable agreement with experimental data. The calculated reaction enthalpy is -24.6 kcal mol(-1) and the reaction rate calculated at the highest CASPT2 level, of k = 6.9 x 10(-18) cm(3) molecule(-1) s(-1). Both results can be regarded also as accurate predictions of the methodology employed in this article.     

Multichannel RRKM study on the mechanism and kinetics of the CH3Cl + OH reaction

The kinetics and mechanism of the reaction of OH with CH₃Cl have been theoretically studied. The potential energy surface for each possible pathway has been investigated by the G2MP2 method. The rate constants for channels leading to several products have been calculated by multichannel‐Rice‐Ramsperger‐Kassel‐Marcus (RRKM) theory over a temperature range 200–2000 K. The results show the major channel is hydrogen abstraction mechanism. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

A theoretical characterization of the photoisomerization channels of 1,2-cyclononadienes on both singlet and triplet potential-energy surfaces

The ground-, (1)(pipi*)-, and (3)(pipi*)-state potential-energy surfaces of 1,2-cyclononadiene and isomeric C(9)H(14) species, as well as 1-methyl-1,2-cyclononadiene and isomeric C(10)H(16) species were all mapped using CASSCF and the 6-31G(d) basis set. Theoretical results were found to be in good agreement with the available experimental observations for both 1,2-cyclononadiene and 1-methyl-1,2-cyclononadiene isomerization reactions under singlet and triplet direct or sensitized irradiation. Extremely efficient decay occurs from the first singlet excited state to the ground state through at least three different conical intersections (surface crossings). The first of these crossing points is accessed by a one-bond ring closure. From this conical intersection point (CI-A or CI-C), some possible subsequent ground-state reaction paths have been identified: 1) intramolecular C--H bond insertion to form the bicyclic photoproduct and 2) intramolecular C--H bond insertion to form tricyclic photoproducts. An excited state [1,3]-sigmatropic shift leads to the second conical intersection (CI-B or CI-E), which can give a three-bond cyclononyne species. Besides these, in the singlet photochemical reactions of 1-methyl-1,2-cyclononadiene, excited-state, one allenic C--H bond insertion leads to a third conical intersection (CI-D). Possible ground-state reaction pathways from this structure lead to the formation of a diene photoproduct or to transannular insertion photoproducts. Moreover, in the case of triplet 1,2-cyclononadiene and 1-methyl-1,2-cyclononadiene photoisomerization reactions, both chemical reactions will adopt a 1,3-biradical (T(1)/S(0)-1, T(1)/S(0)-2, and T(1)/S(0)-3), which may undergo intersystem crossings leading to the formation of tricyclic or bicyclic photoproducts. The results obtained allow a number of predictions to be made.     

Ab Initio Study on the Mechanism of the Atmospheric Reaction OH+O3→HO2+O2

The atmospheric reaction (1) OH + O₃→HO₂ + O₂ was investigated theoretically by using MP2, QCISD, QCISD(T), and CCSD(T) methods with various basis sets. At the highest level of theory, namely, QCISD, the reaction is direct, with only one transition state between reactants and products. However, at the MP2 level, the reaction proceeds through a two‐step mechanism and shows two transition states, TS1 and TS2 , separated by an intermediate, Int . The different methodologies employed in this paper consistently predict the barrier height of reaction (1) to be within the range 2.16–5.11 kcal mol^?1, somewhat higher than the experimental value of 2.0 kcal mol^?1.     

A computational study of thiolate and selenolate oxidation by hydrogen peroxide.

Ab initio molecular orbital calculations have been used to study the effects of the molecular environment on the oxidation of thiolate and selenolate by hydrogen peroxide. The reaction was first examined in vacuo at the QCISD(T)/6-311+G(2df,2pd)//MP2/6-311+G(d,p) level of theory. It was found for both thiolate and selenolate that a reactant aggregate is formed, which has a dissociation rate constant comparable to the activation rate constant (about 10(-3) s(-1) for thiolate and 10(-1) s(-1) for selenolate). Using the polarizable continuum model (PCM) it was then found that the dissociation barrier energy decreases dramatically in water giving a dissociation rate constant of the order of 10(9) s(-1). In this case, the predicted overall rate constant of the thiolate reaction was about 10.2 mol(-1) dm3 s(-1), which is in good agreement with the experimental rate constant of cysteine oxidation in aqueous solution. The calculated rate constant for the selenolate reaction was somewhat higher (about 35.4 mol(-1) dm3 s(-1)). However, this value is several orders of magnitude smaller than the experimental value reported for the oxidation of selenocysteine in glutathione peroxidase. By considering the effect of the PCM dielectric constant on the reaction rate constant it was concluded that the high reactivity of the selenocysteine in glutathione peroxidase, as compared with cysteine, could be mainly due to the molecular environment of the selenocysteine residue.     

Definitive ab initio studies of model SN2 reactions CH(3)X+F- (X=F,Cl, CN,OH, SH,NH(2), PH(2))

The energetics of the stationary points of the gas-phase reactions CH(3)X+F(-)-->CH(3)F+X(-) (X=F, Cl, CN, OH, SH, NH(2) and PH(2)) have been definitively computed using focal point analyses. These analyses entailed extrapolation to the one-particle limit for the Hartree-Fock and MP2 energies using basis sets of up to aug-cc-pV5Z quality, inclusion of higher-order electron correlation [CCSD and CCSD(T)] with basis sets of aug-cc-pVTZ quality, and addition of auxiliary terms for core correlation and scalar relativistic effects. The final net activation barriers for the forward reactions are: E (b/F,F)=-0.8, E (b/F, Cl)=-12.2, E (b/F,OH)=+13.6, E b/F,OH=+16.1, E b/F,SH=+2.8, Eb/F, NH=+32.8, and E b/F,PH =+19.7 kcal x mol(-1). For the reverse reactions E b/F,F= -0.8, Eb/Cl,F =+18.3, E b/CN,F=+12.2, E b/OH,F =-1.8, E b/SH,F =+13.2, E b/NH(2),=-1.5, and E b/PH(2) =+9.6 kcal x mol(-1). The change in energetics between the CCSD(T)/aug-cc-pVTZ reference prediction and the final extrapolated focal point value is generally 0.5-1.0 kcal mol(-1). The inclusion of a tight d function in the basis sets for second-row atoms, that is, utilizing the aug-cc-pV(X+d)Z series, appears to change the relative energies by only 0.2 kcal x mol(-1). Additionally, several decomposition schemes have been utilized to partition the ion-molecule complexation energies, namely the Morokuma-Kitaura (MK), reduced variational space (RVS), and symmetry adapted perturbation theory (SAPT) techniques. The reactant complexes fall into two groups, mostly electrostatic complexes (FCH(3).F(-) and ClCH(3).F(-)), and those with substantial covalent character (NCCH(3).F(-), CH(3)OH.F(-), CH(3)SH.F(-), CH(3)NH(2).F(-) and CH(3)PH(2).F(-)). All of the product complexes are of the form FCH(3).X(-) and are primarily electrostatic.     

Gas‐phase reactions of SiHn+ (n = 1,2) with NF3: A computational investigation on the detailed mechanistic aspects

The mechanism of the gas‐phase reactions of SiH_n⁺ (n = 1,2) with NF₃ were investigated by ab initio calculations at the MP2 and CAS‐MCSCF level of theory. In the reaction of SiH⁺, the kinetically relevant intermediates are the two isomeric forms of fluorine‐coordinated intermediate HSi‐F‐NF₂⁺. These species arise from the exoergic attack of SiH⁺ to one of the F atoms of NF₃ and undergo two competitive processes, namely an isomerization and subsequent dissociation into SiF⁺ + HNF₂, and a singlet‐triplet crossing so to form the spin‐forbidden products HSiF⁺ + NF₂. The reaction of SiH₂⁺ with NF₃ involves instead the concomitant formation of the nitrogen‐coordinated complex H₂Si‐NF₃⁺ and of the fluorine‐coordinated complex H₂Si‐F‐NF₂⁺. The latter isomer directly dissociates into NF₂⁺ + H₂SiF, whereas the former species preferably undergoes the passage through a conical intersection point so to form a H₂SiF‐NF₂⁺ isomer, which eventually dissociates into H₂SiF⁺ and NF₂. © 2012 Wiley Periodicals, Inc.     

An Ab Initio Study on the Mechanism of the Atmospheric Reaction NH2+O3→H2NO+O2

The atmospheric reaction NH₂+O₃→H₂NO+O₂ has been investigated theoretically by using MP2, QCISD, QCISD(T), CCSD(T), CASSCF, and CASPT2 methods with various basis sets. At the MP2 level of theory, the hypersurface of the potential energy (HPES) shows a two step reaction mechanism. Therefore, the mechanism proceeds along two transition states (TS1 and TS2), separated by an intermediate designated as Int. However, when the single‐reference higher correlated QCISD and the multiconfigurational CASSCF methodologies have been employed, the minimum structure Int and TS2 are not found on the HPES, which thus confirms a direct reaction mechanism. Single‐reference high correlated and multiconfigurational methods consistently predict the barrier height of the reaction to be within the range of 3.9 to 6.6 kcal mol^?1, which is somewhat higher than the experimental value. ¹ The calculated reaction enthalpy is ?67.7 kcal mol^?1.