A molecular perspective of water at metal interfaces

Predicting heterogeneous ice nucleation with a data-driven approach

Water in nature predominantly freezes with the help of foreign materials through a process known as heterogeneous ice nucleation. Although this effect was exploited more than seven decades ago in Vonnegut's pioneering cloud seeding experiments, it remains unclear what makes a material a good ice former. Here, we show through a machine learning analysis of nucleation simulations on a database of diverse model substrates that a set of physical descriptors for heterogeneous ice nucleation can be identified. Our results reveal that, beyond Vonnegut's connection with the lattice match to ice, three new microscopic factors help to predict the ice nucleating ability. These are: local ordering induced in liquid water, density reduction of liquid water near the surface and corrugation of the adsorption energy landscape felt by water. With this we take a step towards quantitative understanding of heterogeneous ice nucleation and the in silico design of materials to control ice formation.

Computational Prediction of Novel Ice Phases: A Perspective

Although computational prediction of new ice phases is a niche field in water science, the scientific subject itself is representative of two important areas in physical chemistry, namely, statistical thermodynamics and molecular simulations. The prediction of a variety of novel ice phases has also attracted general public interest since the 1980s. In particular, the prediction of low-dimensional ice phases has gained momentum since the confirmation of a number of low-dimensional "computer ice" phases in the laboratory over the past decade. In this Perspective, the research advancements in computational prediction of novel ice phases over the past few years are reviewed. Particular attention is placed on new ice phases whose physical properties or dimensional structures are distinctly different from conventional bulk ices. Specific topics include the (i) formation of superionic ices, (ii) electrofreezing of water under high pressure and in a high external electric field, (iii) prediction of low-density porous ice at strongly negative pressure, (iv) ab initio computational study of two-dimensional (2D) ice under nanoscale confinement, and (v) 2D ices formed on a solid surface near ambient temperature without nanoscale confinement. Clearly, the formation of most of these novel ice phases demands certain extreme conditions. Ongoing challenges and new opportunities for predicting new ice phases from either classical molecular dynamics simulation or high-level ab initio computation are discussed.

Uniform, Anticorrosive, and Antiabrasive Coatings on Metallic Surfaces for Cation-Metal and Cation-π Interactions

Metals are widely used, from daily life to modern industry. Great efforts have been made to protect the metals with various coatings. However, the well-known conventional electrochemical corrosion induced by cations and the ubiquitous nature of the coffee-ring effect make these processes very difficult. Here, a scheme by two bridges of cations and ethylenediamine (EDA) is proposed to overcome the coffee-ring effect and electrochemical corrosion and experimentally achieve uniform, anticorrosive, and antiabrasive coatings on metallic surfaces. Anticorrosive capability reaches about 26 times higher than that without cation-controlled coatings at 12 h in extremely acidic, high-temperature, and high-humidity conditions and still enhances to 2.7 times over a week. Antiabrasive capability also reaches 2.5 times. Theoretical calculations show that the suspended materials are uniformly adsorbed on the surface mediated by complexed cations through strong cation-metal and cation-π interactions. Notably, the well-known conventional electrochemical corrosion induced by cations is avoided by EDA to control cations solubility in different coating processes. These findings provide a new efficient, cost-effective, facile, and scalable method to fabricate protective coatings on metallic materials and a methodology to study metallic nanostructures in solutions, benefitting practical applications including coatings, printing, dyeing, electrochemical protection, and biosensors.

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.

Interactions of Oxide Surfaces with Water Revealed with Solid-State NMR Spectroscopy

Hydrous materials are ubiquitous in the natural environment and efforts have previously been made to investigate the structures and dynamics of hydrated surfaces for their key roles in various chemical and physical applications, with the help of theoretical modeling and microscopy techniques. However, an overall atomic-scale understanding of the water-solid interface, including the effect of water on surface ions, is still lacking. Herein, we employ ceria nanorods with different amounts of water as an example and demonstrate a new approach to explore the water-surface interactions by using solid-state NMR in combination with density functional theory. NMR shifts and relaxation time analysis provide detailed information on the local structure of oxygen ions and the nature of water motion on the surface: the amount of molecularly adsorbed water decreases rapidly with increasing temperature (from room temperature to 150 °C), whereas hydroxyl groups are stable up to 150 °C, and dynamic water molecules are found to instantaneously coordinate to the surface oxygen ions. The applicability of dynamic nuclear polarization for selective detection of surface oxygen species is also compared to conventional NMR with surface selective isotopic-labeling: the optimal method depends on the feasibility of enrichment and the concentration of protons in the sample. These results provide new insight into the interfacial structure of hydrated oxide nanostructures, which is important to improve performance for various applications.

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.

