Van der Waals supercritical fluid: Exact formulas for special lines

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.

Boiling line and near-critical maxima of propane-nitrogen mixtures

It is well known that some thermodynamic quantities demonstrate maxima in the vicinity of a critical point. The lines of these maxima in the density-temperature or pressure-temperature planes are called "Widom lines." The behavior of Widom lines of one-component fluids has already been well studied in a number of papers by different authors. However, up to now the understanding of Widom lines in binary mixtures is still lacking. In this paper we study the boiling curve and the near-critical maxima of mixtures of nitrogen and propane by means of molecular dynamics simulation. We calculate the boiling curves and estimate the critical temperatures in a set of concentrations from pure nitrogen to pure propane. The influence of the composition of the mixture on the Widom lines of the system is evaluated. We find that the mixture of propane and nitrogen behaves as a type I mixture in the van Konynenburg-Scott classification, i.e., when the concentration is changed, the critical point and the corresponding Widom lines continuously shift in the density-temperature plane.

The Line of Ideal Isothermal Compressibility

We have considered the line along which the values of the isothermal compressibility of a system are the same as they would be for an ideal gas. It was called the κ_T line. Various substances and models have been studied with the use of the multi-parameter equations of states, implemented in REFPROP (Lemmon, E. W.; Bell, I. H.; Huber, M. L.; McLinden, M. O. NIST standard reference database 23: Reference fluid thermodynamic and transport properties-REFPROP, version 10.0 (National Institute of Standards and Technology, 2018)), and numerical simulations, respectively. It is shown that in the reduced variables, all κ_T lines can be described by the same relation following from the van der Waals equation. This similarity appears to be valid even for the systems with non-straight Zeno lines. Deviations within ∼5% are observed only near the binodals.

Frenkel Line in Nitrogen Terminates at the Triple Point

Recent studies on supercritical nitrogen revealed clear changes in structural markers and dynamical properties when the coordination number approaches its maximum value. The line in P and T space where these changes occur is referred to as the Frenkel line. Here, we qualitatively reproduce such changes in the supercritical regime using the popular "optimized potential for liquid simulation" (OPLS) classical force field for molecular dynamics. Unfortunately, at 160 K, OPLS nitrogen predicts sublimation rather than producing a liquid phase; therefore, we developed our own force field to achieve quantitative agreement with experimental data. We confirm the asymptotic behavior of the coordination number on crossing the Frenkel line and note an associated change in the diffusion constant, consistent with the non-rigid to rigid liquid-like character of the "transition". The simulations allow us to track the Frenkel line to subcritical temperatures and demonstrate that it terminates at the triple point. This establishes the experimentally measurable changes, which could unequivocally determine the Frenkel line in other systems.

Thermodynamic heterogeneity and crossover in the supercritical state of matter

A hallmark of a thermodynamic phase transition is the qualitative change of system thermodynamic properties such as energy and heat capacity. On the other hand, no phase transition is thought to operate in the supercritical state of matter and, for this reason, it was believed that supercritical thermodynamic properties vary smoothly and without any qualitative changes. Here, we perform extensive molecular dynamics simulations in a wide temperature range and find that a deeply supercritical state is thermodynamically heterogeneous, as witnessed by different temperature dependence of energy, heat capacity and its derivatives at low and high temperature. The evidence comes from three different methods of analysis, two of which are model-independent. We propose a new definition of the relative width of the thermodynamic crossover and calculate it to be in the fairly narrow relative range of 13%-20%. On the basis of our results, we relate the crossover to the supercritical Frenkel line.

Investigation concerning the uniqueness of separatrix lines separating liquidlike from gaslike regimes deep in the supercritical phase of water with a focus on Widom line concepts

