Decoupling of Translational Diffusion from the Viscosity of Supercooled Water: Role of Translational Jump Diffusion

Dual Nature of Large and Anisotropic Glass-Forming Molecules in Terms of Debye-Stokes-Einstein Relation Revealed

The fundamental Debye-Stokes-Einstein (DSE) relation between rotational relaxation times and shear viscosity attracts longstanding research interest as one of the most important characteristics of many glass-forming liquids. Here, we provide strong evidence, missing so far, for the relevance of anisotropy for DSE-related behavior. Dielectric spectroscopy and shear viscosity measurements were employed to get insight into the decoupling between reorientation relaxation times and viscosity for anisotropic glass-formers with dipole moments oriented parallel or perpendicular to the long molecular axis. We found that in the case of large and anisotropic molecules, the breakdown of DSE relation depends on the component of anisotropic rotation contributing to the dielectric response. Specifically, for glass-formers with dipole moment perpendicular to the long molecular axis, the DSE relation was found to be valid throughout the supercooled regime. Contrary, a departure from the DSE predictions in the intermediate supercooled regime, was observed for glass-formers where only short-axes contributions were relevant in dielectric response. MD simulations revealed differences in the mechanism of short and long axes reorientations suggesting that for anisotropic objects, not the reorientation mechanism itself, but the aspect of anisotropic motion, is the key to understanding the behavior of these systems in the context of DSE relation.

Stokes-Einstein Relation in Different Models of Water

The purpose of this paper is to discuss to which extent a microscopic version of the Stokes-Einstein (SE) relation without the hydrodynamic radius applies to liquid water. We demonstrate that the self-diffusion and shear viscosity data for five popular water models, recently reported by Ando [J. Chem. Phys. 159, 101102 (2023)], are in excellent agreement with the SE relation. The agreement with experimental results is also quite impressive. The limitations on the applicability of the SE relation are briefly discussed.

Breaking the Brownian barrier: models and manifestations of molecular diffusion in complex fluids

Over a century ago, Einstein formulated a precise mathematical model for describing Brownian motion. While this model adequately explains the diffusion of micron-sized particles in fluids, its limitations become apparent when applied to molecular self-diffusion in fluids. The foundational principles of Gaussianity and Markovianity, central to the Brownian diffusion paradigm, are insufficient for describing molecular diffusion, particularly in complex fluids characterized by intricate intermolecular interactions and hindered relaxation processes. This perspective delves into the nuanced behavior observed in diverse complex fluids, including molecular self-assembly systems, deep eutectic solvents, and ionic liquids, with a specific focus on modeling self-diffusion within these media. We explore the possibility of extending diffusion models to incorporate non-Gaussian and non-Markovian effects by augmenting the Brownian model using non-local diffusion equations. Furthermore, we validate the applicability of these models by utilizing them to describe results from quasielastic neutron scattering and MD simulations.

Computation of Overhauser dynamic nuclear polarization processes reveals fundamental correlation between water dynamics, structure, and solvent restructuring entropy

Hydration water dynamics, structure, and thermodynamics are crucially important to understand and predict water-mediated properties at molecular interfaces. Yet experimentally and directly quantifying water behavior locally near interfaces at the sub-nanometer scale is challenging, especially at interfaces submerged in biological solutions. Overhauser dynamic nuclear polarization (ODNP) experiments measure equilibrium hydration water dynamics within 8-15 angstroms of a nitroxide spin probe on instantaneous timescales (10 picoseconds to nanoseconds), making ODNP a powerful tool for probing local water dynamics in the vicinity of the spin probe. As with other spectroscopic techniques, concurrent computational analysis is necessary to gain access to detailed molecular level information about the dynamic, structural, and thermodynamic properties of water from experimental ODNP data. We chose a model system that can systematically tune the dynamics of water, a water-glycerol mixture with compositions ranging from 0 to 0.3 mole fraction glycerol. We demonstrate the ability of molecular dynamics (MD) simulations to compute ODNP spectroscopic quantities, and show that translational, rotational, and hydrogen bonding dynamics of hydration water align strongly with spectroscopic ODNP parameters. Moreover, MD simulations show tight correlations between the dynamic properties of water that ODNP captures and the structural and thermodynamic behavior of water. Hence, experimental ODNP readouts of varying water dynamics suggest changes in local structural and thermodynamic hydration water properties.

