Transport coefficients of the Lennard-Jones model fluid. I. Viscosity

Predicting viscosities and thermal conductivities from dilute gas to dense liquid: Deriving fundamental transfer lengths for momentum and energy exchange in revised Enskog theory

Chapman-Enskog theory has long provided an accurate description of the transport properties of dilute gas mixtures. At elevated densities, revised Enskog theory (RET) provides a framework for describing the departure of the transport properties from their dilute-gas values. Various methods of adapting RET for the description of real fluids have been proposed in the literature. The methods have in common that they incorporate one or more length scales to describe molecular interactions. With few exceptions, the required length scales have been estimated from experimental transport property data. In this work, we introduce two transfer lengths that describe the residual transport of momentum and energy. We derive a model called the exchange-weighted closest approach (EWCA), which links the transfer lengths to the intermolecular potential. Combining the EWCA model with Mie potentials fitted to experimental equilibrium properties yields accurate predictions for several real fluids, including a binary mixture. At higher temperatures, the theory is accurate at surprisingly high densities, even up to the liquid-solid transition of argon. We demonstrate how the transfer lengths can be computed from experimental data or correlations for the transport properties. The transfer lengths obtained in this manner are in good agreement with those obtained from the EWCA model paired with an accurate ab initio potential for argon. The results suggest that kinetic theory, after further developments, can become a predictive theory also for liquids.

Models to predict configurational adiabats of Lennard-Jones fluids and their transport coefficients

A comparison is made between three simple approximate formulas for the configurational adiabat (i.e., constant excess entropy, sex) lines in a Lennard-Jones (LJ) fluid, one of which is an analytic formula based on a harmonic approximation, which was derived by Heyes et al. [J. Chem. Phys. 159, 224504 (2023)] (analytic isomorph line, AIL). Another is where the density is normalized by the freezing density at that temperature (freezing isomorph line, FIL). It is found that the AIL formula and the average of the freezing density and the melting density ("FMIL") are configurational adiabats at all densities essentially down to the liquid-vapor binodal. The FIL approximation departs from a configurational adiabat in the vicinity of the liquid-vapor binodal close to the freezing line. The self-diffusion coefficient, D, shear viscosity, ηs, and thermal conductivity, λ, in macroscopic reduced units are essentially constant along the AIL and FMIL at all fluid densities and temperatures, but departures from this trend are found along the FIL at high liquid state densities near the liquid-vapor binodal. This supports growing evidence that for simple model systems with no or few internal degrees of freedom, isodynes are lines of constant excess entropy. It is shown that for the LJ fluid, ηs and D can be predicted accurately by an essentially analytic procedure from the high temperature limiting inverse power fluid values (apart from at very low densities), and this is demonstrated quite well also for the experimental argon viscosity.

Computing Accurate True Self-Diffusion Coefficients and Shear Viscosities Using the OrthoBoXY Approach

In a recent paper [Busch, J.; Paschek, D. J. Phys. Chem. B 2023, 127, 7983-7987], we have shown that for molecular dynamics (MD) simulations using orthorhombic periodic boundary conditions with "magic" box length ratios of Lz/Lx = Lz/Ly = 2.7933596497, the self-diffusion coefficients Dx and Dy in x- and y-directions are independent of the system size. They both represent the true self-diffusion coefficient D0 = (Dx + Dy)/2, while the shear viscosity can be calculated from diffusion coefficients in x-, y-, and z-directions, using η = kBT·8.1711245653/[3πLz(Dx + Dy - 2Dz)]. In this contribution, we test this "OrthoBoXY" approach by its application to a variety of different systems: liquid water, dimethyl ether, methanol, triglyme, water/methanol mixtures, water/triglyme mixtures, and imidazolium-based ionic liquids. The chosen systems range from small-sized molecular liquids to complex mixtures and ionic liquids, while spanning a viscosity range of almost 3 orders of magnitude. We assess the efficiency of the method for computing true self-diffusion and viscosity data and provide simple formulas for estimating the required MD simulation lengths and sizes for delivering reliable data with targeted uncertainty levels. Our analysis of the system size dependence of statistical uncertainties for both the viscosity and the self-diffusion coefficient leads us to the conclusion that it is preferable to extend the simulation length instead of increasing the system size. MD simulations consisting of 768 molecules or ion pairs seem to be perfectly adequate.

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.

