Structure and dynamics of water near hydrophilic surfaces

Multimodal confined water dynamics in reverse osmosis polyamide membranes

While polyamide (PA) membranes are widespread in water purification and desalination by reverse osmosis, a molecular-level understanding of the dynamics of both confined water and polymer matrix remains elusive. Despite the dense hierarchical structure of PA membranes formed by interfacial polymerization, previous studies suggest that water diffusion remains largely unchanged with respect to bulk water. Here, we employ neutron spectroscopy to investigate PA membranes under precise hydration conditions, and a series of isotopic contrasts, to elucidate water transport and polymer relaxation, spanning ps-ns timescales, and Å-nm lengthscales. We experimentally resolve, for the first time, the multimodal diffusive nature of water in PA membranes: in addition to (slowed down) translational jump-diffusion, we observe a long-range and a localized mode, whose geometry and timescales we quantify. The PA matrix is also found to exhibit rotational relaxations commensurate with the nanoscale confinement observed in water diffusion. This comprehensive ‘diffusion map’ can anchor molecular and nanoscale simulations, and enable the predictive design of PA membranes with tuneable performance.

Polymeric membranes are extensively used in water desalination, but the effect of membrane nanostructure on water transport is still elusive. The authors, using quasi-elastic neutron scattering and contrast variation techniques, provide detailed insight into the dynamics of the polymer network and confined water across a wide range of length and timescales.

Simulations of the IR and Raman spectra of water confined in amorphous silica slit pores

Water in nano-scale confining environments is a key element in many biological, material, and geological systems. The structure and dynamics of the liquid can be dramatically modified under these conditions. Probing these changes can be challenging, but vibrational spectroscopy has emerged as a powerful tool for investigating their behavior. A critical, evolving component of this approach is a detailed understanding of the connection between spectroscopic features and molecular-level details. In this paper, this issue is addressed by using molecular dynamics simulations to simulate the linear infrared (IR) and Raman spectra for isotopically dilute HOD in D₂O confined in hydroxylated amorphous silica slit pores. The effect of slit-pore width and hydroxyl density on the silica surface on the vibrational spectra is also investigated. The primary effect of confinement is a blueshift in the frequency of OH groups donating a hydrogen bond to the silica surface. This appears as a slight shift in the total (measurable) spectra but is clearly seen in the distance-based IR and Raman spectra. Analysis indicates that these changes upon confinement are associated with the weaker hydrogen-bond accepting properties of silica oxygens compared to water molecules.

Quantification of spatio-temporal scales of dynamical heterogeneity of water near lipid membranes above supercooling

A hydrated 1,2-dimyristoyl-sn-glycero-3-phosphorylcholine (DMPC) lipid membrane is investigated using an all atom molecular dynamics simulation at 308 K to determine the physical sources of universal slow relaxations of hydration layers and length-scale of the spatially heterogeneous dynamics. Continuously residing interface water (IW) molecules hydrogen bonded to different moieties of lipid heads in the membrane are identified. The non-Gaussian parameters of all classes of IW molecules show a cross-over from cage vibration to translational diffusion. A significant non-Gaussianity is observed for the IW molecules exhibiting large length correlations in translational van Hove functions. Two time-scales for the ballistic motions and hopping transitions are obtained from the self intermediate scattering functions of the IW molecules with an additional long relaxation, which disappears for bulk water. The long relaxation time-scales for the IW molecules obtained from the self intermediate scattering functions are in good accordance with the hydrogen bond relaxation time-scales irrespective of the nature of the chemical confinement and the confinement lifetime. Employing a block analysis approach, the length-scale of dynamical heterogeneities is captured from a transition from non-Gaussianity to Gaussianity in van Hove correlation functions of the IW molecules. The heterogeneity length-scale is comparable to the wave-length of the small and weak undulations of the membrane calculated by Fourier transforms of lipid tilts. This opens up a new avenue towards a possible correlation between heterogeneity length-scale and membrane curvature more significant for rippled membranes. Thus, our analyses provide a measure towards the spatio-temporal scale of dynamical heterogeneity of confined water near membranes.

Phase Transition in Monolayer Water Confined in Janus Nanopore

The ubiquitous nature of water invariably leads to a variety of physical scenarios that can result in many intriguing properties. We investigate the thermodynamics and associated phase transitions for a water monolayer confined within a quasi-two-dimensional nanopore. An asymmetric nanopore constructed by combining a hydrophilic (hexagonal boron nitride) sheet and a hydrophobic (graphene) sheet leads to an ordered water structure at much higher temperatures compared to a symmetric hydrophobic nanopore consisting of two graphene sheets. The discontinuous change in the thermodynamic quantities, potential energy ( U), and entropy ( S) of confined water molecules computed from the all-atom molecular dynamics simulation trajectories, uncovers a first-order phase transition in the temperature range of T = 320-330 K. Structural analysis reveals that water molecules undergo a disorder-to-order phase transformation in this temperature range with a 4-fold symmetric phase persisting at lower temperatures. Our findings predict a novel confinement system which has the melting transition for monolayer water above the room temperature, and provide a microscopic understanding which will have important implications for other nanofludic systems as well.

