Physisorption interaction of H2 with noble-metal surfaces: A new H2-Cu potential

vdW-DF-ahcx: a range-separated van der Waals density functional hybrid

Hybrid density functionals replace a fraction of an underlying generalized-gradient approximation (GGA) exchange description with a Fock-exchange component. Range-separated hybrids (RSHs) also effectively screen the Fock-exchange component and thus open the door for characterizations of metals and adsorption at metal surfaces. The RSHs are traditionally based on a robust GGA, such as PBE (Perdew J Pet al1996Phys. Rev. Lett.773865), for example, as implemented in the HSE design (Heyd Jet al2003J. Chem. Phys.1188207). Here we define an analytical-hole (Henderson T Met al2008J. Chem. Phys.128194105) consistent-exchange RSH extension to the van der Waals density functional (vdW-DF) method (Berland Ket al2015Rep. Prog. Phys.78066501), launching vdW-DF-ahcx. We characterize the GGA-type exchange in the vdW-DF-cx version (Berland K and Hyldgaard P 2014Phys. Rev. B89035412), isolate the short-ranged exchange component, and define the new vdW-DF hybrid. We find that the performance vdW-DF-ahcx compares favorably to (dispersion-corrected) HSE for descriptions of bulk (broad molecular) properties. We also find that it provides accurate descriptions of noble-metal surface properties, including CO adsorption.

Performance of Made Simple Meta-GGA Functionals with rVV10 Nonlocal Correlation for H2 + Cu(111), D2 + Ag(111), H2 + Au(111), and D2 + Pt(111)

Accurately modeling heterogeneous catalysis requires accurate descriptions of rate-controlling elementary reactions of molecules on metal surfaces, but standard density functionals (DFs) are not accurate enough for this. The problem can be solved with the specific reaction parameter approach to density functional theory (SRP-DFT), but the transferability of SRP DFs among chemically related systems is limited. We combine the MS-PBEl, MS-B86bl, and MS-RPBEl semilocal made simple (MS) meta-generalized gradient approximation (GGA) (mGGA) DFs with rVV10 nonlocal correlation, and we evaluate their performance for the hydrogen (H2) + Cu(111), deuterium (D2) + Ag(111), H2 + Au(111), and D2 + Pt(111) gas-surface systems. The three MS mGGA DFs that have been combined with rVV10 nonlocal correlation were not fitted to reproduce particular experiments, nor has the b parameter present in rVV10 been reoptimized. Of the three DFs obtained the MS-PBEl-rVV10 DF yields an excellent description of van der Waals well geometries. The three original MS mGGA DFs gave a highly accurate description of the metals, which was comparable in quality to that obtained with the PBEsol DF. Here, we find that combining the three original MS mGGA DFs with rVV10 nonlocal correlation comes at the cost of a slightly less accurate description of the metal. However, the description of the metal obtained in this way is still better than the descriptions obtained with SRP DFs specifically optimized for individual systems. Using the Born-Oppenheimer static surface (BOSS) model, simulations of molecular beam dissociative chemisorption experiments yield chemical accuracy for the D2 + Ag(111) and D2 + Pt(111) systems. A comparison between calculated and measured E 1/2(ν, J) parameters describing associative desorption suggests chemical accuracy for the associative desorption of H2 from Au(111) as well. Our results suggest that ascending Jacob's ladder to the mGGA rung yields increasingly more accurate results for gas-surface reactions of H2 (D2) interacting with late transition metals.

Designing new SRP density functionals including non-local vdW-DF2 correlation for H2 + Cu(111) and their transferability to H2 + Ag(111), Au(111) and Pt(111)

Specific reaction parameter density functionals (SRP-DFs) that can describe dissociative chemisorption molecular beam experiments of hydrogen (H₂) on cold transition metal surfaces with chemical accuracy have so far been shown to be only transferable among different facets of the same metal, but not among different metals. We design new SRP-DFs that include non-local vdW-DF2 correlation for the H₂ + Cu(111) system, and evaluate their transferability to the highly activated H₂ + Ag(111) and H₂ + Au(111) systems and the non-activated H₂ + Pt(111) system. We design our functionals for the H₂ + Cu(111) system since it is the best studied system both theoretically and experimentally. Here we demonstrate that a SRP-DF fitted to reproduce molecular beam sticking experiments for H₂ + Cu(111) with chemical accuracy can also describe such experiments for H₂ + Pt(111) with chemical accuracy, and vice versa. Chemically accurate functionals have been obtained that perform very well with respect to reported van der Waals well geometries, and which improve the description of the metal over current generalized gradient approximation (GGA) based SRP-DFs. From a systematic comparison of our new SRP-DFs that include non-local correlation to previously developed SRP-DFs, for both activated and non-activated systems, we identify non-local correlation as a key ingredient in the construction of transferable SRP-DFs for H₂ interacting with transition metals. Our results are in excellent agreement with experiment when accurately measured observables are available. It is however clear from our analysis that, except for the H₂ + Cu(111) system, there is a need for more, more varied, and more accurately described experiments in order to further improve the design of SRP-DFs. Additionally, we confirm that, when including non-local correlation, the sticking of H₂ on Cu(111) is still well described quasi-classically.

