Dynamic Crossover in Fluids: From Hard Spheres to Molecules

Entropy Scaling of Viscosity IV-Application to 124 Industrially Important Fluids

In our previous work [Yang X.J. Chem. Eng. Data2021, 66, 1385-1398], a residual entropy scaling (RES) approach was developed to link viscosity to residual entropy using a 4-term power function for 39 refrigerants. In further research [Yang X.Int. J. Thermophys.2022, 43, 183], this RES approach was extended to 124 pure fluids containing fluids from light gases (hydrogen and helium) to dense fluids (e.g., heavy hydrocarbons) and fluids with strong association force (e.g., water). In these previous research studies, the model was developed by manual optimization of the power function. The average absolute relative deviation (AARD) of experimental data from the RES model is approximately 3.36%, which is higher than the 2.74% obtained with the various models in REFPROP 10.0. In the present work, the power function was optimized by iteratively fitting the global (fluid-independent power terms) and local parameters (fluid-specific and group-specific parameters) and screening the experimental data. The resulting equation has only three terms instead of four. Most notably, the AARD of the new RES model is reduced down to 2.76%; this is very close to the various multiparameter models in REFPROP 10.0, while the average relative deviation (ARD) amounts to 0.03%, which is smaller than REFPROP 10.0's 0.7%. A Python package is provided for the use of the developped model.

Predicting viscosities and thermal conductivities from dilute gas to dense liquid: Deriving fundamental transfer lengths for momentum and energy exchange in revised Enskog theory

Chapman-Enskog theory has long provided an accurate description of the transport properties of dilute gas mixtures. At elevated densities, revised Enskog theory (RET) provides a framework for describing the departure of the transport properties from their dilute-gas values. Various methods of adapting RET for the description of real fluids have been proposed in the literature. The methods have in common that they incorporate one or more length scales to describe molecular interactions. With few exceptions, the required length scales have been estimated from experimental transport property data. In this work, we introduce two transfer lengths that describe the residual transport of momentum and energy. We derive a model called the exchange-weighted closest approach (EWCA), which links the transfer lengths to the intermolecular potential. Combining the EWCA model with Mie potentials fitted to experimental equilibrium properties yields accurate predictions for several real fluids, including a binary mixture. At higher temperatures, the theory is accurate at surprisingly high densities, even up to the liquid-solid transition of argon. We demonstrate how the transfer lengths can be computed from experimental data or correlations for the transport properties. The transfer lengths obtained in this manner are in good agreement with those obtained from the EWCA model paired with an accurate ab initio potential for argon. The results suggest that kinetic theory, after further developments, can become a predictive theory also for liquids.

A Picture is Worth a Thousand Timesteps: Excess Entropy Scaling for Rapid Estimation of Diffusion Coefficients in Molecular-Dynamics Simulations of Fluids

In molecular-dynamics simulations of fluids, the Einstein-Helfand (EH) and Green-Kubo (GK) relationships are frequently used to compute a variety of transport coefficients, including diffusion coefficients. These relationships are formally valid in the limit of infinite sampling time: The error in the estimate of a transport coefficient (relative to an infinitely long simulation) asymptotically approaches zero as more dynamics are simulated and recorded. In practice, of course, one can only simulate a finite number of particles for a finite amount of time. In this work, we show that in this pre-asymptotic regime, an approach for estimating diffusion coefficients based upon excess entropy scaling (EES) achieves a significantly lower error than either EH or GK relationships at fixed online sampling time. This approach requires access only to structural information at the level of the radial distribution function (RDF). We further demonstrate that the use of a recently developed RDF mollification scheme significantly reduces the amount of sampling time needed to converge to the long-time value of the diffusion coefficient. We also demonstrate favorable sample-to-sample variances in the diffusion coefficient estimate obtained using EES as compared to those obtained using EH and GK.

Linking excess entropy and acentric factor in spherical fluids

