J. I. Juaristi

Electronic Stopping Power in LiF from First Principles

Using time-dependent density-functional theory we calculate from first principles the rate of energy transfer from a moving proton or antiproton to the electrons of an insulating material, LiF. The behavior of the electronic stopping power versus projectile velocity displays an effective threshold velocity of approximately 0.2 a.u. for the proton, consistent with recent experimental observations, and also for the antiproton. The calculated proton/antiproton stopping-power ratio is approximately 2.4 at velocities slightly above the threshold (v approximately 0.4 a.u.), as compared to the experimental value of 2.1. The projectile energy loss mechanism is observed to be extremely local.

Vibrational deexcitation and rotational excitation of H2 and D2 scattered from Cu(111): Adiabatic versus non-adiabatic dynamics

We have studied survival and rotational excitation probabilities of H(2)(v(i) = 1, J(i) = 1) and D(2)(v(i) = 1, J(i) = 2) upon scattering from Cu(111) using six-dimensional (6D) adiabatic (quantum and quasi-classical) and non-adiabatic (quasi-classical) dynamics. Non-adiabatic dynamics, based on a friction model, has been used to analyze the role of electron-hole pair excitations. Comparison between adiabatic and non-adiabatic calculations reveals a smaller influence of non-adiabatic effects on the energy dependence of the vibrational deexcitation mechanism than previously suggested by low-dimensional dynamics calculations. Specifically, we show that 6D adiabatic dynamics can account for the increase of vibrational deexcitation as a function of the incidence energy, as well as for the isotope effect observed experimentally in the energy dependence for H(2)(D(2))/Cu(100). Furthermore, a detailed analysis, based on classical trajectories, reveals that in trajectories leading to vibrational deexcitation, the minimum classical turning point is close to the top site, reflecting the multidimensionally of this mechanism. On this site, the reaction path curvature favors vibrational inelastic scattering. Finally, we show that the probability for a molecule to get close to the top site is higher for H(2) than for D(2), which explains the isotope effect found experimentally.

Dissipative effects in the dynamics of N2 on tungsten surfaces

The role of electron-hole pair excitations in the dynamics of N(2) on W(100) and W(110) is evaluated using a theoretical model that accounts for the six-dimensionality of the problem in the whole calculation. The six-dimensional potential energy surface is determined in each case from an extensive grid of energies calculated with density functional theory. Dissipative effects due to electron-hole pair excitations are introduced in the classical dynamics equations through a friction force. Corresponding electron friction coefficients are calculated for each atom in the molecule with density functional theory in a local density approximation. Our results show that electronic friction plays a very minor role in the dissociative dynamics of N(2) in both tungsten faces. A similar conclusion is reached when we calculate the energy lost by the reflecting molecules.

Electronic Friction-Based Vibrational Lifetimes of Molecular Adsorbates : Beyond the Independent-Atom Approximation

We assess the accuracy of vibrational damping rates of diatomic adsorbates on metal surfaces as calculated within the local-density friction approximation (LDFA). An atoms-in-molecules (AIM) type charge partitioning scheme accounts for intramolecular contributions and overcomes the systematic underestimation of the nonadiabatic losses obtained within the prevalent independent-atom approximation. The quantitative agreement obtained with theoretical and experimental benchmark data suggests the LDFA-AIM scheme as an efficient and reliable approach to account for electronic dissipation in ab initio molecular dynamics simulations of surface chemical reactions.

Role of physisorption states in molecular scattering: A semilocal density-functional theory study on O2/Ag (111)

We simulate the scattering of O2 from Ag(111) with classical dynamics simulations performed on a six-dimensional potential energy surface calculated within semilocal density-functional theory. The enigmatic experimental trends that originally required the conjecture of two types of repulsive walls, arising from a physisorption and chemisorption part of the interaction potential, are fully reproduced. Given the inadequate description of the physisorption properties in semilocal density-functional theory, our work casts severe doubts on the prevalent notion to use molecular scattering data as indirect evidence for the existence of such states.

Charge state dependence of the energy loss of slow nitrogen ions reflected from an aluminum surface under grazing incidence

The experimentally observed monotonic increase of the energy loss with charge q for N qa ions impinging on an Al(1 1 1) surface under grazing angle of incidence is explained by a model that accounts for the eAect of K- and L-shell vacancies in the stopping power. Our model allows us to estimate the characteristic time scales (and distances from the surface) for the neutralization and relaxation of multicharged ions. We use a transport cross section formulation of the electronic stopping of ions in an electron gas, as well as a self-consistent calculation of the scattering potential within density functional theory. ” 1999 Elsevier Science B.V. All rights reserved.

Nonlinear screening effects in the interaction of slow multicharged ions with metal surfaces

Abstract Density functional theory is used to treat the problem of screening of highly charged ions in an electron gas when inner shell vacancies are present. Results are presented for N and Ne ions in various electronic configurations. These include information about screening charge densities, induced densities of states in the conduction band, bound energy levels and KLL Auger lines.

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.

Femtosecond-laser-driven molecular dynamics on surfaces: Photodesorption of molecular oxygen from Ag(110)

We simulate the femtosecond laser induced desorption dynamics of a diatomic molecule from a metal surface by including the effect of the electron and phonon excitations created by the laser pulse. Following previous models, the laser induced surface excitation is treated through the two temperature model, while the multidimensional dynamics of the molecule is described by a classical Langevin equation, in which the friction and random forces account for the action of the heated electrons. In this work, we propose the additional use of the generalized Langevin oscillator model to also include the effect of the energy exchange between the molecule and the heated surface lattice in the desorption dynamics. The model is applied to study the laser induced desorption of O