Supercooling and freezing processes in nanoconfined water by time-resolved optical Kerr effect spectroscopy

Using heterodyne-detected optical Kerr effect (HD-OKE) measurements, we investigate the vibrational dynamics and the structural relaxation of water nanoconfined in Vycor porous silica samples (pore size ≃ 4 nm) at different levels of hydration and temperatures. At low levels of hydration corresponding to two complete superficial water layers, no freezing occurs and the water remains mobile at all the investigated temperatures with dynamic features similar, but not equal to, the bulk water. The fully hydrated sample shows the formation of ice at about 248 K. This process does not involve all the contained water; a part of it remains in a supercooled phase. The structural relaxation times measured from the decay of the time-dependent HD-OKE signal shows the temperature dependence largely affected by the hydration level; the low frequency (ν < 500 cm(-1)) vibrational spectra obtained by the Fourier transforms of the HD-OKE signal appear less affected by confinement.

Interaction of neurotransmitters with a phospholipid bilayer: a molecular dynamics study

We have performed a series of molecular dynamics simulations to study the interactions between the neurotransmitters (NTs) γ-aminobutyrate (GABA), glycine (GLY), acetylcholine (ACH) and glutamate (GLU) as well as the amidated/acetylated γ-aminobutyrate (GABA(neu)) and the osmolyte molecule glycerol (GOL) with a dipalmitoylphosphatidylcholine (DPPC) bilayer. In agreement with previously published experimental data, we found the lowest membrane affinity for the charged molecules and a moderate affinity for zwitterionic and polar molecules. The affinity can be ranked as follows: ACH-GLU<<GABA<GLY<<GABA(neu)<<GOL. The latter three penetrated the bilayer at most with the deepest location being close to the glycerol backbone of the phospholipids. Even at that position, these solutes were noticeably hydrated and carried ∼30-80% of the bulk water along. The mobility of hydration water at the solute is also affected by the penetration into the bilayer. Two time scales of exchanging water molecules could be determined. In the bulk phase, the hydration layer contains ∼20% slow exchanging water molecules which increases 2-3 times as the solutes entered the bilayer. Our results indicate that there is no intermediate exchange of slow moving water molecules from the solutes to the lipid atoms and vice versa. Instead, the exchange relies on the reservoir of unbounded ("free") water molecules in the interfacial bilayer region. Results from the equilibrium simulations are in good agreement with the results from umbrella sampling simulations, which were conducted for the four naturally occurring NTs. Free energy profiles for ACH and GLU show a minimum of ∼2-3 kJ/mol close to the bilayer interface, while for GABA and GLY, a minimum of respectively ∼2 kJ/mol and ∼5 kJ/mol is observed when these NTs are located in the vicinity of the lipid glycerol backbone. The most important interaction of NTs with the bilayer is the charged amino group of NTs with the lipid phosphate group.

Glassy nature of water in an ultraconfining disordered material: the case of calcium-silicate-hydrate

We present the structural and dynamic nature of water ultraconfined in the quasi-two-dimensional nanopores of the highly disordered calcium-silicate-hydrate (C-S-H), the major binding phase in cement. Our approach is based on classical molecular simulations. We demonstrate that the C-S-H nanopore space is hydrophilic, particularly because of the nonbridging oxygen atoms on the disordered silicate chains which serve as hydrogen-bond acceptor sites, directionally orienting the hydrogen atoms of the interfacial water molecules toward the calcium-silicate layers. The water in this interlayer space adopts a unique multirange structure: a distorted tetrahedral coordination at short range up to 2.7 Å, a disordered structure similar to that of dense fluids and supercooled phases at intermediate range up to 4.2 Å, and persisting spatial correlations through dipole-dipole interactions up to 10 Å. A three-stage dynamics governs the mean square displacement (MSD) of water molecules, with a clear cage stage characteristic of the dynamics in supercooled liquids and glasses, consistent with its intermediate-range structure identified here. At the intermediate time scales corresponding to the β-relaxation of glassy materials, coincident with the cage stage in MSD, the non-Gaussian parameter indicates a significant heterogeneity in the translational dynamics. This dynamic heterogeneity is induced primarily because of the heterogeneity in the distribution of hydrogen bond strengths. The strongly attractive interactions of water molecules with the calcium silicate walls serve to constrain their motion. Our findings have important implications on describing the cohesion and mechanical behavior of cement from its setting to its aging.

A molecular dynamics simulation of the adsorption of water molecules surrounding an Au nanoparticle