Efficient single-run implementation of generalized Einstein relation to compute transport coefficients: A binary-based time sampling

A robust and simple implementation of the generalized Einstein formulation using single equilibrium molecular dynamics simulation is introduced to compute diffusion and shear viscosity. The unique features underlying this framework are as follows: (1) The use of a simple binary-based method to sample time-dependent transport coefficients results in a uniform distribution of data on a logarithmic time scale. Although we sample "on-the-fly," the algorithm is readily applicable for post-processing analysis. Overlapping same-length segments are not sampled as they indicate strong correlations. (2) Transport coefficients are estimated using a power law fitting function, a generalization of the standard linear relation, that accurately describes the long-time plateau. (3) The use of a generalized least squares (GLS) fitting estimator to explicitly consider correlations between fitted data points results in a reliable estimate of the statistical uncertainties in a single run. (4) The covariance matrix for the GLS method is estimated analytically using the Wiener process statistics and computed variances. (5) We provide a Python script to perform the fits and automate the procedure to determine the optimal fitting domain. The framework is applied to two fluids, binary hard sphere and a Lennard-Jones near the triple point, and the validity of the single-run estimates is verified against multiple independent runs. The approach should be applicable to other transport coefficients since the diffusive limit is universal to all of them. Given its rigor and simplicity, this methodology can be readily incorporated into standard molecular dynamics packages using on-the-fly or post-processing analysis.

Reassessing the transport properties of fluids: A symbolic regression approach

The viscosity and thermal conductivity coefficients of the Lennard-Jones fluid are extracted through symbolic regression (SR) techniques from data derived from simulations at the atomic scale. This data-oriented approach provides closed form relations that achieve fine accuracy when compared to well-established theoretical, empirical, or approximate equations, fully transparent, with small complexity and high interpretability. The novelty is further outlined by suggesting analytical expressions for estimating fluid transport properties across the whole phase space, from a dilute gas to a dense liquid, by considering only two macroscopic properties (density and temperature). In such expressions, the underlying physical mechanisms are reflected, while, at the same time, it can be a computationally efficient alternative to costly in time and size first principle and/or molecular dynamics simulations.

Vibrational model for thermal conductivity of Lennard-Jones fluids: Applicability domain and accuracy level

Exact mechanisms of thermal conductivity in liquids are not well understood, despite a rich research history. A vibrational model of energy transfer in dense simple liquids with soft pairwise interactions seems adequate to partially fill this gap. The purpose of the present paper is to define its applicability domain and to demonstrate how well it works within the identified applicability domain in the important case of the Lennard-Jones model system. The existing results from molecular dynamics simulations are used for this purpose. Additionally, we show that a freezing density scaling approach represents a very powerful tool to estimate the thermal conductivity coefficient across essentially the entire gas-liquid region of the phase diagram, including metastable regions. A simple practical expression serving this purpose is proposed.

Pseudo hard-sphere viscosities from equilibrium Molecular Dynamics

Transport coefficients like shear, bulk and longitudinal viscosities are sensitive to the intermolecular interaction potential and finite size effects when are numerically determined. For the hard-sphere (HS) fluid, such transport properties are determined almost exclusively with computer simulations. However, their systematic determination and analysis throughout shear stress correlation functions and the Green-Kubo formalism can not be done due to discontinuous nature of the interaction potential. Here, we use the pseudo hard-sphere (PHS) potential to determine pressure correlation functions as a function of volume fraction in order to compute mentioned viscosities. Simulation results are compared to available event-driven molecular dynamics of the HS fluid and also used to propose empirical corrections for the Chapman-Enskog zero density limit of shear viscosity. Moreover, we show that PHS potential is a reliable representation of the HS fluid and can be used to compute transport coefficients. The molecular simulation results of the present work are valuable for further exploration of HS-type fluids or extend the approach to compute transport properties of hard-colloid suspensions.

Calculation of the Transport and Relaxation Properties of the Ar···HCl van der Waals Complex Using a New Potential Energy Surface: Comparison of the Classical and Full Quantum Mechanical Kinetic Theory Results with Molecular Dynamics Simulations

