Roth R, Evans R, Dietrich S
H.H. Wills Physics Laboratory, University of Bristol, Bristol BS8 1TL, United Kingdom.
Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Topics. 2000 Oct;62(4 Pt B):5360-77. doi: 10.1103/physreve.62.5360.
We present a versatile density functional approach (DFT) for calculating the depletion potential in general fluid mixtures. For the standard situation of a single big particle immersed in a sea of small particles near a fixed object, the system is regarded as an inhomogeneous binary mixture of big and small particles in the external field of the fixed object, and the limit of vanishing density of the big species, rho(b)-->0, is taken explicitly. In this limit our approach requires only the equilibrium density profile of a one-component fluid of small particles in the field of the fixed object, and a knowledge of the density independent weight functions which characterize the mixture functional. Thus, for a big particle near a planar wall or a cylinder or another fixed big particle, the relevant density profiles are functions of a single variable, which avoids the numerical complications inherent in brute force DFT. We implement our approach for additive hard-sphere mixtures, comparing our results with computer simulations for the depletion potential of a big sphere of radius R(b) in a sea of small spheres of radius R(s) near (i) a planar hard wall, and (ii) another big sphere. In both cases our results are accurate for size ratios s=R(s)/R(b) as small as 0.1, and for packing fractions of the small spheres eta(s) as large as 0.3; these are the most extreme situations for which reliable simulation data are currently available. Our approach satisfies several consistency requirements, and the resulting depletion potentials incorporate the correct damped oscillatory decay at large separations of the big particles or of the big particle and the wall. By investigating the depletion potential for high size asymmetries we assess the regime of validity of the well-known Derjaguin approximation for hard-sphere mixtures and argue that this fails, even for very small size ratios s, for all but the smallest values of eta(s) where the depletion potential reduces to the Asakura-Oosawa potential. We provide an accurate parametrization of the depletion potential in hard-sphere fluids, which should be useful for effective Hamiltonian studies of phase behavior and colloid structure. Our results for the depletion potential in a hard-sphere system, with a size ratio s=0.0755 chosen to mimic a recent experiment on a colloid-colloid mixture, are compared with the experimental data. Although there is good overall agreement, in particular for the dependence of the oscillations on eta(s), there are some significant differences at high values of eta(s).
我们提出了一种通用的密度泛函方法(DFT),用于计算一般流体混合物中的排空势。对于单个大粒子浸没在固定物体附近小粒子海洋中的标准情况,该系统被视为在固定物体外场中大小粒子的非均匀二元混合物,并明确考虑了大粒子密度消失的极限情况,即ρ(b)→0。在此极限下,我们的方法仅需要固定物体场中小粒子单组分流体的平衡密度分布,以及表征混合物泛函的与密度无关的权重函数。因此,对于靠近平面壁、圆柱体或另一个固定大粒子的大粒子,相关密度分布是单个变量的函数,这避免了蛮力DFT中固有的数值复杂性。我们将该方法应用于加和硬球混合物,将我们的结果与计算机模拟结果进行比较,模拟的是半径为R(b)的大球体在半径为R(s)的小球体海洋中靠近(i)平面硬壁和(ii)另一个大球体时的排空势。在这两种情况下,对于小尺寸比s = R(s)/R(b)低至0.1以及小球体的填充分数η(s)高达0.3时,我们的结果都是准确的;这些是目前可获得可靠模拟数据的最极端情况。我们的方法满足几个一致性要求,并且由此产生的排空势在大粒子之间或大粒子与壁之间的大间距处包含正确的阻尼振荡衰减。通过研究高尺寸不对称性下的排空势,我们评估了硬球混合物中著名的Derjaguin近似的有效性范围,并认为即使对于非常小的尺寸比s,除了η(s)非常小的值(此时排空势简化为朝仓 - 大泽势)之外,该近似都不成立。我们提供了硬球流体中排空势的精确参数化,这对于相行为和胶体结构的有效哈密顿量研究应该是有用的。我们将硬球系统中排空势的结果(尺寸比s = 0.0755,选择该值以模拟最近关于胶体 - 胶体混合物的实验)与实验数据进行了比较。虽然总体上有很好的一致性,特别是对于振荡对η(s)的依赖性,但在η(s)较高值时存在一些显著差异。