Interfacial water: A first principles molecular dynamics study of a nanoscale water film on salt

Crumbling crystals: on the dissolution mechanism of NaCl in water

Dissolution of ionic salts in water is ubiquitous, particularly for NaCl. However, an atomistic scale understanding of the process remains elusive. Simulations lend themselves conveniently to studying dissolution since they provide the spatio-temporal resolution that can be difficult to obtain experimentally. Nevertheless, the complexity of various inter- and intra-molecular interactions require careful treatment and long time scale simulations, both of which are typically hindered by computational expense. Here, we use advances in machine learning potential methodology to resolve at an ab initio level of theory the dissolution mechanism of NaCl in water. The picture that emerges is that of a steady ion-wise unwrapping of the crystal preceding its rapid disintegration, reminiscent of crumbling. The onset of crumbling can be explained by a strong increase in the ratio of the surface area to volume of the crystal. Overall, dissolution comprises a series of highly dynamical microscopic sub-processes, resulting in an inherently stochastic mechanism. These atomistic level insights contribute to the general understanding of dissolution mechanisms in other crystals, and the methodology is primed for more complex systems of recent interest such as water/salt interfaces under flow and salt crystals under confinement.

Determination of the Thickness of Interfacial Water by Time-Resolved Sum-Frequency Generation Vibrational Spectroscopy

The physics and chemistry of a charged interface are governed by the structure of the electrical double layer (EDL). Determination of the interfacial water thickness (diw) of the charged interface is crucial to quantitatively describe the EDL structure, but it can be utilized with very scarce experimental methods. Here, we propose and verify that the vibrational relaxation time (T1) of the OH stretching mode at 3200 cm-1, obtained by time-resolved sum frequency generation vibrational spectroscopy with ssp polarizations, provides an effective tool to determine diw. By investigating the T1 values at the SiO2/NaCl solution interface, we established a time-space (T1-diw) relationship. We find that water has a T1 lifetime of ≥0.5 ps for diw ≤ 3 Å, while it displays bulk-like dynamics with T1 ≤ 0.2 ps for diw ≥ 9 Å. T1 decreases as diw increases from ∼3 Å to 9 Å. The hydration water at the DPPG lipid bilayer and LK15β protein interfaces has a thickness of ≥9 Å and shows a bulk-like feature. The time-space relationship will provide a novel tool to pattern the interfacial topography and heterogeneity in Ångstrom-depth resolution by imaging the T1 values.

How Far Is "Bulk Water" from Interfaces? Depends on the Nature of the Surface and What We Measure

Using systematic molecular dynamics (MD) simulations, we revisit the question: At what distance from an interface do the properties of "bulk water" get recovered? We have considered three different kinds of interfaces: nonpolar (hydrophobic; isooctane-water interface), charged (negative; AOT bilayer), and polar (zwitterionic; POPC bilayer). In order to interrogate the extent of perturbation of the interfacial water molecules as a function of the distance from the interface, we utilize a diverse range of structural and dynamical parameters. To capture the structural perturbations, we look into local density (translational order), local tetrahedral order parameter, and dipolar orientation of the water molecules. We also explore the anisotropic diffusion of the water molecules in the direction perpendicular to the interface as well as the planar diffusion parallel to the interface in a distance dependent manner. In addition, the orientational time correlation functions have been computed to understand the extent of slowdown in the rotational dynamics. As expected, the electrostatic field emanating from the charged AOT interface seems to have the highest long-range effect on the orientational order and dynamics of the water molecules, whereas specific interactions like hydrogen bonding and electrostatic interaction lead to significant trapping and kinetic slowdown for both AOT and POPC (zwitterionic) very close to the interface. Our analysis highlights that not only the length-scale of perturbation depends on the nature of the interfaces and specific interactions but also the type of water property that we measure/calculate. Different water properties seem to have widely different length-scale of perturbation. Orientational order parameters seem to be perturbed to a much longer length-scale as compared to translational order parameters. The global orientational order of water can be perturbed even up to ∼4-5 nm near the negatively charged AOT surface in the absence of any extra electrolyte. This observation has significant implication toward the interpretation of experimental measurements as well since different spectroscopic techniques would probe different parameters or water properties with possible mutual disagreement and inconsistency between different types of measurements. Thus, our study provides a broader and unifying perspective toward the aspect of "context dependent" structural and dynamical perturbation of "interfacial water".

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.

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.

