Rickayzen G, Heyes D M
Division of Chemistry, School of Biomedical and Molecular Sciences, University of Surrey, Guildford GU2 7XH, United Kingdom.
J Chem Phys. 2007 Mar 21;126(11):114504. doi: 10.1063/1.2647150.
The authors investigate the behavior of a model fluid for which the interaction energy between molecules at a separation r is of the form 4epsilon[(sigma/r)2n-(sigma/r)n], where epsilon and sigma are constants and n is a large integer. The particular properties they study are the pressure p, the mean square force F2, the elastic shear modulus at infinite frequency Ginfinity, the bulk modulus at infinite frequency Kinfinity, and the potential energy per molecule u. They show that if n is sufficiently large it is possible to derive the properties of the system in terms of two parameters, the values of the cavity function and of its derivative at the position r=sigma. As an example they examine in detail the cases with n=144 and n=72 for three different temperatures and they test the theory by comparison with a computer simulation of the system. They use the simulated pressure and the average mean square force to determine the two parameters and use these values to evaluate other properties; it is found that the theory produces results which agree with computer simulation to within approximately 3%. It is also shown that the model, when the parameter n is large, is equivalent to Baxter's sticky-sphere model with the strength of the adhesion determined by the value of n and the temperature. They use Baxter's solution of the Percus-Yevick equations for the sticky-sphere model to determine the cavity function and from that the values of the same properties. In this second approach there are no free parameters to determine from simulation; all properties are completely determined by the theory. The results obtained agree with computer simulation only to within approximately 6%. This suggests that for this model one needs a better approximation to the cavity function than that provided by the Percus-Yevick solution. Nevertheless, the model looks promising for the study of (typically small) colloidal liquids where the range of attraction is short but finite when compared to its diameter, in contrast to Baxter's sticky-sphere limit where the attractive interaction range is taken to be infinitely narrow. The continuous function approach developed here enables important physical properties such as the infinite shear modulus to be computed, which are finite in experimental systems but are undefined in the sticky-sphere model.
作者研究了一种模型流体的行为,对于该流体,分子间在间距为r时的相互作用能具有4ε[(σ/r)^2n - (σ/r)^n]的形式,其中ε和σ为常数,n为大整数。他们所研究的特定性质包括压力p、均方力F²、无限频率下的弹性剪切模量G∞、无限频率下的体积模量K∞以及每个分子的势能u。他们表明,如果n足够大,就有可能根据两个参数来推导系统的性质,这两个参数即腔函数及其在r = σ位置处的导数的值。作为一个例子,他们详细研究了n = 144和n = 72在三种不同温度下的情况,并通过与该系统的计算机模拟进行比较来检验该理论。他们使用模拟压力和平均均方力来确定这两个参数,并使用这些值来评估其他性质;结果发现该理论得出的结果与计算机模拟结果在约3%的范围内相符。还表明,当参数n很大时,该模型等同于巴克斯特的粘性球模型,其粘附强度由n的值和温度决定。他们使用巴克斯特对粘性球模型的珀库斯 - 耶维克方程的解来确定腔函数,并由此确定相同性质的值。在第二种方法中,无需从模拟中确定自由参数;所有性质完全由该理论确定。所得到的结果与计算机模拟结果仅在约6%的范围内相符。这表明对于该模型,需要一个比珀库斯 - 耶维克解所提供的更好的腔函数近似。尽管如此,与巴克斯特的粘性球极限情况(其中吸引相互作用范围被视为无限窄)相比,该模型对于研究(通常较小的)胶体液体看起来很有前景,在胶体液体中,吸引力范围与其直径相比短但有限。这里开发的连续函数方法能够计算重要的物理性质,如无限剪切模量,这些性质在实验系统中是有限的,但在粘性球模型中是未定义的。