Communication: Simple liquids' high-density viscosity

Invariant dynamics in a united-atom model of an ionic liquid

We study a united-atom model of the ionic liquid 1-butyl-1-methylpyrrolidinium bis(trifluoromethyl)sulfonylamide to determine to what extent there exist curves in the phase diagram along which the microscopic dynamics are invariant when expressed in dimensionless, or reduced, form. The initial identification of these curves, termed isodynes, is made by noting that contours of reduced shear viscosity and reduced self-diffusion coefficient coincide to a good approximation. Choosing specifically the contours of reduced viscosity as nominal isodynes, further simulations were carried out for state points on these, and other aspects of dynamics were investigated to study their degree of invariance. These include the mean-squared displacement, shear-stress autocorrelation function, and various rotational correlation functions. These were invariant to a good approximation, with the main exception being rotations of the anion about its long axis. The dynamical features that are invariant have in common that they are aspects that would be relevant for a coarse-grained description of the system; specifically, removing the most microscopic degrees of freedom in principle leads to a simplification of the potential energy landscape, which allows for the existence of isodynes.

Diffusion mobility increases linearly on liquid binodals above triple point

Self-diffusion in fluids has been thoroughly studied numerically, but even for simple liquids just a few scaling relationships are known. Relations between diffusion, excitation spectra, and character of the interparticle interactions remain poorly understood. Here, we show that diffusion mobility of particles in simple fluids increases linearly on the liquid branch of the liquid-gas binodal, from the triple point almost up to the critical point. With molecular dynamics simulations, we considered bulk systems of particles interacting via a generalised Lennard-Jones potential, as well as ethane. Using a two-oscillator model for the analysis of excitations, we observed that the mobility (inverse diffusion) coefficient on the liquid-gas binodal increases linearly above the triple point until the dispersion of high-frequency spectra has a solid-like (oscillating) shape. In terms of a separate mode analysis (of longitudinal and transverse modes), this corresponds to crossed modes in the intermediate range of wavenumbers q, between the hydrodynamic regime (small q) and the regime of individual particle motion (large q). The results should be interesting for a broad community in physics and chemistry of fluids, since self-diffusion is among the most fundamental transport phenomena, important for prospective chemical technologies, micro-, nanofluidics, and biotechnologies.

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.

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.

Generalized hydrodynamics of the Lennard-Jones liquid in view of hidden scale invariance

In recent years lines along which structure and dynamics are invariant to a good approximation, so-called isomorphs, have been identified in the thermodynamic phase diagrams of several model liquids and solids. This paper reports computer simulation data of the transverse and longitudinal collective dynamics at different length scales along an isomorph of the Lennard-Jones system. Our findings are compared to corresponding results along an isotherm and an isochore. Confirming the theoretical prediction, the reduced-unit dynamics of the transverse momentum density is invariant to a good approximation along the isomorph on all time and length scales. Likewise, the wave-vector dependent shear-stress autocorrelation function is found to be isomorph invariant (with minor deviations at very short times). A similar invariance is not seen along the isotherm or the isochore. Using a spatially nonlocal hydrodynamic model for the transverse momentum-density time-autocorrelation function, the macroscopic shear viscosity and its wave dependence are determined, demonstrating that the shear viscosity is isomorphic invariant on all length scales studied. This analysis implies the existence of a length scale that is isomorph invariant in reduced units, i.e., which characterizes each isomorph. The transverse sound-wave velocity, the Maxwell relaxation time, and the rigidity shear modulus are also isomorph invariant. In contrast to the isomorph invariance of all aspects of the transverse dynamics, the reduced-unit dynamics of the mass density is not invariant on length scales longer than the interparticle distance. By fitting to a generalized hydrodynamic model, we extract values for the wave-vector-dependent thermal diffusion coefficient, sound attenuation coefficient, and adiabatic sound velocity. The isomorph variation of these quantities in reduced units on long length scales can be eliminated by scaling with the density-scaling exponent, a fundamental quantity in the isomorph theory framework; this is an empirical observation that remains to be explained theoretically.

