Association reaction between SiH3 and H2O2: a computational study of the reaction mechanism and kinetics

The association reaction between silyl radical (SiH₃) and H₂O₂ has been studied in detail using high-level composite ab initio CBS-QB3 and G4MP2 methods. The global hybrid meta-GGA M06 and M06-2X density functionals in conjunction with 6-311++G(d,p) basis set have also been applied. To understand the kinetics, variational transition-state theory calculation is performed on the first association step, and successive unimolecular reactions are subjected to Rice–Ramsperger–Kassel–Marcus calculations to predict the reaction rate constants and product branching ratios. The bimolecular rate constant for SiH₃–H₂O₂ association in the temperature range 250–600 K, k(T) = 6.89 × 10^?13 T ^?0.163exp(?0.22/RT) cm³ molecule^?1 s^?1 agrees well with the current literature. The OH production channel, which was experimentally found to be a minor one, is confirmed by the rate constants and branching ratios. Also, the correlation between our theoretical work and experimental literature is established. The production of SiO via secondary reactions is calculated to be one of the major reaction channels from highly stabilized adducts. The H-loss pathway, i.e., SiH₂(OH)₂ + H, is the major decomposition channel followed by secondary dissociation leading to SiO.     

Quantum chemical and theoretical kinetics study of the O(3P) + C2H2 reaction: a multistate process

The potential energy surfaces of the two lowest-lying triplet electronic surfaces 3A' and 3A' for the O(3P) + C2H2 reaction were theoretically reinvestigated, using various quantum chemical methods including CCSD(T), QCISD, CBS-QCI/APNO, CBS-QB3, G2M(CC,MP2), DFT-B3LYP and CASSCF. An efficient reaction pathway on the electronically excited 3A' surface resulting in H(2S) + HCCO(A2A') was newly identified and is predicted to play an important role at higher temperatures. The primary product distribution for the multistate multiwell reaction was then determined by RRKM statistical rate theory and weak-collision master equation analysis using the exact stochastic simulation method. Allowing for nonstatistical behavior of the internal rotation mode of the initial 3A' adducts, our computed primary-product distributions agree well with the available experimental results, i.e., ca. 80% H(2S) + HCCO(X2A' + A2A') and 20% CH2(X3B1) + CO(X1sigma+) independent of temperature and pressure over the wide 300-2000 K and 0-10 atm ranges. The thermal rate coefficient k(O + C2H2) at 200-2000 K was computed using multistate transition state theory: k(T) = 6.14 x 10(-15)T (1.28) exp(-1244 K/T) cm3 molecule(-1) s(-1); this expression, obtained after reducing the CBS-QCI/APNO ab initio entrance barriers by 0.5 kcal/mol, quasi-perfectly matches the experimental k(T) data over the entire 200-2000 K range, spanning 3 orders of magnitude.     

Calculations on the competition between association and reaction for C3H(+) + H2

A potential energy surface has been calculated for the competing associative and reactive ion-molecule processes involving the reactants C3H(+) + H2. Our ab initio results show that the linear ion C3H+ and H2 can directly access the deep potential well of the propargyl ion H2CCCH+, which is calculated to lie 390 kJ mol-1 below the zero-point energy of the reactants. Isomerization between the propargyl ion and the lower energy, cyclic C3H3+ ion, calculated to lie 501 kJ mol-1 below the zero-point energy of reactants, can subsequently occur via two pathways. One of these pathways involves a transition state lying 22 kJ mol-1 below the energy of the reactants while the other, which occurs at much lower energies, involves two transition states and an intermediate. The dissociation of c-C3H3+ into c-C3H2(+) + H is calculated to occur directly, without any intermediate potential energy maximum, but the energy of the products lies 7.3 kJ mol-1 above the energy of the reactants. Using the minimum energy potential pathway and properties of the stationary point structures determined via ab initio methods, we have calculated both the association rate coefficient to produce C3H3+ as a function of density and the branching ratio between the propargyl and cyclic structures of the ion. Our results are in good agreement with some experimental results and in conflict with others. Specifically, we agree with the 1:1 branching ratio measured for the propargyl and cyclic isomers of C3H3+ at 80 and 300 K and we agree with the rate coefficient for radiative association measured at 80 K. We cannot reproduce reported measurements that the reactive channel (C3H2(+) + H) is the dominant channel at 80 K and at low gas densities, or that the association channel at high densities saturates at an effective rate coefficient well below the Langevin value -2x10(-11) cm3 s-1 at 300K and 1x10(-10) cm3 s-1 at 80K.     

