平行板间方势阱偶极流体的结构性质研究

在密度泛函理论(DFT)框架下, 应用改进的基本度量理论(MFMT)表达硬球作用对自由能泛函的贡献, 根据统计力学理论结合加权密度近似(WDA)表达偶极作用对自由能泛函的贡献,得到了方势阱偶极流体在平行板间的密度分布表达式, 计算了偶极流体在两平行板间的密度分布, 并探讨了方势阱深度和宽度对体系密度分布的影响. 此外, 通过体系密度分布, 进一步分析了方势阱宽度和深度以及板间尺度与溶剂化力的关系.     

平行板间方势阱偶极流体的结构性质研究

在密度泛函理论(DFT)框架下, 应用改进的基本度量理论(MFMT)表达硬球作用对自由能泛函的贡献, 根据统计力学理论结合加权密度近似(WDA)表达偶极作用对自由能泛函的贡献,得到了方势阱偶极流体在平行板间的密度分布表达式, 计算了偶极流体在两平行板间的密度分布, 并探讨了方势阱深度和宽度对体系密度分布的影响. 此外, 通过体系密度分布, 进一步分析了方势阱宽度和深度以及板间尺度与溶剂化力的关系.     

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

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

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

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

约束条件下的硬球流体

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

非迭代冻结密度近似方法在计算氢键相互作用的合理性研究

为计算相互作用较弱的分子碎片之间的耦合能,Harris从密度泛函理论出发,提出了一种简化方法,即冻结密度近似(FDA)方法.对该方法在描述分子间氢键作用的合理性进行了验证.对水分子间的HO┉H氢键、甲酰胺与水分子间的NH┉O氢键、二氟甲烷和水分子间的OH┉F氢键,以及DNA中的碱基(AT,GC)之间的N—H┉O,N—H┉N等类型的氢键的计算表明:若电子交换关联采用非定域自旋密度近似,FDA的计算结果同其他abinitio方法的计算结果以及实验结果都符合得很好.FDA在计算过程中既不需要求解泊松方程,也不需要进行反复的自洽迭代,所以运算速度较快,在研究生物大分子体系中的氢键相互作用方面具有一定的使用价值 关键词： 冻结密度近似 氢键 密度泛函     

3C-SiC(001)-(2×1)表面原子与电子结构研究

采用广义梯度近似的密度泛函理论方法计算了3C-SiC（001）-（2×1）表面的原子及电子结构．计算结果表明,3C-SiC（001）-（2×1）表面为非对称性的Si二聚体模型,其二聚体的Si原子间键长为0.232 nm．电子结构的计算结果表明,在费米能级处有明显的态密度,因此3C-SiC（001）-（2×1）表面呈金属性．在带隙附近存在四个表面态带,一个位于费米能级附近,一个位于费米能级以上5 eV处,另外两个位于费米能级以下的价带中． 关键词： 碳化硅 密度泛函理论计算 原子结构 电子结构     

氮气在MCM-41中毛细凝聚的模拟和理论研究

采用巨正则系综 Monte Carlo 方法(GCMC)以及基于扩展基本度量理论和 MBWR 状态方程的密度泛函理论,研究了 77.4 K 时氮气在不同孔径的 MCM-41 分子筛中吸附的密度分布和吸附等温线.提出的平均场权重密度泛函理论,克服了平均场近似预测主体相热力学性质时的偏差,获得的密度分布和吸附等温线与 GCMC 分子模拟结果有着很好的一致性.     

半导体性单壁碳纳米管与其石墨衬底表面相互作用的第一性原理计算

利用基于密度泛函理论的计算方法,研究了处于石墨（HOPG）表面的半导体性单壁碳纳米管,发现碳纳米管下面的石墨受碳管的作用向下凹陷,而纳米管本身虽然保持其形状不变,但它的电子态受石墨衬底影响,造成导带下移,禁带宽度明显减小. 关键词： 单壁碳纳米管 密度泛函理论 局域密度近似方法     

含时密度泛函理论及应用的最新发展

密度泛函理论在材料计算研究领域得到了广泛的应用,然而它无法处理含时问题和材料的激发态性质。Runge-Gross定理奠定了含时密度泛函理论的基础,为研究这两类问题提供了有效的手段。经过三十多年的发展,含时密度泛函已被应用到量子化学、材料计算等多个领域,人们也更加了解其优势和不足。目前,含时密度泛函理论和方法仍在迅速发展。本文简要回顾含时密度泛函方法的发展历史,介绍近年来含时密度泛函在理论和应用方面的一些重要进展,总结当前在含时密度泛函领域存在的重要难题以及面临的挑战,展望其发展方向和趋势。     

Mean Spherical Approximation-Based Partitioned Density Functional Theory

