第四讲(08-2a).pptVIP

分子动力学方法（MD）是按系统内禀力动力学规律来计算并确定状态的改变。 3、Gear方法（预纠法） 预纠法的步骤如下： 第一步：根据当前t时刻的位置坐标、速度（动量）、加速度等，预测出下一个时刻的位置坐标、速度（动量）、加速度等值。 第二步：根据预测的位置坐标，计算出力（因为力是坐标的函数），从而得到一个新的加速度，评价预测的加速度和计算得到两个加速度之间的偏差。 第三步：利用此评价偏差，对原第一步预测得到的时刻位置坐标、速度（动量）、加速度等值进行修正。 Gear方法是在t时刻处按泰勒级数展开，得到的 时刻位置坐标、速度（动量）、加速度等的预测值。即 这是第一步，得到了预测值。然后根据 ， 计算出作用力 ，再得到加速度： 评价预测的加速度 的偏差，即： 最后对预测值进行修正。即 由于物质系统的复杂性以及原子间相互作用类型的不同,很 难得到满足各种不同体系和物质的一般性而又精度较高的势函 数.所以针对不同的物质体系人们陆续发展了大量的经验和半经 验的势函数。 半经验势： 形式经过一定的理论分析而得到, 但其中一些参数则需要 根据宏观实验参数用经验方法来确定,这些宏观实验参数主要有 弹性常数、平衡点阵常数以及内聚能、空位形成能和层错能等。 经验势： 为了拟合的方便, 在选择势函数的形式时,并不一定要求有 确切的理论依据,而是出于经验的估计和拟合方便的需要,相对自 由地选择势函数形式。　 一、对势 对势认为原子之间的相互作用是两两原子之间的作用， 与其他原子的位置无关。 1、间断对势（又称硬球模型） 形式如下： 2 、 Lennard-Jones( LJ)势 在实际应用中常采用形式： C反映了相互作用的强度; 反映了原子的大小。 第一项描述了范德瓦耳斯力。第二项是排斥力，其来源之一是 原子核之间的库仑斥力，来源之二是电子之间由于泡利不相容 原理产生的排斥力。 Lennard-Jones势常用来模拟惰性气体的凝聚态，例如对 于象CO2这类分子晶体可以采用原子-原子之间的势能加和来表 示分子之间的作用势。 为了提高计算精度，有时还需对Lennard-Jones势作进一 步修饰。 3、Morse势 Morse势较适宜与双原子分子。也常用来构造多体势中的对势 部分。 4、Born - Mayer势 第一项代表离子之间的库仑势；第二项代表排斥势；第三项代 表偶极子之间的作用势；第四项代表偶极子与四极子之间的作 用势。 Born - Mayer势适用于碱金属卤化物和碱土金属卤化物 等各种离子型化合物。 二、多体势 对势模型在解决分子晶体和离子型化合物问题中取得 了成功，但对于含有共价键的过渡金属却遇到了很多困难。 1、嵌入原子势（EAM势） 思想是：某原子的原子核除了受到其他原子核的排斥作用外， 还受到该原子的核外电子和其他原子的背景电子的静电作用，其 表达式为： 其中 是原子核之间的两体对势， 是把原子i埋入电子 密度为 的位置所需要的能量。 嵌入原子势很好的描述了金属原子之间的相互作用，是描述金 属材料的一种常用势函数。 2、修正的嵌入原子势（MEAM势） 以上讨论的势都是中心对称的，对于由共价键结合的 有机分子和半导体材料是不适宜的。为了更好的描述各种含有 共价键作用的材料，发展了很多考虑角度效应的多体势，其中 修正的嵌入原子势就是比较常用的一种，它的形式与嵌入原子 势一样，只是在背景电子密度函数 中考虑了角度效应，其形 式为： 其中 是分电子密度，分别代表核外s、 p、d 、f 轨道的电子密度。 是分电子密度 的权重因子。 修正的嵌入原子势能精密的描述共价键结合原子间的作 用，但参数拟合复杂，计算量大，不易收敛。 势参数的确定 原子间相互作用的势函数形式多种多样，待势函数的形式 选定后，对于具体研究对象的势参数又不同，所以势参数的确 定在分子动力学模拟中至关重要。 势参数的确定通常有三种方法： A、通过实验值（例如：晶格常数、弹性系数、振动谱等）拟合势参数。 B、通过蒙特卡罗方法确定势参数。 C、基于量子力学知识得到的各种微观信息来确定势参数。 截断半径法： 其特点是：对某一分子，预先选定一个截断半径rc, 只计算以截断半径rc为球体内的分子间的相