The role of van der Waals forces in water adsorption on metals

Hybrid Quantum Mechanical, Molecular Mechanical, and Machine Learning Potential for Computing Aqueous-Phase Adsorption Free Energies on Metal Surfaces

Performing reliable computer simulations of elementary processes occurring at metal-water interfaces is pivotal for novel catalyst design in sustainable energy applications. Computational catalyst design hinges on the ability to reliably and efficiently compute the potential energy surface (PES) of the system. Due to the large system sizes needed for studying processes at liquid water-metal interfaces, these systems can currently not be described using density functional theory (DFT). In this work, we used a hybrid quantum mechanical, molecular mechanical, and machine learning potential for studying the adsorption behavior of phenol, atomic hydrogen, 2-butanol, and 2-butanone on the (0001) facet of Ru under reducing conditions when Ru is not oxidized. Specifically, we describe the adsorbate and the surrounding metal atoms at the DFT level of theory. Here, we also considered the electrostatic field effect of the water molecules on adsorbate-metal interactions. Next, for the water-water and water-adsorbate interactions, we used established classical force fields. Finally, for the water-Ru surface interaction, for which no reliable force fields have been published, we used Behler-Parrinello high-dimensional neural network potentials (HDNNPs). Employing this setup, we used our explicit solvation for metal surface (eSMS) approach to compute the aqueous-phase effect on the low-coverage adsorption of selected molecules and atoms on the (0001) facet of Ru. In agreement with previous experimental and computational studies of oxygenated molecules over transition metal facets, we found that liquid water destabilizes the tested adsorbates on Ru(0001). Interestingly, our findings indicate that adsorbates on Ru are less affected by the presence of an aqueous phase than on other transition metals (e.g., Pt), highlighting the necessity of experimental investigations of Ru-based catalytic systems in liquid water.

Confinement Effects of Two-Dimensional Surfaces on Water Adsorption and Dissociation over Pt(111)

It has been established that the confined space created by stacking a two dimensional (2D) surface atop a metal catalyst serves as a nano-reactor. According to recent research, when a graphene (Gr) overlayer encloses a catalyst from above, the activation barrier for the water dissociation reaction, a process with major industrial significance, decreases. In order to investigate how the effect of confinement varies among different two-dimensional (2D) materials, we study the adsorption and dissociation barriers of water molecule on (111) under graphene, hexagonal boron nitride (h-BN), and heptazine-based graphitic carbon nitride (g-C3N4) layers using density functional theory calculations. Our findings reveal that the strength of adsorption does not decrease consistently with a reduction in the height of the 2D overlayer. Furthermore, a smaller barrier is not always the consequence of poorer adsorption of the reactant. We also examine the effect of confinement on the shape of the reaction path, on the frequencies of vibrational modes, and on the rate constants derived using the harmonic transition state theory. Overall, all three 2D surfaces cause a decrease in barrier height and a weakening of adsorption, though to differing degrees due to a mix of mechanical, geometric and electronic variables.

Benchmarking water adsorption on metal surfaces with ab initio molecular dynamics

Solid-water interfaces are ubiquitous in nature and technology. In particular, technologies evolving in the green transition, such as electrocatalysis, heavily rely on the junction of an electrolyte and an electrode as a central part of the device. For the understanding of atomic-scale processes taking place at the electrolyte-electrode interface, density functional theory (DFT) has become the de facto standard. The validation of DFT's ability to simulate the interfacial solid/water interaction is crucial, and ideal simulation setups need to be identified in order to prevent avoidable systematic errors. Here, we develop a rigorous sampling protocol for benchmarking the adsorption/desorption energetics of water on metallic surfaces against experimental temperature programmed desorption, single crystal adsorption calorimetry, and thermal energy atom scattering. We screened DFT's quality on a series of transition metal surfaces, applying three of the most common exchange-correlation approximations: PBE-D3, RPBE-D3, and BEEF-vdW. We find that all three xc-functionals reflect the pseudo-zeroth order desorption of water rooted in the combination of attractive adsorbate-adsorbate interactions and their saturation at low and intermediate coverages, respectively. However, both RPBE-D3 and BEEF-vdW lead to more accurate water adsorption strengths, while PBE-D3 clearly overbinds near-surface water. We relate the variations in binding strength to specific variations in water-metal and water-water interactions, highlighting the structural consequences inherent in an uninformed choice of simulation parameters. Our study gives atomistic insight into water's complex adsorption equilibrium. Furthermore, it represents a guideline for future DFT-based simulations of solvated solid interfaces by providing an assessment of systematic errors in specific setups.

Experimental Study on the Preparation of Lignin-Based Activated Carbon and the Adsorption Performance for Phenol

Biomass waste and wastewater are important wastes in the process of industrialization, which need to be effectively treated and utilized. In this work, an innovative method of collaborative treatment of biomass waste and phenol-containing wastewater is proposed. Biomass waste was used to produce activated carbon (AC), and then AC was used for phenol removal in wastewater treatment. Two kinds of typical biomass waste material, namely, coconut shell and lignin, were used. Physical activation (steam activation) and chemical activation methods were compared. Results show that steam activation is an effective method for coconut shell AC production. The largest Brunauer-Emmett-Teller (BET) surface area was 1065 m2/g at 800 °C. Chemical activation could produce AC samples with higher BET specific surface area. The lignin AC with K2CO3 activation has the largest BET surface of 1723.8 m2/g at 800 °C. FTIR results indicated that K2CO3 activation could greatly enhance the formation of surface oxygen-containing functional groups. Both coconut shell AC and lignin AC samples show excellent performance for phenol removal. The highest phenol removal efficiency for coconut shell AC and lignin AC are 96.87% and 98.22%, respectively. Adsorption kinetic analysis show that the pseudo-first-order kinetic model is able to describe the adsorption characteristics of phenol in wastewater treatment. Recycling properties show that regeneration of lignin AC could maintain high adsorption performance for phenol.