Mass-zero constrained molecular dynamics for electrode charges in simulations of electrochemical systems

Classical molecular dynamics simulations have recently become a standard tool for the study of electrochemical systems. State-of-the-art approaches represent the electrodes as perfect conductors, modeling their responses to the charge distribution of electrolytes via the so-called fluctuating charge model. These fluctuating charges are additional degrees of freedom that, in a Born-Oppenheimer spirit, adapt instantaneously to changes in the environment to keep each electrode at a constant potential. Here, we show that this model can be treated in the framework of constrained molecular dynamics, leading to a symplectic and time-reversible algorithm for the evolution of all the degrees of freedom of the system. The computational cost and the accuracy of the new method are similar to current alternative implementations of the model. The advantage lies in the accuracy and long term stability guaranteed by the formal properties of the algorithm and in the possibility to systematically introduce additional kinematic conditions of arbitrary number and form. We illustrate the performance of the constrained dynamics approach by enforcing the electroneutrality of the electrodes in a simple capacitor consisting of two graphite electrodes separated by a slab of liquid water.

A Helical Monolayer Ice

Understanding water wetting layers on solid surfaces is essential for many natural and industrial processes. Here we find a helical ice monolayer with every six water molecules helically arranged along the normal of the basal plane by performing an intensive structural search based on ab initio calculations. The helical ice is more stable than all previous models of monolayer and bilayer ices in a wide range of water densities both in vacuum and on weakly interacting substrates due to a stronger network of hydrogen bonds enabled by the helical geometry. More compelling is the fact that this model adequately explains a recent experimental ice monolayer grown on graphite in terms of the lattice parameter, water density, and Moiré pattern. The helical character in the new ice model echoes previously reported helical motifs in one-dimensional ice structures and suggests an unexpected capability of hydrogen bonds in driving the surface reconstruction of ice structures.

Solid wetting-layers in inorganic nano-reactors: the water in imogolite nanotube case

By combined use of wide-angle X-ray scattering, thermo-gravimetric analysis, inelastic neutron scattering, density functional theory and density functional theory molecular dynamics simulations, we investigate the structure, dynamics and stability of the water wetting-layer in single-walled aluminogermanate imogolite nanotubes (SW Ge-INTs): an archetypal system for synthetically controllable and monodisperse nano-reactors. We demonstrate that the water wetting-layer is strongly bound and solid-like up to 300 K under atmospheric pressure, with dynamics markedly different from that of bulk water. Atomic-scale characterisation of the wetting-layer reveals organisation of the H₂O molecules in a curved triangular sublattice stabilised by the formation of three H-bonds to the nanotube's inner surface, with covalent interactions sufficiently strong to promote energetically favourable decoupling of the H₂O molecules in the adlayer. The evidenced changes in the local composition, structure, electrostatics and dynamics of the Ge-INT's inner surface upon the formation of the solid wetting-layer demonstrate solvent-mediated functionalisation of the nanotube's cavity at room temperature and pressure, suggesting new strategies for the design of nano-rectors towards potential control of chemical reactivity in nano-confined volumes.

Ionization of Water as an Effect of Quantum Delocalization at Aqueous Electrode Interfaces

The enhanced probability of water dissociation at the aqueous electrode interfaces is predicted by path-integral ab initio molecular dynamics. The ionization process is observed at the aqueous platinum interface when nuclear quantum effects are introduced in the statistical sampling, while minor effects have been observed at the gold interface. We characterize the dissociation mechanism of the formed water ions. In spite of the fact that the concentration and lifetime of the ions might be challenging to experimentally detect, they may serve as a guide to future experiments. Our observation might have a significant impact on the understanding of electrochemical processes occurring at the metal electrode surface.

What is the trigger for the hydrogen evolution reaction? - towards electrocatalysis beyond the Sabatier principle

The mechanism of the hydrogen evolution reaction, although intensively studied for more than a century, remains a fundamental scientific challenge. Many important questions are still open, making it elusive to establish rational principles for electrocatalyst design. In this work, a comprehensive investigation was conducted to identify which dynamic phenomena at the electrified interface are prerequisite for the formation of molecular hydrogen. In fact, what we observe as an onset of the macroscopic faradaic current originates from dynamic structural changes in the double layer, which are entropic in nature. Based on careful analysis of the activation process, an electrocatalytic descriptor is introduced, evaluated and experimentally confirmed. The catalytic activity descriptor is named as the potential of minimum entropy. The experimentally verified catalytic descriptor reveals significant potential to yield innovative insights for the design of catalytically active materials and interfaces.

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.

