Widom line and dynamical crossovers as routes to understand supercritical water

Critical Threshold for Bubble-like Nucleation during Pseudoboiling at Supercritical Pressures

Supercritical pseudoboiling was proposed in the 1950s-1960s. Recently, evaporation-like and boiling-like heat transfer have been directly observed in macroscopic scales, and the contribution of pseudoboiling to the total heat transfer rate has been quantitatively characterized experimentally. Here, we explore the critical threshold to generate a bubble-like nucleus at supercritical pressure at the atomic scale, characterized by the total energy (Te = Ke + Pe, where Ke and Pe are kinetic energy and potential energy, respectively). Molecular dynamics simulations are performed, including an argon fluid box heated by a solid wall having its temperature above the fluid temperature. The fluid pressure is controlled by a movable piston wall opposite the heating wall. The effects of pressure, nonuniform heating, and surface wettability on pseudoboiling are investigated. It is found that the criterion Te > 0 should be satisfied for subcritical boiling, matching that reported previously. The criterion for supercritical pseudoboiling was newly obtained such that Te > 0.012 eV at 8 MPa for argon, but the threshold increases as pressure increases. Nonuniform heating and surface wettability do not affect the critical threshold of Te for bubble-like nucleation but affect the location of the initially generated bubble-like nucleus and the stabilized pseudofilm or pseudonucleate heat transfer modes, where the former is similar to (vapor) film boiling and the latter is similar to nucleate boiling at subcritical pressure. Because pseudoboiling occurs without surface tension at supercritical pressure, we observe that the bubble-like structure may not display a perfectly smooth gas-liquid interface but may display an irregular pattern instead. Our work explains pseudoboiling from the viewpoint of the competition between kinetic energy and potential energy and presents a link regarding boiling in the two domains of subcritical pressure and supercritical pressure.

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.

Crossover from gas-like to liquid-like molecular diffusion in a simple supercritical fluid

According to textbooks, no physical observable can be discerned allowing to distinguish a liquid from a gas beyond the critical point. Yet, several proposals have been put forward challenging this view and various transition boundaries between a gas-like and a liquid-like behaviour, including the so-called Widom and Frenkel lines, and percolation line, have been suggested to delineate the supercritical state space. Here we report observation of a crossover from gas-like (Gaussian) to liquid-like (Lorentzian) self-dynamic structure factor by incoherent quasi-elastic neutron scattering measurements on supercritical fluid methane as a function of pressure, along the 200 K isotherm. The molecular self-diffusion coefficient was derived from the best Gaussian (at low pressures) or Lorentzian (at high pressures) fits to the neutron spectra. The Gaussian-to-Lorentzian crossover is progressive and takes place at about the Widom line intercept (59 bar). At considerably higher pressures, a liquid-like jump diffusion mechanism properly describes the supercritical fluid on both sides of the Frenkel line. The present observation of a gas-like to liquid-like crossover in the self dynamics of a simple supercritical fluid confirms emerging views on the unexpectedly complex physics of the supercritical state, and could have planet-wide implications and possible industrial applications in green chemistry.

Thermodynamic crossovers in supercritical fluids

Can liquid-like and gas-like states be distinguished beyond the critical point, where the liquid-gas phase transition no longer exists and conventionally only a single supercritical fluid phase is defined? Recent experiments and simulations report strong evidence of dynamical crossovers above the critical temperature and pressure. Despite using different criteria, many existing theoretical explanations consider a single crossover line separating liquid-like and gas-like states in the supercritical fluid phase. We argue that such a single-line scenario is inconsistent with the supercritical behavior of the Ising model, which has two crossover lines due to its symmetry, violating the universality principle of critical phenomena. To reconcile the inconsistency, we define two thermodynamic crossover lines in supercritical fluids as boundaries of liquid-like, indistinguishable, and gas-like states. Near the critical point, the two crossover lines follow critical scalings with exponents of the Ising universality class, supported by calculations of theoretical models and analyses of experimental data from the standard database. The upper line agrees with crossovers independently estimated from the inelastic X-ray scattering data of supercritical argon, and from the small-angle neutron scattering data of supercritical carbon dioxide. The lower line is verified by the equation of states for the compressibility factor. This work provides a fundamental framework for understanding supercritical physics in general phase transitions.

