Structure and dynamics of supercooled water in neutral confinements

Transport Mechanism in Nafion Revealed by Detailed Comparison of 1H and 17O Nuclear Magnetic Resonance Diffusion Coefficients

We use 1H and 17O NMR static field gradient diffusometry to measure self-diffusion coefficients of protons (DH) and oxygens (DO) in Nafion 212 with various hydration levels (λ = 4-18). For all samples and both nuclei, we obtain activation energies (Ea) of ≈0.19 eV. Analyzing the hydration-level dependence of DH and DO, we find DO/DH ≈ 1 at λ ≈ 18, resembling the situation in bulk water, while oxygen diffusion becomes faster than proton diffusion when the water content is decreased, leading to DO/DH ≈ 1.2 at λ ≈ 4. A comparison with literature data for acidic bulk solutions implies that faster oxygen than proton diffusion results from the existence of the polymer framework. To rationalize the observed ratios DO/DH ≥ 1, we consider a bimodal dynamical model in which the interactions of H+(H2O)m ions with neighboring SO3- groups lead to slower water dynamics in the vicinity of the polymer framework than in the center of the water nanochannels.

Three-Dimensional Confinement of Water: H2O Exhibits Long-Range (>50 nm) Structure while D2O Does Not

Water is the liquid of life thanks to its three-dimensional adaptive hydrogen (H)-bond network. Confinement of this network may lead to dramatic structural changes influencing chemical and physical transformations. Although confinement effects occur on a <1 nm length scale, the upper length scale limit is unknown. Here, we investigate 3D-confinement over lengths scales ranging from 58-140 nm. By confining water in zwitterionic liposomes of different sizes and measuring the change in H-bond network conformation using second harmonic scattering (SHS), we determined long-range confinement effects in light and heavy water. D₂O displays no detectable 3D-confinement effects <58 nm (<3 × 10⁶ D₂O molecules). H₂O is distinctly different. The vesicle enclosed inner H-bond network has a different conformation compared to the outside network and the SHS response scales with the volume of the confining space. H₂O displays confinement effects over distances >100 nm (>2 × 10⁷ H₂O molecules).

Structural and Dynamical Properties of Liquids in Confinements: A Review of Molecular Dynamics Simulation Studies

Molecular dynamics (MD) simulations are a powerful tool for detailed studies of altered properties of liquids in confinement, in particular, of changed structures and dynamics. They allow, on one hand, for perfect control and systematic variation of the geometries and interactions inherent in confinement situations and, on the other hand, for type-selective and position-resolved analyses of a huge variety of structural and dynamical parameters. Here, we review MD simulation studies on various types of liquids and confinements. The main focus is confined aqueous systems, but also ionic liquids and polymer and silica melts are discussed. Results for confinements featuring different interactions, sizes, shapes, and rigidity will be presented. Special attention will be given to situations in which the confined liquid and the confining matrix consist of the same type of particles and, hence, disparate liquid-matrix interactions are absent. Findings for the magnitude and the range of wall effects on molecular positions and orientations and on molecular dynamics, including vibrational motion and structural relaxation, are reviewed. Moreover, their dependence on the parameters of the confinement and their relevance to theoretical approaches to the glass transition are addressed.

Rotational dynamics of proteins in nanochannels: role of solvent's local viscosity

Viscosity variation of solvent in local regions near a solid surface, be it a biological surface of a protein or an engineered surface of a nanoconfinement, is a direct consequence of intermolecular interactions between the solid body and the solvent. The current coarse-grained molecular dynamics study takes advantage of this phenomenon to investigate the anomaly in a solvated protein's rotational dynamics confined using a representative solid matrix. The concept of persistence time, the characteristic time of structural reordering in liquids, is used to compute the solvent's local viscosity. With an increase in the degree of confinement, the confining matrix significantly influences the solvent molecule's local viscosity present in the protein hydration layer through intermolecular interactions. This effect contributes to the enhanced drag force on protein motion, causing a reduction in the rotational diffusion coefficient. Simulation results suggest that the direct matrix-protein non-bonded interaction is responsible for the occasional jump and discontinuity in orientational motion when the protein is in very tight confinement.

