微孔中简单流体粘度的分子动力学模拟及关联模型

用分子动力学模拟计算了微孔介质中流体氩在不同温度、不同密度和不同孔径下的剪切粘度.并根据Chapman-Enskog关于硬球流体传递性质的理论以及Heyes的关于Lennard-Jones流体粘度的表达式,提出了两个描述微孔介质中流体粘度的模型,该模型可以计算微孔中流体氩在不同状态下的粘度值.通过与计算机模拟值的比较,证明这两个微孔流体粘度模型是可用的.     

Mode-coupling theory for multiple-time correlation functions of tagged particle densities and dynamical filters designed for glassy systems

The theoretical framework for higher-order correlation functions involving multiple times and multiple points in a classical, many-body system developed by Van Zon and Schofield [Phys. Rev. E 2002, 65, 011106] is extended here to include tagged particle densities. Such densities have found an intriguing application as proposed measures of dynamical heterogeneities in structural glasses. The theoretical formalism is based upon projection operator techniques which are used to isolate the slow time evolution of dynamical variables by expanding the slowly evolving component of arbitrary variables in an infinite basis composed of the products of slow variables of the system. The resulting formally exact mode-coupling expressions for multiple-point and multiple-time correlation functions are made tractable by applying the so-called N-ordering method. This theory is used to derive for moderate densities the leading mode coupling expressions for indicators of relaxation type and domain relaxation, which use dynamical filters that lead to multiple-time correlations of a tagged particle density. The mode coupling expressions for higher order correlation functions are also successfully tested against simulations of a hard sphere fluid at relatively low density.     

Transport coefficients of square-well fluids

We analyze the analytical form of the velocity time correlation function of a hard sphere system obtained by employing generalized Langevin equation for a square-well fluid. The self-diffusion coefficient and shear viscosity have been calculated using this analytical form of velocity tcf for a square-well fluid. The addition of an attractive square-well potential in place of hard sphere leads to a substantial influence on transport coefficients. Unlike harmonic model diffusion coefficient no longer vanishes. A breakdown of the Stokes–Einstein relation is observed at low densities for a square-well fluid.     

Structure of polymer fluids in colloidal dispersions

We have calculated some structural properties of a fluid of hard sphere polymer chains about a variable sized central hard sphere with the Monte Carlo method in the canonical ensemble. We have additionally calculated these structural properties with an integral equation based on density functional theory. The integral equation theory gives good agreement with the simulations at all but the highest densities.     

A diagrammatic formulation of the kinetic theory of fluctuations in equilibrium classical fluids. IV. The short time behavior of the memory function

Using a recently developed diagrammatic formulation of the kinetic theory of fluctuations in liquids, we investigate the short time behavior of the memory function for density fluctuations in a classical atomic fluid. At short times, the memory function has a large contribution that is generated by the repulsive part of the interatomic potential. We introduce a small parameter that is a measure of the softness of the repulsive part of the potential. The diagrams in the memory function that contribute to lowest order in that small parameter are identified and summed to give an explicit expression for the dominant contribution to the memory function at short times. The result leads to a theory for fluids with continuous potentials that is similar to the Enskog theory for hard sphere fluids.     

The effect of rotational and translational energy exchange on tracer diffusion in rough hard sphere fluids

A study is presented of tracer diffusion in a rough hard sphere fluid. Unlike smooth hard spheres, collisions between rough hard spheres can exchange rotational and translational energy and momentum. It is expected that as tracer particles become larger, their diffusion constants will tend toward the Stokes-Einstein hydrodynamic result. It has already been shown that in this limit, smooth hard spheres adopt "slip" boundary conditions. The current results show that rough hard spheres adopt boundary conditions proportional to the degree of translational-rotational energy exchange. Spheres for which this exchange is the largest adopt "stick" boundary conditions while those with more intermediate exchange adopt values between the "slip" and "stick" limits. This dependence is found to be almost linear. As well, changes in the diffusion constants as a function of this exchange are examined and it is found that the dependence is stronger than that suggested by the low-density, Boltzmann result. Compared with smooth hard spheres, real molecules undergo inelastic collisions and have attractive wells. Rough hard spheres model the effect of inelasticity and show that even without the presence of attractive forces, the boundary conditions for large particles can deviate from "slip" and approach "stick."     

Isotropic-nematic transition of hard rods immersed in random sphere matrices

