Quantum chemical prediction of pathways and rate constants for reaction of cyanomethylene radical with NO

High-level ab initio calculations have been performed to study the mechanism and kinetics of the reaction of the cyanomethylene radical (HCCN) with the NO. The species involved have been optimized at the B3LYP/6-311++G(3df,2p) level, and their corresponding single-point energies are improved by the CCSD(T)/aug-cc-PVQZ//B3LYP/6-311++G(3df,2p) approach. From the calculated potential energy surface, we have predicted the favorable pathways for the formation of several isomers of a HCCN-NO complex. Barrierless formation of HCN + NCO (P1) is also possible. Formation of HCNO + CN (P3) is endoergic but may become significant at high temperatures. To rationalize the scenario of our calculated results, we also employ the Fukui functions and hard-and-soft acid-and-base (HSAB) theory to seek possible clues. The predicted total rate coefficient, k(total), at He pressure 760 Torr can be represented with the equation k(total) = 1.40 × 10(-7) T(-2.01) exp(3.15 kcal mol(-1)/RT) at T = 298-3000 K in units of cm(3) molecule(-1) s(-1). The predicted total rate coefficients at some available conditions (He pressures of 6, 18, and 30 Torr in the temperature of 298 K) are in reasonable agreement with experimental observation. In addition, the rate constants for key individual product channels are provided in different temperature and pressure conditions.     

The reaction rate for dissociative adsorption of N2 on stepped Ru(0001): six-dimensional quantum calculations

Quantum-mechanical calculations of the reaction rate for dissociative adsorption of N2 on stepped Ru(0001) are presented. Converged six-dimensional quantum calculations for this heavy-atom reaction have been performed using the multiconfiguration time-dependent Hartree method. A potential-energy surface for the transition-state region is constructed from density-functional theory calculations using Shepard interpolation. The quantum results are in very good agreement with the results of the harmonic transition-state theory. In contrast to the findings of previous model calculations on similar systems, the tunneling effect is found to be small.     

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.     

Hydrogen adsorption on graphite (0001) surface: a combined spectroscopy-density-functional-theory study

The adsorption of H/D atoms on the graphite (0001) surface is investigated by means of both high-resolution electron-energy loss spectroscopy (HREELS) and periodic first-principle density-functional theory. The two methods converge towards two modes of adsorption: adsorption in clusters of about four hydrogen atoms and adsorption in pairs of atoms on contiguous carbon sites. The desorption energies estimated from the calculated dissociation energies range from 8 to 185 kJ mol(-1) leading to an estimated surface coverage at saturations of 30-44 at. %. These results are compared with previous thermal desorption spectroscopy results. New HREEL signal assignments are proposed based on quantum calculations.     

Water adsorption on the hydroxylated H-(1x1) O-ZnO(0001) surface

The adsorption of water multilayers on a well defined single crystal, hydroxyl-terminated ZnO-surface, H(1x1)-O-ZnO(0001) surface has been investigated using infrared (IR) spectroscopy, helium atom scattering (HAS) and X-ray photoelectron spectroscopy (XPS). The results reveal the formation of well ordered mono-, bi- and multilayers of D2O and H2O on this substrate. On the bare hydroxyl-covered H(1x1) surface the OH-stretch vibration could be clearly identified in the IR-spectra. The water adsorption and desorption kinetics on this hydroxylated surface were studied by monitoring the reflectivity of the surface for helium atoms. The analysis of the data yielded activation energies for desorption of H2O from the H(1x1) O-ZnO surface of 55.2 kJ mol-1. The results reveal the formation of ordered mono- and bilayers. Further exposure to water at 113 K results in the formation of amorphous 3-D islands.     

Quantum chemical and conventional TST calculations of rate constants for the OH + alkane reaction

