Water Structuring Induces Nonuniversal Hydration Repulsion between Polar Surfaces: Quantitative Comparison between Molecular Simulations, Theory, and Experiments.

Polar surfaces in water typically repel each other at close separations, even if they are charge-neutral. This so-called hydration repulsion balances the van der Waals attraction and gives rise to a stable nanometric water layer between the polar surfaces. The resulting hydration water layer is crucial for the properties of concentrated suspensions of lipid membranes and hydrophilic particles in biology and technology, but its origin is unclear. It has been suggested that surface-induced molecular water structuring is responsible for the hydration repulsion, but a quantitative proof of this water-structuring hypothesis is missing. To gain an understanding of the mechanism causing hydration repulsion, we perform molecular simulations of different planar polar surfaces in water. Our simulated hydration forces between phospholipid bilayers agree perfectly with experiments, validating the simulation model and methods. For the comparison with theory, it is important to split the simulated total surface interaction force into a direct contribution from surface-surface molecular interactions and an indirect water-mediated contribution. We find the indirect hydration force and the structural water-ordering profiles from the simulations to be in perfect agreement with the predictions from theoretical models that account for the surface-induced water ordering, which strongly supports the water-structuring hypothesis for the hydration force. However, the comparison between the simulations for polar surfaces with different headgroup architectures reveals significantly different decay lengths of the indirect water-mediated hydration-force, which for laterally homogeneous water structuring would imply different bulk-water properties. We conclude that laterally inhomogeneous water ordering, induced by laterally inhomogeneous surface structures, shapes the hydration repulsion between polar surfaces in a decisive manner. Thus, the indirect water-mediated part of the hydration repulsion is caused by surface-induced water structuring but is surface-specific and thus nonuniversal.

Disentanglement of Surface and Confinement Effects for Diene Metathesis in Mesoporous Confinement.

We study the effects of a planar interface and confinement on a generic catalytically activated ring-closing polymerization reaction near an unstructured catalyst. For this, we employ a coarse-grained polymer model using grand-canonical molecular dynamics simulations with a Monte Carlo reaction scheme. Inspired by recent experiments in the group of M. Buchmeiser that demonstrated an increase in ring-closing selectivity under confinement, we show that both the interface effects, i.e., placing the catalyst near a planar wall, and the confinement effects, i.e., locating the catalyst within a pore, lead to an increase of selectivity. We furthermore demonstrate that curvature effects for cylindrical mesopores (2 nm < d < 12.3 nm) influence the distribution of the chain ends, leading to a further increase in selectivity. This leads us to speculate that specially corrugated surfaces might also help to enhance catalytically activated polymerization processes.

Lipid bilayer properties potentially contributed to the evolutionary disappearance of betaine lipids in seed plants.

BACKGROUND: Many organisms rely on mineral nutrients taken directly from the soil or aquatic environment, and therefore, developed mechanisms to cope with the limitation of a given essential nutrient. For example, photosynthetic cells have well-defined responses to phosphate limitation, including the replacement of cellular membrane phospholipids with non-phosphorous lipids. Under phosphate starvation, phospholipids in extraplastidial membranes are replaced by betaine lipids in microalgae. In higher plants, the synthesis of betaine lipid is lost, driving plants to other strategies to cope with phosphate starvation where they replace their phospholipids by glycolipids. RESULTS: The aim of this work was to evaluate to what extent betaine lipids and PC lipids share physicochemical properties and could substitute for each other. By neutron diffraction experiments and dynamic molecular simulation of two synthetic lipids, the dipalmitoylphosphatidylcholine (DPPC) and the dipalmitoyl-diacylglyceryl-N,N,N-trimethylhomoserine (DP-DGTS), we found that DP-DGTS bilayers are thicker than DPPC bilayers and therefore are more rigid. Furthermore, DP-DGTS bilayers are more repulsive, especially at long range, maybe due to unexpected unscreened electrostatic contribution. Finally, DP-DGTS bilayers could coexist in the gel and fluid phases. CONCLUSION: The different properties and hydration responses of PC and DGTS provide an explanation for the diversity of betaine lipids observed in marine organisms and for their disappearance in seed plants.

