On the Different Faces of the Supercritical Phase of Water at a Near-Critical Temperature: Pressure-Induced Structural Transitions Ranging from a Gaslike Fluid to a Plastic Crystal Polymorph

Substantial breakdown of the hydrogen-bonding network, local density inhomogeneities and fluid-liquid structural transitions in supercritical octanol-1: A molecular dynamics investigation

Molecular dynamics simulations have been employed to explore the hydrogen-bonding structure and dynamics in supercritical octanol-1 at a near-critical temperature and up to high densities and pressures. A substantial breakdown of the hydrogen-bonding network when going from ambient-liquid to supercritical conditions is revealed. The fraction of the non-hydrogen bonded molecules significantly increases in supercritical octanol-1, and a substantial decrease in the intermittent hydrogen-bond lifetime is observed. This behavior is also reflected on the maximum local density augmentation, which is comparable to the values obtained for non-polar and non-hydrogen bonded fluids. The existence of a structural transition from an inhomogeneous fluid phase to a soft-liquid one at densities higher than 2.0 ρc is also revealed. At higher densities, a significant change in the reorientational relaxation process is observed, reflected on the significant increase in the ratio of the Legendre reorientational times τ1R/τ2R. The latter becomes much higher than the value predicted by the Debye model of diffusive reorientation and the corresponding ratio for ambient liquid octanol-1. The non-polar tail of octanol-1 under supercritical conditions reorients more slowly in comparison with the polar tail. Interestingly, the opposite behavior is observed for the ambient liquid, further verifying the strong effect of the breakdown of the hydrogen bonding network on the properties of supercritical octanol-1. In accordance with the above-mentioned findings, the static dielectric constant of supercritical octanol-1 is very low even at high densities and pressures, comparable to the values obtained for non-polar and non-hydrogen bonded fluids.

Widom Line in Supercritical Water in Terms of Changes in Local Structure: Theoretical Perspective

Performing molecular dynamics simulations with the TIP4P/2005 water model along 9 isobars (from 175 to 375 bar) in the temperature range between 300 and 1100 K, we have found that the loci of the extrema in the rate of change of specific structural properties can be used to define purely structure-based Widom lines. We have examined several parameters that describe the local structure of water, such as the tetrahedral arrangement, nearest neighbor distance, local density around the molecules, and the size of the largest dense domain. The last two parameters were determined using the Voronoi polyhedral and density-based spatial clustering of applications with noise methods, respectively. By analyzing the moments of the associated distributions, we show that along a given isobar, the temperature at which we observe a maximum in the fluctuation, the rate of change of the average values, or in the skewness values unambiguously determines the Widom line that is in agreement with the experimentally detected, thermodynamic response function-based ones.

Electrolyte Permeability in Plastic Ice VII

Deep brines in water-rich planets form when electrolytes diffuse from the rocky interior through layers of thick dense ice such as ice VII and the hypothesized plastic ice VII. We perform classical molecular dynamics simulations of Li⁺, Na⁺, and K⁺ alkali ions and F^- and Cl^- halide ions in plastic ice VII at conditions similar to water-rich super-Earths, icy moons, and ocean worlds. We find that plastic ice VII is permeable to electrolytes on geological timescales. Diffusion occurs via jumps between adjacent voids in the bcc crystal structure and is governed by molecular rotations. An exception to this mechanism is Na⁺ which, at variance with other ions, can substitute water molecules on lattice positions. The bulk modulus of pristine plastic ice VII is dependent on the pace of molecular rotations: when the rotations are slow, the bulk modulus is 1 order of magnitude lower compared to the bulk modulus at conditions of fast rotations, hence providing direct evidence of the role of molecular rotations in determining elastic properties. Electrolytes affect the bulk modulus only at high-concentration conditions and slow molecular rotations. Our results show that plastic ice VII may facilitate the development of brines in water-rich planets and ocean worlds, with a clear significance for their potential to support exobiology and for the chemical evolution of their aqueous reservoirs.

Supercritical fluids behave as complex networks

Supercritical fluids play a key role in environmental, geological, and celestial processes, and are of great importance to many scientific and engineering applications. They exhibit strong variations in thermodynamic response functions, which has been hypothesized to stem from the microstructural behavior. However, a direct connection between thermodynamic conditions and the microstructural behavior, as described by molecular clusters, remains an outstanding issue. By utilizing a first-principles-based criterion and self-similarity analysis, we identify energetically localized molecular clusters whose size distribution and connectivity exhibit self-similarity in the extended supercritical phase space. We show that the structural response of these clusters follows a complex network behavior whose dynamics arises from the energetics of isotropic molecular interactions. Furthermore, we demonstrate that a hidden variable network model can accurately describe the structural and dynamical response of supercritical fluids. These results highlight the need for constitutive models and provide a basis to relate the fluid microstructure to thermodynamic response functions.

