M. Blanco-Rey

Methane dissociation and methyl diffusion on PdO{100}

The CH(4)-->CH(3)+H dissociation on the defect-free Pd-2O terminated PdO{100} surface is investigated by first principles calculations. CH(3) chemisorbs preferentially at O sites, where the adsorption is 0.68 eV more stable than at Pd sites. Analysis of the electronic structure of adsorption geometries reveals that the enhanced stability is due to the strongly covalent character of the bonding between surface O and CH(3). However, transition state searches suggest that the Pd atom is the most likely active site for CH(4) dissociation. An alternative route toward the most stable chemisorbed state is proposed that involves diffusion of CH(3) on the surface.

Influence of the van der Waals interaction in the dissociation dynamics of N2 on W(110) from first principles

Using ab initio molecular dynamics (AIMD) calculations, we investigate the role of the van der Waals (vdW) interaction in the dissociative adsorption of N2 on W(110). Hitherto, existing classical dynamics calculations performed on six-dimensional potential energy surfaces based on density functional theory (DFT), and the semi-local PW91 and RPBE [Hammer et al. Phys. Rev. B 59, 7413 (1999)] exchange-correlation functionals were unable to fully describe the dependence of the initial sticking coefficient on the molecular beam incidence conditions as found in experiments. N2 dissociation on W(110) was shown to be very sensitive not only to short molecule-surface distances but also to large distances where the vdW interaction, not included in semilocal-DFT, should dominate. In this work, we perform a systematic study on the dissociative adsorption using a selection of existing non-local functionals that include the vdW interaction (vdW-functionals). Clearly, the inclusion of the non-local correlation term contributes in all cases to correct the unrealistic energy barriers that were identified in the RPBE at large molecule-surface distances. Among the tested vdW-functionals, the original vdW-DF by Dion et al. [Phys. Rev. Lett. 92, 246401 (2004)] and the ulterior vdW-DF2 give also an adequate description of the N2 adsorption energy and energy barrier at the transition state, i.e., of the properties that are commonly used to verify the quality of any exchange-correlation functional. However, the results of our AIMD calculations, which are performed at different incidence conditions and hence extensively probe the multi-configurational potential energy surface of the system, do not seem as satisfactory as the preliminary static analysis suggested. When comparing the obtained dissociation probabilities with existing experimental data, none of the used vdW-functionals seems to provide altogether an adequate description of the N2/W(110) interaction at short and large distances.

Vibrational lifetimes of hydrogen on lead films: An ab initio molecular dynamics with electronic friction (AIMDEF) study

Using density functional theory and Ab Initio Molecular Dynamics with Electronic Friction (AIMDEF), we study the adsorption and dissipative vibrational dynamics of hydrogen atoms chemisorbed on free-standing lead films of increasing thickness. Lead films are known for their oscillatory behaviour of certain properties with increasing thickness, e.g., energy and electron spillout change in discontinuous manner, due to quantum size effects [G. Materzanini, P. Saalfrank, and P. J. D. Lindan, Phys. Rev. B 63, 235405 (2001)]. Here, we demonstrate that oscillatory features arise also for hydrogen when chemisorbed on lead films. Besides stationary properties of the adsorbate, we concentrate on finite vibrational lifetimes of H-surface vibrations. As shown by AIMDEF, the damping via vibration-electron hole pair coupling dominates clearly over the vibration-phonon channel, in particular for high-frequency modes. Vibrational relaxation times are a characteristic function of layer thickness due to the oscillating behaviour of the embedding surface electronic density. Implications derived from AIMDEF for frictional many-atom dynamics, and physisorbed species will also be given.

Surface stress in d-band metal surfaces

Using calculations within the density functional approach, we make a survey of the surface stresses of relevant d-band metal surfaces. Unlike the case for surface energies, where volcano-shaped dependences have been found across the d period, the surface stress values show a non-trivial dependence on the metal species and surface termination, which cannot be anticipated from homogeneous electron gas models. A qualitative interpretation as regards the origin of this dependence can be given on the basis of the decomposition of the surface stress into repulsive and attractive electronic components.

Diffusion of Hydrogen in Pd Assisted by Inelastic Ballistic Hot Electrons

Sykes et al. [Proc. Natl. Acad. Sci. U.S.A. 102, 17907 (2005)] have reported how electrons injected from a scanning tunneling microscope modify the diffusion rates of H buried beneath Pd(111). A key point in that experiment is the symmetry between positive and negative voltages for H extraction, which is difficult to explain in view of the large asymmetry in Pd between the electron and hole densities of states. Combining concepts from the theory of ballistic electron microscopy and electron-phonon scattering we show that H diffusion is driven by the s-band electrons only, which explains the observed symmetry.

Determination of the surface structure of CeO2(111) by low-energy electron diffraction