Reactions of OH with methane, ethane, propane, i-butane, and n-butane have been modeled using ab initio (MP2) and hybrid DFT (BHandHLYP) methods, and the 6-311G(d,p) basis set. Furthermore, single-point calculations at the CCSD(T) level were carried out at the optimized geometries. The rate constants have been calculated using the conventional transition-state theory (CTST). Arrhenius equations are proposed in the temperature range of 250–650 K. Hindered Internal Rotation partition functions calculations were explicitly carried out and included in the total partition functions. These corrections showed to be relevant in the determination of the pre-exponential parameters, although not so important as in the NO₃ + alkane reactions [G. Bravo-Pérez, J.R. Alvarez-Idaboy, A. Cruz-Torres, M.E. Ruíz, J. Phys. Chem. A 106 (2002) 4645]. The explicit participation of the tunnel effect has been taken into account. The calculated rate coefficients provide a very good agreement with the experimental data. The best agreement for the overall alkane + OH reactions seemed to occur when the BHandHLYP geometries and partition functions are used. For propane and i-butane, in addition to the respective secondary and tertiary H-abstraction channels, the primary one has been considered. These pathways are confirmed to be significant in spite of the large differences in activation energies between primary and secondary or primary and tertiary channels, respectively of propane and i-butane reactions and should not be disregarded.     

XUV free-electron laser desorption of NO from graphite (0001)

We report results of femtosecond laser induced desorption of NO from highly oriented pyrolytic graphite using XUV photon energies of hν = 38 eV and hν = 57 eV. Femtosecond pulses with a pulse energy of up to 40 μJ and about 30 fs duration generated at FLASH are applied. The desorbed molecules are detected with rovibrational state selectivity by (1 + 1) REMPI in the A(2)Σ(+) ← X(2)Π γ-bands around λ = 226 nm. A nonlinear desorption yield of neutral NO is observed with an exponent of m = 1.4 ± 0.2. At a fluence of about 4 mJ/cm(2) a desorption cross section of σ(1) = (1.1 ± 0.4) × 10(-17) cm(2) is observed, accompanied with a lower one of σ(2) = (2.6 ± 0.3) × 10(-19) cm(2) observable at higher total fluence. A nonthermal rovibrational population distribution is observed with an average rotational energy of = 38.6 meV (311 cm(-1)), a vibrational energy of = 136 meV (1097 cm(-1)) and an electronic energy of = 3.9 meV (31 cm(-1)).     

A DFT study of reaction pathways of NH(3) decomposition on InN (0001) surface

Reaction pathways for complete decomposition of ammonia on the InN (0001) surface are investigated using first principles calculations. We show that while the initial NH(3) decomposition on this surface can proceed by H dissociation, its further decomposition is most favorable by H transfer. The calculated low diffusion barriers for the decomposed species on the surface imply that the metal-organic chemical vapor deposition growth of InN is a reaction-limited process rather than diffusion-limited at low adsorbate coverage.     

Effect of carbon surface modification with dimethylamine on reactive adsorption of NO(x)

The wood-based activated carbon, either as received or oxidized with nitric acid, was exposed to dimethylamine vapors. This modification was expected to introduce nitrogen groups. Then, the modified samples were used as adsorbents of NO(2) under dynamic conditions. Both NO(2) breakthrough curves and the NO concentration curves were recorded. The samples before and after exposure to NO(2) were characterized using adsorption of nitrogen, elemental analysis, potentiometric titration, FTIR, and thermal analysis. Modifications with amines resulted in an increase in NO(2) adsorption and in a decrease in NO emission. The effects were more visible when oxidation was used as a pretreatment of the carbon surface. This process increased the incorporation of nitrogen to the carbon matrix via acid-based reactions resulting in the formation of amides and amine carboxylic salts. Besides this, dimethylamine was strongly adsorbed on the carbon surface via hydrogen bonding with oxygen-containing groups. When the samples were exposed to nitrogen dioxide, there was an indication that nitramine and nitrosoamine were formed in the reactions of NO(2) with either amides or amines. In the reactions of amines with NO, nitrosoamines are the likely products. As a next step, the surface of the carbon matrix is reoxidized by NO(2), which is accompanied by the release of NO.     

Quantum chemical study of adsorption and dissociation of H2S on the gallium-rich GaAs (001)-4 x 2 surface