Liquid-Liquid Crossover in Water Model: Local Structure vs Kinetics of Hydrogen Bonds

In equilibrium and supercooled liquids, polymorphism is manifested by thermodynamic regions defined in the phase diagram, which are predominantly of different short- and medium-range order (local structure). It is found that on the phase diagram of the water model, the thermodynamic region corresponding to the equilibrium liquid phase is divided by a line of the smooth liquid-liquid crossover. In the case of the water model TIP4P/2005, this crossover is revealed by various local order parameters and corresponds to pressures on the order of 3150 ± 350 atm at ambient temperature. In the vicinity of the crossover, the dynamics of water molecules change significantly, which is reflected, in particular, in the fact that the self-diffusion coefficient reaches its maximum values. In addition, changes in the structure also manifest themselves in changes in the kinetics of hydrogen bonding, which are captured by values of such quantities as the average lifetime of hydrogen bonding, the average lifetimes of different local coordination numbers, and the frequencies of changes in different local coordination numbers. An interpretation of the hydrogen bond kinetics in terms of the free energy landscape concept in the space of possible coordination numbers is proposed.

Complex phase diagram and supercritical matter

The supercritical region is often described as uniform with no definite transitions. The distinct behaviors of the matter therein, e.g., as liquidlike and gaslike, however, suggest "supercritical boundaries." Here we provide a mathematical description of these phenomena by revisiting the Yang-Lee theory and introducing a complex phase diagram, specifically a four-dimensional (4D) one with complex T and p. While the traditional 2D phase diagram with real temperature T and pressure p values (the physical plane) lacks Lee-Yang (LY) zeros beyond the critical point, preventing the occurrence of criticality, the off-plane zeros in this 4D scenario still induce critical anomalies in various physical properties. This relationship is evidenced by the correlation between the Widom line and LY edges in van der Waals, 2D Ising model, and water. The diverged supercritical boundaries manifest the high-dimensional feature of the phase diagram: e.g., when LY zeros of complex T or p are projected onto the physical plane, boundaries defined by isobaric heat capacity C_{p} or isothermal compression coefficient K_{T} emanates. These results demonstrate the incipient phase transition nature of the supercritical matter.

Diffusion of hydrocarbons diluted in supercritical carbon dioxide

Mutual diffusion of six hydrocarbons (methane, ethane, isobutane, benzene, toluene or naphthalene) diluted in supercritical carbon dioxide ([Formula: see text]) is studied by molecular dynamics simulation near the Widom line, i.e., in the temperature range from 290 to 345 K along the isobar 9 MPa. The [Formula: see text] + aromatics mixtures are additionally sampled at 10 and 12 MPa and an experimental database with Fick diffusion coefficient data for those systems is provided. Taylor dispersion experiments of [Formula: see text] with benzene, toluene, n-dodecane and 1,2,3,4-tetrahydronaphthalene are conducted along the [Formula: see text] 10 MPa isobar. Maxwell-Stefan and Fick diffusion coefficients are analyzed, together with the thermodynamic factor that relates them. It is found that the peculiar behavior of the Fick diffusion coefficient of some [Formula: see text] mixtures in the extended critical region is a consequence of the thermodynamic factor minimum due to pronounced clustering on the molecular scale. Further, the strong dependence of the Fick diffusion coefficient on the molecular mass of the solute as well as the breakdown of the Stokes-Einstein relation near the Widom line are confirmed. Eleven correlations for the prediction of the Fick diffusion coefficient of [Formula: see text] mixtures are assessed. An alternative two-step approach for the prediction of the infinite dilution Fick diffusion coefficient of supercritical [Formula: see text] mixtures is proposed. It requires only the state point in terms of temperature and pressure (or density) as well as the molecular solute mass as input parameters. First, entropy scaling is applied to estimate the self-diffusion coefficient of [Formula: see text]. Subsequently, this coefficient is used to determine the infinite dilution Fick diffusion coefficient of the mixture, based on the finding that these two diffusion coefficients exhibit a linear relationship, where the slope depends only on the molecular solute mass.

