The kinetics of the gas phase decomposition reactions of the hexafluoro-t-butoxy radical

Hexafluoro-t-butoxy radicals have been generated by reacting fluorine with hexafluoro-2-methyl isopropanol: Over the temperature range of 406–600 K the hexafluoro-t-butoxy radical decomposes exclusively by loss of a CF₃ radical [reaction (-2)] rather than by loss of a CH₃ radical [reaction (-1)]: (1) The limits of detectability of the product CF₃COCF₃, by gas-chromatographic analysis, place a lower limit on the ratio k_?2/k_-1 of ～80. The implications of this finding in relation to the reverse radical addition reactions to the carbonyl group are briefly discussed. A thermochemical kinetic calculation reveals a discrepancy in the kinetics and thermodynamics of the decomposition and formation reactions of the related t-butoxy radical:      

Ab initio chemical kinetics for the hydrolysis of N2O4 isomers in the gas phase

The mechanism and kinetics for the gas-phase hydrolysis of N(2)O(4) isomers have been investigated at the CCSD(T)/6-311++G(3df,2p)//B3LYP/6-311++G(3df,2p) level of theory in conjunction with statistical rate constant calculations. Calculated results show that the contribution from the commonly assumed redox reaction of sym-N(2)O(4) to the homogeneous gas-phase hydrolysis of NO(2) can be unequivocally ruled out due to the high barrier (37.6 kcal/mol) involved; instead, t-ONONO(2) directly formed by the association of 2NO(2), was found to play the key role in the hydrolysis process. The kinetics for the hydrolysis reaction, 2NO(2) + H(2)O ? HONO + HNO(3) (A) can be quatitatively interpreted by the two step mechanism: 2NO(2) → t-ONONO(2), t-ONONO(2) + H(2)O → HONO + HNO(3). The predicted total forward and reverse rate constants for reaction (A), k(tf) = 5.36 × 10(-50)T(3.95) exp(1825/T) cm(6) molecule(-2) s(-1) and k(tr) = 3.31 × 10(-19)T(2.478) exp(-3199/T) cm(3) molecule(-1) s(-1), respectively, in the temperature range 200-2500 K, are in good agreement with the available experimental data.     

Ab initio chemical kinetics for the unimolecular decomposition of Si2H5 radical and related reverse bimolecular reactions

For plasma enhanced and catalytic chemical vapor deposition (PECVD and Cat‐CVD) processes using small silanes as precursors, disilanyl radical (Si₂H₅) is a potential reactive intermediate involved in various chemical reactions. For modeling and optimization of homogeneous a‐Si:H film growth on large‐area substrates, we have investigated the kinetics and mechanisms for the thermal decomposition of Si₂H₅ producing smaller silicon hydrides including SiH, SiH₂, SiH_3, and Si₂H₄, and the related reverse reactions involving these species by using ab initio molecular‐orbital calculations. The results show that the lowest energy path is the production of SiH + SiH₄ that proceeds via a transition state with a barrier of 33.4 kcal/mol relative to Si₂H₅. Additionally, the dissociation energies for breaking the Si? Si and H? SiH₂ bonds were predicted to be 53.4 and 61.4 kcal/mol, respectively. To validate the predicted enthalpies of reaction, we have evaluated the enthalpies of formation for SiH, SiH₂, HSiSiH₂, and Si₂H₄(C_2h) at 0 K by using the isodesmic reactions, such as ²HSiSiH₂ + ¹C₂H₆→¹Si₂H₆ + ²HCCH₂ and ¹Si₂H₄(C_2h) + ¹C₂H₆ → ¹Si₂H₆ + ¹C₂H₄. The results of SiH (87.2 kcal/mol), SiH₂ (64.9 kcal/mol), HSiSiH₂ (98.0 kcal/mol), and Si₂H₄ (68.9 kcal/mol) agree reasonably well previous published data. Furthermore, the rate constants for the decomposition of Si₂H₅ and the related bimolecular reverse reactions have been predicted and tabulated for different T, P‐conditions with variational Rice–Ramsperger–Kassel–Marcus (RRKM) theory by solving the master equation. The result indicates that the formation of SiH + SiH₄ product pair is most favored in the decomposition as well as in the bimolecular reactions of SiH₂ + SiH₃, HSiSiH₂ + H₂, and Si₂H₄(C_2h) + H under T, P‐conditions typically used in PECVD and Cat‐CVD. © 2013 Wiley Periodicals, Inc.     

