Ab Initio Study on the Oxidation of NCN by O (3P): Prediction of the Total Rate Constant and Product Branching Ratios

The reaction of NCN with O is relevant to the formation of prompt NO according to the new mechanism, CH+N2-->cyclic-C(H)NN- -->HNCN-->H+NCN. The reaction has been investigated by ab initio molecular orbital and transition state theory calculations. The mechanisms for formation of possible product channels involved in the singlet and triplet potential energy surfaces have been predicted at the highest level of the modified GAUSSIAN-2 (G2M) method, G2M (CC1). The barrierless association/dissociation processes on the singlet surface were also examined with the third-order Rayleigh-Schrödinger perturbation (CASPT3) and the multireference configuration interaction methods including Davidson's correction for higher excitations (MRCI+Q) at the CASPT3(6,6)/6-311+G(3df)//UB3LYP/6-311G(d) and MRCI+Q(6,6)/6-311+G(3df)//UB3LYP/6-311G(d) levels. The rate constants for the low-energy channels producing CO+N2, CN+NO, and N(4S)+NCO have been calculated in the temperature range of 200-3000 K. The results show that the formation of CN+NO is dominant and its branching ratio is over 99% in the whole temperature range; no pressure dependence was noted at pressures below 100 atm. The total rate constant can be expressed by: kt=4.23x10(-11) T0.15 exp(17/T) cm3 molecule(-1) s(-1).

Adsorption Configurations and Energetics of BClx (x = 0−3) on TiO2 Anatase (101) and Rutile (110) Surfaces

This study investigates the adsorption and reactions of boron trichloride and its fragments (BClx) on the TiO2 anatase (101) and rutile (110) surfaces by first-principles calculations. The results show that the possible absorbates on the TiO2 anatase and rutile surfaces are very similar. The single- and double-site adsorption configurations are found for both anatase and rutile surfaces. The particular adsorbate feature on the anatase surface is its in-plane double-site adsorption by Ti and O from its sawtooth surface. The potential energy surface shows that BCl3 can be adsorbed on the O site for both the anantase and rutile surfaces and the most of the BClx reaction on both anatase and rutile surfaces are endothermic, except for the dissociative reaction on the rutile surface. The energy levels of the BClx reactions between the anatase and rutile surfaces show that the rutile surface has lower energy levels than those of anatase surface. This result reveals that the BClx dissociative adsorption more easily occurs on rutile surface than on anatase surface.

Experimental and Computational Studies of the Phenyl Radical Reaction with Allene