Radical-molecule reaction C3H+H2O: a mechanistic study

Despite the importance of the C(3)H radical in both combustion and interstellar space, the reactions of C(3)H toward stable molecules have never been studied. In this paper, we report our detailed mechanistic study on the radical-molecule reaction C(3)H+H(2)O at the Becke's three parameter Lee-Yang-Parr-B3LYP6-311G(d,p) and coupled cluster with single, double, and triple excitations-CCSD(T)6-311G(2d,p) (single-point) levels. It is shown that the C(3)H+H(2)O reaction initially favors formation of the carbene-insertion intermediates HCCCHOH (1a,1b) rather than the direct H- or OH-abstraction process. Subsequently, the isomers (1a,1b) can undergo a direct H- extrusion to form the well-known product propynal HCCCHO (P(5)). Highly competitively, (1a,1b) can take the successive 1,4- and 1,2-H-shift interconversion to isomer H(2)CCCHO(2a,2b) and then to isomer H(2)CCHCO(3a,3b), which can finally take a direct C-C bond cleavage to give product C(2)H(3) and CO (P(1)). The other products are kinetically much less feasible. With the overall entrance barrier 10.6 kcal/mol, the title reaction can be important in postburning processes. Particularly, our calculations suggest that the title reaction may play a role in the formation of the intriguing interstellar molecule, propynal HCCCHO. The calculated results will also be useful for the analogous C(3)H reactions such as with ammonia and alkanes.     

Analogies between point mechanics and chemical reaction kinetics

A few rather interesting analogies are shown between classical point mechanics and chemical reaction kinetics. The approach used is that of the differential geometry.

---

. — .

     

Crossed beams study of the reaction 1CH2+C2H2-->C3H3+H

The reaction of electronically excited singlet methylene (1CH2) with acetylene (C2H2) was studied using the method of crossed molecular beams at a mean collision energy of 3.0 kcal/mol. The angular and velocity distributions of the propargyl radical (C3H3) products were measured using single photon ionization (9.6 eV) at the advanced light source. The measured distributions indicate that the mechanism involves formation of a long-lived C3H4 complex followed by simple C-H bond fission producing C3H3+H. This work, which is the first crossed beams study of a reaction involving an electronically excited polyatomic molecule, demonstrates the feasibility of crossed molecular beam studies of reactions involving 1CH2.     

NH2(X^2B1)与C2H4反应机理的从头算研究

用从头算MP2方法,在6-311G^**基组下,对NH2X^2B1)与C2H4的加成和氢迁移反应机理进行了研究,优化得到反应的过渡态,并通过振动分析和内禀反应坐标(IRC)加以证实,计算了两个反应的能垒和1500K～2000K温度范围内的速率常数。结果表明：在1500K～2000K温度范围内加成反应是NH2(XX^2B1)与C2H4的反应的主要通道,同时报道了两个反应沿反应路径变化信息。     

Thermal decomposition of ethanol. III. A computational study of the kinetics and mechanism for the CH3+C2H5OH reaction

The mechanism for the CH3+C2H5OH reaction has been investigated by the modified Gaussian-2 method based on the geometric parameters of the stationary points optimized at the B3LYP/6-311+G(d,p) level of theory. Five transition states have been identified for the production of CH4+CH3CHOH (TS1), CH4+CH3CH2O (TS2), CH4+CH2CH2OH (TS3), CH3OH+CH3CH2 (TS4), and CH3CH2OCH3+H (TS5) with the corresponding barriers 12.0, 13.2, 16.0, 44.7, and 49.9 kcal/mol, respectively. The predicted rate constants and branching ratios for the three lower-energy H-abstraction reactions were calculated using the conventional and variational transition state theory with quantum-mechanical tunneling corrections for the temperature range 300-3000 K. The predicted total rate constant, kt=8.36 x 10(-76) T(20.00) exp(5258/T) cm3 mol(-1) s(-1) (300-600 K) and 6.10 x 10(-25) T(4.10)exp(-4058/T) cm3 mol(-1) s(-1) (600-3000 K), agrees closely with existing experimental data in the temperature range 403-523 K. Similarly, the predicted rate constants for CH3+CH3CD2OH and CD3+C2H5OD are also in reasonable agreement with available low temperature kinetic data.     

Mechanism and kinetics of the reaction NO3 + C2H4