Implications of the fractional charge of hydroxide at the electrochemical interface

Rational design of materials that efficiently convert electrical energy into chemical bonds will ultimately depend on a thorough understanding of the electrochemical interface at the atomic level. Towards this goal, the use of density functional theory (DFT) at the generalized gradient approximation (GGA) level has been applied widely in the past 15 years. In the calculation of electrochemical reaction energetics using GGA-DFT, it is frequently implicitly assumed that ions in the Helmholtz plane have unit charge. However, the ion charge is observed to be fractional near the interface through both a capacitor model and through Bader charge partitioning. In this work, we show that this spurious charge transfer can be effectively mitigated by continuum charging of the electrolyte. We then show that, similar to hydronium, the observed fractional charge of hydroxide is not due to a GGA level self-interaction error, as the partial charge is observed even when using hybrid level exchange-correlation functionals.

Water Structure and Properties at Hydrophilic and Hydrophobic Surfaces

The properties of water on both molecular and macroscopic surfaces critically influence a wide range of physical behaviors, with applications spanning from membrane science to catalysis to protein engineering. Yet, our current understanding of water interfacing molecular and material surfaces is incomplete, in part because measurement of water structure and molecular-scale properties challenges even the most advanced experimental characterization techniques and computational approaches. This review highlights progress in the ongoing development of tools working to answer fundamental questions on the principles that govern the interactions between water and surfaces. One outstanding and critical question is what universal molecular signatures capture the hydrophobicity of different surfaces in an operationally meaningful way, since traditional macroscopic hydrophobicity measures like contact angles fail to capture even basic properties of molecular or extended surfaces with any heterogeneity at the nanometer length scale. Resolving this grand challenge will require close interactions between state-of-the-art experiments, simulations, and theory, spanning research groups and using agreed-upon model systems, to synthesize an integrated knowledge of solvation water structure, dynamics, and thermodynamics.

Rich topologies of monolayer ices via unconventional electrowetting

Accurate manipulation of a substance on the nanoscale and ultimately down to the level of a single atom or molecule is an ongoing subject of frontier research. Herein, we show that topologies of water monolayers on substrates, in the complete wetting condition, can be manipulated into rich forms of ordered structures via electrowetting. Notably, two new topologies of monolayer ices were identified from our molecular dynamics simulations: one stable below room temperature and the other one having the ability to be stable at room temperature. Moreover, the wettability of the substrate can be tuned from superhydrophobic to superhydrophilic by uniformly changing the charge of each atomic site of the dipole or quadrupole distributed in an orderly manner on the model substrate. At a certain threshold value of the atomic charge, water droplets on the substrate can spread out spontaneously, achieving a complete electrowetting. Importantly, unlike the conventional electrowetting, which involves application of a uniform external electric field, we proposed non-conventional electrowetting, for the first time, by invoking the electric field of dipoles and quadrupoles embedded in the substrate. Moreover, different topologies of water monolayers can be achieved by using the non-conventional electrowetting. A major advantage of the non-conventional electrowetting is that the contact-angle saturation, a long-standing and known limitation in the field of electrowetting, can be overcome by tuning uniformly the lattice atomic charge at the surface, thereby offering a new way to mitigate the contact-angle saturation for various electrowetting applications.

Hydration-Effect-Promoting Ni-Fe Oxyhydroxide Catalysts for Neutral Water Oxidation

Oxygen evolution reaction (OER) catalysts that function efficiently in pH-neutral electrolyte are of interest for biohybrid fuel and chemical production. The low concentration of reactant in neutral electrolyte mandates that OER catalysts provide both the water adsorption and dissociation steps. Here it is shown, using density functional theory simulations, that the addition of hydrated metal cations into a Ni-Fe framework contributes water adsorption functionality proximate to the active sites. Hydration-effect-promoting (HEP) metal cations such as Mg²⁺ and hydration-effect-limiting Ba²⁺ into Ni-Fe frameworks using a room-temperature sol-gel process are incorporated. The Ni-Fe-Mg catalysts exhibit an overpotential of 310 mV at 10 mA cm^-2 in pH-neutral electrolytes and thus outperform iridium oxide (IrO₂ ) electrocatalyst by a margin of 40 mV. The catalysts are stable over 900 h of continuous operation. Experimental studies and computational simulations reveal that HEP catalysts favor the molecular adsorption of water and its dissociation in pH-neutral electrolyte, indicating a strategy to enhance OER catalytic activity.