MD simulations of charged binary mixtures reveal a generic relation between high- and low-temperature behavior

Experimental studies of the glassy slowdown in molecular liquids indicate that the high-temperature activation energy E_∞ of glass-forming liquids is directly related to their glass transition temperature T_g. To further investigate such a possible relation between high- and low-temperature dynamics in glass-forming liquids, we analyze the glassy dynamics of binary mixtures using molecular dynamics simulations. We consider a binary mixture of charged Lennard-Jones particles and vary the partial charges of the particles and, thus, the high-temperature activation energy and the glass transition temperature of the system. Based on previous results, we introduce a phenomenological model describing relaxation times over the whole temperature regime from high temperatures to temperatures well inside the supercooled regime. By investigating the dynamics of both particle species on molecular and diffusive length scales along isochoric and isobaric pathways, we find a quadratic charge dependence of both E_∞ and T_g, resulting in an approximately constant ratio of both quantities independent of the underlying observable, the thermodynamic ensemble, and the particle species, and this result is robust against the actual definition of T_g. This generic relation between the activation energy and the glass transition temperature indicates that high-temperature dynamics and the glassy slowdown are related phenomena, and the knowledge of E_∞ may allow us to approximately predict T_g.

Replica exchange MD simulations of two-dimensional water in graphene nanocapillaries: rhombic versus square structures, proton ordering, and phase transitions

The hydrogen bond patterns, proton ordering, and phase transitions of monolayer ice in two-dimensional hydrophobic confinement are fundamentally different from those found for bulk ice. To investigate the behavior of quasi-2D ice, we perform molecular dynamics simulations of water confined between fixed graphene plates at a distance of 0.65 nm. While experimental results are still limited and theoretical investigations are often based on a single, often empirically based force field model, this work presents a systematic study modeling the water-graphene interaction by effective Lennard-Jones potentials previously derived from high-level ab initio CCSD(T) calculations of water adsorbed on graphene [Phys. Chem. Chem. Phys., 2013, 15, 4995]. For the water-water interaction different water force fields, i.e. SPCE, TIP3P, TIP4P, TIP4P/ICE, and TIP5P, are used. The water occupancy of the graphene capillary at a pressure of 1000 MPa is determined to be between 13.5 and 13.9 water molecules per square nanometer, depending on the choice of the water force field. Based on these densities, we explore the structure and dynamics of quasi-2D water for temperatures ranging from 200 K to about 600 K for each of the five force fields. To ensure complete sampling of the configurational space and to overcome the barriers separating metastable structures, these simulations are based on the replica exchange molecular dynamics technique. We report different tetragonal hydrogen bond patterns, which are classified as nearly square or as rhombic. While many of these arrangements are essentially flat, in some cases puckered arrangements are found, too. Also the proton ordering of the quasi-2D water structures is considered, allowing us to identify them as ferroelectric, ferrielectric or antiferroelectric. For temperatures between 200 K and 400 K we find several second-order phase transitions from one ice structure to another, changing in many cases both the arrangements of the oxygen atoms and the proton ordering. For temperatures between 400 K and 600 K there are melting-like transitions from a monolayer of ice to a monolayer of liquid water. These first-order phase transitions have a latent heat between 3.4 and 4.0 kJ mol^-1. Both the values of the transition temperatures and of the latent heats display considerable model dependence for the five different water models investigated here.

Effects of confinement on supercooled tetrahedral liquids