Using replica density functional theory and Monte Carlo computer simulations we investigate a system of annealed hard spherocylinders adsorbed in a matrix of quenched hard spheres. Theoretical predictions for the partition coefficient, defined as the ratio of density of rods in the matrix and that in a reservoir, agree well with simulation results. Theory predicts the isotropic-nematic transition to remain first order upon increasing sphere packing fraction, and to shift towards lower rod densities. This scenario is consistent with our simulation results that clearly show a jump in the nematic order parameter upon increasing the rod density at constant matrix packing fraction, corresponding to the isotropic-nematic transition, even for sphere matrix packing fractions < or approximately equal to 0.3.     

Momentum and stress relaxation in fluids illustrating caging

The self diffusion coefficient, shear viscosity, and velocity time correlation function are calculated for a hard sphere fluid under a severe assumption, namely, the friction arises from uncorrelated binary collisions and from correlated backscattering (caging) collisions as represented in the memory function. Relaxation of the memory function from its zerotime caging value is described as a diffusion process. Derived diffusion coefficients and the shear viscosities, relative to their Enskog values decrease and increase with density, respectively, in a monotonic and gradual fashion in contrast with simulation values that show a precipitous change near the fluid-solid transition. In the present pair diffusion model, the velocity time correlation function vanishes at the proper time but its tail is overly damped relative to the simulation data. A weak breakdown of the Stokes-Einstein relation is also predicted.     

Free energy density functional for adsorption of fluids in nanopores

A classical free energy density functional, which is isomorphic to a usual effective hard sphere model + mean field approximation for tail contribution, is proposed for treatment of real fluids in inhomogeneous states. In the framework of the classical density functional theory (DFT), the present functional is applied to two representative model fluids, namely, a Lennard-Jones fluid and a hard core attractive Yukawa fluid, subject to influence of various external fields. A comprehensive comparison with simulation results and a detailed analysis show that the present functional holds simultaneously all of the desirable properties inherent in an excellent functional, such as high accuracy, computational simplicity, consistency with a hard wall sum rule, nonrecourse to use of adjustable parameter(s) and weighted densities, reproduction of bulk second-order direct correlation function (DCF) in bulk limit, and applicability to subcritical fluid phenomena.     

Structure of spherical electric double layers: a density functional approach

A density functional theory is presented for the structure of spherical electric double layers within the restricted primitive model, where the macroion is considered as a hard sphere having uniform surface charge density, the small ions as charged hard spheres, and the solvent is taken as a dielectric continuum. The theory is partially perturbative as the hard-sphere contribution to the one-particle correlation function is evaluated using suitably averaged weighted density and the ionic part is obtained through a second-order functional Taylor expansion around the uniform fluid. The theory is in quantitative agreement with Monte Carlo simulation for the density profiles and the zeta potentials over a wide range of macroion sizes and electrolyte concentrations. The theory is able to provide interesting insights about the layering and the charge inversion phenomena occurring at the interface.     

A fundamental measure theory for the sticky hard sphere fluid

We construct a density functional theory (DFT) for the sticky hard sphere (SHS) fluid which, like Rosenfeld's fundamental measure theory (FMT) for the hard sphere fluid [Y. Rosenfeld, Phys. Rev. Lett. 63, 980 (1989)], is based on a set of weighted densities and an exact result from scaled particle theory (SPT). It is demonstrated that the excess free energy density of the inhomogeneous SHS fluid Φ(SHS) is uniquely defined when (a) it is solely a function of the weighted densities from Kierlik and Rosinberg's version of FMT [E. Kierlik and M. L. Rosinberg, Phys. Rev. A 42, 3382 (1990)], (b) it satisfies the SPT differential equation, and (c) it yields any given direct correlation function (DCF) from the class of generalized Percus-Yevick closures introduced by Gazzillo and Giacometti [J. Chem. Phys. 120, 4742 (2004)]. The resulting DFT is shown to be in very good agreement with simulation data. In particular, this FMT yields the correct contact value of the density profiles with no adjustable parameters. Rather than requiring higher order DCFs, such as perturbative DFTs, our SHS FMT produces them. Interestingly, although equivalent to Kierlik and Rosinberg's FMT in the case of hard spheres, the set of weighted densities used for Rosenfeld's original FMT is insufficient for constructing a DFT which yields the SHS DCF.     

流体熵相关性质的Monte Carlo模拟新方法