The reaction of NO(3) radical with C(2)H(4) was characterized using the B3LYP, MP2, B97-1, CCSD(T), and CBS-QB3 methods in combination with various basis sets, followed by statistical kinetic analyses and direct dynamics trajectory calculations to predict product distributions and thermal rate constants. The results show that the first step of the reaction is electrophilic addition of an O atom from NO(3) to an olefinic C atom from C(2)H(4) to form an open-chain adduct. A concerted addition reaction mechanism forming a five-membered ring intermediate was investigated, but is not supported by the highly accurate CCSD(T) level of theory. Master-equation calculations for tropospheric conditions predict that the collisionally stabilized NO(3)-C(2)H(4) free-radical adduct constitutes 80-90% of the reaction yield and the remaining products consist mostly of NO(2) and oxirane; the other products are produced in very minor yields. By empirically reducing the barrier height for the initial addition step by 1 kcal mol(-1) from that predicted at the CBS-QB3 level of theory and treating the torsional modes explicitly as one-dimensional hindered internal rotations (instead of harmonic oscillators), the computed thermal rate constants (including quantum tunneling) can be brought into very good agreement with the experimental data for the overall reaction rate constant.     

Group theory and reaction mechanisms: Permutation theoretic prediction and computational support for pseudorotation modes in C2H and C5H rearrangements

The McIver–Stanton rules concerning the symmetry of transition states have a counterpart in rules concerning the permutation symmetry of single steps in degenerate rearrangements, derivable with the aid of Longuet–Higgins group theory. The generalized rules are illustrated by the widely studied PX₅ polytopal rearrangements. The analysis leads to prediction of hitherto unexplored “pseudorotation” pathways for rearrangements in ethyl and homotetrahedryl cations. CNDO computations of system energies, gradients, and curvatures at critical points on the C₂H and C₅H surfaces indicate that symmetry-breaking in keeping with the permutation-theoretic predictions is a key feature of the low-energy rearrangements of these systems. In particular, computation indicates that the C_2v “classical” homotetrahedral cation corresponds to an energy maximum rather than an energy minimum, or a transition state.     

Antioxidant activity of trans-resveratrol toward hydroxyl and hydroperoxyl radicals: a quantum chemical and computational kinetics study

In this work, we have carried out a systematic study of the antioxidant activity of trans-resveratrol toward hydroxyl ((?)OH) and hydroperoxyl ((?)OOH) radicals in aqueous simulated media using density functional quantum chemistry and computational kinetics methods. All possible mechanisms have been considered: hydrogen atom transfer (HAT), proton-coupled electron transfer (PCET), sequential electron proton transfer (SEPT), and radical adduct formation (RAF). Rate constants have been calculated using conventional transition state theory in conjunction with the Collins-Kimball theory. Branching ratios for the different paths contributing to the overall reaction, at 298 K, are reported. For the global reactivity of trans-resveratrol toward (?)OH radicals, in water at physiological pH, the main mechanism of reaction is proposed to be the sequential electron proton transfer (SEPT). However, we show that trans-resveratrol always reacts with (?)OH radicals at a rate that is diffusion-controlled, independent of the reaction pathway. This explains why trans-resveratrol is an excellent but very unselective (?)OH radical scavenger that provides antioxidant protection to the cell. Reaction between trans-resveratrol and the hydroperoxyl radical occurs only by phenolic hydrogen abstraction. The total rate coefficient is predicted to be 1.42 × 10(5) M(-1) s(-1), which is much smaller than the ones for reactions of trans-resveratrol with (?)OH radicals, but still important. Since the (?)OOH half-life time is several orders larger than the one of the (?)OH radical, it should contribute significantly to trans-resveratrol oxidation in aqueous biological media. Thus, trans-resveratrol may act as an efficient (?)OOH, and also presumably (?)OOR, radical scavenger.     

A systematic computational study of the reactions of HO2 with RO2: the HO2 + C2H5O2 reaction

The reaction of HO2 with C2H5O2 has been studied using the density functional theory (B3LYP) and the coupled-cluster theory [CCSD(T)]. The reaction proceeds on the triplet potential energy surface via hydrogen abstraction to form ethyl hydroperoxide and oxygen. On the singlet potential energy surface, the addition-elimination mechanism is revealed. Variational transition state theory is used to calculate the temperature-dependent rate constants in the range 200-1000 K. At low temperatures (e.g., below 300 K), the reaction takes place predominantly on the triplet surface. The calculated low-temperature rate constants are in good agreement with the experimental data. As the temperature increases, the singlet reaction mechanism plays more and more important role, with the formation of OH radical predominantly. The isotope effect of the reaction (DO2 + C2D5O2 vs HO2 + C2H5O2) is negligible. In addition, the triplet abstraction energetic routes for the reactions of HO2 with 11 alkylperoxy radicals (CnHmO2) are studied. It is shown that the room-temperature rate constants have good linear correlation with the activation energies for the hydrogen abstraction.     

