新起点
周期性边界条件
2021-01-27 11:50:58

周期性边界条件(PBCs)是分子模拟与数学模型中的常用近似方法之一,其思想是用一个称为元胞的周期性盒子来描述宏观的体系。在一个元胞周围有紧密堆积的完全相同的元胞。例如在二维体系中,一个方形元胞四周有相同的元胞(称为"镜像"),若一个粒子从右边边界穿出,操作上左边边界对应位置有相同粒子进入,由此保证了元胞中粒子数目即物理量如动量与能量的守恒。

尽管计算机的计算能力发展迅速,但用其模拟阿佛加德罗常数数量级的粒子仍然是不现实的。当今计算机能模拟数百万乃至千万个粒子的体系,要用这有限个粒子反映宏观相的热力学性质,只能借助于周期性边界来实现,使用有限个粒子的体系近似地描述宏观体系。分子动力学模拟中,常使用PBC描述气体、液体、溶液及晶态的行为。

PBC算法主要体现在两部分。其一:对体系内粒子坐标的约束。x_size是一个方向上的元胞边长,则粒子坐标区间为。例如,超方形晶格的周期性边界条件对应着高维下的简单立方堆积。使用最密堆积(英语:Close-packing of equal spheres)构建高维下的周期性边界条件元胞可获得理论最优效率。例如,四维的最密堆积是正十六胞体堆砌,八维下则是E8晶格(英语:E8 lattice)。高维周期性边界条件的实现等价于信息论中基于高维密堆积的纠错码算法。

PBC的引入不影响体系的动量守恒与能量守恒,但是角动量不再守恒。

PBC作为一种人为的近似,不可避免地对模拟引入偏差,影响计算结果的有效范围。即使对于最简单的兰纳-琼斯流体,PBC的引入对NVE系综下的模拟的取样也有轻微影响,从而对物理量的统计值引入系统误差,故而引入PBC的NVE系综被建议称作NVEPG系综。径向分布函数g()在大于半个盒子边长时是不准确的,等于盒子边长时更是有人为的尖峰。在固体体系中,涉及振动的研究,声音,冲击波或声子的波长范围受体系大小的限制。对于长程力例如电磁相互作用的计算,需要使用特殊的方法(如Edwald加和,PPPM方法等)。相行为研究中,连续相的取向由于界面能的影响常常沿着盒子边长,体系受限甚至可能得出错误的结论,因此即使引入PBC,体系的大小也应当足够大以尽可能减小引入近似造成的非物理干扰。

相关:

  • 分子建模
  • 网站公告: