初级碰撞原子能量对bcc

1、初级碰撞原子能量对bcc-Fe辐照损伤缺陷数的影响摘要：材料处于辐照条件下会产生缺陷，这些微观缺陷积累会在宏观上有所体 现。本文选用不同能量的原子向bcc-Fe表面进行碰撞，会产生位移级联现象。本 文利用基于LAMMPS程序的分子动力学模拟方法，研究不同能量的初级碰撞原子 对bcc-Fe缺陷数量的影响。结果表明，辐照损伤会对材料产生一定的缺陷；本文 主要研究了在辐照损伤过程中，初级碰撞原子能量分别为700ev、1Kev、2Kev对 材料缺陷数量的影响。初级碰撞原子能量对最大缺陷数目产生了一定的影响，当 初级碰撞原子能量为700ev时，最大缺陷数目（间隙原子数目=空位数目）为45 个；初级碰撞原

2、子能量为1Kev时，最大缺陷数目（间隙原子数目=空位数目）为 75个；当初级碰撞原子能量为2Kev时，最大缺陷数目（间隙原子数目=空位数目） 为155个；随着初级碰撞原子能量的增加，缺陷数量随之增加，但初级碰撞原子 能量对体系达到稳定状态时的缺陷数目影响不大。关键词：分子动力学、bcc-Fe、初级碰撞原子、缺陷引言随着社会的飞速发展，能源的需求量逐渐增加，但地球上的能源已经出现了 短缺现象，以及生态破坏、环境污染等问题的出现，开发新能源已成为新的发展 趋势，而新能源高利用率尤为重要。核能的开发以及利用已经引起了世界各国专 家的关注，尤其是处在恶劣环境下的材料所受到的影响，所以材料的抗辐照性能

3、决定了聚变堆的使用寿命。在高温高压等恶劣的辐照环境下，选用抗辐照、使用寿命长的材料尤为重要， A508-III是国产反应堆容器用钢，具有抗脆化性1。辐照脆化能够减少材料的使 用寿命，影响脆化程度的因素有很多，如温度、材料中参杂杂质等。近年来，对 于辐照材料的研究得到了广泛地发展。王2等人研究了空位浓度对bcc-Fe辐照损 伤的影响；马3等人用分子动力学模拟了 3C-SiC辐照诱发缺陷的演化。本文利用分子动力学模拟方法，基于LAMMPS软件进行模拟，研究了不同能 量的初级碰撞原子对bcc-Fe材料缺陷数量的影响，从而得到初级碰撞原子能量与 缺陷数量之间的关系。辐照损伤简介点阵原子在受到具有能量的

4、粒子撞击时，核反应将产生嬗变元素，所产生的 缺陷和嬗变元素会引起材料宏观物理性质的变化，这种改变叫做辐照效应。辐照 效应包括两个方面，辐照损伤和辐照改性。辐照损伤是粒子撞击材料产生的缺陷 使材料性能发生改变。在粒子撞击材料到材料达到稳定状态的过程中，可分为两 个阶段：1.初级过程：具有能量的粒子撞击材料时，入射原子会与材料原子产生 弹性碰撞，入射原子会把其中一部分能量传递给被撞击的材料原子。这些被撞击 的材料原子叫做初级碰撞原子，简PKA（primary knock-on atom）。如图1所示。如 果材料原子被带有能量的入射粒子击出，则被击出的材料原子叫做初级离位原子。 2.次级过程：带有能

5、量的初级离位原子可以使其他材料原子离位，从而形成耳机 碰撞原子，而带有能量的耳机碰撞原子又可以使三级离位原子击出，如此延续， 就形成了级联碰撞过程。模型及模拟方法本文利用分子动力学方法模拟bcc-Fe在不同初级能量原子的碰撞下产生位移 级联过程，并运用可视化软件OVITO4直接观察其演化过程。本文采用F-S势函 数5，对bcc-Fe在三维方向上进行堆垛，能量为700ev和1Kev模拟空间为 25a0 x25a0 x25a0, 2Kev 模拟空间为 40a0 x40a0 x40a0 其中 a0=0.287nm 为晶格常数， 初级碰撞原子向Z轴负方向进行碰撞，采用周期性边界条件，模拟体系温度为 3

6、00k，时间步长选取1fs，初级碰撞原子能量分别选为700eV，1KeV，2KeV，分 别对三个模型进行模拟，采用npt系综对体系进行弛豫，使得体系达到平衡状态。模拟结果与讨论3.1初级碰撞原子能量为700eV在体系模型已经建立完全之后，对体系进行能量最小化弛豫，使体系达到稳 定状态，图2和图3分别表示当初始碰撞原子能量为700ev时，体系温度为300k 时初始平衡状态，各个原子会在其平衡位置振动，不会脱离其平衡位置。当给定 初级碰撞原子700ev能量时，体系会发生位移级联现象而产生间隙原子和空位， 图4表示体系在整个位移级联过程中，产生最大的空位和间隙原子时的状态，从 图中可以看出存在许多间