Ion Distribution and Hydration Structure at Solid-Liquid Interface between NaCl Crystal and Its Solution

The interface structure between NaCl crystal and its solution has been investigated at the saturated concentration of 298 K by molecular dynamics simulations. We have found that there are many fine structures at this complex interface. Near the surface of crystal, most of Na⁺ only coordinate with water molecules, while almost all Cl^- coordinate with Na⁺ in addition to water molecules. An ion coordinating with more water molecules is farther away from the epitaxial position of lattice. As approaching to the interface, the first hydration shell of ions has the tendency of being ordered, while the orientation of dipole of water molecules in the first hydration shell becomes more disordered than that in the solution. Generally, the first hydration shell of Na⁺ is less affected by nearest Cl^-, whereas the first hydration shell of Cl^- is significantly affected by nearest Na⁺.

Perspective: Structure and dynamics of water at surfaces probed by scanning tunneling microscopy and spectroscopy

The detailed and precise understanding of water-solid interaction largely relies on the development of atomic-scale experimental techniques, among which scanning tunneling microscopy (STM) has proven to be a noteworthy example. In this perspective, we review the recent advances of STM techniques in imaging, spectroscopy, and manipulation of water molecules. We discuss how those newly developed techniques are applied to probe the structure and dynamics of water at solid surfaces with single-molecule and even submolecular resolution, paying particular attention to the ability of accessing the degree of freedom of hydrogen. In the end, we present an outlook on the directions of future STM studies of water-solid interfaces as well as the challenges faced by this field. Some new scanning probe techniques beyond STM are also envisaged.

Wetting of water on hexagonal boron nitride@Rh(111): a QM/MM model based on atomic charges derived for nano-structured substrates

The wetting of water on corrugated and flat hexagonal boron nitride (h-BN) monolayers on Rh(111) is studied within a hybrid quantum mechanics/molecular mechanics (QM/MM) approach. Water is treated by QM methods, whereas the interactions between liquid and substrate are described at the MM level. The electrostatic properties of the substrate are reproduced by assigning specifically generated partial charges to each MM atom. We propose a method to determine restrained electrostatic potential (RESP) charges that can be applied to periodic systems. Our approach is based on the Gaussian and plane waves algorithm and allows an easy tuning of charges for nano-structured substrates. We have successfully applied it to reproduce the electrostatic potential of the corrugated and flat h-BN/Rh(111) known as nanomesh. Molecular dynamics simulations of water films in contact with these substrates are performed and structural and dynamic properties of the interfaces are analyzed. Based on this analysis and on the interaction energies between water film and substrate, we found that the corrugated nanomesh is slightly more hydrophilic. On a macroscopic scale, we expect a smaller contact angle for a droplet on the corrugated surface.

Dispersion corrected RPBE studies of liquid water

The structure of liquid water has been addressed by ab initio molecular dynamics simulations based on density functional theory. Exchange-correlation effects have been described by the popular PBE and RPBE functionals within the generalized gradient approximation as these functionals also yield satisfactory results for metals which is important to model electrochemical interfaces from first principles. In addition, dispersive interactions are included by using dispersion-corrected schemes. It turns out that the dispersion-corrected RPBE functional reproduces liquid water properties quite well in contrast to the PBE functional. This is caused by the replacement of the over-estimated directional hydrogen-bonding in the PBE functional by non-directional dispersive interactions.

Ab initio Kinetic Monte Carlo simulations of dissolution at the NaCl-water interface