Accessing the usefulness of atomic adsorption configurations in predicting the adsorption properties of molecules with machine learning

We present a systematic study into the effect of adding atomic adsorption configurations into the training and validation dataset for a neural network's predictions of the adsorption energies of small molecules on single metal and bimetallic, single crystal surfaces. Specifically, we examine the efficacy of models trained with and without H and X atomic adsorption configurations, where X is C, N, or O, to predict XHn adsorption energies. In addition, we compare our machine learning models to traditional simple scaling relationships. We find that models trained with the atomic adsorption configurations outperform models trained with only molecular adsorption configurations, with as much as a 0.37 eV decrease in the MAE. We find that models trained with the atomic adsorption configurations slightly outperform traditional scaling relationships. In general, these results suggest it may be possible to vastly reduce the number of adsorption configurations one needs for training and validation datasets by supplementing said data with the adsorption configurations of composite atoms or smaller molecular fragments.

Interaction of water and carbon monoxide with MnO(001) thin films on Au(111)

Atomically defined MnO(001) thin films were grown on an Au(111) substrate, and their interaction with water (D2O) was investigated by infrared reflection absorption spectroscopy (IRAS) and thermal desorption spectroscopy (TDS). Carbon monoxide adsorption experiments were performed to probe surface atoms and defects on oxide films. Next, water interaction was investigated from which an associative binding pathway and a dissociative binding pathway were revealed, where the water molecules adsorb at terraces and water dissociation takes place at oxygen vacancies mediated by nearby Mn2+ sites. The IRAS data are supported by TDS experiments, which also manifest the importance of defects in the adsorption characteristics of MnO.

Stay Hydrated! Impact of Solvation Phenomena on the CO2 Reduction Reaction at Pb(100) and Ag(100) surfaces

Herein, a comprehensive computational study of the impact of solvation on the reduction reaction of CO2 to formic acid (HCOOH) and carbon monoxide on Pb(100) and Ag(100) surfaces is presented. Results further the understanding of how solvation phenomena influence the adsorption energies of reaction intermediates. We applied an explicit solvation scheme harnessing a combined density functional theory (DFT)/microkinetic modeling approach for the CO2 reduction reaction. This approach reveals high selectivities for CO formation at Ag and HCOOH formation on Pb, resolving the prior disparity between ab initio calculations and experimental observations. Furthermore, the detailed analysis of adsorption energies of relevant reaction intermediates shows that the total number of hydrogen bonds formed by HCOO plays a primary role for the adsorption strength of intermediates and the electrocatalytic activity. Results emphasize the importance of explicit solvation for adsorption and electrochemical reaction phenomena on metal surfaces.

3He spin-echo scattering indicates hindered diffusion of isolated water molecules on graphene-covered Ir(111)

The dynamics of water diffusion on carbon surfaces are of interest in fields as diverse as furthering the use of graphene as an industrial-coating technology and understanding the catalytic role of carbon-based dust grains in the interstellar medium. The early stages of water-ice growth and the mobility of water adsorbates are inherently dependent on the microscopic mechanisms that facilitate water diffusion. Here, we use 3He spin-echo quasi-inelastic scattering to probe the microscopic mechanisms responsible for the diffusion of isolated water molecules on graphene-covered and bare Ir(111). The scattering of He atoms provides a non-invasive and highly surface-sensitive means to measure the rate at which absorbates move around on a substrate at very low coverage. Our results provide an approximate upper limit on the diffusion coefficient for water molecules on GrIr(111) of < 10 - 12  m2/s, an order of magnitude lower than the coefficient that describes the diffusion of water molecules on the bare Ir(111) surface. We attribute the hindered diffusion of water molecules on the GrIr(111) surface to water trapping at specific areas of the corrugated moiré superstructure. Lower mobility of water molecules on a surface is expected to lead to a lower ice nucleation rate and may enhance the macroscopic anti-icing properties of a surface.

Preferred surface orientation for CO oxidation on SnO2 surfaces

In the present study, we perform a comparative study on the oxidation mechanism of CO gas molecules on SnO2 (110), (101), and (100) surfaces. The optimized adsorption configurations show that the adsorption of CO molecules could occur similarly on the three SnO2 surfaces via two adsorption modes, physisorption of CO on the Sn5c site that is considered as the first step for CO oxidation, followed by CO chemisorption on the O2c site resulting in the formation of CO2 species. Based on the calculated adsorption energies and CO molecule diffusion on SnO2 surfaces, CO molecule adsorption on the (101) surface exhibits the highest adsorption energy and the lowest reaction barrier for CO oxidation compared to the widely considered (110) surface or the (100) surface. These findings are expected to have a major impact on improving sensing properties toward toxic gas by means of surface-orientation engineering.

Orientational Water Bonding on Pt(111): Beyond the Frontier Orbital Principle

