Simulations of water nano-confined between corrugated planes

Water confined to nanoscale widths in two dimensions between ideal planar walls has been the subject of ample study, aiming at understanding the intrinsic response of water to confinement, avoiding the consideration of the chemistry of actual confining materials. In this work, we study the response of such nanoconfined water to the imposition of a periodicity in the confinement by means of computer simulations, both using empirical potentials and from first-principles. For that we propose a periodic confining potential emulating the atomistic oscillation of the confining walls, which allows varying the lattice parameter and amplitude of the oscillation. We do it for a triangular lattice, with several values of the lattice parameter: one which is ideal for commensuration with layers of Ih ice and other values that would correspond to more realistic substrates. For the former, the phase diagram shows an overall rise of the melting temperature. The liquid maintains a bi-layer triangular structure, however, despite the fact that it is not favoured by the external periodicity. The first-principles liquid is significantly affected by the modulation in its layering and stacking even at relatively small amplitudes of the confinement modulation. Beyond some critical modulation amplitude, the hexatic phase present in flat confinement is replaced by a trilayer crystalline phase unlike any of the phases encountered for flat confinement. For more realistic lattice parameters, the liquid does not display higher tendency to freeze, but it clearly shows inhomogeneous behaviour as the strength of the rugosity increases. In spite of this expected inhomogeneity, the structural and dynamical response of the liquid is surprisingly insensitive to the external modulation. Although the first-principles calculations give a more triangular liquid than the one observed with empirical potentials (TIP4P/2005), both agree remarkably well for the main conclusions of the study.

Studies on electrostatic interactions within model nano-confined aqueous environments of different chemical nature

We study the potential of mean force for pairs of parallel flat surfaces with attractive electrostatic interactions by employing model systems functionalized with different charged, hydrophobic and hydrophilic groups. We study the way in which the local environment (hydrophobic or hydrophilic moieties) modulates the interaction between the attractive charged groups on the plates by removing or attracting nearby water and thus screening or not the electrostatic interaction. To explicitly account for the role of the solvent and the local hydrophobicity, we also perform studies in vacuo. Additionally, the results are compared to that for non-charged plates in order to single out and rationalize the non-additivity of the different non-covalent interactions. Our simulations demonstrate that the presence of neighboring hydrophobic groups promote water removal in the vicinity of the charged groups, thus enhancing charge attraction upon self-assembly. This role of the local hydrophobicity modulating electrostatic interactions is consistent with recent qualitative descriptions in the protein binding context.

Effect of material flexibility on the thermodynamics and kinetics of hydrophobically induced evaporation of water

The evaporation of water induced by confinement between hydrophobic surfaces has received much attention due to its suggested functional role in numerous biophysical phenomena and its importance as a general mechanism of hydrophobic self-assembly. Although much progress has been made in understanding the basic physics of hydrophobically induced evaporation, a comprehensive understanding of the substrate material features (e.g., geometry, chemistry, and mechanical properties) that promote or inhibit such transitions remains lacking. In particular, comparatively little research has explored the relationship between water's phase behavior in hydrophobic confinement and the mechanical properties of the confining material. Here, we report the results of extensive molecular simulations characterizing the rates, free energy barriers, and mechanism of water evaporation when confined between model hydrophobic materials with tunable flexibility. A single-order-of-magnitude reduction in the material's modulus results in up to a nine-orders-of-magnitude increase in the evaporation rate, with the corresponding characteristic time decreasing from tens of seconds to tens of nanoseconds. Such a modulus reduction results in a 24-orders-of-magnitude decrease in the reverse rate of condensation, with time scales increasing from nanoseconds to tens of millions of years. Free energy calculations provide the barriers to evaporation and confirm our previous theoretical predictions that making the material more flexible stabilizes the confined vapor with respect to liquid. The mechanism of evaporation involves surface bubbles growing/coalescing to form a subcritical gap-spanning tube, which then must grow to cross the barrier.

Surface dissimilarity affects critical distance of influence for confined water

In this study, the properties of nano-confined water, such as density, orientation etc., are monitored across varying confinement spacing to determine the critical distance of influence between dissimilar surfaces.

A nonmonotonic dependence of the contact angles on the surface polarity for a model solid surface