H(2)S adsorption and dissociation on the gallium-rich GaAs(001)-4 x 2 surface is investigated using hybrid density functional theory. Starting from chemisorbed H(2)S on the GaAs(001)-4 x 2 surface, two possible reaction routes have been proposed. We find that H(2)S adsorbs molecularly onto GaAs(001)-4 x 2 via the formation of a dative bond, and this process is exothermic with adsorption energy of 6.6 kcal/mol. For the first reaction route, one of the H atoms from the chemisorbed H(2)S is transferred to a second-layer As atom and the dissociated SH is inserted into the Ga-As bond with an activation barrier of 8.2 kcal/mol, which is found to be 29.3 kcal/mol more stable than the reactants. For the second case, the dissociated species may insert themselves into the Ga-Ga dimer resulting in the Ga-H-Ga and Ga-HS-Ga bridge-bonded states, which are found to be 29.8 and 22.2 kcal/mol more stable than the reactants, respectively. However, the calculations also show that the activation barrier (16.1 kcal/mol) for chemisorbed H(2)S dissociation through the second route is higher than the transfer of one H atom into a second-layer As atom. As a result, we conclude that sulfur insertion into the Ga-As bond is more kinetically favorable.     

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

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

Theoretical study of hydrogen dissociative adsorption on the Cu(110) surface

We have calculated the six-dimensional (6D) potential energy surface for H2 in front of a frozen Cu(110) surface using density functional theory for 22 H2-surface configurations and the corrugation reducing procedure to interpolate between them. We carry out classical trajectory calculations on the dissociative adsorption process and find excellent agreement with measurements. We find that it is of prominent importance to account for the rovibrational state distribution in the incident H2 beam. A straightforward analysis leads to the conclusion that the motion along the surface does not play an appreciable role in the dynamics whereas the dynamical role of molecular rotation is crucial. The latter fact precludes any interpretation of dissociation in terms of a static concept such as "barrier distributions."     

Harmonic force field between the (0001) surface of graphite and adsorbed methane

Empirical potentials were used to approximate the interaction between the (0001) surface of graphite and a single methane molecule. The frequency of the vibration perpendicular to the surface was obtained by assuming that the force field under consideration is harmonic, and also by solving the one-dimensional Schrödinger equation of the vibration. It was found that the two approaches yielded nearly identical excitation energies for the fundamental transition.     

The dissociative adsorption of silane and disilane on Si(100)-(2 x 1)

We investigate the dissociative adsorption of silane and disilane on Si(100)-(2 x 1) using pseudopotential planewave density functional theory calculations. These are important steps in the growth of silicon films. Although silane has been studied computationally in some detail previously, we find physisorbed precursor states for the intradimer and interdimer channels. The silane energetics calculated here are in good agreement with experimental data and previous theoretical estimates and provide us with a useful reference point for our disilane calculations. Disilane has not been studied as intensively as silane. We investigate both silicon-silicon bond cleavage and silicon-hydrogen bond cleavage mechanisms, and for each we investigate intradimer, interdimer, and inter-row channels. As in the case of silane, we also find precursor states in the adsorption path in agreement with molecular beam experiments. The qualitative picture that emerges is that adsorption takes place through a weakly bound precursor state with a transition state to chemisorption that is low lying in energy relative to the gas phase. This is in good agreement with experimental data. However, the calculated energetics are only in fair agreement with experiments, with our transition state to chemisorption being about 0.02 eV above the gas phase while experimentally it is estimated to be approximately 0.28 eV below the gas phase. This suggests that accurate theoretical characterization of these weakly bound precursor states and the adsorption barriers requires further computational work.     

Mechanisms of H2 dissociative adsorption on the Pt(211) stepped surface

We utilize classical trajectory calculations to study the reaction dynamics of the dissociative adsorption of H2 on the stepped Pt(211) surface. The potential-energy surface has been obtained through an accurate interpolation of density-functional theory data at the generalized gradient approximation level, using the corrugation reduction procedure. New techniques for visualizing the collective dynamics of trajectories are introduced to elucidate the reaction mechanisms involved. Reaction exhibits a nonmonotonic dependence on collision energy, first decreasing with energy, and then increasing. A strong component of direct nonactivated reaction exists at the top edge of the step over the entire range of energies. The inverse relationship between reaction and collision energy at low energies is attributed to trapping in weak chemisorption wells. These wells also influence the direct reaction at the step, leading to a strong asymmetric dependence on incidence angle. Reaction on the terrace is activated, and only contributes significantly at high energies. Agreement with experiments on Pt(533) [A. T. Gee, B. E. Hayden, C. Mormiche, and T. S. Nunney, J. Chem. Phys. 112, 7660 (2000); Surf. Sci. 512, 165 (2002)] is good, and we are able to suggest new interpretations of the experimental data.     

