Size dependence of hydrophobic hydration at electrified gold/water interfaces

The role of the water contact layer on hydration and transport at solid/liquid interfaces

Understanding the structure in the nanoscopic region of water that is in direct contact with solid surfaces, so-called contact layer, is key to quantifying macroscopic properties that are of interest to e.g. catalysis, ice nucleation, nanofluidics, gas adsorption, and sensing. We explore the structure of the water contact layer on various technologically relevant solid surfaces, namely graphene, MoS[Formula: see text], Au(111), Au(100), Pt(111), and Pt(100), which have been previously hampered by time and length scale limitations of ab initio approaches or force field inaccuracies, by means of molecular dynamics simulations based on ab initio machine learning potentials built using an active learning scheme. Our results reveal that the in-plane intermolecular correlations of the water contact layer vary greatly among different systems: Whereas the contact layer on graphene and on Au(111) is predominantly homogeneous and isotropic, it is inhomogeneous and anisotropic on MoS[Formula: see text], on Au(100), and on the Pt surfaces, where it additionally forms two distinct sublayers. We apply hydrodynamics and the theory of the hydrophobic effect, to relate the energy corrugation and the characteristic length-scales of the contact layer with wetting, slippage, the hydration of small hydrophobic solutes and diffusio-osmotic transport. Thus, this work provides a microscopic picture of the water contact layer and links it to macroscopic properties of liquid/solid interfaces that are measured experimentally and that are relevant to wetting, hydrophobic solvation, nanofluidics, and osmotic transport.

Slip Opacity and Fast Osmotic Transport of Hydrophobes at Aqueous Interfaces with Two-Dimensional Materials

We investigate the interfacial transport of water and hydrophobic solutes on van der Waals bilayers and heterostructures formed by stacking graphene, hBN, and MoS2 using extensive ab initio molecular dynamics simulations. We compute water slippage and the diffusio-osmotic transport coefficient of hydrophobic particles at the interface by combining hydrodynamics and the theory of the hydrophobic effect. We find that slippage is dominated by the layer that is in direct contact with water and only marginally altered by the second layer, leading to a so-called "slip opacity". The screening of the lateral forces, where the liquid does not feel the forces coming from the second nearest layer, is one of the factors leading to the "slip opacity" in our systems. The diffusio-osmotic transport of small hydrophobes (with a radius below 2.5 Å) is also affected by the slip opacity, being dramatically enhanced by slippage. Furthermore, the direction of diffusio-osmotic flow is controlled by the solute size, with the flow in the opposite direction of the concentration gradient for smaller hydrophobes, and vice versa for larger ones. We connect our findings to the wetting properties of two-dimensional materials, and we propose that slippage and wetting can be controlled separately: whereas the slippage is mostly determined by the layer in closer proximity to water, wetting can be finely tuned by stacking different two-dimensional materials. Our study advances the computational design of two-dimensional materials and van der Waals heterostructures, enabling precise control over wetting and slippage properties for applications in coatings and water purification membranes.

Tuning hydrogen bond network connectivity in the electric double layer with cations

Hydrogen bond (H-bond) network connectivity in electric double layers (EDLs) is of paramount importance for interfacial HER/HOR electrocatalytic processes. However, it remains unclear whether the cation-specific effect on H-bond network connectivity in EDLs exists. Herein, we report simulation evidence from ab initio molecular dynamics that cations at Pt(111)/water interfaces can tune the structure and the connectivity of H-bond networks in EDLs. As the surface charge density σ becomes more negative, we show that the connectivity of the H-bond networks in EDLs of the Na+ and Ca2+ systems decreases markedly; in stark contrast, the connectivity of the H-bond networks in EDLs of the Mg2+ system increases slightly. Further analysis revealed that the interplay between the hydration of cations and the interfacial water structure plays a key role in the connectivity of H-bond networks in EDLs. These findings highlight the key roles of cations in EDLs and electrocatalysis.

Tuning Acid-Base Chemistry at an Electrified Gold/Water Interface

Acid-base reactions are ubiquitous in solution chemistry, as well as in electrochemistry. However, macroscopic concepts derived in solutions, such as pKa and pH, differ significantly at electrified metal-aqueous interfaces due to specific solvation and applied voltage. Here, we measure the pKa values of an amino acid, glycine, at a gold/water interface under a varying applied voltage by means of spectroscopic titration. With the help of simulations, we propose a general model to understand potential-dependent shifts in pKa values in terms of local hydrophobicity and electric fields. These parameters can be tuned by adjusting the metal surface and applied voltage, respectively, offering promising, but still unexplored, paths to regulate reactivity. Our results change the focus with respect to common interpretations based on, for example, apparent local pH effects and open interesting perspectives for electrochemical reaction steering.