Unraveling the dynamic slowdown in supercooled water: The role of dynamic disorder in jump motions

When a liquid is rapidly cooled below its melting point without inducing crystallization, its dynamics slow down significantly without noticeable structural changes. Elucidating the origin of this slowdown has been a long-standing challenge. Here, we report a theoretical investigation into the mechanism of the dynamic slowdown in supercooled water, a ubiquitous yet extraordinary substance characterized by various anomalous properties arising from local density fluctuations. Using molecular dynamics simulations, we found that the jump dynamics, which are elementary structural change processes, deviate from Poisson statistics with decreasing temperature. This deviation is attributed to slow variables competing with the jump motions, i.e., dynamic disorder. The present analysis of the dynamic disorder showed that the primary slow variable is the displacement of the fourth nearest oxygen atom of a jumping molecule, which occurs in an environment created by the fluctuations of molecules outside the first hydration shell. As the temperature decreases, the jump dynamics become slow and intermittent. These intermittent dynamics are attributed to the prolonged trapping of jumping molecules within extended and stable low-density domains. As the temperature continues to decrease, the number of slow variables increases due to the increased cooperative motions. Consequently, the jump dynamics proceed in a higher-dimensional space consisting of multiple slow variables, becoming slower and more intermittent. It is then conceivable that with further decreasing temperature, the slowing and intermittency of the jump dynamics intensify, eventually culminating in a glass transition.

Enhanced fluidity of water in superhydrophobic nanotubes: estimating viscosity using jump-corrected confined Stokes-Einstein approach

Accurately predicting the viscosity of water confined within nanotubes is vital for various technological applications. Traditional methods have failed in this regard, necessitating a novel approach. We introduced the jump-corrected confined Stokes-Einstein (JCSE) method and now employ the same to estimate the viscosity and diffusion in superhydrophobic nanotubes. Our study covers a temperature range of 230-300 K and considers three nanotube diameters. Results show that water inside superhydrophobic nanotubes exhibits a significantly lower viscosity and higher diffusion than those inside hydrophobic nanotubes. Narrower nanotubes and lower temperatures accentuate these effects. Furthermore, water inside superhydrophobic nanotubes display a lower viscosity than bulk water, with the difference increasing at lower temperatures. This reduction is attributed to weaker water-water interactions caused by a lower water density in the interfacial region. These findings highlight the importance of interfacial water density and its influence on nanotube viscosity, shedding light on nanoscale fluid dynamics and opening avenues for diverse applications.

Anomalous lateral diffusion of lipids during the fluid/gel phase transition of a lipid membrane

A lipid membrane undergoes a phase transition from fluid to gel phase upon changing external thermodynamic conditions, such as decreasing temperature and increasing pressure. Extremophilic organisms face the challenge of preventing this deleterious phase transition. The main focus of their adaptive strategy is to facilitate effective temperature sensing through sensor proteins, relying on the drastic changes in packing density and membrane fluidity during the phase transition. Although the changes in packing density parameters due to the fluid/gel phase transition are studied in detail, the impact on membrane fluidity is less explored in the literature. Understanding the lateral diffusive dynamics of lipids in response to temperature, particularly during the fluid/gel phase transition, is albeit crucial. Here we have simulated the phase transition of a single component lipid membrane composed of dipalmitoylphosphatidylcholine (DPPC) lipids using a coarse-grained (CG) model and studied the changes of the structural and dynamical properties. It is observed that near the phase transition point, both fluid and gel phase domains coexist together. The dynamics remains highly non-Gaussian for a long time even when the mean square displacement reaches the Fickian regime at a much earlier time. This Fickian yet non-Gaussian diffusion (FnGD) is a characteristic of a highly heterogeneous system, previously observed for the lateral diffusion of lipids in raft mimetic membranes having liquid-ordered and liquid-disordered phases co-existing together. We have analyzed the molecular trajectories and calculated the jump-diffusion of the lipids, stemming from sudden jump translations, using a translational jump-diffusion (TJD) approach. An overwhelming contribution of the jump-diffusion of the lipids is observed suggesting anomalous diffusion of lipids during fluid/gel phase transition of the membrane. These results are important in unravelling the intricate nature of lipid diffusion during the phase transition of the membrane and open up a new possibility of investigating the most significant change of membrane properties during phase transition, which can be effectively sensed by proteins.

