Molecular Dynamics Study of the Lennard−Jones Fluid Viscosity: Application to Real Fluids

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.

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.

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.

Predicting Polymer Brush Behavior in Solvents Using the Steepest-Entropy-Ascent Quantum Thermodynamic Framework

The steepest-entropy-ascent quantum thermodynamic (SEAQT) framework is utilized to study the effects of temperature on polymer brushes. The brushes are represented by a discrete energy spectrum, and energy degeneracies obtained through the replica-exchange Wang-Landau algorithm. The SEAQT equation of motion is applied to the density of states to establish a unique kinetic path from an initial thermodynamic state to a stable equilibrium state. The kinetic path describes the brush's evolution in state space, as it interacts with a thermal reservoir. The predicted occupation probabilities along the kinetic path are used to determine the expected thermodynamic and structural properties. The polymer density profile of a polystyrene brush in cyclohexane solvent is predicted using the equation of motion, and it agrees qualitatively with the experimental density profiles. The Flory-Huggins parameter chosen to describe brush-solvent interactions affects the solvent distribution in the brush but has a minimal impact on the polymer density profile. Three types of nonequilibrium kinetic paths with differing amounts of entropy production are considered: a heating path, a cooling path, and a heating-cooling path. Properties such as tortuosity, radius of gyration, brush density, solvent density, and brush chain conformations are calculated for each path.

The receding contact line cools down during dynamic wetting

When a contact line (CL)-where a liquid-vapor interface meets a substrate-is put into motion, it is well known that the contact angle differs between advancing and receding CLs. Using non-equilibrium molecular dynamics simulations, we reveal another intriguing distinction between advancing and receding CLs: while temperature increases at an advancing CL-as expected from viscous dissipation, we show that temperature can drop at a receding CL. Detailed quantitative analysis based on the macroscopic energy balance around the dynamic CL showed that the internal energy change of the fluid due to the change of the potential field along the pathline out of the solid-liquid interface induced a remarkable temperature drop around the receding CL, in a manner similar to latent heat upon phase changes. This result provides new insights for modeling the dynamic CL, and the framework for heat transport analysis introduced here can be applied to a wide range of nanofluidic systems.

Microscopic Origins of the Viscosity of a Lennard-Jones Liquid

Unlike crystalline solids or ideal gases, transport properties remain difficult to describe from a microscopic point of view in liquids, whose dynamics result from complex energetic and entropic contributions at the atomic scale. Two scenarios are generally proposed: one represents the dynamics in a fluid as a series of energy-barrier crossings, leading to Arrhenius-like laws, while the other assumes that atoms rearrange themselves by collisions, as exemplified by the free volume model. To assess the validity of these two views, we computed, using molecular dynamics simulations, the transport properties of the Lennard-Jones fluid and tested to what extent the Arrhenius equation and the free volume model describe the temperature dependence of the viscosity and of the diffusion coefficient at fixed pressure. Although both models reproduce the simulation results over a wide range of pressure and temperature covering the liquid and supercritical states of the Lennard-Jones fluid, we found that the parameters of the free volume model can be estimated directly from local structural parameters, also obtained in the simulations. This consistency of the results gives more credibility to the free volume description of transport properties in liquids.

Predicting thermodiffusion in simple binary fluid mixtures

The predictive capabilities of some existing theoretical models to quantify thermodiffusion have been investigated in this work. To do so, the tests have been performed on two model fluids, the hard-sphere and the Lennard-Jones (including spheres and dimers) ones, exploring different mixtures and thermodynamic conditions thanks to extensive molecular simulations. It has been confirmed that the thermal diffusion factor should be expressed as the sum of one term related to the isotope effect and one term related to the "chemical" effects and that a kinetic term is required to quantify thermodiffusion from the gas state to the liquid state. In addition, regarding the isotope effects, it has been obtained that none of the available theoretical models are able to yield a reasonable prediction relatively to the molecular simulations results and that the moment of inertia contribution is one order of magnitude smaller than the mass contribution in the liquid state. Finally, concerning the chemical effects, it has been shown the Shukla and Firoozabadi model, complemented with a kinetic term, is probably the most reasonable option to estimate the chemical contribution to the thermal diffusion factor, even if it fails in capturing the effect of the asymmetry in size and in shape between the species. Overall, this works confirms that there is still a lack of a generic model able to predict accurately thermal diffusion factors, or equivalently Soret coefficient, in simple binary mixtures from the gas state to the liquid state.

Enhanced local viscosity around colloidal nanoparticles probed by equilibrium molecular dynamics simulations