Ab initio thermochemistry and kinetics for carbon-centered radical addition and beta-scission reactions

A quantitative comparison of ab initio calculated rate coefficients using five computational methods and five different approaches of treating hindered internal rotation and tunneling with experimental values of rate coefficients for nine carbon-centered radical additions/beta scissions at 300, 600, and 1000 K is performed. The high-accuracy compound methods, CBS-QB3 and G3B3, and the density functionals, MPW1PW91, BB1K, and BMK, have been evaluated using the following approaches: (i) the harmonic oscillator approximation; (ii) the hindered internal rotor approximation for the internal rotation about the forming/breaking bond in the transition state and product; and the hindered internal rotation approximation combined with (iii) Wigner, (iv) Skodje and Truhlar, and (v) Eckart zero-curvature tunneling corrections. The density functional theory (DFT) based values for beta-scission rate coefficients deviate significantly from the experimental ones at 300 K, and the DFT methods do not accurately predict the equilibrium coefficient. The hindered rotor approximation offers a significant improvement in the agreement with experimental rate coefficients as compared to the harmonic oscillator treatment, especially at higher temperatures. Tunneling correction factors are smaller than 1.40 at 300 K and 1.03 at 1000 K. For both the CBS-QB3 method, including the hindered rotor treatment but excluding tunneling corrections, and the G3B3 method, including hindered rotor and Eckart tunneling corrections, a mean factor of deviation with experimentally observed values of 3 is found.     

Ab initio chemical kinetics for reactions of ClO with Cl2O2 isomers

The mechanisms for the reactions of ClO with ClOClO, ClOOCl, and ClClO(2) have been investigated at the CCSD(T)/6-311+G(3df)//PW91PW91∕6-311+G(3df) level of theory. The rate constants for their low energy channels have been calculated by statistical theory. The results show that the main products for the reaction of ClO with ClOClO are ClOCl + ClOO, which can be produced readily by ClO abstracting the terminal O atom from ClOClO. This process occurs without an intrinsic barrier, with the predicted rate constant: k (ClO + ClOClO) = 7.26 × 10(-10) T(-0.15) × exp (-40/T) cm(3)molecule(-1)s(-1) for 200-1500 K. For the reactions of ClO + ClOOCl and ClClO(2), the lowest abstraction barriers are 7.2 and 7.3 kcal/mol, respectively, suggesting that these two reactions are kinetically unimportant in the Earth's stratosphere as their rate constants are less than 10(-14) cm(3)molecule(-1)s(-1) below 700 K. At T = 200-1500 K, the computed rate constants can be represented by k (ClO+ ClOOCl) = 1.11 × 10 (-14) T (0.87) exp (-3576/T) and k (ClO+ ClClO(2)) = 4.61 × 10(-14) T(0.53) exp (-3588/T) cm(3)molecule(-1)s(-1). For these systems, no experimental or theoretical kinetic data are available for comparison.     

Ab initio studies of aspartic acid conformers in gas phase and in solution

Systematic and extensive conformational searches of aspartic acid in gas phase and in solution have been performed. For the gaseous aspartic acid, a total of 1296 trial canonical structures and 216 trial zwitterionic structures were generated by allowing for all combinations of internal single-bond rotamers. All the trial structures were optimized at the B3LYP/6-311G* level and then subjected to further optimization at the B3LYP/6-311++G** level. A total of 139 canonical conformers were found, but no stable zwitterionic structure was found. The rotational constants, dipole moments, zero-point vibrational energies, harmonic frequencies, and vertical ionization energies of the canonical conformers were determined. Single-point energies were also calculated at the MP2/6-311++G** and CCSD/6-311++G** levels. The equilibrium distributions of the gaseous conformers at various temperatures were calculated. The proton affinity and gas phase basicity were calculated and the results are in excellent agreement with the experiments. The conformations in the solution were studied with different solvation models. The 216 trial zwitterionic structures were first optimized at the B3LYP/6-311G* level using the Onsager self-consistent reaction field model (SCRF) and then optimized at the B3LYP/6-311++G** level using the conductorlike polarized continuum model (CPCM) SCRF theory. A total of 22 zwitterions conformers were found. The gaseous canonical conformers were combined with the CPCM model and optimized at the B3LYP/6-311++G** level. The solvated zwitterionic and canonical structures were further examined by the discrete/SCRF model with one and two water molecules. The incremental solvation of the canonical and zwitterionic structures with up to six water molecules in gas phase was systematically examined. The studies show that combining aspartic acid with at least six water molecules in the gas phase or two water molecules and a SCRF solution model is required to provide qualitatively correct results in the solution.     

