纳米单晶氩扩散性质的分子动力学模拟

采用分子动力学(MD)方法对纳米单晶氩的扩散性质进行了模拟研究。采用L-J12-6势函数描述氩原子势,半步长的Leap-frog算法求解运动方程,Berendsen热浴法调温,P-R方法调压,在充分弛豫后进行扩散模拟的时间为100 ps,计算得到了10个恒定温度下原子均方位移(MSD)与时间(t)的关系曲线。根据MSD-t曲线特点,分为初始阶段和稳态阶段2个阶段进行数据处理。利用Einstein关系求出了稳态阶段的扩散系数D,由Arrhenius表达式求出了纳米单晶氩的稳态扩散常数D0为0.100 004 509 m2/s、稳态扩散激活能Q为4.848 81×10-21J。