Noncanonical Relationship between Heterogeneity and the Stokes-Einstein Breakdown in Deep Eutectic Solvents

The relationship between Stokes-Einstein breakdown (SEB) and dynamical heterogeneity (DH) is of paramount importance in the physical chemistry of complex fluids. In this work, we employ neutron scattering to probe the DH and SEB in a series of deep eutectic solvents (DESs) composed of acetamide and lithium salts. Quasielastic neutron scattering experiments reveal SEB in the jump diffusion of acetamide, represented by a fractional Stokes-Einstein relationship. Among these DESs, lithium perchlorate exhibits the most pronounced SEB while lithium bromide displays the weakest. Concurrently, elastic incoherent neutron scans identify that bromide DES is the most heterogeneous and perchlorate is the least. For the first time, our study unveils a counterintuitive incommensurate relationship between DH and SEB. Further, it reveals the intricate contrasting nature of the SEB-DH relationship when investigated in proximity to the glass-transition temperature and further away from it.

Does Viscosity Decoupling Guarantee Dynamic Heterogeneity? A Clue through an Excitation and Emission Wavelength-Dependent Time-Resolved Fluorescence Anisotropy Study

Traditionally, deviation from Stokes-Einstein-Debye (SED) relation in terms of viscosity dependence of medium dynamics, i.e., τ x ∝ ( η T ) p with p ≠ 1, is taken as a signature of dynamic heterogeneity. However, it does not guarantee medium heterogeneity, as the decoupling may also originate from the deviation of the basic assumption of SED. Here, we developed a method to find a stronger relation between viscosity decoupling (p ≠ 1) and dynamic heterogeneity in terms of rotational motion. Our approach exploited the fact that in heterogeneous media, a solvatochromic probe will be solvated to a different extent at different microdomains (subpopulations), and photoselection of these subpopulations can be achieved by excitation or emission wavelength-dependent measurements. We hypothesized that the dynamics of a homogeneous system might show viscosity decoupling, but the extent of decoupling at different excitations (or at different emissions) should not be different. On the other hand, in a heterogeneous medium, this extent of viscosity decoupling (p-value) should be different at different excitations (or at different emissions). As proof of concept, we investigated three versatile solvent media: squalane (viscous molecular liquid), 1-ethyle-3-methylimidazolium ethyl sulfate ionic liquid (IL), and [0.78 acetamide + 0.22 LiNO3] deep eutectic solvent (DES). We found that squalane is homogeneous, although it shows fractional viscosity dependence (p ≠ 1). Interestingly, mild heterogeneity in IL and significant heterogeneity in the DES were observed. Overall, we conclude that the difference in the p-value as a function of excitation (or emission) wavelength-dependent might be a superior way for the detection of dynamic heterogeneity.

Absorption of pressurized methane in normal and supercooled p-xylene revealed via high-resolution neutron imaging

Supercooling of liquids leads to peculiarities which are scarcely studied under high-pressure conditions. Here, we report the surface tension, solubility, diffusivity, and partial molar volume for normal and supercooled liquid solutions of methane with p-xylene. Liquid bodies of perdeuterated p-xylene (p-C₈D₁₀), and, for comparison, o-xylene (o-C₈D₁₀), were exposed to pressurized methane (CH₄, up to 101 bar) at temperatures ranging 7.0-30.0 °C and observed at high spatial resolution (pixel size 20.3 μm) using a non-tactile neutron imaging method. Supercooling led to the increase of diffusivity and partial molar volume of methane. Solubility and surface tension were insensitive to supercooling, the latter substantially depended on methane pressure. Overall, neutron imaging enabled to reveal and quantify multiple phenomena occurring in supercooled liquid p-xylene solutions of methane under pressures relevant to the freeze-out in the production of liquefied natural gas.

Difference in "Supercooling" Affinity between (Acetamide + Na/KSCN) Deep Eutectics: Reflections in the Simulated Anomalous Motions of the Constituents and Solution Microheterogeneity Features