Introduced by Pitzer in 1955, the acentric factor (ω) has been used to evaluate a molecule's deviation from the corresponding state principle. Pitzer devised ω based on a concept called perfect liquid (or centric fluid), a hypothetical species perfectly adhering to this principle. However, its physical significance remains unclear. This work attempts to clarify the centric fluid from an excess entropy perspective. We observe that the excess entropy per particle of centric fluids approximates -kB at their critical points, akin to the communal entropy of an ideal gas in classical cell theory. We devise an excess entropy dissection and apply it to model fluids (square-well, Lennard-Jones, Mie n-6, and the two-body ab initio models) to interpret this similarity. The dissection method identifies both centricity-independent and centricity-dependent entropic features. Regardless of the acentric factor, the attractive interaction contribution to the excess entropy peaks at the density where local density is most enhanced due to the competition between the local attraction and critical fluctuations. However, only in centric fluids does the entropic contribution from the local attractive potential become comparable to that of the hard sphere exclusion, making the centric fluid more structured than acentric ones. These findings elucidate the physical significance of the centric fluid as a system of particles where the repulsive and attractive contributions to the excess entropy become equal at its gas-liquid criticality. We expect these findings to offer a way to find suitable intermolecular potentials and assess the physical adequacy of equations of state.

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.

Influence of repulsion on entropy scaling and density scaling of monatomic fluids

Entropy scaling is applied to the shear viscosity, self-diffusion coefficient, and thermal conductivity of simple monatomic fluids. An extensive molecular dynamics simulation series is performed to obtain these transport properties and the residual entropy of three potential model classes with variable repulsive exponents: n, 6 Mie (n = 9, 12, 15, and 18), Buckingham's exponential-six (α = 12, 14, 18, and 30), and Tang-Toennies (αT = 4.051, 4.275, and 4.600). A wide range of liquid and supercritical gas- and liquid-like states is covered with a total of 1120 state points. Comparisons to equations of state, literature data, and transport property correlations are made. Although the absolute transport property values within a given potential model class may strongly depend on the repulsive exponent, it is found that the repulsive steepness plays a negligible role when entropy scaling is applied. Hence, the plus-scaled transport properties of n, 6 Mie, exponential-six, and Tang-Toennies fluids lie basically on one master curve, which closely corresponds with entropy scaling correlations for the Lennard-Jones fluid. This trend is confirmed by literature data of n, 6 Mie, and exponential-six fluids. Furthermore, entropy scaling holds for state points where the Pearson correlation coefficient R is well below 0.9. The condition R > 0.9 for strongly correlating liquids is thus not necessary for the successful application of entropy scaling, pointing out that isomorph theory may be a part of a more general framework that is behind the success of entropy scaling. Density scaling reveals a strong influence of the repulsive exponent on this particular approach.

Universal scaling of transverse sound speed and its isomorphic property in Yukawa fluids

Molecular dynamical simulations are performed to investigate the scaling of the transverse sound speed in two-dimensional (2D) and 3D Yukawa fluids. From the calculated diagnostics of the radial distribution function, the mean-squared displacement, and the Pearson correlation coefficient, the approximate isomorphic curves for 2D and 3D liquidlike Yukawa systems are obtained. It is found that the structure and dynamics of 2D and 3D liquidlike Yukawa systems exhibit the isomorphic property under the conditions of the same relative coupling parameter Γ/Γ_{m}=const. It is demonstrated that the reduced transverse sound speed, i.e., the ratio of the transverse sound speed to the thermal speed, is an isomorph invariant, which is a quasiuniversal function of Γ/Γ_{m}. The obtained isomorph invariant of the reduced transverse sound speed can be useful to estimate the transverse sound speed, or determine the coupling strength, with applications to dusty (complex) plasma or colloidal systems.

Scaling Properties of Liquid Dynamics Predicted from a Single Configuration: Small Rigid Molecules

