可变元胞分子动力学方法对晶体结构稳定性的研究

本文利用可变元胞的分子动力学方法（MD），研究了离子束辐照条件下的晶体结构稳定性问题，计算证实了铜、银、铝等金属在常压下，面心结构是稳定的。在三维分子动力学程序中，使用了跳点法差分格式，并对原格式作了修正。计算表明，这样的差分格式较常用的Euler-Cauchy法和四阶Runge-Kutta法节省机时，且占用较小的贮存单元。