Transport properties of Mie(14,7) fluids: molecular dynamics simulation and theory

An extensive computer simulation study is presented for the self-diffusion coefficient, the shear viscosity, and the thermal conductivity of Mie(14,7) fluids. The time-correlation function formalism of Green-Kubo is utilized in conjunction with molecular dynamics (MD) simulations. In addition to molecular simulations, the results of a recent study [A. Eskandari Nasrabad, J. Chem. Phys. 128, 154514 (2008)] for the mean free volume are applied to calculate the self-diffusion coefficients within a free volume theory framework. A detailed comparison between the MD simulation and free volume theory results for the diffusion coefficient is given. The density fluctuation theory of shear viscosity is used to compute the shear viscosity and the results are compared to those from MD simulations. The density and temperature dependences of different time-correlation functions and transport coefficients are studied and discussed.     

Scaling law of shear viscosity in atomic liquid and liquid mixtures

A scaling law relating the shear viscosity of one and two component liquid mixtures to their excess thermodynamic entropies defined through pair correlation functions is derived by approximating the mode coupling theory expressions of frictions and then combining with the Stokes-Einstein relation. Molecular dynamics simulation has been performed to generate the data of shear viscosity for one and two component liquid mixtures to test the derived scaling law. The derived scaling laws yield numerical results of shear viscosity for one component and two component liquid mixtures, which are in excellent agreement with the molecular dynamics simulation results for a wide range of density and interaction potential.     

Modified free volume theory of self-diffusion and molecular theory of shear viscosity of liquid carbon dioxide

In previous work on the density fluctuation theory of transport coefficients of liquids, it was necessary to use empirical self-diffusion coefficients to calculate the transport coefficients (e.g., shear viscosity of carbon dioxide). In this work, the necessity of empirical input of the self-diffusion coefficients in the calculation of shear viscosity is removed, and the theory is thus made a self-contained molecular theory of transport coefficients of liquids, albeit it contains an empirical parameter in the subcritical regime. The required self-diffusion coefficients of liquid carbon dioxide are calculated by using the modified free volume theory for which the generic van der Waals equation of state and Monte Carlo simulations are combined to accurately compute the mean free volume by means of statistical mechanics. They have been computed as a function of density along four different isotherms and isobars. A Lennard-Jones site-site interaction potential was used to model the molecular carbon dioxide interaction. The density and temperature dependence of the theoretical self-diffusion coefficients are shown to be in excellent agreement with experimental data when the minimum critical free volume is identified with the molecular volume. The self-diffusion coefficients thus computed are then used to compute the density and temperature dependence of the shear viscosity of liquid carbon dioxide by employing the density fluctuation theory formula for shear viscosity as reported in an earlier paper (J. Chem. Phys. 2000, 112, 7118). The theoretical shear viscosity is shown to be robust and yields excellent density and temperature dependence for carbon dioxide. The pair correlation function appearing in the theory has been computed by Monte Carlo simulations.     

Temperature dependence of fluid transport in nanopores

Understanding the temperature-dependent nanofluidic transport behavior is critical for developing thermomechanical nanodevices. By using non-equilibrium molecular dynamics simulations, the thermally responsive transport resistance of liquids in model carbon nanotubes is explored as a function of the nanopore size, the transport rate, and the liquid properties. Both the effective shear stress and the nominal viscosity decrease with the increase of temperature, and the temperature effect is coupled with other non-thermal factors. The molecular-level mechanisms are revealed through the study of the radial density profile and hydrogen bonding of confined liquid molecules. The findings are verified qualitatively with an experiment on nanoporous carbon.     

Effect of wall roughness on fluid transport resistance in nanopores

Using non-equilibrium molecular dynamics simulations, we investigate the effect of wall roughness on the transport resistance of water molecules inside modified carbon nanotubes. The effective shear stress, which characterizes the strong interaction between liquid molecules and solid wall, is a quantity that dominates the nanofluidic transport resistance. Both the effective shear stress and nominal viscosity arise with the increase of the amplitude or the decrease of the wavelength of roughness. The effect of roughness is also relatively more prominent in smaller nanotubes. The molecular mechanism is elucidated through the study of the radial density profile, hydrogen bonding, and velocity field of the confined water molecules.     

Generic van der Waals equation of state, modified free volume theory of diffusion, and viscosity of simple liquids