The supercritical phase of fluids has long been known to feature significantly different liquidlike and gaslike regimes. However, it is textbook knowledge that the supercritical state is a homogeneous fluid phase where properties change continuously. Nevertheless, there has been an increasing amount of evidence published that suggests that there might exist a unique line that rigorously separates different regimes in supercritical phases, particularly in the case of water. Here, we use the quasiexact IAPWS95 equation of state to rigorously assess the macroscopic thermodynamic properties of supercritical water without invoking any water model or related approximations. We focus on how these properties change deep in the supercritical phase, in particular if they allow one to introduce a unique "thermodynamic separatrix." Our rigorous thermodynamic analysis, which relies exclusively on accurate experimental data, makes clear that there is no unique separatrix in real supercritical water-such as the recently much-invoked "Widom line." A comparison to the van der Waals equation of state reproduces qualitatively all our findings for real water, thereby suggesting that our analysis should be transferable to other fluids and critical points. Topological analysis of the H-bond network structure of supercritical water, as obtained from molecular-dynamics simulations using a standard water model, demonstrates that also the percolation line does not provide a meaningful separatrix to rigorously distinguish liquidlike from gaslike regimes.

Dynamics, thermodynamics and structure of liquids and supercritical fluids: crossover at the Frenkel line

We review recent work aimed at understanding dynamical and thermodynamic properties of liquids and supercritical fluids. The focus of our discussion is on solid-like transverse collective modes, whose evolution in the supercritical fluids enables one to discuss the main properties of the Frenkel line separating rigid liquid-like and non-rigid gas-like supercritical states. We subsequently present recent experimental evidence of the Frenkel line showing that structural and dynamical crossovers are seen at a pressure and temperature corresponding to the line as predicted by theory and modelling. Finally, we link dynamical and thermodynamic properties of liquids and supercritical fluids by the new calculation of liquid energy governed by the evolution of solid-like transverse modes. The disappearance of those modes at high temperature results in the observed decrease of heat capacity.

Widom line, dynamical crossover, and percolation transition of supercritical oxygen via molecular dynamics simulations

Supercritical oxygen, a cryogenic fluid, is widely used as an oxidizer in jet propulsion systems and is therefore of paramount importance in gaining physical insights into processes such as transcritical and supercritical vaporization. It is well established in the scientific literature that the supercritical state is not homogeneous but, in fact, can be demarcated into regions with liquid-like and vapor-like properties, separated by the "Widom line." In this study, we identified the Widom line for oxygen, constituted by the loci of the extrema of thermodynamic response functions (heat capacity, volumetric thermal expansion coefficient, and isothermal compressibility) in the supercritical region, via atomistic molecular dynamics simulations. We found that the Widom lines derived from these response functions all coincide near the critical point until about 25 bars and 15-20 K, beyond which the isothermal compressibility line begins to deviate. We also obtained the crossover from liquid-like to vapor-like behavior of the translational diffusion coefficient, shear viscosity, and rotational relaxation time of supercritical oxygen. While the crossover of the translational diffusion coefficient and shear viscosity coincided with the Widom lines, the rotational relaxation time showed a crossover that was largely independent of the Widom line. Further, we characterized the clustering behavior and percolation transition of supercritical oxygen molecules, identified the percolation threshold based on the fractal dimension of the largest cluster and the probability of finding a cluster that spans the system in all three dimensions, and found that the locus of the percolation threshold also coincided with the isothermal compressibility Widom line. It is therefore clear that supercritical oxygen is far more complex than originally perceived and that the Widom line, dynamical crossovers, and percolation transitions serve as useful routes to better our understanding of the supercritical state.

Classical and quantum many-body effects on the critical properties and thermodynamic regularities of silicon

Using molecular simulation, we determine the critical properties of Si as well as the loci for several remarkable thermodynamic contours spanning the supercritical region of the phase diagram. We consider a classical three-body potential as well as a quantum (tight-binding) many-body model, and determine the loci for the ideality contours, including the Zeno line and the H line of ideal enthalpy. The two strategies (classical or quantum) lead to strongly asymmetric binodals and to critical properties in good agreement with each other. The Zeno and H lines are found to remain linear over a wide temperature interval, despite the changes in electronic structure undergone by the fluid along these contours. We also show that the classical and quantum model yield markedly different results for the parameters defining the H line, the exponents for the power-laws underlying the line of minima for the isothermal enthalpy and for the density required to achieve ideal behavior, most notably for the enthalpy.

