分子动力学概述.docxVIP

分子动力学 分子动力学方法是一种计算机模拟实验方法，是研究凝聚态系统的有力工具。该技术不仅可以得到原子的运动轨迹，还可以观察到原子运动过程中各种微观细节。它是对理论计算和实验的有力补充。 分子动力学总是假定原子的运动服从某种确定的描述，这种描叙可以牛顿方程、拉格朗日方程或哈密顿方程所确定的描述，也就是说原子的运动和确定的轨迹联系在一起。在忽略核子的量子效应和Born-Oppenheimer绝热近似下，分子动力学的这一种假设是可行的⑴。所谓绝热近似也就是要求在分子动力学过程中的每一瞬间电子都处于原子结构的基态。要进行分子动力学模拟就必须知道原子间的相互作用势。在分子动力学模拟中，我们一般采用经验势来代替原子间的相互作用势，如Lennard-Jones势、Mores势、EAM原子嵌入势、F-S多体势。然而采用经验势必然丢失了局域电子结构之间存在的强相关作用信息，即不能得到原子动力学过程中的电子性质［1］。 事实上，分子动力学就是模拟原子系统的趋衡过程。实际上，分子动力学方法就是确定某一描述与初始条件、边值关系的数值解。我们假定系统经过M步长之后达到稳定，而这一稳定状态正是我们所求的。 1、分子动力学的算法分析 首先，我们假定我们研究的系统服从Newton方程所确定的描述，即： j:(t)=-F(t)(1)m 式中r(t)表征原子在t时刻的位置矢量 F(t)表征原子在t时刻所受到的力，