Deep depression of freezing points of ionic amide deep eutectic solvents (DESs) is known to exhibit a significant dependence on the identity of ions present in those systems and the nature of the functional group attached to the host amide. This deep depression of the freezing point is sometimes termed as "supercooling". For (acetamide + electrolyte) DESs, experiments have revealed signatures of ion-dependent spatiotemporal heterogeneity features. The focus of this work is to provide microscopic explanations of these experimentally observed macroscopic system properties in terms of particle jumps and insights about the origin of the cation dependence. For this purpose, extensive molecular dynamics simulations have been performed employing (acetamide + Na/KSCN) deep eutectics as representative ionic systems at 303, 318, 333, and 348 K. The individual translational motions of acetamide and the ions are followed, and their connections to solution heterogeneity are explored. The center-of-mass motion for Na⁺ has been found to be more anomalous than that for K⁺. This difference corroborates well with experimental reports on heterogeneous relaxations in these systems. Simulated viscosity coefficients and dynamic heterogeneity features also reflect this difference. Moreover, simulated reorientational relaxations of acetamide molecules in these ionic DESs suggest that a Na⁺-containing DES is more heterogeneous than the corresponding K⁺-containing system. Estimated void and neck distributions for acetamide molecules differ as the alkali metal ions differ. In brief, this study provides a detailed microscopic view of the cation dependence of the microheterogeneous relaxation dynamics of these DESs reported repeatedly by different experiments.

Self-Diffusion in Confined Water: A Comparison between the Dynamics of Supercooled Water in Hydrophobic Carbon Nanotubes and Hydrophilic Porous Silica

Confined liquids are model systems for the study of the metastable supercooled state, especially for bulk water, in which the onset of crystallization below 230 K hinders the application of experimental techniques. Nevertheless, in addition to suppressing crystallization, confinement at the nanoscale drastically alters the properties of water. Evidently, the behavior of confined water depends critically on the nature of the confining environment and the interactions of confined water molecules with the confining matrix. A comparative study of the dynamics of water under hydrophobic and hydrophilic confinement could therefore help to clarify the underlying interactions. As we demonstrate in this work using a few representative results from the relevant literature, the accurate assessment of the translational mobility of water molecules, especially in the supercooled state, can unmistakably distinguish between the hydrophilic and hydrophobic nature of the confining environments. Among the numerous experimental methods currently available, we selected nuclear magnetic resonance (NMR) in a field gradient, which directly measures the macroscopic translational self-diffusion coefficient, and quasi-elastic neutron scattering (QENS), which can determine the microscopic translational dynamics of the water molecules. Dielectric relaxation, which probes the re-orientational degrees of freedom, are also discussed.

Non-monotonic composition dependence of the breakdown of Stokes-Einstein relation for water in aqueous solutions of ethanol and 1-propanol: explanation using translational jump-diffusion approach

A series of experimental and simulation studies examined the validity of the Stokes-Einstein relationship (SER) of water in binary water/alcohol mixtures of different mixture compositions. These studies revealed a strong non-monotonic composition dependence of the SER with maxima at the specific alcohol mole fraction where the non-idealities of the thermodynamic and transport properties are observed. The translational jump-diffusion (TJD) approach elucidated the breakdown of the SER in pure supercooled water as caused by the jump translation of molecules. The breakdown of SER in the supercooled water/methanol binary mixture was successfully explained using the same TJD approach. To further generalize the picture, here we focus on the non-monotonic composition dependence of SER breakdown of water in two water/alcohol mixtures (water/ethanol and water/propanol) for a broad temperature range. In agreement with previous studies, maximum breakdown of SER is observed for the mixture with alcohol mole fraction x = 0.2. Diffusion of the water molecules at the maximum SER breakdown point is largely contributed by jump-diffusion. The residual-diffusion, obtained by subtracting the jump-diffusion from the total diffusion, approximately follows the SER for different compositions and temperatures. We also performed hydrogen (H-)bond dynamics and observed that the contribution of jump-diffusion is proportional to the total free energy of activation of breaking all H-bonds that exist around a molecule. This study, therefore, suggests that the more a molecule is trapped by H-bonding, the more likely it is to diffuse through the jump-diffusion mechanism, eventually leading to an increasing degree of SER breakdown.

Shear viscosity and Stokes-Einstein violation in supercooled light and heavy water