Ab initio study of the decomposition of 2,5-dimethylfuran

The initial steps in the thermal decomposition of 2,5-dimethylfuran are identified as scission of the C-H bond in the methyl side chain and formation of β- and α-carbenes via 3,2-H and 2,3-methyl shifts, respectively. A variety of channels are explored which prise the aromatic ring open and lead to a number of intermediates whose basic properties are essentially unknown. Once the furan ring is opened demethylation to yield highly unsaturated species such as allenylketenes appears to be a feature of this chemistry. The energetics of H abstraction by the hydroxyl radical (and other abstracting species) from a number of mono- and disubstituted methyl furans has been studied. H-atom addition to 2,5-dimethylfuran followed by methyl elimination is shown to be the most important route to formation of the less reactive 2-methylfuran. Identification of 2-ethenylfuran as an C(6)H(6)O intermediate in 2,5-dimethylfuran flames is probably not correct and is more likely the isomeric 2,5-dimethylene-2,5-dihydrofuran for which credible formation channels exist.     

Ab initio study of 1,4-pentadienyl electrocyclic reactions

The thermochemistry and transition states of the electrocyclic ring closures of the resonance-stabilized 1,4-pentadienyl radical to cyclopenten-3-yl, cyclobut-2-enylmethyl, and 2-vinylcyclopropyl are investigated at Hartree-Fock and coupled-cluster levels of theory. The CCSD(T)//QCISD/cc-pVDZ calculations predict activation barriers of 130, 169, and 236 kJ/mol, respectively, and DeltaH values of -60, 115, and 155 kJ/mol. Experimental evidence for the appearance of vinylcyclopropyl following photolytic generation of pentadienyl is more likely the result of a distinct electrocyclic reaction than quenching of a two-step mechanism for formation of cyclopentenyl. Higher energy pathways for formation of polycyclic structures are also briefly examined.     

Ab initio studies on the pyrolysis mechanism of 2-haloacetic acids in the gas phase

The dehydrohalogenation mechanism of 2-haloacetic acids (XCH2CO2H, X=F, Cl and Br) has been studied theoretically by HF/3-21G and AM1 methods. The results indicate that these reactions are most probably proceeded in terms of a polar five-membered cyclic transition state in the gas phase. Their microscopic processes are beleived to be a stepwise reaction and the rate-determining step is the first one. By comparing the energy barriers of different 2-haloacetic acids, it can be realized that 2-fluoroacetic acid is easier to react than 2-chloroacetic and 2-bromoacetic acids.     

Ab initio study on tautomerism of 2-thiouracil in the gas phase and in solution

Ab initio geometry optimizations were carried out at the HF/3-21G and HF/6-31+G^** levels for the six tautomeric forms of 2-thiouracil (2TU, 2TU1, 2TU2, 2TU3, 2TU4, 2TU5) in the gas phase and in solution. To obtain a more definitive estimate of the relative stabilities for 2-thiouracil tautomers in the gas phase, single-point MP2/6-31+G^** calculations were performed on the HF/6-31+G^** optimized geometries. The tautomeric equilibria in 1,4-dioxane (=2.21), acetonitrile (=38), and in water (=78.54) were studied using the self-consistent reaction field (SCRF) theory. The calculated relative free energies indicated that 2TU is the energetically preferred tautomer in the gas phase and in solution. The stability order of 2-thiouracil tautomers depends on the level of theory and the dielectric constant of the solvent. The obtained results are compared with the available experimental data.     

Ab initio chemical kinetics for reactions of H atoms with SiHx (x = 1–3) radicals and related unimolecular decomposition processes