The intermolecular potential energy surface (PES) of the Ar···HCl complex was calculated at the RCCSD(T)/aug-cc-pvQz-BF level of theory. The obtained potential was expanded in terms of Legendre polynomials and fitted to a mathematical model. The fitting results are highly correlated with the ab initio PES data with SD = 5.9 × 10^-3 cm^-1 and average absolute deviation (AAD) = 4.0 × 10^-6 cm^-1. The interaction second virial coefficients (B₁₂) in the temperature range of 190-480 K were calculated by considering classical and first quantum corrections and compared with the available experimental data. A reasonable agreement with the experimental and calculated B₁₂ was obtained. The PES was also used to obtain the rovibrational energy levels, and the spectroscopic rovibrational constants were obtained. It was found that the D₀ values differ ∼2.25 cm^-1 from the experimental values of the ground rovibrational state. Furthermore, the obtained potential was used to calculate the transport and relaxation properties using full quantum close-coupling (CC) formalism and the classical kinetic theory methods based on the Mason-Monchik approximation (MMA). It was found that the deviation between MMA and CC calculations is increased with increasing the temperature due to the higher influence of the rotational degrees of freedom on the transport properties. Also, the contribution of the inelastic (off-diagonal) transitions for diffusion coefficient is higher than the viscosity. Furthermore, the classical molecular dynamics simulations were performed using LJ(12,6) and Vashishta models, to calculate the interaction diffusion and viscosity coefficients, and compared with the results of the full quantum CC calculations. The obtained results confirm that the Vashishta model is better fitted to the ab initio potentials and is more accurate than LJ(12,6) in calculation of the diffusion coefficients.

Modeling the viscoelastic relaxation dynamics of soft particles via molecular dynamics simulation-informed multi-dimensional transition-state theory

Viscoelastic soft colloidal particles have been widely explored in mechanical, chemical, pharmaceutical and other engineering applications due to their unique combination of viscosity and elasticity. The characteristic viscoelastic relaxation time shows an Arrhenius-type (or super-Arrhenius due to temperature-dependent transition attempts) thermally-activated behavior, but a holistic explanation from the relevant transition-state theory remains elusive. In this paper, the viscoelastic relaxation times of Lennard-Jones soft colloidal particle systems, including a single particle type system and a binary particle mixture based on the Kob-Andersen model, are determined using molecular dynamics (MD) simulations as the benchmark. First, the particle systems show a non-Maxwellian behavior after comparing the MD-predicted viscoelastic relaxation time and dynamic moduli (storage and loss modulus) to the classic Maxwell viscoelastic model and the recent particle local connectivity theory. Surprisingly, neither the Maxwell relaxation time τ_Maxwell (obtained from the static shear viscosity η and the high-frequency shear modulus G_∞) nor the particle local connectivity lifetime τ_LC can capture the super-Arrhenius temperature-dependent behavior in the MD-predicted relaxation time τ_MD. Then, the particle dissociation and association transition kinetics, fractal dimensions of the particle systems, and neighbor particle structure (obtained from the radial distribution functions) are shown to collectively determine the viscoelastic relaxation time. These factors are embedded into a new multi-dimensional transition kinetics model to directly estimate the viscoelastic relaxation time τ_Model, which is found to agree with the MD-predicted τ_MD remarkably well. This work highlights the microscopic origin of viscoelastic relaxation dynamics of soft colloidal particles, and theoretically connects rheological dynamics and transition kinetics in soft matters.

Intrinsic viscuit probability distribution functions for transport coefficients of liquids and solids

A reformulation of the Green–Kubo expressions for the transport coefficients of liquids in terms of a probability distribution function (PDF) of short trajectory contributions, which were named “viscuits,” has been explored in a number of recent publications. The viscuit PDF, P, is asymmetric on the two sides of the distribution. It is shown here using equilibrium 3D and 2D molecular dynamics simulations that the viscuit PDF of a range of simple molecular single component and mixture liquid and solid systems can be expressed in terms of the same intrinsic PDF ( P0), which is derived from P with the viscuit normalized by the standard deviation separately on each side of the distribution. P0 is symmetric between the two sides and can be represented for not very small viscuit values by the same gamma distribution formulated in terms of a single disposable parameter. P0 tends to an exponential in the large viscuit wings. Scattergrams of the viscuits and their associated single trajectory correlation functions are shown to distinguish effectively between liquids, solids, and glassy systems. The so-called viscuit square root method for obtaining the transport coefficients is shown to be a useful probe of small and statistically zero self-diffusion coefficients of molecules in the liquid and solid states, respectively. The results of this work suggest that the transport coefficients have a common underlying physical origin, reflecting at a coarse-grained level the traversal statistics of the system through its high-dimensioned potential energy landscape.

