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.     

Computational study on the mechanisms and energetics of trimethylindium reactions with H2O and H2S

The reactions of trimethylindium (TMIn) with H2O and H2S are relevant to the chemical vapor deposition of indium oxide and indium sulfide thin films. 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)/[6-31G(d,p)+Lanl2dz]//B3LYP/[6-31G(d,p)+Lanl2dz] and CCSD(T)/[6-31G(d,p)+Lanl2dz] //MP2/[6-31G(d,p)+Lanl2dz] methods. The results of both methods are in good agreement for the optimized geometries and relative energies. When TMIn reacts with H2O and H2S, initial molecular complexes [(CH3)3In:OH2 (R1)] and [(CH3)3In:SH2 (R2)] are formed with 12.6 and 3.9 kcal/mol binding energies. Elimination of a CH4 molecule from each complex occurs with a similar energy barrier at TS1 (19.9 kcal/mol) and at TS3 (22.1 kcal/mol), respectively, giving stable intermediates (CH3)2InOH and (CH3)2InSH. The elimination of the second CH4 molecule from these intermediate products, however, has to overcome very high and much different barriers of 66.1 and 53.2 kcal/mol, respectively. In the case of DMIn with H2O and H2S reactions, formation of both InO and InS is exothermic by 3.1 and 30.8 kcal/mol respectively. On the basis of the predicted heats of formation of R1 and R2 at 0 K and -20.1 and 43.6 kcal/mol, the heats of formation of (CH3)2InOH, (CH3)2InSH, CH3InO, CH3InS, InO, and InS are estimated to be -20.6, 31.8, and 29.0 and 48.4, 35.5, and 58.5 kcal/mol, respectively. The values for InO and InS are in good agreement with available experimental data. A similar study on the reactions of (CH3)2In with H2O and H2S has been carried out; in these reactions CH3InOH and CH3InSH were found to be the key intermediate products.     

Potential energy surface and thermochemistry for the direct gas phase reaction of germane and water

Gas phase reaction between germane GeH₄ and water H₂O was investigated at CCSD(T)/[aug-cc-pVTZ-pp for Ge + Lanl2dz for H and O]//MP2/6-31G(d,p) level. Only the hydrogen elimination channels are monitored. Within the energy range of 100 kcal/mol, we located nine equilibrium and six transition states on the potential energy surface (PES) of the Ge–O–H systems. GeH₄ reacts with H₂O exothermically (by 2.37 kcal/mol) without a barrier to form a non-planar complex GeH₄·H₂O which isomerizes to GeH₃OH·H₂ and H₂GeOH₂·H₂ with a barrier of 44.34 kcal/mol and 53.75 kcal/mol respectively. The first step of hydrogen elimination gives two non-planar species, GeH₃OH and H₂GeOH₂ but germinol GeH₃OH is found to be more stable. Further thermal decomposition reactions of GeH₃OH involving hydrogen elimination have been studied extensively using the same method. The final hydrogen elimination step gives HGeOH which can exist in cis and trans forms. As the trans form is more stable, only the trans form is considered on the potential energy surface (PES) of the reaction. The important thermochemical parameters (∆_rE_tot + ZPE), ∆_rH and ∆_rG for the H₂ elimination pathways are predicted accurately.     

Computational methods to calculate accurate activation and reaction energies of 1,3-dipolar cycloadditions of 24 1,3-dipoles

Theoretical calculations were performed on the 1,3-dipolar cycloaddition reactions of 24 1,3-dipoles with ethylene and acetylene. The 24 1,3-dipoles are of the formula X≡Y(+)-Z(-) (where X is HC or N, Y is N, and Z is CH(2), NH, or O) or X═Y(+)-Z(-) (where X and Z are CH(2), NH, or O and Y is NH, O, or S). The high-accuracy G3B3 method was employed as the reference. CBS-QB3, CCSD(T)//B3LYP, SCS-MP2//B3LYP, B3LYP, M06-2X, and B97-D methods were benchmarked to assess their accuracies and to determine an accurate method that is practical for large systems. Several basis sets were also evaluated. Compared to the G3B3 method, CBS-QB3 and CCSD(T)/maug-cc-pV(T+d)Z//B3LYP methods give similar results for both activation and reaction enthalpies (mean average deviation, MAD, < 1.5 kcal/mol). SCS-MP2//B3LYP and M06-2X give small errors for the activation enthalpies (MAD < 1.5 kcal/mol), while B3LYP has MAD = 2.3 kcal/mol. SCS-MP2//B3LYP and B3LYP give the reasonable reaction enthalpies (MAD < 5.0 kcal/mol). The B3LYP functional also gives good results for most 1,3-dipoles (MAD = 1.9 kcal/mol for 17 common 1,3-dipoles), but the activation and reaction enthalpies for ozone and sulfur dioxide are difficult to calculate by any of the density functional methods.     