Higher hydrogen uptake capacity of C2H4Ti+ than C2H4Ti: a quantum chemical study

Using density functionals theory, we show that gravimetric hydrogen uptake of C₂H₄Ti complex and its cation, C₂H₄Ti⁺, differ by about 2 wt%. Six and five hydrogen molecules are found to be adsorbed on C₂H₄Ti⁺ and C₂H₄Ti complexes thereby showing a hydrogen-uptake capacity of 13.74 and 11.72 wt%, respectively. All hydrogen molecules are adsorbed in molecular form on C₂H₄Ti⁺ ion with an increase in metal bond strength, whereas in some cases, the hydrogen molecules are found to be dissociated on C₂H₄Ti neutral complex. The uptake capacity of neutral C₂H₄Ti complex shown in this work is in excellent agreement with that reported experimentally, Phillips and Shivaram (Phys Rev Lett 100:105505, 2008). The H₂ adsorption energy and its dependence on exchange and correlation functions in density functionals method were illustrated. Even after the adsorption of maximum number of hydrogen molecules on C₂H₄Ti and C₂H₄Ti⁺ complexes, Ti and Ti⁺ remain strongly bound to C₂H₄ substrate.     

Addition reaction of adamantylideneadamantane with Br2 and 2Br2: a computational study

Ab initio calculations were carried out for the reaction of adamantylideneadamantane (Ad=Ad) with Br2 and 2Br2. Geometries of the reactants, transition states, intermediates, and products were optimized at HF and B3LYP levels of theory using the 6-31G(d) basis set. Energies were also obtained using single point calculations at the MP2/6-31G(d)//HF/6-31G(d), MP2/6-31G(d)//B3LYP/6-31G(d), and B3LYP/6-31+G(d)//B3LYP/6-31G(d) levels of theory. Intrinsic reaction coordinate (IRC) calculations were performed to characterize the transition states on the potential energy surface. Only one pathway was found for the reaction of Ad=Ad with one Br2 producing a bromonium/bromide ion pair. Three mechanisms for the reaction of Ad=Ad with 2Br2 were found, leading to three different structural forms of the bromonium/Br3- ion pair. Activation energies, free energies, and enthalpies of activation along with the relative stability of products for each reaction pathway were calculated. The reaction of Ad=Ad with 2Br2 was strongly favored over the reaction with only one Br2. According to B3LYP/6-31G(d) and single point calculations at MP2, the most stable bromonium/Br3- ion pair would form spontaneously. The most stable of the three bromonium/Br3- ion pairs has a structure very similar to the observed X-ray structure. Free energies of activation and relative stabilities of reactants and products in CCl4 and CH2ClCH2Cl were also calculated with PCM using the united atom (UA0) cavity model and, in general, results similar to the gas phase were obtained. An optimized structure for the trans-1,2-dibromo product was also found at all levels of theory both in gas phase and in solution, but no transition state leading to the trans-1,2-dibromo product was obtained.     

The mean association number: Ultraviolet spectroscopic study of PbCl2−H2O−C2H5OH

The concept of mean asociation number (MAN) has been proposed to describe the association state of an electrolyte. Both the mean spherical approximation (MSA) and the generalized mean spherical approximation (GMSA) are used to calculate the MAN. The theoretical results obtained are compared to the experiment results from an ultraviolet spectroscopic study of PbCl₂–H₂O–C₂H₅OH solutions and also to the predictions of the Bjerrum theory as given by Friedman.     

Theoretical study of HCN(+) + C2H2 reaction