Liquid-liquid critical point in a simple analytical model of water

A statistical model for a simple three-dimensional Mercedes-Benz model of water was used to study phase diagrams. This model on a simple level describes the thermal and volumetric properties of waterlike molecules. A molecule is presented as a soft sphere with four directions in which hydrogen bonds can be formed. Two neighboring waters can interact through a van der Waals interaction or an orientation-dependent hydrogen-bonding interaction. For pure water, we explored properties such as molar volume, density, heat capacity, thermal expansion coefficient, and isothermal compressibility and found that the volumetric and thermal properties follow the same trends with temperature as in real water and are in good general agreement with Monte Carlo simulations. The model exhibits also two critical points for liquid-gas transition and transition between low-density and high-density fluid. Coexistence curves and a Widom line for the maximum and minimum in thermal expansion coefficient divides the phase space of the model into three parts: in one part we have gas region, in the second a high-density liquid, and the third region contains low-density liquid.

Crossover of collective modes and positive sound dispersion in supercritical state

Supercritical state has been viewed as an intermediate state between gases and liquids with largely unknown physical properties. Here, we address the important ability of supercritical fluids to sustain collective excitations. We directly study propagating modes on the basis of correlation functions calculated in molecular dynamics simulations and find that the supercritical system sustains propagating solid-like transverse modes below the Frenkel line but not above where there is one longitudinal mode only. Important thermodynamic implications of this finding are discussed. We directly detect positive sound dispersion (PSD) below the Frenkel line where transverse modes are operative and quantitatively explain its magnitude on the basis of transverse and longitudinal velocities. PSD disappears above the Frenkel line which therefore demarcates the supercritical phase diagram into two areas where PSD does and does not operate.

Impact of Friedel oscillations on vapor-liquid equilibria and supercritical properties in two and three dimensions

We determine the impact of the Friedel oscillations on the phase behavior, critical properties, and thermodynamic contours in films [two dimensions (2D)] and bulk phases [three dimensions (3D)]. Using expanded Wang-Landau simulations, we calculate the grand-canonical partition function and, in turn, the thermodynamic properties of systems modeled with a linear combination of the Lennard-Jones and Dzugutov potentials, weighted by a parameter X (0<X<1). Varying X allows us to control the height of the first Friedel oscillation and to provide a complete characterization of the effect of the metal-like character in the potential on the thermodynamic properties over a wide range of conditions. For 3D systems, we are able to show that the critical parameters exhibit a linear dependence on X and that the loci for the thermodynamic state points, for which the system shows the same compressibility factor or enthalpy as an ideal gas, are two straight lines spanning the subcritical and supercritical regions of the phase diagram for all X values. Reducing the dimensionality to 2D results in a loss of impact of the Friedel oscillation on the critical properties, as evidenced by the virtually constant critical density across the range of X values. Furthermore, our results establish that the straightness of the two ideality lines is retained in 2D and is independent from the height of the first Friedel oscillation in the potential.

Collective modes and thermodynamics of the liquid state

