Steering and ro-vibrational effects in the dissociative adsorption and associative desorption of H_2/Pd(100)

The interaction of hydrogen with many transition metal surfaces is characterized by a coexistence of activated with non-activated paths to adsorption with a broad distribution of barrier heights. By performing six-dimensional quantum dynamical calculations using a potential energy surface derived from ab initio calculations for the system H_2/Pd(100) we show that these features of the potential energy surface lead to strong steering effects in the dissociative adsorption and associative desorption dynamics. In particular, we focus on the coupling of the translational, rotational and vibrational degrees of freedom of the hydrogen molecule in the reaction dynamics.Comment: 8 pages, 5 figures, subm. to the Proceedings of ISSP-18, June 1996, Poland, to appear in Prog. Surf. Sc

Mechanism of Poisoning the Catalytic Activity of Pd(100) by a Sulfur Adlayer

The modification of the potential-energy surface (PES) of H_2 dissociation over Pd(100) as induced by the presence of a (2x2) S adlayer is investigated by density-functional theory and the linear augmented plane wave method. It is shown that the poisoning effect of S originates from the formation of energy barriers hindering the dissociation of H_2. The barriers are in the entrance channel of the PES and their magnitude strongly depends on the lateral distance of the H_2 molecule from the S adatoms.Comment: RevTeX, 14 pages, 3 figure

Theory of self-diffusion in GaAs

Ab initio molecular dynamics simulations are employed to investigate the dominant migration mechanism of the gallium vacancy in gaas as well as to assess its free energy of formation and the rate constant of gallium self-diffusion. our analysis suggests that the vacancy migrates by second nearest neighbour hops. the calculated self-diffusion constant is in good agreement with the experimental value obtained in ^69 GaAs/ ^71 GaAs isotope heterostructures and at significant variance with that obtained earlier from interdiffusion experiments in GaAlAs/GaAs-heterostructures.Comment: 15 pages, 4 figures. Z. Phys. Chem, in prin

Effect of a humid environment on the surface structure of RuO2(110)

Combining density-functional theory and thermodynamics we compute the phase diagram of surface structures of rutile RuO2 (110) in equilibrium with water vapor in the complete range of experimentally accessible gas phase conditions. Through the formation of hydroxyl or water-like groups, already lowest concentrations of hydrogen in the gas phase are sufficient to stabilize an oxygen-rich polar oxide termination even at very low oxygen pressure.Comment: 7 pages including 5 figures, Submitted to Phys. Rev. B., Related publications can be found at http://www.fhi-berlin.mpg.de/th/paper.htm

Bridging the length and time scales: from ab initio electronic structure calculations to macroscopic proportions

Density functional theory (DFT) primarily provides a good description of the electronic structure. Thus, DFT primarily deals with length scales as those of a chemical bond, i.e. 10^-10 meter, and with time scales of the order of atomic vibrations, i.e. 10^-13 seconds. However, several interesting phenomena happen and/or become observable on different scales, namely meso- or macroscopic lengths and on time scales of seconds or even minutes. To bridge the gap between 10^-13 seconds and a second or between 10^-10 meter and 10 and more nano meters is one of the important challenges we are facing today. In this paper we show how we are overcoming these time and size problems for the example of crystal growth and the evolution of nano-scale structures. The key is a kinetic Monte Carlo approach with detailed input from DFT calculations of the relevant atomistic processes.Comment: 13 pages, 5 figures, to be published in Comments on Condens. Matt. Phys. (1998). Other related publications can be found at http://www.rz-berlin.mpg.de/th/paper.htm

Impact of vibrational entropy on the stability of unsolvated peptide helices with increasing length

Helices are a key folding motif in protein structure. The question which factors determine helix stability for a given polypeptide or protein is an ongoing challenge. Here we use van der Waals corrected density-functional theory to address a part of this question in a bottom-up approach. We show how intrinsic helical structure is stabilized with length and temperature for a series of experimentally well studied unsolvated alanine based polypeptides, Ac-Alan-LysH+. By exploring extensively the conformational space of these molecules, we find that helices emerge as the preferred structure in the length range n=4-8 not just due to enthalpic factors (hydrogen bonds and their cooperativity, van der Waals dispersion interactions, electrostatics), but importantly also by a vibrational entropic stabilization over competing conformers at room temperature. The stabilization is shown to be due to softer low-frequency vibrational modes in helical conformers than in more compact ones. This observation is corroborated by including anharmonic effects explicitly through \emph{ab initio} molecular dynamics, and generalized by testing different terminations and considering larger helical peptide models

A Model for the Thermal Expansion of Ag(111) and other Metal Surfaces

We develop a model to study the thermal expansion of surfaces, wherein phonon frequencies are obtained from ab initio total energy calculations. Anharmonic effects are treated exactly in the direction normal to the surface, and within a quasiharmonic approximation in the plane of the surface. We apply this model to the Ag(111) and Al(111) surfaces, and find that our calculations reproduce the experimental observation of a large and anomalous increase in the surface thermal expansion of Ag(111) at high temperatures [P. Statiris, H.C. Lu and T. Gustafsson, Phys. Rev. Lett. 72, 3574 (1994)]. Surprisingly, we find that this increase can be attributed to a rapid softening of the in-plane phonon frequencies, rather than due to the anharmonicity of the out-of-plane surface phonon modes. This provides evidence for a new mechanism for the enhancement of surface anharmonicity. A comparison with Al(111) shows that the two surfaces behave quite differently, with no evidence for such anomalous behavior on Al(111).Comment: 17 pages, 4 figures, to appear in Z. Chem. Phy

Multi-Lattice Kinetic Monte Carlo Simulations from First-Principles: Reduction of the Pd(100) Surface Oxide by CO

We present a multi-lattice kinetic Monte Carlo (kMC) approach that efficiently describes the atomistic dynamics of morphological transitions between commensurate structures at crystal surfaces. As an example we study the reduction of a $(\sqrt{5}\times \sqrt{5})R27^{\circ}$ PdO(101) overlayer on Pd(100) in a CO atmosphere. Extensive density-functional theory calculations are used to establish an atomistic pathway for the oxide reduction process. First-principles multi-lattice kMC simulations on the basis of this pathway fully reproduce the experimental temperature dependence of the reduction rate [Fernandes et al., Surf. Sci. 2014, 621, 31-39] and highlight the crucial role of elementary processes special to the boundary between oxide and metal domains.Comment: 19 pages, 10 figure