The shear viscosity formula derived by the density fluctuation theory in previous papers is computed for argon, krypton, and methane by using the self-diffusion coefficients derived in the modified free volume theory with the help of the generic van der Waals equation of state. In the temperature regime near or above the critical temperature, the density dependence of the shear viscosity can be accounted for by ab initio calculations with the self-diffusion coefficients provided by the modified free volume theory if the minimum (critical) free volume is set equal to the molecular volume and the volume overlap parameter (alpha) is taken about unity in the expression for the self-diffusion coefficient. In the subcritical temperature regime, if the density fluctuation range parameter is chosen appropriately at a temperature, then the resulting expression for the shear viscosity can well account for its density and temperature dependence over the ranges of density and temperature experimentally studied. In the sense that once the density fluctuation range is fixed at a temperature, the theory can account for the experimental data at other subcritical temperatures on the basis of the intermolecular force only; the theory is predictive even in the subcritical regime of temperature. Theory is successfully tested in comparison with experimental data for self-diffusion coefficients and shear viscosity for argon, krypton, and methane.     

Vapor-liquid phase coexistence and transport properties of two-dimensional oligomers

Grand-canonical transition-matrix Monte Carlo and histogram reweighting techniques are used herein to study the vapor-liquid coexistence properties of two-dimensional (2D) flexible oligomers with varying chain lengths (m = 1-8). The phase diagrams of the various 2D oligomers follow the correspondence state (CS) principle, akin to the behavior observed for bulk oligomers. The 2D critical density is not influenced by the oligomer chain length, which contrasts with the observation for the bulk oligomers. Line tension, calculated using Binder's formalism, in the reduced plot is found to be independent of chain length in contrast to the 3D behavior. The dynamical properties of 2D fluids are evaluated using molecular dynamics simulations, and the velocity and pressure autocorrelation functions are investigated using Green-Kubo (GK) relations to yield the diffusion and viscosity. The viscosity determined from 2D non-equilibrium molecular dynamics simulation is compared with the viscosity estimated from the GK relations. The GK relations prove to be reliable and efficient for the calculation of 2D transport properties. Normal diffusive regions are identified in dense oligomeric fluid systems. The influence of molecular size on the diffusivity and viscosity is found to be diminished at specific CS points for the 2D oligomers considered herein. In contrast, the viscosity and diffusion of the 3D bulk fluid, at a reduced temperature and density, are strongly dependent on the molecular size at the same CS points. Furthermore, the viscosity increases and the diffusion decreases multifold in the 2D system relative to those in the 3D system, at the CS points.     

Structure and transport properties of polymer grafted nanoparticles

We perform molecular dynamics simulations on a bead-spring model of pure polymer grafted nanoparticles (PGNs) and of a blend of PGNs with a polymer melt to investigate the correlation between PGN design parameters (such as particle core concentration, polymer grafting density, and polymer length) and properties, such as microstructure, particle mobility, and viscous response. Constant strain-rate simulations were carried out to calculate viscosities and a constant-stress ensemble was used to calculate yield stresses. The PGN systems are found to have less structural order, lower viscosity, and faster diffusivity with increasing length of the grafted chains for a given core concentration or grafting density. Decreasing grafting density causes depletion effects associated with the chains leading to close contacts between some particle cores. All systems were found to shear thin, with the pure PGN systems shear thinning more than the blend; also, the pure systems exhibited a clear yielding behavior that was absent in the blend. Regarding the mechanism of shear thinning at the high shear rates examined, it was found that the shear-induced decrease of Brownian stresses and increase in chain alignment, both correlate with the reduction of viscosity in the system with the latter being more dominant. A coupling between Brownian stresses and chain alignment was also observed wherein the non-equilibrium particle distribution itself promotes chain alignment in the direction of shear.     

Determination of the distance-dependent viscosity of mixtures in parallel slabs using non-equilibrium molecular dynamics

We generalize a technique for determination of the shear viscosity of mixtures in planar slabs using non-equilibrium computer simulations by applying an external force parallel to the surface generating Poiseuille flow. The distance-dependent viscosity of the mixture, given as a function of the distance from the surface, is determined by analysis of the resulting velocity profiles of all species. We present results for a highly non-ideal water + methanol mixture in the whole concentration range between rutile (TiO(2)) walls. The bulk results are compared to the existing equilibrium molecular dynamics and experimental data while the inhomogeneous viscosity profiles at the interface are interpreted using the structural data and information on hydrogen bonding.     

Transport Coefficients Of Ar—Kr Mixtures by Molecular Dynamics Computer Simulation

Abstract

Equilibrium molecular dynamics computer simulations have been used to determine the transport coefficients of model Ar—Kr mixtures, which are represented by Lennard-Jones pair potentials with Lorentz—Berthelot rules for the cross-species interactions. The component self-diffusion and mutual-diffusion coefficients are calculated from time correlation functions and mean square displacements. Time correlation functions are used to evaluate the shear and bulk viscosity, thermal conductivity and the thermal diffusion coefficient (Soret/Dufour coefficient). In the case of the thermal transport coefficients, the partial enthalpy of the two species is required at each state point to define the heat flux rigorously. We obtain this and the partial volume (and species resolved chemical potential) using particle-exchange (and particle insertion) techniques implemented in separate [NPT] simulations at the same state point.

