Relation between the two-body entropy and the relaxation time in supercooled water

Understanding pectin cross-linking in plant cell walls

Pectin is a major component of plant cells walls. The extent to which pectin chains crosslink with one another determines crucial properties including cell wall strength, porosity, and the ability of small, biologically significant molecules to access the cell. Despite its importance, significant gaps remain in our comprehension, at the molecular level, of how pectin cross-links influence the mechanical and physical properties of cell walls. This study employs a multidisciplinary approach, combining molecular dynamics simulations, experimental investigations, and mathematical modelling, to elucidate the mechanism of pectin cross-linking and its effect on cell wall porosity. The computational aspects of this work challenge the prevailing egg-box model, favoring instead a zipper model for pectin cross-linking, whilst our experimental work highlights the significant impact of pectin cross-linking on cell wall porosity. This work advances our fundamental understanding of the biochemistry underpinning the structure and function of the plant cell wall. This knowledge has important implications for agricultural biotechnology, informing us about the chemical properties of plant pectins that are best suited for improving crop resilience and amenability to biofuel extraction by modifying the cell wall.

Effect of confinement on water properties in super-hydrophilic pores using MD simulations with the mW model

CONTEXT

We explore the influence of strongly hydrophilic confinement on various properties of water, such as density, enthalpy, potential energy, radial distribution function, entropy, specific heat capacity, structural dynamics, and transition temperatures (freezing and melting temperatures), using monatomic water (mW) model. The properties of water are found to be dependent on confinement and the wall-fluid surface interaction. Hysteresis loops are observed for density, enthalpy, potential energy, and entropy around the transition temperatures, while the size of hysteresis loops varies with confinement and surface interaction. In smaller pore sizes (H ≤ 20), the solid phase displays a higher density compared to the liquid phase, which is unconventional behavior compared to bulk water systems due to the pronounced hydrophilic properties of the confinement surface. Specific heat capacity exhibits more oscillations in the confined system compared to bulk water, stemming from uneven enthalpy differences across equal temperature intervals. During phase transformation in both heating and quenching processes, there is an abrupt change observed in specific heat capacity. Confinement exerts a notable impact on entropy in the solid phase, but its influence is negligible in the liquid phase. At lower pore sizes (H < 25 Å), there is more fluctuation in freezing temperature for all wall-fluid interactions, which diminishes beyond pore sizes of H > 25 Å. Similarly, more oscillatory behavior is observed in melting temperatures at lower pore sizes (H < 40 Å), which diminishes at higher pore sizes (H > 40 Å). During the quenching process, a sudden jump in the in-plane orientational and tetrahedral order parameters indicates the formation of an ordered phase, specifically a diamond crystalline structure. The percentages of different crystalline structures (cubic diamond, hexagonal diamond, and 2D hexagonal) vary with both the confinement size and the wall-fluid interaction strength.

METHODS

Cooling and heating simulations are conducted with the mW water model using LAMMPS for different nanoscale confinement separation sizes ranging from 8.5 to 70 Å within the temperature range of 100-350 K. The water is modeled using two-body and three-body interaction potential (Stillinger-Weber potential) and the confinement is introduced using LJ 9-3 water-wall interaction potential. Entropy is calculated using RDF data obtained from the simulation experiments for each temperature point with increments or decrements of 2.5 K. The transition temperatures are estimated using the specific heat capacity analysis.

Understanding dynamics in coarse-grained models. IV. Connection of fine-grained and coarse-grained dynamics with the Stokes-Einstein and Stokes-Einstein-Debye relations

