A first-principles study of surface and subsurface H on and in Ni(1 1 1): diffusional properties and coverage-dependent behavior

Periodic, self-consistent, density functional theory (GGA-PW91) calculations are performed for both surface and subsurface atomic hydrogen on and in Ni(1 1 1). At a low coverage (θ=0.25 ML), the binding energies (BEs) of a hydrogen atom in surface fcc, subsurface octahedral (first layer), and subsurface octahedral (second layer) sites are −2.89, −2.18, and −2.11 eV, respectively. The activation energy barriers for hydrogen diffusion from the surface to the first subsurface layer and from the first to the second subsurface layer are estimated to be 0.88 and 0.52 eV, respectively. In the entire coverage range studied, hydrogen occupies surface fcc and hcp sites and subsurface octahedral sites. In addition, the magnitude of the BE per hydrogen atom and the magnetization of the nickel slabs both decrease as hydrogen coverage increases. Vibrational frequencies of hydrogen at various surface and subsurface sites are calculated and are in reasonable agreement with experimental data. A phase stability calculation with a 2 × 2 surface unit cell shows that a p(2 × 2)-2H overlayer structure (θ=0.5 ML) and a p(1 × 1)-1H structure (θ=1.0 ML) are stable at low hydrogen pressures, in agreement with numerous experimental results. A very large increase in pressure is required to populate subsurface sites. After such an increase occurs, the first subsurface layer is filled completely.