S.O.S: Shape, orientation, and size tune solvation in electrocatalysis

Current models to understand the reactivity of metal/aqueous interfaces in electrochemistry, e.g., volcano plots, are based on the adsorption free energies of reactants and products, which are often small hydrophobic molecules (such as in CO2 and N2 reduction). Calculations played a major role in the quantification and comprehension of these free energies in terms of the interactions that the reactive species form with the surface. However, solvation free energies also come into play in two ways: (i) by modulating the adsorption free energy together with solute-surface interactions, as the solute has to penetrate the water adlayer in contact with the surface and get partially desolvated (which costs free energy); (ii) by regulating transport across the interface, i.e., the free energy profile from the bulk to the interface, which is strongly non-monotonic due to the unique nature of metal/aqueous interfaces. Here, we use constant potential molecular dynamics to study the solvation contributions, and we uncover huge effects of the shape and orientation (on top of the already known size effect) of small hydrophobic and amphiphilic solutes on their adsorption free energy. We propose a minimal theoretical model, the S.O.S. model, that accounts for size, orientation, and shape effects. These novel aspects are rationalized by recasting the concepts at the base of the Lum-Chandler-Weeks theory of hydrophobic solvation (for small solutes in the so-called volume-dominated regime) into a layer-by-layer form, where the properties of each interfacial region close to the metal are explicitly taken into account.

Solvent-pair surfactants enabled assembly of clusters and copolymers towards programmed mesoporous metal oxides

Organic-inorganic molecular assembly has led to numerous nano/mesostructured materials with fantastic properties, but it is dependent on and limited to the direct interaction between host organic structure-directing molecules and guest inorganic species. Here, we report a "solvent-pair surfactants" enabled assembly (SPEA) method to achieve a general synthesis of mesostructured materials requiring no direct host-guest interaction. Taking the synthesis of mesoporous metal oxides as an example, the dimethylformamide/water solvent pairs behave as surfactants and induce the formation of mesostructured polyoxometalates/copolymers nanocomposites, which can be converted into metal oxides. This SPEA method enables the synthesis of functional ordered mesoporous metal oxides with different pore sizes, structures, compositions and tailored pore-wall microenvironments that are difficult to access via conventional direct organic-inorganic assembly. Typically, nitrogen-doped mesoporous ε-WO3 with high specific surface area, uniform mesopores and stable framework is obtained and exhibits great application potentials such as gas sensing.

Cavity formation at metal-water interfaces

The free energy cost of forming a cavity in a solvent is a fundamental concept in rationalizing the solvation of molecules and ions. A detailed understanding of the factors governing cavity formation in bulk solutions has inter alia enabled the formulation of models that account for this contribution in coarse-grained implicit solvation methods. Here, we employ classical molecular dynamics simulations and multistate Bennett acceptance ratio free energy sampling to systematically study cavity formation at a wide range of metal-water interfaces. We demonstrate that the obtained size- and position-dependence of cavitation energies can be fully rationalized by a geometric Gibbs model, which considers that the creation of the metal-cavity interface necessarily involves the removal of interfacial solvent. This so-called competitive adsorption effect introduces a substrate dependence to the interfacial cavity formation energy that is missed in existing bulk cavitation models. Using expressions from scaled particle theory, this substrate dependence is quantitatively reproduced by the Gibbs model through simple linear relations with the adsorption energy of a single water molecule. Besides providing a better general understanding of interfacial solvation, this paves the way for the derivation and efficient parametrization of more accurate interface-aware implicit solvation models needed for reliable high-throughput calculations toward improved electrocatalysts.

On the Trail of Molecular Hydrophilicity and Hydrophobicity at Aqueous Interfaces

Uncovering microscopic hydrophilicity and hydrophobicity at heterogeneous aqueous interfaces is essential as it dictates physico/chemical properties such as wetting, the electrical double layer, and reactivity. Several molecular and spectroscopic descriptors were proposed, but a major limitation is the lack of connections between them. Here, we combine density functional theory-based MD simulations (DFT-MD) and SFG spectroscopy to explore how interfacial water responds in contact with self-assembled monolayers (SAM) of tunable hydrophilicity. We introduce a microscopic metric to track the transition from hydrophobic to hydrophilic interfaces. This metric combines the H/V descriptor, a structural descriptor based on the preferential orientation within the water network in the topmost binding interfacial layer (BIL) and spectroscopic fingerprints of H-bonded and dangling OH groups of water carried by BIL-resolved SFG spectra. This metric builds a bridge between molecular descriptors of hydrophilicity/hydrophobicity and spectroscopically measured quantities and provides a recipe to quantitatively or qualitatively interpret experimental SFG signals.