Applying an excess entropy scaling formalism to the coarse-grained (CG) dynamics of liquids, we discovered that missing rotational motions during the CG process are responsible for artificially accelerated CG dynamics. In the context of the dynamic representability between the fine-grained (FG) and CG dynamics, this work introduces the well-known Stokes-Einstein and Stokes-Einstein-Debye relations to unravel the rotational dynamics underlying FG trajectories, thereby allowing for an indirect evaluation of the effective rotations based only on the translational information at the reduced CG resolution. Since the representability issue in CG modeling limits a direct evaluation of the shear stress appearing in the Stokes-Einstein and Stokes-Einstein-Debye relations, we introduce a translational relaxation time as a proxy to employ these relations, and we demonstrate that these relations hold for the ambient conditions studied in our series of work. Additional theoretical links to our previous work are also established. First, we demonstrate that the effective hard sphere radius determined by the classical perturbation theory can approximate the complex hydrodynamic radius value reasonably well. Furthermore, we present a simple derivation of an excess entropy scaling relationship for viscosity by estimating the elliptical integral of molecules. In turn, since the translational and rotational motions at the FG level are correlated to each other, we conclude that the "entropy-free" CG diffusion only depends on the shape of the reference molecule. Our results and analyses impart an alternative way of recovering the FG diffusion from the CG description by coupling the translational and rotational motions at the hydrodynamic level.

Surface premelting and melting of colloidal glasses

The nature of liquid-to-glass transition is a major puzzle in science. A similar challenge exists in glass-to-liquid transition, i.e., glass melting, especially for the poorly investigated surface effects. Here, we assemble colloidal glasses by vapor deposition and melt them by tuning particle attractions. The structural and dynamic parameters saturate at different depths, which define a surface liquid layer and an intermediate glassy layer. The power-law growth of both layers and melting front behaviors at different heating rates are similar to crystal premelting and melting, suggesting that premelting and melting can be generalized to amorphous solids. The measured single-particle kinetics reveal various features and confirm theoretical predictions for glass surface layer.

Effects of structural inhomogeneity on equilibration processes in Langevin dynamics

In recent decades, computer experiments have led to an accurate and fundamental understanding of atomic and molecular mechanisms in fluids, such as different kinds of relaxation processes toward steady physical states. In this paper, we investigate how exactly the configuration of initial states in a molecular-dynamics simulation can affect the rates of decay toward equilibrium for the widely known Langevin canonical ensemble. For this purpose, we derive an original expression relating the system relaxation time τ_{sys} and the radial distribution function g(r) in the near-zero and high-density limit. We found that, for an initial state which is slightly marginally inhomogeneous in the number density of atoms, the system relaxation time τ_{sys} is much longer than that for the homogeneous case and an increasing function of the Langevin coupling constant, γ. We also found, during structural equilibration, g(r) at large distances approaches 1 from above for the inhomogeneous case and from below for the macroscopically homogeneous one.

Structure and dynamics of nanoconfined water and aqueous solutions

This review is devoted to discussing recent progress on the structure, thermodynamic, reactivity, and dynamics of water and aqueous systems confined within different types of nanopores, synthetic and biological. Currently, this is a branch of water science that has attracted enormous attention of researchers from different fields interested to extend the understanding of the anomalous properties of bulk water to the nanoscopic domain. From a fundamental perspective, the interactions of water and solutes with a confining surface dramatically modify the liquid's structure and, consequently, both its thermodynamical and dynamical behaviors, breaking the validity of the classical thermodynamic and phenomenological description of the transport properties of aqueous systems. Additionally, man-made nanopores and porous materials have emerged as promising solutions to challenging problems such as water purification, biosensing, nanofluidic logic and gating, and energy storage and conversion, while aquaporin, ion channels, and nuclear pore complex nanopores regulate many biological functions such as the conduction of water, the generation of action potentials, and the storage of genetic material. In this work, the more recent experimental and molecular simulations advances in this exciting and rapidly evolving field will be reported and critically discussed.

Role of local order in anomalous ion diffusion: Interrogation through tetrahedral entropy of aqueous solvation shells

