Freezing density scaling of transport coefficients in the Weeks-Chandler-Andersen fluid

It is shown that the transport coefficients (self-diffusion, shear viscosity, and thermal conductivity) of the Weeks-Chandler-Andersen (WCA) fluid along isotherms exhibit a freezing density scaling (FDS). The functional form of this FDS is essentially the same or closely related to those in the Lennard-Jones fluid, hard-sphere fluid, and some liquefied noble gases. This proves that this FDS represents a quasi-universal corresponding state principle for simple classical fluids with steep interactions. Some related aspects, such as a Stokes-Einstein relation without a hydrodynamic diameter and gas-to-liquid dynamical crossover, are briefly discussed. Simple fitting formulas for the transport coefficients of the dense WCA fluid are suggested.

System Size Dependence of the Diffusion Coefficients in MD Simulations: A Simple Correction Formula for Pure Dense Fluids

A practical correction formula relating the self-diffusion coefficient of dense liquids from molecular dynamics (MD) simulations with periodic boundary conditions to the self-diffusion coefficient in the thermodynamic limit is discussed. This formula applies to pure dense fluids and has a very simple form D = D0(1 - γN-1/3), where D0 is the self-diffusion coefficient in the thermodynamic limit and N is the number of particles in the simulation. The numerical factor γ is dependent on the geometry of the simulation cell. Remarkably, γ ≃ 1.0 is the most popular cubic geometry. The success of this formula is supported by results from MD simulations, including very recent simulations with a "magic" simulation geometry.

Comparing four hard-sphere approximations for the low-temperature WCA melting line

By combining interface-pinning simulations with numerical integration of the Clausius–Clapeyron equation, we accurately determine the melting-line coexistence pressure and fluid/crystal densities of the Weeks–Chandler–Andersen system, covering four decades of temperature. The data are used for comparing the melting-line predictions of the Boltzmann, Andersen–Weeks–Chandler, Barker–Henderson, and Stillinger hard-sphere approximations. The Andersen–Weeks–Chandler and Barker–Henderson theories give the most accurate predictions, and they both work excellently in the zero-temperature limit for which analytical expressions are derived here.

Extreme case of density scaling: The Weeks-Chandler-Andersen system at low temperatures

This paper studies numerically the Weeks-Chandler-Andersen system, which is shown to obey hidden scale invariance with a density-scaling exponent that varies from below 5 to above 500. This unprecedented variation makes it advantageous to use the fourth-order Runge-Kutta algorithm for tracing out isomorphs. Good isomorph invariance of structure and dynamics is observed over more than three orders of magnitude temperature variation. For all state points studied, the virial potential-energy correlation coefficient and the density-scaling exponent are controlled mainly by the temperature. Based on the assumption of statistically independent pair interactions, a mean-field theory is developed that rationalizes this finding and provides an excellent fit to data at low temperatures.

The EXP pair-potential system. II. Fluid phase isomorphs

This paper continues the investigation of the exponentially repulsive EXP pair-potential system of Paper I [A. K. Bacher et al., J. Chem. Phys. 149, 114501 (2018)] with a focus on isomorphs in the low-temperature gas and liquid phases. As expected from the EXP system's strong virial potential-energy correlations, the reduced-unit structure and dynamics are isomorph invariant to a good approximation. Three methods for generating isomorphs are compared: the small-step method that is exact in the limit of small density changes and two versions of the direct-isomorph-check method that allows for much larger density changes. Results from the latter two approximate methods are compared to those of the small-step method for each of the three isomorphs generated by 230 one percent density changes, covering one decade of density variation. Both approximate methods work well.

Simple liquids' quasiuniversality and the hard-sphere paradigm

This topical review discusses the quasiuniversality of simple liquids' structure and dynamics and two possible justifications of it. The traditional one is based on the van der Waals picture of liquids in which the hard-sphere system reflects the basic physics. An alternative explanation argues that all quasiuniversal liquids to a good approximation conform to the same equation of motion, referring to the exponentially repulsive pair-potential system as the basic reference system. The paper, which is aimed at non-experts, ends by listing a number of open problems in the field.

Equilibrium and nonequilibrium dynamics of soft sphere fluids

We use computer simulations to test the freezing-point scaling relationship between equilibrium transport coefficients (self-diffusivity, viscosity) and thermodynamic parameters for soft sphere fluids. The fluid particles interact via the inverse-power potential (IPP), and the particle softness is changed by modifying the exponent of the distance-dependent potential term. In the case of IPP fluids, density and temperature are not independent variables and can be combined to obtain a coupling parameter to define the thermodynamic state of the system. We find that the rescaled coupling parameter, based on its value at the freezing point, can approximately collapse the diffusivity and viscosity data for IPP fluids over a wide range of particle softness. Even though the collapse is far from perfect, the freezing-point scaling relationship provides a convenient and effective way to compare the structure and dynamics of fluid systems with different particle softness. We further show that an alternate scaling relationship based on two-body excess entropy can provide an almost perfect collapse of the diffusivity and viscosity data below the freezing transition. Next, we perform nonequilibrium molecular dynamics simulations to calculate the shear-dependent viscosity and to identify the distinct role of particle softness in underlying structural changes associated with rheological properties. Qualitatively, we find a similar shear-thinning behavior for IPP fluids with different particle softness, though softer particles exhibit stronger shear-thinning tendency. By investigating the distance and angle-dependent pair correlation functions in these systems, we find different structural features in the case of IPP fluids with hard-sphere like and softer particle interactions. Interestingly, shear-thinning in hard-sphere like fluids is accompanied by enhanced translational order, whereas softer fluids exhibit loss of order with shear. Our results provide a systematic evaluation of the role of particle softness in equilibrium and nonequilibrium transport properties and their underlying connection with thermodynamic and structural properties.