We report shear viscosity of heavy water supercooled 33K below its melting point, revealing a 15-fold increase compared to room temperature. We also confirm our previous data for the viscosity of supercooled light water and reach a better accuracy. Our measurements, based on the spontaneous Brownian motion of 350nm spheres, disagree at the lowest temperature with the only other available data, based on Poiseuille flow in a narrow capillary, which may have been biased by electro-osmotic effects. Here we provide a detailed description of the experiment and its analysis. We review the literature data about dynamic properties of water (viscosity, self-diffusion coefficient, and rotational correlation time), discuss their temperature dependence, and compare their decoupling in the two isotopes.

Single molecule demonstration of Debye-Stokes-Einstein breakdown in polystyrene near the glass transition temperature

Rotational-translational decoupling, in which translational motion is apparently enhanced over rotational motion in violation of Stokes-Einstein (SE) and Debye-Stokes-Einstein (DSE) predictions, has been observed in materials near their glass transition temperatures (T_g). This has been posited to result from ensemble averaging in the context of dynamic heterogeneity. In this work, ensemble and single molecule experiments are performed in parallel on a fluorescent probe in high molecular weight polystyrene near its T_g. Ensemble results show decoupling onset at approximately 1.15T_g, increasing to over three orders of magnitude at T_g. Single molecule measurements also show a high degree of decoupling, with typical molecules at T_g showing translational diffusion coefficients nearly 400 times higher than expected from SE/DSE predictions. At the single molecule level, higher degree of breakdown is associated with particularly mobile molecules and anisotropic trajectories, providing support for anomalous diffusion as a critical driver of rotational-translational decoupling and SE/DSE breakdown.

Experiments with high-molecular-weight polystyrene provide insights into the mechanisms behind rotational-translational decoupling in glassy systems. Specifically, particularly mobile molecules exhibiting anisotropic trajectories are found to play a key role in Debye-Stokes-Einstein breakdown.

Water Diffusion Proceeds via a Hydrogen-Bond Jump Exchange Mechanism

The self-diffusion of water molecules plays a key part in a broad range of essential processes in biochemistry, medical imaging, material science, and engineering. However, its molecular mechanism and the role played by the water hydrogen-bond network rearrangements are not known. Here we combine molecular dynamics simulations and analytic modeling to determine the molecular mechanism of water diffusion. We establish a quantitative connection between the water diffusion coefficient and hydrogen-bond jump exchanges, and identify the features that determine the underlying energetic barrier. We thus provide a unified framework to understand the coupling between translational, rotational, and hydrogen-bond dynamics in liquid water. It explains why these different dynamics do not necessarily exhibit identical temperature dependences although they all result from the same hydrogen-bond exchange events. The consequences for the understanding of water diffusion in supercooled conditions and for water transport in complex aqueous systems, including ionic, biological, and confined solutions, are discussed.

Non-Markovian rate theory on a multidimensional reaction surface: Complex interplay between enhanced configuration space and memory

A theory of barrier crossing rate on a multidimensional reaction energy surface is presented. The theory is a generalization of the earlier theoretical schemes to higher dimensions, with the inclusion of non-Markovian friction along both the reactive and the nonreactive coordinates. The theory additionally includes the bilinear coupling between the reactive and the nonreactive modes at the Hamiltonian level. Under suitable conditions, we recover the rate expressions of Langer and Hynes and establish a connection with the rate treatment of Pollak. Within the phenomenology of generalized Langevin equation description, our formulation provides an improvement over the existing ones because we explicitly include both the non-Markovian effects along the reaction coordinate and the bilinear coupling at the Hamiltonian level. At intermediate-to-large friction, an increase in dimensionality by itself tends to reduce the rate, while the inclusion of the memory effects increases the rate. The theory predicts an increase in rate when off-diagonal friction terms are included. We present a model calculation to study isomerization of a stilbene-like molecule using the prescription of Hochstrasser and co-workers on a two-dimensional reaction energy surface, employing Zwanzig–Bixon hydrodynamic theory of frequency-dependent friction. The calculated rate shows a departure from the predictions of Langer’s theory and also from the two-dimensional transition state theory.

Translational Jump-Diffusion of Hydroxide Ion in Anion Exchange Membrane: Deciphering the Nature of Vehicular Diffusion