7、隙原子，通过可视化软件OVITO中W-S缺陷分析的到最 大缺陷数目为45个；随着时间的推移，体系在退火过程中缺陷数量减小，逐渐 趋于稳定状态，图5表示体系退火过程中最终所达到的稳定状态，其空位数和间 隙原子数均为3个。最终状态的缺陷数目要比最大缺陷数目小得多。图3初级碰撞原子能量为700ev，初始平衡状态三维立体图图5初级碰撞原子能量为700ev，体系趋于稳定时x-y平面图3.2初级碰撞原子能量为1KeV当给定初级碰撞原子1Kev能量时，体系会发生位移级联现象而产生间隙原子 和空位，图6表示体系在整个位移级联过程中，产生最大的空位和间隙原子时的 状态，从图中可以看出存在许多间隙原子，通过可视化

8、软件OVITO中W-S缺陷分 析的到最大缺陷数目为75个；随着时间的推移，体系在退火过程中缺陷数量减 小，逐渐趋于稳定状态，图7表示体系退火过程中最终所达到的稳定状态，其空 位数和间隙原子数均为5个。最终状态的缺陷数目要比最大缺陷数目小得多。通 过比较初级碰撞原子能量为700ev和1Kev时，产生最大缺陷数目相差较大，说 明初级碰撞原子能量对其影响较大，比较最终状态缺陷数目，说明初级碰撞原子 能量对其影响不大。图6初级碰撞原子能量为1Kev，缺陷最大时x-y平面图图7初级碰撞原子能量为1000ev，体系趋于稳定时x-y平面图3.3初级碰撞原子能量为2Kev由于初级碰撞原子能量较大，此模型采用4

9、0a0 x40a0 x40a0，在体系模型已经 建立完全之后，对体系进行能量最小化弛豫，使体系达到稳定状态，图8和图9 分别表示当初始碰撞原子能量为2Kev时，体系温度为300k时初始平衡状态，此 时，各个原子会在其平衡位置振动，不会脱离其平衡位置。当给定初级碰撞原子 1Kev能量时，体系会发生位移级联现象而产生间隙原子和空位，图10表示体系 在整个位移级联过程中，产生最大的空位和间隙原子时的状态，从图中可以看出 存在许多间隙原子，通过可视化软件OVITO中W-S缺陷分析的到最大缺陷数目为 155 个；随着时间的推移，体系在退火过程中缺陷数量减小，逐渐趋于稳定状态， 图11表示体系退火过程中最

10、终所达到的稳定状态，其空位数和间隙原子数均为9 个。最终状态的缺陷数目要比最大缺陷数目小得多。通过比较初级碰撞原子能量 为700ev、1Kev、2Kev时，产生最大缺陷数目相差较大，说明初级碰撞原子能量 对其影响较大，比较最终状态缺陷数目，说明初级碰撞原子能量对其影响不大。图9初级碰撞原子能量为2kev，初始平衡状态三维立体图图11初级碰撞原子能量为2kev，体系趋于稳定时x-y平面图结论bcc-fe在受到辐照损伤时，能够产生位移级联从而使体系产生缺陷(空位、间 隙原子)。用过对体系进行模拟研究，在温度为300k，初级碰撞原子能量为 700ev时，利用W-S缺陷分析法得到最大缺陷数目为45个，

11、最终缺陷数目为3 个。初级碰撞原子能量为1Kev时，利用W-S缺陷分析法得到最大缺陷数目为75 个，最终缺陷数目为5个。初级碰撞原子能量为2Kev时，利用W-S缺陷分析法 得到最大缺陷数目为155个，最终缺陷数目为9个。通过比较以上数据得出，最 大缺陷数目相差较大，说明初级碰撞原子能量对其影响较大，且随着初级碰撞原 子能量的增加，最大缺陷数目逐渐增大；比较最终状态缺陷数目，说明初级碰撞 原子能量对其影响不大。参考文献李承亮，张明乾.压水堆核电站反应堆压力容器材料概述J.材料导报,2008,22(9):65-68.王建伟，尚新春，吕国才.bcc-Fe空位浓度对辐照损伤影响的分子动力学模拟J.材 料工程,2011(10):15-18.马小强，袁大庆,夏海鸿，等.3C-SiC辐照诱发缺陷演化及温度效应分子动力学模拟 J.原子能科学技术,2016,50(2) :219-226.Stukowski A. Visualization and analysis of atomisticsimulation data with OVITO- the Open Visualization ToolJ