Repulsive interaction of the helium atom with a metal surface

The repulsive part of the helium scattering potential at a surface is approximately proportional to the surface electron density. The proportionality coefficient is shown to be a well-defined quantity, which can be related to the electron-helium scattering length. The spread in the values of the proportionality constant suggested in the literature is shown to be due to different definitions of the coefficient or due to inadequate calculational methods. The value calculated using the local density approximation with a self-interaction correction is in very good agreement with the electron-scattering-length measurements.Peer reviewe

Bayesian Error Estimation in Density Functional Theory

We present a practical scheme for performing error estimates for Density Functional Theory calculations. The approach which is based on ideas from Bayesian statistics involves creating an ensemble of exchange-correlation functionals by comparing with an experimental database of binding energies for molecules and solids. Fluctuations within the ensemble can then be used to estimate errors relative to experiment on calculated quantities like binding energies, bond lengths, and vibrational frequencies. It is demonstrated that the error bars on energy differences may vary by orders of magnitude for different systems in good agreement with existing experience.Comment: 5 pages, 3 figure

Quantum Motion of Chemisorbed Hydrogen on Ni Surfaces

Quantum mechanical energy levels and wave functions have been calculated for the motion of chemisorbed hydrogen atoms on Ni surfaces. The results show considerable quantum effects for the adatom in both the ground and the excited states. The description of the adparticles as being delocalized along the surface offers a novel interpretation of several phenomena, in particular the vibrational excitations.Peer reviewe

The Path Integral Monte Carlo Calculation of Electronic Forces

We describe a method to evaluate electronic forces by Path Integral Monte Carlo (PIMC). Electronic correlations, as well as thermal effects, are included naturally in this method. For fermions, a restricted approach is used to avoid the ``sign'' problem. The PIMC force estimator is local and has a finite variance. We applied this method to determine the bond length of H$_2$ and the chemical reaction barrier of H+H$_2\longrightarrow$H$_2$+H. At low temperature, good agreement is obtained with ground state calculations. We studied the proton-proton interaction in an electron gas as a simple model for hydrogen impurities in metals. We calculated the force between the two protons at two electronic densities corresponding to Na ($r_s=3.93$) and Al ($r_s=2.07$) using a supercell with 38 electrons. The result is compared to previous calculations. We also studied the effect of temperature on the proton-proton interaction. At very high temperature, our result agrees with the Debye screening of electrons. As temperature decreases, the Debye theory fails both because of the strong degeneracy of electrons and most importantly, the formation of electronic bound states around the protons.Comment: 18 pages, 10 figure