Renormalized charge and dielectric effects in colloidal interactions: a numerical solution of the nonlinear Poisson-Boltzmann equation for unknown boundary conditions.

 The Derjaguin-Landau-Verwey-Overbeek (DLVO) theory, introduced more than 70 years ago, is a hallmark of colloidal particle modeling. For highly charged particles in the dilute regime, it is often supplemented by Alexander's prescription (Alexander et al. in J Chem Phys 80:5776, 1984) for using a renormalized charge. Here, we solve the problem of the interaction between two charged colloids at finite ionic strength, including dielectric mismatch effects, using an efficient numerical scheme to solve the nonlinear Poisson-Boltzmann (NPB) equation with unknown boundary conditions. Our results perfectly match the analytical predictions for the renormalized charge by Trizac and coworkers (Aubouy et al. in J Phys A 36:5835, 2003). Moreover, they allow us to reinterpret previous molecular dynamics (MD) simulation results by Kreer et al. (Phys Rev E 74:021401, 2006), rendering them now in agreement with the expected behavior. We furthermore find that the influence of polarization becomes important only when the Debye layers overlap significantly.

Solvent Effects on Structure and Screening in Confined Electrolytes.

Using classical density functional theory, we investigate the influence of solvent on the structure and ionic screening of electrolytes under slit confinement and in contact with a reservoir. We consider a symmetric electrolyte with implicit and explicit solvent models and find that spatially resolving solvent molecules is essential for the ion structure at confining walls, excess ion adsorption, and the pressure exerted on the walls. Despite this, we observe only moderate differences in the period of oscillations of the pressure with the slit width and virtually coinciding decay lengths as functions of the scaling variable σ_{ion}/λ_{D}, where σ_{ion} is the ion diameter and λ_{D} the Debye length. Moreover, in the electrostatic-dominated regime, this scaling behavior is practically independent of the relative permittivity and its dependence on the ion concentration. In contrast, the crossover to the hard-core-dominated regime depends sensitively on all three factors.

A screening of results on the decay length in concentrated electrolytes.

Screening of electrostatic interactions in room-temperature ionic liquids and concentrated electrolytes has recently attracted much attention as surface force balance experiments have suggested the emergence of unanticipated anomalously large screening lengths at high ion concentrations. Termed underscreening, this effect was ascribed to the bulk properties of concentrated ionic systems. However, underscreening under experimentally relevant conditions is not predicted by classical theories and challenges our understanding of electrostatic correlations. Despite the enormous effort in performing large-scale simulations and new theoretical investigations, the origin of the anomalously long-range screening length remains elusive. This contribution briefly summarises the experimental, analytical and simulation results on ionic screening and the scaling behaviour of screening lengths. We then present an atomistic simulation approach that accounts for the solvent and ion exchange with a reservoir. We find that classical density functional theory (DFT) for concentrated electrolytes under confinement reproduces ion adsorption at charged interfaces surprisingly well. With DFT, we study confined electrolytes using implicit and explicit solvent models and the dependence on the solvent's dielectric properties. Our results demonstrate how the absence vs. presence of solvent particles and their discrete nature affect the short and long-range screening in concentrated ionic systems.

Assessing the Validity of NMR Relaxation Rates Obtained from Coarse-Grained Simulations of PEG-Water Mixtures.

NMR relaxometry is a powerful and well-established experimental approach for characterizing dynamic processes in soft matter systems. All-atom (AA) resolved simulations are typically employed to gain further microscopic insights while reproducing the relaxation rates R1. However, such approaches are limited to time and length scales that prevent to model systems such as long polymer chains or hydrogels. Coarse graining (CG) can overcome this barrier at the cost of losing atomistic details that impede the calculation of NMR relaxation rates. Here, we address this issue by performing a systematic characterization of dipolar relaxation rates R1 on a PEG-H2O mixture at two different levels of details: AA and CG. Remarkably, we show that NMR relaxation rates R1 obtained at the CG level obey the same trends when compared to AA calculations but with a systematic offset. This offset is due to, on the one hand, the lack of an intramonomer component and, on the other hand, the inexact positioning of the spin carriers. We show that the offset can be corrected for quantitatively by reconstructing a posteriori the atomistic details for the CG trajectories.

The Three-Phase Contact Potential Difference Modulates the Water Surface Charge.

The surface charge of an open water surface is crucial for solvation phenomena and interfacial processes in aqueous systems. However, the magnitude of the charge is controversial, and the physical mechanism of charging remains incompletely understood. Here we identify a previously overlooked physical mechanism determining the surface charge of water. Using accurate charge measurements of water microdrops, we demonstrate that the water surface charge originates from the electrostatic effects in the contact line vicinity of three phases, one of which is water. Our experiments, theory, and simulations provide evidence that a junction of two aqueous interfaces (e.g., liquid-solid and liquid-air) develops a pH-dependent contact potential difference Δϕ due to the longitudinal charge redistribution between two contacting interfaces. This universal static charging mechanism may have implications for the origin of electrical potentials in biological, nanofluidic, and electrochemical systems and helps to predict and control the surface charge of water in various experimental environments.

NMR Investigation of Water in Salt Crusts: Insights from Experiments and Molecular Simulations.