Small rigid ions perturb the water structure around them significantly. At constant viscosity, alkali cations (Li⁺, Na⁺, and so on) exhibit an anomalous non-monotonic dependence of diffusivity on ion-size, in stark violation of the Stokes-Einstein expression. Although this is a well-known problem, we find that an entropic view of the problem can be developed, which provides valuable insight. The local entropy experienced by the solute ion is relevant here, which leads to the connection with local viscosity, discussed earlier by many. Due to the strong interactions with ions, the translational and rotational entropy of solvation water decreases sharply; however, an opposite effect comes from the disruption of the tetrahedral network structure of water near the charges. We compute the tetrahedral order of water molecules (q_tet) around the ion and suitably defined tetrahedral entropy [S(q_tet)] that is a contribution to the excess entropy of the system. Our results reveal that although the structural properties of the second shell become nearly identical to the bulk, S(q_tet) of the second shell is found to play an important role in giving rise to the non-monotonic ion-size dependence. The detailed study of the static and dynamic fluctuations in q_tet and the number of hydration water molecules provides interesting insights into correlation between the structure and dynamics; the smallest static fluctuation of q_tet for the first hydration shell water molecules of Li⁺ is indicative of the iceberg picture. The study of fluctuation properties of q_tet and the coordination number also reveals the role of the second hydration layer and could explain the anomalous behavior of the Rb⁺ ion.

Characterization of hydration water in supercooled water-trehalose solutions: The role of the hydrogen bonds network

The structural and dynamical properties of hydration water in aqueous solutions of trehalose are studied with molecular dynamics simulation. We simulate the systems in the supercooled region to investigate how the interaction with the trehalose molecules modifies the hydrogen bond network, the structural relaxation, and the diffusion properties of hydration water. The analysis is performed by considering the radial distribution functions, the residence time of water molecules in the hydration shell, the two body excess entropy, and the hydrogen bond water-water and water-trehalose correlations of the hydration water. The study of the two body excess entropy shows the presence of a fragile to strong crossover in supercooled hydration water also found in the relaxation time of the water-water hydrogen bond correlation function, and this is in agreement with predictions of the mode coupling theory and of previous studies of the oxygen-oxygen density correlators [A. Iorio et al., J. Mol. Liq. 282, 617 (2019); Sci. China: Phys., Mech. Astron. 62, 107011 (2019)]. The water-trehalose hydrogen bond correlation function instead evidences a strong to strong crossover in the relaxation time, and this crossover is related to a trehalose dynamical transition. This signals the role that the strong interplay between the soluted molecules and the surrounding solvent has in determining the dynamical transition common to both components of the system that happens upon cooling and that is similar to the well known protein dynamical transition. We connect our results with the cryoprotecting role of trehalose molecules.

Perspective: Excess-entropy scaling

This article gives an overview of excess-entropy scaling, the 1977 discovery by Rosenfeld that entropy determines properties of liquids like viscosity, diffusion constant, and heat conductivity. We give examples from computer simulations confirming this intriguing connection between dynamics and thermodynamics, counterexamples, and experimental validations. Recent uses in application-related contexts are reviewed, and theories proposed for the origin of excess-entropy scaling are briefly summarized. It is shown that if two thermodynamic state points of a liquid have the same microscopic dynamics, they must have the same excess entropy. In this case, the potential-energy function exhibits a symmetry termed hidden scale invariance, stating that the ordering of the potential energies of configurations is maintained if these are scaled uniformly to a different density. This property leads to the isomorph theory, which provides a general framework for excess-entropy scaling and illuminates, in particular, why this does not apply rigorously and universally. It remains an open question whether all aspects of excess-entropy scaling and related regularities reflect hidden scale invariance in one form or other.

Free energy and entropy of a dipolar liquid by computer simulations

Thermodynamic properties for a system composed of dipolar molecules are computed. Free energy is evaluated by means of the thermodynamic integration technique, and it is also estimated by using a perturbation theory approach, in which every molecule is modeled as a hard sphere within a square well, with an electric dipole at its center. The hard sphere diameter, the range and depth of the well, and the dipole moment have been calculated from properties easily obtained in molecular dynamics simulations. Connection between entropy and dynamical properties is explored in the liquid and supercooled states by using instantaneous normal mode calculations. A model is proposed in order to analyze translation and rotation contributions to entropy separately. Both contributions decrease upon cooling, and a logarithmic correlation between excess entropy associated with translation and the corresponding proportion of imaginary frequency modes is encountered. Rosenfeld scaling law between reduced diffusion and excess entropy is tested, and the origin of its failure at low temperatures is investigated.