Individual adsorption of low volatility pheromones: Amphiphilic molecules on a clean water-air interface

Environmental conditions can alter olfactory scent and chemical communication among biological species. In particular, odorant molecules interact with aerosols. Thermodynamics variables governing the adsorption from air to water surface of bombykol, the most studied pheromone, and of three derivative molecules, bombykal, bombykoic acid, and bombykyle acetate, are computed by steered and un-biased molecular dynamics in order to compare the role of their polar head group on adsorption on aqueous aerosols. When adsorbed, the molecule center of mass stands at about 1.2 Å from the interface and oscillates on the same length scale, trapped in an energy well. Gibbs energy of adsorption and desorption time of bombykol are found to be 9.2 k_BT and 59 µs, respectively. The following ordering between the molecules is observed, reading from the more to the least adsorbed: bombykoic acid > bombykol > bombykoic acetate > bombykal. It originates from a complex interplay of entropy and enthalpy. The entropy and enthalpy of adsorption are discussed in the light of structural arrangement, H-bonding, and hydrophilic tail positioning of the molecules at the interface. Our results show that, when dispersed in the air, pheromones adsorb on aqueous aerosols. However, the individual residence time is quite short on pure water surfaces. Aerosols can, therefore, only have a decisive influence on chemical communication through collective effects or through their chemical composition that is generally more complex than that of a pure water surface.

Spectroscopic Fingerprints of Cavity Formation and Solute Insertion as a Measure of Hydration Entropic Loss and Enthalpic Gain

Hydration free energies are dictated by a subtle balance of hydrophobic and hydrophilic interactions. We present here a spectroscopic approach, which gives direct access to the two main contributions: Using THz-spectroscopy to probe the frequency range of the intermolecular stretch (150-200 cm-1 ) and the hindered rotations (450-600 cm-1 ), the local contributions due to cavity formation and hydrophilic interactions can be traced back. We show that via THz calorimetry these fingerprints can be correlated 1 : 1 with the group specific solvation entropy and enthalpy. This allows to deduce separately the hydrophobic (i.e. cavity formation) and hydrophilic contributions to thermodynamics, as shown for hydrated alcohols as a case study. Accompanying molecular dynamics simulations quantitatively support our experimental results. In the future our approach will allow to dissect hydration contributions in inhomogeneous mixtures and under non-equilibrium conditions.

Binding of polar and hydrophobic molecules at the LiCoO2 (001)-water interface: force field development and molecular dynamics simulations

A classical model in the framework of the INTERFACE force field has been developed for treating the LiCoO₂ (LCO) (001)/water interface. In comparison to ab initio molecular dynamics (MD) simulations based on density functional theory, MD simulations using the classical model lead to generally reliable descriptions of interfacial properties, such as the density distribution of water molecules. Water molecules in close contact with the LCO surface form a strongly adsorbed layer, which leads to a free energy barrier for the adsorption of polar or charged molecules to the LCO surface. Moreover, due to the strong hydrogen bonding interactions with the LCO surface, the first water layer forms an interface that exhibits hydrophobic characters, leading to favorable adsorption of non-polar molecules to the interface. Therefore, despite its highly polar nature, the LCO (001) surface binds not only polar/charged but also non-polar solutes. As an application, the model is used to analyze the adsorption of reduced nicotinamide adenine dinucleotide (NADH) and its molecular components to the LCO (001) surface in water. The results suggest that recently observed redox activity of NADH at the LCO/water interface was due to the co-operativity between the ribose component, which drives binding to the LCO surface, and the nicotinamide moiety, which undergoes oxidation.

Unexpectedly High Capacitance of the Metal Nanoparticle/Water Interface: Molecular-Level Insights into the Electrical Double Layer

The electrical double-layer plays a key role in important interfacial electrochemical processes from catalysis to energy storage and corrosion. Therefore, understanding its structure is crucial for the progress of sustainable technologies. We extract new physico-chemical information on the capacitance and structure of the electrical double-layer of platinum and gold nanoparticles at the molecular level, employing single nanoparticle electrochemistry. The charge storage ability of the solid/liquid interface is larger by one order-of-magnitude than predicted by the traditional mean-field models of the double-layer such as the Gouy-Chapman-Stern model. Performing molecular dynamics simulations, we investigate the possible relationship between the measured high capacitance and adsorption strength of the water adlayer formed at the metal surface. These insights may launch the active tuning of solid-solvent and solvent-solvent interactions as an innovative design strategy to transform energy technologies towards superior performance and sustainability.

