分子动力学模拟方法概述.pdfVIP

《装备制造技术》 年第 期 2007 10 分子动力学模拟方法概述 申海兰,赵靖松 (西安电子科技大学机电工程学院,陕西 西安 ) 710071 摘要:介绍了分子动力学模拟的基本原理及常用的原子间相互作用势,如Lennard-Jones势;论述了几种常用的有限差分算法,如Verlet 算法;说明了分子动力学模拟的几种系综及感兴趣的宏观统计量的提取。 关键词:分子动力学模拟;原子间相互作用势;有限差分算法;系综 中图分类号: 文献标识码: 文章编号: O3 A 1672-545X(2007)10-0029-02 从统计物理学中衍生出来的分子动力学模拟方法(molec- , ),实践证明是一种描述纳米科技 ulardynamicssimulationMDS 2势函数 研究对象的有效方法,得到越来越广泛的重视。所谓分子动力 学模拟,是指对于原子核和电子所构成的多体系统,用计算机 ( )势 2.1Lennard-JonesL-J 模拟原子核的运动过程,从而计算系统的结构和性质,其中每 [5] 最先提出了两体势的解析形式, 将它 Mie Lennard-Jones 一个原子核被视为在全部其他原子核和电子所提供的经验势 应用到铜的自扩散研究计算中。该两体势可表示为 [1] 场作用下按牛顿定律运动 。它被认为是本世纪以来除理论分 σ 12 σ 6 ( ) ( ) #() - () 析和实验观察之外的第三种科学研究手段,称之为“计算机实 r=4ε$ r r % 3 [2] 式中 和 是势参数,表示原子间距。通常为更好地 验”手段 ,在物理学、化学、生物学和材料科学等许多领域中得 ε σ r 拟合元素已有的实验数据,如结合能、晶格常数等, 势则使 到广泛地应用。 L-J [6] 根据模拟对象的不同,将它分为平衡态分子动力学模拟 用一种普适的形式 。 (