For decades, our understanding of water-metal bonding has been dominated by the frontier orbital principle in which globally stable water-metal interactions are ruled by HOMO interacting with metal surfaces. Using density functional theory calculations, herein, we have revealed that the frontier orbital principle cannot be applied to metastable water bonding on Pt(111), where the decisive role of HOMO is replaced by HOMO-1 in terms of the greatest orbital shifts and depopulations as the two different bonding indicators. Unlike the stable water configuration in which both HOMO-1 and HOMO prefer to overlap with metal states through σ-like orbital interactions, metastable configurations exhibit delicate competition or balance between σ-like and π-like orbital interactions exerted by HOMO-1 and HOMO, respectively. These findings have significantly deepened our understanding of orbital roles in water-metal bonding interactions and bridged the gap between theoretical understanding of electrified waters at electrochemical interfaces and water science on metal surfaces.

Ab initio study of water dissociation on a charged Pd(111) surface

The interactions between molecules and electrode surfaces play a key role in electrochemical processes and are a subject of extensive research, both experimental and theoretical. In this paper, we address the water dissociation reaction on a Pd(111) electrode surface, modeled as a slab embedded in an external electric field. We aim at unraveling the relationship between surface charge and zero-point energy in aiding or hindering this reaction. We calculate the energy barriers with dispersion-corrected density-functional theory and an efficient parallel implementation of the nudged-elastic-band method. We show that the lowest dissociation barrier and consequently the highest reaction rate take place when the field reaches a strength where two different geometries of the water molecule in the reactant state are equally stable. The zero-point energy contributions to this reaction, on the other hand, remain nearly constant across a wide range of electric field strengths, despite significant changes in the reactant state. Interestingly, we show that the application of electric fields that induce a negative charge on the surface can make nuclear tunneling more significant for these reactions.

The interaction mechanism of cesium with water on the SrTiO3(100) surface at room temperature

The interaction of water with cesium on the strontium titanate surface SrTiO3(100), was studied, mainly by means of work function measurements and thermal desorption spectroscopy. The catalytic role of cesium with respect to the dissociation of water on surface was investigated, by applying two different adsorption processes at room temperature (RT): (1) The adsorption of water on the cesium covered surface (sequential adsorption), and (2) the co-adsorption process (simultaneous adsorption) on surface. Based on the results and by adopting the Lewis acid–base model, we conclude that during the sequential adsorption the water molecules are mostly adsorbs non-dissociatively on surface, without oxidizing the alkaline overlayer. This seems to be due, first to the strong interaction between the alkaline adatoms and the substrate, and secondly to the limited maximum pre-deposited amount of cesium (≤ 0.45 ML). Instead, water dissociation appears to merely occur on defective sites of the substrate in accordance with previous studies. For a full cesium layer covered surface, the adsorbed water retracts the metallicity of cesium due to electrostatic interactions. In contrast to the sequential adsorption, during the co-adsorption process the oxidation of cesium takes place above a critical coverage of cesium (≥ 0.45 ML). It appears that the co-adsorbed cesium with water modifies the surface potential providing an effective template for cesium oxide, Cs2O development. Based on that, we suggest a catalytic reaction of water dissociation according to the Langmuir–Hinshelwood mechanism. Finally, we propose atomistic adsorption models for both processes of cesium with water adsorption.

Structure and interactions at the Mg(0001)/water interface: An ab initio study

A molecular level understanding of metal/bulk water interface structure is key for a wide range of processes, including aqueous corrosion, which is our focus, but their buried nature makes experimental investigation difficult and we must mainly rely on simulations. We investigate the Mg(0001)/water interface using second generation Car–Parrinello molecular dynamics (MD) to gain structural information, combined with static density functional theory calculations to probe the atomic interactions and electronic structure (e.g., calculating the potential of zero charge). By performing detailed structural analyses of both metal–surface atoms and the near-surface water, we find that, among other insights: (i) water adsorption causes significant surface roughening (the planar distribution for top-layer Mg has two peaks separated by [Formula: see text]), (ii) strongly adsorbed water covers only [Formula: see text] of available surface sites, and (iii) adsorbed water avoids clustering on the surface. Static calculations are used to gain a deeper understanding of the structuring observed in MD. For example, we use an energy decomposition analysis combined with calculated atomic charges to show that adsorbate clustering is unfavorable due to Coulombic repulsion between adsorption site surface atoms. Results are discussed in the context of previous simulations carried out on other metal/water interfaces. The largest differences for the Mg(0001)/water system appear to be the high degree of surface distortion and the minimal difference between the metal work function and metal/water potential of zero charge (at least compared to other interfaces with similar metal–water interaction strengths). The structural information, in this paper, is important for understanding aqueous Mg corrosion, as the Mg(0001)/water interface is the starting point for key reactions. Furthermore, our focus on understanding the driving forces behind this structuring leads to important insights for general metal/water interfaces.

Interactions between glucosides of the tip of the S1 subunit of SARS-CoV-2 spike protein and dry and wet surfaces of CuO and Cu-A model for the surfaces of coinage metals

Despite their importance there is little knowledge at the atomic scale on the interactions between fragments of SARS-CoV-2 and inorganic materials. Such knowledge is important to understand the survival of the virus at surfaces and for the development of antiviral materials. Here is reported a study of the interactions between glucoside monomers of the tip of the S1 subunit of SARS-CoV-2 spike protein with dry and wet surfaces of CuO and Cu, performed with dispersion corrected density functional theory—DFT. The three glucoside monomers that constitute the tip of S1: 6VSB, 6VXX and 6X6P, were adsorbed onto dry and wet CuO(111) and Cu(110) with different orientations and surface alignments.