Excess entropy and Stokes-Einstein relation in simple fluids

The Stokes-Einstein (SE) relation between the self-diffusion and shear viscosity coefficients operates in sufficiently dense liquids not too far from the liquid-solid phase transition. By considering four simple model systems with very different pairwise interaction potentials (Lennard-Jones, Coulomb, Debye-Hückel or screened Coulomb, and the hard sphere limit) we identify where exactly on the respective phase diagrams the SE relation holds. It appears that the reduced excess entropy s_{ex} can be used as a suitable indicator of the validity of the SE relation. In all cases considered the onset of SE relation validity occurs at approximately s_{ex}≲-2. In addition, we demonstrate that the line separating gaslike and liquidlike fluid behaviours on the phase diagram is roughly characterized by s_{ex}≃-1.

Comparison of equilibrium techniques for the viscosity calculation from DPD simulations

Dissipative Particle Dynamics (DPD) is a powerful mesoscopic modelling technique that is routinely used to predict complex fluid morphology and structural properties. While its ability to quickly scan the conformational space is well known, it is unclear if DPD can correctly calculate the viscosity of complex fluids. In this work, we estimate the viscosity of several unentangled polymer solutions using both the Einstein and Green-Kubo formulas. For this purpose, an Einstein relation is derived analogous to the revised Green-Kubo formula suggested by Jung and Schmid, J. Chem. Phys., 2016, 144, 204104. We show that the DPD simulations reproduce the dynamical behaviour predicted by the theory irrespectively of the values of the conservative and friction parameters used and estimate a Schmidt number compatible to that of a fluid system. Moreover, we observe that the Einstein method requires shorter trajectories to achieve the same statistical accuracy as the Green-Kubo formula. This work shows that DPD can confidently be used to calculate the viscosity of complex fluids and that the statistical accuracy of short trajectories can be improved by using our revised Einstein formula.

Dynamic Crossover in Fluids: From Hard Spheres to Molecules

We propose a simple and generic definition of a demarcation reconciling structural and dynamic frameworks when combined with the entropy scaling framework. This crossover line between gas- and liquid-like behaviors is defined as the curve for which an individual property, the contribution to viscosity due to molecules' translation, is exactly equal to a collective property, the contribution to viscosity due to molecular interactions. Such a definition is shown to be consistent with the one based on the minima of the kinematic viscosity. For the hard sphere, this is shown to be an exact solution. For Lennard-Jones spheres and dimers and for some simple real fluids, this relation holds very well. This crossover line passes nearby the critical point, and for all studied fluids, it is well captured by the critical excess entropy curve for atomic fluids, emphasizing the link between transport properties and local structure.

Correlation in thermal fluctuations induced by phase-locked hydrodynamic modes

Thermal fluctuations constitute a fundamental equilibrium phenomenon whose spatial and temporal correlations are governed by the relevant scales of molecular collisions. From the continuum point of view, thermal fluctuations in a fluid can be regarded as comprising a multitude of hydrodynamic modes (HMs) with random phases, each one having one degree of freedom. We show that in a two-dimensional fluid channel with the Navier slip boundary condition, in which the HMs are represented by periodic arrays of vortex and antivortex pairs, periodic modulation of the slip boundary condition can selectively suppress noncommensurate HMs while phase lock the remaining eigenmodes. As a result, thermal fluctuations would exhibit mesoscopic-scale spatial correlations, manifest as a spatially varying diffusion constant when evaluated from the fluctuation-dissipation theorem. Good agreement is shown with the molecular dynamics results. Such manifestation of equilibrium collective motion implies that instead of just being an alternative mathematical basis for expressing thermal fluctuations, in mesoscopic systems the HMs may be manipulated to have physical consequences very different from those expected in bulk fluid.