Structural predictor for nonlinear sheared dynamics in simple glass-forming liquids

Glass-forming liquids subjected to sufficiently strong shear universally exhibit striking nonlinear behavior; for example, a power-law decrease of the viscosity with increasing shear rate. This phenomenon has attracted considerable attention over the years from both fundamental and applicational viewpoints. However, the out-of-equilibrium and nonlinear nature of sheared fluids have made theoretical understanding of this phenomenon very challenging and thus slower to progress. We find here that the structural relaxation time as a function of the two-body excess entropy, calculated for the extensional axis of the shear flow, collapses onto the corresponding equilibrium curve for a wide range of pair potentials ranging from harsh repulsive to soft and finite. This two-body excess entropy collapse provides a powerful approach to predicting the dynamics of nonequilibrium liquids from their equilibrium counterparts. Furthermore, the two-body excess entropy scaling suggests that sheared dynamics is controlled purely by the liquid structure captured in the form of the two-body excess entropy along the extensional direction, shedding light on the perplexing mechanism behind shear thinning.

Slow Dynamics and Structure of Supercooled Water in Confinement

We review our simulation results on properties of supercooled confined water. We consider two situations: water confined in a hydrophilic pore that mimics an MCM-41 environment and water at interface with a protein. The behavior upon cooling of the α relaxation of water in both environments is well interpreted in terms of the Mode Coupling Theory of glassy dynamics. Moreover, we find a crossover from a fragile to a strong regime. We relate this crossover to the crossing of the Widom line emanating from the liquid-liquid critical point, and in confinement we connect this crossover also to a crossover of the two body excess entropy of water upon cooling. Hydration water exhibits a second, distinctly slower relaxation caused by its dynamical coupling with the protein. The crossover upon cooling of this long relaxation is related to the protein dynamics.

Water: A Tale of Two Liquids

Water is the most abundant liquid on earth and also the substance with the largest number of anomalies in its properties. It is a prerequisite for life and as such a most important subject of current research in chemical physics and physical chemistry. In spite of its simplicity as a liquid, it has an enormously rich phase diagram where different types of ices, amorphous phases, and anomalies disclose a path that points to unique thermodynamics of its supercooled liquid state that still hides many unraveled secrets. In this review we describe the behavior of water in the regime from ambient conditions to the deeply supercooled region. The review describes simulations and experiments on this anomalous liquid. Several scenarios have been proposed to explain the anomalous properties that become strongly enhanced in the supercooled region. Among those, the second critical-point scenario has been investigated extensively, and at present most experimental evidence point to this scenario. Starting from very low temperatures, a coexistence line between a high-density amorphous phase and a low-density amorphous phase would continue in a coexistence line between a high-density and a low-density liquid phase terminating in a liquid–liquid critical point, LLCP. On approaching this LLCP from the one-phase region, a crossover in thermodynamics and dynamics can be found. This is discussed based on a picture of a temperature-dependent balance between a high-density liquid and a low-density liquid favored by, respectively, entropy and enthalpy, leading to a consistent picture of the thermodynamics of bulk water. Ice nucleation is also discussed, since this is what severely impedes experimental investigation of the vicinity of the proposed LLCP. Experimental investigation of stretched water, i.e., water at negative pressure, gives access to a different regime of the complex water diagram. Different ways to inhibit crystallization through confinement and aqueous solutions are discussed through results from experiments and simulations using the most sophisticated and advanced techniques. These findings represent tiles of a global picture that still needs to be completed. Some of the possible experimental lines of research that are essential to complete this picture are explored.