There are large differences—of up to 1.3 eV—in binding energies between these monomers and the surfaces. These differences depend on: the type of surface; if the surface is wet or dry; if the glucosidic O-atom points towards or away from the surfaces; and to a smaller extent on the surface alignment of the monomers. All monomers bind strongly to the surfaces via molecular adsorption that does not involve bond breaking in the monomers at this stage. 6VSB has the larger adsorption energies—that reach 2.2 eV—due to its larger dipole moment. Both materials bind the monomers more strongly when their surfaces are dry. At Cu(110) the bonds are on average 1 eV stronger when the surface is dry when compared to wet. The difference between dry and wet CuO(111) is smaller, in the order of 0.2 eV. Overall, it is here shown that the stability of the monomers of the tip of the spike protein of the virus is very different at different surfaces. For a given surface the larger binding energies in dry conditions could explain the differences in the surface stability of the spike protein depending on the presence of moisture.

Ab Initio Simulations of Water/Metal Interfaces

Structures and processes at water/metal interfaces play an important technological role in electrochemical energy conversion and storage, photoconversion, sensors, and corrosion, just to name a few. However, they are also of fundamental significance as a model system for the study of solid-liquid interfaces, which requires combining concepts from the chemistry and physics of crystalline materials and liquids. Particularly interesting is the fact that the water-water and water-metal interactions are of similar strength so that the structures at water/metal interfaces result from a competition between these comparable interactions. Because water is a polar molecule and water and metal surfaces are both polarizable, explicit consideration of the electronic degrees of freedom at water/metal interfaces is mandatory. In principle, ab initio molecular dynamics simulations are thus the method of choice to model water/metal interfaces, but they are computationally still rather demanding. Here, ab initio simulations of water/metal interfaces will be reviewed, starting from static systems such as the adsorption of single water molecules, water clusters, and icelike layers, followed by the properties of liquid water layers at metal surfaces. Technical issues such as the appropriate first-principles description of the water-water and water-metal interactions will be discussed, and electrochemical aspects will be addressed. Finally, more approximate but numerically less demanding approaches to treat water at metal surfaces from first-principles will be briefly discussed.

Static and dynamic water structures at interfaces: A case study with focus on Pt(111)

An accurate atomistic treatment of aqueous solid-liquid interfaces necessitates the explicit description of interfacial water ideally via ab initio molecular dynamics simulations. Many applications, however, still rely on static interfacial water models, e.g., for the computation of (electro)chemical reaction barriers and focus on a single, prototypical structure. In this work, we systematically study the relation between density functional theory-derived static and dynamic interfacial water models with specific focus on the water-Pt(111) interface. We first introduce a general construction protocol for static 2D water layers on any substrate, which we apply to the low index surfaces of Pt. Subsequently, we compare these with structures from a broad selection of reference works based on the Smooth Overlap of Atomic Positions descriptor. The analysis reveals some structural overlap between static and dynamic water ensembles; however, static structures tend to overemphasize the in-plane hydrogen bonding network. This feature is especially pronounced for the widely used low-temperature hexagonal ice-like structure. In addition, a complex relation between structure, work function, and adsorption energy is observed, which suggests that the concentration on single, static water models might introduce systematic biases that are likely reduced by averaging over consistently created structural ensembles, as introduced here.

Dynamics and Surface Propensity of H+ and OH- within Rigid Interfacial Water: Implications for Electrocatalysis

Facile solvent reorganization promoting ion transfer across the solid-liquid interface is considered a prerequisite for efficient electrocatalysis. We provide first-principles insight into this notion by examining water self-ion dynamics at a highly rigid NaCl(100)-water interface. Through extensive density functional theory molecular dynamics simulations, we demonstrate for both acidic and alkaline solutions that Grotthuss dynamics is not impeded by a rigid water structure. Conversely, decreased proton transfer barriers and a striking propensity of H₃O⁺ and OH^- for stationary interfacial water are found. Differences in the ideal hydration structure of the ions, however, distinguish their behavior at the water contact layer. While hydronium can maintain its optimal solvation, the preferentially hypercoordinated hydroxide is repelled from the immediate vicinity of the surface due to interfacial coordination reduction. This has implications for alkaline hydrogen electrosorption in which the formation of undercoordinated OH^- at the surface is proposed to contribute to the observed sluggish kinetics.

Simulation of Metal-Supported Metal-Nanoislands: A Comparison of DFT Methods

We have evaluated various density functional theory (DFT) methods to simulate geometric, energetic, electronic, and hydrogen adsorption properties of metal-nanoparticles supported on metal surfaces. We used Pt and Pd nanoislands on Au(111) as model systems. The evaluated DFT methods include GGA (PW91, PBE, RPBE, revPBE, and PBESol), GGA with van der Waals (vdW) corrected (PBE-D3), GGA with optimized vdW functionals (revPBE-vdW), meta-GGA (SCAN and MS2), and the machine learning-based method BEEF-vdW. The results show that the various DFT methods yield similar geometric and electronic properties for Pt (or Pd) nanoislands on Au(111). The DFT methods also produce similar relative energetics for small Pt (or Pd) clusters with different conformations on Au(111). The results show that a triatomic cluster of Pt on Au(111) is more stable with a linear conformation. In contrast, a triatomic cluster of Pd is more stable with a triangular conformation. For clusters with four or more atoms, Pt and Pd clusters on Au(111) prefer non-linear conformation. We found that the various DFT methods yield different results only for the adsorption energy of hydrogen.