We found an unusual nonmonotonic contact angle dependence of the surface polarity (denoted as q) on a solid surface with specific charge patterns, where the contact angle firstly decreases and then increases as q increases from 0 e to 1.0 e.

Anomalous Dynamics of Water Confined in Protein-Protein and Protein-DNA Interfaces

Confined water often exhibits anomalous properties not observable in the bulk phase. Although water in hydrophobic confinement has been the focus of intense investigation, the behavior of water confined between hydrophilic surfaces, which are more frequently found in biological systems, has not been fully explored. Here, we investigate using molecular dynamics simulations dynamical properties of the water confined in hydrophilic protein-protein and protein-DNA interfaces. We find that the interfacial water exhibits glassy slow relaxations even at 300 K. In particular, the rotational dynamics show a logarithmic decay that was observed in glass-forming liquids at deeply supercooled states. We argue that such slow water dynamics are indeed induced by the hydrophilic binding surfaces, which is in opposition to the picture that the hydration water slaves protein motions. Our results will significantly impact the view on the role of water in biomolecular interactions.

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.

Water Determines the Structure and Dynamics of Proteins

Water is an essential participant in the stability, structure, dynamics, and function of proteins and other biomolecules. Thermodynamically, changes in the aqueous environment affect the stability of biomolecules. Structurally, water participates chemically in the catalytic function of proteins and nucleic acids and physically in the collapse of the protein chain during folding through hydrophobic collapse and mediates binding through the hydrogen bond in complex formation. Water is a partner that slaves the dynamics of proteins, and water interaction with proteins affect their dynamics. Here we provide a review of the experimental and computational advances over the past decade in understanding the role of water in the dynamics, structure, and function of proteins. We focus on the combination of X-ray and neutron crystallography, NMR, terahertz spectroscopy, mass spectroscopy, thermodynamics, and computer simulations to reveal how water assist proteins in their function. The recent advances in computer simulations and the enhanced sensitivity of experimental tools promise major advances in the understanding of protein dynamics, and water surely will be a protagonist.

Confined Water as Model of Supercooled Water

Water in confined geometries has obvious relevance in biology, geology, and other areas where the material properties are strongly dependent on the amount and behavior of water in these types of materials. Another reason to restrict the size of water domains by different types of geometrical confinements has been the possibility to study the structural and dynamical behavior of water in the deeply supercooled regime (e.g., 150-230 K at ambient pressure), where bulk water immediately crystallizes to ice. In this paper we give a short review of studies with this particular goal. However, from these studies it is also clear that the interpretations of the experimental data are far from evident. Therefore, we present three main interpretations to explain the experimental data, and we discuss their advantages and disadvantages. Unfortunately, none of the proposed scenarios is able to predict all the observations for supercooled and glassy bulk water, indicating that either the structural and dynamical alterations of confined water are too severe to make predictions for bulk water or the differences in how the studied water has been prepared (applied cooling rate, resulting density of the water, etc.) are too large for direct and quantitative comparisons.

Two Dimensional Ice from First Principles: Structures and Phase Transitions

Despite relevance to disparate areas such as cloud microphysics and tribology, major gaps in the understanding of the structures and phase transitions of low-dimensional water ice remain. Here, we report a first principles study of confined 2D ice as a function of pressure. We find that at ambient pressure hexagonal and pentagonal monolayer structures are the two lowest enthalpy phases identified. Upon mild compression, the pentagonal structure becomes the most stable and persists up to ∼2  GPa, at which point the square and rhombic phases are stable. The square phase agrees with recent experimental observations of square ice confined within graphene sheets. This work provides a fresh perspective on 2D confined ice, highlighting the sensitivity of the structures observed to both the confining pressure and the width.

Influence of disordered porous media on the anomalous properties of a simple water model

The thermodynamic, dynamic, and structural behavior of a water-like system confined in a matrix is analyzed for increasing confining geometries. The liquid is modeled by a two-dimensional associating lattice gas model that exhibits density and diffusion anomalies, similar to the anomalies present in liquid water. The matrix is a triangular lattice in which fixed obstacles impose restrictions to the occupation of the particles. We show that obstacles shorten all lines, including the phase coexistence, the critical and the anomalous lines. The inclusion of a very dense matrix not only suppresses the anomalies but also the liquid-liquid critical point.

Pathways to dewetting in hydrophobic confinement

