Influence of the repulsive coefficient α and approximate corresponding states in Mie α-6 and exponential α-6 fluids

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.

Beyond Born-Mayer: Improved Models for Short-Range Repulsion in ab Initio Force Fields

Short-range repulsion within intermolecular force fields is conventionally described by either Lennard-Jones (A/r(12)) or Born-Mayer (A exp(-Br)) forms. Despite their widespread use, these simple functional forms are often unable to describe the interaction energy accurately over a broad range of intermolecular distances, thus creating challenges in the development of ab initio force fields and potentially leading to decreased accuracy and transferability. Herein, we derive a novel short-range functional form based on a simple Slater-like model of overlapping atomic densities and an iterated stockholder atom (ISA) partitioning of the molecular electron density. We demonstrate that this Slater-ISA methodology yields a more accurate, transferable, and robust description of the short-range interactions at minimal additional computational cost compared to standard Lennard-Jones or Born-Mayer approaches. Finally, we show how this methodology can be adapted to yield the standard Born-Mayer functional form while still retaining many of the advantages of the Slater-ISA approach.

Calculation of the transport properties of a dilute gas consisting of Lennard-Jones chains

The transport properties in the dilute gas limit have been calculated by the classical-trajectory method for a gas consisting of chain-like molecules. The molecules were modelled as rigid chains consisting of spherical segments that interact through a combination of site-site Lennard-Jones 12-6 potentials. Results are reported for shear viscosity, self-diffusion, and thermal conductivity for chains consisting of 1, 2, 3, 4, 5, 6, 7, 8, 10, 13, and 16 segments in the reduced temperature range of 0.3-50. The results indicate that the transport properties increase with temperature and decrease with chain length. At high temperatures the dependence of the transport properties is governed effectively by the repulsive part of the potential. No simple scaling with chain length has been observed. The higher order correction factors are larger than observed for real molecules so far, reaching asymptotic values of 1.019-1.033 and 1.060-1.072 for viscosity and thermal conductivity, respectively. The dominant contribution comes from the angular momentum coupling. The agreement with molecular dynamics calculations for viscosity is within the estimated accuracy of the two methods for shorter chains. However, for longer chains differences of up to 7% are observed.