This study uses molecular dynamics simulations performed in a parallel computing environment to investigate the adsorption of water molecules surrounding Au nanoparticles of various sizes. An observation of the oxygen and hydrogen atom distributions reveals that the adsorption of the water molecules creates two shell-like formations of water in close vicinity to the Au nanoparticle surface. These shell-like formations are found to be more pronounced around smaller Au nanoparticles. The rearrangement of water molecules in this region reduces the local hydrogen bond strength to below that which is observed in the bulk region. Finally, the simulation results indicate that the absolute value of the interaction energy between the water molecules and the Au nanoparticle is reduced when the water molecules surround a nanoparticle of larger diameter. This observation implies that a stronger adsorption effect exists between smaller Au nanoparticles and water molecules. Hence, the value of the adsorption constant increases for smaller Au nanoparticles.

X-ray Absorption Spectroscopy Measurements of Liquid Water

Recent studies, based on X-ray absorption spectroscopy (XAS) and X-ray Raman scattering (XRS), have shown that the hydrogen bond network in liquid water consists mainly of water molecules with only two strong hydrogen bonds. Since this result is controversial, it is important to demonstrate the reliability of the experimental data, which is the purpose of this paper. Here we compare X-ray absorption spectra of liquid water recorded with five very different techniques sensitive to the local environment of the absorbing molecule. Overall, the spectra obtained with photon detection show a very close similarity and even the observable minor differences can be understood. The comparison demonstrates that XAS and XRS can indeed be applied reliably to study the local bonding of the water molecule and thus to reveal the hydrogen bond situation in bulk water.

The effects of confinement on the behavior of water molecules between parallel Au plates of (001) planes

Molecular dynamics simulation is utilized to investigate the behavior of water molecules confined between two Au plates of (001) planes separated by gaps of 24.48, 16.32, 12.24, 11.22, and 10.20 A. The simulation results indicate that the arrangements of the water molecules are dependent on the gap size. For the largest gap size, adsorption of the Au surface creates two permanent water layers in the vicinity of each Au plate. Furthermore, in this case, the gap size is sufficiently large to permit the formation of a central region within which the water molecules are randomly oriented in a similar manner to bulk water molecules. The results indicate that the orientation of the first water layer directly absorbed by the plate surface does not change as the gap size between the two Au plates is reduced. However, the orientations of the O-H bonds in the second water layer parallel to the surface rearrange to form hydrogen bonds between the water layers as the separation between the plates is decreased. Finally, an inspection of the variation of the self-diffusion coefficients with the gap size suggests that the difference between the dynamic properties of the water molecules in the z direction and the x-y plane decreases as the distance between the two Au plates increases.

Water motion in reverse micelles studied by quasielastic neutron scattering and molecular dynamics simulations

Motion of water molecules in Aerosol OT [sodium bis(2-ethylhexyl) sulfosuccinate, AOT] reverse micelles with water content w(0) ranging from 1 to 5 has been explored both experimentally through quasielastic neutron scattering (QENS) and with molecular dynamics (MD) simulations. The experiments were performed at the energy resolution of 85 microeV over the momentum transfer (Q) range of 0.36-2.53 A(-1) on samples in which the nonpolar phase (isooctane) and the AOT alkyl chains were deuterated, thereby suppressing their contribution to the QENS signal. QENS results were analyzed via a jump-diffusion/isotropic rotation model, which fits the results reasonably well despite the fact that confinement effects are not explicitly taken into account. This analysis indicates that in reverse micelles with low-water content (w(0)=1 and 2.5) translational diffusion rate is too slow to be detected, while for w(0)=5 the diffusion coefficient is much smaller than for bulk water. Rotational diffusion coefficients obtained from this analysis increase with w(0) and are smaller than for bulk water, but rotational mobility is less drastically reduced than translational mobility. Using the Faeder/Ladanyi model [J. Phys. Chem. B 104, 1033 (2000)] of reverse micelle interior, MD simulations were performed to calculate the self-intermediate scattering function F(S)(Q,t) for water hydrogens. Comparison of the time Fourier transform of this F(S)(Q,t) with the QENS dynamic structure factor S(Q,omega), shows good agreement between the model and experiment. Separate intermediate scattering functions F(S) (R)(Q,t) and F(S) (CM)(Q,t) were determined for rotational and translational motion. Consistent with the decoupling approximation used in the analysis of QENS data, the product of F(S) (R)(Q,t) and F(S) (CM)(Q,t) is a good approximation to the total F(S)(Q,t). We find that the decay of F(S) (CM)(Q,t) is nonexponential and our analysis of the MD data indicates that this behavior is due to lower water mobility close to the interface and to confinement-induced restrictions on the range of translational displacements. Rotational relaxation also exhibits nonexponential decay. However, rotational mobility of O-H bond vectors in the interfacial region remains fairly high due to the lower density of water-water hydrogen bonds in the vicinity of the interface.

Direct observation of terahertz surface modes in nanometer-sized liquid water pools

The far-infrared absorption spectrum of nanometer-sized water pools at the core of AOT micelles exhibits a pronounced resonance which is absent in bulk water. The amplitude and spectral position of this resonance are sensitive to the size of the confined water core. This resonance results from size-dependent modifications in the vibrational density of states, and thus has far-reaching implications for chemical processes which involve water sequestered within small cavities. These data represent the first study of the terahertz dielectric properties of confined liquids.