Size-Dependent Order-Disorder Crossover in Hydrophobic Hydration: Comparison between Spherical Solutes and Linear Alcohols

Theory and computer simulation studies have predicted that water molecules around hydrophobic molecules should undergo an order-disorder transition with increasing solute size around a 1 nm length scale. Some theories predict the formation of a clathrate-like ordered structure around smaller hydrophobic solutes (<1 nm) and the formation of disordered vapor-liquid interfaces around larger solutes (>1 nm) and surfaces. Experimental validation of these predictions has often been elusive and contradictory. High-resolution Raman spectroscopy has detected that water around small hydrophobic solutes shows a signature similar to that of bulk water at lower temperature (increased ordering and a stronger hydrogen-bonded network). Similarly, water around larger solutes shows an increasing population of dangling OH bonds very similar to higher temperature bulk water. Thus, the solute size dependence of the structure and dynamics of water around hydrophobic molecules seems to have an analogy with the temperature dependence in bulk water. In this work, using atomistic classical molecular dynamics (MD) simulations, we have systematically investigated this aspect and characterized this interesting analogy. Structural order parameters including the tetrahedral order parameter (Q), hydrogen bond distribution, and vibrational power spectrum highlight this similarity. However, in contrast to the experimental observations, we do not observe any length-dependent crossover for linear hydrophobic alcohols (n-alkanols) using classical MD simulations. This is in agreement with earlier findings that linear alkane chains do not demonstrate the length-dependent order-disorder transition due to the presence of a sub-nanometer length scale along the cross section of the chain. Moreover, the collapsed state of linear hydrocarbon chains is not significantly populated for smaller chains (number of carbons below 20). In the context of our computational results, we raise several pertinent questions related to the sensitivity of various structural and dynamical parameters toward capturing these complex phenomena of hydrophobic hydration.

Evidence for Entropically Controlled Interfacial Hydration in Mesoporous Organosilicas

At aqueous interfaces, the distribution and dynamics of adsorbates are modulated by the behavior of interfacial water. Hydration of a hydrophobic surface can store entropy via the ordering of interfacial water, which contributes to the Gibbs energy of solute binding. However, there is little experimental evidence for the existence of such entropic reservoirs, and virtually no precedent for their rational design in systems involving extended interfaces. In this study, two series of mesoporous silicas were modified in distinct ways: (1) progressively deeper thermal dehydroxylation, via condensation of surface silanols, and (2) increasing incorporation of nonpolar organic linkers into the silica framework. Both approaches result in decreasing average surface polarity, manifested in a blue-shift in the fluorescence of an adsorbed dye. For the inorganic silicas, hydrogen-bonding of water becomes less extensive as the number of surface silanols decreases. Overhauser dynamic nuclear polarization (ODNP) relaxometry indicates enhanced surface water diffusivity, reflecting a loss of enthalpic hydration. In contrast, organosilicas show a monotonic decrease in surface water diffusivity with decreasing polarity, reflecting enhanced hydrophobic hydration. Molecular dynamics simulations predict increased tetrahedrality of interfacial water for the organosilicas, implying increased ordering near the nm-size organic domains (relative to inorganic silicas, which necessarily lack such domains). These findings validate the prediction that hydrophobic hydration at interfaces is controlled by the microscopic length scale of the hydrophobic regions. They further suggest that the hydration thermodynamics of structurally heterogeneous silica surfaces can be tuned to promote adsorption, which in turn tunes the selectivity in catalytic reactions.

The role of hydrophobic hydration in the free energy of chemical reactions at the gold/water interface: Size and position effects

Metal/water interfaces catalyze a large variety of chemical reactions, which often involve small hydrophobic molecules. In the present theoretical study, we show that hydrophobic hydration at the Au(100)/water interface actively contributes to the reaction free energy by up to several hundreds of meV. This occurs either in adsorption/desorption reaction steps, where the vertical distance from the surface changes in going from reactants to products, or in addition and elimination reaction steps, where two small reactants merge into a larger product and vice versa. We find that size and position effects cannot be captured by treating them as independent variables. Instead, their simultaneous evaluation allows us to map the important contributions, and we provide examples of their combinations for which interfacial reactions can be either favored or disfavored. By taking a N₂ and a CO₂ reduction pathway as test cases, we show that explicitly considering hydrophobic effects is important for the selectivity and rate of these relevant interfacial processes.

A molecular perspective on induced charges on a metallic surface

