缔合Lennard-Jones流体在纳米缝隙中的界面性质研究

结合改进的基础度量理论(modified fundamental measure theory, MFMT)和一阶平均球近似理论(first-order mean-spherical approximation theory, FMSA), 建立了二缔合Lennard-Jones(LJ)流体在纳米缝隙中的自由能模型. 应用密度泛函理论(Density functional theory, DFT), 研究了界面张力等温线及汽液平衡时的界面张力, 阐明了外场强度、缔合能量和缝隙宽度对汽液界面张力的影响.     

Yukawa流体的相平衡和界面张力研究

应用二阶微扰理论, Duh-Mier-Y-Teran状态方程和在平均球近似(mean spherical approximation, MSA)的基础上获得的直接相关函数, 建立了适用于均匀流体和非均匀流体的状态方程. 结合此状态方程, 重整化群理论(renormalization group theory, RG)和密度泛函理论(density functional theory, DFT), 分别研究了Yukawa流体的相平衡和界面张力. 结果与分子模拟数据吻合良好.     

AB2型聚合物流体的表面结构性质

在密度泛函理论(DFT)框架下, 应用改进的基本度量理论(MFMT)表达硬球作用对自由能泛函的贡献, 根据统计力学理论结合加权密度近似(WDA)表达聚合作用对自由能泛函的贡献, 建立了描述AB2型聚合物流体的化学势, 得到了聚合物流体在硬球颗粒表面的密度分布表达式, 计算了聚合物流体在硬球颗粒表面附近的密度分布, 并探讨了体积分数、聚合程度和硬球颗粒尺度对体系密度分布的影响. 此外, 通过体系密度分布, 进一步分析了体积分数、聚合程度和硬球颗粒尺度与剩余吸附的关系.     

由L-J流体的FMSA状态方程计算流体的PVT性质

运用Tang等提出的Lennard-Jones (L-J)流体两参数的一阶平均球形近似(FMSA)状态方程, 计算了流体的汽液共存相图和饱和蒸汽压曲线, 以及非饱和区的PVT性质, 并与文献数据进行比较. L-J参数由Tr＜0.95的汽液相共存数据回归得到. 计算结果表明, 对于分子较接近球形的流体, 除临界点附近外, 该方程可以在较大的温度和压力范围内计算真实流体的PVT性质, 结果满意. 对于球形分子, 该方程的精确度随分子尺寸的变大基本保持稳定. 该方程不适用于强极性物质. 在高密度区, 该方程的计算结果明显优于P-R方程. 对于分子偏离球形较远的流体, 该方程的适用性变差, 此时要考虑分子形状的影响, 可采用三参数的FMSA状态方程进行计算.     

平行板间超支化聚合物流体的密度分布和溶剂化力

在密度泛函理论(DFT)框架下,应用改进的基本度量理论(MFMT)和统计力学理论结合加权密度近似(WDA)分别表示硬球作用和聚合作用对自由能泛函的贡献,建立了描述ABg型超支化聚合物流体的化学势,得到了聚合物流体在平行板间的密度分布表达式,计算了聚合物流体在两平行板间的密度分布,并探讨了体积分数、反应程度和单体中B类官能团数对体系密度分布的影响.此外,通过体系密度分布,进一步分析了反应程度、板间宽度与溶剂化力的关系.     

计算机模拟铌原子簇的稳定构型和能量性质

采用密度泛函(DFT)方法结合全局优化"Basin-Hopping"算法研究了铌原子簇: 对于Nb小簇n＝2～6我们用密度泛函方法计算了它们的稳定几何构型和电子结构, 通过拟合计算结果构造铌原子簇势能函数, 并利用该函数和全局优化"Basin-Hopping"算法得到较大铌原子簇(n＝7～20)能量极小的结构. 计算结果表明与实验及其它计算结果相一致.     

基于密度泛函理论研究二元排斥Yukawa流体的表面结构性质

基于自由能密度泛函理论(DFT)考察了二元排斥Yukawa (HCRY)流体在不同外场下的密度分布. 基于微扰理论, 体系的Helmholtz自由能泛函采用硬球排斥部分和长程色散部分贡献之和, 其中Kierlik和Rosinberg的加权密度近似(WDA)被用来计算硬球排斥部分, 而色散部分采用平均场理论(MFT)进行描述. 为了验证DFT计算结果的合理性, 研究中采用巨正则Monte Carlo(GCMC)模拟计算了在不同主体相密度、硬核直径和位能参数比的条件下二元HCRY混合流体的密度分布. 结果表明, 该DFT计算结果与GCMC模拟值吻合良好.     

CuC、CuN分子基态的结构与分析势能函数

在Cu 的有效核势近似下, 运用密度泛函(B3LYP)方法, 对Cu采用基集合LANL2DZ, 但对其价电子基组的p轨道函数部分做了必要的修改, 而对C、N采用基集合6-311+G(d), 对CuC和CuN分子的微观结构进行了理论计算. 优化并计算了两分子基态的能量, 平衡结构和谐振频率. 根据原子分子反应静力学原理, 导出了CuC和CuN分子基态的合理离解极限和离解能. 应用密度泛函(B3LYP)方法扫描了CuC和CuN分子基态的势能曲线, 并采用最小二乘法拟合了两分子基态的Murrell-Sorbie势能函数及其在平衡位置附近的Dunham展开式. 同时根据Herzberg 和Dunham的公式, 计算了CuC和CuN分子基态的光谱参数.     

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

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