We have used ab initio molecular dynamics (AIMD) simulations to study the interaction of water with the NaCl surface. As expected, we find that water forms several ordered hydration layers, with the first hydration layer having water molecules aligned so that oxygen atoms are on average situated above Na sites. In an attempt to understand the dissolution of NaCl in water, we have then combined AIMD with constrained barrier searches, to calculate the dissolution energetics of Na(+) and Cl(-) ions from terraces, steps, corners and kinks of the (100) surface. We find that the barrier heights show a systematic reduction from the most stable flat terrace sites, through steps to the smallest barriers for corner and kink sites. Generally, the barriers for removal of Na(+) ions are slightly lower than for Cl(-) ions. Finally, we use our calculated barriers in a Kinetic Monte Carlo as a first order model of the dissolution process.

An unconventional bilayer ice structure on a NaCl(001) film

Water-solid interactions are of broad importance both in nature and technology. The hexagonal bilayer model based on the Bernal-Fowler-Pauling ice rules has been widely adopted to describe water structuring at interfaces. Using a cryogenic scanning tunnelling microscope, here we report a new type of two-dimensional ice-like bilayer structure built from cyclic water tetramers on an insulating NaCl(001) film, which is completely beyond this conventional bilayer picture. A novel bridging mechanism allows the interconnection of water tetramers to form chains, flakes and eventually a two-dimensional extended ice bilayer containing a regular array of Bjerrum D-type defects. Ab initio density functional theory calculations substantiate this bridging growth mode and reveal a striking proton-disordered ice structure. The formation of the periodic Bjerrum defects with unusually high density may have a crucial role as H donor sites in directing multilayer ice growth and in catalysing heterogeneous chemical reactions on water-coated salt surfaces.

Dissolution of NaCl nanocrystals: an ab initio molecular dynamics study

NaCl nanocrystal dissolution was investigated in atomistic detail revealing a difference in the solvation of two different ionic species.

Ab initio molecular dynamics simulation on the formation process of He@C₆₀ synthesized by explosion

The applications of endohedral non-metallic fullerenes are limited by their low production rate. Recently, an explosive method developed in our group shows promise to prepare He@C₆₀ at fairly high yield, but the mechanism of He inserting into C₆₀ cage at explosive conditions was not clear. Here, ab initio molecular dynamics analysis has been used to simulate the collision between C₆₀ molecules at high-temperature and high-pressure induced by explosion. The results show that defects formed on the fullerene cage by collidsion can effectively decrease the reaction barrier for the insertion of He into C₆₀, and the self-healing capability of the defects was also observed.

Salt effects on water/hydrophobic liquid interfaces: a molecular dynamics study

We investigate structural and dynamical properties of NaCl and KCl solution/n-decane interfaces. We find that K(+) and Cl(-) ions localize close to the interface and modify the orientations of both first-layer and second-layer water molecules in the interfacial region. Na(+) and Cl(-) ions are also localized; however, they do not affect the orientation significantly. NaCl salts cause an increase of the surface tension and the residence time of interfacial water molecules. They also decrease the residence time of interfacial n-decane molecules. KCl salt has similar effects, although to a smaller extent. These findings may help us to better understand biological interfaces in physiological solutions.

Molecular simulation of water confined in nanoporous silica

This paper reports on a molecular simulation study of the thermodynamics, structure and dynamics of water confined at ambient temperature in hydroxylated silica nanopores of a width H = 10 and 20 Å. The adsorption isotherms for water in these nanopores resemble those observed for experimental samples; the adsorbed amount increases continuously in the multilayer adsorption regime until a jump occurs due to capillary condensation of the fluid within the pore. Strong layering of water in the vicinity of the silica surfaces is observed as marked density oscillations are observed up to 8 Å from the surface in the density profiles for confined water. Our results indicate that water molecules within the first adsorbed layer tend to adopt a H-down orientation with respect to the silica substrate. For all pore sizes and adsorbed amounts, the self-diffusivity of confined water is lower than the bulk, due to the hydrophilic interaction between the water molecules and the hydroxylated silica surface. Our results also suggest that the self-diffusivity of confined water is sensitive to the adsorbed amount.