We use molecular dynamics simulations to ascertain the effects of geometrical restriction on glass-forming tetrahedral liquids. Striving for a broad approach, we study families of waterlike and silicalike liquids, for which we systematically scale the partial charges and, hence, the relevance of the tetrahedral networks. The confined liquids and the confining matrices consist of the same type of particles to avoid disruptive interactions and distorted structures at the interfaces. Spatially resolved analyses show that these neutral confinements still impose static mobility gradients and density correlations on the liquids. We quantify the increasing degree and range of the altered properties upon cooling. For both families of models, common relations describe the confinement effects of all systems with tetrahedral order, while deviations occur for systems with lower polarities and different structures. The observations are rationalized by considering the fact that a pinned wall imprints a static energy landscape to a neighboring liquid. We explore the properties of this landscape based on changes in vibrational motion and structural relaxation and find that typical barrier heights amount to two to three times the activation energy of bulk dynamics. Combining the present and previous results, we predict the evolution of confinement effects down to the glass transition temperature for liquids without fragile-to-strong crossover. In addition, it is found for silicalike liquids that the temperature dependence of dynamic and static correlation lengths from confinement studies is not affected when cooling through fragile-to-strong transitions of the bulk materials, casting doubt on the relevance of these length scales for the glassy slowdown.

DMA study of water's glass transition in nanoscale confinement

Dynamic mechanical analysis (DMA) measurements of water confined in nanoporous silica have been performed as a function of temperature and frequency for different pore sizes (2.5-10 nm) at heating and cooling. Most of the data show three processes, P1, P2 and P3, where P1 and P2 depend on measurement frequency and P3 does not. The characteristic shift of P3 with pore size shows that this process corresponds to freezing/melting of "internal water", i.e. in the core of the pores. Thermal expansion data indicate - in agreement with e.g. [A. Taschin, P. Bartolini and R. Torre, Meas. Sci. Technol., 2017, 28, 014009] - that in all our nanoporous systems about 2 layers of water remain liquid much below the freezing point. Dynamic elastic measurements show clear signatures of glass freezing of this supercooled water in the vicinity of P1. Extrapolating the DMA data to the timescale (103 s) of adiabatic calorimetry unveils a systematic behaviour: P1(T) shows a clear size dependence for a broad range of pore diameters, i.e. 2.5 nm ≤ d ≤ 52 nm, implying (together with the corresponding activation energy 0.5 eV) that P1 corresponds to the glass-liquid transition of a few layers of supercooled water at Tg(d). An extrapolation of Tg(d) to d → ∞ yields Tg(∞) ≈ 136 K, the traditional value for bulk water. The small (liquid like) value of Young's modulus in a temperature region above P1 is most naturally explained assuming that the supercooled water in this range is still liquid, implying that Tg values of 160 K or even 210 K - as suggested by various authors - are unlikely.

Molecular Dynamics Simulations of Water, Silica, and Aqueous Mixtures in Bulk and Confinement

Aqueous systems are omnipresent in nature and technology. They show complex behaviors, which often originate in the existence of hydrogen-bond networks. Prominent examples are the anomalies of water and the non-ideal behaviors of aqueous solutions. The phenomenology becomes even richer when aqueous liquids are subject to confinement. To this day, many properties of water and its mixtures, in particular, under confinement, are not understood. In recent years, molecular dynamics simulations developed into a powerful tool to improve our knowledge in this field. Here, our simulation results for water and aqueous mixtures in the bulk and in various confinements are reviewed and some new simulation data are added to improve our knowledge about the role of interfaces. Moreover, findings for water are compared with results for silica, exploiting that both systems form tetrahedral networks.

Properties of Hydrogen-Bonded Liquids at Interfaces

Effects of interfaces on hydrogen-bonded liquids play major roles in nature and technology. Despite their importance, a fundamental understanding of these effects is still lacking. In large parts, this shortcoming is due to the high complexity of these systems, leading to an interference of various interactions and effects. Therefore, it is advisable to take gradual approaches, which start from well designed and defined model systems and systematically increase the level of intricacy towards more complex mimetics. Moreover, it is necessary to combine insights from a multitude of methods, in particular, to link novel preparation strategies and comprehensive experimental characterization with inventive computational and theoretical modeling. Such concerted approach was taken by a group of preparative, experimentally, and theoretically working scientists in the framework of Research Unit FOR 1583 funded by the Deutsche Forschungsgemeinschaft (German Research Foundation). This special issue summarizes the outcome of this collaborative research. In this introductory article, we give an overview of the covered topics and the main results of the whole consortium. The following contributions are review articles or original works of individual research projects.

2H NMR Studies on the Dynamics of Pure and Mixed Hydrogen-Bonded Liquids in Confinement