Isomorphs are curves in the thermodynamic phase diagram along which structure and dynamics are invariant to a good approximation. There are two main ways to trace out isomorphs, the configurational-adiabat method and the direct-isomorph-check method. Recently a new method based on the scaling properties of forces was introduced and shown to work very well for atomic systems [T. B. Schrøder, Phys. Rev. Lett. 2022, 129, 245501]. A unique feature of this method is that it only requires a single equilibrium configuration for tracing out an isomorph. We here test generalizations of this method to molecular systems and compare to simulations of three simple molecular models: the asymmetric dumbbell model of two Lennard-Jones spheres, the symmetric inverse-power-law dumbbell model, and the Lewis-Wahnström o-terphenyl model. We introduce and test two force-based and one torque-based methods, all of which require just a single configuration for tracing out an isomorph. Overall, the method based on requiring invariant center-of-mass reduced forces works best.

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.

Excess Entropy Scaling in Active-Matter Systems

Active-matter systems feature discrete particles that can convert stored or ambient free energy into motion. To realize the engineering potential of active matter, there is a strong need for predictive and theoretically grounded techniques for describing transport in these systems. In this work, we perform molecular-dynamics (MD) simulations of a model active-matter system, in which we vary the total fraction of active particles (0.01 ≤ ϕ ≤ 0.5) as well as the degree of activity of the active particles. These simulations reveal a fascinating array of transport phenomena, including activity-enhanced diffusion coefficients. By adapting an existing result for binary (inactive) fluids, we demonstrate the existence of an excess entropy scaling relation in an active system. This relationship is well supported by our MD results and establishes a new connection between transport (dynamics) and structure (statics) in active matter, a promising step for predictive and generalizable models of other transport phenomena in such systems.

Freezing density scaling of fluid transport properties: Application to liquified noble gases

A freezing density scaling of transport properties of the Lennard-Jones fluid is rationalized in terms of the Rosenfeld's excess entropy scaling and isomorph theory of Roskilde-simple systems. Then, it is demonstrated that the freezing density scaling operates reasonably well for viscosity and thermal conductivity coefficients of liquid argon, krypton, and xenon. Quasi-universality of the reduced transport coefficients at their minima and at freezing conditions is discussed. The magnitude of the thermal conductivity coefficient at the freezing point is shown to agree remarkably well with the prediction of the vibrational model of thermal transport in dense fluids.

Predicting thermodiffusion in simple binary fluid mixtures

The predictive capabilities of some existing theoretical models to quantify thermodiffusion have been investigated in this work. To do so, the tests have been performed on two model fluids, the hard-sphere and the Lennard-Jones (including spheres and dimers) ones, exploring different mixtures and thermodynamic conditions thanks to extensive molecular simulations. It has been confirmed that the thermal diffusion factor should be expressed as the sum of one term related to the isotope effect and one term related to the "chemical" effects and that a kinetic term is required to quantify thermodiffusion from the gas state to the liquid state. In addition, regarding the isotope effects, it has been obtained that none of the available theoretical models are able to yield a reasonable prediction relatively to the molecular simulations results and that the moment of inertia contribution is one order of magnitude smaller than the mass contribution in the liquid state. Finally, concerning the chemical effects, it has been shown the Shukla and Firoozabadi model, complemented with a kinetic term, is probably the most reasonable option to estimate the chemical contribution to the thermal diffusion factor, even if it fails in capturing the effect of the asymmetry in size and in shape between the species. Overall, this works confirms that there is still a lack of a generic model able to predict accurately thermal diffusion factors, or equivalently Soret coefficient, in simple binary mixtures from the gas state to the liquid state.

Freezing Temperature and Density Scaling of Transport Coefficients

It is demonstrated that the freezing density scaling of transport coefficients in fluids, similar to the freezing temperature scaling, originates from the quasi-universal excess entropy scaling approach proposed by Rosenfeld. The freezing density scaling has a considerably wider applicability domain on the phase diagram of Lennard-Jones and related systems. As an illustration of its predictive power, we show that it reproduces with an excellent accuracy the shear viscosity coefficients of saturated liquid argon, krypton, xenon, and methane.

Gas–liquid crossover in the Lennard-Jones system

It is demonstrated that the crossover between gas- and liquid-like regions on the phase diagram of the Lennard-Jones system occurs at a fixed value of the density divided by its value at the freezing point, ρ/ ρfr ≃ 0.35. This definition is consistent with other definitions proposed recently. As a result, a very simple practical expression for the gas-to-liquid crossover line emerges.