Density scaling and quasiuniversality of flow-event statistics for athermal plastic flows

Athermal steady-state plastic flows were simulated for the Kob-Andersen binary Lennard-Jones system and its repulsive version in which the sign of the attractive terms is changed to a plus. Properties evaluated include the distributions of energy drops, stress drops, and strain intervals between the flow events. We show that simulations at a single density in conjunction with an equilibrium-liquid simulation at the same density allow one to predict the plastic flow-event statistics at other densities. This is done by applying the recently established "hidden scale invariance" of simple liquids to the glass phase. The resulting scaling of flow-event properties reveals quasiuniversality, i.e., that the probability distributions of energy drops, stress drops, and strain intervals in properly reduced units are virtually independent of the microscopic pair potentials.

NVU perspective on simple liquids' quasiuniversality

The last half-century of research into the structure, dynamics, and thermodynamics of simple liquids has revealed a number of approximate universalities. This paper argues that simple liquids' reduced-coordinate constant-potential-energy hypersurfaces constitute a quasiuniversal family of compact Riemannian manifolds parametrized by a single number; from this follows the quasiuniversalities.

Thermodynamic properties of non-conformal soft-sphere fluids with effective hard-sphere diameters

In this work we study a set of soft-sphere systems characterised by a well-defined variation of their softness. These systems represent an extension of the repulsive Lennard-Jones potential widely used in statistical mechanics of fluids. This type of soft spheres is of interest because they represent quite accurately the effective intermolecular repulsion in fluid substances and also because they exhibit interesting properties. The thermodynamics of the soft-sphere fluids is obtained via an effective hard-sphere diameter approach that leads to a compact and accurate equation of state. The virial coefficients of soft spheres are shown to follow quite simple relationships that are incorporated into the equation of state. The approach followed exhibits the rescaling of the density that produces a unique equation for all systems and temperatures. The scaling is carried through to the level of the structure of the fluids.

Pressure-energy correlations in liquids. V. Isomorphs in generalized Lennard-Jones systems

This series of papers is devoted to identifying and explaining the properties of strongly correlating liquids, i.e., liquids with more than 90% correlation between their virial W and potential energy U fluctuations in the NVT ensemble. Paper IV [N. Gnan et al., J. Chem. Phys. 131, 234504 (2009)] showed that strongly correlating liquids have "isomorphs," which are curves in the phase diagram along which structure, dynamics, and some thermodynamic properties are invariant in reduced units. In the present paper, using the fact that reduced-unit radial distribution functions are isomorph invariant, we derive an expression for the shapes of isomorphs in the WU phase diagram of generalized Lennard-Jones systems of one or more types of particles. The isomorph shape depends only on the Lennard-Jones exponents; thus all isomorphs of standard Lennard-Jones systems (with exponents 12 and 6) can be scaled onto a single curve. Two applications are given. One tests the prediction that the solid-liquid coexistence curve follows an isomorph by comparing to recent simulations by Ahmed and Sadus [J. Chem. Phys. 131, 174504 (2009)]. Excellent agreement is found on the liquid side of the coexistence curve, whereas the agreement is less convincing on the solid side. A second application is the derivation of an approximate equation of state for generalized Lennard-Jones systems by combining the isomorph theory with the Rosenfeld-Tarazona expression for the temperature dependence of the potential energy on isochores. It is shown that the new equation of state agrees well with simulations.

Pressure-energy correlations in liquids. IV. "Isomorphs" in liquid phase diagrams

This paper is the fourth in a series devoted to identifying and explaining the properties of strongly correlating liquids, i.e., liquids where virial and potential energy correlate better than 90% in their thermal equilibrium fluctuations in the NVT ensemble. For such liquids we here introduce the concept of "isomorphic" curves in the phase diagram. A number of thermodynamic, static, and dynamic isomorph invariants are identified. These include the excess entropy, the isochoric specific heat, reduced-unit static and dynamic correlation functions, as well as reduced-unit transport coefficients. The dynamic invariants apply for both Newtonian and Brownian dynamics. It is shown that after a jump between isomorphic state points the system is instantaneously in thermal equilibrium; consequences of this for generic aging experiments are discussed. Selected isomorph predictions are validated by computer simulations of the Kob-Andersen binary Lennard-Jones mixture, which is a strongly correlating liquid. The final section of the paper relates the isomorph concept to phenomenological melting rules, Rosenfeld's excess entropy scaling, Young and Andersen's approximate scaling principle, and the two-order parameter maps of Debenedetti and co-workers. This section also shows how the existence of isomorphs implies an "isomorph filter" for theories for the non-Arrhenius temperature dependence of viscous liquids' relaxation time, and it explains isochronal superposition for strongly correlating viscous liquids.