The evaporation of water from bare soil is often accompanied by the formation of a layer of crystallized salt, a process that must be understood in order to address the issue of soil salinization. Here, we use nuclear magnetic relaxation dispersion measurements to better understand the dynamic properties of water within two types of salt crusts: sodium chloride (NaCl) and sodium sulfate (Na2SO4). Our experimental results display a stronger dispersion of the relaxation time T1 with frequency for the case of sodium sulfate as compared to sodium chloride salt crusts. To gain insight into these results, we perform molecular dynamics simulations of salt solutions confined within slit nanopores made of either NaCl or Na2SO4. We find a strong dependence of the value of the relaxation time T1 on pore size and salt concentration. Our simulations reveal the complex interplay between the adsorption of ions at the solid surface, the structure of water near the interface, and the dispersion of T1 at low frequency, which we attribute to adsorption-desorption events.

Transport of thin water films: From thermally activated random walks to hydrodynamics.

Under ambient atmospheric conditions, a thin film of water wets many solid surfaces, including insulators, ice, and salt. The film thickness as well as its transport behavior sensitively depend on the surrounding humidity. Understanding this intricate interplay is of the highest relevance for water transport through porous media, particularly in the context of soil salinization induced by evaporation. Here, we use molecular simulations to evaluate the transport properties of thin water films on prototypical salt and soil interfaces, namely NaCl and silica solid surfaces. Our results show two distinct regimes for water transport: at low water coverage, the film permeance scales linearly with the adsorbed amount, in agreement with the activated random walk model. For thicker water films, the permeance scales as the adsorbed amount to the power of 3, in line with the Stokes equation. By comparing results obtained for silica and NaCl surfaces, we find that, at low water coverage, water permeance at the silica surface is considerably lower than at the NaCl surface, which we attribute to difference in hydrogen bonding. We also investigate the effect of atomic surface defects on the transport properties. Finally, in the context of water transport through the porous material, we determine the humidity-dependent crossover between a vapor-dominated and a thin film-dominated transport regimes depending on the pore size.

Dielectric constant of aqueous solutions of proteins and organic polymers from molecular dynamics simulations.

The dielectric constant of water/oligomer mixtures, spanning the range from pure water to pure oligomeric melts, is investigated using molecular dynamics (MD) simulations. As prototypical water-soluble organic substances, we consider neutral poly-glycine, poly-ethylene glycol, and charged monomeric propionic acid. As the water content is reduced, the dielectric constant decreases but does not follow an ideal mixing behavior. The deviations from ideal mixing originate primarily in the non-linear relation between the oligomer mass fraction and collective polarization effects. We find that the dielectric constant is dominated by water polarization, even if the oligomer mass fraction exceeds 50%. By a double extrapolation of the MD simulation results to the limit of vanishing water fraction and to the limit of infinite oligomeric chain length, we estimate the orientational contribution to the dielectric constant of the pure polymeric melts. By this procedure, we obtain ɛ = 17 ± 2 for polyglycine and ɛ = 1 ± 0.3 for polyethylene glycol. The large difference is rationalized by polarization correlations of glycine units. Interestingly, we find constant temperature simulations to outperform replica exchange simulations in terms of equilibration speed.

Electronic screening using a virtual Thomas-Fermi fluid for predicting wetting and phase transitions of ionic liquids at metal surfaces.

Of relevance to energy storage, electrochemistry and catalysis, ionic and dipolar liquids display unexpected behaviours-especially in confinement. Beyond adsorption, over-screening and crowding effects, experiments have highlighted novel phenomena, such as unconventional screening and the impact of the electronic nature-metallic versus insulating-of the confining surface. Such behaviours, which challenge existing frameworks, highlight the need for tools to fully embrace the properties of confined liquids. Here we introduce a novel approach that involves electronic screening while capturing molecular aspects of interfacial fluids. Although available strategies consider perfect metal or insulator surfaces, we build on the Thomas-Fermi formalism to develop an effective approach that deals with any imperfect metal between these asymptotes. Our approach describes electrostatic interactions within the metal through a 'virtual' Thomas-Fermi fluid of charged particles, whose Debye length sets the screening length λ. We show that this method captures the electrostatic interaction decay and electrochemical behaviour on varying λ. By applying this strategy to an ionic liquid, we unveil a wetting transition on switching from insulating to metallic conditions.

Water at charged interfaces.

