Díez A, Largo J, Solana J R
Departamento de Física Aplicada, Universidad de Cantabria, E-39005 Santander, Spain.
J Phys Chem B. 2007 Aug 30;111(34):10194-201. doi: 10.1021/jp072059z. Epub 2007 Aug 7.
Monte Carlo simulations have been performed to determine the excess energy and the equation of state of fcc solids with Sutherland potentials for wide ranges of temperatures, densities, and effective potential ranges. The same quantities have been determined within a perturbative scheme by means of two procedures: (i) Monte Carlo simulations performed on the reference hard-sphere system and (ii) second-order Barker-Henderson perturbation theory. The aim was twofold: on the one hand, to test the capability of the "exact" MC-perturbation theory of reproducing the direct MC simulations and, on the other hand, the reliability of the Barker-Henderson perturbation theory, as compared with direct MC simulations and MC-perturbation theory, to determine the thermodynamic properties of these solids depending on temperature, density, and potential range. We have found that the simulation data for the excess energy obtained from the two procedures are in close agreement with each other. For the equation of state, the results from the MC-perturbation procedure also agree well with the direct MC simulations except for very low temperatures and extremely short-ranged potentials. Regarding the Barker-Henderson perturbation theory, we have found that in general the second-order approximation does not provide significant improvement over the first-order one.
已经进行了蒙特卡罗模拟,以确定具有萨瑟兰势的面心立方固体在很宽的温度、密度和有效势范围下的过剩能量和状态方程。通过两种方法在微扰方案中确定了相同的量:(i) 在参考硬球系统上进行的蒙特卡罗模拟,以及(ii) 二阶巴克-亨德森微扰理论。目的有两个:一方面,测试“精确”的蒙特卡罗微扰理论再现直接蒙特卡罗模拟的能力;另一方面,与直接蒙特卡罗模拟和蒙特卡罗微扰理论相比,测试巴克-亨德森微扰理论确定这些固体依赖于温度、密度和势范围的热力学性质的可靠性。我们发现,从这两种方法获得的过剩能量的模拟数据彼此非常吻合。对于状态方程,除了非常低的温度和极短程的势之外,蒙特卡罗微扰方法的结果也与直接蒙特卡罗模拟非常吻合。关于巴克-亨德森微扰理论,我们发现一般来说,二阶近似并不比一阶近似有显著改进。