Adiabatic and nonadiabatic dynamics in the CH3(CD3)+HCl reaction

The scattering dynamics leading to the formation of Cl (2P(3/2)) and Cl* (2P(1/2)) products of the CH(3)+HCl reaction (at a mean collision energy =22.3 kcal mol(-1)) and the Cl (2P(3/2)) products of the CD(3)+HCl reaction (at =19.4 kcal mol(-1)) have been investigated by using photodissociation of CH(3)I and CD(3)I as sources of translationally hot methyl radicals and velocity map imaging of the Cl atom products. Image analysis with a Legendre moment fitting procedure demonstrates that, in all three reactions, the Cl/Cl* products are mostly forward scattered with respect to the HCl in the center-of-mass (c.m.) frame but with a backward scattered component. The distributions of the fraction of the available energy released as translation peak at f(t)=0.31-0.33 for all the reactions, with average values that lie in the range =0.42-0.47. The detailed analysis indicates the importance of collision energy in facilitating the nonadiabatic transitions that lead to Cl* production. The similarities between the c.m.-frame scattering and kinetic energy release distributions for Cl and Cl* channels suggest that the nonadiabatic transitions to a low-lying excited potential energy surface (PES) correlating to Cl* products occur after passage through the transition state region on the ground-state PES. Branching fractions for Cl* are determined to be 0.14+/-0.02 for the CH(3)+HCl reaction and 0.20+/-0.03 for the CD(3)+HCl reaction. The difference cannot be accounted for by changes in collision energy, mass effects, or vibrational excitation of the photolytically generated methyl radical reagents and instead suggests that the low-frequency bending modes of the CD(3)H or CH(4) coproduct are important mediators of the nonadiabatic couplings occurring in this reaction system.     

Gaseous reaction mechanism of C2F radical with water

The kinetic properties of the carbon-fluorine radicals are little understood except those of CFn (n =1-3). In this article, a detailed mechanistic study was reported on the gas-phase reaction between the simplest pi-bonded C2F radical and water as the first attempt to understand the chemical reactivity of the C2F radical. Various reaction channels are considered. The most kinetically competitive channel is the quasi-direct hydrogen-abstraction route forming P5 HCCF + OH. At the CCSD(T)/6-311+G(2d,2p)//B3LYP/6-311G(d,p)+ZPVE, CCSD(T)/6-311+G(3df,2p)//QCISD/6-311G(d,p)+ZPVE and Gaussian-3//B3LYP/6-31G(d) levels, the overall H-abstraction barriers (4.5, 4.7, and 4.2 kcal/mol) for the C2F + H2O reaction are comparable to the corresponding values (5.5, 3.7, and 5.7 kcal/mol) for the analogous C2H + H2O reaction. This suggests that C2F is a reactive radical like the extensively studied C2H, in contrast to the situation of the CF and CF2 radicals that have much lower reactivity than the corresponding hydrocarbon species. Thus, the C2F radical is expected to play an important role in the combustion processes of the carbon-fluorine chemistry. Furthermore, addition of a second H2O can catalyze the reaction with the H-abstraction barrier significantly reduced to a marginally zero value (0.5 kcal/mol). This is also indicative of the potential relevance of the title reactions in the low-temperature atmospheric chemistry.     

Dynamics of the gas-phase reactions of chloride ion with fluoromethane: high excess translational activation energy for an endothermic S(N)2 reaction.