Transport properties of Lennard-Jones fluids: Freezing density scaling along isotherms

It is demonstrated that properly reduced transport coefficients (self-diffusion, shear viscosity, and thermal conductivity) of Lennard-Jones fluids along isotherms exhibit quasi-universal scaling on the density divided by its value at the freezing point. Moreover, this scaling is closely related to the density scaling of transport coefficients of hard-sphere fluids. The Stokes-Einstein relation without the hydrodynamic diameter is valid in the dense fluid regime. The lower density boundary of its validity can serve as a practical demarcation line between gaslike and liquidlike regimes.

Viscuit and the fluctuation theorem investigation of shear viscosity by molecular dynamics simulations: The information and the noise

The shear viscosity, η, of model liquids and solids is investigated within the framework of the viscuit and Fluctuation Theorem (FT) probability distribution function (PDF) theories, following Heyes et al. [J. Chem. Phys. 152, 194504 (2020)] using equilibrium molecular dynamics (MD) simulations on Lennard-Jones and Weeks-Chandler-Andersen model systems. The viscosity can be obtained in equilibrium MD simulation from the first moment of the viscuit PDF, which is shown for finite simulation lengths to give a less noisy plateau region than the Green-Kubo method. Two other formulas for the shear viscosity in terms of the viscuit and PDF analysis are also derived. A separation of the time-dependent average negative and positive viscuits extrapolated from the noise dominated region to zero time provides another route to η. The third method involves the relative number of positive and negative viscuits and their PDF standard deviations on the two sides for an equilibrium system. For the FT and finite shear rates, accurate analytic expressions for the relative number of positive to negative block average shear stresses is derived assuming a shifted Gaussian PDF, which is shown to agree well with non-equilibrium molecular dynamics simulations. A similar treatment of the positive and negative block average contributions to the viscosity is also shown to match the simulation data very well.

Thermodynamic or density scaling of the thermal conductivity of liquids

Thermodynamic or density scaling is applied to thermal conductivity (λ) data from the literature for the model Lennard-Jones (12-6) fluid; the noble gases neon to xenon; nitrogen, ethene, and carbon dioxide as examples of linear molecules; the quasi-spherical molecules methane and carbon tetrachloride; the flexible chain molecules n-hexane and n-octane; the planar toluene and m-xylene; the cyclic methylcyclohexane; the polar R132a and chlorobenzene; and ammonia and methanol as H-bonded fluids. Only data expressed as Rosenfeld reduced properties could be scaled successfully. Two different methods were used to obtain the scaling parameter γ, one based on polynomial fits to the group (TV^γ) and the other based on the Avramov equation. The two methods agree well, except for λ of CCl₄. γ for the thermal conductivity is similar to those for the viscosity and self-diffusion coefficient for the smaller molecules. It is significantly larger for the Lennard-Jones fluid, possibly due to a different dependence on packing fraction, and much larger for polyatomic molecules where heat transfer through internal modes may have an additional effect. Methanol and ammonia, where energy can be transmitted through intermolecular hydrogen bonding, could not be scaled. This work is intended as a practical attempt to examine thermodynamic scaling of the thermal conductivity of real fluids. The divergence of the scaling parameters for different properties is unexpected, suggesting that refinement of theory is required to rationalize this result. For the Lennard-Jones fluid, the Ohtori-Iishi version of the Stokes-Einstein-Sutherland relation applies at high densities in the liquid and supercritical region.

The Lennard-Jones potential: when (not) to use it

The Lennard-Jones 12-6 potential (LJ) is arguably the most widely used pair potential in molecular simulations. In fact, it is so popular that the question is rarely asked whether it is fit for purpose. In this paper, we argue that, whilst the LJ potential was designed for noble gases such as argon, it is often used for systems where it is not expected to be particularly realistic. Under those circumstances, the disadvantages of the LJ potential become relevant: most important among these is that in simulations the LJ potential is always modified such that it has a finite range. More seriously, there is by now a whole family of different potentials that are all called Lennard-Jones 12-6, and that are all different - and that may have very different macroscopic properties. In this paper, we consider alternatives to the LJ 12-6 potential that could be employed under conditions where the LJ potential is only used as a typical short-ranged potential with attraction. We construct a class of potentials that are, in many respects LJ-like but that are by construction finite ranged, vanishing quadratically at the cut-off distance, and that are designed to be computationally cheap. Below, we present this potential and report numerical data for its thermodynamic and transport properties, for the most important cases (cut-off distance rc = 2σ ("LJ-like") and rc = 1.2σ (a typical "colloidal" potential)).

