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.

Treatment of nontraditional source waters (e.g., produced water, municipal and industrial wastewaters, agricultural runoff) offers exciting opportunities to expand water and energy resources via water reuse and resource recovery. While conventional polymer membranes perform water/ion separations well, they do not provide solute-specific separation, a key component for these treatment opportunities. Herein, we discuss the selectivity limitations plaguing all conventional membranes, which include poor removal of small, neutral solutes and insufficient discrimination between ions of the same valence. Moreover, we present synthetic approaches for solute-tailored selectivity including the incorporation of single-digit nanopores and solute-selective ligands into membranes. Recent progress in these areas highlights the need for fundamental studies to rationally design membranes with selective moieties achieving desired separations.

We introduce a powerful, widely applicable approach to characterizing polymer conformational distributions, specifically the end-to-end distance distributions, P(R ee ), accessed through double electron−electron resonance (DEER) spectroscopy in conjunction with molecular dynamics (MD) simulations. The technique is demonstrated on one of the most widely used synthetic, disordered, water-soluble polymers: poly(ethylene oxide) (PEO). Despite its widespread importance, no systematic experimental characterization of PEO's R ee conformational landscape exists. The evaluation of P(R ee ) is particularly important for short polymers or (bio)polymers with sequence complexities that deviate from simple polymer physics scaling laws valid for long chains. In this study, we characterize the R ee landscape by measuring P(R ee ) for low molecular weight (MW: 0.22−2.6 kDa) dilute PEO chains. We use DEER with end-conjugated spin probes to resolve R ee populations from ∼2−9 nm and compare them with full distributions from MD. The P( R ee )'s from DEER and MD show remarkably good agreement, particularly at longer chain lengths where populations in the DEER-unresolvable range (<1.5 nm) are low. Both the P(R ee ) and the root-mean-square R ̃ee indicate that aqueous PEO is a semiflexible polymer in a good solvent, with the latter scaling linearly with molecular weight up to its persistence length (l p ∼ 0.48 nm), and rapidly transitioning to excluded volume scaling above l p . The R ̃ee scaling is quantitatively consistent with that from experimental scattering data on high MW (>10 kDa) PEO and the P(R ee )'s crossover to the theoretical distribution for an excluded volume chain.

Performance of membranes for water purification is highly influenced by the interactions of solvated species with membrane surfaces, including surface adsorption of solutes upon fouling. Current efforts toward fouling-resistant membranes often pursue surface hydrophilization, frequently motivated by macroscopic measures of hydrophilicity, because hydrophobicity is thought to increase solute–surface affinity. While this heuristic has driven diverse membrane functionalization strategies, here we build on advances in the theory of hydrophobicity to critically examine the relevance of macroscopic characterizations of solute–surface affinity. Specifically, we use molecular simulations to quantify the affinities to model hydroxyl- and methyl-functionalized surfaces of small, chemically diverse, charge-neutral solutes represented in produced water. We show that surface affinities correlate poorly with two conventional measures of solute hydrophobicity, gas-phase water solubility and oil–water partitioning. Moreover, we find that all solutes show attraction to the hydrophobic surface and most to the hydrophilic one, in contrast to macroscopically based hydrophobicity heuristics. We explain these results by decomposing affinities into direct solute interaction energies (which dominate on hydroxyl surfaces) and water restructuring penalties (which dominate on methyl surfaces). Finally, we use an inverse design algorithm to show how heterogeneous surfaces, with multiple functional groups, can be patterned to manipulate solute affinity and selectivity. These findings, importantly based on a range of solute and surface chemistries, illustrate that conventional macroscopic hydrophobicity metrics can fail to predict solute–surface affinity, and that molecular-scale surface chemical patterning significantly influences affinity—suggesting design opportunities for water purification membranes and other engineered interfaces involving aqueous solute–surface interactions.

We use molecular dynamics simulations to investigate the effect of polypeptoid sequence on the structure and dynamics of its hydration waters. Polypeptoids provide an excellent platform to study small-molecule hydration in disordered polymers, as they can be precisely synthesized with a variety of sidechain chemistries. We examine water behavior near a set of peptoid oligomers in which the number and placement of nonpolar versus polar sidechains are systematically varied.To do this, we leverage a new computational workflow enabling accurate sampling of polypeptoid conformations. We find that the hydration waters are less dense, are more tetrahedral, and have slower dynamics compared to bulk water. The magnitude of these shifts increases with the number of nonpolar groups. We also find that shifts in the water structure and dynamics are strongly correlated, suggesting that experimental insight into the dynamics of hydration water obtained by Overhauser dynamic nuclear polarization (ODNP) also contains information about water structural properties. We then demonstrate the ability of ODNP to probe site-specific dynamics of hydration water near these model peptoid systems.