Strongly interacting, dynamically disordered and with no small parameter, liquids took a theoretical status between gases and solids with the historical tradition of hydrodynamic description as the starting point. We review different approaches to liquids as well as recent experimental and theoretical work, and propose that liquids do not need classifying in terms of their proximity to gases and solids or any categorizing for that matter. Instead, they are a unique system in their own class with a notably mixed dynamical state in contrast to pure dynamical states of solids and gases. We start with explaining how the first-principles approach to liquids is an intractable, exponentially complex problem of coupled non-linear oscillators with bifurcations. This is followed by a reduction of the problem based on liquid relaxation time τ representing non-perturbative treatment of strong interactions. On the basis of τ, solid-like high-frequency modes are predicted and we review related recent experiments. We demonstrate how the propagation of these modes can be derived by generalizing either hydrodynamic or elasticity equations. We comment on the historical trend to approach liquids using hydrodynamics and compare it to an alternative solid-like approach. We subsequently discuss how collective modes evolve with temperature and how this evolution affects liquid energy and heat capacity as well as other properties such as fast sound. Here, our emphasis is on understanding experimental data in real, rather than model, liquids. Highlighting the dominant role of solid-like high-frequency modes for liquid energy and heat capacity, we review a wide range of liquids: subcritical low-viscous liquids, supercritical state with two different dynamical and thermodynamic regimes separated by the Frenkel line, highly-viscous liquids in the glass transformation range and liquid-glass transition. We subsequently discuss the fairly recent area of liquid-liquid phase transitions, the area where the solid-like properties of liquids have become further apparent. We then discuss gas-like and solid-like approaches to quantum liquids and theoretical issues that are similar to the classical case. Finally, we summarize the emergent view of liquids as a unique system with a mixed dynamical state, and list several areas where interesting insights may appear and continue the extraordinary liquid story.

Dynamical crossover line in supercritical water

Dynamical crossover in water is studied by means of computer simulation. The crossover temperature is calculated from the behavior of velocity autocorrelation functions. The results are compared with experimental data. It is shown that the qualitative behavior of the dynamical crossover line is similar to the melting curve behavior. Importantly, the crossover line belongs to experimentally achievable (P, T) region which stimulates the experimental investigation in this field.

Reply to "Comment on 'Dynamic transition of supercritical hydrogen: Defining the boundary between interior and atmosphere in gas giants' "

We have recently proposed the dynamic transition of molecular hydrogen in gas giants around 10 GPa based on the recent understanding of a supercritical state related to the Frenkel line. In the preceding Comment, Bryk makes several remarks with regard to the Frenkel line and its relationship to the speed of sound. We disagree with this discussion and respond to it in this Reply.

Thermodynamic properties of supercritical carbon dioxide: Widom and Frenkel lines

Supercritical fluids are widely used in a number of important technological applications, yet the theoretical progress in the field has been rather moderate. Fairly recently, a new understanding of the liquidlike and gaslike properties of supercritical fluids has come to the fore, particularly with the advent of the Widom and Frenkel lines that aim to demarcate different physical properties on the phase diagram. Here, we report the results of a computational study of supercritical carbon dioxide, one of the most important fluids in the chemical industry. We study the response functions of CO_{2} in the supercritical state and calculate the locations of their maxima (Widom lines). We also report the preliminary calculations of the Frenkel line, the line of crossover of microscopic dynamics of particles. Our insights are relevant to physical processes in the atmosphere of Venus and its evolution.

How dimensionality changes the anomalous behavior and melting scenario of a core-softened potential system?

We present a computer simulation study of the phase diagram and anomalous behavior of two-dimensional (2D) and three-dimensional (3D) classical particles repelling each other through an isotropic core-softened potential. As in the analogous three-dimensional case, in 2D a reentrant-melting transition occurs upon compression under not too high pressure, along with a spectrum of thermodynamic and dynamic anomalies in the fluid phase. However, in two dimensions the order of the region of anomalous diffusion and the region of structural anomaly is inverted in comparison with the 3D case, where there exists a water-like sequence of anomalies, and has a silica-like sequence. In the low density part of the 2D phase diagram, melting is a continuous two-stage transition, with an intermediate hexatic phase. All available evidence supports the Kosterlitz-Thouless-Halperin-Nelson-Young (KTHNY) scenario for this melting transition. On the other hand, at high density part of the phase diagram one first-order transition takes place.

True Widom line for a square-well system

