Accurate thermodynamics of simple fluids and chain fluids based on first-order perturbation theory and second virial coefficients: uv-theory

Simulation and Data-Driven Modeling of the Transport Properties of the Mie Fluid

This work reports the computation and modeling of the self-diffusivity (D*), shear viscosity (η*), and thermal conductivity (κ*) of the Mie fluid. The transport properties were computed using equilibrium molecular dynamics simulations for the Mie fluid with repulsive exponents (λr) ranging from 7 to 34 and at a fixed attractive exponent (λa) of 6 over the whole fluid density (ρ*) range and over a wide temperature (T*) range. The computed database consists of 17,212, 14,288, and 13,099 data points for self-diffusivity, shear viscosity, and thermal conductivity, respectively. The database is successfully validated against published simulation data. The above-mentioned transport properties are correlated using artificial neural networks (ANNs). Two modeling approaches were tested: a semiempirical formulation based on entropy scaling and an empirical formulation based on density and temperature as input variables. For the former, it was found that a unique formulation based on entropy scaling does not yield satisfactory results over the entire density range due to a divergent and incorrect scaling of the transport properties at low densities. For the latter empirical modeling approach, it was found that regularizing the data, e.g., modeling ρ*D* instead of D*, ln η* instead of η*, and ln κ* instead of κ*, as well as using the inverse of the temperature as an input feature, helps to ease the interpolation efforts of the artificial neural networks. The trained ANNs can model seen and unseen data over a wide range of density and temperature. Ultimately, the ANNs can be used alongside equations of state to regress effective force field parameters from volumetric and transport data.

Local equilibrium in liquid phase shock waves

We have assessed the assumption of local thermodynamic equilibrium in a shock wave by comparing local thermodynamic data generated with nonequilibrium molecular dynamics (NEMD) simulations with results from corresponding equilibrium simulations. The shock had a Mach number approximately equal to 2 in a Lennard-Jones spline liquid. We found that the local equilibrium assumption holds perfectly well behind the wave front, and is a very good approximation in the front itself. This was supported by calculations of the excess entropy production in the shock front with four different methods that use the local equilibrium assumption in different ways. Two of the methods assume local equilibrium between excess thermodynamic variables by treating the shock as an interface in Gibbs's sense. The other two methods are based on the local equilibrium assumption in a continuous description of the shock front. We show for the shock studied in this work that all four methods give excess entropy productions that are in excellent agreement, with an average variance of 3.5% for the nonequilibrium molecular dynamics (NEMD) simulations. In addition, we solved the Navier-Stokes (N-S) equations numerically for the same shock wave using an equilibrium equation of state (EoS) based on a recently developed perturbation theory. The results for the density, pressure, and temperature profiles agree well with the profiles from the NEMD simulations. For instance, the shock waves generated in the two simulations travel with almost the same speed; the average absolute Mach-number deviation of the N-S simulations relative to NEMD is 2.6% in the investigated time interval.

Equation of state for confined fluids

Fluids confined in small volumes behave differently than fluids in bulk systems. For bulk systems, a compact summary of the system’s thermodynamic properties is provided by equations of state. However, there is currently a lack of successful methods to predict the thermodynamic properties of confined fluids by use of equations of state, since their thermodynamic state depends on additional parameters introduced by the enclosing surface. In this work, we present a consistent thermodynamic framework that represents an equation of state for pure, confined fluids. The total system is decomposed into a bulk phase in equilibrium with a surface phase. The equation of state is based on an existing, accurate description of the bulk fluid and uses Gibbs’ framework for surface excess properties to consistently incorporate contributions from the surface. We apply the equation of state to a Lennard-Jones spline fluid confined by a spherical surface with a Weeks–Chandler–Andersen wall-potential. The pressure and internal energy predicted from the equation of state are in good agreement with the properties obtained directly from molecular dynamics simulations. We find that when the location of the dividing surface is chosen appropriately, the properties of highly curved surfaces can be predicted from those of a planar surface. The choice of the dividing surface affects the magnitude of the surface excess properties and its curvature dependence, but the properties of the total system remain unchanged. The framework can predict the properties of confined systems with a wide range of geometries, sizes, interparticle interactions, and wall–particle interactions, and it is independent of ensemble. A targeted area of use is the prediction of thermodynamic properties in porous media, for which a possible application of the framework is elaborated.

Perturbation theories for fluids with short-ranged attractive forces: A case study of the Lennard-Jones spline fluid

It is generally not straightforward to apply molecular-thermodynamic theories to fluids with short-ranged attractive forces between their constituent molecules (or particles). This especially applies to perturbation theories, which, for short-ranged attractive fluids, typically must be extended to high order or may not converge at all. Here, we show that a recent first-order perturbation theory, the uv-theory, holds promise for describing such fluids. As a case study, we apply the uv-theory to a fluid with pair interactions defined by the Lennard-Jones spline potential, which is a short-ranged version of the LJ potential that is known to provide a challenge for equation-of-state development. The results of the uv-theory are compared to those of third-order Barker–Henderson and fourth-order Weeks–Chandler–Andersen perturbation theories, which are implemented using Monte Carlo simulation results for the respective perturbation terms. Theoretical predictions are compared to an extensive dataset of molecular simulation results from this (and previous) work, including vapor–liquid equilibria, first- and second-order derivative properties, the critical region, and metastable states. The uv-theory proves superior for all properties examined. An especially accurate description of metastable vapor and liquid states is obtained, which might prove valuable for future applications of the equation-of-state model to inhomogeneous phases or nucleation processes. Although the uv-theory is analytic, it accurately describes molecular simulation results for both the critical point and the binodal up to at least 99% of the critical temperature. This suggests that the difficulties typically encountered in describing the vapor–liquid critical region are only to a small extent caused by non-analyticity.