Liquid Water: A Single Approach to Its Two Continuous Phase Transitions

In this paper, we show that both continuous phase transitions of liquid water, the liquid-gas and the liquid-liquid, can be articulated within a single thermodynamic analytical formalism. This result follows from a combination of the two-liquid model (TLM), recently confirmed for water, with the idea of a thermal-dependent excluded volume, v_e, concept introduced by van der Waals, in his famous state equation. Starting from the fundamentals of thermodynamics, it will be shown that the TLM naturally leads to the idea of an extensive thermal-dependent v_e that acts as a parameter of the sample thermodynamic potentials. This procedure effectively separates the thermodynamics of the system into two parts: the first concerns the clusters' thermodynamics, taken as wandering particles, and the second concerns the thermal behavior of its internal structure (geometry and number of particles). From this result, we demonstrate that the condition of mechanical instability leads to not one but two critical points, each happening in one of the above-described parts of the system.

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.

Numerical Investigation of Heat Transfer Characteristics of scCO2 Flowing in a Vertically-Upward Tube with High Mass Flux

In this work, the heat transfer characteristics of supercritical pressure CO₂ in vertical heating tube with 10 mm inner diameter under high mass flux were investigated by using an SST k-ω turbulent model. The influences of inlet temperature, heat flux, mass flux, buoyancy and flow acceleration on the heat transfer of supercritical pressure CO₂ were discussed. Our results show that the buoyancy and flow acceleration effect based on single phase fluid assumption fail to explain the current simulation results. Here, supercritical pseudo-boiling theory is introduced to deal with heat transfer of scCO₂. scCO₂ is treated to have a heterogeneous structure consisting of vapor-like fluid and liquid-like fluid. A physical model of scCO₂ heat transfer in vertical heating tube was established containing a gas-like layer near the wall and a liquid-like fluid layer. Detailed distribution of thermophysical properties and turbulence in radial direction show that scCO₂ heat transfer is greatly affected by the thickness of gas-like film, thermal properties of gas-like film and turbulent kinetic energy in the near-wall region. Buoyancy parameters Bu < 10⁻⁵, Bu* < 5.6 × 10⁻⁷ and flow acceleration parameter K_v < 3 × 10⁻⁶ in this paper, which indicate that buoyancy effect and flow acceleration effect has no influence on heat transfer of scCO₂ under high mass fluxes. This work successfully explains the heat transfer mechanism of supercritical fluid under high mass flux.

Universal Two-Component Dynamics in Supercritical Fluids

Despite the technological importance of supercritical fluids, controversy remains about the details of their microscopic dynamics. In this work, we study four supercritical fluid systems—water, Si, Te, and Lennard-Jones fluid—via classical molecular dynamics simulations. A universal two-component behavior is observed in the intermolecular dynamics of these systems, and the changing ratio between the two components leads to a crossover from liquidlike to gaslike dynamics, most rapidly around the Widom line. We find evidence to connect the liquidlike component dominating at lower temperatures with intermolecular bonding and the component prominent at higher temperatures with free-particle, gaslike dynamics. The ratio between the components can be used to describe important properties of the fluid, such as its self-diffusion coefficient, in the transition region. Our results provide an insight into the fundamental mechanism controlling the dynamics of supercritical fluids and highlight the role of spatiotemporally inhomogeneous dynamics even in thermodynamic states where no large-scale fluctuations exist in the fluid.

Advances in the study of supercooled water

