计算机模拟在材料研究中的应用 - 钱学森空间技术实验室.pdf

钱学森空间技术实验室研究通讯 计算机模拟在材料研究中的应用 ——非晶合金尺寸效应的分子动力学模拟研究 张 琪 （钱学森空间技术实验室空间技术与应用基础研究部，北京100094） 非晶合金又被称作金属玻璃，在结构上具有短程有序、长程无序的特点。与 传统的氧化物玻璃相比，由于非晶合金原子间是通过金属键结合，而不是共价键， 所以保留了与金属相关的特性。与传统的金属材料相比，由于非晶态结构没有位 错和晶界等缺陷，使得非晶合金具有高强度、高弹性极限、良好的耐腐蚀性能、 良好的耐磨性等优异性能。此外，非晶合金还具有热塑性，这使得非晶合金具有 良好的可加工性，可被加工制备成各种形状复杂的零件。 目前，铁磁非晶合金由于具有优异的磁学性能，已逐步取代传统的硅钢片， 被大量用作变压器的磁芯。在刀具、手机、医疗植入器件中也越来越多地看到非 晶合金的身影。在军事应用中，锆基非晶合金由于具有良好的自锐性，可被应用 于穿甲弹弹芯；利用高硬度的特性，非晶合金也可用作穿甲防护材料。在航天航 空应用中，非晶合金由于其物理性能和力学性能在-200 ℃～400 ℃的温度范围内 受温度变化影响很小，可以很好地适应空间温度的变化。利用非晶合金高比强度、 高弹性极限以及低热膨胀系数等特性可将非晶合金制备成盘压杆，用作卫星探测 器的撑杆；利用非晶合金良好的可加工性以及良好的耐磨、耐腐蚀性，可将非晶 合金制备成空间活动部件用的齿轮等零件；利用非晶合金的高弹性极限可将非晶 合金应用于传感器的部件。 作为一种新兴的材料，有关非晶合金的研究还处于发展阶段，关于其结构和 性能的研究还存在很多有争议的话题，尺寸效应就是其中之一。目前的研究中既 有材料强度随尺寸减小而增加的，也有随尺寸减小而降低的，还有几乎不随尺寸 减小发生变化的。这种争议既可能来源于非晶合金本身的性质，也可能来源于制 备方法、实验条件等外加因素的影响。 小尺寸材料的实验研究很大程度上受到制备方法和检测手段的限制，而计算 机模拟实验则可以避免这些限制。通过计算机模拟实验，可以模拟制备出几何形 状完美的纳米线材料，也可以避免离子辐照对样品表面结构的破坏。此外，计算 机模拟实验还可以直接从原子层面进行微观结构的分析。利用量子力学、统计物 理、固体物理等学科的基本原理，计算机模拟实验已经发展出第一性原理计算、 蒙特卡洛方法、分子动力学方法等一系列模拟方法。其中，第一性原理计算方法 以量子力学为基础，虽然相对比较准确，但精确求解薛定谔方程需要耗费大量的 机时，其应用仅仅限制在规模小、结构简单的体系，还难以用到包含上万个原子 的体系，也难以进行长时间的模拟。分子动力学方法，利用给定的经验势函数， 通过求解体系中所有粒子的牛顿运动方程，得到各粒子的运动轨迹，虽然相比第 1 钱学森空间技术实验室研究通讯 一性原理准确性稍差，但是可以应用于上十万、上百万甚至上千万个原子的庞大 体系，随着经验势函数的不断发展，其准确性和可靠性也在不断得到改善。 我的工作采用分子动力学模拟方法对非晶合金的尺寸效应进行了一系列研 究，最近的研究工作发表在Scientific Reports 上，其它相关工作发表在Journal of Chemical Physics、Applied Physics Letters 、Acta Materialia 杂志。 首先采用分子动力学方法分别模拟物理实验中的聚焦粒子束切割（focused ion beam，FIB ）、热压印及铸造方法制备Cu64Zr36 非晶纳米线。聚焦离子束切 割方法的模拟通过在室温下从块体非晶中删除纳米线结构以外的原子实现；热压 印方法的模拟通过在过冷液相温区从块体非晶中删除纳米线结构以外的原子，然 后将其冷却至室温实现；铸造方法的模拟通过施加排斥力，限制液态纳米线边界 的移动，将液态纳米线冷却至室温实现。然后通过调整Cu64Zr36 块体非晶及非晶 纳米线中的原子坐标位置实现单轴拉伸和压缩实验的模拟。应变速率为 108 s-1， 通过控制调整原子坐标位置的时间间隔实现。分子动力学模拟中的应变速率远高