The kinetics for the gas-phase reaction of phenyl radicals with allene has been measured by cavity ring-down spectrometry (CRDS), and the mechanism and initial product branching have been elucidated with the help of quantum-chemical calculations. The absolute rate constant measured by the CRDS technique can be expressed by the following Arrhenius equation: kallene (T=301-421 K)=(4.07+/-0.38)x10(11) exp[-(1865+/-85)/T] cm3 mol(-1) s(-1). Theoretical calculations, employing high level G2M energetic and IRCMax(RCCSD(T)//B3LYP-DFT) molecular parameters, indicate that under our experimental conditions the most preferable reaction channel is the addition of phenyl radicals to the terminal carbon atoms in allene. Predicted total rate constants agree with the experimental values within 40%. Calculated total and branching rate constants are provided for high-T kinetic modeling.

Computational Study of Reaction Pathways for the Formation of Indium Nitride from Trimethylindium with HN3: Comparison of the Reaction with NH3 and That on TiO2 Rutile (110) Surface

The reactions of trimethylindium (TMIn) with HN3 and NH3 are relevant to the chemical vapor deposition of indium nitride thin film. The mechanisms and energetics of these reactions in the gas phase have been investigated by density functional theory and ab initio calculations using the CCSD(T)/Lanl2dz//B3LYP/Lanl2dz and CCSD(T)/Lanl2dz//MP2/Lanl2dz methods. The results of both methods are in good agreement for the optimized geometries and relative energies. These results suggest that the reaction with HN3 forms a new stable product, dimethylindiumnitride, CH3-In=N-CH3 via another stable In(CH3)2N3 (dimethylindium azide, DMInA) intermediate. DMInA may undergo unimolecular decomposition to form CH3InNCH3 by two main possible pathways: (1) a stepwise decomposition process through N2 elimination followed by CH3 migration from In to the remaining N atom and (2) a concerted process involving the concurrent CH3 migration and N2 elimination directly giving N2+CH3InNCH3. The reaction of TMIn with NH3 forms a most stable product DMInNH2 following the initial association and CH4-elimination reaction. The required energy barrier for the elimination of the second CH4 molecule from DMInNH2 is 74.2 kcal/mol. Using these reactions, we predict the heats of formation at 0 K for all the products and finally for InN which is 123+/-1 kcal/mol predicted by the two methods. The gas-phase reaction of HN3 with TMIn is compared with that occurring on rutile TiO2 (110). The most noticeable difference is the high endothermicity of the gas-phase reaction for InN production (53 kcal/mol) and the contrasting large exothermicity (195 kcal/mol) released by the low-barrier Langmuir-Hinshelwood type processes following the adsorption of TMIn and HN3 on the surface producing a horizontally adsorbed InN(a), Ti-NIn-O(a), and other products, CH4(g)+N2(g)+2CH3O(a) [J. Phys. Chem. B 2006, 110, 2263].

The mechanisms of the reactions of W and W+ with COx (x=1, 2): a computational study

The mechanisms of the reactions of W and W+ with COx (x=1, 2) were studied at the CCSD(T)/[SDD+6-311G(d)]//B3LYP/[SDD+6-31G(d)] level of theory. It was shown that the gas-phase reaction of W with CO2 proceeds with a negligible barrier via an insertion pathway, W(7S)+CO2(1A1)-->W(eta2-OCO)(6A')-->OW(eta1-CO)(1A)-->WO (3Sigma+)+CO(1Sigma). This oxidation process is calculated to be exothermic by 32.4 kcal/mol. Possible intermediates of this reaction are the W(eta2-OCO) and OWCO complexes, among which the latter is 37.4 kcal/mol more stable and lies 39.7 and 7.3 kcal/mol lower than the reactants, W(7S)+CO2(1A1), and the products, WO (3Sigma+)+CO(1Sigma), respectively. The barrier separating W(eta2-OCO) from OWCO is 8.0 kcal/mol (relative to the W(eta2-OCO) complex), which may be characterized as a W+delta-(CO2)-delta charge-transfer complex. Ionization of W does not change the character of the reaction of W with CO2: the reaction of W+ with CO2, like its neutral analog, proceeds via an insertion pathway and leads to oxidation of the W-center. The overall reaction W+(6D) + CO2(1A1)-->W(eta1-OCO)+(6A)-->OW(eta1-CO)+(4A)-->WO+(4Sigma+)+CO(1Sigma) is calculated to be exothermic by 25.4 kcal/mol. The cationic reaction proceeds with a somewhat large (9.9 kcal/mol) barrier and produces two intermediates, W(eta1-OCO)+(6A) and OW(eta1-CO)+(4A). Intermediate W(eta1-OCO)+(6A) is 20.0 kcal/mol less stable than OW(eta1-CO)+(4A), and separated from the latter by a 35.2 kcal/mol barrier. Complex W(eta1-OCO)+(6A) is characterized as an ion-molecular complex type of W+-(CO2). Gas-phase reactions of M=W/W+ with CO lead to the formation of a W-carbonyl complex M(eta1-CO) for both M=W and W+. The C-O insertion product, OMC, lies by 5.2 and 69.3 kcal/mol higher than the corresponding M(eta1-CO) isomer, for M=W and W+, respectively, and is separated from the latter by a large energy barrier.

Computational Study on the Kinetics and Mechanisms for the Unimolecular Decomposition of Formic and Oxalic Acids

The kinetics and mechanisms for the unimolecular decomposition reactions of formic acid and oxalic acid have been studied computationally by the high-level G2M(CC1) method and microcanonical RRKM theory. There are two reaction pathways in the decomposition of formic acid: The dehydration process starting from the Z conformer is found to be the dominant, whereas the decarboxylation reaction starting from the E conformer is less competitive. The predicted rate constants for the dehydration and decarboxylation reactions are in good agreement with the experimental data. The calculated CO/CO2 ratio, 13.6-13.9 between 1300 and 2000 K, is in close agreement with the ratio of 10 measured experimentally by Hsu et al. (In The 19th International Symposium on Combustion; The Combustion Institute: Pittsburgh, PA, 1983; p 89). For oxalic acid, its isomer with two intramolecular hydrogen bonds is the most stable structure, similar to earlier reports. Two primary decomposition channels of oxalic acid producing CO2+HOCOH with barriers of 33-36 kcal/mol and CO2+CO+H2O with a barrier of 39 kcal/mol were found. At high temperatures, the latter process becomes more competitive. The rate constant predicted for the formation of CO2 and HOCOH (the precursor of HCOOH) agrees well with available experimental data. The mechanism for the isomerization of HOCOH to HCOOH is also discussed.

Photodissociation dynamics of nitrobenzene and o-nitrotoluene

Photodissociation of nitrobenzene at 193, 248, and 266 nm and o-nitrotoluene at 193 and 248 nm was investigated separately using multimass ion imaging techniques. Fragments corresponding to NO and NO(2) elimination from both nitrobenzene and o-nitrotoluene were observed. The translational energy distributions for the NO elimination channel show bimodal distributions, indicating two dissociation mechanisms involved in the dissociation process. The branching ratios between NO and NO(2) elimination channels were determined to be NONO(2)=0.32+/-0.12 (193 nm), 0.26+/-0.12 (248 nm), and 0.4+/-0.12(266 nm) for nitrobenzene and 0.42+/-0.12(193 nm) and 0.3+/-0.12 (248 nm) for o-nitrotoluene. Additional dissociation channels, O atom elimination from nitrobenzene, and OH elimination from o-nitrotoluene, were observed. New dissociation mechanisms were proposed, and the results are compared with potential energy surfaces obtained from ab initio calculations. Observed absorption bands of photodissociation are assigned by the assistance of the ab initio calculations for the relative energies of the triplet excited states and the vertical excitation energies of the singlet and triplet excited states of nitrobenzene and o-nitrotoluene. Finally, the dissociation rates and lifetimes of photodissociation of nitrobenzene and o-nitrotoluene were predicted and compared to experimental results.

Mechanisms of the Reactions of W and W+ with NOx (x = 1, 2): A Computational Study

The mechanisms of the reactions of W and W+ with NOx (x=1, 2) were studied at the CCSD(T)/[SDD+6-311G(d)]//B3LYP/[SDD+6-31G(d)] level of theory. It was shown that the insertion pathway of the reaction W(7S)+NO2(2A1) is a multistate process, which involves several lower lying electronic states of numerous intermediates and transition states, and leads to oxidation, WO(3Sigma)+NO(2Pi), and/or nitration, WN(4Sigma)+O2(3Sigmag-), of the W-center. Oxidation products WO(3Sigma)+NO(2Pi) lie 87.6 kcal/mol below the reactants, while the nitration channel is only 31.0 kcal/mol exothermic. Furthermore, it was shown that nitration of W with NO2 is kinetically less favorable than its oxidation. The addition-dissociation pathway of the reaction W(7S)+NO2(2A1) proceeds via the octet (ground) state potential energy surface of the reaction, requires 3.3 kcal/mol barrier, and leads exclusively to oxidation products. Calculations show that oxidation of the W+ cation by NO2 is a barrierless process in the gas phase, proceeds exclusively via the insertion pathway, and is exothermic by 82.9 kcal/mol. The nitration of W+ by NO2 is only 14.1 kcal/mol exothermic and could be accessible only under high-temperature conditions. Reactions of M=W/W+ with NO are also barrierless processes in the gas phase and lead to the N-O insertion product NMO, which are 105.4 and 77.4 kcal/mol lower than the reactants for W and W+, respectively.

Experimental and theoretical studies of rate coefficients for the reaction O(3P)+CH3OH at high temperatures

Rate coefficients of the reaction O((3)P) + CH(3)OH in the temperature range of 835-1777 K were determined using a diaphragmless shock tube. O atoms were generated by photolysis of SO(2) with a KrF excimer laser at 248 nm or an ArF excimer laser at 193 nm; their concentrations were monitored via atomic resonance absorption excited by emission from a microwave-discharged mixture of O(2) and He. The rate coefficients determined for the temperature range can be represented by the Arrhenius equation, k(T) = (2.29 +/- 0.18) x 10(-10) exp[-(4210 +/- 100)T] cm(3) molecule(-1) s(-1); unless otherwise noted, all the listed errors represent one standard deviation in fitting. Combination of these and previous data at lower temperature shows a non-Arrhenius behavior described as the three-parameter equation, k(T) = (2.74 +/- 0.07) x 10(-18)T(2.25 +/- 0.13) exp[-(1500 +/- 90)T] cm(3)molecule(-1) s(-1). Theoretical calculations at the Becke-3-Lee-Yang-Parr (B3LYP)6-311 + G(3df,2p) level locate three transition states. Based on the energies computed with coupled clusters singles, doubles (triples) [CCSD(T)]/6-311 + G(3df,2p)B3LYP6-311 + G(3df,2p), the rate coefficients predicted with canonical variational transition state theory with small curvature tunneling corrections agree satisfactorily with the experimental observations. The branching ratios of two accessible reaction channels forming OH + CH(2)OH (1a) and OH + CH(3)O (1b) are predicted to vary strongly with temperature. At 300 K, reaction (1a) dominates, whereas reaction (1b) becomes more important than reaction (1a) above 1700 K.

Computational Studies on the Kinetics and Mechanisms for NH3 Reactions with ClOx (x = 0−4) Radicals

Kinetics and mechanisms for NH3 reactions with ClOx (x = 0-4) radicals have been investigated at the G2M level of theory in conjunction with statistical theory calculations. The geometric parameters of the species and stationary points involved in the reactions have been optimized at the B3LYP/6-311+G(3df,2p) level of theory. Their energetics have been further refined with the G2M method. The results show that the H-abstraction process is the most favorable channel in each reaction and the barriers predicted in decreasing order are OClO > ClO > Cl > ClO3 > ClO4. All reactions were found to occur by hydrogen-bonding complexes; the rate constants for these complex metathetical processes have been calculated in the temperature range 200-2000 K by the microcanonical VTST and/or RRKM theory (for ClO4 + NH3) with Eckart tunneling and multiple reflection corrections. The predicted rate constants are in good agreement with the available experimental data.

Experimental and theoretical investigation of rate coefficients of the reaction S(P3)+OCS in the temperature range of 298–985K

The reaction S(3P)+OCS in Ar was investigated over the pressure range of 50-710 Torr and the temperature range of 298-985 K with the laser photolysis technique. S atoms were generated by photolysis of OCS with light at 248 nm from a KrF excimer laser; their concentration was monitored via resonance fluorescence excited by atomic emission of S produced from microwave-discharged SO2. At pressures less than 250 Torr, our measurements give k(298 K)=(2.7+/-0.5)x10(-15) cm3 molecule-1 s-1, in satisfactory agreement with a previous report by Klemm and Davis [J. Phys. Chem. 78, 1137 (1974)]. New data determined for 407-985 K connect rate coefficients reported previously for T>or=860 and T<or=478 K and show a non-Arrhenius behavior. Combining our results with data reported at high temperatures, we derived an expression k(T)=(6.1+/-0.3)x10(-18) T1.97+/-0.24 exp[-(1560+/-170)/T] cm3 molecule-1 s-1 for 298<or=TK<or=1680. At 298 K and P>or=500 Torr, the reaction rate was enhanced. Theoretical calculations at the G2M(CC2) level, using geometries optimized with the B3LYP6-311+G(3df) method, yield energies of transition states and products relative to those of the reactants. Rate coefficients predicted with multichannel Rice-Ramsperger-Kassel-Marcus (RRKM) calculations agree satisfactorily with experimental observations. According to our calculations, the singlet channel involving formation of SSCO followed by direct dissociation into S2(a 1Deltag)+CO dominates below 2000 K; SSCO is formed via intersystem crossing from the triplet surface. At low temperature and under high pressure the stabilization of OCS2, formed via isomerization of SSCO, becomes important; its formation and further reaction with S atoms partially account for the observed increase in the rate coefficient under such conditions.

Quantum Chemical Prediction of Reaction Pathways and Rate Constants for Dissociative Adsorption of COx and NOx on the Graphite (0001) Surface

We present predictions of reaction rate constants for dissociative adsorption reactions of CO(x) (x = 1, 2) and NO(x) (x = 1, 2) molecules on the basal graphite (0001) surface based on potential energy surfaces (PES) obtained by the integrated ONIOM(B3LYP:DFTB-D) quantum chemical hybrid approach with dispersion-augmented density functional tight binding (DFTB-D) as low level method. Following an a priori methodology developed in a previous investigation of water dissociative adsorption reactions on graphite, we used a C(94)H(24) dicircumcoronene graphene slab as model system for the graphite surface in finite-size molecular structure investigations, and single adsorbate molecules reacting with the pristine graphene sheet. By employing the ONIOM PES information in RRKM theory we predict reaction rate constants in the temperature range between 1,000 and 5,000 K. We find that among CO(x) and NO(x) adsorbate species, the dissociative adsorption reactions of CO(2) and both radical species NO and NO(2) are likely candidates as a cause for high temperature oxidation and erosion of graphite (0001) surfaces, whereas reaction with CO is not likely to lead to long-lived surface defects. High temperature quantum chemical molecular dynamics simulations (QM/MD) at T = 5,000 K using on-the-fly DFTB-D energies and gradients confirm the results of our PES study.

Characterization of O2–CeO2 Interactions Using In Situ Raman Spectroscopy and First-Principle Calculations

Interactions between O(2) and CeO(2) are examined experimentally using in situ Raman spectroscopy and theoretically using density-functional slab-model calculations. Two distinct oxygen bands appear at 825 and 1131 cm(-1), corresponding to peroxo- and superoxo-like species, respectively, when partially reduced CeO(2) is exposed to 10 % O(2). Periodic density-functional theory (DFT) calculations aid the interpretation of spectroscopic observations and provide energetic and geometric information for the dioxygen species adsorbed on CeO(2). The O(2) adsorption energies on unreduced CeO(2) surfaces are endothermic (0.91<DeltaE(ads)<0.98 eV), while those on reduced surfaces are exothermic (-4. 0<DeltaE(ads)<-0.9 eV), depending on other relevant surface processes such as chemisorption and diffusion into the bulk. Partial reduction of surface Ce(4+) to Ce(3+) (together with formation of oxygen vacancies) alters geometrical parameters and, accordingly, leads to a shift in the vibrational frequencies of adsorbed oxygen species compared to those on unreduced CeO(2). Moreover, the location of oxygen vacancies affects the formation and subsequent dissociation of oxygen species on the surfaces. DFT predictions of the energetics support the experimental observation that the reduced surfaces are energetically more favorable than the unreduced surfaces for oxygen adsorption and reduction.

A Computational Study on the Kinetics and Mechanism for the Unimolecular Decomposition of o-Nitrotoluene

The kinetics and mechanism for the unimolecular decomposition of o-nitrotoluene (o-CH(3)C(6)H(4)NO(2)) have been studied computationally at the G2M(RCC, MP2)//B3LYP/6-311G(d, p) level of theory in conjunction with rate constant predictions with RRKM and TST calculations. The results of the calculations reveal 10 decomposition channels for o-nitrotoluene and its six isomeric intermediates, among them four channels give major products: CH(3)C(6)H(4) + NO(2), C(6)H(4)C(H)ON (anthranil) + H(2)O, CH(3)C(6)H(4)O (o-methyl phenoxy) + NO, and C(6)H(4)C(H(2))NO + OH. The predicted rate constants in the 500-2000 K temperature range indicate that anthranil production, taking place initially by intramolecular H-abstraction from the CH(3) group by NO(2) followed by five-membered ring formation and dehydration, dominates at temperatures below 1000 K, whereas NO(2) elimination becomes predominant above 1100 K and CH(3)C(6)H(4)O formation by the nitro-nitrite isomerization/decomposition process accounts for only 5-11% of the total product yield in the middle temperature range 800-1300 K. The branching ratio for CH(2)C(6)H(4)NO formation by the decomposition process of CH(2)C(6)H(4)N(O)OH is negligible. The predicted high-pressure-limit rate constants with the rate expression of 4.10 x 10(17) exp[-37000/T] s(-1) for the NO(2) elimination channel and 9.09 x 10(12) exp[-25800/T] s(-1) for the H(2)O elimination channel generally agree reasonably with available experimental data. The predicted high-pressure-limit rate constants for the NO and OH elimination channels are represented as 1.49 x 10(14) exp[-30000/T] and 1.31 x 10(15) exp[-38000/T] s(-1), respectively.

Experimental and Theoretical Studies of the Phenyl Radical Reaction with Propene

The kinetics for the reaction of C6H5 with propene has been measured by cavity ring-down spectrometry (CRDS) at temperatures 296-496 K under an Ar pressure of 40 Torr. The total rate constant can be given by the following Arrhenius expression (in units of cm3 mol(-1) s(-1)): k(C6H5 + C3H6) = 10(11.93+/-0.06) exp[-(1512 +/- 51)/T]. Density functional and higher level of theory calculations (up to the G2M level) have been carried out to provide additional insights about the mechanism of this reaction, and we also performed transition state theory (TST) calculation for the rate constant prediction. Our theoretical kinetic calculations predict that the C6H5 addition to the terminal =CH2 site in propene is dominant at the temperature range of our CRDS measurements. However, the H-abstraction channel forming benzene and the allyl radical becomes increasingly important at higher temperatures. The total high-pressure limiting rate constant calculated on the basis of the G2M reaction barriers is in reasonable agreement with the experimental values.

Theoretical study on the kinetics and mechanism for the reaction of FCO with NO

The radical reaction mechanism of FCO + NO on the ground electronic state energy surface has been studied at the G2M level of theory based on the geometric parameters optimized at the B3LYP/6-311+G(d) level of theory. The two kinds of reaction pathways include the direct fluorine abstraction channel producing CO + FNO and the association channel forming the FC(O)NO complex. The former has a distinct barrier of 8.9 kcal mol(-1), while the latter is a barrierless association process. The rate constant of this reaction system in the temperature range 200-3000 K has been calculated by the microcanonical VTST/RRKM theory. The theoretical result shows that the predicted total rate constants exhibit a negative-temperature dependence and positive-pressure effect at lower temperatures. Under the experimental conditions, the predicted values are in good agreement with the experimental results. In addition, the predicted branching ratios clearly indicate that the dominant product channel is the formation of FC(O)NO at low temperatures and FNO + CO at high temperatures (>500 K).

Reactions of Hydrazoic Acid and Trimethylindium on TiO2 Rutile (110) Surface: A Computational Study on the Formation of the First Monolayer InN

This article reports the result of a computational study on the reaction of hydrazoic acid and trimethylindium (TMIn), coadsorbed on TiO2 rutile (110) surface. The adsorption geometries and energies of possible adsorbates including HN3-In(CH3)3(a) and its derivatives, HN3-In(CH3)2(a), N3-In(CH3)2(a), N3-In(CH3)(a), and N-In(a), have been predicted by first-principles calculations based on the density functional theory (DFT) and the pseudopotential method. The mechanisms of these surface reactions have also been explicitly elucidated with the computed potential energy surfaces. Starting from the interaction of three stable HN3 adsorbates, HN3-Ob(a), H(N2)N-Ob(a), and Ti-NN(H)N-Ob(a), where Ob is the bridged O site on the surface, with two stable intermediates from the adsorption and dissociative adsorption of TMIn, (H3C)3In-Ob(a) and (H3C)2In-Ob(a)+H3C-Ob(a), InN products can be formed exothermically via four reaction paths following the initial barrierless In-atom association with the N atom directly bonded to H, by CH4 elimination (with approximately 40 kcal/mol barriers), the InN-N bond breaking and the final CH3 elimination or migration (with <20 kcal/mol barriers). These Langmuir-Hinshelwood processes producing the two most stable InN(a) side-on adsorptions confirm that HN3 and TMIn are indeed very efficient precursors for the deposition of InN films on TiO2 nanoparticles. The result of similar calculations for the reactions occurring by the Rideal-Eley mechanism involving HN3(a)+TMIn(g) and HN3(g)+TMIn(a) indicates that they are energetically less favored and produce the less stable InN(a) with end-on configurations.

Reactions of Trimethylindium on TiO2 Nanoparticles: Experimental and Computational Study

This article reports the results of an experimental and computational study on the reaction of trimethylindium, In(CH3)(3), adsorbed on TiO2 nanoparticle films. Experimentally, Fourier transform infrared (FTIR) spectra have been measured by varying In(CH3)(3) dosing pressure, UV irradiation time in the absence and presence of oxygen, and surface annealing temperature on both "clean" and HO-covered TiO2 nanoparticle films. Computationally, adsorption energies, molecular structures, and vibrational frequencies of possible adsorbates have been predicted by first-principles calculations based on the density functional theory (DFT) and the pseudopotential method. Three important reactions involving CH3 elimination, CH4 elimination, and CH3 migration from the adsorbed trimethylindium have been elucidated in detail. CH(3 migration is the only exothermic process with the lowest reaction barrier. On the basis of experimental and computational results, the two sharpest peaks at 2979 and 2925 cm(-1), detected in the dosage and UV irradiation experiments in the absence of oxygen, are attributable to the asymmetric and symmetric C-H vibrations of methyl groups in In(CH3)3(a) and its derivatives, (H3C)2In(a), H3CIn(a), and H3CO(a). In the UV irradiation experiment in the presence of oxygen, the methyl groups attached to the In atom were quickly oxidized to the methoxy with the C-H vibrations at 2925 and 2822 cm(-1) and to the carboxyl group with vibrations at 2888 cm(-1) (vs(CH)), 1577 cm(-1) (va(OCO)), 1380 cm(-1) (delta(CH)), and 1355 cm(-1) (vs(OCO)). Finally, from the computed energies with vibrational analysis, the adsorbed structure of the carboxyl group was confirmed to involve two oxygen atoms doubly adsorbed on two surface Ti atoms.

Water Clusters on Graphite: Methodology for Quantum Chemical A Priori Prediction of Reaction Rate Constants

The properties, interactions, and reactions of cyclic water clusters (H(2)O)(n=1-5) on model systems for a graphite surface have been studied using pure B3LYP, dispersion-augmented density functional tight binding (DFTB-D), and integrated ONIOM(B3LYP:DFTB-D) methods. Coronene C(24)H(12) as well as polycircumcoronenes C(96)H(24) and C(216)H(36) in monolayer, bilayer, and trilayer arrangements were used as model systems to simulate ABA bulk graphite. Structures, binding energies, and vibrational frequencies of water clusters on mono- and bilayer graphite models have been calculated, and structural changes and frequency shifts due to the water cluster-graphite interactions are discussed. ONIOM(B3LYP:DFTB-D) with coronene and water in the high level and C(96)H(24) in the low level mimics the effect of extended graphite pi-conjugation on the water-graphite interaction very reasonably and suggests that water clusters only weakly interact with graphite surfaces, as suggested by the fact that water is an excellent graphite lubricant. We use the ONIOM(B3LYP:DFTB-D) method to predict rate constants for model pathways of water dissociative adsorption on graphite. Quantum chemical molecular dynamics (QM/MD) simulations of water clusters and water addition products on the C(96)H(24) graphite model are presented using the DFTB-D method. A three-stage strategy is devised for a priori investigations of high temperature corrosion processes of graphite surfaces due to interaction with water molecules and fragments.

Computational Study on the Mechanism and Rate Constant for the C6H5 + C6H5NO Reaction

The reaction mechanism of C6H5 + C6H5NO involving four product channels on the doublet-state potential energy surface has been studied at the B3LYP/6-31+G(d, p) level of theory. The first reaction channel occurs by barrierless association forming (C6H5)2NO (biphenyl nitroxide), which can undergo isomerization and decomposition. The second channel takes place by substitution reaction producing C12H10 (biphenyl) and NO. The third and fourth channels involve direct hydrogen abstraction reactions producing C6H4NO + C6H6 and C6H5NOH + C6H4, respectively. Bimolecular rate constants of the above four product channels have been calculated in the temperature range 300-2000 K by the microcanonical Rice-Ramsperger-Kassel-Marcus theory and/or variational transition-state theory. The result shows the dominant reactions are channel 1 at lower temperatures (T < 800 K) and channel 3 at higher temperatures (T > 800 K). The total rate constant at 7 Torr He is predicted to be k(t) = 3.94 x 10(21) T(-3.09) exp(-699/T) for 300-500 K, 2.09 x 10(20) T(-3.56) exp(2315/T) for 500-1000 K, and 1.51 x 10(2) T(3.30) exp(-3043/T) for 1000-2000 K (in units of cm3 mol(-1) s(-1)), agreeing reasonably with the experimental data within their reported errors. The heats of formation of key products including biphenyl nitroxide, hydroxyl phenyl amino radical, and N-hydroxyl carbazole have been estimated.

Ab initio Studies of ClOx 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..

Ab initio kinetics of the reaction of HCO with NO: Abstraction versus association/elimination mechanism

The kinetics and mechanism for the reaction of HCO with NO occurring by both singlet and triplet electronic state potential-energy surfaces (PESs) have been studied at the modified Gaussian-2 level of theory based on the geometric parameters optimized by the Becke-3 Lee-Yang-Parr/6-311G(d,p) method. There are two major reaction channels on both singlet and triplet PESs studied: one is direct H abstraction producing CO+HNO and the other is association forming a stable HC(O)NO (nitrosoformaldehyde) molecule. The dominant reaction is predicted to be the direct H abstraction occurring primarily by the lowest-energy path via a loose hydrogen-bonding singlet molecular complex, ON...HCO, with a 2.9-kcal/mol binding energy and a small decomposition barrier (1.9 kcal/mol). The commonly assumed HC(O)NO intermediate, predicted to lie below the reactants by 27.7 kcal/mol, has a high HNO-elimination barrier (34.5 kcal/mol). Bimolecular rate constants for the formation of the singlet products and their branching ratios have been calculated in the temperature range of 200-3000 K. The rate constant for the disproportionation process producing HNO+CO, found to be affected strongly by multiple reflections above the well of the complex at low temperature, is predicted to be k(HNO)=3.08 x 10(-12) T(0.10) exp(242T) for 200-500 K, and 1.72 x 10(-16) T(1.47) exp(888T) for 500-3000 K in units of cm(3) molecule(-1) s(-1). The high- and low-pressure rate constants for the association process forming HC(O)NO can be represented by k(infinity)=4.42 x 10(-11) T(0.25) exp(-28T) cm(3) molecule(-1) s(-1) (200-3000 K) and k(0)=7.30x10(-16) T(-5.75) exp(-719T) (200-1000 K) and 1.82 x 10(2) T(-11.92) exp(1846T) (1000-3000 K) cm(6) molecule(-2) s(-1) for N(2)-buffer gas. The absolute values of total rate constant, predicted to be weakly dependent negatively on temperature but positively on pressure, are in close agreement with most experimental data within their reported errors.

Reaction mechanism and kinetics of the NCN+NO reaction: Comparison of theory and experiment

The rate constants for the NCN + NO reaction have been measured by laser photolysis/laser-induced fluorescence technique in the temperature range of 254-353 K in the presence of He (40-600 Torr) and N2 (30-528 Torr) buffer gases. The NCN radical was produced from the photodissociation of NCN3 at 193 nm and monitored with a dye laser at 329.01 nm. The reaction was found to be strongly positive-pressure dependent with negative-temperature dependence, as was reported previously. The experimental data could be reasonably accounted for by dual-channel Rice-Ramsperger-Kassel-Marcus calculations based on the predicted potential-energy surface using the modified Gaussian-2 method. The reaction is predicted to occur via weak intermediates, cis- and trans-NCNNO, in the 2A" state which crosses with the 2A' state containing more stable cis- and trans-NCNNO isomers. The high barriers for the fragmentation of these isomers and their trapping in the 2A' state by collisional stabilization give rise to the observed positive-pressure dependence and negative-temperature effect. The predicted energy barrier for the fragmentation of the cis-NCNNO (2A') to CN + N2O also allows us to quantitatively account for the rate constant previously measured for the reverse process CN + N2O --> NCN + NO.

Computational Study on the Kinetics and Mechanism for the Unimolecular Decomposition of C6H5NO2and the Related C6H5+ NO2and C6H5O + NO Reactions†

The kinetics and mechanisms for the unimolecular dissociation of nitrobenzene and related association reactions C(6)H(5) + NO(2) and C(6)H(5)O + NO have been studied computationally at the G2M(RCC, MP2) level of theory in conjunction with rate constant prediction with multichannel RRKM calculations. Formation of C(6)H(5) + NO(2) was found to be dominant above 850 K with its branching ratio > 0.78, whereas the formation of C(6)H(5)O + NO via the C(6)H(5)ONO intermediate was found to be competitive at lower temperatures, with its branching ratio increasing from 0.22 at 850 K to 0.97 at 500 K. The third energetically accessible channel producing C(6)H(4) + HONO was found to be uncompetitive throughout the temperature range investigated, 500-2000 K. The predicted rate constants for C(6)H(5)NO(2) --> C(6)H(5) + NO(2) and C(6)H(5)O + NO --> C(6)H(5)ONO under varying experimental conditions were found to be in good agreement with all existing experimental data. For C(6)H(5) + NO(2), the combination processes producing C(6)H(5)ONO and C(6)H(5)NO(2) are dominant at low temperature and high pressure, while the disproportionation process giving C(6)H(5)O + NO via C(6)H(5)ONO becomes competitive at low pressure and dominant at temperatures above 1000 K.