Hydrogen atoms and SiH_x (x = 1–3) radicals coexist during the chemical vapor deposition (CVD) of hydrogenated amorphous silicon (a‐Si:H) thin films for Si‐solar cell fabrication, a technology necessitated recently by the need for energy and material conservation. The kinetics and mechanisms for H‐atom reactions with SiH_x radicals and the thermal decomposition of their intermediates have been investigated by using a high high‐level ab initio molecular‐orbital CCSD (Coupled Cluster with Single and Double)(T)/CBS (complete basis set extrapolation) method. These reactions occurring primarily by association producing excited intermediates, ¹SiH₂, ³SiH₂, SiH₃, and SiH₄, with no intrinsic barriers were computed to have 75.6, 55.0, 68.5, and 90.2 kcal/mol association energies for x = 1–3, respectively, based on the computed heats of formation of these radicals. The excited intermediates can further fragment by H₂ elimination with 62.5, 44.3, 47.5, and 56.7 kcal/mol barriers giving ¹Si, ³Si, SiH, and ¹SiH₂ from the above respective intermediates. The predicted heats of reaction and enthalpies of formation of the radicals at 0 K, including the latter evaluated by the isodesmic reactions, SiH_x + CH₄ = SiH₄ + CH_x, are in good agreement with available experimental data within reported errors. Furthermore, the rate constants for the forward and unimolecular reactions have been predicted with tunneling corrections using transition state theory (for direct abstraction) and variational Rice–Ramsperger–Kassel–Marcus theory (for association/decomposition) by solving the master equation covering the P,T‐conditions commonly employed used in industrial CVD processes. The predicted results compare well experimental and/or computational data available in the literature. © 2013 Wiley Periodicals, Inc.     

Ab initio study on the thermolysis of β-hydroxyolefins in gas phase

Thermolysis studies of β-hydroxyolefins in gas phase were realized using ab initio MP2 and DFT methods at the 6-31G^* levels to explore the possibility of determining a possible concerted process with a six-membered cyclic transition state (TS). Vibrational frequency calculations were carried out in order to confirm the stationary states, including TS structures. IRC calculations have been performed in all cases in order to verify that localized TS structures connect with the corresponding minimum stationary points associated with the reactant and products. With the aim of corroborating the postulated mechanism in the experimental study, we present a theoretical study in order to calculate the rate constants and the activation parameters. The results obtained are in accordance with the experimental conclusions.     

Ab initio molecular orbital studies of nonidentity allyl transfer reactions

Ab initio molecular orbital (MO) calculations are carried out on the nonidentity allyl transfer processes, X^? + CH₂CHCH₂Y ? CH₂CHCH₂ X + Y^?, with X^? = H, F, and Cl and Y = H, NH₂, OH, F, PH₂, SH, and Cl. The Marcus equation applies well to the allyl transfer reactions. The transition state (TS) position along the reaction coordinate and the TS structure are strongly influenced by the thermodynamic driving force, whereas the TS looseness is originated from the intrinsic barrier. The intrinsic barrier, ΔE, looseness, %L?, and absolute asymmetry, %AS?, are well correlated with the percentage bond elongation, %CY? = [(d ? d)/d] × 100 and/or %CX?. The %CY? and the bond orders indicate that a stronger nucleophile and/or a stronger nucleofuge (or a better leaving group) leads to an earlier TS on the reaction coordinate with a lesser degree of bond making as well as bond breaking. These are consistent with the Bell-Evans-Polanyi principle and the Leffler-Hammond postulate. © 1995 by John Wiley & Sons, Inc.     

Ab initio study of C + H3+ reactions

The reaction C + H3+ --> CH(+) + H2 is frequently used in models of dense interstellar cloud chemistry with the assumption that it is fast, i.e. there are no potential energy barriers inhibiting it. Ab initio molecular orbital study of the triplet CH3+ potential energy surface (triplet because the reactant carbon atom is a ground state triplet) supports this hypothesis. The reaction product is 3 pi CH+; the reaction is to exothermic even though the product is not in its electronic ground state. No path has been found on the potential energy surface for C + H3+ --> CH2(+) + H reaction.     

Ab initio chemical kinetics for the OH+HNCN reaction

The kinetics and mechanism of the reaction of the cyanomidyl radical (HNCN) with the hydroxyl radical (OH) have been investigated by ab initio calculations with rate constants prediction. The single and triplet potential energy surfaces of this reaction have been calculated by single-point calculations at the CCSD(T)/6-311+G(3df,2p) level based on geometries optimized at the B3LYP/6-311+G(3df,2p) and CCSD/6-311++G(d,p) levels. The rate constants for various product channels in the temperature range of 300-3000 K are predicted by variational transition-state and Rice-Ramsperger-Kassel-Marcus (RRKM) theories. The predicted total rate constants can be represented by the expressions ktotal=2.66 x 10(+2)xT-4.50 exp(-239/T) in which T=300-1000 K and 1.38x10(-20)xT2.78 exp(1578/T) cm3 molecule(-1) s(-1) where T=1000-3000 K. The branching ratios of primary channels are predicted: k1 for forming singlet HON(H)CN accounts for 0.32-0.28, and k4 for forming singlet HONCNH accounts for 0.68-0.17 in the temperature range of 300-800 K. k2+k7 for producing H2O+NCN accounts for 0.55-0.99 in the high-temperature range of 800-3000 K. The branching ratios of k3 for producing HCN+HNO, k6 for producing H2N+NCO, k8 for forming 3HN(OH)CN, k9 for producing CNOH+3NH, and k5+k10 for producing NH2+NCO are negligible. The rate constants for key individual product channels are provided in a table for different temperature and pressure conditions.     