Earlier, ab initio and reactive force-field-based molecular dynamics (MD) simulation studies suggested an overwhelming contribution of the vehicular diffusion in the total diffusion of hydroxide ions rather than structural diffusion. But does the vehicular diffusion occur via small-step displacement? This question is important to have an understanding of the real characteristics of vehicular diffusion. To answer this question, we perform a classical molecular dynamics simulation of a system containing a hydroxide ion exchange membrane polymer and hydroxide ion at different hydration levels and temperatures using the same molecular force field (Dubey, V. Chem. Phys. Lett. 2020, 755, 137802), which successfully captured the microscopic structure and dynamics of the system. We use the translational jump-diffusion approach, used previously in supercooled water for understanding the origin of breakdown of the Stokes-Einstein relation, to calculate the jump-diffusion coefficient of hydroxide ion and water in the anion exchange membrane. We have seen a significant role of hydration level and temperature in the mechanism of vehicular diffusion. In overhydrated membrane, both hydroxide ions and water molecules diffuse via both small- and large-step displacement. With decreasing hydration level and temperature, the diffusion is increasingly governed by the jump-diffusion mechanism. The larger contribution of jump-diffusion comes from the stronger caging of the diffusing species by the solvent at lower hydration levels and temperature. These results, therefore, suggest that the hydration level and temperature of the hydroxide ion exchange membrane determine the detailed mechanism of the vehicular diffusion of hydroxide ion, especially whether the diffusion follows hydrodynamics or not.

Exploring the validity of the Stokes-Einstein relation in supercooled water using nanomolecular probes

The breakdown of Stokes-Einstein relation in liquid water is one of the many anomalies that take place upon cooling and indicates the decoupling of diffusion and viscosity. It is hypothesized that these anomalies manifest due to the appearance of nanometer-scale spatial fluctuations, which become increasingly pronounced in the supercooled regime. Here, we explore the validity of the Stokes-Einstein relation in supercooled water using nanomolecular probes. We capture the diffusive dynamics of the probes using dynamic light scattering and target dynamics at different length scales by varying the probe size, from ≈100 nm silica spheres to molecular-sized polyhydroxylated fullerenes (≈1 nm). We find that all the studied probes, independent of size, display similar diffusive dynamics with an Arrhenius activation energy of ≈23 kJ mol^-1. Analysis of the diffusion coefficient further indicates that the probes, independent of their size, experience similar dynamic environment, which coincides with the macroscopic viscosity, while single water molecules effectively experience a comparatively lower viscosity. Finally, we conclude that our results indicate that the Stokes-Einstein relation is preserved for diffusion of probes in supercooled water T ≥ 260 K with size as small as ≈1 nm.

Breakdown of the Stokes-Einstein relation in supercooled water: the jump-diffusion perspective

Although water is the most ubiquitous liquid it shows many thermodynamic and dynamic anomalies. Some of the anomalies further intensify in the supercooled regime. While many experimental and theoretical studies have focused on the thermodynamic anomalies of supercooled water, fewer studies explored the dynamical anomalies very extensively. This is due to the intricacy of the experimental measurement of the dynamical properties of supercooled water. Violation of the Stokes-Einstein relation (SER), an important relation connecting the diffusion of particles with the viscosity of the medium, is one of the major dynamical anomalies. In absence of experimentally measured viscosity, researchers used to check the validity of SER indirectly using average translational relaxation time or α-relaxation time. Very recently, the viscosity of supercooled water was accurately measured at a wide range of temperatures and pressures. This allowed direct verification of the SER at different temperature-pressure thermodynamic state points. An increasing breakdown of the SER was observed with decreasing temperature. Increasing pressure reduces the extent of breakdown. Although some well-known theories explained the above breakdown, a detailed molecular mechanism was still elusive. Recently, a translational jump-diffusion (TJD) approach has been able to quantitatively explain the breakdown of the SER in pure supercooled water and an aqueous solution of methanol. The objective of this article is to present a detailed and state-of-the-art analysis of the past and present works on the breakdown of SER in supercooled water with a specific focus on the new TJD approach for explaining the breakdown of the SER.

Dehydration induced dynamical heterogeneity and ordering mechanism of lipid bilayers

