关于分子动力学模拟中扩散系数的系统尺寸依赖性：纯稠密流体的简单修正公式

讨论了一个实用的修正公式，该公式将具有周期性边界条件的分子动力学模拟中稠密液体的自扩散系数与热力学极限中的自扩散系数联系起来。该公式适用于纯稠密流体，形式非常简单：$D=D_0(1-\gamma N^{-1/3})$，其中 $D_0$ 是热力学极限中的自扩散系数，$N$ 是模拟中的粒子数。数值因子 $\gamma$ 取决于模拟单元的几何形状。值得注意的是，对于最流行的立方几何形状，$\gamma\simeq 1.0$。分子动力学模拟的结果支持了该公式的成功，包括最近使用“魔法”模拟几何形状的模拟结果。