The viscoelasticity of the fluids is characterised by the relaxation times for bulk and shear stress relaxation. The results are for dense liquids close to the triple point temperature and density. Agreement with experiment and previous simulation is particularly good for the density of the mixtures, the shear modulus, shear viscosity, shear stress relaxation time and thermal conductivity. As for the single component noble gas fluids (simulated and experiment) there is a significant qualitative difference in the temperature and, for mixtures, composition dependence of the bulk viscosity.     

Estimation and analysis of the rheological properties of a perfluoropolyether through molecular dynamics simulation

Equilibrium and non-equilibrium molecular dynamics simulations of a perfluoropolyether C₈F₁₈O₄ are reported using an atomistic interaction potential. The bulk rheological properties of the perfluoropolyether are investigated through molecular dynamics simulations as a function of both temperature and shear rate. The effect of molecular structure on viscosity is explored in detail. The rotational relaxation time is reported as a function of temperature. Structural properties, including the mean-square end-to-end chain length, the mean-square radius of gyration of chains, and the distribution functions of bond lengths, bond angles, and bond torsional angles are collected and analyzed as functions of shear rate. After an initial plateau, both mean-square end-to-end chain length and mean-square radius of gyration decrease monotonically with increasing shear rate. The behaviors of the rheological and structural properties are explained through an analysis of the individual contributions due to bond stretching, bond bending, and bond torsion, as well as both intramolecular and intermolecular non-bonded interactions. A further analysis is possible through a meticulous breakdown of each contribution into a specific type of mode; e.g., the total bond stretching is comprised of CC, CO, and CF bond stretching terms. In this way, one can relate the shear viscosity to the specific chemical structure of C₈F₁₈O₄.     

Comparison between theoretical values and simulation results of viscosity for the dissipative particle dynamics method

In order to investigate the validity of the dissipative particle dynamics method, which is a mesoscopic simulation technique, we have derived an expression for viscosity from the equation of motion of dissipative particles. In the concrete, we have shown the Fokker-Planck equation in phase space, and macroscopic conservation equations such as the equation of continuity and the equation of momentum conservation. The basic equations of the single-particle and pair distribution functions have been derived using the Fokker-Planck equation. The solutions of these distribution functions have approximately been solved by the perturbation method under the assumption of molecular chaos. The expressions of the viscosity due to momentum and dissipative forces have been obtained using the approximate solutions of the distribution functions. Also, we have conducted nonequilibrium dynamics simulations to investigate the influence of the parameters, which have appeared in defining the equation of motion in the dissipative particle dynamics method. The theoretical values of the viscosity due to dissipative forces in the Hoogerbrugge-Koelman theory are in good agreement with the simulation results obtained by the nonequilibrium dynamics method, except in the range of small number densities. There are restriction conditions for taking appropriate values of the number density, number of particles, time interval, shear rate, etc., to obtain physically reasonable results by means of dissipative particle dynamics simulations.     

Transport and Helfand moments in the Lennard-Jones fluid. I. Shear viscosity

The authors propose a new method, the Helfand-moment method, to compute the shear viscosity by equilibrium molecular dynamics in periodic systems. In this method, the shear viscosity is written as an Einstein-type relation in terms of the variance of the so-called Helfand moment. This quantity is modified in order to satisfy systems with periodic boundary conditions usually considered in molecular dynamics. They calculate the shear viscosity in the Lennard-Jones fluid near the triple point thanks to this new technique. They show that the results of the Helfand-moment method are in excellent agreement with the results of the standard Green-Kubo method.     

Friction between brush layers of charged and neutral bottle-brush macromolecules. molecular dynamics simulations

Using molecular dynamics simulations, we study the lubricating properties of neutral and charged bottle-brush coatings as a function of the compression and shear stresses and brush grafting density. Our simulations have shown that in charged bottle-brush systems under shear there is a layer with excess counterions located in the middle between brush-bearing surfaces. The main deformation mode of the charged bottle-brush layers is associated with the backbone deformation, resulting in the backbone deformation ratio, α, and shear viscosity, η, being universal functions of the Weissenberg number. In the case of neutral bottle-brush systems, in addition to the backbone deformation there is also side chain deformation. The coupling between backbone and side chain deformation violates universality in the deformation ratio, α, dependence on the Weissenberg number and results in scaling exponents varying with the compression stress and brush grafting density. The existence of different length scales controlling deformation of neutral bottle brushes manifests itself in the shear viscosity, η, dependence on the shear rate, ?γ. Shear viscosity, η, as a function of the shear rate, ?γ, has two plateaus and two shear thinning regimes. The low shear rate plateau and shear thinning regime correspond to the backbone deformation, while the second plateau and shear thinning regime at moderate shear rates are due to side chain deformation. For both systems the value of the friction coefficient increases with increasing shear rate. The values of the friction coefficient for charged bottle-brush systems are about ten times smaller than corresponding values for neutral systems at the same shear rate.     