In this review, we report recent progress in the field of supercooled water. Due to its uniqueness, water presents numerous anomalies with respect to most simple liquids, showing polyamorphism both in the liquid and in the glassy state. We first describe the thermodynamic scenarios hypothesized for the supercooled region and in particular among them the liquid-liquid critical point scenario that has so far received more experimental evidence. We then review the most recent structural indicators, the two-state model picture of water, and the importance of cooperative effects related to the fact that water is a hydrogen-bonded network liquid. We show throughout the review that water's peculiar properties come into play also when water is in solution, confined, and close to biological molecules. Concerning dynamics, upon mild supercooling water behaves as a fragile glass former following the mode coupling theory, and it turns into a strong glass former upon further cooling. Connections between the slow dynamics and the thermodynamics are discussed. The translational relaxation times of density fluctuations show in fact the fragile-to-strong crossover connected to the thermodynamics arising from the existence of two liquids. When considering also rotations, additional crossovers come to play. Mobility-viscosity decoupling is also discussed in supercooled water and aqueous solutions. Finally, the polyamorphism of glassy water is considered through experimental and simulation results both in bulk and in salty aqueous solutions. Grains and grain boundaries are also discussed.

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

The present study reports a systematic analysis of a wide variety of structural, thermodynamic, and dynamic properties of supercritical water along the near-critical isotherm of T = 1.03T_c and up to extreme pressures, using molecular dynamics and Monte Carlo simulations. The methodology employed provides solid evidence about the existence of a structural transition from a liquidlike fluid to a compressed, tightly packed liquid, in the density and pressure region around 3.4ρ_c and 1.17 GPa, introducing an alternative approach to locate the crossing of the Frenkel line. Around 8.5 GPa another transition to a face-centered-cubic plastic crystal polymorph with density 5.178ρ_c is also observed, further confirmed by Gibbs free energy calculations using the two-phase thermodynamic model. The isobaric heat capacity maximum, closely related to the crossing of the Widom line, has also been observed around 0.8ρ_c, where the local density augmentation is also maximized. Another structural transition has been observed at 0.2ρ_c, related to the transformation of the fluid to a dilute gas at lower densities. These findings indicate that a near-critical isotherm can be divided into different domains where supercritical water exhibits distinct behavior, ranging from a gaslike one to a plastic crystal one.

Towards Understanding the Structure of Subcritical and Transcritical Liquid–Gas Interfaces Using a Tabulated Real Fluid Modeling Approach

A fundamental understanding and simulation of fuel atomization, phase transition, and mixing are among the topics researchers have struggled with for decades. One of the reasons for this is that the accurate, robust, and efficient simulation of fuel jets remains a challenge. In this paper, a tabulated multi-component real-fluid model (RFM) is proposed to overcome most of the limitations and to make real-fluid simulations affordable. Essentially, a fully compressible two-phase flow and a diffuse interface approach are used for the RFM model, which were implemented in the CONVERGE solver. PISO and SIMPLE numerical schemes were modified to account for a highly coupled real-fluid tabulation approach. These new RFM model and numerical schemes were applied to the simulation of different fundamental 1-D, 2-D, and 3-D test cases to better understand the structure of subcritical and transcritical liquid–gas interfaces and to reveal the hydro-thermodynamic characteristics of multicomponent jet mixing. The simulation of a classical cryogenic injection of liquid nitrogen coaxially with a hot hydrogen jet is performed using thermodynamic tables generated by two different equations of state: Peng–Robinson (PR) and Soave–Redlich–Kwong (SRK). The numerical results are finally compared with available experimental data and published numerical studies with satisfactory agreement.

Phase distribution including a bubblelike region in supercritical fluid

Pseudoboiling in supercritical fluid (SF) has been paid great attention in recent years. Available works mainly focus on thermodynamics analysis. Fewer studies were reported on the spatial time phase distribution. Here, SF is investigated in a multiphase fluid framework using molecular dynamics (MD) simulations. A simulation box contains 10 976 argon atoms, with periodic boundary conditions applied on all the box surfaces. Pressure and temperature are well controlled. Based on MD simulation results, an onset pseudoboiling temperature T^{-} and a termination pseudoboiling temperature T^{+} are defined using the neighboring molecules method, the radial distribution function method, and the two-body excess entropy method. The two transition temperatures divide the whole phase diagram into three regimes of liquidlike, two-phase-like (TPL), and gaslike, and the MD determined T^{-} and T^{+} well matched the thermodynamics-determined values. In the TPL regime, nanovoids are observed to have two distinct characteristics: (1) Particles are sparsely distributed to have gas density inside the void, but are densely populated to have liquid density outside the void. (2) Voids have a curved interface. These characteristics are very similar to bubble characteristics in subcritical pressure. Hence, voids in the supercritical state are called "bubblelike" in this paper. Nonlinear dynamics demonstrates chaotic behavior in the TPL regime, similar to the two-phase regime in the subcritical domain. The above findings give strong evidence that SF in the TPL regime consists of a mixture of bubblelike voids and surrounding liquids. Our work highlights the multiphase feature of a SF, hence, the well-established multiphase theory in subcritical pressures can be introduced to handle the complex SF.

