液态金属大系统急冷过程中微结构转变机制的模拟研究.pdfVIP

第30卷 第2朔 稀有金■材料与工程 V01．30．No．2 RARE 2001年4月 METALMATERIALSANDENGINEERING Apfn200l 液态金属大系统急冷过程中微结构 转变机制的模拟研究 郑采星1 刘让苏1 董科军1 彭 平1 刘海蓉1徐仲榆1谢泉2 (1湖南大学长抄410082) (2长抄交通学院基础科学部长抄410076) 摘要 首次对于由50000个Al原子组成的液态金属系统在由急冷形成非晶态结构的过程中所出现的t观结构组 卷的变化进行了分子动力学模拟研究，探讨各种有关t型、原子团结构组态的转变机制．褥到一■十分清晰的关于 金属原于成固及其演变的躅景．同时也对模拭计算结果的赣观机理给出了相应的讨论，这对于搽人理解非晶卷结构 的形成机制及其徽观过程．将有重要的启示作用。 美t词 藏态金属 般结构转变机制 计算机模撤 中圈法分类号：TGl46．2+1文■标识码：A 文章■号：1002．185x(2001)02—0089-04 众所周知，由液态金属急冷形成的非晶态金属具 开放实验室银河I型大型计算机，采用分子动力学方 有许多优良的性能，主要是由其特殊的微观结构组态 法首次对50000个Al原子的液态金属大系统在急冷 决定的。因而对其非晶形成过程进行跟踪研究，以弄 过程中的结构组态变化进行模拟研究。模拟计算的条 清其微观结构组态的转变特性，对于深人理解其徽观 件为：将50000个Al原子置于一立方体盒中，按周 结构与宏观性能的关系，将具有重要的理论和实践意 期性边界条件运行，其原子间相互作用势仍采用由 义。然而，在目前的实验条件下，要实现跟踪研究是 wang“。1等人所发展的扩展非局域模型庸势理论中 难以完成的。随着计算机技术的飞速发展，我们已将 的双体有效势函数： 分子动力学方法运用于对由500个原子组成的液态 r f O、∞r 1 yo)=(z知^J f l—f詈II如F(口)sin(q)向f 金属小系统的微观结构组态的瞬变过程，进行物理图 L ㈡… J 象十分清晰的模拟研究，并已取得了一些重要的研究 (1) 结果[|“]。但由于在模拟计算中利用了周期性边界条 式中zⅢ和F(口)分别表示有效离子价和归一化能量 件，对于由500个原子组成的小系统而言，处于或接 波数特性函数，其定义详见文献[6，7]。双体势的截 近于边界的原子将达到近300个，约为总原子数的 止距离为20u(原子单位)。运行的时间步长为10_15 60％。这将严重影响其模拟研究的结果。为此，本研 s。急冷速度为3．4×10”K／s。模拟计算从943K(Al 究则在原有工作的基础上采用并行算法，对于由 的熔点为933K)开始，首先让系统等沮运行750步使 50 之处在乎衡态(按系统的能量在运行中的变化程度来 000个原子组成的液态金属大系统(其中处于或接 近于边界的原子约为8000个，约为总原子数的 判断其平衡程度时，该系统在运行100步左右时，实 16％，其边界效应的影响将大为减弱)的急冷凝固过 际上就已经达到了平衡态)。然后再让系统以3．4× 1013 883，833， 程进行