The ubiquity of aqueous solutions in contact with charged surfaces and the realization that the molecular-level details of water-surface interactions often determine interfacial functions and properties relevant in many natural processes have led to intensive research. Even so, many open questions remain regarding the molecular picture of the interfacial organization and preferential alignment of water molecules, as well as the structure of water molecules and ion distributions at different charged interfaces. While water, solutes and charge are present in each of these systems, the substrate can range from living tissues to metals. This diversity in substrates has led to different communities considering each of these types of aqueous interface. In this Review, by considering water in contact with metals, oxides and biomembranes, we show the essential similarity of these disparate systems. While in each case the classical mean-field theories can explain many macroscopic and mesoscopic observations, it soon becomes apparent that such theories fail to explain phenomena for which molecular properties are relevant, such as interfacial chemical conversion. We highlight the current knowledge and limitations in our understanding and end with a view towards future opportunities in the field.

Universal and Nonuniversal Aspects of Electrostatics in Aqueous Nanoconfinement.

Dielectric water properties, which significantly change in confinement, determine electrostatic interactions and thereby influence all molecular forces and chemical reactions. We present comparative simulations of water between graphene sheets, decanol monolayers, and phospholipid and glycolipid bilayers. Generally, dielectric profiles strongly differ in perpendicular and parallel surface directions and for large surface separation decay to the bulk value 1-2 nm away from the surface. Polar surface groups enhance the local interfacial dielectric response and for phospholipid bilayers induce a giant parallel contribution. A mapping on a box model with asymptotically determined effective water layer widths demonstrates that the perpendicular effective dielectric constant for all systems decreases for confinement below a nanometer, while the parallel one stays rather constant. The confinement-dependent perpendicular effective dielectric constant for graphene is in agreement with experimental data only if the effective water layer width is suitably adjusted. The interactions between two charges at small separation depend on the product of parallel and perpendicular effective water dielectric components; for large separation the interactions depend on the confining medium. For metallic confining media the interactions at large separation decay exponentially with a decay length that depends on the ratio of the effective parallel and perpendicular water dielectric components.

Giant Axial Dielectric Response in Water-Filled Nanotubes and Effective Electrostatic Ion-Ion Interactions from a Tensorial Dielectric Model.

Molecular dynamics simulations in conjunction with effective medium theory are used to investigate dielectric effects in water-filled nanotubes. The resulting effective axial dielectric constant shows a divergent increase for small nanotube radii that depends on the nanotube length, while the effective radial dielectric constant decreases significantly for thin nanotubes. By solving Poisson's equation for an anisotropic dielectric medium in cylindrical geometry, we show that the axial ion-ion interaction depends for small separations primarily on the radial dielectric constant, not on the axial one. This means that electrostatic ion-ion interactions in thin water-filled nanotubes are on the linear dielectric level significantly enhanced due to water confinement effects at small separations, while at large separations the outside medium dominates. If the outside medium is metallic, then the ion-ion interaction decays exponentially for large ion separation.

Memory-kernel extraction for different molecular solutes in solvents of varying viscosity in confinement.

The friction coefficient of molecular solutes depends on the solute, on the solvent, and on the solute-solvent interactions, but is typically assumed to not depend on an externally applied force that acts on the solute. In this paper we compute the friction memory function from molecular dynamics simulations and show that the friction coefficients of harmonically confined methane, water, Na^{+}, an artificial Na^{-} ion, and glycerol in water in fact increase with confinement strength. The results show that the friction increase with confinement strength is a fundamental effect that occurs for hydrophobic, hydrophilic, as well as charged molecules. We demonstrate that a parameter-free extraction of the running integral over the memory function yields the most robust results when compared to methods based on parametrization or Fourier transforms. In all systems, this friction increase is accompanied by a slowdown of the solvent dynamics in the first hydration shell of the solutes. By simulations of a confined glycerol molecule in water-glycerol mixtures, we furthermore demonstrate that the friction dependence on the confining potential is magnified in more viscous solvents, which suggests that this effect plays an important role for larger molecules in highly viscous solutions like polymer melts, in line with dynamic scaling arguments.

Influence of polar co-solutes and salt on the hydration of lipid membranes.

The influence of the co-solutes TMAO, urea, and NaCl on the hydration repulsion between lipid membranes is investigated in a combined experimental/simulation approach. Pressure-hydration curves obtained via sorption experiments reveal that the repulsion significantly increases when the membranes are loaded with co-solutes, most strongly for TMAO. As a result, the co-solutes retain additional water molecules and therefore provide membranes with a fluid and more physiological environment. The experimental data are quantitatively reproduced in complementary solvent-explicit atomistic molecular dynamics simulations, which yield the chemical potential of water. Simulation analysis reveals that the additional repulsion arises from the osmotic pressure generated by the co-solutes, an effect which is maximal for TMAO, due to its unfavorable interactions with the lipid headgroup layer and its extraordinarily high osmotic coefficient.