Interfacial Hydrogen-Bonding Dynamics in Surface-Facilitated Dehydrogenation of Water on TiO2(110)

Molecular-level understanding of the dehydrogenation of interfacial water molecules on metal oxides and their interactive nature relies on the ability to track the motion of light and small hydrogen atoms, which is known to be difficult. Here, we report precise measurements of the surface-facilitated water dehydrogenation process at terminal Ti sites of TiO₂(110) using scanning tunneling microscopy. Our measured hydrogen-bond dynamics of H₂O and D₂O reveal that the vibrational and electronic excitations dominate the sequential transfer of two H (D) atoms from a H₂O (D₂O) molecule to adjacent surface oxygen sites, manifesting the active participation of the oxide surface in the dehydrogenation processes. Our results show that, at the stoichiometric Ti_5c sites, individual H₂O molecules are energetically less stable than the dissociative form, where a barrier is expected to be as small as approximately 70-120 meV on the basis of our experimental and theoretical results. Moreover, our results reveal that interfacial hydrogen bonds can effectively assist H atom transfer and exchange across the surface. The revealed quantitative hydrogen-bond dynamics provide a new atomistic mechanism for water interactions on metal oxides in general.

Nanoscopic diffusion of water on a topological insulator

The microscopic motion of water is a central question, but gaining experimental information about the interfacial dynamics of water in fields such as catalysis, biophysics and nanotribology is challenging due to its ultrafast motion, and the complex interplay of inter-molecular and molecule-surface interactions. Here we present an experimental and computational study of the nanoscale-nanosecond motion of water at the surface of a topological insulator (TI), Bi[Formula: see text]Te[Formula: see text]. Understanding the chemistry and motion of molecules on TI surfaces, while considered a key to design and manufacturing for future applications, has hitherto been hardly addressed experimentally. By combining helium spin-echo spectroscopy and density functional theory calculations, we are able to obtain a general insight into the diffusion of water on Bi[Formula: see text]Te[Formula: see text]. Instead of Brownian motion, we find an activated jump diffusion mechanism. Signatures of correlated motion suggest unusual repulsive interactions between the water molecules. From the lineshape broadening we determine the diffusion coefficient, the diffusion energy and the pre-exponential factor.

A Comprehensive Review on Water Diffusion in Polymers Focusing on the Polymer-Metal Interface Combination

Water diffusion in polymers is relevant to a broad range of physicochemical phenomena and technological processes. Although many fields contributed to rapid progress in the fundamental knowledge of water-polymer interactions, detailed understandings come mainly from interpreting numerous experiments. These studies showed that a remarkably rich variety of diffusion forms between water and even seemingly simple polymers. In this review, focusing on the gravimetric and capacitance method, we discuss contradictions and problems existing for water diffusion in polymers in detail from perspectives of experiments and models, focusing on the analysis of error derived from widely used methods, especially for the Brasher-Kingsbury equation. We also provide a perspective on outstanding problems, challenges, and open questions, including water clusters, relaxation, and electrochemical reactions at the metal/polymer interface, as well as expanding the theoretical prospective.

Charge fluctuations from molecular simulations in the constant-potential ensemble

Statistical mechanics of constant-potential molecular simulations yields a new fluctuation–dissipation relation for the differential capacitance.

Dependence of water adsorption on the surface structure of silicon wafers aged under different environmental conditions

Most materials exposed to ambient air can adsorb water molecules and the adsorption capability strongly depends on the surface property. The water contact angle has been widely used as a measure for surface wettability; however, a question can still be asked whether the water contact angle can be used as an adequate sole predictor for water adsorption on the surface in humid air. In this paper, HF-etched silicon wafers were aged (oxidized) under different environmental conditions at room temperature to grow surface layers with varying water contact angles from ∼0° (fully hydrophilic) to ∼83° (highly hydrophobic), and water adsorption as a function of relative humidity (RH) was studied on such surfaces. The thickness and structure of the adsorbed water layer were found to depend on not only the surface wettability on each surface, but also the history of surface oxidation conditions. In particular, the silicon wafer surface oxidized in liquid water uptakes significantly more water from humid air than the fully-hydroxylated native oxide surface (SiO_x/OH), even though its water contact angle is higher than that on the SiO_x/OH surface. This could be attributed to the formation of a gel-like structure during oxidation in liquid water.

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.

Understanding supported noble metal catalysts using first-principles calculations