Guided ion beam tandem mass spectrometry techniques are used to examine the competing product channels in the reaction of Cl(-) with CH(3)F in the center-of-mass collision energy range 0.05-27 eV. Four anionic reaction products are detected: F(-), CH(2)Cl(-), FCl(-), and CHCl(-). The endothermic S(N)2 reaction Cl(-) + CH(3)F --> CH(3)Cl + F(-) has an energy threshold of E(0) = 181 +/- 14 kJ/mol, exhibiting a 52 +/- 16 kJ/mol effective barrier in excess of the reaction endothermicity. The potential energy of the S(N)2 transition state is well below the energy of the products. Dynamical impedances to the activation of the S(N)2 reaction are discussed, including angular momentum constraints, orientational effects, and the inefficiency of translational energy in promoting the reaction. The fluorine abstraction reaction to form CH(3) + FCl(-) exhibits a 146 +/- 33 kJ/mol effective barrier above the reaction endothermicity. Direct proton transfer to form HCl is highly inefficient, but HF elimination is observed above 268 +/- 95 kJ/mol. Potential energy surfaces for the reactions are calculated using the CCSD(T)/aug-cc-pVDZ and HF/6-31+G(d) methods and used to interpret the dynamics.     

Reliable predictions of the thermochemistry of boron-nitrogen hydrogen storage compounds: BxNxHy, x = 2, 3

Thermochemical data calculated using ab initio molecular orbital theory are reported for 16 BxNxHy compounds with x = 2, 3 and y > or = 2x. Accurate gas-phase heats of formation were obtained using coupled cluster with single and double excitations and perturbative triples (CCSD(T)) valence electron calculations extrapolated to the complete basis set (CBS) limit with additional corrections including core/valence, scalar relativistic, and spin-orbit corrections to predict the atomization energies and scaled harmonic frequencies to correct for zero point and thermal energies and estimate entropies. Computationally cheaper calculations were also performed using the G3MP2 and G3B3 variants of the Gaussian 03 method, as well as density functional theory (DFT) using the B3LYP functional. The G3MP2 heats of formation are too positive by up to approximately 6 kcal/mol as compared with CCSD(T)/CBS values. The more expensive G3B3 method predicts heats of formation that are too negative as compared with the CCSD(T)/CBS values by up to 3-4 kcal/mol. DFT using the B3LYP functional and 6-311+G** basis set predict isodesmic reaction energies to within a few kcal/mol compared with the CCSD(T)/CBS method so isodesmic reactions involving BN compounds and the analogous hydrocarbons can be used to estimate heats of formation. Heats of formation of c-B3N3H12 and c-B3N3H6 are -95.5 and -115.5 kcal/mol at 298 K, respectively, using our best calculated CCSD(T)/CBS approach. The experimental value for c-B3N3H6 appears to be approximately 7 kcal/mol too negative. Enthalpies, entropies, and free energies are calculated for many dehydrocoupling and dehydrogenation reactions that convert BNH6 to alicyclic and cyclic oligomers and H2(g). Generally, the reactions are highly exothermic and exergonic as well because of the release of 1 or more equivalents of H2(g). For c-B3N3H12 and c-B3N3H6, available experimental data for sublimation and vaporization lead to estimates of their condensed phase 298 K heats of formation: DeltaHf degrees [c-B3N3H12(s)] = -124 kcal/mol and DeltaHf degrees [c-B3N3H6(l)] = -123 kcal/mol. The reaction thermochemistries for the dehydrocoupling of BNH6(s) to c-B3N3H12(s) and the dehydrogenation of c-B3N3H12(s) to c-B3N3H6(l) are much less exothermic compared with the gas-phase reactions due to intermolecular forces which decrease in the order BNH6 > cyclo-B3N3H12 > cyclo-B3N3H6. The condensed phase reaction free energies are less negative compared with the gas-phase reactions but are still too favorable for BNH6 to be regenerated from either c-B3N3H12 or c-B3N3H6 by just an overpressure of H2.     

Energetics and mechanism of the decomposition of trifluoromethanol

