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.

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.

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.

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.

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.

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.

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.

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.

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.

Simulations of Nanoseparated Charged Surfaces Reveal Charge-Induced Water Reorientation and Nonadditivity of Hydration and Mean-Field Electrostatic Repulsion.

We perform atomistic simulations of nanometer-separated charged surfaces in the presence of monovalent counterions at fixed water chemical potential. The counterion density profiles are well described by a modified Poisson-Boltzmann (MPB) approach that accounts for nonelectrostatic ion-surface interactions, while the effects of smeared-out surface-charge distributions and dielectric profiles are found to be relatively unimportant. The simulated surface interactions are for weakly charged surfaces well described by the additive contributions of hydration and MPB repulsions, but already for a moderate surface charge density of σ = -0.77 e/nm2 this additivity breaks down. This we rationalize by a combination of different effects, namely, counterion correlations as well as the surface charge-induced reorientation of hydration water, which modifies the effective water dielectric constant as well as the hydration repulsion.

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.

Dispersion truncation affects the phase behavior of bulk and confined fluids: Coexistence, adsorption, and criticality.

We present molecular simulations of bulk and confined Lennard-Jones fluids to assess the effect of dispersion truncation through a simple spherical cutoff. The latter is well corrected on a mean field level for bulk fluids if the cutoff distance is larger than about three molecular diameters. In confinement, however, there is no general analytical treatment, and thus, the truncated and shifted Lennard-Jones potential has to be employed, with drastic consequences on the bulk critical temperature, vapor/liquid coexistence pressure, and surface tension. We show using grand-canonical Monte-Carlo simulations of nitrogen adsorption in amorphous silica nanopores that the choice of the cutoff significantly modifies the pressure at which capillary condensation occurs and compute the capillary critical temperature in terms of a first order transition between an adsorbed film and filled pores.

Hydration Friction in Nanoconfinement: From Bulk via Interfacial to Dry Friction.

The viscous properties of nanoscopically confined water are important when hydrated surfaces in close contact are sheared against each other. Numerous experiments have probed the friction between atomically flat hydrated surfaces in the subnanometer separation regime and suggested an increased water viscosity, but the value of the effective viscosity of ultraconfined water, the mechanism of hydration layer friction, and the crossover to the dry friction limit are unclear. We study the shear friction between polar surfaces by extensive nonequilibrium molecular dynamics simulations in the linear-response regime at low shearing velocity, which is the relevant regime for typical biological applications. With decreasing water film thickness we find three consecutive friction regimes: For thick films friction is governed by bulk water viscosity. At separations of about a nanometer the highly viscous interfacial water layers dominate and increase the surface friction, while at the transition to the dry friction limit interfacial slip sets in. Based on our simulation results, we construct a confinement-dependent friction model which accounts for the additive friction contributions from bulklike water, interfacial water layers, and interfacial slip and which is valid for arbitrary water film thickness.

Water Dielectric Effects in Planar Confinement.

We investigate the dielectric profile of water confined between two planar polar walls using atomistic molecular dynamics simulations. For a water slab thickness below 1 nm the dielectric response is highly asymmetric: while the parallel component slightly increases compared to bulk, the perpendicular one decreases drastically due to anticorrelated polarization of neighboring water molecules. We demonstrate the importance of the dielectric contribution due to flexible polar headgroups and derive an effective dielectric tensorial box model suitable for coarse-grained electrostatic modeling.

Water-Mediated Interactions between Hydrophilic and Hydrophobic Surfaces.

All surfaces in water experience at short separations hydration repulsion or hydrophobic attraction, depending on the surface polarity. These interactions dominate the more long-ranged electrostatic and van der Waals interactions and are ubiquitous in biological and colloidal systems. Despite their importance in all scenarios where the surface separation is in the nanometer range, the origin of these hydration interactions is still unclear. Using atomistic solvent-explicit molecular dynamics simulations, we analyze the interaction free energies of charge-neutral model surfaces with different elastic and water-binding properties. The surface polarity is shown to be the most important parameter that not only determines the hydration properties and thereby the water contact angle of a single surface but also the surface-surface interaction and whether two surfaces attract or repel. Elastic properties of the surfaces are less important. On the basis of surface contact angles and surface-surface binding affinities, we construct a universal interaction diagram featuring three different interaction regimes-hydration repulsion, cavitation-induced attraction-and for intermediate surface polarities-dry adhesion. On the basis of scaling arguments and perturbation theory, we establish simple combination rules that predict the interaction behavior for combinations of dissimilar surfaces.

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.

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.