Heterogeneous catalysis on supported and nonsupported nanoparticles is of fundamental importance in the energy and chemical conversion industries. Rather than laboratory analysis, first-principles calculations give us an atomic-level understanding of the structure and reactivity of nanoparticles and supports, greatly reducing the efforts of screening and design. However, unlike catalysis on low index single crystalline surfaces, nanoparticle catalysis relies on the tandem properties of a support material as well as the metal cluster itself, often with charge transfer processes being of key importance. In this perspective, we examine current state-of-the-art quantum-chemical research for the modeling of reactions that utilize small transition metal clusters on metal oxide supports. This should provide readers with useful insights when dealing with chemical reactions on such systems, before discussing the possibilities and challenges in the field.

Theoretical investigations of a platinum-water interface using quantum-mechanics-molecular-mechanics based molecular dynamics simulations

Pt-water interfaces have been of immense interest in the field of energy storage and conversion. Studying this interface using both experimental and theoretical tools is challenging. On the theoretical front, typically one uses classical molecular dynamics (MD) simulations to handle large system sizes or time scales while for a more accurate quantum mechanical description Born Oppenheimer MD (BOMD) is typically used. The latter is limited to smaller system sizes and time-scales. In this study using quantum-mechanics-molecular-mechanics (QMMM), we have performed atomistic MD simulations to have a microscopic understanding of the structure of the Pt-water interface using a system size that is much larger than that accessible when using BOMD simulations. In contrast to recent reports using BOMD simulations, our study reveals that the water molecules typically form two distinct layers above the Pt-surface before they form bulk like structures. Further, we also find that a significant fraction of the water molecules at the interface are pointed towards the surface thereby disrupting the H-bond network. Consistent with this observation, the layer resolved oxygen-oxygen radial distribution function for the water molecules belonging to the solvating water layer shows a high density liquid like behaviour even though the overall water behaves like a low density liquid. A charge transfer analysis reveals that this solvating water layer donates electrons to the Pt atoms in contact with it thereby resulting in the formation of an interface dipole that is pointing towards the surface. Our results suggest that, using QMMM-MD, on one hand it is possible to study more realistic models of solid-liquid interfaces that are inaccessible with BOMD, while on the other hand one also has access to information about such systems that are not obtained from conventional classical MD simulations.

Different platinum crystal surfaces show very distinct protein denaturation capabilities

Different platinum (Pt) surfaces of nanocrystals usually exhibit significant distinctions with regard to various biological, physical, and chemical characteristics, such as bio-recognition, surface wetting, and catalytic activities. In this study, we report for the first time that two shape-controlled Pt nanocrystals with the most common low-index surfaces, Pt(100) and Pt(111), show very dissimilar protein denaturation capabilities based on all-atom molecular dynamics simulations employing the widely used model protein, villin headpiece (HP35). We demonstrate that HP35 is well preserved on the Pt(100) crystal surface, whereas it is severely disrupted on the Pt(111) crystal surface. This surprising difference originates from the distinct water behavior in the first solvation shell (FSS) of the two Pt crystal surfaces. Within the FSS of the Pt(100) crystal surface, water molecules form a very compact and stable monolayer through a highly uniform rhombic hydrogen-bond network. This water monolayer prefers the adsorption of acidic residues (such as Glu and Asp) and acts as a shield to prevent other residues from directly coming into contact with the metal surface. On the other hand, the hydrogen bond network in the water monolayer in the FSS of the Pt(111) crystal surface is very sparse and quite defective, which makes it more vulnerable to the penetration of various residues, particularly those with planar side chains such as Phe, Trp and Arg due to strong dispersion interactions, leading to subsequent protein unfolding. The binding free energy calculations for some key amino acids on the two different crystal surfaces further uncover the molecular origin behind their distinct protein denaturation capability. Our study reveals the vital importance of interfacial water in determining the structure of proteins when binding to different metal crystal surfaces. The discovered molecular mechanisms may be helpful for the future development of a bio-assisted programmable synthetic strategy of sophisticated Pt nanostructures for biomedical applications.

Hydrogen bonding structure of confined water templated by a metal-organic framework with open metal sites

Water in confinement exhibits properties significantly different from bulk water due to frustration in the hydrogen-bond network induced by interactions with the substrate. Here, we combine infrared spectroscopy and many-body molecular dynamics simulations to probe the structure and dynamics of confined water as a function of relative humidity within a metal-organic framework containing cylindrical pores lined with ordered cobalt open coordination sites. Building upon the agreement between experimental and theoretical spectra, we demonstrate that water at low relative humidity binds initially to open metal sites and subsequently forms disconnected one-dimensional chains of hydrogen-bonded water molecules bridging between cobalt atoms. With increasing relative humidity, these water chains nucleate pore filling, and water molecules occupy the entire pore interior before the relative humidity reaches 30%. Systematic analysis of rotational and translational dynamics indicates heterogeneity in this pore-confined water, with water molecules displaying variable mobility as a function of distance from the interface.

