FeCr合金级联碰撞的分子动力学模拟 .docxVIP

FeCr合金级联碰撞的分子动力学模拟郑颖，杨新华（华中科技大学力学系，武汉 430074）摘要：采用分子动力学模拟的方法，在中子辐射条件下模拟了含不同Cr浓度的FeCr合金的级联碰撞过程。结果表明，在纯铁和不同Cr浓度下的FeCr合金中平均残留的弗兰克尔缺陷对数近似的相同；空位团簇所占残留空位的比例在某一范围内波动，与Cr浓度无关；在Cr浓度低于5%时自间隙团簇所占自间隙原子的比例稍高Cr浓度低于7%时的情况。研究了关于残留的自间隙原子中Cr原子所占比例的情况，数据表明，不同Cr浓度的FeCr合金级联碰撞后残留的自间隙原子中，Cr原子所占的比例是模拟前FeCr合金中Cr的浓度的2-3倍左右，说明Cr原子在自间隙缺陷中有聚集趋势。关键词：FeCr合金，级联碰撞，分子动力学引言核聚变反应堆由于运行温度高、产生的中子通量大等严酷的工况对结构的材料提出了苛刻的要求。实验表明含2%~6%Cr的铁素体/马氏体（Ferritic/Martensitic, F/M）钢在强辐射下具有几何稳定性、较低的辐射肿胀和优良的耐腐蚀性能等特性，因此被认为是未来聚变示范堆、ADS和四代堆的重要候选结构材料[1]。关于核聚变堆材料辐射损伤的问题一直是研究的重点。而对由级联碰撞导致的初始损伤状态的研究是中子辐射损伤模拟的出发点。早在1995年Phythian等[2]利用分子动力学模拟方法模拟了纯Fe和纯Cu在能量处于60eV-10keV范围的入射粒子以及不同温度100K-900K条件下对级联碰撞过程的影响。在以往的FeCr合金相关研究中，Malerba等[3]运用分子动力学模拟方法研究了含10%Cr的FeCr合金的与纯铁的级联碰撞过程。Vortler等[4]模拟了在同一温度和不同入射粒子能量下两种分别含5%Cr和15%Cr的FeCr合金级联碰撞过程的差异。为了对比FeCr合金中不同Cr浓度对级联碰撞过程的影响，本文采用分子动力学的方法针对辐射损伤的级联碰撞过程进行了相关研究。1分子动力学模拟方法分子动力学模拟中至关重要的一步是原子间的力的计算，这有赖与原子间的作用势，可以说分子动力学模拟的准确度取决于原子间的作用势准确度。本文采用Caro等[5]文中所给的浓度依赖嵌入原子势法(CD-EAM)。同嵌入原子势法(EAM)相似，具有N个原子的FeCr合金系统总势能表示为： (1)其中表示为原子类型（Fe或Cr），是原子类型的嵌入函数，是原子和间的距离。是所有其他原子的核外电子在第个原子处产生的电子云密度之和，表达式为：(2)表示对势，对于FeCr合金来说有：、和。表示局部浓度，表达式为： (3)是原子处电子云密度之和（）中Fe原子产生的部分。是四阶的Rcdlich-Kister多项式，当时，。为了更准确地表示高能原子间短程相互作用，本文采用ZBL势函数[6]与FeCr合金CD-EAM势函数相结合的方式来描述级联碰撞阶段原子间的相互作用。模拟过程采用Verlet算法进行动力学积分和Berendsen热浴法[7]进行温度控制。2FeCr合金模型和模拟步骤FeCr合金是体心立方结构，晶格常数取，模拟的几何尺寸为（共250000个原子）。在级联碰撞模拟中分别取了6中不同Cr浓度的FeCr合金：Fe-2%Cr、Fe-5%Cr、Fe-7%Cr、Fe-10%Cr、Fe-15%Cr和Fe-20%Cr。结果与纯铁的级联碰撞结果作为对比。以Fe-2%Cr和Fe-15%Cr合金为例，模型如图1所示。建立模型的方法是：先建立纯Fe 体系模型，再按相应比例对原子类型进行随机替换，其中淡蓝色的是Fe 原子，深蓝色的是Cr 原子。(a) Fe-2%Cr (b)Fe-15%Cr图1Fe-2%Cr和Fe-15%Cr合金三维模型Fig.1 3D models of Fe-2%Cr and Fe-15%Cr alloy级联碰撞模拟采用的是分子动力学软件LAMMPS。模拟中边界条件为周期性边界条件，模拟温度设为600K，采用三维周期性边界条件。弛豫阶段，整个区域采用NVE 系综，并施加Benerdsen热浴，弛豫时的步长为4fs，时长为12ps；级联碰撞阶段，整个模拟过程采用NVE 系综[8-10]。为了让碰撞区域在中心附近，避免碰撞区域与边界相互作用，在弛豫后稳定体系的中心区域选择一个Fe 原子，并赋予5keV 的能量，作为初始离位原子(PKA)。同时为避免沟渠效应，选择沿135方向入射[8-10]。时间步长的变化由PKA入射时的0.01fs到级联碰撞模拟末期时的1fs。并且对纯Fe和各Cr浓度的FeCr合金模拟了8次级联碰撞过程。本文采用WS原胞法[11]（Wigner-Seitz cell method）来鉴别缺陷，具体方法如下：将弛豫足够时间的模拟盒子中的原子构型作为参照