The thermal instability of alpha-fluoroalcohols is generally attributed to a unimolecular 1,2-elimination of HF, but the barrier to intramolecular HF elimination from CF3OH is predicted to be 45.1 +/- 2 kcal/mol. The thermochemical parameters of trifluoromethanol were calculated using coupled-cluster theory (CCSD(T)) extrapolated to the complete basis set limit. High barriers of 42.9, 43.1, and 45.0 kcal/mol were predicted for the unimolecular decompositions of CH2FOH, CHF2OH, and CF3OH, respectively. These barriers are lowered substantially if cyclic H-bonded dimers of CF3OH with complexation energies of approximately 5 kcal/mol are involved. A six-membered ring dimer has an energy barrier of 28.7 kcal/mol and an eight-membered dimer has an energy barrier of 32.9 kcal/mol. Complexes of CF3OH with HF lead to strong H-bonded dimers, trimers and tetramers with complexation energies of approximately 6, 11, and 16 kcal/mol, respectively. The dimer, CH3OH:HF, and the trimers, CF3OH:2HF and (CH3OH)2:HF, have decomposition energy barriers of 26.7, 20.3, and 22.8 kcal/mol, respectively. The tetramer (CH3OH:HF)2 gives rise to elimination of two HF molecules with a barrier of 32.5 kcal/mol. Either CF3OH or HF can act as catalysts for HF-elimination via an H-transfer relay. Because HF is one of the decomposition products, the decomposition reactions become autocatalytic. If the energies due to complexation for the CF3OH-HF adducts are not dissipated, the effective barriers to HF elimination are lowered from approximately 20 to approximately 9 kcal/mol, which reconciles the computational results with the experimentally observed stabilities.     

Carboranes and baskets from reaction of B4H10 with allene

The reaction of the boron hydride B4H10 with allene was studied at the CCSD(T)/6-311+G(d)//MP2/6-31G(d) level. The mechanism is surprisingly complex with 44 transition states and several branching points located. The four carboranes and one basket that have been observed experimentally are all connected by pathways that have very similar free energies of activation. In addition, two new structures, a basket (2,4-(CH2CH2CH2)B4H8, 5a) and a "classical" structure (1,4-(Me2C)bisdiborane, 7), which might be obtained from the B4H10 + C3H4 reaction under the right conditions (hot/cold, quenched, etc.) have been identified. The first branch point in the reaction is the competition between H2 elimination from B4H10 (DeltaG(298 K) = 32.2 kcal/mol) and the hydroboration of allene by B4H10 (DeltaG(298 K) = 31.3 kcal/mol). The next branch point in the hydroboration mechanism controls the formation of 2,4-(MeCHCH2)B4H8 (1) (DeltaG(298 K) = 31.5 kcal/mol) and arachno-1,2/arachno-1,3-Me2-1-CB4H7 (8 and 8a) (DeltaG(298 K) = 34.3 kcal/mol). Another branch point in the H2-elimination mechanism controls the formation of 1-Me-2,5-micro-CH2-1-CB4H7 (29) (DeltaG(298 K) = 0.1 kcal/mol) and 2,5-micro-CHMe-1-CB4H7 (25/26) (DeltaG(298 K) = 7.3 kcal/mol). Formation of 2-Me-2,3-C2B4H7, a carborane observed in the reaction of methylacetylene with B4H10, is calculated to be blocked by a high barrier for H2 elimination. All free energies are relative to B4H10 + allene. An interesting reaction step discovered is the "reverse hydroboration step" in which a hydrogen atom is transferred from carbon back to boron, which allows a CH hydrogen to shuttle between the terminal and central carbon of allene.     

Kinetics of the multichannel reaction of methanethiyl radical (CH3S*) with 3O2