Screening nature of the van der Waals density functional method: a review and analysis of the many-body physics foundation

We review the screening nature and many-body physics foundation of the van der Waals density functional (vdW-DF) method [Berland Ket al2015Rep. Prog. Phys.78066501], a systematic approach to construct truly nonlocal exchange-correlation energy density functionals. To that end we define and focus on a class of consistent vdW-DF versions that adhere to the Lindhard screening logic of the full method formulation. The consistent-exchange vdW-DF-cx version [Berland K and Hyldgaard P 2014Phys. Rev. B89035412] and its spin extension [Thonhauser Tet al2015Phys. Rev. Lett.115136402] represent the first examples of this class; in general, consistent vdW-DFs reflect a concerted expansion of a formal recast of the adiabatic-connection formula [Hyldgaard Pet al2014Phys. Rev. B90075148], an exponential summation of contributions to the local-field response, and the Dyson equation. We argue that the screening emphasis is essential because the exchange-correlation energy reflects an effective electrodynamics set by a long-range interaction. Two consequences are that (1) there are, in principle, no wiggle room in how one balances exchange and correlation, for example, in vdW-DF-cx, and that (2) consistent vdW-DFs have a formal structure that allows them to incorporate vertex-correction effects, at least in the case of levels that experience recoil-less interactions (for example, near the Fermi surface). We explore the extent to which the strictly nonempirical vdW-DF-cx formulation can serve as a systematic extension of the constraint-based semilocal functionals. For validation, we provide a complete survey of vdW-DF-cx performance for broad molecular processes, for the full set of 55 benchmarks in GMTKN55 [Goerigk Let al2017Phys. Chem. Chem. Phys.1932184] and comparing to the quantum-chemistry calculations that are summarized in that paper. We also provide new vdW-DF-cx results for metal surface energies and work functions that we compare to experiment. Finally, we use the screening insight to separate the vdW-DF nonlocal-correlation term into pure-vdW-interaction and local-field-susceptibility effects and present tools to compute and map the binding signatures.

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.

van der Waals forces in density functional theory: a review of the vdW-DF method

A density functional theory (DFT) that accounts for van der Waals (vdW) interactions in condensed matter, materials physics, chemistry, and biology is reviewed. The insights that led to the construction of the Rutgers-Chalmers van der Waals density functional (vdW-DF) are presented with the aim of giving a historical perspective, while also emphasizing more recent efforts which have sought to improve its accuracy. In addition to technical details, we discuss a range of recent applications that illustrate the necessity of including dispersion interactions in DFT. This review highlights the value of the vdW-DF method as a general-purpose method, not only for dispersion bound systems, but also in densely packed systems where these types of interactions are traditionally thought to be negligible.

Desorption of H2, HD, and D2 through temporary negative ion formation

Physically adsorbed H(2), HD, and D(2) desorbs via the low-energy 2Σ(u)+ electron scattering resonance with characteristic electron energy dependent cross sections as we show by desorption and electron energy-loss measurements. Only H(2) and HD desorb, for energetic reasons, via excitation of the j = 0 → 2 rotational mode and its subsequent decay into free translational motion, while all three isotopes desorb by direct transitions to unbound continuum states, with a probability close to unity and an unexpectedly weak mass dependence ~m(-1/2).

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.

Electron-induced desorption of physisorbed H2 via resonance vibrational excitation

We show by electron energy-loss measurements that desorption of physisorbed H(2) and D(2) induced by low-energy electrons takes place with large cross sections, predominantly via resonance excitation of the molecule-surface vibrational mode. The observed H(2), D(2) cross-section ratio supports a picture where rotation-translation conversion of the resonance excited j = 0 --> 2 rotational transition contributes to the desorption of H(2), while this channel is energetically closed for D(2).

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.

Selective adsorption, bound states, and potential parameters for He, Ne, and Ar interacting with a Cu(110) surface

Using nozzle beams of He, Ne, and Ar, we have measured diffractive selective adsorption resonances from a Cu(110) surface kept at 20 K. Bound state energies of the atom-surface potentials have been determined from plots of the measured resonance energies versus incident angle and their fits to calculated kinematical dispersion relations. For 3He and 4He we have found a unique level assignment that is compatible with a single gas-surface potential curve with a well depth of 6.05 meV of the He-Cu(110) potential. This value is about 10% larger than the prediction of 5.55 meV from the current physisorption theory. The Ne and Ar data reveal a large number of closely spaced levels with level separations and estimated van der Waals coefficients that are compatible with available theoretical data.