Nanofluids-dispersions of nanometer-sized particles in a liquid medium-have been proposed for a wide variety of thermal management applications. It is known that a solid-like nanolayer of liquid of typical thicknesses of 0.5-1 nm surrounding the colloidal nanoparticles can act as a thermal bridge between the nanoparticle and the bulk liquid. Yet, its effect on the nanofluid viscosity has not been elucidated so far. In this article, we compute the local viscosity of the nanolayer using equilibrium molecular dynamics based on the Green-Kubo formula. We first assess the validity of the method to predict the viscosity locally. We apply this methodology to the calculation of the local viscosity in the immediate vicinity of a metallic nanoparticle for a wide range of solid-liquid interaction strength, where a nanolayer of thickness 1 nm is observed as a result of the interaction with the nanoparticle. The viscosity of the nanolayer, which is found to be higher than its corresponding bulk value, is directly dependent on the solid-liquid interaction strength. We discuss the origin of this viscosity enhancement and show that the liquid density increment alone cannot explain the values of the viscosity observed. Rather, we suggest that the solid-like structure of the distribution of the liquid atoms in the vicinity of the nanoparticle contributes to the nanolayer viscosity enhancement. Finally, we observe a failure of the Stokes-Einstein relation between viscosity and diffusion close to the wall, depending on the liquid-solid interaction strength, which we rationalize in terms of the hydrodynamic slip.

Residual entropy model for predicting the viscosities of dense fluid mixtures

In this work, we have investigated the mono-variant relationship between the reduced viscosity and residual entropy in pure fluids and in binary mixtures of hydrocarbons and hydrocarbons with dissolved carbon dioxide. The mixtures considered were octane + dodecane, decane + carbon dioxide, and 1,3-dimethylbenzene (m-xylene) + carbon dioxide. The reduced viscosity was calculated according to the definition of Bell, while the residual entropy was calculated from accurate multi-parameter Helmholtz-energy equations of state and, for mixtures, the multi-fluid Helmholtz energy approximation. The mono-variant dependence of reduced viscosity upon residual molar entropy was observed for the pure fluids investigated, and by incorporating two scaling factors (one for reduced viscosity and the other for residual molar entropy), the data were represented by a single universal curve. To apply this method to mixtures, the scaling factors were determined from a mole-fraction weighted sum of the pure-component values. This simple model was found to work well for the systems investigated. The average absolute relative deviation (AARD) was observed to be between 1% and 2% for pure components and a mixture of similar hydrocarbons. Larger deviations, with AARDs of up to 15%, were observed for the asymmetric mixtures, but this compares favorably with other methods for predicting the viscosity of such systems. We conclude that the residual-entropy concept can be used to estimate the viscosity of mixtures of similar molecules with high reliability and that it offers a useful engineering approximation even for asymmetric mixtures.

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.

Cross-over in the dynamics of polymer confined between two liquids of different viscosity

Using molecular dynamics simulations, we analysed the polymer dynamics of chains of different molecular weights entrapped at the interface between two immiscible liquids. We showed that on increasing the viscosity of one of the two liquids the dynamic behaviour of the chain changes from a Zimm-like dynamics typical of dilute polymer solutions to a Rouse-like dynamics where hydrodynamic interactions are screened. We observed that when the polymer is in contact with a high viscosity liquid, the number of solvent molecules close to the polymer beads is reduced and ascribed the screening effect to this reduced number of polymer-solvent contacts. For the longest chain simulated, we calculated the distribution of loop length and compared the results with the theoretical distribution developed for solid/liquid interfaces. We showed that the polymer tends to form loops (although flat against the interface) and that the theory works reasonably well also for liquid/liquid interfaces.

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.

Microscopic flow around a diffusing particle

We report here on the computation of the microscopic flow induced by the motion of a small tagged particle in a fluid from molecular dynamic simulations. It is found that the hydrodynamical Stokes solution with slip boundary conditions is recovered at only a few diameters away from the tagged particle. However, fluctuations of the diffusing particle itself induce a renormalization of the bath viscosity and, more strikingly, an apparent violation of the non-penetrability of the particles in the laboratory frame. The expected zero normal velocity at contact is satisfied only in the particle frame, or for heavy particles. Further evidence of this generalized boundary condition is given by the evaluation of the flow in a granular gas using data from particle tracking experiments.

Bubble evolution and properties in homogeneous nucleation simulations