The properties of water under confinement are significantly altered with respect to the bulk phase. Here the authors use infrared spectroscopy and many-body molecular dynamics simulations to show the structure and dynamics of confined water as a function of relative humidity within a metal-organic framework.

Phase Diagram of Nanoscale Water on Solid Surfaces with Various Wettabilities

Understanding structural and dynamic properties of water in contact with solid surfaces is essential for diverse fields, including environmental sciences, nanofluidics, lubrication, and electrochemistry. Despite tremendous efforts, how interfacial water phase behaviors correlate with a surface's wettability remains elusive. Here, we investigate the structure and dynamics of nanoscale water droplets or adlayers on solid surfaces with wettabilities spanning from strongly hydrophobic to strongly hydrophilic using extensive molecular dynamics simulations. It is shown that liquid water drops on solid surfaces with contact angles greater than 42.6° transform into drops of ordinary hexagonal ice (I_h) upon cooling. In contrast, water forms a liquid disc on a completely wetted surface with a zero contact angle, which freezes into a hexagonal bilayer ice disc at low temperatures. Unexpectedly, on surfaces with a mild contact angle in the range of 21.9°-29.2°, the originally stable liquid drop at room temperature further wets the surface upon cooling and eventually transforms into a bilayer ice disc. These results establish a phase diagram of nanoscale water at the wettability versus temperature plane, which may expand our knowledge of water-surface interactions as well as enrich the complexity of water behaviors at interfaces.

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.

Water nanostructure formation on oxide probed in situ by optical resonances

The dynamic characterization of water multilayers on oxide surfaces is hard to achieve by currently available techniques. Despite this, there is an increasing interest in the evolution of water nanostructures on oxides to fully understand the complex dynamics of ice nucleation and growth in natural and artificial environments. Here, we report the in situ detection of the dynamic evolution of nanoscale water layers on an amorphous oxide surface probed by optical resonances. In the water nanolayer growth process, we find an initial nanocluster morphology that turns into a planar layer beyond a critical thickness. In the reverse process, the planar water film converts to nanoclusters, accompanied by a transition from a planar amorphous layer to crystalline nanoclusters. Our results are explained by a simple thermodynamic model as well as kinetic considerations. Our work represents an approach to reveal the nanostructure and dynamics at the water-oxide interface using resonant light probing.

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.

Molecular Dynamics Simulation of the Influence of Nanoscale Structure on Water Wetting and Condensation

Recent advances in the microfabrication technology have made it possible to control surface properties at micro- and nanoscale levels. Functional surfaces drastically change wettability and condensation processes that are essential for controlling of heat transfer. However, the direct observation of condensation on micro- and nanostructure surfaces is difficult, and further understanding of the effects of the microstructure on the phase change is required. In this research, the contact angle of droplets with a wall surface and the initial condensation process were analyzed using a molecular dynamics simulation to investigate the impact of nanoscale structures and their adhesion force on condensation. The results demonstrated the dependence of the contact angle of the droplets and condensation dynamics on the wall structure and attractive force of the wall surface. Condensed water droplets were adsorbed into the nanostructures and formed a water film in case of a hydrophilic surface.

Replica exchange MD simulations of two-dimensional water in graphene nanocapillaries: rhombic versus square structures, proton ordering, and phase transitions