Liquid water can become metastable with respect to its vapor in hydrophobic confinement. The resulting dewetting transitions are often impeded by large kinetic barriers. According to macroscopic theory, such barriers arise from the free energy required to nucleate a critical vapor tube that spans the region between two hydrophobic surfaces--tubes with smaller radii collapse, whereas larger ones grow to dry the entire confined region. Using extensive molecular simulations of water between two nanoscopic hydrophobic surfaces, in conjunction with advanced sampling techniques, here we show that for intersurface separations that thermodynamically favor dewetting, the barrier to dewetting does not correspond to the formation of a (classical) critical vapor tube. Instead, it corresponds to an abrupt transition from an isolated cavity adjacent to one of the confining surfaces to a gap-spanning vapor tube that is already larger than the critical vapor tube anticipated by macroscopic theory. Correspondingly, the barrier to dewetting is also smaller than the classical expectation. We show that the peculiar nature of water density fluctuations adjacent to extended hydrophobic surfaces--namely, the enhanced likelihood of observing low-density fluctuations relative to Gaussian statistics--facilitates this nonclassical behavior. By stabilizing isolated cavities relative to vapor tubes, enhanced water density fluctuations thus stabilize novel pathways, which circumvent the classical barriers and offer diminished resistance to dewetting. Our results thus suggest a key role for fluctuations in speeding up the kinetics of numerous phenomena ranging from Cassie-Wenzel transitions on superhydrophobic surfaces, to hydrophobically driven biomolecular folding and assembly.

The role of material flexibility on the drying transition of water between hydrophobic objects: a thermodynamic analysis

Liquid water confined between hydrophobic objects of sufficient size becomes metastable with respect to its vapor at separations smaller than a critical drying distance. Macroscopic thermodynamic arguments predicting this distance have been restricted to the limit of perfectly rigid confining materials. However, no material is perfectly rigid and it is of interest to account for this fact in the thermodynamic analysis. We present a theory that combines the current macroscopic theory with the thermodynamics of elasticity to derive an expression for the critical drying distance for liquids confined between flexible materials. The resulting expression is the sum of the well-known drying distance for perfectly rigid confining materials and a new term that accounts for flexibility. Thermodynamic arguments show that this new term is necessarily positive, meaning that flexibility increases the critical drying distance. To study the expected magnitude and scaling behavior of the flexible term, we consider the specific case of water and present an example of drying between thin square elastic plates that are simply supported along two opposite edges and free at the remaining two. We find that the flexible term can be the same order of magnitude or greater than the rigid solution for materials of biological interest at ambient conditions. In addition, we find that when the rigid solution scales with the characteristic size of the immersed objects, the flexible term is independent of size and vice versa. Thus, the scaling behavior of the overall drying distance will depend on the relative weights of the rigid and flexible contributions.

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.

Communication: Chemical functionality of interfacial water enveloping nanoscale structural defects in proteins

Building upon a non-Debye multiscale treatment of water dielectrics, this work reveals the biochemical role of interfacial water enveloping nanoscale structural defects in soluble proteins, asserting its role as a chemical base. This quasi-reactant status is already implied by the significant concentration of structural defects in the vicinity of an enzymatically active site, delineating their role as promoters or enhancers of catalytic activity.

Role of arginine in mediating protein-carbon nanotube interactions

Arginine-rich proteins (e.g., lysozyme) or poly-L-arginine peptides have been suggested as solvating and dispersing agents for single-wall carbon nanotubes (CNTs) in water. In addition, protein structure-function in porous and hydrophobic materials is of broad interest. The amino acid residue, arginine (Arg(+)), has been implicated as an important mediator of protein/peptide-CNT interactions. To understand the structural and thermodynamic aspects of this interaction at the molecular level, we employ molecular dynamics (MD) simulations of the protein lysozyme in the interior of a CNT, as well as of free solutions of Arg(+) in the presence of a CNT. To dissect the Arg(+)-CNT interaction further, we also perform simulations of aqueous solutions of the guanidinium ion (Gdm(+)) and the norvaline (Nva) residue in the presence of a CNT. We show that the interactions of lysozyme with the CNT are mediated by the surface Arg(+) residues. The strong interaction of Arg(+) residue with the CNT is primarily driven by the favorable interactions of the Gdm(+) group with the CNT wall. The Gdm(+) group is not as well-hydrated on its flat sides, which binds to the CNT wall. This is consistent with a similar binding of Gdm(+) ions to a hydrophobic polymer. In contrast, the Nva residue, which lacks the Gdm(+) group, binds to the CNT weakly. We present details of the free energy of binding, molecular structure, and dynamics of these solutes on the CNT surface. Our results highlight the important role of Arg(+) residues in protein-CNT or protein-carbon-based material interactions. Such interactions could be manipulated precisely through protein engineering, thereby offering control over protein orientation and structure on CNTs, graphene, or other hydrophobic interfaces.

