Anisotropy of the physisorption interaction between H2 and metal surfaces

Hindered rotational physisorption states of H2 on Ag(111) surfaces

We have investigated the physisorption states of H2 on Ag(111) surfaces. To clarify the accurate adsorption properties of H2 on Ag(111), we performed first-principles calculations based on spin-polarized density functional theory (DFT) with the semiempirical DFT-D2 method and the newly-developed exchange functional with the non-local correlation functional vdW-DF2 (rev-vdW-DF2). We constructed exhaustive potential energy surfaces, and revealed that non-negligible out-of-plane potential anisotropy with a perpendicular orientation preference exists even for H2 physisorption on planar Ag(111), as predicted by previous results of resonance-enhanced multiphoton ionization spectroscopy and temperature-programmed desorption experiments. Therefore, the molecular rotational ground states of ortho-H2 split into two energy levels in the anisotropic potential. The obtained adsorption energy and the number of bound states, including the zero-point energies and the rotational energy shift, agree with diffractive and rotationally mediated selective adsorption scattering resonance measurements. The origin of the potential anisotropy on Ag(111) is a combination of the London dispersion interaction and the virtual transition of the metal electron to the unoccupied molecular state.

Effects of rotational-symmetry breaking on physisorption of ortho- and para-H2 on Ag(111)

Quantum-state-selective thermal desorption of H2 weakly physisorbed on Ag(111) demonstrates significantly different desorption features between the nuclear-spin modifications. An energy shift due to the rotational-symmetry breaking induced by an anisotropic interaction affects not only the enthalpy but also the entropy of adsorption. The preexponential factor for desorption of the ortho-H2 is about three times as large as that of the para-H2. The entropy difference indicates a perpendicular orientation preference of anisotropic physisorption potential, which also suggests the importance of partial hybridization interaction for weak physisorption.

Benchmarking van der Waals density functionals with experimental data: potential-energy curves for H2 molecules on Cu(111), (100) and (110) surfaces

Detailed physisorption data from experiment for the H(2) molecule on low-index Cu surfaces challenge theory. Recently, density functional theory (DFT) has been developed to account for nonlocal correlation effects, including van der Waals (dispersion) forces. We show that the functional vdW-DF2 gives a potential-energy curve, potential-well energy levels and difference in lateral corrugation promisingly close to the results obtained by resonant elastic backscattering-diffraction experiments. The backscattering barrier is sensitive to the choice of exchange functional approximation. Further, the DFT-D3 and TS-vdW corrections to traditional DFT formulations are also benchmarked, and deviations are analyzed.

Dissociative chemisorption of H2 on the Cu(110) surface: a quantum and quasiclassical dynamical study

Six-dimensional quantum dynamical and quasiclassical trajectory (QCT) calculations are reported for the reaction and vibrationally inelastic scattering of (v = 0,1,j = 0) H(2) scattering from Cu(110), and for the reaction and rovibrationally elastic and inelastic scattering of (v = 1,j = 1) H(2) scattering from Cu(110). The dynamics results were obtained using a potential energy surface obtained with density functional theory using the PW91 functional. The reaction probabilities computed with quantum dynamics for (v = 0,1,j = 0) were in excellent agreement with the QCT results obtained earlier for these states, thereby validating the QCT approach to sticking of hydrogen on Cu(110). The vibrational de-excitation probability P(v=1,j = 0 --> v = 0) computed with the QCT method is in remarkably good agreement with the quantum dynamical results for normal incidence energies E(n) between 0.2 and 0.6 eV. The QCT result for the vibrational excitation probability P(v = 0,j = 0 --> v = 1) is likewise accurate for E(n) between 0.8 and 1 eV, but the QCT method overestimates vibrational excitation for lower E(n). The QCT method gives probabilities for rovibrationally (in)elastic scattering, P(v = 1,j = 1 --> v('),j(')), which are in remarkably good agreement with quantum dynamical results. The rotationally averaged, initial vibrational state-selective reaction probability obtained with QCT agrees well with the initial vibrational state-selective reaction probability extracted from molecular beam experiments for v = 1, for the range of collision energies for which the v=1 contribution to the measured total sticking probability dominates. The quantum dynamical probabilities for rovibrationally elastic scattering of (v = 1,j = 1) H(2) from Cu(110) are in good agreement with experiment for E(n) between 0.08 and 0.25 eV.