提出了用于计算实际体系熵相关性质的Monte Carlo 多级取样分子模拟方法.应用这一方法,对硬球流体的化学势及Helmholtz 自由能进行了估算,得到了满意的结果.计算化学势时,不存在通常试验粒子方法所遇到的高密度问题.该方法特别适合规律性的系统研究,较之普通模拟方法要有效得多.模拟得到的硬球体系无限稀释组份的超额化学势与对比直径的关系,在相变区域为一条双凹曲线;无论是在相变区还是在单相区,Carnahan-Starling 公式对这一关系的描述均有较大偏差.     

流体熵相关性质的Monte Carlo模拟新方法

提出了用于计算实际体系熵相关性质的Monte Carlo多级取样分子模拟方法. 应用这一方法, 对硬球流体的化学势及Helmholtz自由能进行了估算, 得到了满意的结果. 计算化学势时, 不存在通常试验粒子方法所遇到的高密度问题. 该方法特别适合规律性的系统研究, 较之普通模拟方法要有效得多. 模拟得到的硬球体系无限稀释组份的超额化学势与对比直径的关系, 在相变区域为一条双凹曲线; 无论是在相变区还是在单相区, Carnahan-Starling公式对这一关系的描述均有较大偏差.     

Rough hard spheres treatment of tracer diffusion of each component in two-component liquid systems

The rough hard spheres treatment of diffusion is extended to cover tracer diffusion of both species in five two-component organic liquid systems. Self-consistent results are obtained for the factor expressing coupling of translational-rotational motion and reasonable agreement is found where comparisons can be made with mutual diffusion data.     

Theory of transport in liquid metals: III. Calculation of shear viscosity coefficients of binary alloys

A theoretical procedure for the calculation of the shear viscosity coefficients of liquid binary alloys is proposed based on the simulation of an alloy by a hypothetical single-component hard-sphere liquid and treatment of the latter according to a corrected Enskog approach. The Enskog result is corrected according to a corresponding states principles developed by the authors in a previous work. The pair correlation function at contact necessary for the Enskog calculation is obtained from the solution of the Carnahan and Starling equation of state for single component hard-sphere fluids. The proposed theory is tested numerically by comparing its viscosity results for fifteen alloys with experimental data and with the corresponding results of the Eyring approach and the Tham and Gubbins theory. The latter theory is implemented numerically with pair-correlation values at contact satisfying exactly the Carnahan and Starling equation of state for binary hard-sphere fluids. These values result from analytic and closed-form expressions for the pair-correlation function developed presently through a correction of the corresponding Lebowitz expressions.     

一个新的桥泛函及其在非均一流体密度泛函理论中的应用

基于对OZ 方程的渐近行为与Taylor级数展开的分析，提出了一个新的桥泛函，桥泛函被表达为间接相关函数的函数，Taylor级数展开的重整化导致了一个可调参数，通过将所提出的桥泛函与一个最近提出的密度泛函理论方法学，以及单个硬墙的sum 规则结合，可以确定可调参数.所提出的桥泛函能预言如下非均一流体的密度分布：硬球流体接近一个硬墙与在球形空隙内，Lennard-Jones 流体与缔合硬球流体在两个硬墙之内.理论预言与文献所报导的模拟数据符合很好.     

Viscosity and thermal conductivity of some refrigerants at moderate and high densities on the basis of Rainwater–Friend theory and the corresponding states principle

The initial density dependence of viscosity and thermal conductivity was formulated on the basis of Rainwater–Friend (RF) theory. In this work, we have first focused on the calculation of viscosity and thermal conductivity of moderately dense argon by using RF theory and an accurate ab initio potential function. This theory which was originally presented for spherical potentials have been adapted for calculation of viscosity and translational contribution of thermal conductivity of some refrigerants by introducing the corresponding states correlations for the second transport virial coefficients. Then the internal states contribution for the thermal conductivity has been determined based on the Mason–Monchick and modified Enskog theories. So, we have calculated the viscosity and thermal conductivity of some refrigerants, R32, R14, R12, R13, R22, R134a, R143a, R125, R123, R142b, at moderate densities up to about 2 mol dm⁻³. At high densities, beyond the validity range of RF theory, we have applied correlation expressions for the viscosity and thermal conductivity residual functions to calculate the viscosity and thermal conductivity of supercritical refrigerants and then compared with the available experimental data. In conclusion, we have shown that the RF theory in conjunction with the corresponding states residual functions present the reliable model for calculation of viscosity and thermal conductivity of refrigerants over a comparatively wide temperature and pressure range up to 65 MPa within the experimental errors.