Heat capacity of square-well fluids of variable width

We have obtained by Monte Carlo NVT simulations the constant-volume excess heat capacity of square-well fluids for several temperatures, densities and potential widths. Heat capacity is a thermodynamic property much more sensitive to the accuracy of a theory than other thermodynamic quantities, such as the compressibility factor. This is illustrated by comparing the reported simulation data for the heat capacity with the theoretical predictions given by the Barker-Henderson perturbation theory as well as with those given by a non-perturbative theoretical model based on Baxter's solution of the Percus-Yevick integral equation for sticky hard spheres. Both theories give accurate predictions for the equation of state. By contrast, it is found that the Barker-Henderson theory strongly underestimates the excess heat capacity for low to moderate temperatures, whereas a much better agreement between theory and simulation is achieved with the non-perturbative theoretical model, particularly for small well widths, although the accuracy of the latter worsens for high densities and low temperatures, as the well width increases.     

方阱流体的配位数模型

基于格子气体理论，提出了方阱流体及其混合物的新配位数模型和构型能表达式。新模型与Ｌｅｅ等人的ＭｏｎｔｅＣａｒｌｏ计算机模拟结果进行了比较，取得了较好的计算结果。将新模型应用于分子尺寸不同的混合物时，计算结果尤为满意。     

The role of higher-order terms in perturbation approaches to the monomer and bonding contributions in a SAFT-type equation of state for square-well chain fluids

A perturbation theory for square-well chain fluids is developed within the scheme of the (generalised) Wertheim thermodynamic perturbation theory. The theory is based on the Pavlyukhin parametrisations [Y. T. Pavlyukhin, J. Struct. Chem. 53, 476 (2012)] of their simulation data for the first four perturbation terms in the high temperature expansion of the Helmholtz free energy of square-well monomer fluids combined with a second-order perturbation theory for the contact value of the radial distribution function of the square-well monomer fluid that enters into bonding contribution. To obtain the latter perturbation terms, we have performed computer simulations in the hard-sphere reference system. The importance of the perturbation terms beyond the second-order one for the monomer fluid and of the approximations of different orders in the bonding contribution for the chain fluids in the predicted equation of state, excess energy and liquid–vapour coexistence densities is analysed.     

地层油气的相平衡热力学研究

地层油气的相平衡热力学研究王利生（北京理工大学化工学院北京１０００８１）关键词混合规则，状态方程，油藏流体，相平衡１链状分子的普遍化范德华配分函数链状分子的转动和振动运动会受阻于其它分子的运动，因而其配分函数中的转动和振动贡献与温度、密度均有关，被分...     

Thermodynamic properties of hard ellipsoid square-well fluid

The perturbation theory with non-spherical reference system is used for molecular fluid with angle-dependent square-well type potential. Simple analytic expressions are given for the thermodynamic properties such as the equation of state, excess free energy per particle, internal energy and internal heat capacity. The effects of anisotropy on the thermodynamic properties are discussed. The anisotropy effects increase with increase of density and decrease of temperature and depends on the anisotropy parameterx ₀.     

Gibbs ensemble simulation of the vapour-liquid equilibrium of square well spherocylinders

Gibbs ensemble Monte Carlo simulations have been performed for systems of square-well spherocylinders of different length-to-breadth ratio. The results are used to test a recent perturbation theory proposed for this kind of system. In addition, the results are compared to similar simulations performed for a Kihara fluid of elongated molecules. An unexpected good agreement is found for the coexistence thermodynamic and structural properties of both model fluids, hence suggesting that the hard spherocylinder plus square-well interaction should be considered as a reference potential for a perturbative treatment of more complex fluid models.     

A self-consistent approach for modelling the interfacial properties and phase diagrams of Yukawa,Lennard-Jones and square-well fluids

A self-consistent density-functional approach is presented for describing the phase behaviour and interfacial tensions of van der Waals fluids represented by the hard-core Yukawa (HCY), Lennard-Jones (LJ) and square-well (SW) potentials. The excess Helmholtz energy functional is formulated in terms of a modified fundamental measure theory (MFMT) for the short-ranged repulsion and a density-gradient expansion for the van der Waals attractions. Analytical expressions for the direct correlation functions of uniform fluids are utilized to take into account the effect of van der Waals’ attraction on intermolecular correlations. For bulk phases, the density functional theory is reduced to an equation of state (EOS) that provides accurate saturation pressures and vapour–liquid phase diagrams. Near the critical region, the long-range fluctuations can be corrected by using the renormalization group (RG) theory. With the same set of molecular parameters, the theory also yields satisfactory surface tensions and interfacial density profiles at all relevant temperatures.     

Quantum corrections to the thermodynamic properties of a two-dimensional fluid with hard-core potential