Understanding the response of the surface of metallic solids to external electric field sources is crucial to characterize electrode-electrolyte interfaces. Continuum electrostatics offer a simple description of the induced charge density at the electrode surface. However, such a simple description does not take into account features related to the atomic structure of the solid and to the molecular nature of the solvent and of the dissolved ions. In order to illustrate such effects and assess the ability of continuum electrostatics to describe the induced charge distribution, we investigate the behavior of a gold electrode interacting with sodium or chloride ions fixed at various positions, in a vacuum or in water, using all-atom constant-potential classical molecular dynamics simulations. Our analysis highlights important similarities between the two approaches, especially under vacuum conditions and when the ion is sufficiently far from the surface, as well as some limitations of the continuum description, namely, neglecting the charges induced by the adsorbed solvent molecules and the screening effect of the solvent when the ion is close to the surface. While the detailed features of the charge distribution are system-specific, we expect some of our generic conclusions on the induced charge density to hold for other ions, solvents, and electrode surfaces. Beyond this particular case, the present study also illustrates the relevance of such molecular simulations to serve as a reference for the design of improved implicit solvent models of electrode-electrolyte interfaces.

Stripping away ion hydration shells in electrical double-layer formation: Water networks matter

For centuries the double layer at the solid/electrolyte interface has been a central concept in electrochemistry. Today, it is still crucial for virtually all renewable energy storage and conversion technologies. Here, the double-layer formation is probed by THz spectroscopy with ultrabright synchrotron light as a source. Our results capture the molecular details of double-layer formation at positively/negatively charged Au electrodes for an NaCl electrolyte. We reveal a contrasting response applying positive versus negative bias, which is dictated by the interfacial water network and rationalized by accompanying molecular dynamics simulations and electronic-structure calculations. While Na⁺ is directly attracted toward the negatively charged electrode, stripping of the Cl⁻ hydration shell is observed only at larger potential values.

The double layer at the solid/electrolyte interface is a key concept in electrochemistry. Here, we present an experimental study combined with simulations, which provides a molecular picture of the double-layer formation under applied voltage. By THz spectroscopy we are able to follow the stripping away of the cation/anion hydration shells for an NaCl electrolyte at the Au surface when decreasing/increasing the bias potential. While Na⁺ is attracted toward the electrode at the smallest applied negative potentials, stripping of the Cl⁻ hydration shell is observed only at higher potential values. These phenomena are directly measured by THz spectroscopy with ultrabright synchrotron light as a source and rationalized by accompanying molecular dynamics simulations and electronic-structure calculations.

Cation enrichment in the ion atmosphere is promoted by local hydration of DNA

Electrostatic interactions are central to the structure and function of nucleic acids, including their folding, condensation, and interaction with proteins and other charged molecules. These interactions are profoundly affected by ions surrounding nucleic acids, the constituents of the so-called ion atmosphere. Here, we report precise Fourier Transform-Terahertz/Far-Infrared (FT-THz/FIR) measurements in the frequency range 30-500 cm^-1 for a 24-bp DNA solvated in a series of alkali halide (NaCl, NaF, KCl, CsCl, and CsF) electrolyte solutions which are sensitive to changes in the ion atmosphere. Cation excess in the ion atmosphere is detected experimentally by observation of cation modes of Na⁺, K⁺, and Cs⁺ in the frequency range between 70-90 cm^-1. Based on MD simulations, we propose that the magnitude of cation excess (which is salt specific) depends on the ability of the electrolyte to perturb the water network at the DNA interface: In the NaF atmosphere, the ions reduce the strength of interactions between water and the DNA more than in case of a NaCl electrolyte. Here, we explicitly take into account the solvent contribution to the chemical potential in the ion atmosphere: A decrease in the number of bound water molecules in the hydration layer of DNA is correlated with enhanced density fluctuations, which decrease the free energy cost of ion-hydration, thus promoting further ion accumulation within the DNA atmosphere. We propose that taking into account the local solvation is crucial for understanding the ion atmosphere.

Effect of the metallicity on the capacitance of gold-aqueous sodium chloride interfaces

Electrochemistry experiments have established that the capacitance of electrode-electrolyte interfaces is much larger for good metals, such as gold and platinum, than for carbon-based materials. Despite the development of elaborate electrode interaction potentials, to date molecular dynamics simulations are not able to capture this effect. Here, we show that changing the width of the Gaussian charge distribution used to represent the atomic charges in gold is an effective way to tune its metallicity. Larger Gaussian widths lead to a capacitance of aqueous solutions (pure water and 1 M NaCl) in good agreement with recent ab initio molecular dynamics results. For pure water, the increase in the capacitance is not accompanied by structural changes, while in the presence of salt, the Na⁺ cations tend to adsorb significantly on the surface. For a strongly metallic gold electrode, these ions can even form inner sphere complexes on hollow sites of the surface.