Understanding the influence of dehydration on the membrane structure is crucial to control membrane functionality related to domain formation and cell fusion under anhydrobiosis conditions. To this end, we perform all-atom molecular dynamic simulations of 1,2-dimyristoyl-sn-glycero-3-phosphocholine dimyristoylphosphatidylcholine lipid membranes at different hydration levels at 308 K. As dehydration increases, the lipid area per head group decreases with an increase in bilayer thickness and lipid order parameters indicating bilayer ordering. Concurrently, translational and rotational dynamics of interfacial water (IW) molecules near membranes slow down. On the onset of bilayer ordering, the IW molecules exhibit prominent features of dynamical heterogeneity evident from non-Gaussian parameters and one-dimensional van Hove correlation functions. At a fully hydrated state, diffusion constants (D) of the IW follow a scaling relation, D∼τ_α ^-1, where the α relaxation time (τ_α) is obtained from self-intermediate scattering functions. However, upon dehydration, the relation breaks and the D of the IW follows a power law behavior as D∼τ_α ^-0.57, showing the signature of glass dynamics. τ_α and hydrogen bond lifetime calculated from intermittent hydrogen bond auto-correlation functions undergo a similar crossover in association with bilayer ordering on dehydration. The bilayer ordering is accompanied with an increase in fraction of caged lipids spanned over the bilayer surface and a decrease in fraction of mobile lipids due to the non-diffusive dynamics. Our analyses reveal that the microscopic mechanism of lipid ordering by dehydration is governed by dynamical heterogeneity. The fundamental understanding from this study can be applied to complex bio-membranes to trap functionally relevant gel-like domains at room temperature.

Supercooled liquid-like dynamics in water near a fully hydrated titania surface: Decoupling of rotational and translational diffusion

We report an ab initio molecular dynamics (MD) simulation investigating the effect of a fully hydrated surface of TiO₂ on the water dynamics. It is found that the universal relation between the rotational and translational diffusion characteristics of bulk water is broken in the water layers near the surface with the rotational diffusion demonstrating progressive retardation relative to the translational diffusion when approaching the surface. This kind of rotation-translation decoupling has so far only been observed in the supercooled liquids approaching glass transition, and its observation in water at a normal liquid temperature is of conceptual interest. This finding is also of interest for the application-significant studies of the water interaction with fully hydrated nanoparticles. We note that this is the first observation of rotation-translation decoupling in an ab initio MD simulation of water.

Breakdown of the Stokes-Einstein Relation in Supercooled Water/Methanol Binary Mixtures: Explanation Using the Translational Jump-Diffusion Approach

A recent experiment has directly checked the validity of the Stokes-Einstein (SE) relation for pure water, pure methanol, and their binary mixtures of three different compositions at different temperatures. The effect of composition on the nature of breakdown of the SE relation is interesting. While in the majority of the systems, an increasing SE breakdown is observed with decreasing temperature, the breakdown is already significant at higher temperatures for the equimolar mixture. Violations of the SE relation in pure supercooled water at different temperatures and pressures have been previously explained using the translational jump-diffusion (TJD) approach, which provides a fundamental molecular basis, by directly connecting the SE breakdown with jump-diffusion of the molecules. We have used the same TJD approach for explaining the SE breakdown for the methanol/water binary mixtures of compositions studied in the experiment over a wide range of temperatures between 220 K and 300 K. We have understood that the jump-diffusion is the key responsible factor for the SE breakdown. The maximum jump-diffusion contribution gives rise to the early SE breakdown observed for the equimolar mixture observed in the experiment. This study, therefore, provides molecular insight into the SE breakdown for the supercooled water/methanol binary mixture, as found in the experiment.

Breakdown of the Stokes-Einstein Equation for Solutions of Water in Oil Reverse Micelles