Dynamical crossover and its connection to the Widom line in supercooled TIP4P/Ice water

We perform molecular dynamics simulations with the TIP4P/Ice water model to characterize the relationship between dynamics and thermodynamics of liquid water in the supercooled region. We calculate the relevant properties of the phase diagram, and we find that TIP4P/Ice presents a retracing line of density maxima, similar to what was previously found for atomistic water models and models of other tetrahedral liquids. For this model, a liquid-liquid critical point between a high-density liquid and a low-density liquid was recently found. We compute the lines of the maxima of isothermal compressibility and the minima of the coefficient of thermal expansion in the one phase region, and we show that these lines point to the liquid-liquid critical point while collapsing on the Widom line. This line is the line of the maxima of correlation length that emanates from a second order critical point in the one phase region. Supercooled water was found to follow mode coupling theory and to undergo a transition from a fragile to a strong behavior right at the crossing of the Widom line. We find here that this phenomenology also happens for TIP4P/Ice. Our results appear, therefore, to be a general characteristic of supercooled water, which does not depend on the interaction potential used, and they reinforce the idea that the dynamical crossover from a region where the relaxation mechanism is dominated by cage relaxation to a region where cages are frozen and hopping dominates is correlated in water to a phase transition between a high-density liquid and a low-density liquid.

Diffusion of the carbon dioxide-ethanol mixture in the extended critical region

The effect of traces of ethanol in supercritical carbon dioxide on the mixture's thermodynamic properties is studied by molecular simulations and Taylor dispersion measurements. This mixture is investigated along the isobar p = 10 MPa in the temperature range between T = 304 and 343 K. Along this path, the mixture undergoes two transitions: First, the Widom line is crossed, marking the transition from liquid-like to gas-like conditions. A second transition occurs from the supercritical gas-like domain to a subcritical gas. The Widom line crossover entails inflection points for most of the studied properties, i.e. density, enthalpy, shear viscosity, Maxwell-Stefan and intradiffusion coefficients. On the other hand, the transition between the super- and subcritical regions is found to be generally smooth, an observation that is qualitatively confirmed by experimental Taylor dispersion measurements. Dedicated atomistic simulations show the presence of microheterogeneities due to ethanol self-association along the investigated path, which lead to the mixture's anomalous behavior in its extended critical region.

Thermodynamics and Dynamics of Supercritical Water Pseudo-Boiling

Supercritical fluid pseudo-boiling (PB), recently brought to the attention of the scientific community, is the phenomenon occurring when fluid changes its structure from liquid-like (LL) to gas-like (GL) states across the Widom line. This work provides the first quantitative analysis on the thermodynamics and the dynamics of water's PB, since the understanding of this phase transition is mandatory for the successful implementation of technologies using supercritical water (scH2O) for environmental, energy, and nanomaterial applications. The study combines computational techniques with in situ neutron imaging measurements. The results demonstrate that, during isobaric heating close to the critical point, while water density drops by a factor of three in the PB transitional region, the system needs >16 times less energy to increase its temperature by 1 K than to change its structure from LL to GL phase. Above the PB-Widom line, the structure of LL water consists mainly of tetramers and trimers, while below the line mostly dimers and monomers form in the GL phase. At atomic level, the PB dynamics are similar to those of the subcritical water vaporization. This fundamental knowledge has great impact on water science, as it helps to establish the structure-properties relationship of scH2O.