Hydration and Nanoconfined Water: Insights from Computer Simulations

The comprehension of the structure and behavior of water at interfaces and under nanoconfinement represents an issue of major concern in several central research areas like hydration, reaction dynamics and biology. From one side, water is known to play a dominant role in the structuring, the dynamics and the functionality of biological molecules, governing main processes like protein folding, protein binding and biological function. In turn, the same principles that rule biological organization at the molecular level are also operative for materials science processes that take place within a water environment, being responsible for the self-assembly of molecular structures to create synthetic supramolecular nanometrically-sized materials. Thus, the understanding of the principles of water hydration, including the development of a theory of hydrophobicity at the nanoscale, is imperative both from a fundamental and an applied standpoint. In this work we present some molecular dynamics studies of the structure and dynamics of water at different interfaces or confinement conditions, ranging from simple model hydrophobic interfaces with different geometrical constraints (in order to single out curvature effects), to self-assembled monolayers, proteins and phospholipid membranes. The tendency of the water molecules to sacrifice the lowest hydrogen bond (HB) coordination as possible at extended interfaces is revealed. This fact makes the first hydration layers to be highly oriented, in some situations even resembling the structure of hexagonal ice. A similar trend to maximize the number of HBs is shown to hold in cavity filling, with small subnanometric hydrophobic cavities remaining empty while larger cavities display an alternation of filled and dry states with a significant inner HB network. We also study interfaces with complex chemical and geometrical nature in order to determine how different conditions affect the local hydration properties. Thus, we show some results for protein hydration and, particularly, some preliminary studies on membrane hydration. Finally, calculations of a local hydrophobicity measure of relevance for binding and self-assembly are also presented. We then conclude with a few words of further emphasis on the relevance of this kind of knowledge to biology and to the design of new materials by highlighting the context-dependent and non-additive nature of different non-covalent interactions in an aqueous nanoenvironment, an issue that is usually greatly overlooked.

Solvent cavitation under solvophobic confinement

The stability of liquids under solvophobic confinement can tip in favor of the vapor phase, nucleating a liquid-to-vapor phase transition that induces attractive forces between confining surfaces. In the case of water adjacent to hydrophobic surfaces, experimental and theoretical evidence support confinement-mediated evaporation stabilization of biomolecular and colloidal assemblies. The macroscopic thermodynamic theory of cavitation under confinement establishes the connection between the size of the confining surfaces, interfacial free energies, and bulk solvent pressure with the critical evaporation separation and interfacial forces. While molecular simulations have confirmed the broad theoretical trends, a quantitative comparison based on independent measurements of the interfacial free energies and liquid-vapor coexistence properties has, to the best of our knowledge, not yet been performed. To overcome the challenges of simulating a large number of systems to validate scaling predictions for a three-dimensional fluid, we simulate both the forces and liquid-vapor coexistence properties of a two-dimensional Lennard-Jones fluid confined between solvophobic plates over a range of plate sizes and reservoir pressures. Our simulations quantitatively agree with theoretical predictions for solvent-mediated forces and critical evaporation separations once the length dependence of the solvation free energy of an individual confining plate is taken into account. The effective solid-liquid line tension length dependence results from molecular scale correlations for solvating microscopic plates and asymptotically decays to the macroscopic value for plates longer than 150 solvent diameters. The success of the macroscopic thermodynamic theory at describing two-dimensional liquids suggests application to surfactant monolayers to experimentally confirm confinement-mediated cavitation.

Water Dynamics in Metal-Organic Frameworks: Effects of Heterogeneous Confinement Predicted by Computational Spectroscopy

