Dynamic Crossover in Fluids: From Hard Spheres to Molecules

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.