We determine the atomic structure of the (111) surface of an epitaxial ceria film using low-energy electron diffraction (LEED). The 3-fold-symmetric LEED patterns are consistent with a bulk-like termination of the (111) surface. By comparing the experimental dependence of diffraction intensity on electron energy (LEED-I(V) data) with simulations of dynamic scattering from different surface structures, we find that the CeO2(111) surface is terminated by a plane of oxygen atoms. We also find that the bond lengths in the top few surface layers of CeO2(111) are mostly undistorted from their bulk values, in general agreement with theoretical predictions. However, the topmost oxygen layer is further from the underlying cerium layer than the true bulk termination, an expansion that differs from theoretical predictions.

On the tautomerisation of porphycene on copper (111)

We use density-functional theory (DFT) to analyse the interaction of trans- and cis-porphycene with Cu(111) and their interconversion by intramolecular H-transfer. This tautomerisation reaction is characterised by small values for the reaction energy and barrier, on the order of ∼0.1 eV, where the trans configuration is thermodynamically more stable upon adsorption according to the experiments [J. N. Ladenthin et al., ACS Nano 9, 7287-7295 (2015)]. To gain even a qualitatively correct description of this reaction at the DFT level, an accurate treatment of dispersion interactions and a careful choice of the exchange contribution are required in order to predict the subtle energetics. Analysis of the electronic structure shows that adsorption is contributed by a van der Waals (vdW) interaction, mainly responsible for stabilising the polyaromatic fragments, and by a significant charge redistribution localised between Cu and the unsaturated N atoms of the molecule central cavity. We find that different vdW functionals can produce qualitatively different electronic structures, while yielding small trans vs. cis energy differences. Unlike other functionals surveyed here, vdW-DF with PBE exchange satisfactorily reproduces not only the experimental energetics but also the scanning tunneling microscopy images. This gives us confidence that this functional achieves a reliable balance between the two mechanisms contributing to the adsorption of porphycene.

Quantitative LEED analysis using a simultaneous optimization algorithm

The performance of a combinatorial simultaneous optimization (SO) algorithm is tested using experimental LEED I (E) data from Cu(100) and Fe0.57Al0.47(100) surfaces. SO optimizes structures taking advantage of the experimental database at two levels: (i) commensurate subsets of the database with the number of unknown parameters are chosen to find local solutions using Broyden’s method and (ii) these partial structural solutions are used to build a Markov chain over the whole database. This procedure is of global character, the same as simulated annealing or genetic algorithm methods, but displays a very competitive scaling law because after the first iteration candidates are not chosen by a blind/random pick; they are already solutions to the problem with a restricted experimental database. (Some figures in this article are in colour only in the electronic version)

Diffusion of hydrogen interstitials in the near-surface region of Pd(111) under the influence of surface coverage and external static electric fields

Past scanning tunneling microscopy (STM) experiments of H manipulation on Pd(111), at low temperature, have shown that it is possible to induce diffusion of surface species as well as of those deeply buried under the surface. Several questions remain open regarding the role of subsurface site occupancies. In the present work, the interaction potential of H atoms with Pd(111) under various H coverage conditions is determined by means of density functional theory calculations in order to provide an answer to two of these questions: (i) whether subsurface sites are the final locations for the H impurities that attempt to emerge from bulk regions, and (ii) whether penetration of the surface is a competing route of on-surface diffusion during depletion of surface H on densely covered Pd(111). We find that a high H coverage has the effect of blocking resurfacing of H atoms travelling from below, which would otherwise reach the surface fcc sites, but it hardly alters deeper diffusion energy barriers. Penetration is unlikely and restricted to high occupancies of hcp hollows. In agreement with experiments, the Pd lattice expands vertically as a consequence of H atoms being blocked at subsurface sites, and surface H enhances this expansion. STM tip effects are included in the calculations self-consistently as an external static electric field. The main contribution to the induced surface electric dipoles originates from the Pd substrate polarisability. We find that the electric field has a non-negligible effect on the H-Pd potential in the vicinity of the topmost Pd atomic layer, yet typical STM intensities of 1-2 VÅ(-1) are insufficient to invert the stabilities of the surface and subsurface equilibrium sites.

A FORTRAN-90 Low-Energy Electron Diffraction program (LEED90 v1.1)

We describe a FORTRAN-90 program to compute low-energy electron diffraction I(V) curves. Plane-waves and layer doubling are used to compute the inter-layer multiple-scattering, while the intra-layer multiple-scattering is computed in the standard way expanding the wavefield on a basis of spherical waves. The program is kept as general as possible, in order to allow testing different parts of multiple-scattering calculations. In particular, it can handle non-diagonal t-matrices describing the scattering of non-spherical potentials, anisotropic vibrations, anharmonicity, etc. The program does not use old FORTRAN flavours, and has been written keeping in mind the advantage for parallelism brought forward by FORTRAN-90.