这组 MATLAB 例程模拟具有周期性边界条件的三维硬球堆积。 它初始化一个均匀晶格，然后使用 Metropolis 算法对粒子位置进行多次迭代退火。 然后计算 $g_2(r)$，或成对径向分布函数，在集合中的许多中心上取平均值。 对于大的初始粒子密度，结果 $g_2$ 与 Percus-Yevick 近似给出的解析近似匹配，给定足够高的粒子密度的任意精度。
2022-02-26 14:34:44
4KB
1