We analyze the properties of naturally formed nanobubbles in Lennard-Jones molecular dynamics simulations of liquid-to-vapor nucleation in the boiling and the cavitation regimes. The large computational volumes provide a realistic environment at unchanging average temperature and liquid pressure, which allows us to accurately measure properties of bubbles from their inception as stable, critically sized bubbles, to their continued growth into the constant speed regime. Bubble gas densities are up to 50% lower than the equilibrium vapor densities at the liquid temperature, yet quite close to the gas equilibrium density at the lower gas temperatures measured in the simulations: The latent heat of transformation results in bubble gas temperatures up to 25% below those of the surrounding bulk liquid. In the case of rapid bubble growth-typical for the cavitation regime-compression of the liquid outside the bubble leads to local temperature increases of up to 5%, likely significant enough to alter the surface tension as well as the local viscosity. The liquid-vapor bubble interface is thinner than expected from planar coexistence simulations by up to 50%. Bubbles near the critical size are extremely nonspherical, yet they quickly become spherical as they grow. The Rayleigh-Plesset description of bubble-growth gives good agreement in the cavitation regime.

Local shear viscosity of strongly inhomogeneous dense fluids: from the hard-sphere to the Lennard-Jones fluids

This work aims at providing a tractable approach to model the local shear viscosity of strongly inhomogeneous dense fluids composed of spherical molecules, in which the density variations occur on molecular distance. The proposed scheme, which relies on the local density average model, has been applied to the quasi-hard-sphere, the Week-Chandler-Andersen and the Lennard-Jones fluids. A weight function has been developed to deal with the hard-sphere fluid given the specificities of momentum exchange. To extend the approach to the smoothly repulsive potential, we have taken into account that the non-local contributions to the viscosity due to the interactions of particles separated by a given distance are temperature dependent. Then, using a simple perturbation scheme, the approach is extended to the Lennard-Jones fluids. It is shown that the viscosity profiles of inhomogeneous dense fluids deduced from this approach are consistent with those directly computed by non-equilibrium molecular dynamics simulations.

Thermal fluctuations of hydrodynamic flows in nanochannels

Flows at the nanoscale are subject to thermal fluctuations. In this work, we explore the consequences for a fluid confined within a channel of nanometric size. First, the phenomenon is illustrated on the basis of molecular dynamics simulations. The center of mass of the confined fluid is shown to perform a stochastic, non-Markovian motion, whose diffusion coefficient satisfies Einstein's relation, and which can be further characterized by the fluctuation relation. Next, we develop an analytical description of the thermally induced fluid motion. We compute the time- and space-dependent velocity correlation function, and characterize its dependence on the nanopore shape, size, and boundary slip at the surface. The experimental implications for mass and charge transports are discussed for two situations. For a particle confined within the nanopore, we show that the fluid fluctuating motion results in an enhanced diffusion. The second situation involves a charged nanopore in which fluid motion within the double layer induces a fluctuating electric current. We compute the corresponding contribution to the current power spectrum.

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.

Local viscosity of a fluid confined in a narrow pore

In this paper, molecular dynamics simulations of a simple Lennard-Jones fluid confined in narrow slit pores and undergoing shear have been performed. The aim is to investigate the effects of density inhomogeneities at the fluid-solid interfaces on the shear viscosity profiles. It has been found that the local viscosity was varying strongly with the distance from the solid walls for both dilute and dense fluid states with oscillations correlated to the density ones. To describe the computed viscosity profiles, we propose a scheme that uses the local average density model, combined with an adequate weight function, for the configurational viscosity and a semiempirical model for the translational viscosity. It is shown that the proposed approach is able to provide viscosity profiles in good agreement with those coming from simulations for different pore widths and for different fluid states (dilute to dense).

Nanoscale simple-fluid behavior under steady shear

In this study, we use two nonequilibrium molecular dynamics algorithms, boundary-driven shear and homogeneous shear, to explore the rheology and flow properties of a simple fluid undergoing steady simple shear. The two distinct algorithms are designed to elucidate the influences of nanoscale confinement. The results of rheological material functions, i.e., viscosity and normal pressure differences, show consistent Newtonian behaviors at low shear rates from both systems. The comparison validates that confinements of the order of 10 nm are not strong enough to deviate the simple fluid behaviors from the continuum hydrodynamics. The non-Newtonian phenomena of the simple fluid are further investigated by the homogeneous shear simulations with much higher shear rates. We observe the "string phase" at high shear rates by applying both profile-biased and profile-unbiased thermostats. Contrary to other findings where the string phase is found to be an artifact of the thermostats, we perform a thorough analysis of the fluid microstructures formed due to shear, which shows that it is possible to have a string phase and second shear thinning for dense simple fluids.