2H NMR is used to ascertain dynamical behaviors of pure and mixed hydrogen-bonded liquids in bulk and in confinement. Detailed comparisons of previous and new results in broad dynamic and temperature ranges reveal that confinement effects differ for various liquids and confinements. For water, molecular reorientation strongly depends on the confinement size, with much slower and less fragile structural relaxation under more severe geometrical restriction. Moreover, a dynamical crossover occurs when a fraction of solid water forms so that the dynamics of the fraction of liquid water becomes even more restricted and, as a consequence, changes from bulk-like to interface-dominated. For glycerol, by contrast, confinement has weak effects on the reorientation dynamics. Mixed hydrogen-bonded liquids show even more complex dynamical behaviors. For aqueous solutions, the temperature dependence of the structural relaxation becomes discontinuous when the concentration changes due to a freezing of water fractions. This tendency for partial crystallization is enhanced rather than reduced by confinement, because different liquid-matrix interactions of the molecular species induce micro-phase segregation, which facilitates ice formation in water-rich regions. In addition, dynamical couplings at solvent-protein interfaces are discussed. It is shown that, on the one hand, solvent dynamics are substantially slowed down at protein surfaces and, on the other hand, protein dynamics significantly depend on the composition and, thus, the viscosity of the solvent. Furthermore, a protein dynamical transition occurs when the amplitude of water-coupled restricted backbone dynamics vanishes upon cooling.

Ab-Initio Molecular Dynamics Simulations and Calculations of Spectroscopic Parameters in Hydrogen-Bonding Liquids in Confinement (Project 8)

We investigate the effect of several nanoscale confinements on structural and dynamical properties of liquid water and binary aqueous mixtures. By means of molecular dynamics simulations based on density functional theory and atomistic force fields. Our main focus is on the dependence on the structure and the hydrogen-bonding-network of the liquids near the confinement interface at atomistic resolution. As a complementary aspect, spatially resolved profiles of the proton NMR chemical shift values are used to quantify the local strength of the hydrogen-bond-network.

Phase separation of triethylamine and water in native and organically modified silica nanopores

A mixture of triethylamine and water is a lower critical solution temperature system that demixes (separates into individual phases) on heating. Differential scanning calorimetry has been applied to study the process of demixing in native and organically modified silica nanopores whose size varied from 4 to 30 nm. It has been found that in both types of nanopores, the temperature and enthalpy of demixing decrease significantly with decreasing the pore size. Isoconversional kinetic analysis has been utilized to determine the activation energy and pre-exponential factor of the process. It has been demonstrated that the depression of the transition temperature upon nanoconfinement is associated with acceleration of the process due to lowering of the activation energy. Nanoconfinement has also been found to lower the pre-exponential factor of the process that has been linked to a decrease in the molecular mobility.

Continuous melting through a hexatic phase in confined bilayer water

Liquid water is not only of obvious importance but also extremely intriguing, displaying many anomalies that still challenge our understanding of such an a priori simple system. The same is true when looking at nanoconfined water: The liquid between constituents in a cell is confined to such dimensions, and there is already evidence that such water can behave very differently from its bulk counterpart. A striking finding has been reported from computer simulations for two-dimensionally confined water: The liquid displays continuous or discontinuous melting depending on its density. In order to understand this behavior, we have analyzed the melting exhibited by a bilayer of nanoconfined water by means of molecular dynamics simulations. At high density we observe the continuous melting to be related to the phase change of the oxygens only, with the hydrogens remaining liquidlike throughout. Moreover, we find an intermediate hexatic phase for the oxygens between the liquid and a triangular solid ice phase, following the Kosterlitz-Thouless-Halperin-Nelson-Young theory for two-dimensional melting. The liquid itself tends to maintain the local structure of the triangular ice, with its two layers being strongly correlated yet with very slow exchange of matter. The decoupling in the behavior of the oxygens and hydrogens gives rise to a regime in which the complexity of water seems to disappear, resulting in what resembles a simple monoatomic liquid. This intrinsic tendency of our simulated water may be useful for understanding novel behaviors in other confined and interfacial water systems.