Effects of applied voltage on water at a gold electrode interface from ab initio molecular dynamics

Electrode–water interfaces under voltage bias demonstrate anomalous electrostatic and structural properties that are influential in their catalytic and technological applications. Mean-field and empirical models of the electrical double layer (EDL) that forms in response to an applied potential do not capture the heterogeneity that polarizable, liquid-phase water molecules engender. To illustrate the inhomogeneous nature of the electrochemical interface, Born–Oppenheimer ab initio molecular dynamics calculations of electrified Au(111) slabs interfaced with liquid water were performed using a combined explicit–implicit solvent approach. The excess charges localized on the model electrode were held constant and the electrode potentials were computed at frequent simulation times. The electrode potential in each trajectory fluctuated with changes in the atomic structure, and the trajectory-averaged potentials converged and yielded a physically reasonable differential capacitance for the system. The effects of the average applied voltages, both positive and negative, on the structural, hydrogen bonding, dynamical, and vibrational properties of water were characterized and compared to literature where applicable. Controlled-potential simulations of the interfacial solvent dynamics provide a framework for further investigation of more complex or reactive species in the EDL and broadly for understanding electrochemical interfaces in situ.

Ab initio molecular dynamics of an aqueous electrode interface reveal the electrostatic, structural, and dynamic effects of quantifiable voltage biases on water.

A first-principles analysis of the charge transfer in magnesium corrosion

Magnesium is the lightest structural engineering material and bears high potential to manufacture automotive components, medical implants and energy storage systems. However, the practical use of untreated magnesium alloys is restricted as they are prone to corrosion. An essential prerequisite for the control or prevention of the degradation process is a deeper understanding of the underlying corrosion mechanisms. Prior investigations of the formation of gaseous hydrogen during the corrosion of magnesium indicated that the predominant mechanism for this process follows the Volmer-Heyrovský rather than the previously assumed Volmer-Tafel pathway. However, the energetic and electronic states of both reaction paths as well as the charge state of dissolved magnesium have not been fully unraveled yet. In this study, density functional theory calculations were employed to determine these parameters for the Volmer, Tafel and Heyrovský steps to gain a comprehensive understanding of the major corrosion mechanisms responsible for the degradation of magnesium.

Conductance and configuration of molecular gold-water-gold junctions under electric fields

Water molecules can mediate charge transfer in biological and chemical reactions by forming electronic coupling pathways. Understanding the mechanism requires a molecular-level electrical characterization of water. Here, we describe the measurement of single water molecular conductance at room temperature, characterize the structure of water molecules using infrared spectroscopy, and perform theoretical studies to assist in the interpretation of the experimental data. The study reveals two distinct states of water, corresponding to a parallel and perpendicular orientation of the molecules. Water molecules switch from parallel to perpendicular orientations on applying an electric field, producing switching from high to low conductance states, thus enabling the determination of single water molecular dipole moments. The work further shows that water-water interactions affect the atomic scale configuration and conductance of water molecules. These findings demonstrate the importance of the discrete nature of water molecules in electron transfer and set limits on water-mediated electron transfer rates.

Ten Facets, One Force Field: The GAL19 Force Field for Water-Noble Metal Interfaces

Understanding the structure of the water/metal interfaces plays an important role in many areas ranging from surface chemistry to environmental processes. The size, required phase-space sampling, and the slow diffusion of molecules at the water/metal interfaces motivate the development of accurate force fields. We develop and parametrize GAL19, a novel force field, to describe the interaction of water with two facets (111 and 100) of five metals (Pt, Pd, Au, Ag, Cu). To increase transferability compared to its predecessor GAL17, the water-metal interaction is described as a sum of pairwise terms. The interaction energy has three contributions: (i) physisorption is described via a Tang and Toennies potential, (ii) chemisorption and surface corrugation rely on an attractive Gaussian term, and (iii) the angular dependence is explicitly included as a truncated Fourier series. Thirteen parameters are used for each metal surface and were fitted on 250 water adsorption energies computed at the PBE+dDsC level. The performance of GAL19 was evaluated on a set of more than 600 DFT adsorption energies for each surface, leading to an average root-mean-square deviation of only 1 kcal/mol, correctly reproducing the adsorption trends: strong on Pt and Pd but weaker on Ag, Au, and Cu. This force field was then used to simulate the water/metal interface for all ten surfaces for 1 ns. Structural analyses reveal similar tendencies for all surfaces: a first, dense water layer that is mostly adsorbed on the metal top sites and a second layer up to around 6 Å, which is less structured. On Pt and Pd, the first layer is strongly organized with water lying flat on the surface. The pairwise additive functional form allows one to simulate the water adsorption on alloys, which is demonstrated at the example of Ag/Cu and Au/Pt alloys. The water/Ag-Cu interface is predicted to be disordered with water mostly adsorbed on Cu which should exacerbate the Ag reactivity. On the contrary, incorporating Pt into Au materials leads to a structuring of the water interface. Our promising results make GAL19 an ideal candidate to get representative sampling of complex metal/water interfaces as a first step toward accurate estimation of free energies of reactions in solution at the metal interface.