Twist viscosities and flow alignment of biaxial nematic liquid crystal phases of a soft ellipsoid-string fluid studied by molecular dynamics simulation

We have calculated the twist viscosity and the alignment angle between the director and the stream lines in shear flow of a liquid crystal model system, which forms biaxial nematic liquid crystals, as functions of the density, from the Green-Kubo relations by equilibrium molecular dynamics simulation and by a nonequilibrium molecular dynamics algorithm, where a torque conjugate to the director angular velocity is applied to rotate the director. The model system consists of a soft ellipsoid-string fluid where the ellipsoids interact according a repulsive version of the Gay-Berne potential. Four different length-to-width-to-breadth ratios have been studied. On compression, this system forms discotic or calamitic uniaxial nematic phases depending on the dimensions of the molecules, and on further compression a biaxial nematic phase is formed. In the uniaxial nematic phase there is one twist viscosity and one alignment angle. In the biaxial nematic phase there are three twist viscosities and three alignment angles corresponding to the rotation around the various directors and the different alignments of the directors relative to the stream lines, respectively. It is found that the smallest twist viscosity arises by rotation around the director formed by the long axes, the second smallest one arises by rotation around the director formed by the normals of the broadsides, and the largest one by rotation around the remaining director. The first twist viscosity is rather independent of the density whereas the last two ones increase strongly with density. One finds that there is one stable director alignment relative to the streamlines, namely where the director formed by the long axes is almost parallel to the stream lines and where the director formed by the normals of the broadsides is almost parallel to the shear plane. The relative magnitudes of the components of the twist viscosities span a fairly wide interval so this model should be useful for parameterisation experimental data.     

From Effective Hamiltonians to Fluctuation and Dissipation

Irreversible quantum dynamics is reviewed by means of projected resolvents and effective Hamiltonians. It is shown that for a particular partition the dynamical equations are formally identical to the generalized Langevin equation of non-equilibrium statistical mechanics. Three model applications underline the fundamental role played by fluctuations to study line profiles in spectroscopy and to discuss the relationship between fluctuation and dissipation. A new expression of the parameter q characterizing Fano profiles is presented.Dedicated to Jean-Paul Malrieu who advocated the analyses of molecular wave functions in terms of atomic observables and of their fluctuations.     

System size and control parameter effects in reverse perturbation nonequilibrium molecular dynamics

The issue of system size effects in the reverse perturbation nonequilibrium molecular dynamics method for determining transport coefficients of fluids is examined for the case of the Lennard-Jones model. It is found that when adequate precautions are observed in obtaining linear temperature or momentum profiles, a 250 atom system is adequate for determining the thermal conductivity and the shear viscosity. Also, a means of determining the uncertainties in the transport coefficients is described. The conclusion is that this method is computationally competitive with other simulation methods for estimating transport coefficients.     

Some aspects of theory and simulations of polymers in steady flows

We present a study of various properties of bead spring chains in steady flows. The Langevin equation of the normal modes of the chain is solved by Fourier transformation. From the resulting power spectrum, the autocorrelation functions of all configuration-dependent quantities can be calculated. In equilibrium, the influence of the bead masses on the short-time dynamics is discussed. The influence of different flow fields (shear, elongational and Kramers potential flow) on the mean-square chain dimension is calculated. A comparison with results obtained from non-equilibrium molecular dynamics and Monte Carlo calculations is made. Finally, the influence of shear flow on the configurational and rheological properties of cyclic polymers and on the excluded volume behavior of linear chains is examined.     

Transient-time correlation function applied to mixed shear and elongational flows

The transient-time correlation function (TTCF) method is used to calculate the nonlinear response of a homogeneous atomic fluid close to equilibrium. The TTCF response of the pressure tensor subjected to a time-independent planar mixed flow of shear and elongation is compared to directly averaged non-equilibrium molecular dynamics (NEMD) simulations. We discuss the consequence of noise in simulations with a small rate of deformation. The generalized viscosity for planar mixed flow is also calculated with TTCF. We find that for small rates of deformation, TTCF is far more efficient than direct averages of NEMD simulations. Therefore, TTCF can be applied to fluids with deformation rates which are much smaller than those commonly used in NEMD simulations. Ultimately, TTCF applied to molecular systems is amenable to direct comparison between NEMD simulations and experiments and so in principle can be used to study the rheology of polymer melts in industrial processes.