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

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

密度泛函与分子模拟计算介孔孔径分布比较

用巨正则系综Monte Carlo模拟(GCMC)方法和密度泛函理论（ DFT）结合统计积分方程（SIE）计算了介孔材料的孔径分布.为比较这两种方法,以77 K氮气在介孔活性碳微球中的吸附数据为依据,求出其孔径分布.在GCMC模拟和DFT计算中,流体分子模型化为单点的Lerrnard-Jones球;流体分子与吸附剂材料之间的作用采用平均场理论中的10-4-3模型.在DFT方法中,自由能采用Tarazona 提出的加权近似密度泛函方法(weighted density approximation,WDA)求解.结果表明,对于孔径大于1.125 nm的介孔材料,GCMC和DFT两种方法都可以用来研究介孔材料的孔径分布;对于小于1.125 nm的介孔材料,不能用DFT方法计算孔径分布（DFT方法本身的近似产生了误差）,只能用分子模拟方法.     

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

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

平行硬壁间荷电硬球流体的密度泛函理论研究

基于经典流体的密度泛函理论并结合改进的基本度量理论, 研究了受限于平行硬壁间的荷电硬球流体的平衡密度分布. 通过对比有、无外电场时分子数密度分布的相应变化情况, 分析了荷电硬球分子间的库仑排斥势能与场致势能间的竞争作用, 研究了电场对于受限荷电硬球流体系统聚集态结构的影响.     

氢气在碳纳米管内吸附的密度泛函理论研究

本文应用密度泛函理论研究纯流体氢气在单壁碳纳米管内吸附过程,采用硬球状态方程改进的基本测量理论表征硬球的斥力作用,离子吸引项的贡献则用微扰理论描述.在温度为300k,氢气本体对比密度范围为0.2～0.7条件下,计算了三种不同尺寸的碳纳米管氢气吸附的密度分布,其密度泛函计算结果与计算机分子模拟数据完全一致.     

约束条件下的硬球流体

利用密度泛函理论和分子动力学方法,对处于两平行硬墙之间的硬球流体的密度分布进行了计算通过比较两种方法的结果,发现在墙之间距离较大时,Rosenfeld密度泛函理论的结果与分子动力学模拟的结果符合很好;当两堵墙间的距离很小时,这两个结果之间存在明显的不一致.另外,还研究了约束条件下密度分布的结构.     

甲烷晶体的晶格能和弹性性质: 不同方法及泛函的评估

通过对甲烷晶体进行结构、晶格能和弹性特性的研究, 评估了不包含和包含色散能量修正的密度泛函理论的性能. 我们分别利用不包含色散能量修正的密度泛函理论(DFT) (包含不同的标准泛函和杂化泛函)和包含色散能量修正的密度泛函理论(DFT-D)计算了甲烷晶体特性, 并与实验作对比. 尽管DFT-D 与传统密度泛函理论及杂化密度泛函理论相比, 修正了甲烷晶体中的范德华(vdW)相互作用, 但是一些修正方案过分修正了这种相互作用. 因此, 人们在使用DFT-D方法时务必谨慎.     

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

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

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

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

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

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

Further test of third order + second-order perturbation DFT approach: hard core repulsive Yukawa fluid subjected to diverse external fields

Grand canonical Monte Carlo simulation is used to investigate density profiles of hard-core repulsive Yukawa (HCRY) model fluid under the influence of various external fields and radial distribution function (RDF) of the bulk HCRY system. The aim of these extensive simulations is to provide exact data for purely repulsive interaction potential against which the validity of a third order + second-order perturbation DFT approach can be tested. It is found that a semiempirical parametrized bridge function due to Malijevsky and Labik performs very well for the RDF of the bulk HCRY fluid. Incorporation of a bulk second-order direct correlation function (DCF) of the HCRY fluid based on the Malijevsky-Labik bridge function into the third order + second-order perturbation DFT approach yields the resulting theoretical predictions for the density profiles of inhomogeneous HCRY fluid that are in a very good agreement with the simulation data, an exception being somewhat larger deviations appearing for the structure of the fluid around the center of a hard spherical cavity. Both theory and simulation predict layering transition and gas-liquid coexistence phenomena occurring with the HCRY model fluid under confined conditions. For the case of an inverse sixth-power repulsive potential under the influence of a flat stationary wall defined by an inverse twelfth-power repulsive potential, the present third order + second-order perturbation DFT approach is found to be superior to several existing weighted density approximations (WDA) and partitioned WDA.     

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 of inhomogeneous attractive and repulsive hard-core yukawa fluid: grand canonical Monte Carlo simulation and density functional theory study