A detailed theoretical investigation for the ion-molecule reaction of HCN (+) with C 2H 2 is performed at the B3LYP/6-311G(d,p) and CCSD(T)/6-311++G(3df,2pd) (single-point) levels. Possible energetically allowed reaction pathways leading to various low-lying dissociation products are probed. It is shown that eight dissociation products P 1 (H 2C 3N (+)+H), P 2 (CN+C 2H 3 (+)), P 3 (HC 3N (+)+H 2), P 4 (HCCCNH (+)+H), P 5 (H 2NCCC (+)+H), P 6 (HCNCCH (+)+H), P 7 (C 2H 2 (+)+HCN), and P 8 (C 2H 2 (+)+HNC) are both thermodynamically and kinetically accessible. Among the eight dissociation products, P 1 is the most abundant product. P 7 and P 3 are the second and third feasible products but much less competitive than P 1 , followed by the almost negligible product P 2 . Other products, P 4 (HCCCNH (+)+H), P 5 (HCNCCH (+)+H), P 6 (H 2NCCC (+)+H), and P 8 (C 2H 2 (+)+HNC) may become feasible at high temperatures. Because the intermediates and transition states involved in the reaction HCN (+) + C 2H 2 are all lower than the reactant in energy, the title reaction is expected to be rapid, as is consistent with the measured large rate constant at room temperature. The present calculation results may provide a useful guide for understanding the mechanism of HCN (+) toward other pi-bonded molecules.     

Isotope effect in Ni + C2H4 (C2D4) association reaction

The association reactions of atomic nickel with ethene and fully deuterated ethene in carbon dioxide buffer gas at 295 K have been investigated in the pressure range 5–100 torr, using a laser photolysis-laser fluorescence technique. By comparison with results of ab initio quantum chemistry calculations for the complex Ni[C₂H₄], the data are shown to be consistent with reaction on both ground and excited state potential energy surfaces. A simple rate equations treatment is described which shows the form of the pressure dependence of the second-order recombination rate coefficient in this case. Under conditions which are expected to hold for the Ni + C₂H₄ (C₂D₄) reaction, the pressure dependence has the standard Lindemann-Hinshelwood form, with the limiting high pressure rate constant given by an apparent value which reflects the degree to which the participating electronic states are coupled by nonadiabatic transitions. The limiting high pressure behavior of the recombination rate coefficient for Ni + C₂H₄ is not strongly affected by deuterium isotope substitution. However, the effect on the low pressure rate constant is large and consistent with RRKM unimolecular reaction theory. This validates the use of RRKM calculations for estimating the binding energy of the complex from kinetic data. The binding energy of Ni[C₂H₄] is estimated to be 35.2 ± 4.2 kcal mol^?1. © 1994 John Wiley & Sons, Inc.     

A shock-tube study of the kinetics of reaction of hydroxyl radicals with H2, CO,CH4, CF3H,C2H4, and C2H6

Concentration-time profiles have been measured for hydroxyl radicals generated by the shock-tube decomposition of hydrogen peroxide in the presence of a variety of additives. At temperatures close to 1300°K the rate constants for the reaction are found to be in the ratio 0.18:0.19:0.59:1.00:2.33:2.88 for the additives CO:CF₃H:H₂:CH₄:C₂H₄:C₂H₆, respectively.     

Radical reaction C3H+NO: a mechanistic study

Although a number of hydrocarbon radicals including the heavier C(3)-radicals C(3)H(3) and C(3)H(5) have been experimentally shown to deplete NO effectively, no theoretical or experimental attempts have been made on the reactivity of the simplest C(3)-radical towards NO. In this article, we report our detailed mechanistic study on the C(3)H+NO reaction at the Gussian-3//B3LYP/6-31G(d) level by constructing the singlet and triplet electronic state [H,C(3),N,O] potential energy surfaces (PESs). The l-C(3)H+NO reaction is shown to barrierlessly form the entrance isomer HCCCNO followed by the direct O-elimination leading to HCCCN+(3)O on triplet PES, or by successive O-transfer, N-insertion, and CN bond-rupture to generate the product (1)HCCN+CO on singlet PES. The possible singlet-triplet intersystem crossings are also discussed. Thus, the novel reaction l-C(3)H+NO can proceed effectively even at low temperatures and is expected to play an important role in both combustion and interstellar processes. For the c-C(3)H+NO reaction, the initially formed H-cCCC-NO can most favorably isomerize to HCCCNO, and further evolution follows that of the l-C(3)H+NO reaction. Quantitatively, the c-C(3)H+NO reaction can take place barrierlessly on singlet PES, yet it faces a small barrier 2.7 kcal/mol on triplet PES. The results will enrich our understanding of the chemistry of the simplest C(3)-radical in both combustion and interstellar processes, which to date have received little attention despite their importance and available abundant studies on its structural and spectroscopic properties.     

The reaction of triplet nitrenes with oxygen: a computational study

Triplet carbenes react much more rapidly with oxygen than do triplet nitrenes. This trend is explained by DFT and MO calculations. [reaction: see text]