An entropy scaling demarcation of gas- and liquid-like fluid behaviors

In this work, we propose a generic and simple definition of a line separating gas-like and liquid-like fluid behaviors from the standpoint of shear viscosity. This definition is valid even for fluids such as the hard sphere and the inverse power law that exhibit a unique fluid phase. We argue that this line is defined by the location of the minimum of the macroscopically scaled viscosity when plotted as a function of the excess entropy, which differs from the popular Widom lines. For hard sphere, Lennard-Jones, and inverse-power-law fluids, such a line is located at an excess entropy approximately equal to -2/3 times Boltzmann's constant and corresponds to points in the thermodynamic phase diagram for which the kinetic contribution to viscosity is approximately half of the total viscosity. For flexible Lennard-Jones chains, the excess entropy at the minimum is a linear function of the chain length. This definition opens a straightforward route to classify the dynamical behavior of fluids from a single thermodynamic quantity obtainable from high-accuracy thermodynamic models.

Density dependence of the dynamical processes governing the velocity autocorrelation function of a quantum fluid

We present an exponential mode analysis of the dynamical processes determining the time behavior of the Kubo velocity autocorrelation function (KVAF) of fluid para-H_{2}, as obtained by ring polymer molecular dynamics simulations at various fluid densities. The mechanisms contributing to the decay of the KVAF are thoroughly characterized at a slightly supercritical temperature, in a density interval ranging from the critical point to the fluid-solid transition. We show that the quantum nature of the system does not influence the specific phenomena and decay channels through which a loss in velocity correlation takes place, since these are the same as found in classical fluids. Similarly, a dynamical crossover is observed with increasing density, signaling the onset of a transverse-like dynamics like in classical systems. We also investigate the effect of density on the processes contributing to the most relevant property of a quantum fluid, namely, the large values of the total and zero-point kinetic energy arising through the Heisenberg uncertainty principle.

A multi-scale model for fluid transport through a bio-inspired passive valve

Tillandsia landbeckii is a rootless plant thriving in the hyper-arid Atacama Desert of Chile. These plants use unique cellulose-based microscopic structures called trichomes to collect fresh water from coastal fog. The trichomes rely on a passive mechanism to maintain an asymmetrical transport of water: they allow for the fast absorption of liquid water deposited by sporadic fog events while preventing evaporation during extended drought periods. Inspired by the trichome's design, we study fluid transport through a micrometric valve. Combining Grand Canonical Monte Carlo with Non-Equilibrium Molecular Dynamics simulations, we first analyze the adsorption and transport of a fluid through a single nanopore at different chemical potentials. We then scale up the atomic results using a lattice approach, and simulate the transport at the micrometric scale. Results obtained for a model Lennard-Jones fluid and TIP4P/2005 water were compared, allowing us to identify the key physical parameters for achieving a passive hydraulic valve. Our results show that the difference in transport properties of water vapor and liquid water within the cellulose layer is the basis for the ability of the Tillandsia trichome to function as a water valve. Finally, we predict a critical pore dimension above which the cellulose layer can form an efficient valve.

Dynamical Origin of the Total and Zero-Point Kinetic Energy in a Quantum Fluid

By applying an exponential mode analysis to ring polymer molecular dynamics simulations of dense fluid parahydrogen, we find that the dynamical processes establishing the time behavior of the Kubo velocity autocorrelation function display the same nature as those already observed in high-density classical fluids. This result permits us to demonstrate that the exponential mode decomposition is a unique tool to identify which dynamical processes lead to one of the most notable properties of quantum fluids: the large value of the mean kinetic energy per particle and the importance of the zero-temperature quantum effects in determining it.

Modified Entropy Scaling of the Transport Properties of the Lennard-Jones Fluid

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.