Previous literature claims that the density functional theory for non-uniform non-hard sphere interaction potential fluid can be improved on by treating the tail part by the third order functional perturbation expansion approximation (FPEA) with the symmetrical and intuitive consideration-based simple function C0^(3)(r1, r2, r3) =(∫dr4a(r4-r1)a(r4-r2)a(r4-r3) as the uniform third order direct correlation function (DCF) for the tail part,here kernel function a(r) = (6/πσ^3)Heaviside(σ/2 - r). The present contribution concludes that for the mean spherical approximation-based second order DCF, the terms higher than second order in the FPEA of the tail part of the non-uniform first order DCF are exactly zero. The reason for the partial success of the previous a kernel function-based third order FPEA for the tail part is due to the adjustable parameter ξ and the short range of the a kernel function.Improvement over the previous theories is proposed and tested.     

Bridge density functional approximation for non-uniform hard core repulsive Yukawa fluid

In this work, a bridge density functional approximation （BDFA） （J. Chem. Phys. 112, 8079 （2000）） for a nonuniform hard-sphere fluid is extended to a non-uniform hard-core repulsive Yukawa （HCRY） fluid. It is found that the choice of a bulk bridge functional approximation is crucial for both a uniform HCRY fluid and a non-uniform HCRY fluid. A new bridge functional approximation is proposed, which can accurately predict the radial distribution function of the bulk HCRY fluid. With the new bridge functional approximation and its associated bulk second order direct correlation function as input, the BDFA can be used to well calculate the density profile of the HCRY fluid subjected to the influence of varying external fields, and the theoretical predictions are in good agreement with the corresponding simulation data. The calculated results indicate that the present BDFA captures quantitatively the phenomena such as the coexistence of solid-like high density phase and low density gas phase, and the adsorption properties of the HCRY fluid, which qualitatively differ from those of the fluids combining both hard-core repulsion and an attractive tail.     

Density functional approximation for van der Waals fluids： based on hard sphere density functional approximation

A universal theoretical approach is proposed which enables all hard sphere density functional approximations (DFAs) applicable to van der Waals fluids. The resultant DFA obtained by combining the universal theoretical approach with any hard sphere DFAs only needs as input a second-order direct correlation function (DCF) of a coexistence bulk fluid, and is applicable in both supercritical and subcritical temperature regions. The associated effective hard sphere density can be specified by a hard wall sum rule. It is indicated that the value of the effective hard sphere density so determined can be universal, i.e. can be applied to any external potentials different from the single hard wall. As an illustrating example, the universal theoretical approach is combined with a hard sphere bridge DFA to predict the density profile of a hard core attractive Yukawa model fluid influenced by diverse external fields; agreement between the present formalism's predictions and the corresponding simulation data is good or at least comparable to several previous DFT approaches. The primary advantage of the present theoretical approach combined with other hard sphere DFAs is discussed.     

Augmented Kierlik-Rosinberg Fundamental Measure Functional andExtension of Fundamental Measure Functional to Inhomogeneous Non-hard Sphere Fluids

From point of view of weighted density procedure, it isguessed that a Percus-Yevick (PY) compressibility excess free energydensity, appearing in the Kierlik--Rosinberg type fundamentalmeasure functional (KR-FMF) and expressed in terms of scaledparticle variables, can be substituted by a corresponding expressiondictated by a more accurate Mansoori-Carnahan-Starling-Leland(MCSL) equation of state, while retaining the original weightingfunctions; it is numerically indicated that the resultantundesirable non-self-consistency between the PY type weightingfunction and MCSL type excess free energy density had no bad effecton the performance of the resultant augmented KR-FMF which, on theone hand, preserves the exact low-density limit of the originalKR-FMF and holds a high degree of pressure self-consistency, on theother hand, improves significantly, as expected, the predictions ofdensity profile of hard sphere fluid at single hard wall contactlocation and its vicinity, and of the bulk hard sphere second orderdirect correlation function (DCF), obtained from functionaldifferentiation. The FMF is made applicable to inhomogeneousnon-hard sphere fluids by supplementing a functional perturbationexpansion approximation truncated at the lowest order with summationof higher order terms beyond the lowest term calculated by the FMFfor an effective hard sphere fluid; the resultant extended FMF onlyneeds second order DCF and pressure of the fluid considered atcoexistence state as inputs, consequently is applicable whether theconsidered temperature is above critical point or below criticalpoint. The extended MCSL-augmented KR-FMF is found to be endowedwith an excellent performance for predictions of density profile andsurface tension by comparing the present predictions of these twoquantities with available computer simulation data for inhomogeneoushard core attractive Yukawa fluid and Lennard-Jones fluid.     

Statistical Mechanics Approach for Uniform and Non-uniform Fluid with Hard Core and Interaction Tail

One recently proposed self-consistent hard sphere bridge functional was combined with an exponential function exp(-cr) and a re-normalized indirect correlation function to construct the bridge function for fluid with hard core and interaction tail. In the present approach, the adjustable parameter α was determined by the thermodynamic consistency realized on the compressibility modulus, the re-normalization of the indirect correlation function was realized by a modified Mayer function with the interaction potential replaced by the perturbative part of the interaction potential. As an example, the present bridge function was combined with the Ornstein-Zernike (OZ) equation to predict structure and thermodynamics properties in very good agreement with the simulation data available for Lennard-Jones (L J). Based on the universality principle of the free energy density functional and the test particle trick, the numerical solution of the OZ equation was employed to construct the first order direct correlation function of the non-uniform fluid as a functional of the density distribution by means of the indirect correlation function. In the framework of the density functional theory, the numerically obtained functional predicted the density distribution of LJ fluid confined in two planar hard walls that is in good agreement with the simulation data.     

Inhomogeneous 2D Lennard-Jones Fluid: Theory and Computer Simulation

In this work, thermodynamical properties of a two-dimensional (2D) Lennard-Jones (LJ) fluid are studied. Here, to increase the accuracy of our theoretical calculations, the correlation functions in three-particle level (triplet) are applied. To obtain the triplet correlation functions, the Attard's source particle method is extended to 2D systems. In the Attard's procedure, the inhomogeneous Ornstein-Zernike (OZ) equation is solved using the Treizenberg-Zwanzwig (TZ) expression and a closure relation like the hypernetted-chain (HNC) approximation. In the present work, we also have performed the Monte Carlo (MC) simulation. The theoretical results are in fairly agreement with the MC simulation. Also, our results show that the approach proposed here is suitable to study the 2D LJ fluid.     

Inhomogeneous 2D Lennard-Jones Fluid: Theory and Computer Simulation

In this work, thermodynamical properties of a two-dimensional (2D) Lennard-Jones (LJ) fluid are studied. Here, to increase the accuracy of our theoretical calculations, the correlation functions in three-particle level (triplet) are applied. To obtain the triplet correlation functions, the Attard's source particle method is extended to 2D systems. In the Attard's procedure, the inhomogeneous Ornstein-Zernike (OZ) equation is solved using the Treizenberg-Zwanzwig (TZ) expression and a closure relation like the hypernetted-chain (HNC) approximation. In the present work, we also have performed the Monte Carlo (MC) simulation. The theoretical results are in fairly agreement with the MC simulation. Also, our results show that the approach proposed here is suitable to study the 2D LJ fluid.     

LJ模型的热力学特性、相平衡及用于真实分子的研究

结合描述硬球固体Helmholtz自由能的自由体积方法与描述硬球固体径向分布函数的拟合的分析表达式与一阶热力学摄动理论,用于描述Lennard Jones(LJ)固体的Helmholtz自由能.按照一个修正的WCA方法将LJ势分为短程排斥部分与长程吸引部分,将文献中一个用于求取液相的等价的硬球直径的简单的迭代法扩展到固相,用于求取固相的等价的硬球直径.在固体Helmholtz自由能的计算中,使用200壳层,以便获得精确的结果.体相LJ液体的热力学特性由一个最近提出的状态方程求取.该方法很好地描述了LJ固体的过量Helmholtz自由能与状态方程,满意地描述了Lennard Jones模型的相平衡;通过选取合适的LJ势参数,能很好地描述了真实分子的融化曲线.     

A classical density functional approach to depletion interaction of Lennard-Jones binary mixtures

In this article, we apply classical density functional theory to investigate the characteristics of depletion interaction in Lennard-Jones (LJ) binary fluid mixtures. First, to confirm the validity of our adopted density functional formalism, we calculate the radial distribution functions using a theoretical approach and compare them with results obtained by molecular dynamics simulation. Then, this approach is applied to two colloids immersed in LJ solvent systems. We investigate the variation of depletion interaction with respect to the distance of two colloids in LJ binary systems. We find that depletion interaction may be attractive or repulsive, mostly depending on the bulk density of the solvent and the temperature of the binary system. For high bulk densities, the repulsive barrier of depletion force is remarkable when the total excluded volume of colloids touches each other and reaches a maximum. The height of the repulsive barrier is related to the parameters of the LJ potential and bulk density. Moreover, the depletion force may exhibit attractive wells if the bulk density of the solvent is low. The attractive well tends to appear when the surface–surface distance of colloids is half of the size of the polymer and deepens with temperature lowering in a fixed bulk density. In contrast with the hard-sphere system, no oscillation of depletion potential around zero is observed.