Model calculations for copper clusters on Au(111) surfaces

Using the semi-empirical embedded-atom method, the structure of small copper clusters on Au(111) surfaces has been investigated both by static and dynamic calculations. By varying the size of roughly circular clusters, the edge energy per atom is obtained; it agrees quite well with estimates based on experimental results. Small three-dimensional clusters tend to have the shape of a pyramid, whose sides are oriented in the directions of small surface energy. The presence of a cluster is found to distort the underlying lattice of adsorbed copper atoms. (C) 2002 Published by Elsevier Science B.V.

Prediction of the formation and stabilities of energetic salts and ionic liquids based on ab initio electronic structure calculations

A computational approach to predict the thermodynamics for forming a variety of imidazolium-based salts and ionic liquids from typical starting materials is described. The gas-phase proton and methyl cation acidities of several protonating and methylating agents, as well as the proton and methyl cation affinities of many important methyl-, nitro-, and cyano- substituted imidazoles, have been calculated reliably by using the computationally feasible DFT (B3LYP) and MP2 (extrapolated to the complete basis set limit) methods. These accurately calculated proton and methyl cation affinities of neutrals and anions are used in conjunction with an empirical approach based on molecular volumes to estimate the lattice enthalpies and entropies of ionic liquids, organic solids, and organic liquids. These quantities were used to construct a thermodynamic cycle for salt formation to reliably predict the ability to synthesize a variety of salts including ones with potentially high energetic densities. An adjustment of the gas phase thermodynamic cycle to account for solid- and liquid-phase chemistries provides the best overall assessment of salt formation and stability. This has been applied to imidazoles (the cation to be formed) with alkyl, nitro, and cyano substituents. The proton and methyl cation donors studied were as follows: HCl, HBr, HI, (HO)(2)SO2, HSO3CF3 (TfOH), and HSO3(C6H4)CH3 (TsOH); CH3Cl, CH3Br, CH3I, (CH3O)(2)SO2, CH3SO3CF3 (TfOCH3) and CH3SO3(C6H4)CH3 (TsOCH3). As substitution of the cation with electron-withdrawing groups increases, the triflate reagents appear to be the best overall choice as protonating and methylating agents. Even stronger alkylating agents should be considered to enhance the chances of synthetic success. When using the enthalpies of reaction for the gas-phase reactants (eq 6) to form a salt, a cutoff value of - 13 kcal mol(-1) or lower (more negative) should be used as the minimum value for predicting whether a salt can be synthesized.

NO Oxidation on Platinum Group Metals Oxides:First Principles Calculations Combined with Microkinetic Analysis

By combining density functional theory calculation and microkinetic analysis, NO oxidation on the platinum group metal oxides (PtO(2), IrO(2), OsO(2)) is investigated, aiming at shedding light on the activities of metal oxides and exploring the activity variations of metal oxides compared to their corresponding metals. A microkinetic model, taking into account the possible low diffusion of surface species on metal oxide surfaces, is proposed for NO oxidation. The resultant turnover frequencies of NO oxidation show that under the typical experimental condition, T = 600 K, p(O2) = 0.1 atm, p(NO) = 3 x 10(-4) atm, p(NO2) = 1.7 x 10(-4) atm; (i) IrO(2)(110) exhibits higher activity than PtO(2)(110) and OsO(2)(110), and (ii) compared to the corresponding metallic Pt, Ir, and Os, the activity of PtO(2) to catalyze NO oxidation is lower, but interestingly IrO(2) and OsO(2) exhibit higher activities. The reasons for the activity differences between the metals and oxides are addressed. Moreover, other possible reaction pathways of NO oxidation on PtO(2)(110), involving O(2) molecule (NO + O(2) -> OONO) and lattice bridge-O(2c), are also found to give low activities. The origin of the Pt catalyst deactivation is also discussed.

SrTiO3(001)(2x1) reconstructions: First-principles calculations of surface energy and atomic structure compared with scanning tunneling microscopy images

(1x1) and (2x1) reconstructions of the (001) SrTiO3 surface were studied using the first-principles full-potential linear muffin-tin orbital method. Surface energies were calculated as a function of TiO2 chemical potential, oxygen partial pressure and temperature. The (1x1) unreconstructed surfaces were found to be energetically stable for many of the conditions considered. Under conditions of very low oxygen partial pressure the (2x1) Ti2O3 reconstruction [Martin R. Castell, Surf. Sci. 505, 1 (2002)] is stable. The question as to why STM images of the (1x1) surfaces have not been obtained was addressed by calculating charge densities for each surface. These suggest that the (2x1) reconstructions would be easier to image than the (1x1) surfaces. The possibility that the presence of oxygen vacancies would destabilise the (1x1) surfaces was also investigated. If the (1x1) surfaces are unstable then there exists the further possibility that the (2x1) DL-TiO2 reconstruction [Natasha Erdman Nature (London) 419, 55 (2002)] is stable in a TiO2-rich environment and for p(O2)>10(-18) atm.

Free energy and molecular dynamics calculations for the cubic-tetragonal phase transition in zirconia

The high-temperature cubic-tetragonal phase transition of pure stoichiometric zirconia is studied by molecular dynamics (MD) simulations and within the framework of the Landau theory of phase transformations. The interatomic forces are calculated using an empirical, self-consistent, orthogonal tight-binding model, which includes atomic polarizabilities up to the quadrupolar level. A first set of standard MD calculations shows that, on increasing temperature, one particular vibrational frequency softens. The temperature evolution of the free-energy surfaces around the phase transition is then studied with a second set of calculations. These combine the thermodynamic integration technique with constrained MD simulations. The results seem to support the thesis of a second-order phase transition but with unusual, very anharmonic behavior above the transition temperature.

Density functional calculations for polymers and clusters - progress and limitations

The applicability of density functional (DF) methods has progressed greatly since the first workshop of this series ten years ago. Applications that show both the successes and the limitations can be found in the fields of: (a) the structures of the isomers of atomic clusters. and (b) the structure of organic molecules and polymers, and their reactions with additional molecules. We shall review some of the results and the lessons to be learned from them. (C) 2001 Elsevier Science B.V. All rights reserved.