-
氢(H)同位素滞留问题是聚变堆第一壁材料设计的关键, 而深入理解H在缺陷(如空位)处的非均匀形核长大过程有助于揭示H起泡及滞留的机制. 针对第一壁材料钨(W)中空位捕获H及其解离的动力学过程开展研究, 通过耦合捕获和解离两过程, 构建新物理模型, 避免了原单一过程的物理模型需准确记录相应事件首次发生时间的不足, 另外新模型可同时提取解离系数和有效捕获半径等动力学参数. 通过分子动力学模拟发现新模型能较好地描述W中空位-H复合体对H捕获和解离的动力学过程, 根据空位-H复合体随时间的演化曲线, 提取了有效捕获半径和解离系数等动力学参数. 一方面能为动力学蒙特卡罗和速率理论等长时空尺度方法提供输入参数, 另一方面促进了分子动力学的发展, 进而实现了以较低计算资源获得更可靠的计算结果.Tungsten (W) alloys and W-based alloys are the primary candidate materials for plasma-facing components in future fusion reactors (e.g. ITER and CFETR). One of the critical issues still to be clarified in the design of the fusion reactor materials is the retention of hydrogen (H) isotopes in W, when the plasma-facing materials are supposed to sustain high-flux plasma and high-energy neutron. The dynamical behaviours of H in W with radiation defects (e.g. vacancy) are of serious concerns for understanding the mechanism of H capture, retention and permeation in W. In this work, a new model to extract the effective capture radius (ECR) and dissociation coefficient simultaneously is presented through coupling the trapping process and detrapping process of H in W vacancy. In the new model, the quantity ratio of vacancy to H atom in vacancy-H complex (VHx+1) in the molecular dynamics (MD) simulations is described as a function of time, while the exact occurrence time of corresponding event is not required. This new model, combined with extensive MD calculations, enables the simultaneous determining of the ECR and dissociation coefficient of H in W vacancy. It is found that the parameters are dependent not only on the event type but also on temperature. The dissociation energy of H from vacancy-H complex decreases gradually with the increase of the trapped number of H atoms in the vacancy-H complex. It is also found that the common assumption (i.e. the ECR is equal to one lattice constant and the pre-exponential factor is equal to 1013 s–1) in the long-term simulation methods (e.g. kinetic Monte Carlo and rate theory) is not always valid, since these calculated dynamical parameters are dispersive. The new model to obtain more reliable results with lower cost of computing resources can be easily extended into the other similar kinetic processes (e.g. H/He trapping and detrapping processes in other materials systems). These calculated dynamical parameters should be potentially helpful in supplying the initial input parameters for the long-term simulation methods.
-
Keywords:
- capture radius /
- dissociation coefficient /
- molecular dynamics /
- vacancy
1. 引 言
金属钨(W)以其具有优异性能(高熔点、高热导、高硬度、抗腐蚀和低溅射等)而被广泛应用于国防军事工业及核技术等领域, W基材料被认为是最具潜力的面对等离子体材料(plasma-facing material, PFM)[1,2]. 而作为PFM的W需承受等离子体辐照、高热负荷和中子辐照, 会引起复杂且不同尺度的损伤[3–6], 导致氢(H)同位素(包括氘(D)和氚(T))的滞留, 影响PFM的服役及燃料的自持, H同位素滞留问题是聚变堆第一壁材料设计的关键. W中H滞留受多种因素影响, 如: H的引入方式、温度、压力、材料的成分及结构等. 在聚变堆偏离器区域, W将受到高束流密度(1022—1024 m–2·s–1)的低能(通常 < 200 eV) D/T辐照[7], 该辐照也称为“亚临界”辐照, 即入射粒子能量小, 不能直接使W阵点离位, 且对应射程只有几纳米[8]. 然而H滞留深度可超过微米[9,10], 这是因为H在W晶格中扩散势垒低[11–13], 且不能通过自捕获机制形成H团簇, 因而可快速迁移至表面层(指射程区域)以下的扩散层. 然而H在W晶格中的平衡溶解度极低[12], 表明扩散层的H滞留主要受到材料本征缺陷或中子辐照诱导缺陷的影响. 缺陷对H的持续捕获可能导致H泡(或称H团簇)的形核及长大, 因而对空位[14–16]、位错[17,18]和晶界[19,20]等缺陷处H泡非均匀形核机制的研究至关重要. 虽然W中空位的平衡浓度低, 但是中子辐照引起离位损伤使材料中存在大量的过饱和空位[21], 这些空位的大小、分布及其演化直接影响W的性能及寿命, 因而深入理解H在空位处非均匀形核长大过程有助于揭示W中H起泡及滞留的机制.
W中H演化微观过程的实验研究极为困难, 需借助多尺度方法研究. 第一性原理[22]可计算小尺寸(指约100原子数的体系)缺陷结构的形成能, 也可选定部分路径计算迁移势垒. 分子动力学(molecular dynamics, MD)[23]可对大尺寸超胞(如百万原子数的体系)开展计算, 且能研究动力学事件(一般指小于纳秒的过程), 获取动力学参数. 动力学蒙特卡罗和速率理论可对更大尺寸体系研究更长时间的演化, 然而其结果的可靠性取决于输入参数的准确性. 这些输入参数包括热力学参数和动力学参数, 可通过实验或计算(第一性原理/MD)获得. 然而对于动力学参数, 实验上往往通过预设物理模型间接外推, 难以准确直接测量, 又由于温度的影响, 动力学转变机制复杂多变, 所以需借助MD开展研究.
MD经过数十年的发展, 模拟技术已比较成熟, 当前MD模拟发展趋势主要是以更高的效率、模拟更大的体系、实现更长的演化时间、取得更精确的结果[24]. 本课题组前期已自主开发了GPU并行加速的分子动力学软件(MDPSCU), 实现了MD模拟从计算技术和算法上的发展[25]; 动力学过程的模拟是MD模拟的特色, 而提取动力学参数的物理模型还不完善, 需进一步发展. 比如, 前期在计算有效捕获半径(ECR)时, 通过在一个超胞中引入多个H, 从而加速了超胞中空位捕获H事件的发生, 这也是提高MD模拟效率的重要方法[26]; 另外, 前期针对多路径转变(如双空位-H复合体(V2Hx+1)的H解离(V2Hx+1 → V2Hx + H)或空位解离(V2Hx+1 → VHx+1 + V), 计算各路径的动力学参数时, 重构了物理模型, 给出了新公式, 以获取更精确的计算结果[27]. 然而前期推导的物理模型[26–28]均需准确记录相应动力学事件首次发生的时间, 否则会发生计算精度不足的问题, 为了解决这一问题, 本文通过耦合捕获与解离过程, 推导出新物理模型, 并针对W中空位捕获H及其解离的微观动力学过程开展MD研究, 验证了新模型的可靠性, 且新模型可同时提取有效捕获半径和解离系数等动力学参数. 因而本文目的不仅是为动力学蒙特卡罗和速率理论长时空尺度方法提供输入参数, 同时也是要在方法上实现以较低计算资源获得可靠的物理参数, 这也是MD的重要发展方向.
2. 物理模型与计算细节
2.1 物理模型
针对W中空位捕获H和解离H的两个单一动力学过程, 前期分别推导了提取有效捕获半径和解离系数的物理模型[26]. 其中有效捕获半径用以表征空位及空位-H复合体(VHx,
x=0,1,⋯ )捕获溶质H原子的能力, 可等效为VHx与溶质H的反应距离[29], 因而捕获速率可表示为捕获H后的复合体(VHx+1)的浓度(CVHx+1 )随时间(t)的变化率, 即dCVHx+1dt=4πRc,x(DH+DVHx)CVHxCH, (1) 式中, Rc, x即VHx对溶质H的有效捕获半径; DA (A = H或VHx)表示A的扩散系数, 而溶质H的扩散系数一般要比W空位大得多(即DH
≫ DV)[13,26], 且空位-H复合体的可动性比空位更低[30], 因而DVHx 在计算过程中可忽略. 另外根据体积为V的超胞中A的数量(NA)与浓度(CA)的关系, NA = VCA, 则(1)式可变换为VdNVHx+1dt=4πRc,xDHNVHxNH. (2) 在MD模拟初始时, 超胞中只创建一个VHx和M个溶质H原子, 如前所述, 这里设置M个H是为了加速捕获事件的发生, 但高浓度的溶质H将影响其扩散行为[31], 因而M不宜太大, 前期研究[26]发现对于2000个W原子的体系, M取值为10时, 计算结果较为可靠, 受浓度影响可忽略. 对于捕获过程的MD模拟, 只关注捕获事件(VHx + H → VHx+1), 那么对于超胞, 要么未发生捕获(存在一个VHx), 要么已发生捕获(存在一个VHx+1), 对于已发生捕获的超胞将停止运行, 且认为该超胞中在后续时间都存在已捕获溶质H的VHx+1. 进行Nb次模拟, 即针对Nb个超胞开展MD模拟, 每个超胞设置不同随机数种子(溶质H的位置也随机设置). 假定Nt, x+1表示t时刻已发生捕获事件的超胞数, 可相当于含有VHx+1的数量, 则(2)式可变换为
VdNt,x+1dt=4πMRc,xDH(Nb−Nt,x+1). (3) 对(3)式进行积分, 可得
ln(Nb−Nt,x+1)=−4MπRc,xDHVt+ln(Nb−N0,x+1). (4) 因而VHx对溶质H的有效捕获半径(Rc, x)可通过ln(Nb–Nt, x+1)-t曲线的斜率推导获得. 类似上述推导过程, 可推导出空位-H复合体(VHx+1)解离H过程的物理模型[26–28], 假定该动力学过程为一级反应, 则
dCVHx+1dt=−kdiss,x+1CVHx+1, (5) 式中kdiss, x+1为解离系数, 表征上述事件发生的快慢. 同样进行Nb次模拟, 每个超胞初始时仅含有一个VHx+1, 这儿Nt, x+1可表示为t时刻未曾发生解离事件的超胞数, 则
lnNt,x+1N0,x+1=−kdiss,x+1t. (6) 因而VHx+1的解离系数(kdiss, x+1)可通过ln(Nt, x+1)-t曲线的斜率推导获得. 基于上述两单一过程的物理模型, 前期对W中空位/双空位-H复合体(V1/2Hx)捕获H及其解离动力学过程进行了研究[26,27].
上述物理模型均要求准确记录已发生相应事件的时间, 然而MD模拟过程中一般间隔一定时间(或步数)输出超胞的构型(包括原子坐标、速度和受力等信息), 通过分析构型中是否存在空位-H复合体(VHx+1)来判定相应事件是否已发生. 该处理方法可能存在误差, 比如, 对于捕获过程, 在某间隔时间内, 某超胞可能连续发生捕获和解离的事件, 从而在输出构型时, 并未鉴别出该超胞已发生过捕获事件. 解决方法之一是缩短输出构型的时间间隔, 然而该方法将增大数据的存储负担; 第二种方法是在线分析法, 该方法不需输出构型, 而是直接在MD模拟过程中一边计算一边分析, 然而该方法将增加计算时间, 尤其是对于GPU并行计算, 往往存在不同内存之间的数据拷贝, 影响并行效率, 另外该方法还不利于非预设的数据处理. 事实上最好的方法是通过改进物理模型来防止上述误差的产生. 记录的参数在新模型中应只与当前构型有关, 而不是还与之前的构型有关. 可耦合捕获与解离两基本动力学过程, 构建新模型以避免上述不足. 下面以解离过程为例推导, 假定MD模拟初始时, 超胞中只含有VHx+1, MD运行过程中存在H的解离, 同时也可能存在解离H的再捕获, 那么根据(1)和(5) 式可知
dCVHx+1dt=4πRc,xDHCVHxCH−kdiss,x+1CVHx+1. (7) (7)式已忽略了相比溶质H扩散系数(DH)较小的
DVHx+1 , 乘以超胞体积(V)可得dnVHx+1dt=4πRc,xDHVnVHxnH−kdiss,x+1nVHx+1, (8) 式中
nVHx 和nH特指单个超胞中VHx和溶质H的数量, 若进行Nb次独立的MD模拟, 则dNVHx+1dt=4πRc,xDHV(Nb−NVHx+1)−kdiss,x+1NVHx+1, (9) 式中
NVHx+1=∑nVHx+1 , 对(9)式进行数学转换和积分可得yVHx+1(t)=[yVHx+1(0)−q]exp(−pt)+q. (10) (10)式定义了
yVHx+1(t)=NVHx+1(t)Nb, q=4πRc,xDH4πRc,xDH+Vkdiss,x+1, p=4πRc,xDHV+kdiss,x+1. 因而基于(10)式, 相应过程的解离系数(甚至包括其逆过程的有效捕获半径)可通过拟合
yVHx+1(t) -t曲线推导获得, 这里VHx+1的比率yVHx+1(t)=NVHx+1(t)/Nb 仅与t时刻的输出构型有关, 因而新模型原则上可获得更精确的动力学参数. 根据y-t曲线拟合的参数p与q, 解离系数(kdiss,x+1)和有效捕获半径(Rc,x)可得Rc,x=qpV4πDH, kdiss,x+1=p(1−q). 2.2 计算细节
本文针对W中空位捕获H及其解离过程开展MD模拟, 并基于新物理模型拟合方程, 提取有效捕获半径和解离系数等动力学参数, 这儿选用Bonny等[32]拟合的W-H嵌入原子势函数(EAM1), 时间步长选用0.2 fs, 使用速度-Verlet算法[33]来求解原子的运动方程. W一般为体心立方(BCC)结构[34], 超胞大小为10a × 10a × 10a, 其中a为BCC-W的晶胞参数(a = 3.14 Å[35]), x, y, z取向分别为[100], [010]和[001], 均使用周期性边界条件, 完美超胞含有2000个W原子, 超胞尺寸足以消除点缺陷之间的相互作用, 可保证计算的可靠性. 对于解离过程, 在超胞中心删除一个W原子以创建一个空位, 然后在空位中随机插入若干个H原子, 以构建空位-H复合体(VHx+1), 如图1所示, VH2位于超胞的中心, 该图原子的显示使用了OVITO软件[36].
MD模拟的基本过程与前期研究类似[37,38]: 首先是将含有VHx+1的初始超胞进行原子位置局部最优化处理, 然后按麦克斯韦-玻尔兹曼分布设定原子速度热化超胞, 将超胞温度稳定至所需温度, 再进行22 ps的自由演化, 此后进行正式的模拟且采集有效的数据, 在该阶段采用微正则(NVE)系综, 运行100万步(即200 ps), 超胞的瞬间构型每隔2万步(即4 ps)输出一次. 在运行结束后, 对输出的构型进行原子位置局部最优化处理, 然后通过对比优化后的构型和完美点阵的Wigner-Seitz胞, 判定超胞中是否存在空位-H复合体(VHx+1). 对于捕获过程, 主要的不同是初始超胞的构建, 除了设置VHx, 还随机设置10个溶质H原子.
3. 结果与讨论
3.1 解 离
对于解离过程, 初始超胞中只设置空位-H复合体(VHx+1), 在MD模拟过程, 会出现H的解离, 当然也可能会出现解离后的溶质H被VHx重新捕获的情况. 但根据新物理模型, 只要知道当前构型下空位-H复合体(VHx+1)的数量即可推导出相应过程的动力学参数, 而不必准确记录相应事件的首次发生时间. 图2所示是不同类型的事件在不同温度下空位-H复合体(VHx+1)的比率(
yVHx+1(t) )随时间(t)的变化, 由图2可以发现MD模拟获得各输出时间点的VHx+1的比率(图2中点表示)能被新物理模型(图2中曲线表示)较好地拟合, 表明该模型可用于描述捕获和解离的耦合过程.图2(a)—(f)主要用于研究解离过程, 即初始时超胞中只设置了空位-H复合体(VHx+1). 根据拟合参数p与q以及(12b)式可获得解离系数(kdiss, x+1). 另外若解离是热激活过程, 则解离系数可用Arrhenius关系来描述温度(T )的影响, 即
kdiss,x+1=vdiss,x+1exp(−Ediss,x+1kBT), (13) 式中, νdiss,x+1是前因子, 一般认为与温度无关; Ediss,x+1是解离能, 表示空位-H复合体(VHx+1)的H解离形成溶质H的势垒; kB是玻尔兹曼常数. 如图3所示, lnkdiss,x+1与1/T呈现较好的线性关系, 表明解离过程是个热激活过程, 因而可根据曲线截距和斜率推导出前因子和解离能.
新物理模型推导出的解离能在0.58—1.9 eV之间, 与第一性原理计算结果[30]接近. 如图4(a)所示, 新模型推导的解离能(Enew)与前期单一过程物理模型给出的解离能(Eold) [26]比较接近, 表明两种模型在计算解离能时都比较合理, 另外从图4(a)还可发现空位-H复合体(VHx+1)的H解离能随捕获H的数量增加而减小, 该趋势也与前期结果[26]及第一性原理计算结果[30]一致. 图4(b)是解离前因子的比较, 可以发现二者的接近程度不如解离能, 这主要是由于经过了指数的数学变换. 事实上对lnν, 前后模型给出的数值也比较接近. 另外由图4(b)还可以发现, 前因子粗略地随捕获H数量的增加而增加, 即粗略地随解离能的减小而增加. 前因子的变化范围在6—70 ps–1, 表明在一些计算中前系数被认为是常数(1013 s–1)的假定是不合理的.
3.2 捕 获
如前文所述, 在研究空位-H复合体(VHx+1)解离过程时也可能存在解离H再被捕获的情况, 而且根据(12a)式, 有效捕获半径(Rc,x)理应也能推导获得, 事实上通过解离过程拟合获得的Rc,x大多在0.5—4 Å, 与前期[26]结果相近. 然而应指出的是, 对于有效捕获半径的计算, MD模拟过程中应该保证有足够数量的捕获事件发生以满足统计学规律. 如图2(a)所示的2250 K时VH → V+H事件, 含有VH超胞的比率 (yVH)几乎随时间(t)线性下降, 表明存在较少的捕获事件, 因而该y-t曲线不宜用来推导有效捕获半径, 类似的情况还有1500 K时的VH2 → VH+H事件(图2(b))、1000 K时的VH3 →VH2+H事件(图2(c))、1000 K时的VH4 → VH3+H事件(图2(d))、500 K时的VH5 →VH4+H事件(图2(e))和700 K时的VH6 → VH5 +H事件(图2(f))等. 另外对于VH6 → VH5+H事件, 由于解离能较低, 不易发生捕获事件, 因而对于800 K和900 K也不宜用于计算有效捕获半径. 事实上, 前期事实上, 前期[26]在研究高捕获状态的空位-H复合体(如VH7)的解离过程时就发现, 在热化阶段绝大多数超胞中的(VHx+1)已发生解离, 因而在后期MD演化过程中很难发生捕获事件.
因而针对捕获过程, 在超胞中设置了VHx和溶质H原子, 增加捕获事件的发生. 同样在MD模拟过程中, 溶质H原子可能会被VHx捕获, 捕获的H也可能再解离, 那么空位-H复合体的比率(
yVHx+1(t) )随时间(t)的变化关系仍可用(10)式来描述. 图2(g)表示的是V+H → VH事件, 图2(h)表示的是VH+H → VH2事件, 可以发现MD模拟获得比率y(t)能被新物理模型较好地拟合, 同样表明新模型的可靠性, 根据拟合参数且运用(12a)式可推导出有效捕获半径. 计算结果如图5所示, 新模型获得的有效捕获半径与单一过程物理模型[26]给出的结果比较接近. 偏离较大的主要是高捕获状态的空位-H复合体(VHx), 这是因为这些复合体的解离能较低且存在较复杂的捕获与解离过程. 另外由图5还可发现有效捕获半径在0.5—4 Å, 表明在一些计算中有效捕获半径被认为是一个晶格常数(3.14 Å)的假定可能并不合理.4. 结 论
由于单一过程提取动力学参数的物理模型需准确记录相应动力学事件首次发生的时间, 为了弥补这一不足, 本文通过耦合捕获与解离过程, 推导了新物理模型, 新模型通过不同时刻的瞬间构型, 获得空位-H复合体(VHx+1)的比率与时间的关系, 即可同时提取有效捕获半径和解离系数等动力学参数. 针对W空位捕获H及其解离的微观动力学过程开展了MD研究, 验证了新模型的可靠性. 提取了解离系数和有效捕获半径等动力学参数, 发现空位-H复合体(VHx+1)的H解离能随捕获H数量的增加而减小, 且解离前因子及有效捕获半径的数值较分散, 不应假定为常数. 本文获取的动力学参数可为动力学蒙特卡罗和速率理论等长时空尺度方法提供输入参数, 同时本文推导的物理模型可实现以较低计算资源获得可靠的物理参数, 这也是MD的重要发展方向.
[1] Lu G H, Zhou H B, Becquart C S 2014 Nucl. Fusion 54 086001
Google Scholar
[2] 李建刚, 宋云涛, 刘永, 汪小琳, 万元熙 2016 聚变工程试验堆装置主机设计 (北京: 科学出版社) 第250−252页
Li J G, Song Y T, Liu Y, Wang X L, Wan Y X 2016 Main Machine Design of Fusion Engineering Test Reactor (Beijing: Science Press) pp250−252 (in Chinese)
[3] Jia Y Z, Liu W, Xu B, Luo G N, Li C, Fu B Q, de Temmerman G 2015 J. Nucl. Mater. 463 312
Google Scholar
[4] Yi X, Jenkins M L, Hattar K, Edmondson P D, Roberts S G 2015 Acta Mater. 92 163
Google Scholar
[5] Hu X, Koyanagi T, Fukuda M, Kumar N A P K, Snead L L, Wirth B D, Katoh Y 2016 J. Nucl. Mater. 480 235
Google Scholar
[6] Fu B Q, Lai W S, Yuan Y, Xu H Y, Li C, Jia Y Z, Liu W 2013 Chinese Phys. B 22 126601
Google Scholar
[7] Ye M Y 2005 Plasma Sci. Technol. 7 2828
Google Scholar
[8] Fu B Q, Qiu M J, Zhai L, Yang A L, Hou Q 2019 Nucl. Instrum. Meth. B 450 220
Google Scholar
[9] Oyaidzu M, Hayashi T, Alimov V K 2016 Nucl. Mater. Energy 9 93
Google Scholar
[10] Jia Y Z, Liu W, Xu B, Qu S L, Shi L Q, Morgan T W 2017 Nucl. Fusion 57 034003
Google Scholar
[11] Heinola K, Ahlgren T 2010 Phys. Rev. B 81 073409
Google Scholar
[12] Frauenfelder R 1969 J. Vac. Sci. Technol. 6 388
Google Scholar
[13] Qiu M J, Zhai L, Cui J C, Fu B Q, Li M, Hou Q 2018 Chin. Phys. B 27 073103
Google Scholar
[14] Liu Y L, Zhang Y, Luo G N, Lu G H 2009 J. Nucl. Mater. 390/391 1032
Google Scholar
[15] Jiang B, Wan F R, Geng W T 2010 Phys. Rev. B 81 134112
[16] Liu Y L, Gao A Y, Lu W, Zhou H B, Zhang Y 2012 Chin. Phys. Lett. 29 077101
[17] Guo W, Ge L, Yuan Y, Cheng L, Wang S, Zhang X, Lu G H 2019 Nucl. Fusion 59 026005
Google Scholar
[18] Terentyev D, Dubinko V, Bakaev A, Zayachuk Y, van Renterghem W, Grigorev P 2014 Nucl. Fusion 54 042004
Google Scholar
[19] Zhou H B, Liu Y L, Jin S, Zhang Y, Luo G N, Lu G H 2010 Nucl. Fusion 50 025016
Google Scholar
[20] Wang X X, Niu L L, Wang S 2017 J. Nucl. Mater. 487 158
Google Scholar
[21] Fu B Q, Xu B, Lai W S, Yuan Y, Xu H Y, Li C, Jia Y Z, Liu W 2013 J. Nucl. Mater. 441 24
Google Scholar
[22] 王欣欣, 张颖, 周洪波, 王金龙 2014 物理学报 63 046103
Google Scholar
Wang X X, Zhang Y, Zhou H B, Wang J L 2014 Acta Phys. Sin. 63 046103
Google Scholar
[23] 汪俊, 张宝玲, 周宇璐, 侯氢 2011 物理学报 60 106601
Google Scholar
Wang J, Zhang B L, Zhou Y L, Hou Q 2011 Acta Phys. Sin. 60 106601
Google Scholar
[24] 严六明, 朱素华 2013 分子动力学模拟的理论与实践 (北京: 科学出版社) 第6−12页
Yan L M, Zhu S H 2013 Theory and Practice of Molecular Dynamics Simulation (Beijing: Science Press) pp6−12 (in Chinese)
[25] Hou Q, Li M, Zhou Y, Cui J, Cui Z, Wang J 2013 Comput. Phys. Commun. 184 2091
Google Scholar
[26] Fu B Q, Qiu M J, Cui J C, Li M, Hou Q 2018 J. Nucl. Mater. 508 278
Google Scholar
[27] Fu B Q, Qiu M J, Cui J C, Li M, Wang J, Hou Q 2019 Nucl. Instrum. Meth. B 452 21
Google Scholar
[28] Zhou Y L, Wang J, Hou Q, Deng A H 2014 J. Nucl. Mater. 446 49
Google Scholar
[29] Chandrasekhar S 1943 Rev. Mod. Phys. 15 1
Google Scholar
[30] Fernandez N, Ferro Y, Kato D 2015 Acta Mater. 94 307
Google Scholar
[31] Ahlgren T, Bukonte L 2016 J. Nucl. Mater. 479 195
Google Scholar
[32] Bonny G, Grigorev P, Terentyev D 2014 J. Phys.-Condens. Mat. 26 485001
Google Scholar
[33] Swope W C, Andersen H C, Berens P H, Wilson K R 1982 J. Chem. Phys. 76 637
Google Scholar
[34] Fu B Q, Liu W, Li Z L 2009 Appl. Surf. Sci. 255 8511
Google Scholar
[35] Marinica M C, Ventelon L, Gilbert M R, Proville L, Dudarev S L, Marian J, Bencteux G, Willaime F 2013 J. Phys.-Condens. Mat. 25 395502
Google Scholar
[36] Stukowski A 2010 Model. Simul. Mater. Sci. 18 015012
Google Scholar
[37] Cheng G J, Fu B Q, Hou Q, Zhou X S, Wang J 2016 Chinese Phys. B 25 076602
Google Scholar
[38] Fu B Q, Fitzgerald S P, Hou Q, Wang J, Li M 2017 Nucl. Instrum. Meth. B 393 169
Google Scholar
-
-
[1] Lu G H, Zhou H B, Becquart C S 2014 Nucl. Fusion 54 086001
Google Scholar
[2] 李建刚, 宋云涛, 刘永, 汪小琳, 万元熙 2016 聚变工程试验堆装置主机设计 (北京: 科学出版社) 第250−252页
Li J G, Song Y T, Liu Y, Wang X L, Wan Y X 2016 Main Machine Design of Fusion Engineering Test Reactor (Beijing: Science Press) pp250−252 (in Chinese)
[3] Jia Y Z, Liu W, Xu B, Luo G N, Li C, Fu B Q, de Temmerman G 2015 J. Nucl. Mater. 463 312
Google Scholar
[4] Yi X, Jenkins M L, Hattar K, Edmondson P D, Roberts S G 2015 Acta Mater. 92 163
Google Scholar
[5] Hu X, Koyanagi T, Fukuda M, Kumar N A P K, Snead L L, Wirth B D, Katoh Y 2016 J. Nucl. Mater. 480 235
Google Scholar
[6] Fu B Q, Lai W S, Yuan Y, Xu H Y, Li C, Jia Y Z, Liu W 2013 Chinese Phys. B 22 126601
Google Scholar
[7] Ye M Y 2005 Plasma Sci. Technol. 7 2828
Google Scholar
[8] Fu B Q, Qiu M J, Zhai L, Yang A L, Hou Q 2019 Nucl. Instrum. Meth. B 450 220
Google Scholar
[9] Oyaidzu M, Hayashi T, Alimov V K 2016 Nucl. Mater. Energy 9 93
Google Scholar
[10] Jia Y Z, Liu W, Xu B, Qu S L, Shi L Q, Morgan T W 2017 Nucl. Fusion 57 034003
Google Scholar
[11] Heinola K, Ahlgren T 2010 Phys. Rev. B 81 073409
Google Scholar
[12] Frauenfelder R 1969 J. Vac. Sci. Technol. 6 388
Google Scholar
[13] Qiu M J, Zhai L, Cui J C, Fu B Q, Li M, Hou Q 2018 Chin. Phys. B 27 073103
Google Scholar
[14] Liu Y L, Zhang Y, Luo G N, Lu G H 2009 J. Nucl. Mater. 390/391 1032
Google Scholar
[15] Jiang B, Wan F R, Geng W T 2010 Phys. Rev. B 81 134112
[16] Liu Y L, Gao A Y, Lu W, Zhou H B, Zhang Y 2012 Chin. Phys. Lett. 29 077101
[17] Guo W, Ge L, Yuan Y, Cheng L, Wang S, Zhang X, Lu G H 2019 Nucl. Fusion 59 026005
Google Scholar
[18] Terentyev D, Dubinko V, Bakaev A, Zayachuk Y, van Renterghem W, Grigorev P 2014 Nucl. Fusion 54 042004
Google Scholar
[19] Zhou H B, Liu Y L, Jin S, Zhang Y, Luo G N, Lu G H 2010 Nucl. Fusion 50 025016
Google Scholar
[20] Wang X X, Niu L L, Wang S 2017 J. Nucl. Mater. 487 158
Google Scholar
[21] Fu B Q, Xu B, Lai W S, Yuan Y, Xu H Y, Li C, Jia Y Z, Liu W 2013 J. Nucl. Mater. 441 24
Google Scholar
[22] 王欣欣, 张颖, 周洪波, 王金龙 2014 物理学报 63 046103
Google Scholar
Wang X X, Zhang Y, Zhou H B, Wang J L 2014 Acta Phys. Sin. 63 046103
Google Scholar
[23] 汪俊, 张宝玲, 周宇璐, 侯氢 2011 物理学报 60 106601
Google Scholar
Wang J, Zhang B L, Zhou Y L, Hou Q 2011 Acta Phys. Sin. 60 106601
Google Scholar
[24] 严六明, 朱素华 2013 分子动力学模拟的理论与实践 (北京: 科学出版社) 第6−12页
Yan L M, Zhu S H 2013 Theory and Practice of Molecular Dynamics Simulation (Beijing: Science Press) pp6−12 (in Chinese)
[25] Hou Q, Li M, Zhou Y, Cui J, Cui Z, Wang J 2013 Comput. Phys. Commun. 184 2091
Google Scholar
[26] Fu B Q, Qiu M J, Cui J C, Li M, Hou Q 2018 J. Nucl. Mater. 508 278
Google Scholar
[27] Fu B Q, Qiu M J, Cui J C, Li M, Wang J, Hou Q 2019 Nucl. Instrum. Meth. B 452 21
Google Scholar
[28] Zhou Y L, Wang J, Hou Q, Deng A H 2014 J. Nucl. Mater. 446 49
Google Scholar
[29] Chandrasekhar S 1943 Rev. Mod. Phys. 15 1
Google Scholar
[30] Fernandez N, Ferro Y, Kato D 2015 Acta Mater. 94 307
Google Scholar
[31] Ahlgren T, Bukonte L 2016 J. Nucl. Mater. 479 195
Google Scholar
[32] Bonny G, Grigorev P, Terentyev D 2014 J. Phys.-Condens. Mat. 26 485001
Google Scholar
[33] Swope W C, Andersen H C, Berens P H, Wilson K R 1982 J. Chem. Phys. 76 637
Google Scholar
[34] Fu B Q, Liu W, Li Z L 2009 Appl. Surf. Sci. 255 8511
Google Scholar
[35] Marinica M C, Ventelon L, Gilbert M R, Proville L, Dudarev S L, Marian J, Bencteux G, Willaime F 2013 J. Phys.-Condens. Mat. 25 395502
Google Scholar
[36] Stukowski A 2010 Model. Simul. Mater. Sci. 18 015012
Google Scholar
[37] Cheng G J, Fu B Q, Hou Q, Zhou X S, Wang J 2016 Chinese Phys. B 25 076602
Google Scholar
[38] Fu B Q, Fitzgerald S P, Hou Q, Wang J, Li M 2017 Nucl. Instrum. Meth. B 393 169
Google Scholar
期刊类型引用(2)
1. 包宏伟,李燕,马飞. 金属钨辐照缺陷与氢/氦作用的计算模拟研究进展. 稀有金属材料与工程. 2022(03): 1100-1110 . 百度学术
2. 董向成,王国伟,陈建宏. 闪电等离子体连续光谱轮廓形成机理研究(英文). 中国光学. 2021(05): 1243-1250 . 百度学术
其他类型引用(0)
计量
- 文章访问数: 9433
- PDF下载量: 122
- 被引次数: 2