Enhancement and maximum in the isobaric specific-heat capacity measurements of deeply supercooled water using ultrafast calorimetry

Knowledge of the temperature dependence of the isobaric specific heat (C_p) upon deep supercooling can give insights regarding the anomalous properties of water. If a maximum in C_p exists at a specific temperature, as in the isothermal compressibility, it would further validate the liquid-liquid critical point model that can explain the anomalous increase in thermodynamic response functions. The challenge is that the relevant temperature range falls in the region where ice crystallization becomes rapid, which has previously excluded experiments. Here, we have utilized a methodology of ultrafast calorimetry by determining the temperature jump from femtosecond X-ray pulses after heating with an infrared laser pulse and with a sufficiently long time delay between the pulses to allow measurements at constant pressure. Evaporative cooling of ∼15-µm diameter droplets in vacuum enabled us to reach a temperature down to ∼228 K with a small fraction of the droplets remaining unfrozen. We observed a sharp increase in C_p, from 88 J/mol/K at 244 K to about 218 J/mol/K at 229 K where a maximum is seen. The C_p maximum is at a similar temperature as the maxima of the isothermal compressibility and correlation length. From the C_p measurement, we estimated the excess entropy and self-diffusion coefficient of water and these properties decrease rapidly below 235 K.

Two-Component Dynamics and the Liquidlike to Gaslike Crossover in Supercritical Water

Molecular-scale dynamics in sub- to supercritical water is studied with inelastic x-ray scattering and molecular dynamics simulations. The obtained longitudinal current correlation spectra can be decomposed into two main components: a low-frequency (LF), gaslike component and a high-frequency (HF) component arising from the O-O stretching mode between hydrogen-bonded molecules, reminiscent of the longitudinal acoustic mode in ambient water. With increasing temperature, the hydrogen-bond network diminishes and the spectral weight shifts from HF to LF, leading to a transition from liquid- to gaslike dynamics with rapid changes around the Widom line.

Could Chemical Reaction at the Molecular Level Show Distinction between Two Liquid-Water States? Study of the Excited-State Water-Catalyzed Proton Transfer Reaction Provides a Clue

The two liquid-water states, which lead to some anomalies when temperature crosses over 50 ± 10 °C at the atmospheric pressure, have been continuously catching popular attention. In this study, using the excited-state proton transfer (ESPT) catalyzed by water molecules as a prototypical reaction, we demonstrate that the kinetics of ESPT indeed is influenced by the two liquid-water states. In the water-catalyzed ESPT of 3-cyano-7-azaindole (3CAI), a repetitive and comprehensive temperature-dependent study of ESPT in H₂O from 0 to 90 °C shows anomalous behavior. The plot of the logarithm of ESPT rate constant as a function of inverse of absolute temperature deviates from a straight line. The convex-Arrhenius behavior manifests the activation free energy for water-assisted ESPT being dependent on temperature and hence the liquid water structure. To simplify the discussion, the plot is well fitted by using two straight lines that are crossed over in the vicinity of 40 °C. The free energy difference between water-solvated 3CAI and the 1:1 H₂O:3CAI complex is deduced to be 2.29 ± 0.04 and 1.96 ± 0.04 kcal·mol^-1 in the regions of 0-40 and 40-90 °C water, respectively, which also results in different frequency factors, i.e., the proton transfer/tunneling rates of (5.83 ± 0.36) × 10¹⁰ and (3.48 ± 0.27) × 10¹⁰ s^-1, respectively. In a qualitative manner, the results are then rationalized by the different types of H-bonding configuration as proposed for two liquid-water phases, rendering experimental evidence to support the different water phases in ambient temperatures at 1 bar.

Hydrogen Dynamics in Supercritical Water Probed by Neutron Scattering and Computer Simulations

