Three‐body interactions in fluids from molecular simulation: Vapor–liquid phase coexistence of argon

Influence of molecular parameters on the representativeness of interfacial properties of simple fluids

New parameterizations for the Lennard-Jones 12/6 potential capable of reproducing the vapor pressure and surface tension with sufficient precision, but not the liquid-vapor equilibrium densities for the case of simple fluids that include Ar, Kr, Xe, Ne, and CH4 are presented in this work. These results are compared with those derived from the family of Mie(n, 6) potentials, which adequately reproduce the coexistence curve and the vapor pressure, leaving aside the surface tension. In addition, a detailed analysis is presented on different parameterizations and methodologies, which have been developed in recent decades to estimate the interfacial properties of interest here for simple fluids, such as argon, which is a molecule that is, in principle, "simple" to study but that clearly reveals the enormous discrepancy between the results reported in the literature throughout these years. These facts undoubtedly reveal one of the fundamental problems in the context of molecular thermodynamics of fluids: reproducing different thermodynamic properties with sufficient precision from a single set of free parameters for some interaction potential. In order to show the scope of the parameterizations presented for the Lennard-Jones model, they were successfully applied to the case of binary mixtures, which included Ar-Kr, Ar-CH4, and Xe-Kr. Finally, and with the aim of showing a possible solution to the problem posed in this research, results of the same interfacial properties above mentioned for argon and methane were presented in this work by using a set of molecular interactions, called ANC2s, whose flexibility allowed to reproduce the experimental evidence with just one parameterization. The results reported in this work were generated using molecular dynamics simulations.

Molecular simulation of orthobaric isochoric heat capacities near the critical point

A molecular simulation strategy is investigated for detecting the divergence of the isochoric heat capacity (C_{V}) on the vapor and liquid coexistence branches of a fluid near the critical point. The procedure is applied to the empirical Lennard-Jones potential and accurate state-of-the-art ab initio two-body and two-body + three-body potentials for argon. Simulations with the Lennard-Jones potential predict the divergence of C_{V}, and the phenomenon is also observed for both two-body and two-body + three body potentials. The potentials also correctly predict the crossover between vapor and liquid C_{V} values and the subcritical liquid C_{V} minimum, which marks the commencement C_{V} divergence. The effect of three-body interactions is to delay the onset of divergence to higher subcritical temperatures.

Temperature and composition dependence of the Soret coefficient in Lennard-Jones mixtures presenting evaporation/condensation phase transition

The thermodiffusive behavior of a Lennard-Jones binary mixture has been studied by using nonequilibrium molecular dynamics. In particular, the dependence of the Soret coefficient, S(T), on the temperature and composition has been investigated, exploring a wide range of temperatures from 1000 K to the condensation temperature of the mixture. In a previous paper the dependence of S(T) on the temperature and the composition was studied for Lennard-Jones binary mixtures presenting mixing/demixing (consolute) phase transition, and the results allowed the formulation of a very simple expression with the computed values of S(T) in the one phase region outside the critical region closely fitted by the function [T - T(c)(x(1))](-1), with T(c)(x(1)) the demixing temperature of the mixture under study. The results of the present work show that the same expression of S(T) can be found for the one phase region outside the evaporation/condensation region but now with T(c) representing the condensation temperature of the mixture under study.

Hard-wall potential function for transport properties of alkali metal vapors

This study demonstrates that the transport properties of alkali metals are determined principally by the repulsive wall of the pair interaction potential function. The (hard-wall) Lennard-Jones (LJ) (15-6) effective pair potential function is used to calculate the transport collision integrals. Accordingly, reduced collision integrals of K, Rb, and Cs metal vapors are obtained from the Chapman-Enskog solution of the Boltzmann equation. The law of corresponding states based on the experimental transport reduced collision integral is used to verify the validity of a LJ(15-6) hybrid potential in describing the transport properties. LJ(8.5-4) potential function and a simple thermodynamic argument with the input PVT data of liquid metals provide the required molecular potential parameters. Values of the predicted viscosity of monatomic alkali metal vapor are in agreement with typical experimental data with average absolute deviations of 2.97% for K in the range of 700-1500 K, 1.69% for Rb, and 1.75% for Cs in the range of 700-2000 K. In the same way, the values of predicted thermal conductivity are in agreement with experiment within 2.78%, 3.25%, and 3.63% for K, Rb, and Cs, respectively. The LJ(15-6) hybrid potential with a hard-wall repulsion character conclusively predicts the best transport properties of the three alkali metal vapors.

Analytic dependence of the pressure and energy of an atomic fluid under shear

Nonequilibrium molecular dynamics simulations are reported at different strain rates (gamma;) for a shearing atomic fluid interacting via accurate two- and three-body potentials. We report that the hydrostatic pressure has a strain-rate dependence of gamma;(2), in contrast to the gamma;(3/2) dependence predicted by mode-coupling theory. Our results indicate that the pressure and energy of real fluids may display an analytic dependence on the strain rate. This is in contrast to previous work using either Lennard-Jones or Weeks-Chandler-Anderson potentials that had shown a gamma;(3/2) dependence of pressure and energy.