The CH3S* + O2 reaction system is considered an important process in atmospheric chemistry and in combustion as a pathway for the exothermic conversion of methane-thiyl radical, CH3S*. Several density functional and ab initio computational methods are used in this study to determine thermochemical parameters, reaction paths, and kinetic barriers in the CH3S* + O2 reaction system. The data are also used to evaluate feasibility of the DFT methods for higher molecular weight oxy-sulfur hydrocarbons, where sulfur presents added complexity from its many valence states. The methods include: B3LYP/6-311++G(d,p), B3LYP/6-311++G(3df,2p), CCSD(T)/6-311G(d,p)//MP2/6-31G(d,p), B3P86/6-311G(2d,2p)//B3P86/6-31G(d), B3PW91/6-311++G(3df,2p), G3MP2, and CBS-QB3. The well depth for the CH3S* + 3O2 reaction to the syn-CH3SOO* adduct is found to be 9.7 kcal/mol. Low barrier exit channels from the syn-CH3SOO* adduct include: CH2S + HO2, (TS6, E(a) is 12.5 kcal/mol), CH3 + SO2 via CH3SO2 (TS2', E(a) is 17.8) and CH3SO + O (TS17, E(a) is 24.7) where the activation energy is relative to the syn-CH3SOO* stabilized adduct. The transition state (TS5) for formation of the CH3SOO adduct from CH3S* + O2 and the reverse dissociation of CH3SOO to CH3S* + O2 is relatively tight compared to typical association and simple bond dissociation reactions; this is a result of the very weak interaction. Reverse reaction is the dominant dissociation path due to enthalpy and entropy considerations. The rate constants from the chemical activation reaction and from the stabilized adduct to these products are estimated as functions of temperature and pressure. Our forward rate constant and CH3S loss profile are in agreement with the experiments under similar conditions. Of the methods above, the G3MP2 and CBS-QB3 composite methods are recommended for thermochemical determinations on these carbon-sulfur-oxygen systems, when they are feasible.     

Establishment of the C(2)H(5)+O(2) reaction mechanism: a combustion archetype

The celebrated C(2)H(5)+O(2) reaction is an archetype for hydrocarbon combustion, and the critical step in the process is the concerted elimination of HO(2) from the ethylperoxy intermediate (C(2)H(5)O(2)). Master equation kinetic models fitted to measured reaction rates place the concerted elimination barrier 3.0 kcal mol(-1) below the C(2)H(5)+O(2) reactants, whereas the best previous electronic structure computations yield a barrier more than 2.0 kcal mol(-1) higher. We resolve this discrepancy here by means of the most rigorous computations to date, using focal point methods to converge on the ab initio limit. Explicit computations were executed with basis sets as large as cc-pV5Z and correlation treatments as extensive as coupled cluster through full triples with a perturbative inclusion of quadruple excitations [CCSDT(Q)]. The final predicted barrier is -3.0 kcal mol(-1), bringing the concerted elimination mechanism into precise agreement with experiment. This work demonstrates that higher correlation treatments such as CCSDT(Q) are not only feasible on systems of chemical interest but are necessary to supply accuracy beyond 0.5 kcal mol(-1), which is not obtained with the "gold standard" CCSD(T) method. Finally, we compute the enthalpy of formation of C(2)H(5)O(2) to be Delta(f)H degrees (298 K)=-5.3+/-0.5 kcal mol(-1) and Delta(f)H degrees (0 K)=-1.5+/-0.5 kcal mol(-1).     

Exploring the dynamics of reaction N((2)D)+C2H4 with crossed molecular-beam experiments and quantum-chemical calculations

We conducted the title reaction using a crossed molecular-beam apparatus, quantum-chemical calculations, and RRKM calculations. Synchrotron radiation from an undulator served to ionize selectively reaction products by advantage of negligibly small dissociative ionization. We observed two products with gross formula C(2)H(3)N and C(2)H(2)N associated with loss of one and two hydrogen atoms, respectively. Measurements of kinetic-energy distributions, angular distributions, low-resolution photoionization spectra, and branching ratios of the two products were carried out. Furthermore, we evaluated total branching ratios of various exit channels using RRKM calculations based on the potential-energy surface of reaction N((2)D)+C(2)H(4) established with the method CCSD(T)/6-311+G(3df,2p)//B3LYP/6-311G(d,p)+ZPE[B3LYP/6-311G(d,p)]. The combination of experimental and computational results allows us to reveal the reaction dynamics. The N((2)D) atom adds to the C=C π-bond of ethene (C(2)H(4)) to form a cyclic complex c-CH(2)(N)CH(2) that directly ejects a hydrogen atom or rearranges to other intermediates followed by elimination of a hydrogen atom to produce C(2)H(3)N; c-CH(2)(N)CH+H is the dominant product channel. Subsequently, most C(2)H(3)N radicals, notably c-CH(2)(N)CH, further decompose to CH(2)CN+H. This work provides results and explanations different from the previous work of Balucani et al. [J. Phys. Chem. A, 2000, 104, 5655], indicating that selective photoionization with synchrotron radiation as an ionization source is a good choice in chemical dynamics research.     

HNNC radical and its role in the CH+N2 reaction