Origins of fast diffusion of water dimers on surfaces

The diffusion of water molecules and clusters across the surfaces of materials is important to a wide range of processes. Interestingly, experiments have shown that on certain substrates, water dimers can diffuse more rapidly than water monomers. Whilst explanations for anomalously fast diffusion have been presented for specific systems, the general underlying physical principles are not yet established. We investigate this through a systematic ab initio study of water monomer and dimer diffusion on a range of surfaces. Calculations reveal different mechanisms for fast water dimer diffusion, which is found to be more widespread than previously anticipated. The key factors affecting diffusion are the balance of water-water versus water-surface bonding and the ease with which hydrogen-bond exchange can occur (either through a classical over-the-barrier process or through quantum-mechanical tunnelling). We anticipate that the insights gained will be useful for understanding future experiments on the diffusion and clustering of hydrogen-bonded adsorbates.

The experimental observation that water dimers diffuse more rapidly than monomers across materials’ surfaces is yet to be clarified. Here the authors show by ab initio calculations classical and quantum mechanical mechanisms for faster water dimer diffusion on a broad range of metal and non-metal surfaces.

Catalytic Effect of Hydrogen Bond on Oxhydryl Dehydrogenation in Methanol Steam Reforming on Ni(111)

Dehydrogenation of H₃COH and H₂O are key steps of methanol steam reforming on transition metal surfaces. Oxhydryl dehydrogenation reactions of H_xCOH (x = 0-3) and OH on Ni (111) were investigated by DFT calculations with the OptB88-vdW functional. The transition states were searched by the climbing image nudged elastic band method and the dimer method. The activation energies for the dehydrogenation of individual H_xCOH* are 68 to 91 kJ/mol, and reduced to 12-17 kJ/mol by neighboring OH*. Bader charge analysis showed the catalysis role of OH* can be attributed to the effect of hydrogen bond (H-bond) in maintaining the charge of oxhydryl H in the reaction path. The mechanism of H-bond catalysis was further demonstrated by the study of OH* and N* assisted dehydrogenation of OH*. Due to the universality of H-bond, the H-bond catalysis shown here, is of broad implication for studies of reaction kinetics.

Density Functional Theory Based Micro- and Macro-Kinetic Studies of Ni-Catalyzed Methanol Steam Reforming

The intrinsic mechanism of Ni-catalyzed methanol steam reforming (MSR) is examined by considering 54 elementary reaction steps involved in MSR over Ni(111). Density functional theory computations and transition state theory analyses are performed on the elementary reaction network. A microkinetic model is constructed by combining the quantum chemical results with a continuous stirring tank reactor model. MSR rates deduced from the microkinetic model agree with the available experimental data. The microkinetic model is used to identify the main reaction pathway, the rate determining step, and the coverages of surface species. An analytical expression of MSR rate is derived based on the dominant reaction pathway and the coverages of surface species. The analytical rate equation is easy to use and should be very helpful for the design and optimization of the operating conditions of MSR.

Atomic and molecular adsorption on single platinum atom at the graphene edge: A density functional theory study

We present a density functional theory study of atomic and molecular adsorption on a single Pt atom deposited at the edges of graphene. We investigate geometric and electronic structures of atoms (H, C, N, and O) and molecules (O₂, CO, OH, NO, H₂O, and OOH) on a variety of Pt deposited graphene edges and compare the adsorption states with those on a Pt(111) surface and on a Pt single atom. Furthermore, using the calculated adsorption energy and simple kinetic models, the catalytic activities of a Pt single-atom catalyst for the oxygen reduction reaction and CO oxidation are discussed.

DFT + U study of H2O adsorption and dissociation on stoichiometric and nonstoichiometric CuO(1 1 1) surfaces

Surface interaction through adsorption and dissociation between H₂O and metal oxides plays an important role in many industrial as well as fundamental processes. To gain further insights on the interaction, this study performs dispersion-corrected Hubbard-corrected density functional theory calculations in H₂O adsorption and dissociation on stoichiometric and nonstoichiometric CuO(1 1 1) surfaces. The nonstoichiometric surfaces consist of oxygen vacancy defect and oxygen-preadsorbed surfaces. This study finds that H₂O is chemically adsorbed on the top of Cu_sub and Cu_sub-Cu_sub bridge due to the interaction of its p  orbital with d orbital of Cu. The adsorption is found to be the strongest on the surface with the oxygen vacancy defect, followed by the stoichiometric surface, and the oxygen-preadsorbed surface. The oxygen vacancy increases the reactivity for H₂O adsorption and reduces the reaction energy required for H₂O dissociation on the surface. However, the surface modification by the oxygen-preadsorbed significantly reduces the barrier energy for H₂O dissociation when compared with the other surfaces.

Adsorption processes on a Pd monolayer-modified Pt(111) electrode

