Stepwise addition reactions in ammonia synthesis: A first principles study

Catalytic ammonia synthesis is believed to proceed via dissociation of N-2 and H-2 with subsequent stepwise addition reactions from an adsorbed nitrogen atom to NH3. The first step, N-2 dissociation, has been thoroughly studied. However, little is known about the microscopic details of the stepwise addition reactions. To shed light on these stepwise addition reactions, density functional theory calculations with the generalized gradient approximation are employed to investigate NHx (x=1,3) formation on Ru(0001). Transition states and reaction barriers are determined in each elementary step. It is found that the reaction barriers for stepwise addition reactions are rather high, for example, the barrier for NH hydrogenation is calculated to be 1.28 eV, which is comparable with that of N-2 dissociation. In addition, one of the stepwise addition reactions on a stepped surface is also considered. The reaction barrier is found to be much higher than that of N-2 dissociation on the same stepped surface, which indicates the importance of stepwise addition reactions in ammonia synthesis. (C) 2001 American Institute of Physics.

CO oxidation on Pd(100) and Pd(111): A comparative study of reaction pathways and reactivity at low and medium coverages

We have performed density functional theory calculations with the generalized gradient approximation to investigate CO oxidation on a close-packed transition metal surface, Pd(lll), and a more open surface, Pd(100), aiming to shed light on surface structure effects on reaction pathways and reactivity, an important issue in catalysis. Reaction pathways on both surfaces at two different coverages have been studied. It is found that the reaction pathways on both surfaces possess crucial common features despite the fact that they have different surface symmetries. Having determined reaction barriers in these systems, we find that the reaction on Pd(lll) is strongly coverage dependent. Surface coverages, however, have little effect on the reaction on Pd(100). Calculations also reveal that the low coverage reactions are structure sensitive while the medium coverage reactions are not. Detailed discussions on these results are given.

Reduced-size polarized basis sets for calculations of molecular electric properties. III. Second-row atoms

Reduced-size polarized (ZmPolX) basis sets are developed for the second-row atoms X = Si, P, S, and Cl. The generation of these basis sets follows from a simple physical model of the polarization effect of the external electric field which leads to highly compact polarization functions to be added to the chosen initial basis set. The performance of the ZmPolX sets has been investigated in calculations of molecular dipole moments and polarizabilities. Only a small deterioration of the quality of the calculated molecular electric properties has been found. Simultaneously the size of the present reduced-size ZmPolX basis sets is about one-third smaller than that of the usual polarized (PolX) sets. This reduction considerably widens the range of applications of the ZmPolX sets in calculations of molecular dipole moments, dipole polarizabilities, and related properties.

Effect of quantization of vibrations on the structural properties of crystals

We study the structural effects produced by the quantization of vibrational degrees of freedom in periodic crystals at zero temperature. To this end we introduce a methodology based on mapping a suitable subspace of the vibrational manifold and solving the Schrödinger equation in it. A number of increasingly accurate approximations ranging from the quasiharmonic approximation (QHA) to the vibrational self-consistent field (VSCF) method and the exact solution are described. A thorough analysis of the approximations is presented for model monatomic and hydrogen-bonded chains, and results are presented for a linear H-F chain where the potential-energy surface is obtained via first-principles electronic structure calculations. We focus on quantum nuclear effects on the lattice constant and show that the VSCF is an excellent approximation, meaning that correlation between modes is not extremely important. The QHA is excellent for covalently bonded mildly anharmonic systems, but it fails for hydrogen-bonded ones. In the latter, the zero-point energy exhibits a nonanalytic behavior at the lattice constant where the H atoms center, which leads to a spurious secondary minimum in the quantum-corrected energy curve. An inexpensive anharmonic approximation of noninteracting modes appears to produce rather good results for hydrogen-bonded chains for small system sizes. However, it converges to the incorrect QHA results for increasing size. Isotope effects are studied for the first-principles H-F chain. We show how the lattice constant and the H-F distance increase with decreasing mass and how the QHA proves to be insufficient to reproduce this behavior.