In this work, an investigation of supercritical water is presented combining inelastic and deep inelastic neutron scattering experiments and molecular dynamics simulations based on a machine-learned potential of ab initio quality. The local hydrogen dynamics is investigated at 250 bar and in the temperature range of 553-823 K, covering the evolution from subcritical liquid to supercritical gas-like water. The evolution of libration, bending, and stretching motions in the vibrational density of states is studied, analyzing the spectral features by a mode decomposition. Moreover, the hydrogen nuclear momentum distribution is measured, and its anisotropy is probed experimentally. It is shown that hydrogen bonds survive up to the higher temperatures investigated, and we discuss our results in the framework of the coupling between intramolecular modes and intermolecular librations. Results show that the local potential affecting hydrogen becomes less anisotropic within the molecular plane in the supercritical phase, and we attribute this result to the presence of more distorted hydrogen bonds.

Supercritical Water is not Hydrogen Bonded

Thinking about water is inextricably linked to hydrogen bonds, which are highly directional in character and determine the unique structure of water, in particular its tetrahedral H-bond network. Here, we assess if this common connotation also holds for supercritical water. We employ extensive ab initio molecular dynamics simulations to systematically monitor the evolution of the H-bond network mode of water from room temperature, where it is the hallmark of its fluctuating three-dimensional network structure, to supercritical conditions. Our simulations reveal that the oscillation period required for H-bond vibrations to occur exceeds the lifetime of H-bonds in supercritical water by far. Instead, the corresponding low-frequency intermolecular vibrations of water pairs as seen in supercritical water are found to be well represented by isotropic van-der-Waals interactions only. Based on these findings, we conclude that water in its supercritical phase is not a H-bonded fluid.

Pronounced structural crossover in water at supercritical pressures

There have been ample studies of the many phases of H2O in both its solid and low pressure liquid states, and the transitions between them. Using molecular dynamics simulations we address the hitherto unexplored deeply supercritical pressures, where no qualitative transitions are thought to take place and where all properties are expected to vary smoothly. On the basis of these simulations we predict that water at supercritical pressures undergoes a structural crossover across the Frenkel line at pressures as high as 45 times the critical pressure. This provides a new insight into the water phase diagram and establishes a link between the structural and dynamical properties of supercritical water. Specifically, the crossover is demonstrated by a sharp and pronounced at low pressures, and smooth at high pressures, signified by changes in the pair distribution functions and local coordination which coincide with the dynamical transition (the loss of all oscillatory molecular motion) at the Frenkel line on the phase diagram.

Universality, Scaling, and Collapse in Supercritical Fluids

Supercritical fluid (SCF) is known to exhibit salient dynamic and thermodynamic crossovers and an inhomogeneous molecular distribution. However, the question as to what basic physics underlies these microscopic and macroscopic anomalies remains open. Here, using an order parameter extracted by machine learning, the fraction of gas-like (or liquid-like) molecules, we find simplicity and universality in SCF: First, all isotherms of a given fluid collapse onto a single master curve described by a scaling relation. The observed power law holds from the high-temperature and -pressure regime down to the critical point where it diverges. Second, phase diagrams of different compounds collapse onto their master curves by the same scaling exponent, thereby demonstrating a putative law of corresponding supercritical states in simple fluids. The reported results support a model of the SCF as a mixture of two interchangeable microstates, whose spatiotemporal dynamics gives rise to unique macroscopic properties.

Thermodynamics, dynamics, and structure of supercritical water at extreme conditions

Molecular dynamics (MD) simulations to understand the thermodynamic, dynamic, and structural changes in supercritical water across the Frenkel line and the melting line have been performed.

Supercritical water anomalies in the vicinity of the Widom line

Supercritical water is used in a variety of chemical and industrial applications. As a consequence, a detailed knowledge of the structure-properties correlations is of uttermost importance. Although supercritical water was considered as a homogeneous fluid, recent studies revealed an anomalous behaviour due to nanoscale density fluctuations (inhomogeneity). The inhomogeneity is clearly demarked through the Widom line (maxima in response factions) and drastically affect the properties. In the current study the physical properties of supercritical water have been determined by classical molecular dynamics simulations using a variety of polarized and polarizable interatomic potentials. Their validity which was not available at supercritical conditions has been assessed based on the ability to reproduce experimental data. Overall, the polarized TIP4P/2005 model accurately predicted the properties of water in both liquid-like and gas-like regions. All interatomic potentials captured the anomalous behaviour providing a direct evidence of molecular-scale inhomogeneity.