Evidence of stable high-temperature D(x)-CO intermediates on the Ru(0001) surface

We demonstrate the formation of complexes involving attractive interactions between D and CO on Ru(0001) that are stable at significantly higher temperatures than have previously been reported for such intermediate species on this surface. These complexes are evident by the appearance of new desorption features upon heating of the sample. They decompose in stages as the sample temperature is increased, with the most stable component desorbing at >500 K. The D:CO ratio remaining on the surface during the final stages of desorption tends towards 1:1. The new features are populated during normally incident molecular beam dosing of D(2) on to CO pre-covered Ru(0001) surfaces (180 K) when the CO coverage exceeds 50% of the saturation value. The amount of complex formed decreases somewhat with increasing CO pre-coverage. It is almost absent in the case of dosing on to the fully saturated surface. The results are interpreted in terms of both local and long-range rearrangements of the overlayer that give rise to the observed CO coverage dependence and limit the amount of complex that can be formed.     

First-principle study of adsorption of hydrogen on Ti-doped Mg(0001) surface

Ab initio density functional theory (DFT) calculations are performed to study the adsorption of H2 molecules on a Ti-doped Mg(0001) surface. We find that two hydrogen molecules are able to dissociate on top of the Ti atom with very small activation barriers (0.103 and 0.145 eV for the first and second H2 molecules, respectively). Additionally, a molecular adsorption state of H2 above the Ti atom is observed for the first time and is attributed to the polarization of the H2 molecule by the Ti cation. Our results parallel recent findings for H2 adsorption on Ti-doped carbon nanotubes or fullerenes. They provide new insight into the preliminary stages of hydrogen adsorption onto Ti-incorporated Mg surfaces.     

Quantum chemical simulations of water adsorption on a diamond (100) surface with vacancy defects

Results from quantum chemical calculations of the structural, electronic, and energy characteristics of the chemisorption of water on a diamond C(100)-(2 × 1) surface with a vacancy defect are presented. The metastable state of the surface with an adsorbed H₂O molecule and possible configurations of the surface with adsorbed -H and -OH water dissociation fragments are described. It is shown that the presence of a vacancy on the surface decreases the activation energy of the dissociative adsorption of a water molecule.     

Reactivity of V2O3(0001) surfaces: molecular vs dissociative adsorption of water

The adsorption of water on V2O3(0001) surfaces has been investigated by thermal desorption spectroscopy, high-resolution electron energy loss spectroscopy, and X-ray photoelectron spectroscopy with use of synchrotron radiation. The V2O3(0001) surfaces have been generated in epitaxial thin film form on a Rh(111) substrate with three different surface terminations according to the particular preparation conditions. The stable surface in thermodynamic equilibrium with the bulk is formed by a vanadyl (VO) (1x1) surface layer, but an oxygen-rich (radical3xradical3)R30 degrees reconstruction can be prepared under a higher chemical potential of oxygen (microO), whereas a V-terminated surface consisting of a vanadium surface layer requires a low microO, which can be achieved experimentally by the deposition of V atoms onto the (1x1) VO surface. The latter two surfaces have been used to model, in a controlled way, oxygen and vanadium containing defect centres on V2O3. On the (1x1) V=O and (radical3xradical3)R30 degrees surfaces, which expose only oxygen surface sites, the experimental results indicate consistently that the molecular adsorption of water provides the predominant adsorption channel. In contrast, on the V-terminated (1/radical3x1/radical3)R30 degrees surface the dissociation of water and the formation of surface hydroxyl species at 100 K is readily observed. Besides the dissociative adsorption a molecular adsorption channel exists also on the V-terminated V2O3(0001) surface, so that the water monolayer consists of both OH and molecular H2O species. The V surface layer on V2O3 is very reactive and is reoxidised by adsorbed water at 250 K, yielding surface vanadyl species. The results of this study indicate that V surface centres are necessary for the dissociation of water on V2O3 surfaces.