The behavior of water confined in MIL-53(Cr), a flexible metal-organic framework (MOF), is investigated through computational infrared spectroscopy. As the number of molecules adsorbed inside of the pores increases, the water OH stretch band of the linear infrared spectrum grows in intensity and approaches that of bulk water. To assess whether the water confined in MIL-53(Cr) becomes liquid-like, two-dimensional infrared spectra (2DIR) are also calculated. Confinement effects result in distinct chemical environments that appear as specific features in the 2DIR spectra. The evolution of the 2DIR line shape as a function of waiting time is well described in terms of the orientational dynamics of the water molecules, with chemical exchange cross peaks appearing at a time scale similar to the hydrogen bond rearrangement lifetime. The confining environment considerably slows the hydrogen bond dynamics relative to bulk as a result of the competition between water-framework and water-water interactions.

Role of Charge Transfer in Water Diffusivity in Aqueous Ionic Solutions

We performed molecular dynamics simulations on four types of systems containing ion and solvating water. Two systems contained a cation (Na(+) or K(+)), and two other systems an anion (Cl(-) or I(-)). Classical molecular dynamics simulations were performed using three different force fields: a fixed charge force field, a polarizable force field that includes explicit polarization, and also a recently developed force field that includes polarization and charge transfer. These simulations were then compared to first-principles molecular dynamics simulations. While the first-principles simulations showed that the anions accelerated water translational diffusion, the cations slowed it down. In simulations with the classical force fields, only the force field that incorporates explicit charge transfer reproduced this ion-specific behavior. Additional simulations performed to understand the effect of charge transfer demonstrated that two competitive factors determine the behavior of water translational diffusion: the ions diminished charge accelerates water, while the net charge acquired by water either accelerates or slows down its dynamics. Our results show that charge transfer plays a crucial role in governing the water dynamics in aqueous ionic solutions.

Water promotes the sealing of nanoscale packing defects in folding proteins

A net dipole moment is shown to arise from a non-Debye component of water polarization created by nanoscale packing defects on the protein surface. Accordingly, the protein electrostatic field exerts a torque on the induced dipole, locally impeding the nucleation of ice at the protein-water interface. We evaluate the solvent orientation steering (SOS) as the reversible work needed to align the induced dipoles with the Debye electrostatic field and computed the SOS for the variable interface of a folding protein. The minimization of the SOS is shown to drive protein folding as evidenced by the entrainment of the total free energy by the SOS energy along trajectories that approach a Debye limit state where no torque arises. This result suggests that the minimization of anomalous water polarization at the interface promotes the sealing of packing defects, thereby maintaining structural integrity and committing the protein chain to fold.

Water under inner pressure: a dielectric spectroscopy study

Water is the most studied substance on Earth. However, it is not completely understood why its structural and dynamical properties give rise to some anomalous behaviors. Some of them emerge when experiments at low temperatures and/or high pressures are performed. Here we report dielectric measurements on cold water under macroscopically constrained conditions, i.e., water in a large container at constant volume that cannot freeze below the melting point. The inner pressure in these conditions shifts the α relaxation peak to similar frequencies as seen in ice Ih. At 267 K we observe a peculiar response possibly due to the Grotthuss mechanism. At 251 K (the triple point) ice III forms.

The interfacial-organized monolayer water film (MWF) induced "two-step" aggregation of nanographene: both in stacking and sliding assembly pathways

A computational investigation was carried out to understand the aggregation of nanoscale graphene with two typical pathways of stacking assembly and sliding assembly in water. The interfacial-organized monolayer water film (MWF) induced "two-step" aggregation of nanographene in both stacking and sliding assembly pathways was reported for the first time. By means of potential mean forces (PMFs) calculation, no energy barrier was observed during the sliding assembly of two graphene nanosheets, while the PMF profiles could be impacted by the contact forms of nanographene and the MWF within the interplate of two graphene nanosheets. To explore the potential physical basis of the "hindering role" of self-organized interfacial water, the dynamical and structural properties as well as the status of hydrogen bonds (H-bonds) for interfacial water were investigated. We found that the compact, ordered structure and abundant H-bonds of the MWF could be taken as the fundamental aspects of the "hindering role" of interfacial water for the hydrophobic assembly of nanographene. These findings are displaying a potential to further understand the hydrophobic assembly which mostly dominate the behaviors of nanomaterials, proteins etc. in aqueous solutions.