The hydrogen bond patterns, proton ordering, and phase transitions of monolayer ice in two-dimensional hydrophobic confinement are fundamentally different from those found for bulk ice. To investigate the behavior of quasi-2D ice, we perform molecular dynamics simulations of water confined between fixed graphene plates at a distance of 0.65 nm. While experimental results are still limited and theoretical investigations are often based on a single, often empirically based force field model, this work presents a systematic study modeling the water-graphene interaction by effective Lennard-Jones potentials previously derived from high-level ab initio CCSD(T) calculations of water adsorbed on graphene [Phys. Chem. Chem. Phys., 2013, 15, 4995]. For the water-water interaction different water force fields, i.e. SPCE, TIP3P, TIP4P, TIP4P/ICE, and TIP5P, are used. The water occupancy of the graphene capillary at a pressure of 1000 MPa is determined to be between 13.5 and 13.9 water molecules per square nanometer, depending on the choice of the water force field. Based on these densities, we explore the structure and dynamics of quasi-2D water for temperatures ranging from 200 K to about 600 K for each of the five force fields. To ensure complete sampling of the configurational space and to overcome the barriers separating metastable structures, these simulations are based on the replica exchange molecular dynamics technique. We report different tetragonal hydrogen bond patterns, which are classified as nearly square or as rhombic. While many of these arrangements are essentially flat, in some cases puckered arrangements are found, too. Also the proton ordering of the quasi-2D water structures is considered, allowing us to identify them as ferroelectric, ferrielectric or antiferroelectric. For temperatures between 200 K and 400 K we find several second-order phase transitions from one ice structure to another, changing in many cases both the arrangements of the oxygen atoms and the proton ordering. For temperatures between 400 K and 600 K there are melting-like transitions from a monolayer of ice to a monolayer of liquid water. These first-order phase transitions have a latent heat between 3.4 and 4.0 kJ mol^-1. Both the values of the transition temperatures and of the latent heats display considerable model dependence for the five different water models investigated here.

Direct observation of 2-dimensional ices on different surfaces near room temperature without confinement

Water-solid interfaces play important roles in a wide range of fields, including atmospheric science, geochemistry, electrochemistry, and food science. Herein, we report simulation evidence of 2-dimensional (2D) ice formation on various surfaces and the dependence of the 2D crystalline structure on the hydrophobicity and morphology of the underlying surface. Contrary to the prevailing view that nanoscale confinement is necessary for the 2D liquid-to-bilayer ice transition, we find that the liquid-to-bilayer hexagonal ice (BHI) transition can occur either on a model smooth surface or on model fcc-crystal surfaces with indices of (100), (110), and (111) near room temperature. We identify a critical parameter that characterizes the water-surface interaction, above which the BHI can form on the surface. This critical parameter increases as the temperature increases. Even at temperatures above the freezing temperature of bulk ice (Ih ), we find that BHI can also form on a superhydrophilic surface due to the strong water-surface interaction. The tendency toward the formation of BHI without confinement reflects a proper water-surface interaction that can compensate for the entropy loss during the freezing transition. Furthermore, phase diagrams of 2D ice formation are described on the plane of the adsorption energy versus the fcc lattice constant (Eads-afcc), where 4 monolayer square-like ices are also identified on the fcc model surfaces with distinct water-surface interactions.

Thickness and Structure of Adsorbed Water Layer and Effects on Adhesion and Friction at Nanoasperity Contact

Most inorganic material surfaces exposed to ambient air can adsorb water, and hydrogen bonding interactions among adsorbed water molecules vary depending on, not only intrinsic properties of material surfaces, but also extrinsic working conditions. When dimensions of solid objects shrink to micro- and nano-scales, the ratio of surface area to volume increases greatly and the contribution of water condensation on interfacial forces, such as adhesion (Fa) and friction (Ft), becomes significant. This paper reviews the structural evolution of the adsorbed water layer on solid surfaces and its effect on Fa and Ft at nanoasperity contact for sphere-on-flat geometry. The details of the underlying mechanisms governing water adsorption behaviors vary depending on the atomic structure of the substrate, surface hydrophilicity and atmospheric conditions. The solid surfaces reviewed in this paper include metal/metallic oxides, silicon/silicon oxides, fluorides, and two-dimensional materials. The mechanism by which water condensation influences Fa is discussed based on the competition among capillary force, van der Waals force and the rupture force of solid-like water bridge. The condensed meniscus and the molecular configuration of the water bridge are influenced by surface roughness, surface hydrophilicity, temperature, sliding velocity, which in turn affect the kinetics of water condensation and interfacial Ft. Taking the effects of the thickness and structure of adsorbed water into account is important to obtain a full understanding of the interfacial forces at nanoasperity contact under ambient conditions.

Insight into hydrogen production through molecular simulation of an electrode-ionomer electrolyte system

In this work, we examine metal electrode-ionomer electrolyte systems at high voltage (negative surface charge) and at high pH to assess factors that influence hydrogen production efficiency. We simulate the hydrogen evolution electrode interface investigated experimentally in the work of Bates et al. [J. Phys. Chem. C 119, 5467 (2015)] using a combination of first principles calculations and classical molecular dynamics. With this detailed molecular information, we explore the hypotheses posed in the work of Bates et al. In particular, we examine the response of the system to increased bias voltage and oxide coverage in terms of the potential profile, changes in solvation and species concentrations away from the electrode, surface concentrations, and orientation of water at reactive surface sites. We discuss this response in the context of hydrogen production.