Molecular rotations trigger a glass-to-plastic fcc heterogeneous crystallization in high-pressure water

We report a molecular dynamics study of the heterogeneous crystallization of high-pressure glassy water using (plastic) ice VII as a substrate. We focus on the thermodynamic conditions P ∈ [6-8] GPa and T ∈ [100-500] K, at which (plastic) ice VII and glassy water are supposed to coexist in several (exo)planets and icy moons. We find that (plastic) ice VII undergoes a martensitic phase transition to a (plastic) fcc crystal. Depending on the molecular rotational lifetime τ, we identify three rotational regimes: for τ > 20 ps, crystallization does not occur; for τ ∼ 15 ps, we observe a very sluggish crystallization and the formation of a considerable amount of icosahedral environments trapped in a highly defective crystal or in the residual glassy matrix; and for τ < 10 ps, crystallization takes place smoothly, resulting in an almost defect-free plastic fcc solid. The presence of icosahedral environments at intermediate τ is of particular interest as it shows that such a geometry, otherwise ephemeral at lower pressures, is, indeed, present in water. We justify the presence of icosahedral structures based on geometrical arguments. Our results represent the first study of heterogeneous crystallization occurring at thermodynamic conditions of relevance for planetary science and unveil the role of molecular rotations in achieving it. Our findings (i) show that the stability of plastic ice VII, widely reported in the literature, should be reconsidered in favor of plastic fcc, (ii) provide a rationale for the role of molecular rotations in achieving heterogeneous crystallization, and (iii) represent the first evidence of long-living icosahedral structures in water. Therefore, our work pushes forward our understanding of the properties of water.

Fingerprints of the Crossing of the Frenkel and Melting Line on the Properties of High-Pressure Supercritical Water

Using molecular dynamics simulations in combination with the two-phase thermodynamic model, we reveal novel characteristic fingerprints of the crossing of the Frenkel and melting line on the properties of high-pressure water at a near-critical temperature (1.03T_c). The crossing of the Frenkel line at about 1.17 GPa is characterized by a crossover in the rotational and translational entropy ratio S^rot/S^trans, indicating a change in the coupling between translational and rotational motions which is also reflected in the shape of the rotational density of states. The observed isosbestic points in the translational and rotational density of states are also blue-shifted at density and pressure conditions higher than the ones corresponding to the Frenkel line. The first-order phase transition from a rigid liquid to a face-centered cubic plastic crystal phase at about 8.5 GPa is reflected in the discontinuous changes in the translational and rotational entropy, particularly in the significant increase of the ratio S^rot/S^trans. A noticeable discontinuous increase of the dielectric constant has also been revealed when crossing this melting line, which is attributed to the different arrangement of the water molecules in the plastic crystal phase. The reorientational dynamics in the plastic crystal phase is faster in comparison with the "rigid" liquid-like phase, but it remains unchanged upon a further pressure increase in the range of 8.5-11 GPa.

The Phonon Theory of Liquids and Biological Fluids: Developments and Applications

Among the three basic states of matter (solid, liquid, and gas), the liquid state has always eluded general theoretical approaches for describing liquid energy and heat capacity. In this Viewpoint, we derive the phonon theory of liquids and biological fluids stemming from Frenkel's microscopic picture of the liquid state. Specifically, the theory predicts the existence of phonon gaps in vibrational spectra of liquids and a thermodynamic boundary in the supercritical state. Direct experimental evidence reaffirming these theoretical predictions was achieved through a combination of techniques using static compression X-ray diffraction and inelastic X-ray scattering on deeply supercritical argon in a diamond anvil cell. Furthermore, these findings inspired and then led to the discovery of phonon gaps in liquid crystals (mesogens), block copolymers, and biological membranes. Importantly, phonon gaps define viscoelastic crossovers in cellular membranes responsible for lipid self-diffusion, lateral molecular-level stress propagation, and passive transmembrane transport of small molecules and solutes. Finally, molecular interactions mediated by external stimuli result in synaptic activity controlling biological membranes' plasticity resulting in learning and memory. Therefore, we also discuss learning and memory effects─equally important for neuroscience as well as for the development of neuromorphic devices─facilitated in biological membranes by external stimuli.

Temperature- and pressure-dependence of the hydrogen bond network in plastic ice VII

We model, via classical molecular dynamics simulations, the plastic phase of ice VII across a wide range of the phase diagram of interest for planetary investigations. Although structural and dynamical properties of plastic ice VII are mostly independent on the thermodynamic conditions, the hydrogen bond network (HBN) acquires a diverse spectrum of topologies distinctly different from that of liquid water and of ice VII simulated at the same pressure. We observe that the HBN topology of plastic ice carries some degree of similarity with the crystal phase, stronger at thermodynamic conditions proximal to ice VII, and gradually lessening upon approaching the liquid state. Our results enrich our understanding of the properties of water at high pressure and high temperature, and may help in rationalizing the geology of