The role of chemical order in the temperature and composition dependence of the viscosity of liquid alloys

The influence of the chemical order on the viscosity of liquid alloys is investigated by numerical simulation of molecular dynamics. The temperature and composition dependence is discussed in the case of two contrasting alloys: K-Cs and Li-Bi. These two mixtures have different chemical orders, the first one being random and the second one having strong heterocoordination tendencies. In the case of K-Cs, the behavior of the mixture vs temperature is similar to a pure system and its viscosity varies monotonically with the composition. It is not the case for Li-Bi due to its marked chemical order and the heterocoordination tendency is accompanied by a maximum of the viscosity of the mixture when the composition is changed. For the first time, estimates of the temperature dependence of the viscosity of three representative Li-Bi alloys are given.

Probing the link between residual entropy and viscosity of molecular fluids and model potentials

This work investigates the link between residual entropy and viscosity based on wide-ranging, highly accurate experimental and simulation data. This link was originally postulated by Rosenfeld in 1977 [Rosenfeld Y (1977) Phys Rev A 15:2545-2549], and it is shown that this scaling results in an approximately monovariate relationship between residual entropy and reduced viscosity for a wide range of molecular fluids [argon, methane, [Formula: see text], [Formula: see text], refrigerant R-134a (1,1,1,2-tetrafluoroethane), refrigerant R-125 (pentafluoroethane), methanol, and water] and a range of model potentials (hard sphere, inverse power, Lennard-Jones, and Weeks-Chandler-Andersen). While the proposed "universal" correlation of Rosenfeld is shown to be far from universal, when used with the appropriate density scaling for molecular fluids, the viscosity of nonassociating molecular fluids can be mapped onto the model potentials. This mapping results in a length scale that is proportional to the cube root of experimentally measurable liquid volume values.

Shear Viscosity Computed from the Finite-Size Effects of Self-Diffusivity in Equilibrium Molecular Dynamics

A method is proposed for calculating the shear viscosity of a liquid from finite-size effects of self-diffusion coefficients in Molecular Dynamics simulations. This method uses the difference in the self-diffusivities, computed from at least two system sizes, and an analytic equation to calculate the shear viscosity. To enable the efficient use of this method, a set of guidelines is developed. The most efficient number of system sizes is two and the large system is at least four times the small system. The number of independent simulations for each system size should be assigned in such a way that 50%–70% of the total available computational resources are allocated to the large system. We verified the method for 250 binary and 26 ternary Lennard-Jones systems, pure water, and an ionic liquid ([Bmim][Tf₂N]). The computed shear viscosities are in good agreement with viscosities obtained from equilibrium Molecular Dynamics simulations for all liquid systems far from the critical point. Our results indicate that the proposed method is suitable for multicomponent mixtures and highly viscous liquids. This may enable the systematic screening of the viscosities of ionic liquids and deep eutectic solvents.

Self-diffusion coefficient of the square-well fluid from molecular dynamics simulations within the constant force approach

We present a systematic study of the self-diffusion coefficient for a fluid of particles interacting via the square-well pair potential by means of molecular dynamics simulations in the canonical (N, V, T) ensemble. The discrete nature of the interaction potential is modeled by the constant force approximation, and the self-diffusion coefficient is determined for several fluid densities at supercritical thermodynamic states. The dependence of the self-diffusion coefficient on the potential range λ is analyzed in the range of 1.1 ≤ λ ≤ 1.5. The obtained simulation results are in agreement with the self-diffusion coefficient predicted by the Enskog method. Additionally, we show that the diffusion coefficient is very sensitive to the potential range λ. Our results for the self-diffusion coefficient times density extrapolate well to the values in the zero-density limit obtained from the Chapman-Enskog theory for dilute gases. The constant force approximation used in this work to model the discrete pair potentials has shown to be an excellent scheme to compute the transport properties of square-well fluids using molecular dynamics simulations. Finally, the simulation results presented here are useful for improving theoretical approaches, such as the Enskog method.

Viscosity of Lennard-Jones mixtures: A systematic study and empirical law

