Diffusivity anomaly in modified Stillinger-Weber liquids

Unraveling thermodynamic anomalies of water: A molecular simulation approach to probe the two-state theory with atomistic and coarse-grained water models

Thermodynamic and dynamic anomalies of water play a crucial role in supporting life on our planet. The two-state theory attributes these anomalies to a dynamic equilibrium between locally favored tetrahedral structures (LFTSs) and disordered normal liquid structures. This theory provides a straightforward, phenomenological explanation for water's unique thermodynamic and dynamic characteristics. To validate this two-state feature, it is critical to unequivocally identify these structural motifs in a dynamically fluctuating disordered liquid. In this study, we employ a recently introduced structural parameter (θavg) that characterizes the local angular order within the first coordination shell to identify these LFTSs through molecular dynamics simulations. We employ both realistic water models with a liquid-liquid critical point (LLCP) and a coarse-grained water model without an LLCP to study water's anomalies in low-pressure regions below 2 kbar. The two-state theory consistently describes water's thermodynamic anomalies in these models, both with and without an LLCP. This suggests that the anomalies predominantly result from the two-state features rather than criticality, particularly within experimentally accessible temperature-pressure regions.

Thermodynamic anomalies in silicon and the relationship to the phase diagram

The evolution of thermodynamic anomalies are investigated in the pressure-temperature (pT) plane for silicon using the well-established Stillinger-Weber potential. Anomalies are observed in the density, compressibility and heat capacity. The relationships between them and with the liquid stability limit are investigated and related to the known thermodynamic constraints. The investigations are extended into the deeply supercooled regime using replica exchange techniques. Thermodynamic arguments are presented to justify the extension to low temperature, although a region of phase space is found to remain inaccessible due to unsuppressible crystallisation. The locus corresponding to the temperature of minimum compressibility is shown to display a characteristic 'S'-shape in thepTprojection which appears correlated with the underlying crystalline phase diagram. The progression of the anomalies is compared to the known underlying phase diagrams for both the crystal/liquid and amorphous/liquid states. The locations of the anomalies are also compared to those obtained from previous simulation work and (limited) experimental observations.

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.

Water-like anomalies as a function of tetrahedrality

Water is the most common and yet least understood material on Earth. Despite its simplicity, water tends to form tetrahedral order locally by directional hydrogen bonding. This structuring is known to be responsible for a vast array of unusual properties, e.g., the density maximum at 4 °C, which play a fundamental role in countless natural and technological processes, with the Earth’s climate being one of the most important examples. By systematically tuning the degree of tetrahedrality, we succeed in continuously interpolating between water-like behavior and simple liquid-like behavior. Our approach reveals what physical factors make water so anomalous and special even compared with other tetrahedral liquids.

Tetrahedral interactions describe the behavior of the most abundant and technologically important materials on Earth, such as water, silicon, carbon, germanium, and countless others. Despite their differences, these materials share unique common physical behaviors, such as liquid anomalies, open crystalline structures, and extremely poor glass-forming ability at ambient pressure. To reveal the physical origin of these anomalies and their link to the shape of the phase diagram, we systematically study the properties of the Stillinger–Weber potential as a function of the strength of the tetrahedral interaction λ. We uncover a unique transition to a reentrant spinodal line at low values of λ, accompanied with a change in the dynamical behavior, from non-Arrhenius to Arrhenius. We then show that a two-state model can provide a comprehensive understanding on how the thermodynamic and dynamic anomalies of this important class of materials depend on the strength of the tetrahedral interaction. Our work establishes a deep link between the shape of the phase diagram and the thermodynamic and dynamic properties through local structural ordering in liquids and hints at why water is so special among all substances.

Triplet correlations dominate the transition from simple to tetrahedral liquids

The total, triplet, and pair contributions to the entropy with increasing tetrahedrality are mapped out for the Stillinger-Weber liquids to demonstrate the qualitative and quantitative differences between triplet-dominated, tetrahedral liquids and pair-dominated, simple liquids with regard to supercooling and crystallization. The heat capacity anomaly of tetrahedral liquids originates in local ordering due to both pair and triplet correlations. The results suggest that structural correlations can be directly related to thermodynamic anomalies, phase changes, and self-assembly in other atomic and colloidal fluids.