Rosenfeld proposed two different scaling approaches to model the transport properties of fluids, separated by 22 years, one valid in the dilute gas, and another in the liquid phase. In this work, we demonstrate that these two limiting cases can be connected through the use of a novel approach to scaling transport properties and a bridging function. This approach, which is empirical and not derived from theory, is used to generate reference correlations for the transport properties of the Lennard-Jones 12-6 fluid of viscosity, thermal conductivity, and self-diffusion. This approach, with a very simple functional form, allows for the reproduction of the most accurate simulation data to within nearly their statistical uncertainty. The correlations are used to confirm that for the Lennard-Jones fluid the appropriately scaled transport properties are nearly monovariate functions of the excess entropy from low-density gases into the supercooled phase and up to extreme temperatures. This study represents the most comprehensive metastudy of the transport properties of the Lennard-Jones fluid to date.

Although the freezing of liquids and melting of crystals are fundamental for many areas of the sciences, even simple properties like the temperature–pressure relation along the melting line cannot be predicted today. Here we present a theory in which properties of the coexisting crystal and liquid phases at a single thermodynamic state point provide the basis for calculating the pressure, density and entropy of fusion as functions of temperature along the melting line, as well as the variation along this line of the reduced crystalline vibrational mean-square displacement (the Lindemann ratio), and the liquid's diffusion constant and viscosity. The framework developed, which applies for the sizable class of systems characterized by hidden scale invariance, is validated by computer simulations of the standard 12-6 Lennard-Jones system.

RUMD is a general purpose, high-performance molecular dynamics (MD) simulation package running on graphical processing units (GPU's). RUMD addresses the challenge of utilizing the many-core nature of modern GPU hardware when simulating small to medium system sizes (roughly from a few thousand up to hundred thousand particles). It has a performance that is comparable to other GPU-MD codes at large system sizes and substantially better at smaller sizes. RUMD is open-source and consists of a library written in C++ and the CUDA extension to C, an easy-to-use Python interface, and a set of tools for set-up and post-simulation data analysis. The paper describes RUMD's main features, optimizations and performance benchmarks.

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/TRef(ρ) in which T is the temperature, ρ is the density, and TRef(ρ) 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.

Molecular dynamics simulations have been carried out along four Lennard-Jones (LJ) fluid isomorphs close to the freezing line, covering a temperature, T, in the range of 0.8-350 and a number density, ρ, in the range of 1.1-3.0 in LJ units. Analysis of the transport coefficients is via the Green-Kubo time correlation function method. The radial distribution function, percolation threshold connectivity distance, self-diffusion coefficient, and shear viscosity are shown to be invariant along an isomorph to a very good approximation when scaled with Rosenfeld's macroscopic units, although there are some small departures for T ≃ 1 and lower temperatures. The thermal conductivity is shown for the first time also to be isomorph invariant. In contrast, the Einstein and moment-based frequencies, and especially the bulk viscosity, η b , show poor isomorphic collapse at low T but not surprisingly tend to an "inverse power" potential limiting value in the high T limit. In the case of the bulk viscosity, the significant departures from invariance arise from oscillations in the pressure autocorrelation function at intermediate times, which scale for inverse power potential systems but not for the LJ case, at least in part, as the pressure and bulk elastic moduli are not isomorph invariant.