An experimental study is presented for the reverse micellar system of 15% by mass polydisperse hexaethylene glycol monodecylether (C₁₀E₆) in cyclohexane with varying amounts of added water up to 4% by mass. Measurements of viscosity and self-diffusion coefficients were taken as a function of temperature between 10 and 45 °C at varying sample water loads but fixed C₁₀E₆/cyclohexane composition. The results were used to inspect the validity of the Stokes-Einstein equation for this system. Unreasonably small reverse average micelle radii and aggregation numbers were obtained with the Stokes-Einstein equation, but reasonable values for these quantities were obtained using the ratio of surfactant-to-cyclohexane self-diffusion coefficients. While bulk viscosity increased with increasing water load, a concurrent expected decrease of self-diffusion coefficient was only observed for the surfactant and water but not for cyclohexane, which showed independence of water load. Moreover, a spread of self-diffusion coefficients was observed for the protons associated with the ethylene oxide repeat unit in samples with polydisperse C₁₀E₆ but not in a sample with monodisperse C₁₀E₆. These findings were interpreted by the presence of reverse micelle to reverse micelle hopping motions that with higher water load become increasingly selective toward C₁₀E₆ molecules with short ethylene oxide repeat units, while those with long ethylene oxide repeat units remain trapped within the reverse micelle because of the increased hydrogen bonding interactions with the water inside the growing core of the reverse micelle. Despite the observed breakdown of the Stokes-Einstein equation, the temperature dependence of the viscosities and self-diffusion coefficients was found to follow Arrhenius behavior over the investigated range of temperatures.

Cryostabilization of the Cell Membrane of a Psychrotolerant Bacteria via Homeoviscous Adaptation

Homeoviscous adaptation (maintenance of a critical balance between the saturated and unsaturated lipids) of the cell membrane of psychrotolerant bacteria is essential to protect them against freeze-thaw cycle. But how does the homeoviscous adaptation protect the cell membrane during cold stress? In this Letter, we answer this question using a coarse-grained molecular dynamics simulation technique. On the basis of the reported fatty acid profiles of psychrotolerant bacteria at different temperatures, multiple lipid membranes are simulated at a wide range of temperatures between 250 and 300 K. We explicate how the homeoviscous adaptation minimizes the effect of cold stress on the structure and fluidity of the membrane. Partial freezing of the saturated lipid domain occurs with the self-aggregation of saturated and unsaturated lipids near the melting temperature of the unadapted lipid membrane. The gel-like phase provides necessary local packing density that can be sensed by sensor proteins responsible for the homeoviscous adaptation.

Dynamic coupling of a hydration layer to a fluid phospholipid membrane: intermittency and multiple time-scale relaxations

Understanding the coupling of a hydration layer and a lipid membrane is crucial to gaining access to membrane dynamics and understanding its functionality towards various biological processes.

Understanding the Origin of the Breakdown of the Stokes-Einstein Relation in Supercooled Water at Different Temperature-Pressure Conditions

A recent experiment has measured the viscosity of water down to approximately 244 K and up to 300 MPa. The correct viscosity and translational diffusivity data at various temperature-pressure (T-P) state points allowed for checking the validity of the Stokes-Einstein (SE) relation, which accounts for the coupling between translational self-diffusion and medium viscosity. The diffusion-viscosity decoupling increases with decreasing temperature, but the increasing pressure reduces the extent of the decoupling. Earlier simulation studies explained the breakdown of the SE relation in terms of the location of the Widom line, emanating from the liquid-liquid critical point (LLCP). Although these studies made a significant contribution to the current understanding of the above phenomena, a detailed molecular picture is still lacking. Recently, our group has explained the diffusion-viscosity decoupling from a jump-diffusion perspective. The jump-diffusion coefficient, emanating from the jump translation of water molecules, is calculated using a quantitative approach for different temperatures at ambient pressure. It has been observed that jump-diffusion is the key factor for diffusion-viscosity decoupling in supercooled water. The same method is adopted in the present work to estimate the jump-diffusion coefficient for different T-P state points and, thereby, explains the role of jump-diffusion for the different extents of the SE relation breakdown at different pressures. The residual diffusion coefficient, the other component of the total diffusion that originates from small step displacement and that is calculated by subtracting the jump-diffusion coefficient from the total diffusion, is seen to be fairly coupled to the viscosity at the entire range of temperature and pressure. Furthermore, we have calculated the average number of H-bonds per water molecule and the tetrahedral order for different T-P state points and investigated an approximate correlation between the average local structure and the contribution of the jump-diffusion to the total diffusion of water. This study, therefore, puts forward a new perspective for explaining the SE relation breakdown in supercooled water under different pressure conditions.