Two-dimensional infrared spectroscopy from the gas to liquid phase: density dependent J-scrambling, vibrational relaxation, and the onset of liquid character

Ultrafast 2DIR spectra and pump-probe responses of the N₂O ν₃ asymmetric stretch in SF₆ as a function of density from the gas to supercritical phase and liquid are reported. 2DIR spectra unequivocally reveal free rotor character at all densities studied in the gas and supercritical region. Analysis of the 2DIR spectra determines that J-scrambling or rotational relaxation in N₂O is highly efficient, occurring in ∼1.5 to ∼2 collisions with SF₆ at all non-liquid densities. In contrast, N₂O ν₃ vibrational energy relaxation requires ∼15 collisions, and complete vibrational equilibrium occurs on the ∼ns scale at all densities. An independent binary collision model is sufficient to describe these supercritical state point dynamics. The N₂O ν₃ in liquid SF₆ 2DIR spectrum shows no evidence of free rotor character or spectral diffusion. Using these 2DIR results, hindered rotor or liquid-like character is found in gas and all supercritical solutions for SF₆ densities ≥ρ* = 0.3, and increases with SF₆ density. 2DIR spectral analysis offers direct time domain evidence of critical slowing for SF₆ solutions closest to the critical point density. Applications of 2DIR to other high density and supercritical solution dynamics and descriptions are discussed.

Visualization of supercritical water pseudo-boiling at Widom line crossover

Supercritical water is a green solvent used in many technological applications including materials synthesis, nuclear engineering, bioenergy, or waste treatment and it occurs in nature. Despite its relevance in natural systems and technical applications, the supercritical state of water is still not well understood. Recent theories predict that liquid-like (LL) and gas-like (GL) supercritical water are metastable phases, and that the so-called Widom line zone is marking the crossover between LL and GL behavior of water. With neutron imaging techniques, we succeed to monitor density fluctuations of supercritical water while the system evolves rapidly from LL to GL as the Widom line is crossed during isobaric heating. Our observations show that the Widom line of water can be identified experimentally and they are in agreement with the current theory of supercritical fluid pseudo-boiling. This fundamental understanding allows optimizing and developing new technologies using supercritical water as a solvent.

Supercritical water exists in gas- and liquid-like forms, but these have not been distinguished yet at the macroscale. Here the authors investigate supercritical water interacting with microporous carbon by neutron imaging, and observe the coexistence of gas- and liquid-like states upon crossing the Widom line.

Ferromagnetism-induced phase separation in a two-dimensional spin fluid

We study the liquid-gas phase separation observed in a system of repulsive particles dressed with ferromagnetically aligning spins, a so-called "spin fluid." Microcanonical ensemble numerical simulations of finite-size systems reveal that magnetization sets in and induces a liquid-gas phase separation between a disordered gas and a ferromagnetic dense phase at low enough energies and large enough densities. The dynamics after a quench into the coexistence region show that the order parameter associated with the liquid-vapor phase separation follows an algebraic law with an unusual exponent, as it is forced to synchronize with the growth of the magnetization: this suggests that for finite size systems the magnetization sets in along a Curie line, which is also the gas-side spinodal line, and that the coexistence region ends at a tricritical point. This picture is confirmed at the mean-field level with different approximation schemes, namely, a Bethe lattice resolution and a virial expansion complemented by the introduction of a self-consistent Weiss-like molecular field. However, a detailed finite-size scaling analysis shows that in two dimensions the ferromagnetic phase escapes the Berezinskii-Kosterlitz-Thouless scenario and that the long-range order is not destroyed by the unbinding of topological defects. The Curie line thus becomes a magnetic crossover in the thermodynamic limit. Finally, the effects of the magnetic interaction range and those of the interaction softness are characterized within a mean-field semianalytical low-density approach.

Apparent power-law behavior of water's isothermal compressibility and correlation length upon supercooling

Apparent power-law analysis of water's isothermal compressibility and correlation length in the temperature range from 280 K to 229 K.