Specific adsorption of anions is an important aspect in surface electrochemistry for its influence on reaction kinetics in either a promoted or inhibited fashion. Perchloric acid is typically considered as an ideal electrolyte for investigating electrocatalytic reactions due to the lack of specific adsorption of the perchlorate anion on several metal electrodes. In this work, cyclic voltammetry and computational methods are combined to investigate the interfacial processes on a Pd monolayer deposited on Pt(111) single crystal electrode in perchloric acid solution. The “hydrogen region” of this Pd_MLPt(111) surface exhibits two voltammetric peaks: the first “hydrogen peak” at 0.246 V_RHE actually involves the replacement of hydrogen by hydroxyl, and the second “hydrogen peak” H_II at 0.306 V_RHE appears to be the replacement of adsorbed hydroxyl by specific perchlorate adsorption. The two peaks merge into a single peak when a more strongly adsorbed anion, such as sulfate, is involved. Our density functional theory calculations qualitatively support the peak assignment and show that anions generally bind more strongly to the Pd_MLPt(111) surface than to Pt(111).

Specific adsorption of anions is an important aspect in surface electrochemistry for its influence on reaction kinetics in either a promoted or inhibited fashion.

Including dispersion in density functional theory for adsorption on flat oxide surfaces, in metal-organic frameworks and in acidic zeolites

