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

Role of hydroxyl groups in the NH(x) (x = 1-3) adsorption on the TiO2 anatase (101) surface determined by a first-principles study

A spin-polarized density functional theory calculation was carried out to study the adsorption of NH(x) species (x = 1-3) on a TiO2 anatase (101) surface with and without hydroxyl groups by using first-principles calculations. It was found that the present hydroxyl group has the effect of significantly enhancing the adsorption of monodentate adsorbates H2N-Ti(a) compared to that on a bare surface. The nature of the interaction between the adsorbate (NH(x)) and the hydroxylated or bare surface was analyzed by the Mulliken charge and density of states (DOS) calculations. This facilitation of NH2 is caused by the donation of coadsorbed H filling the nonbonding orbital of NH2, resulting in an electron gain in NH2 from the bonding. In addition, the upper valence band, which originally consisted of the mixing of O 2p and Ti 3d orbitals, has been broadened by the two adjacent H 1s and NH2 sigma(y)(b) orbitals joined to the bottom of the original TiO2 valence band. The results are important to understand the OH effect in heterogeneous catalysis.

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].

Adaptive Partitioning in Combined Quantum Mechanical and Molecular Mechanical Calculations of Potential Energy Functions for Multiscale Simulations

In many applications of multilevel/multiscale methods, an active zone must be modeled by a high-level electronic structure method, while a larger environmental zone can be safely modeled by a lower-level electronic structure method, molecular mechanics, or an analytic potential energy function. In some cases though, the active zone must be redefined as a function of simulation time. Examples include a reactive moiety diffusing through a liquid or solid, a dislocation propagating through a material, or solvent molecules in a second coordination sphere (which is environmental) exchanging with solvent molecules in an active first coordination shell. In this article, we present a procedure for combining the levels smoothly and efficiently in such systems in which atoms or groups of atoms move between high-level and low-level zones. The method dynamically partitions the system into the high-level and low-level zones and, unlike previous algorithms, removes all discontinuities in the potential energy and force whenever atoms or groups of atoms cross boundaries and change zones. The new adaptive partitioning (AP) method is compared to Rode's "hot spot" method and Morokuma's "ONIOM-XS" method that were designed for multilevel molecular dynamics (MD) simulations. MD simulations in the microcanonical ensemble show that the AP method conserves both total energy and momentum, while the ONIOM-XS method fails to conserve total energy and the hot spot method fails to conserve both total energy and momentum. Two versions of the AP method are presented, one scaling as O(2N) and one with linear scaling in N, where N is the number of groups in a buffer zone separating the active high-level zone from the environmental low-level zone. The AP method is also extended to systems with multiple high-level zones to allow, for example, the study of ions and counterions in solution using the multilevel approach.