A density functional theory is proposed for an inhomogeneous hard-core Yukawa (HCY) fluid based on Rosenfeld's perturbative method. The excess Helmholtz energy functional is derived from a modified fundamental measure theory for the hard-core repulsion and a quadratic functional Taylor expansion for the long-ranged attractive or repulsive interactions. To test the established theory, grand canonical ensemble Monte Carlo simulations are carried out to simulate the density profiles of attractive and repulsive HCY fluid near a wall. Comparison with the results from the Monte Carlo simulations shows that the present density functional theory gives accurate density profiles for both attractive and repulsive HCY fluid near a wall. Both the present theory and simulations suggest that there is depletion for attractive HCY fluid at low temperature, but no depletion is found for repulsive HCY fluid. The calculated results indicate that the present density functional theory is better than those of the modified version of the Lovett-Mou-Buff-Wertheim and other density functional theories. The present theory is simple in form and computationally efficient. It predicts accurate radial distribution functions of both attractive and repulsive HCY fluid except for the repulsive case at high density, where the theory overestimates the radial distribution function in the vicinity of contact.     

Structures and adsorption of binary hard-core Yukawa mixtures in a slitlike pore: grand canonical Monte Carlo simulation and density-functional study

The grand canonical ensemble Monte Carlo simulation and density-functional theory are applied to calculate the structures, local mole fractions, and adsorption isotherms of binary hard-core Yukawa mixtures in a slitlike pore as well as the radial distribution functions of bulk mixtures. The excess Helmholtz energy functional is a combination of the modified fundamental measure theory of Yu and Wu [J. Chem. Phys. 117, 10156 (2002)] for the hard-core contribution and a corrected mean-field theory for the attractive contribution. A comparison of the theoretical results with the results from the Monte Carlo simulations shows that the corrected theory improves the density profiles of binary hard-core Yukawa mixtures in the vicinity of contact over the original mean-field theory. Both the present corrected theory and the simulations suggest that depletion and desorption occur at low temperature, and the local segregation can be observed in most cases. For binary mixtures in the hard slitlike pore, the present corrected theory predicts more accurate surface excesses than the original one does, while in the case of the attractive pore, no improvement is found in the prediction of a surface excess of the smaller molecule.     

Density functional theory of homopolymer mixtures confined in a slit

A density functional theory (DFT) is developed for polymer mixtures with shorted-ranged attractive interparticle interactions confined in a slit. Different weighting functions are used separately for the repulsive part and the attractive part of the excess free energy functional by applying the weighted density approximation. The predicted results by DFT are in good agreement with the corresponding simulation data indicating the reliability of the theory. Furthermore, the center-of-mass profiles and the end-to-end distance distributions are obtained by the single chain simulation; the predictions also agree well with simulation data. The results reveal that both the attraction of the slit wall and the temperature has stronger effect on longer chains than on shorter ones because the intrasegment correlation of chains increases with increasing chain length.     

Hybrid density functional theory for homopolymer mixtures confined in a selective nanoslit

By integrating polymer density function theory (DFT) and single-chain molecular simulation, a hybrid DFT is developed for homopolymer mixtures confined in a selective nanoslit. Two weighting functions are adopted separately in the polymer DFT for repulsive and attractive contributions to the excess free energy functional. The theoretical results agree well with simulation data for the density profiles, configurations (tail, loop and train), adsorption amounts, layer thicknesses, and partition coefficients. The polymer-slit interaction is found to have a large effect on the density profiles and partition coefficients but is found to have a small effect on the average sizes and percentages of the configurations. Nearly half of the polymer segments form tails, and the other half form trains. In addition, bridges are observed to form for sufficiently long polymer chains. As the length difference between two polymers increases, the effect of chain connectivity becomes increasingly important.