A systematic study of the viscosity of the binary Lennard-Jones (LJ) mixtures is carried out by equilibrium molecular dynamics simulations via the Green-Kubo relation. The effects of mass, size, and energy-parameter asymmetries on the viscosity and the self-diffusion coefficients are examined separately, both in equimolar mixtures and by varying the molar fractions. The systems are mapped into an effective one-component model according to the van der Waals one-fluid (vdW1) model. Furthermore, using an empirical law for pure LJ liquids, similar to the one proposed recently for liquid sodium, it is shown that the viscosity of the mixtures studied here are well-predicted by the combination of vdW1 fluid and empirical law. The Stokes-Einstein relation in the mixtures has also been investigated. A possible simple extension of this relation, from pure liquids to mixtures, has been proposed and tested.

Determining Quiescent Colloidal Suspension Viscosities Using the Green-Kubo Relation and Image-Based Stress Measurements

By combining confocal microscopy and stress assessment from local structural anisotropy, we directly measure stresses in 3D quiescent colloidal liquids. Our noninvasive and nonperturbative method allows us to measure forces ≲50  fN with a small and tunable probing volume, enabling us to resolve the stress fluctuations arising from particle thermal motions. We use the Green-Kubo relation to relate these measured stress fluctuations to the bulk Brownian viscosity at different volume fractions, comparing against simulations and conventional rheometry measurements. We find that the Green-Kubo analysis gives excellent agreement with these prior results, suggesting that similar methods could be applied to investigations of local flow properties in many poorly understood far-from-equilibrium systems, including suspensions that are glassy, strongly sheared, or highly confined.

Reference Correlation for the Viscosity of Carbon Dioxide

A comprehensive database of experimental and computed data for the viscosity of carbon dioxide (CO2) was compiled and a new reference correlation was developed. Literature results based on an ab initio potential energy surface were the foundation of the correlation of the viscosity in the limit of zero density in the temperature range from 100 K to 2000 K. Guided symbolic regression was employed to obtain a new functional form that extrapolates correctly to T → 0 K and to 10 000 K. Coordinated measurements at low density made it possible to implement the temperature dependence of the Rainwater-Friend theory in the linear-in-density viscosity term. The residual viscosity could be formulated with a scaling term ργ /T the significance of which was confirmed by symbolic regression. The final viscosity correlation covers temperatures from 100 K to 2000 K for gaseous CO2, and from 220 K to 700 K with pressures along the melting line up to 8000 MPa for compressed and supercritical liquid states. The data representation is more accurate than with the previous correlations, and the covered pressure and temperature range is significantly extended. The critical enhancement of the viscosity of CO2 is included in the new correlation.

Density of states and dynamical crossover in a dense fluid revealed by exponential mode analysis of the velocity autocorrelation function

Extending a preceding study of the velocity autocorrelation function (VAF) in a simulated Lennard-Jones fluid [Phys. Rev. E 92, 042166 (2015)PLEEE81539-375510.1103/PhysRevE.92.042166] to cover higher-density and lower-temperature states, we show that the recently demonstrated multiexponential expansion method allows for a full account and understanding of the basic dynamical processes encompassed by a fundamental quantity as the VAF. In particular, besides obtaining evidence of a persisting long-time tail, we assign specific and unambiguous physical meanings to groups of exponential modes related to the longitudinal and transverse collective dynamics, respectively. We have made this possible by consistently introducing the interpretation of the VAF frequency spectrum as a global density of states in fluids, generalizing a solid-state concept, and by giving to specific spectral components, obtained through the VAF exponential expansion, the corresponding meaning of partial densities of states relative to specific dynamical processes. The clear identification of a high-frequency oscillation of the VAF with the near-top excitation frequency in the dispersion curve of acoustic waves is a neat example of the power of the method. As for the transverse mode contribution, its analysis turns out to be particularly important, because the multiexponential expansion reveals a transition marking the onset of propagating excitations when the density is increased beyond a threshold value. While this finding agrees with the recent literature debating the issue of dynamical crossover boundaries, such as the one identified with the Frenkel line, we can add detailed information on the modes involved in this specific process in the domains of both time and frequency. This will help obtain a still missing full account of transverse dynamics, in both its nonpropagating and propagating aspects which are linked through dynamical transitions depending on both the thermodynamic states and the excitation wave vectors.