摘要：本文利用python编写了分子动力学模拟程序,并利用该程序对He分子体系进行了详细研究。分别研究了不同初始条件,不同边界条件,截断点位置等的研究。在分子数为500,分子初速度为500的初始条件下,体系最终经过2.25e-11s发展成为麦克斯韦平衡体系。分别研究周期性边界条件与刚性边界对系统发展的影响,研究发现周期性边界条件与刚性边界条件最终对系统的发展没有明显影响。对截断点进行研究后发现,必须选取大于1.1*sigma(sigma为L-J势能中的常数,由实验可测定)作为截断,如果小于该值,体系由于数值误差将偏离物理实际。同时研究发现势函数对于体系的发展也具有重要影响。

关键词：分子动力学 边界条件 verlet算法 he分子体系 

单位：浙江省杭州第二中学; 浙江杭州310053