用密度泛函理论研究Lennard-Jones 流体在狭缝中的相平衡

用改进的基础度量理论(modified fundamental measure theory, MFMT)和密度Taylor展开分别表达过剩自由能中的短程作用和色散作用. 流体分子与狭缝壁之间的相互作用以10-4-3势能函数表达. 由巨势最小原理确定Lennard-Jones (LJ)流体在狭缝中的密度分布和过剩吸附量, 所得结果与分子模拟数据吻合良好. 根据平衡时两相温度, 化学势及巨势相等, 计算了LJ流体在狭缝中的相平衡.     

First-order mean spherical approximation for inhomogeneous fluids

The first-order mean-spherical approximation (FMSA) [Y. Tang, J. Chem. Phys., 118, 4140 (2003)] is extended to the studies of inhomogeneous fluids by combining with Rosenfeld's perturbative method [Y. Rosenfeld, J. Chem. Phys. 98, 8126 (1993)]. In the extension, the key input-direct correlation function of FMSA-is applied to constructing the free energy density functional. Preserving its high fidelity at the bulk limit, the FMSA shows satisfactory performance for Yukawa fluids near hard and attractive walls. The results are better than or comparable to several other theories reported before for the geometry. The FMSA is found, in particular, more satisfactory than the traditional mean-field theory for predicting density profiles around hard walls. The FMSA is also compared with the full MSA for inhomogeneous fluids, showing no appreciable differences. The inhomogeneous FMSA goes successfully through the self-consistency test for reproducing the radial distribution function of the bulk Yukawa fluid. As far as the computation is concerned, the FMSA can be executed much faster than any nonmean-field theories, and the speed is virtually identical to that of the mean-field theory.     

Theoretical study of Sutherland fluids with long-range, short-range, and highly short-range potential parameters

On the basis of the first-order mean spherical approximation (FMSA) theory the behavior of Sutherland fluids with a number of parameters (gamma=3.1-36) is investigated. The investigation includes its modification by the simplified exponent approximation, renormalization group (RG) transformation, and density functional theory (DFT). For long-range parameters, the original FMSA is found sufficiently good to describe the global phase behavior, including inside the critical region. For short-range parameters, the modified FMSA by the simplified exponent approximation outside the critical region and RG transform inside the critical region are applied. For extremely short-range forces, the success is achieved by its combination with the DFT. This work gives a general sense about the capability of a theory for different ranges of potential, as well as for different thermodynamic states.     

First-order mean-spherical approximation for interfacial phenomena: a unified method from bulk-phase equilibria study

The recently proposed first-order mean-spherical approximation (FMSA) [Y. Tang, J. Chem. Phys. 121, 10605 (2004)] for inhomogeneous fluids is extended to the study of interfacial phenomena. Computation is performed for the Lennard-Jones fluid, in which all phase equilibria properties and direct correlation function for density-functional theory are developed consistently and systematically from FMSA. Three functional methods, including fundamental measure theory for the repulsive force, local-density approximation, and square-gradient approximation, are applied in this interfacial investigation. Comparisons with the latest computer simulation data indicate that FMSA is satisfactory in predicting surface tension, density profile, as well as relevant phase equilibria. Furthermore, this work strongly suggests that FMSA is very capable of unifying homogeneous and inhomogeneous fluids, as well as those behaviors outside and inside the critical region within one framework.     

Prediction of phase behavior of nanoconfined Lennard-Jones fluids with density functional theory based on the first-order mean spherical approximation

The recently proposed first-order mean spherical approximation (FMSA) [Y. Tang, J. Chem. Phys. 121, 10605 (2004)] for inhomogeneous fluids is extended to study the phase behavior of nanoconfined Lennard-Jones fluids, which is consistent with the phase equilibria calculation of the corresponding bulk fluid. With a combination of fundamental measure theory, FMSA provides Helmholtz free energy and direct correlation function to formulate density functional theory, which implementation is as easy as the mean-field theory. Following previous success in predicting density profiles inside slit pores, this work is focused specially on the vapor-liquid equilibrium of the Lennard-Jones fluids inside these pores. It is found that outside the critical region FMSA predicts well the equilibrium diagram of slit pores with the sizes of 5.0, 7.5, and 10 molecular diameters by comparing with available computer simulation data. As a quantitative method, FMSA can be treated as an extension from its bulk calculation, while the mean-field theory is only qualitative, as its bulk version.     

Thermodynamic and structural properties of mixed colloids represented by a hard-core two-Yukawa mixture model fluid: Monte Carlo simulations and an analytical theory