Transport properties of Lennard-Jones fluids: Freezing density scaling along isotherms

It is demonstrated that properly reduced transport coefficients (self-diffusion, shear viscosity, and thermal conductivity) of Lennard-Jones fluids along isotherms exhibit quasi-universal scaling on the density divided by its value at the freezing point. Moreover, this scaling is closely related to the density scaling of transport coefficients of hard-sphere fluids. The Stokes-Einstein relation without the hydrodynamic diameter is valid in the dense fluid regime. The lower density boundary of its validity can serve as a practical demarcation line between gaslike and liquidlike regimes.

Isomorph Invariance of Higher-Order Structural Measures in Four Lennard-Jones Systems

In the condensed liquid phase, both single- and multicomponent Lennard-Jones (LJ) systems obey the "hidden-scale-invariance" symmetry to a good approximation. Defining an isomorph as a line of constant excess entropy in the thermodynamic phase diagram, the consequent approximate isomorph invariance of structure and dynamics in appropriate units is well documented. However, although all measures of the structure are predicted to be isomorph invariant, with few exceptions only the radial distribution function (RDF) has been investigated. This paper studies the variation along isomorphs of the nearest-neighbor geometry quantified by the occurrence of Voronoi structures, Frank-Kasper bonds, icosahedral local order, and bond-orientational order. Data are presented for the standard LJ system and for three binary LJ mixtures (Kob-Andersen, Wahnström, NiY2). We find that, while the nearest-neighbor geometry generally varies significantly throughout the phase diagram, good invariance is observed along the isomorphs. We conclude that higher-order structural correlations are no less isomorph invariant than is the RDF.

Vibrational model of thermal conduction for fluids with soft interactions

A vibrational model of heat transfer in simple liquids with soft pairwise interatomic interactions is discussed. A general expression is derived, which involves an averaging over the liquid collective mode excitation spectrum. The model is applied to quantify heat transfer in a dense Lennard-Jones liquid and a strongly coupled one-component plasma. Remarkable agreement with the available numerical results is documented. A similar picture does not apply to the momentum transfer and shear viscosity of liquids.

Excess-entropy scaling in supercooled binary mixtures

Transport coefficients, such as viscosity or diffusion coefficient, show significant dependence on density or temperature near the glass transition. Although several theories have been proposed for explaining this dynamical slowdown, the origin remains to date elusive. We apply here an excess-entropy scaling strategy using molecular dynamics computer simulations and find a quasiuniversal, almost composition-independent, relation for binary mixtures, extending eight orders of magnitude in viscosity or diffusion coefficient. Metallic alloys are also well captured by this relation. The excess-entropy scaling predicts a quasiuniversal breakdown of the Stokes-Einstein relation between viscosity and diffusion coefficient in the supercooled regime. Additionally, we find evidence that quasiuniversality extends beyond binary mixtures, and that the origin is difficult to explain using existing arguments for single-component quasiuniversality.

The role of chemical order in the temperature and composition dependence of the viscosity of liquid alloys

The influence of the chemical order on the viscosity of liquid alloys is investigated by numerical simulation of molecular dynamics. The temperature and composition dependence is discussed in the case of two contrasting alloys: K-Cs and Li-Bi. These two mixtures have different chemical orders, the first one being random and the second one having strong heterocoordination tendencies. In the case of K-Cs, the behavior of the mixture vs temperature is similar to a pure system and its viscosity varies monotonically with the composition. It is not the case for Li-Bi due to its marked chemical order and the heterocoordination tendency is accompanied by a maximum of the viscosity of the mixture when the composition is changed. For the first time, estimates of the temperature dependence of the viscosity of three representative Li-Bi alloys are given.

Onset of transverse (shear) waves in strongly-coupled Yukawa fluids