Upper and lower bounds of 180 degrees unit element of a ridge waveguide: Calculations and measurements

The upper and lower bounds on the actual solution of any microwave structure is of general interest. The purpose of this letter is to compare some calculations using the mode-matching and finite-element methods, with some measurements on a 180 degrees ridge waveguide insert between standard WR62 rectangular waveguides. The work suggests that the MMM produces an upper bound, while the FEM places a lower bound on the measurement. (C) 2001 John Wiley & Sons, Inc.

Crystal structure and ab initio calculations of a cyano-carbamimidic acid ethyl ester

The structure of one tautomer (amine form) of cyano-carbamimidic acid ethyl ester or (amino-ethoxy-methylidene)aminoformonitrile (CAS: 13947-84-7) was determined by single crystal X-ray diffraction. Ab initio quantum chemical calculations at the B3LYP, MP2 and G3 levels were performed to investigate the stability and the formation of the different tautomers and conformers. The calculations indicate that the amine form is the more stable tautomer, showing a high degree of election conjugation. The most stable amine conformer located by the calculations corresponds to the crystallized structure. On the contrary, in the less stable imine form, the conjugation is separated by a N2-C2 single bond. (C) 2007 Elsevier B.V. All rights reserved.

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.

Electron-impact excitation of O II fine-structure levels.

Effective collision strengths for forbidden transitions among the five energetically lowest fine-structure levels of O ii are calculated in the Breit-Pauli approximation using the R-matrix method. Results are presented for the electron temperature range 100-100 000 K. The accuracy of the calculations is evaluated via the use of different types of radial orbital sets and a different configuration expansion basis for the target wavefunctions. A detailed assessment of previous available data is given, and erroneous results are highlighted. Our results reconfirm the validity of the original Seaton and Osterbrock scaling for the optical O ii ratio, a matter of some recent controversy. Finally, we present plasma diagnostic diagrams using the best collision strengths and transition probabilities.

An analysis of geometric uncertainty calculations for prostate radiotherapy in clinical practice.

Margins are used in radiotherapy to assist in the calculation of planning target volumes. These margins can be determined by analysing the geometric uncertainties inherent to the radiotherapy planning and delivery process. An important part of this process is the study of electronic portal images collected throughout the course of treatment. Set-up uncertainties were determined for prostate radiotherapy treatments at our previous site and the new purpose-built centre, with margins determined using a number of different methods. In addition, the potential effect of reducing the action level from 5 mm to 3 mm for changing a patient set-up, based on off-line bony anatomy-based portal image analysis, was studied. Margins generated using different methodologies were comparable. It was found that set-up errors were reduced following relocation to the new centre. Although a significant increase in the number of corrections to a patient's set-up was predicted if the action level was reduced from 5 mm to 3 mm, minimal reduction in patient set-up uncertainties would be seen as a consequence. Prescriptive geometric uncertainty analysis not only supports calculation and justification of the margins used clinically to generate planning target volumes, but may also best be used to monitor trends in clinical practice or audit changes introduced by new equipment, technology or practice. Simulations on existing data showed that a 3 mm rather than a 5 mm action level during off-line, bony anatomy-based portal imaging would have had a minimal benefit for the patients studied in this work.

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.

Inelastic quantum transport in nanostructures: The self-consistent Born approximation and correlated electron-ion dynamics

A dynamical method for inelastic transport simulations in nanostructures is compared to a steady-state method based on nonequilibrium Green's functions. A simplified form of the dynamical method produces, in the steady state in the weak-coupling limit, effective self-energies analogous to those in the Born approximation due to electron-phonon coupling. The two methods are then compared numerically on a resonant system consisting of a linear trimer weakly embedded between metal electrodes. This system exhibits an enhanced heating at high biases and long phonon equilibration times. Despite the differences in their formulation, the static and dynamical methods capture local current-induced heating and inelastic corrections to the current with good agreement over a wide range of conditions, except in the limit of very high vibrational excitations where differences begin to emerge.