Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, Tennessee 37235, USA.
J Chem Phys. 2012 Sep 14;137(10):104104. doi: 10.1063/1.4749381.
We present a Fourier space density functional approach for hard particles with attractive interactions, which is based on a previously developed two-dimensional approach [S. Hlushak, W. Rżysko, and S. Sokołowski, J. Chem. Phys. 131, 094904 (2009)] for hard-sphere chains. The interactions are incorporated by means of a three-dimensional Fourier image of the direct correlation function that is obtained from the first-order mean-spherical approximation. In order to improve the computational efficiency, we make extensive use of fast Fourier transforms for calculating density convolution integrals. A two-dimensional implementation of the new density functional approach, based on the expansion of the functional around the bulk fluid density, is used to study structure and adsorption of two model fluids in narrow cylindrical pores. We also investigate two methods that improve the accuracy of the theory as compared to the conventional DFT approach, which expands the free energy functional around the bulk fluid density: One a variant of the reference fluid density functional theory used by Gillespie et al. [Phys. Rev. E 68, 031503 (2003)], and the second a weighted density approach with energy route thermodynamics. Results from these two methods are compared to the conventional approach and also to the results of Monte Carlo simulations. We find that the method of Gillespie et al. and the weighted density approach with energy route thermodynamics yield significant improvement over the conventional approach.
我们提出了一种用于具有吸引力相互作用的硬粒子的傅里叶空间密度泛函方法,该方法基于先前开发的用于硬球链的二维方法 [S. Hlushak、W. Rżysko 和 S. Sokołowski,J. Chem. Phys. 131, 094904 (2009)]。相互作用通过直接相关函数的三维傅里叶像来引入,该直接相关函数是通过第一性平均球近似获得的。为了提高计算效率,我们广泛使用快速傅里叶变换来计算密度卷积积分。新的密度泛函方法的二维实现基于围绕体相密度展开的函数,用于研究窄圆柱形孔中两种模型流体的结构和吸附。我们还研究了两种方法,这些方法与传统的 DFT 方法(即围绕体相密度展开自由能泛函)相比,可以提高理论的准确性:一种是 Gillespie 等人使用的参考流体密度泛函理论的变体 [Phys. Rev. E 68, 031503 (2003)],另一种是具有能量路径热力学的加权密度方法。这两种方法的结果与传统方法以及蒙特卡罗模拟的结果进行了比较。我们发现,Gillespie 等人的方法和具有能量路径热力学的加权密度方法比传统方法有显著的改进。