A molecular dynamics simulations study on the relations between dynamical heterogeneity, structural relaxation, and self-diffusion in viscous liquids

We perform molecular dynamics simulations for viscous liquids to study the relations between dynamical heterogeneity, structural (α) relaxation, and self-diffusion. For atomistic models of supercooled water, polymer melts, and an ionic liquid, we characterize the space-time characteristics of dynamical heterogeneity by the degree of deviations from Gaussian displacement statistics (α2), the size of clusters comprising highly mobile particles (S(w)), and the length of strings consisting of cooperatively moving particles (L(w)). Comparison of our findings with previous simulation results for a large variety of viscous liquids, ranging from monoatomic liquids to silica melt, reveals a nearly universal decoupling between the time scales of maximum non-Gaussian parameter (τ(α2)) and the time constant of the α relaxation (τ(α)) upon cooling, explicitly, τ(α2) ∝τ(α)(3/4). Such uniform relation was not observed between the peak times of S(w) or L(w) and τ(α). On the other hand, the temperature-dependent time scale of maximum string length (τ(L)) follows the inverse of the self-diffusion coefficient (D) for various systems at sufficiently low temperatures, i.e., τ(L) ∝ D(-1). These observations are discussed in view of a breakdown of the Stokes-Einstein relation for the studied systems. It is found that the degree of deviation from this relation is correlated with the stretching of the α relaxation.

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.

Slow Water Dynamics near a Glass Transition or a Solid Interface: A Common Rationale

Performing molecular dynamics simulations, we investigate the enormous slowdowns of water dynamics when approaching a glass transition or a solid interface. We show that both effects can be described on common grounds within a theoretical framework, which was recently proposed by Schweizer et al. and considers coupled local hopping and elastic distortion. For confined water, we correctly describe the variation of the α-relaxation time, τα, as a function of both temperature and position with respect to the interface. Exploiting our knowledge of a cooperative length scale ξ(T) from the confinement studies, we quantitatively rationalize the glassy slowdown, τα(T), and the Stokes-Einstein breakdown of bulk water. For both confined and bulk liquid, variations of the α-relaxation time are intimately related to changes of the cage-rattling amplitude.

Bulk and interfacial liquid water as a transient network

The special macroscopic properties of liquid water stem from its structure as a complex network of molecules connected by hydrogen bonds. While the dynamics of single molecules within this network has been extensively investigated, only little attention has been paid to the closed loops (meshes) of hydrogen-bonded molecules which determine the network topology. Using molecular dynamics simulations we analyze the size, shape, geometrical arrangement, and dynamical stability of loops containing up to 10 hydrogen bonds. We find that six-membered loops in liquid water even at room temperature retain a striking similarity with the well-known structure of ice. Analyzing the network dynamics we find that rings of more than five hydrogen bonds are stabilized compared to a random collection containing the same number of single bonds. We finally show that in the vicinity of hydrophobic and hydrophilic interfaces loops arrange in a preferred orientation.

Structure and dynamics of TIP3P, TIP4P, and TIP5P water near smooth and atomistic walls of different hydroaffinity

We perform molecular dynamics simulations to observe the structure and dynamics of water using different water models (TIP3P, TIP4P, TIP5P) at ambient conditions, constrained by planar walls, which are either modeled by smooth potentials or regular atomic lattices, imitating the honeycomb-structure of graphene. We implement walls of different hydroaffinity, different lattice constant, and different types of interaction with the water molecules. We find that in the hydrophobic regime the smooth wall generally represents a good abstraction of the atomically rough walls, while in the hydrophilic regime there are noticeable differences in structure and dynamics between all stages of wall roughness. For a small lattice constant however the smooth and the atomically rough wall still share a number of structural and dynamical similarities. Out of the three water models, TIP5P water shows the largest degree of tetrahedral ordering and is often the one that is least perturbed by the presence of the wall.

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.

NMR studies on the temperature-dependent dynamics of confined water

NMR studies of water in nanoscopic confinements of various sizes reveal two dynamical crossovers related to a partial solidification of internal molecules and a glass transition of interfacial molecules, respectively.