铁团簇的半经验原子间相互作用对势的研究.pdfVIP

铁团簇的半经验原子间相互作用对势的研究 马永庆，王逊 大连海事大学 (116026) email:myq@ 摘 要：由量子化学从头计算方法得出不同距离时两个原子间的相互作用能，以势能曲线最 低点与量子化学从头计算的最小值相重合的拟合原则，拟合出原子间相互作用对势。并尝试 用原子间相互作用对势研究铁团簇 （奥氏体、铁素体和马氏体）的原子平均结合能和最近平 衡原子间距。计算结果表明，这种近似方法对理解金属晶体的相变和许多性能可以提供证据。 如此推断，可以计算 Fe-Cr-Mn-C-B 系复杂合金的奥氏体的能量变化，并为解释相变的趋势 提供参考。 关键词：原子间相互作用势；量子化学从头计算；拟合对势；铁团簇 1．引言 由第一性原理量子力学所作的计算可以给出精确的结果。然而，由于需要很大的计算量 来精确求解Schrödinger方程，目前这些方法仅限于研究包含几十个原子的体系。原子间相 互作用势是分子动力学、蒙特卡罗方法模拟计算的基础，也是在原子、分子层次上进行合金 设计的重要部分。原子间相互作用势有很多种，按相互作用的类型划分为两体势（即对势、 偶对势）、三体势、四体势、N 体势、嵌入原子势等。物质（包括各种相，如固相、液相和 气相）的许多物理性质既可以通过实验观察和测定，也可以在理论上计算和预测。另一方面， 经验的原子间相互作用势虽不及上述方法精确，但能处理大得多的体系。对于复杂的材料科 [1] 。本文尝试用原 学问题的计算机模拟，原子间相互作用势及分子动力学方法占有重要地位 子间相互作用对势研究奥氏体、铁素体和马氏体铁团簇的原子平均结合能和最近平衡原子间 距，为进行Fe-Cr-Mn-C-B系实用合金的复杂体系的特性的预测提供参考。 2．原子间相互作用对势 原子间相互作用对势的简单函数是Lennard-Jones 势。尽管Lennard-Jones 型函数仅仅可 以用来描述多体势作用最小的稀有气体，而对于象金属晶体这样多体势作用强的晶体被认为 是不适用的，但是，作为对复杂体系合金的初步探讨，这种方法很可能对理解许多性能和过 程可以提供有益的参考。 [2] Lennard-Jones 势的形式为 ：  σ  m  σ  n  U r ( ) 4ε   −   ， （1） ij    r   r    本课题得到高等学校博士学科点专项科研基金（项目编号：20030151006）资助 交通部应用基础研究项目（项目编号：200332922503）资助 - 1 - 式中m 和n 通常取为 12 和6，ε和σ 将由量子化学从头计算方法拟合确定。求出N 个原子的 总势能 N U (R) ∑U (p R) ， （2） t ij ij i 1, ij 式中R 为原子间最近距离，p ij R 为第i、j 原子间距离。由 dU t (R) 0