The interaction between colloidal particles is well represented by a hard-core two-Yukawa potential. In order to assess the accuracy of theoretical predictions for the thermodynamic and structural properties of mixed colloids, standard Monte Carlo simulations are carried out for the hard-core two-Yukawa mixtures. In the simulations, one range parameter in the two-Yukawa potential is taken as 1.8 or 2.8647, and another is taken as 4, 8, or 13.5485. Both attractive and repulsive dominant cases of the potential outside the hard core are considered. The effects of temperature, density, composition, size and energy parameter ratios on internal energy, compressibility factor, and radial distribution function are investigated extensively. Theoretical calculations are performed in the framework of analytical solution for the Ornstein-Zernike equation with the first-order mean spherical approximation (FMSA). Our analysis shows that the FMSA is very accurate for the prediction of the compressibility factor of the hard-core two-Yukawa mixtures at all conditions studied. The FMSA generally predicts accurate internal energy, but overestimates the internal energy of the systems at lower temperatures. Furthermore, we found that a simplified exponential version of the FMSA predicts fairly good radial distribution function at contact for the mixed two-Yukawa fluids. The comparison of the theoretical compressibility factor with that from the Monte Carlo simulations suggests that the FMSA can be used to investigate the fluid-fluid equilibria of hard-core two-Yukawa mixtures.     

Structure and adsorption of a hard-core multi-Yukawa fluid confined in a slitlike pore: grand canonical Monte Carlo simulation and density functional study

Because of the increasing interest in studying the phenomenon exhibited by charge-stabilized colloidal suspensions in confining geometry, we present a density functional theory (DFT) for a hard-core multi-Yukawa fluid. The excess Helmholtz free-energy functional is constructed by using the modified fundamental measure theory and Rosenfeld's perturbative method, in which the bulk direct correlation function is obtained from the first-order mean spherical approximation. To validate the established theory, grand canonical ensemble Monte Carlo (GCMC) simulations are carried out to determine the density profiles and surface excesses of multi-Yukawa fluid in a slitlike pore. Comparisons of the theoretical results with the GCMC data suggest that the present DFT gives very accurate density profiles and surface excesses of multi-Yukawa fluid in the slitlike pore as well as the radial distribution functions of the bulk fluid. Both the DFT and the GCMC simulations predict the depletion of the multi-Yukawa fluid near a nonattractive wall, while the mean-field theory fails to describe this depletion in some cases. Because the simple form of the direct correlation function is used, the present DFT is computationally as efficient as the mean-field theory, but reproduces the simulation data much better than the mean-field theory.     

Structure and thermodynamics of hard-core Yukawa fluids: thermodynamic perturbation approaches

The thermodynamic perturbation theories, which are based on the power series of a coupling constant (λ-expansion), have been proposed for studying the structural and thermodynamic properties of a hard-core Yukawa (HCY) fluid: one (A1-approximation) is the perturbation theory based on the hard-sphere repulsion as a reference system. The other (A2-approximation) is the perturbation theory based on the reference system which incorporates both the repulsive and short-range attractive interactions. The first-order mean-spherical approximation (FMSA) provided by Tang and Lu [J. Chem. Phys. 99, 9828 (1993)] has been employed for investigating the thermodynamic properties of a HCY fluid using the alternative method via the direct correlation function. The calculated results show that (i) the A1 and A2 approximations are in excellent agreements with previous computer simulation results in the literature and compare with the semi-empirical works of Shukla including the higher-order free energy terms, (ii) the A1 and A2 approximations are better than the FMSA and the mean-spherical approximation, (iii) the A2-approximation compares with the A1-approximation, even though the perturbation effect of an A2-approximation is much smaller than that of an A1-approximation, and that (iv) the FMSA study is particularly of advantage in providing the structure and thermodynamics in a simple and analytic manner.     

Sedimentation equilibrium of colloidal suspensions in a planar pore based on density functional theory and the hard-core attractive Yukawa model

The sedimentation equilibrium of colloidal suspensions modeled by hard-core attractive Yukawa (HCAY) fluids in a planar pore is studied. The density profile of the HCAY fluid in a gravitational field and its distribution between the pore and uniform phases are investigated by a density functional theory (DFT) approach, which results from employing a recently proposed parameter-free version of the Lagrangian theorem-based density functional approximation (Zhou, S. Phys. Lett. A 2003, 319, 279) for hard-sphere fluids to the hard-core part of the HCAY fluid, and the second-order functional perturbation expansion approximation to the tail part as was done in a recent partitioned density functional approximation (Zhou, S. Phys. Rev. E: Stat. Phys., Plasmas, Fluids, Relat. Interdiscip. Top. 2003, 68, 061201). The resultant DFT approach is, thus, the first adjustable parameter-free DFT for HCAY fluids. The validity of the present DFT for HCAY fluids of reduced range parameter z(red) = 1.8 under various external potentials is established in the first of the papers cited previously. The present DFT for HCAY fluids can predict the radial distribution function for the bulk HCAY fluid accurately in the colloidal limit (large value of z(red)), and in the hard-sphere limit, its prediction for the density profile of the hard-sphere fluid in a gravitational field is in very good agreement with the existing simulation data. The dependence of the density profile and distribution coefficient on the magnitude of the interparticle attraction, gravitational field, and degree of confinement is investigated in detail by the present DFT approach. Intuitive and qualitative analyses are also compared with the quantitative DFT calculational results.