A simple theory, based on the physical interpretation of the reciprocal of activity, is developed to evaluate the thermodynamic properties of a two-dimensional fluid in the semi-classical limit. The theory is applied to calculate the quantum corrections to the equation of state and excess free energy of two-dimensional fluids, whose molecules interactvia hard-disc and square-well potential. It is found that the quantum effect increases with the increase of density and decrease of temperature.     

The thermodynamics of molecules with discrete potentials

Fluids formed by molecules interacting with discrete potentials are examined in the context of perturbation theory and the reference hypernetted chain equation (RHNC) solution to the Ornstein—zernike equation. A perturbation theory for discrete-potential fluids (DPT) is presented, which only requires one to know the properties of a square-well fluid of variable range. Several potentials are studied: square-shoulder, a combination of a square-well and square-shoulder, and a discrete representation of a continuous potential model. We have found that the DPT approach reproduces the RHNC predictions in most of the cases.     

A general local composition and coordination number model for square-well fluids with variable well width and diameter ratio

A radial distribution function for attractive hard-core systems is obtained from the equilibrium of a molecular pair between local and bulk environments. With this function, a general model is established for the coordination number (CN) and local composition (LC) of square-well fluids. It meets the low-density, high-density and high-temperature limit conditions, as well as the unlike pair conservation and quasi-chemical equilibrium conditions. It also has some other features that many other models do not have: (1) its CN and LC expressions contain all pair potentials; (2) it yields temperature-dependent CN and LC for closely packed mixtures with different pair potentials; (3) its energy parameter is the difference of the total potentials of one pair in local and bulk environments, not the difference of two pair potentials. This model can accurately predict the CN, LC and compressibility factors of square-well fluids from computer simulation over a wide range of density, well width (λ?=?1–2) and diameter ratio. For the case λ?=?1.5, this model is better than or comparable with semi-empirical models; in other cases, it is far better than semi-empirical models. It does not need any empirical parameter for LC prediction. For the prediction of CN and compressibility factors, it only needs the smoothed radial distribution function of pure hard-sphere fluids. It also gives excellent results for lattice gases and highly nonideal lattice mixtures.     

接枝聚合物对小分子的选择性吸附研究

文章应用密度泛函理论研究接枝于壁面的方阱链对二元小分子混合物的选择性吸附特性. 系统的Helmholtz剩余自由能泛函被表示为硬球排斥和方阱吸引两部分贡献之和,分别由硬球链流体状态方程和变阱宽方阱链流体状态方程的简单加权密度近似来进行计算. 用此理论方法,分别考察了接枝聚合物的结构性质,以及不同温度下接枝分子层对二元方阱流体的选择性吸附性能. 结果表明:分子刷厚度随接枝密度线性增长而随温度非线性增加,并且在高温下趋于饱和;在较低温度下,接枝聚合物刷能表现出很好的选择性吸附能力,当聚合物刷被加热到高于饱和温度时,该能力将大幅度地减弱. 关键词： 密度泛函理论 接枝聚合物 选择性吸附 方阱链     

Molecular dynamic simulations and global equation of state of square-well fluids with the well-widths from λ = 1.1 to 2.1

We present the results of extensive new molecular dynamic (MD) simulations in the one-phase region for square well fluids with well widths λ?=?1.10, 1.15, 1.20, 1.25, 1.375, 1.50, 1.75, 1.90, 2.0, and 2.10. These data have been used in developing a crossover equation of state (CR EOS) for square-well fluids with well widths 1.1?≤?λ?≤?2.1. The CR EOS incorporates non-analytic scaling laws in the critical region, and in the limit of low densities yields the exact second and third virial coefficients. Also in the high-temperature region, it reproduces first-order perturbation theory results. The CR EOS was tested against our new MD simulations, and earlier MD and Monte-Carlo (MC) simulations reported by other authors as well. Excellent agreement between calculated values and simulation data for all SW fluids is observed. In combination with the density-functional theory, the CR EOS is also capable of reproducing surface tension simulations with high accuracy. Application of the CR EOS for solid–liquid equilibrium calculations in combination with the Lennard–Jones and Devonshire cell model for the solid phase, is also discussed.     

The effect of steepness of soft-core square-well potential model on some fluid properties