We examine the performance of nine commonly used methods for including dispersion interactions in density functional theory (DFT): three different parametrizations of damped 1/Rn terms (n = 6, 8, …) added to the DFT energy (Grimme's D2 and D3 parameterizations as well as that of Tkatchenko and Scheffler), three different implementations of the many-body dispersion approach (MBD, MBD/HI and MBD/FI), the density-dependent energy correction, called dDsC, and two "first generation" van der Waals density functionals, revPBE-vdW and optB86b-vdW. As test set we use eight molecule-surface systems for which agreement has been reached between experiment and hybrid QM:QM calculations within chemical accuracy limits (±4.2 kJ mol-1). It includes adsorption of carbon monoxide and dioxide in the Mg2(2,5-dioxido-1,4-benzenedicarboxylate) metal-organic framework (Mg-MOF-74, CPO-27-Mg), adsorption of carbon monoxide as well as of monolayers of methane and ethane on the MgO(001) surface, as well as adsorption of methane, ethane and propane in H-chabazite (H-CHA). D2 with Ne parameters for Mg2+, D2(Ne), MBD/HI and MBD/FI perform best. With the PBE functional, the mean unsigned errors are 6.1, 5.6 and 5.4 kJ mol-1, respectively.

Elementary reaction pathway study and a deduced macrokinetic model for the unified understanding of Ni-catalyzed steam methane reforming

DFT based microkinetics and macrokinetics that give quantitative explanations of the Ni-catalyzed steam methane reforming reactions.

Influence of Van der Waals Interactions on the Solvation Energies of Adsorbates at Pt-Based Electrocatalysts

Solvation can significantly modify the adsorption energy of species at surfaces, thereby influencing the performance of electrocatalysts and liquid-phase catalysts. Thus, it is important to understand adsorbate solvation at the nanoscale. Here we evaluate the effect of van der Waals (vdW) interactions described by different approaches on the solvation energy of *OH adsorbed on near-surface alloys (NSAs) of Pt. Our results show that the studied functionals can be divided into two groups, each with rather similar average *OH solvation energies: (1) PBE and PW91; and (2) vdW functionals, RPBE, PBE-D3 and RPBE-D3. On average, *OH solvation energies are less negative by ∼0.14 eV in group (2) compared to (1), and the values for a given alloy can be extrapolated from one functional to another within the same group. Depending on the desired level of accuracy, these concrete observations and our tabulated values can be used to rapidly incorporate solvation into models for electrocatalysis and liquid-phase catalysis.

Interaction of H2O with the Platinum Pt (001), (011), and (111) Surfaces: A Density Functional Theory Study with Long-Range Dispersion Corrections

Platinum is a noble metal that is widely used for the electrocatalytic production of hydrogen, but the surface reactivity of platinum toward water is not yet fully understood, even though the effect of water adsorption on the surface free energy of Pt is important in the interpretation of the morphology and catalytic properties of this metal. In this study, we have carried out density functional theory calculations with long-range dispersion corrections [DFT-D3-(BJ)] to investigate the interaction of H₂O with the Pt (001), (011), and (111) surfaces. During the adsorption of a single H₂O molecule on various Pt surfaces, it was found that the lowest adsorption energy (E _ads) was obtained for the dissociative adsorption of H₂O on the (001) surface, followed by the (011) and (111) surfaces. When the surface coverage was increased up to a monolayer, we noted an increase in E _ads/H₂O with increasing coverage for the (001) surface, while for the (011) and (111) surfaces, E _ads/H₂O decreased. Considering experimental conditions, we observed that the highest coverage was obtained on the (011) surface, followed by the (111) and (001) surfaces. However, with an increase in temperature, the surface coverage decreased on all the surfaces. Total desorption occurred at temperatures higher than 400 K for the (011) and (111) surfaces, but above 850 K for the (001) surface. From the morphology analysis of the Pt nanoparticle, we noted that, when the temperature increased, only the electrocatalytically active (111) surface remained.

Structure, Dynamics, and Wettability of Water at Metal Interfaces

The water/metal interface often governs important chemophysical processes in various technologies. Therefore, from scientific and engineering perspectives, the detailed molecular-level elucidation of the water/metal interface is of high priority, but the related research is limited. In experiments, the surface-science techniques, which can provide full structural details of the surface, are not easy to directly apply to the interfacial systems under ambient conditions, and the well-defined facets cannot be entirely free from contamination at the contact with water. To answer long-standing debates regarding the wettability, structure, and dynamics of water at metal interfaces, we here develop reliable first-principles-based multiscale simulations. Using the state-of-the-art simulations, we find that the clean metal surfaces are actually superhydrophilic and yield zero contact angles. Furthermore, we disclose an inadequacy of widespread ice-like bilayer model of the water adlayers on metal surfaces from both averaged structural and dynamic points of view. Our findings on the nature of water on metal surfaces provide new molecular level perspectives on the tuning and design of water/metal interfaces that are at the heart of many energy applications.

Imaging and Dynamics of Water Hexamer Confined in Nanopores

Epitaxial two-dimensional (2D) nanostructures with regular patterns show great promise as templates for adsorbate confinement. Prospectively, employing 2D semiconductors with reduced density of states leads to a long excited-state lifetime that allows us to directly image the dynamics of the adsorbate. We show that epitaxial blue phosphorene (blueP) on Au(111) provides such a platform to trap water molecules in the periodic nanopores without formation of strong bonds. The trapped water aggregate is tentatively assigned to a hexamer based on our scanning tunneling microscopy studies and first-principles calculations. Real-space observation of conformational switching of the hexamer induced by inelastic electrons is achieved by using low-temperature scanning tunneling microscopy with molecular resolution. We found a localized interfacial charge rearrangement between the water hexamer and P atoms underneath that is responsible for the reversible desorption and adsorption of water molecules by changing the sample bias polarity from positive to negative, offering a promising strategy for engineering the electronic properties of blueP.

A Water Solvation Shell Can Transform Gold Metastable Nanoparticles in the Fluxional Regime

Solvated gold nanoparticles have been modeled in the fluxional regime by density functional theory including dispersion forces for an extensive set of conventional morphologies. The study of isolated adsorption of one water molecule shows that the most stable adsorption forms are similar (corners and edges) regardless of the nanoparticle shape and size, although the adsorption strength differs significantly (0.15 eV). When a complete and explicit water solvation shell interacts with gold nanoclusters, metastable in vacuum and presenting a predominance of (100) square facets (ino-decahedra Au₅₅ and Au₁₄₇), these nanoparticles are found unstable and transform into the closest morphologies exhibiting mainly (111) triangular facets and symmetries. The corresponding adsorption strength per water molecule becomes independent of shape and size and is enhanced by the formation of two hydrogen bonds on average. For applications in radiotherapy, this study suggests that the shapes of small gold nanoparticles should be homogenized by interacting with the biological environment.

Interconversion of hydrated protons at the interface between liquid water and platinum

The free energy barriers for hydrogen transfer at the H₂O/Pt(111) interface calculated usingab initiomolecular dynamics and umbrella sampling.

Adsorption of Amino Acids on Gold: Assessing the Accuracy of the GolP-CHARMM Force Field and Parametrization of Au-S Bonds

The interaction of amino acids with metal electrodes plays a crucial role in bioelectrochemistry and the emerging field of bionanoelectronics. Here we present benchmark calculations of the adsorption structure and energy of all natural amino acids on Au(111) in vacuum using a van-der-Waals density functional (revPBE-vdW) that showed good performance on the S22 set of weakly bound dimers (mean relative unsigned error (MRUE) wrt CCSD(T)/CBS = 13.3%) and adsorption energies of small organic molecules on Au(111) (MRUE wrt experiment = 11.2%). The vdW-DF results are then used to assess the accuracy of a popular force field for Au-amino acid interactions, GolP-CHARMM, which explicitly describes image charge interactions via rigid-rod dipoles. We find that while the force field underestimates adsorption distances, it does reproduce the binding energy rather well (MRUE wrt revPBE-vdW = 11.3%) with the MRUE decreasing in the order Cys, Met > amines > aliphatic > carboxylic > aromatic. We also present a parametrization of the bonding interaction between sulfur-containing molecules and the Au(111) surface and report force field parameters that are compatible with GolP-CHARMM. We believe the vdW-DF calculations presented herein will provide useful reference data for further force field development, and that the new Au-S bonding parameters will enable improved simulations of proteins immobilized on Au-electrodes via S-linkages.

Competing adsorption mechanisms of pyridine on Cu, Ag, Au, and Pt(110) surfaces

We explore the adsorption of pyridine on Cu, Ag, Au, and Pt(110) surfaces using density functional theory. To account for the van der Waals interaction, we use the optB86b-vdW, optB88-vdW, optPBE-vdW, revPBE-vdW, and rPW86-vdW2 functionals. For comparison, we also run calculations using the generalized gradient approximation-PBE (Perdew-Burke-Ernzerhof) functional. We find the most stable adsorption site to depend on both metal and functional, with two energetically favorable adsorption sites, namely, a vertically oriented site and a flat pyridine site. We calculate that every functional predicts pyridine to lie in the vertical configuration on the coinage metals at a low coverage. On Pt(110), by contrast, we calculate all the functionals-except rPW86-vdW2-to predict pyridine to lie flat at a low coverage. By analyzing these differences for these adsorption configurations, along with various geometric and electronic properties of the adsorbate/substrate system, we access in detail the performance of the 6 functionals we use. We also characterize the nature of the bonding of pyridine on the coinage metals from weak to strong physisorption, depending on the functional used. On Pt(110), we characterize the nature of the bonding of pyridine as ranging from strong physisorption to strong chemisorption depending again on the functional used, illustrating both the importance of the van der Waals interaction to this system and that this system can make a stringent test for computational methods.