A previously unreported channel in the spin-allowed reaction path for the CH+N2 reaction that involves the HNNC radical is presented. The structures and energetics of the HNNC radical and its isomers HCNN and HNCN and the relevant intermediates and transition states that are involved in the proposed mechanism are obtained at the coupled cluster singles and doubles level of theory with noniterative triples correction (CCSD(T)) using a converging series of basis sets aug-cc-pVDZ, aug-cc-pVTZ, and aug-cc-pVQZ. The aug-cc-pVQZ basis is used for all the final single point energy calculations using the CCSD(T)/aug-cc-pVTZ optimized geometries. We find the HNNC radical to have a heat of formation of DeltafH0 (HNNC)=116.5 kcal mol(-1). An assessment of the quality of computed data of the radical species HNCN and HCNN is presented by comparison with the available experimental data. We find that HNNC can convert to HNCN, the highest barrier in this path being 14.5 kcal mol(-1) above the energy of the CH+N2 reactants. Thus, HNNC can play a role in the high-temperature spin-allowed mechanism for the reaction of CH+N2 proposed by Moskaleva, Xia, and Lin (Chem. Phys. Lett. 2000, 331, 269).     

Theoretical study of the reactions of the 1Σ+ ground state of MS+ (M = Sc,Y, and La) with oxygen‐transfer reagent MS+ + CO → ScO+ + CS in the gas phase

The reaction mechanisms of the ¹Σ⁺ ground state of MS⁺ (M = Sc, Y, and La) with oxygen‐transfer reagent MS⁺ + CO → MO⁺ + CS in the gas phase has been proposed and investigated by ab initio methods with the 6‐31G* basis set for nonmetal atoms and the effective core potentials of Lanl2dz for the metal atoms. A carbon migration from oxygen atom to sulfur atom via a four‐center transition state is involved on the reaction potential surface. The activation energies of the reactions are 34.0, 24.1, and 36.7 kcal/mol relative to their corresponding reactants and the reaction heats are 15.7, 18.6, and 18.0 kcal/mol (respectively, for M = Sc, Y, and La) at the MP4 (SDTQ)/6‐31G*//MP2/6‐31G* level plus zero‐point energy, which indicates that the cationic yttrium sulfide is more favorable for this type of reaction. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

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.     

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

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

Reactions of hydrazoic acid on TiO2 nanoparticles: an experimental and computational study

This article reports the results of a computational and experimental study on the reaction of hydrazoic acid, HN3, adsorbed on 15-20 nm TiO2 particle films. Experimentally, FTIR spectra of HN3(a) have been measured by varying HN3 dosage, UV irradiation time and surface annealing temperature. Three sharp peaks, related to v(a)(NNN) of HN3(a) and N3(a) with different configurations in the 2000-2200 cm(-1) region, and a broad band absorption, related to associated and isolated HN(a) and HO(a) adsorptions in the 3000-3800 cm(-1) region, have been detected. Computationally, molecular structures, vibrational frequencies and adsorption energies of possible adsorbates including HN3 and its derivatives, N3, N2, NH, and H, have been predicted by first-principles calculations based on the density functional theory (DFT) and the pseudopotential method. On the basis of the experimental and computational results, the peak appeared at 2075 cm(-1), which increases at a faster rate with HN3 exposure time, is attributed to a stable adsorbate, N3-Ti(a), with the predicted adsorption energy, E(ads) = 13 kcal/mol. The peak at 2118 cm(-1), which survives at the highest surface temperature in the heating experiment, is attributable to the most stable adsorbate, Ti-N2N(H)-O(a) with E(ads) = 36 kcal/mol. The peak at 2170 cm(-1), which vanishes most readily in all of the aforementioned experiments, is related to less stable molecular adsorbates, end-on HN3-Ti(a) with E(ads) = 5 kcal/mol and side-on HN(N2)-Ti(a) with E(ads) = 8 kcal/mol. A potential energy diagram for the formation of various absorbates with their transition states has been established for the HN3/TiO2 system. On the basis of the predicted desorption energies, the four most stable products of the HN3 reaction on TiO2 are H-O(a), 118 kcal/mol; HN-O(a), 85 kcal/mol; Ti-N2N(H)-O(a), 36 kcal/mol; and N3-O(a), 19 kcal/mol.