A simple practical approach to describe transverse (shear) waves in strongly-coupled Yukawa fluids is presented. Theoretical dispersion curves, based on hydrodynamic consideration, are shown to compare favorably with existing numerical results for plasma-related systems in the long-wavelength regime. The existence of a minimum wave number below which shear waves cannot propagate and its magnitude are properly accounted in the approach. The relevance of the approach beyond plasma-related Yukawa fluids is demonstrated by using experimental data on transverse excitations in liquid metals Fe, Cu, and Zn, obtained from inelastic x-ray scattering. Some potentially important relations, scalings, and quasi-universalities are discussed. The results should be interesting for a broad community in chemical physics, materials physics, physics of fluids and glassy state, complex (dusty) plasmas, and soft matter.

Probing the link between residual entropy and viscosity of molecular fluids and model potentials

This work investigates the link between residual entropy and viscosity based on wide-ranging, highly accurate experimental and simulation data. This link was originally postulated by Rosenfeld in 1977 [Rosenfeld Y (1977) Phys Rev A 15:2545-2549], and it is shown that this scaling results in an approximately monovariate relationship between residual entropy and reduced viscosity for a wide range of molecular fluids [argon, methane, [Formula: see text], [Formula: see text], refrigerant R-134a (1,1,1,2-tetrafluoroethane), refrigerant R-125 (pentafluoroethane), methanol, and water] and a range of model potentials (hard sphere, inverse power, Lennard-Jones, and Weeks-Chandler-Andersen). While the proposed "universal" correlation of Rosenfeld is shown to be far from universal, when used with the appropriate density scaling for molecular fluids, the viscosity of nonassociating molecular fluids can be mapped onto the model potentials. This mapping results in a length scale that is proportional to the cube root of experimentally measurable liquid volume values.

Revisiting the Stokes-Einstein relation without a hydrodynamic diameter

We present diffusion coefficient and shear viscosity data for the Lennard-Jones fluid along nine isochores above the critical density, each involving a temperature variation of roughly two orders of magnitude. The data are analyzed with respect to the Stokes-Einstein (SE) relation, which breaks down gradually at high temperatures. This is rationalized in terms of the fact that the reduced diffusion coefficient D ̃ and the reduced viscosity η ̃ are both constant along the system's lines of constant excess entropy (the isomorphs). As a consequence, D ̃ η ̃ is a function of T/T _Ref(ρ) in which T is the temperature, ρ is the density, and T _Ref(ρ) is the temperature as a function of the density along a reference isomorph. This allows one to successfully predict the viscosity from the diffusion coefficient in the studied region of the thermodynamic phase diagram.

Perspective: Excess-entropy scaling

This article gives an overview of excess-entropy scaling, the 1977 discovery by Rosenfeld that entropy determines properties of liquids like viscosity, diffusion constant, and heat conductivity. We give examples from computer simulations confirming this intriguing connection between dynamics and thermodynamics, counterexamples, and experimental validations. Recent uses in application-related contexts are reviewed, and theories proposed for the origin of excess-entropy scaling are briefly summarized. It is shown that if two thermodynamic state points of a liquid have the same microscopic dynamics, they must have the same excess entropy. In this case, the potential-energy function exhibits a symmetry termed hidden scale invariance, stating that the ordering of the potential energies of configurations is maintained if these are scaled uniformly to a different density. This property leads to the isomorph theory, which provides a general framework for excess-entropy scaling and illuminates, in particular, why this does not apply rigorously and universally. It remains an open question whether all aspects of excess-entropy scaling and related regularities reflect hidden scale invariance in one form or other.

Viscosity of Lennard-Jones mixtures: A systematic study and empirical law

A systematic study of the viscosity of the binary Lennard-Jones (LJ) mixtures is carried out by equilibrium molecular dynamics simulations via the Green-Kubo relation. The effects of mass, size, and energy-parameter asymmetries on the viscosity and the self-diffusion coefficients are examined separately, both in equimolar mixtures and by varying the molar fractions. The systems are mapped into an effective one-component model according to the van der Waals one-fluid (vdW1) model. Furthermore, using an empirical law for pure LJ liquids, similar to the one proposed recently for liquid sodium, it is shown that the viscosity of the mixtures studied here are well-predicted by the combination of vdW1 fluid and empirical law. The Stokes-Einstein relation in the mixtures has also been investigated. A possible simple extension of this relation, from pure liquids to mixtures, has been proposed and tested.