It's not just the defects - a curved crystal study of H2O desorption from Ag

We investigate water desorption from hydrophobic surfaces using two curved Ag single crystals centered at (111) and (001) apices. On these types of crystals the step density gradually increases along the curvature, allowing us to probe large ranges of surface structures in between the (001), (111) and (110) planes. Subtle differences in desorption of submonolayer water coverages point toward structure dependencies in water cluster nucleation. The B-type step on hydrophobic Ag binds water structures more strongly than adjacent (111) planes, leading to preferred desorption from steps. This driving force is smaller for A-type steps on (111) terraces. The A'-type step flanked by (001) terraces shows no indication of preferred desorption from steps. Extrapolation to the (311) surface, not contained within either curved surface, demonstrates that both A- and A'-type steps can be regarded chemically identical for water desorption. The different trends in desorption temperature on the two crystals can thus be attributed to stronger water adsorption at (001) planes than at (111) planes and identical to adsorption at the step. These results show that our approach to studying the structure dependence of water desorption is sensitive to variations in desorption energy smaller than 'chemical accuracy', i.e. 1 kcal mol-1.

In situ probing electrified interfacial water structures at atomically flat surfaces

Solid/liquid interfaces are ubiquitous in nature and knowledge of their atomic-level structure is essential in elucidating many phenomena in chemistry, physics, materials science and Earth science¹. In electrochemistry, in particular, the detailed structure of interfacial water, such as the orientation and hydrogen-bonding network in electric double layers under bias potentials, has a significant impact on the electrochemical performances of electrode materials^2-4. To elucidate the structures of electric double layers at electrochemical interfaces, we combine in situ Raman spectroscopy and ab initio molecular dynamics and distinguish two structural transitions of interfacial water at electrified Au single-crystal electrode surfaces. Towards negative potentials, the interfacial water molecules evolve from structurally 'parallel' to 'one-H-down' and then to 'two-H-down'. Concurrently, the number of hydrogen bonds in the interfacial water also undergoes two transitions. Our findings shed light on the fundamental understanding of electric double layers and electrochemical processes at the interfaces.

Strain Relief during Ice Growth on a Hexagonal Template

Heterogeneous ice nucleation at solid surfaces impacts many areas of science, from environmental processes, such as precipitation, to microbiological systems and food processing, but the microscopic mechanisms underpinning nucleation remain unclear. Discussion of ice growth has often focused around the role of the surface in templating the structure of water, forcing the first layer to adopt the registry of the underlying substrate rather than that of ice. To grow a thick ice film, water in the first few ice layers must accommodate this strain, but understanding how this occurs requires detailed molecular-scale information that is lacking. Here we combine scanning tunneling microscopy, low-energy electron diffraction, and work-function measurements with electronic structure calculations to investigate the initial stages of ice growth on a Pt alloy surface, having a lattice spacing 6% larger than ice. Although the first layer of water forms a strictly commensurate hexagonal network, this behavior does not extend to the second layer. Instead, water forms a 2D structure containing extended defect rows made from face-sharing pentamer and octamer rings. The defect rows allow the majority of second-layer water to remain commensurate with the solid surface while compensating lateral strain by increasing the water density close to that of an ice surface. The observation of octamer-pentamer rows in ice films formed on several surfaces suggests that the octamer-pentamer defect motif acts as a flexible strain relief mechanism in thin ice films, providing a mechanism that is not available during the growth of strained films in other materials, such as semiconductors.

Quasi-Liquid Layer on Ice and Its Effect on the Confined Freezing of Porous Materials

Freezing of the water confined in thin pores can be destructive to the porous frame, but the effect of the quasi-liquid layer (QLL) between the confined ice and the pore walls remains still far from being fully understood. In the present study, the physical origins of the intermediate phase of QLL were discussed by thermodynamic analyses. Different interactions on QLL bring different models to estimate its thickness, which generally decays with temperature decreasing. Four representative models of QLL thickness were selected to unveil its effect on the growing rates and extents of ice in a concrete. The engineering consequences of the confined freezing were then discussed in the aspects of effective pore pressures built from the confined ice growth and deformations framed by a poro-elastic model. Overall, thickening QLL depresses ice growing rates and contents and, consequentially, decreases pore pressures and material deformations during freezing. The QLL corrections also narrow the gaps between the predicted and measured freezing deformations. The findings of this study contribute to profound understandings of confined freezing that may bridge over physical principles and engineering observations.