In the present paper we propose a van der Waals-like model that allows a purely analytical study of fluid properties including the equation of state, phase behavior, and supercritical fluctuations. We take a square-well system as an example and calculate its liquid-gas transition line and supercritical fluctuations. Employing this model allows us to calculate not only the thermodynamic response functions (isothermal compressibility βT, isobaric heat capacity CP, density fluctuations ζT, and thermal expansion coefficient αT), but also the correlation length in the fluid ξ. It is shown that the bunch of extrema widens rapidly upon departure from the critical point. It seems that the Widom line defined in this way cannot be considered as a real boundary that divides the supercritical region into gaslike and liquidlike regions. As it has been shown recently, a dynamic line on the phase diagram in the supercritical region, namely, the Frenkel line, can be used for this purpose.

Silicalike sequence of anomalies in core-softened systems

We present a simulation study of density, structural, and diffusion anomalies in a core-softened system, a remarkable model liquid that exhibits anomalous properties seen in tetrahedral liquids such as silica and water. It is widely believed that core-softened potentials demonstrate waterlike sequence of anomalies. Here, we show that the order of the region of anomalous diffusion and the regions of density and structural anomalies are inverted with increasing depth of the attractive part of the potential and have silicalike sequence. We also show that the Widom line slope is negative as in water.

Dynamics and Thermodynamics beyond the critical point

Sudden changes in the dynamical properties of a supercritical fluid model have been found as a function of pressure and temperature (T/T(c) = 2-5 and P/P(c) = 10-10(3)), striving with the notion of a single phase beyond the critical point established by thermodynamics. The sound propagation in the Terahertz frequency region reveals a sharp dynamic crossover between the gas like and the liquid like regimes along several isotherms, which involves, at sufficiently low densities, the interplay between purely acoustic waves and heat waves. Such a crossover allows one to determine a dynamic line in the phase diagram which exhibits a very tight correlation with a number of thermodynamic observables, showing that the supercritical state is remarkably more complex than thought so far.

Widom line and noise-power spectral analysis of a supercritical fluid

We have performed extensive molecular dynamics simulations to study noise-power spectra of density and potential energy fluctuations of a Lennard-Jones model of a fluid in the supercritical region. Emanating from the liquid-vapor critical point, there is a locus of isobaric specific heat maxima, called the Widom line, which is often regarded as an extension of the liquid-vapor coexistence line. Our simulation results show that the noise-power spectrum of the density fluctuations on the Widom line of the liquid-vapor transition exhibits three distinct 1/f^{γ} behaviors with exponents γ=0, 1.2, and 2, depending on the frequency f. We find that the intermediate frequency region with an exponent γ∼ 1 appears as the temperature approaches the Widom temperature from above or below. On the other hand, we do not find three distinct regions of 1/f^{γ} in the power spectrum of the potential energy fluctuations on the Widom line. Furthermore, we find that the power spectra of both the density and potential energy fluctuations at low frequency have a maximum on the Widom line, suggesting that the noise power can provide an alternative signature of the Widom line.

Riemannian geometry study of vapor-liquid phase equilibria and supercritical behavior of the Lennard-Jones fluid

The behavior of thermodynamic response functions and the thermodynamic scalar curvature in the supercritical region have been studied for a Lennard-Jones fluid based on a revised modified Benedict-Webb-Rubin equation of state. Response function extrema are sometimes used to estimate the Widom line, which is characterized by the maxima of the correlation lengths. We calculated the Widom line for the Lennard-Jones fluid without using any response function extrema. Since the volume of the correlation length is proportional to the Riemannian thermodynamic scalar curvature, the locus of the Widom line follows the slope of maximum curvature. We show that the slope of the Widom line follows the slope of the isobaric heat capacity maximum only in the close vicinity of the critical point and that, therefore, the use of response function extrema in this context is problematic. Furthermore, we constructed the vapor-liquid coexistence line for the Lennard-Jones fluid using the fact that the correlation length, and therefore the thermodynamic scalar curvature, must be equal in the two coexisting phases. We compared the resulting phase envelope with those from simulation data where multiple histogram reweighting was used and found striking agreement between the two methods.