Thermal decomposition of ethanol. 4. Ab initio chemical kinetics for reactions of H atoms with CH3CH2O and CH3CHOH radicals

The potential energy surfaces of H-atom reactions with CH(3)CH(2)O and CH(3)CHOH, two major radicals in the decomposition and oxidation of ethanol, have been studied at the CCSD(T)/6-311+G(3df,2p) level of theory with geometric optimization carried out at the BH&HLYP/6-311+G(3df,2p) level. The direct hydrogen abstraction channels and the indirect association/decomposition channels from the chemically activated ethanol molecule have been considered for both reactions. The rate constants for both reactions have been calculated at 100-3000 K and 10(-4) Torr to 10(3) atm Ar pressure by microcanonical VTST/RRKM theory with master equation solution for all accessible product channels. The results show that the major product channel of the CH(3)CH(2)O + H reaction is CH(3) + CH(2)OH under atmospheric pressure conditions. Only at high pressure and low temperature, the rate constant for CH(3)CH(2)OH formation by collisonal deactivation becomes dominant. For CH(3)CHOH + H, there are three major product channels; at high temperatures, CH(3)+CH(2)OH production predominates at low pressures (P < 100 Torr), while the formation of CH(3)CH(2)OH by collisional deactivation becomes competitive at high pressures and low temperatures (T < 500 K). At high temperatures, the direct hydrogen abstraction reaction producing CH(2)CHOH + H(2) becomes dominant. Rate constants for all accessible product channels in both systems have been predicted and tabulated for modeling applications. The predicted value for CH(3)CHOH + H at 295 K and 1 Torr pressure agrees closely with available experimental data. For practical modeling applications, the rate constants for the thermal unimolecular decomposition of ethanol giving key accessible products have been predicted; those for the two major product channels taking place by dehydration and C-C breaking agree closely with available literature data.     

The kinetics and mechanism of the shock induced gas phase decomposition of ethylsilane

The decomposition kinetics of ethylsilane under shock tube conditions (P_T ca. 3100 torr, T ? 1080–1245 K), both in the absence and presence of silylene trapping agents (butadiene and acetylene) are reported. Arrhenius parameters under maximum butadiene inhibition are: log k(C₂H₅SiH₃) = 15.14-64,769 ± 1433 cal/2.303 RT; log k(C₂H₅SiD₃) = 15.29-66,206 ± 1414/2.303 RT. The uninhibited reaction is subject to silylene induced decomposition (63% lowest T -- 24% highest T). Major reaction products are ethylene and hydrogen, consistent with two dominant primary dissociation reactions: C₂H₅SiD₃ → C₂H₅SiD + D₂, ? ? 0.66; C₂H₅SiD₃ → CH₃CH = SiD₂ + HD, ? ? 0.30. Minor products suggest several other less important primary processes: alkane elimination, ? ?0.02, and free-radical production via simple bond fission, ? ?0.02. An upper limit for the activation energy of the decomposition, C₂H₅SiH → C₂H₄ + SiH₂, of E < 30 ± 4 kcal is established, and speculations on the mechanism of this decomposition (concerted or stepwise) with conclusions in favor of the stepwise path are made. Computer modeling studies for the reaction both in the absence and presence of butadiene are shown to be in good agreement with the experimental observations.     

Ab initio computation of low-temperature phase diagrams exhibiting miscibility gaps

A new methodology for the computation of the low-temperature part of phase diagrams without recourse to any experimental information is presented. A central element is a procedure for deciding whether formation of crystalline solid solution phases can take place in the chemical system. Via global exploration of the enthalpy landscapes for many different compositions in the system, candidates for ordered stoichiometric and crystalline solid solution phases are identified. Next, their free enthalpies are computed at ab initio level and a low-temperature phase diagram is derived. As examples, the low-temperature phase diagrams for the ternary alkali halides NaCl/LiCl NaBr/LiBr and NaCl/KCl are presented.