The effect of repulsive steepness of the soft-core square well (SCSW) potential model on the second virial coefficient, critical behaviour (two- phase region and the position of critical point), and coordination number are investigated. The soft-core thermodynamic perturbation theory (TPT) presented by Weeks-Chandler-Anderson (WCA) recently developed by Ben-Amotz and Stell (BAS) has been used for the reference system, and the Barker-Henderson TPT for the perturbed system. The Barker-Henderson macroscopic compressibility approximation has been used for all order perturbation terms in which the second-order one is improved by assuming that the molecules in every two neighbouring shells are correlated upon the original assumption. By using the hard-sphere isothermal compressibility consistency for the radial distribution function (RDF), an analytical closed expression has been derived for the Helmholtz free energy function contained effective hard-sphere diameter. The accuracy of the model has been examined for the hard-core system, and an appropriate range found for the attractive width of the potential well (R), then the effect of steepness parameter on the critical quantities, coordination number, and the inversion temperature of the second virial coefficient, has been investigated qualitatively. The predicted results are in good agreement with the computer simulation data for the critical constants, and coordination number at the limit of the hard-core square-well potential model at least qualitatively, and for the attractive range 1.55 ≤ R ≤ 1.7, quantitatively. It was found that the steepness of the potential model has a marginal effect on the critical behaviour, and also every thermodynamic quantity at low and medium temperatures for which the molecular penetration is negligible, but since the penetration at high temperatures is significant, the role of the steepness of potential on the inversion temperature of the second virial coefficient and coordination number is highlighted.     

Perturbation theory for thermodynamic properties of a ν-dimensional square-well fluid

The Barker-Henderson perturbation theory is used for a ν-dimensional fluid with square-well potential. Analytic expressions are given for the equation of state, excess free energy per particle and internal energy. The numerical results are discussed. A significant feature is the increase of the thermodynamic properties with increasing dimensionality.     

Perturbation theory for the square-well dimer fluid

An analytical equation of state is presented for the square-well dimer fluid of variable well width (1 ≤ λ ≥ 2) based on Barker-Henderson perturbation theory using the recently developed analytical expression for radial distribution function of hard dimers. The integral in the first- and the second-order perturbation terms utilizes the Tang, Y and Lu, B. C.-Y., 1994, J. chem. Phys., 100, 6665 formula for the Hilbert transform. To test the equation of state, NVT and Gibbs ensemble Monte Carlo simulations for square-well dimer fluids are performed for three different well widths (λ = 1.3, 1.5 and 1.8). The prediction of the perturbation theory is also compared with that of thermodynamic perturbation theory in which the equation of state for the square-well dimer is written in terms of that of square-well monomers and the contact value of the radial distribution function.     

用微扰理论推算超临界流体比定容热容

本文利用统计热力学中的BH微扰理论建立数学模型,推算了Ar超临界区域比定容热容,与文献值比较平均偏差为9.54%。引入了修正函数,重新对比热容公式进行推导,再次对Ar超临界区域部分比定容热容进行计算,计算平均偏差0.39%。对CO、CO₂、CH₄等其余11种简单气体的超临界比定容热容进行计算,平均偏差在1%左右。最后给出了计算中各物质势能参数值。     

Perturbation theory for square-well plus sutherland (∞, 6) fluid

The perturbation approach to the theory of fluids, originally introduced by Zwanzig, is applied to the square-well plus Sutherland (∞, 6) fluid. It is shown that the virial coefficients up to the fourth are in close agreement with those for the LJ (12-6) with a hard-sphere cutoff when the bowl width of the square well used is chosen to be equal to 0.25 σ where σ is the hard-sphere diameter. Also, for this suitable choice of the bowl width, the first and the second-order corrections to the free energy have been found to agree reasonably well with those for the square-well fluid of Barker and Henderson.     

方阱链流体在固液界面分布的密度泛函理论研究

应用Yethiraj的加权密度近似泛函理论研究平板狭缝中方阱链流体的密度分布，系统的Helm holtz自由能泛函分为理想气体的贡献利剩余贡献两部分，其中剩余贡献部分分别采用刘洪 来等人建立的基于空穴相关函数的方阱链流体状态方程和Gil-Villegas等人提出的统计缔合 流体理论状态方程(SAFT-VR)结合简单加权密度近似计算.考察了不同链长、温度、系统密度 和壁面吸引强度下平板狭缝中方阱链流体的密度分布，并与Monte Carlo(MC)模拟结果进行 了比较.结果表明采用不同的状态方程对密度分布的计算有明显的影响，对于受限于硬壁狭 缝中的方阱链流体，温度和密度比较高时，两种状态方程计算的结果均与MC模拟符合得比较 好，在低温和低密度下效果变差，SAFT-VR方程的计算结果更接近于MC模拟结果.对于受限于 方阱壁狭缝中的方阱链流体，由于系统密度分布的非均匀性加强，采用两种状态方程计算的 结果均与MC模拟结果有一定偏差，寻找更合适的权重函数是进一步改进的关键. 关键词： 密度泛函理论 非均匀流体 密度分布 固液界面 方阱链