H与C/Be混合层相互作用的分子动力学模拟

核聚变与等离子体物理 ›› 2014, Vol. 34 ›› Issue (2): 111-117.

H与C/Be混合层相互作用的分子动力学模拟

曹小岗1，田树平1，欧巍1，王建强1，张静全2，潘宇东3，苟富均1，陈顺礼1

(1. 四川大学原子核科学技术研究所，成都 610064； 2. 四川大学材料科学与工程学院，成都 610064； 3. 核工业西南物理研究院，成都 610041)

收稿日期:2013-07-07 修回日期:2014-04-03 出版日期:2014-06-15 发布日期:2014-06-16
作者简介:曹小岗(1989-)，男，江西余干人，硕士研究生，从事等离子体与材料表面相互作用研究。
基金资助:
国家磁约束核聚变能发展研究专项(2013GB114003)；科技部863基金资助项目(2011AA050515)；国家自然科学基金资助项目(11275135)

Molecular dynamic simulation of H interaction with C/Be surface

CAO Xiao-gang1, TIAN Shu-ping1, OU Wei1, WANG Jian-qiang1,ZHANG Jing-quan2, PAN Yu-dong3, GOU Fu-jun, CHEN Shun-li1

(1. Institute of Nuclear Science and Technology Sichuan University, Chengdu 610064; 2. College of Materials Science and Engineering, Sichuan University, Chengdu, 610064; 3. Southwestern Institute of Physics, Chengdu 610041)

Received:2013-07-07 Revised:2014-04-03 Online:2014-06-15 Published:2014-06-16

摘要/Abstract

摘要：

利用分子动力学方法研究了H原子与C/Be样品的相互作用过程，当H原子轰击C/Be样品时，发现有一些H原子渗入样品

中并且滞留在样品中，H原子的滞留率随H原子的初始入射能量的升高呈线性增长,有些沉积在样品中H原子与C原子相

互作用形成H−C键。溅射产物以H原子和H₂分子为主。H和H₂的产额率随初始入射能量的变化趋势相反，分析了不同机

制下产物H和H₂的产额率随初始入射能量的关系，且通过分析H原子的入射能量和样品的原子密度的关系来研究轰击后

的样品，发现样品中原子分布变化很小，同时分析了化合物中的化学键分布变化较小，只是其化学键的分布峰向样品

表面移动。

关键词: 等离子体, 分子动力学, 聚变材料, 沉积, 化学键

Abstract:

　　Molecular dynamics method was used to study H interaction with C/Be samples in this paper. It is found that with H atoms bombarding C/Be sample, some H atoms penetrate into the bulk and deposit in the samples. The rate of H retention increases linearly with increasing incident energy. H atoms deposited in the sample react with C atoms forming H-C bonds. The main sputtering product is H atoms and H₂ molecules. With the initial incident energy increase, the trend about yield rate of H and H₂ is opposite. The relationship between yield rate of H and H₂ in different mechanisms and the initial incident energy is analyzed. The modified samples were studied by analyzing incident energy and atom density. The distribution of chemical bonds in the sample was examined. It is found that the change of the distribution of the atoms in the sample is small, the peak of the distribution of chemical bonds only moves to the surface of the sample.

Key words: Plasma, Molecular dynamics, Fusion material, Deposition, Chemistry bond

中图分类号:

O561.4

曹小岗1，田树平1，欧巍1，王建强1，张静全2，潘宇东3，苟富均1，陈顺礼1. H与C/Be混合层相互作用的分子动力学模拟[J]. 核聚变与等离子体物理, 2014, 34(2): 111-117.

CAO Xiao-gang1, TIAN Shu-ping1, OU Wei1, WANG Jian-qiang1,ZHANG Jing-quan2, PAN Yu-dong3, GOU Fu-jun, CHEN Shun-li1. Molecular dynamic simulation of H interaction with C/Be surface[J]. NUCLEAR FUSION AND PLASMA PHYSICS, 2014, 34(2): 111-117.

参考文献

[1]Shibata K. New observational facts about solar flares from Yohkoh studies-evidence of magnetic reconnection and a unified model of flares [J]. Adv. Space Res, 1996, 17(1): 9-18.

[2]Baker D N, Pulkkinen T I, Angelopoulos V, et al. Neutral line model of substorms: past results and present view [J]. J. Geophys. Res., 1996, 101(A6): 12975-13010.

[3]Waelbroeck F L. Current sheets and nonlinear growth of the m=1 kink-tearing mode [J]. Phys. Fluids, 1989, B 1(12): 2372-2380.

[4]Furth P H, Killeen J, Rosenbluth M N. Finite-resistivity instabilities of a sheet pinch [J]. Phys. Fluids, 1963, 6(3): 459-484.

[5]Chacon L, Knoll D A. A 2D high-β Hall MHD implicit nonlinear solver [J]. J. Comput. Phys., 2003, 188(4): 573-592.

[6]Fitzpatrick R. Scaling of forced magnetic reconnection in the Hall-magnetohydrodynamic Taylor problem [J]. Phys. Plasmas, 2004, 11(3): 937-946.

[7]He Zhixiong, Dong J Q, Long Y X, et al. Double tearing mode induced by parallel electron viscosity in tokamak plasmas [J]. Phys. Plasmas, 2010, 17(11): 112102.

[8]Fitzpatrick R, Waelbroeck F L. Two-fluid magnetic island dynamics in slab geometry [J]. Phys. Plasmas, 2005, 12(2): 022307.

[9]Bian N, Vekstein G. On the two-fluid modification of the resistive tearing instability [J]. Phys. Plasmas, 2007, 14(7): 072107.

[10]Chacón L, Knoll D A, Finn J M. An implicit nonlinear reduced resistive MHD solver [J]. J. Comput. Phys., 2002, 178(1): 15-36.

[11]Yabe T, Aoki T. A universal solver for hyperbolic equations by cubic-polynomial interpolation I. One- dimensional solver [J]. Comput. Phys. Commun., 1991, 66(2-3): 219-232.

[12]Yabe T, Ishikawa T, Wang P Y, et al. A universal solver for hyperbolic equations by cubic-polynomial inter- polation II. Two- and three-dimensional solvers [J]. Comput. Phys. Commun., 1991, 66(2?3): 233-242.

[13]Knoll D A, Chacón L. Magnetic reconnection in the two-dimensional Kelvin-Helmholtz instability [J]. Phys. Rev. Lett., 2002, 88(21): 215003.

[1]	何宗钰慧, 高金明, 郝广周, 梁绍勇, 才来中. HL-3装置Bolometer系统金属膜电阻探测器的研制[J]. 核工业西南物理研究院, 2025, 45(1): 13-18.
[2]	张金恒, 苌磊, 杨鑫, 周海山, 罗广南. 外界磁场和天线构型对射频等离子体微放电的影响研究 [J]. 核工业西南物理研究院, 2025, 45(1): 117-124.
[3]	李梦超, 郭勇, 荣垂才, 陈维, 黄骏, 王兴权. 时间分步平均法求解 1D LPIC 中的电流密度[J]. 核工业西南物理研究院, 2024, 44(4): 491-496.
[4]	杨秋磊, 陈伟, 张洁, 程时葵, 梁绍勇, 张轶泼, HL- 实验团队. HL-3 装置欧姆等离子体中高频磁流体不稳定性的实验观测[J]. 核工业西南物理研究院, 2024, 44(3): 328-333.
[5]	孙海岗, 靳琛垚, 张炜, 江堤, 叶孜崇. 鞘层电容对快速扫描双探针伏安曲线的影响[J]. 核工业西南物理研究院, 2024, 44(2): 198-205.
[6]	陈明, 沈飊, Shinichiro Kado, 陈大龙, 张恒, 郭笔豪, 黄耀, 蔡福瑞, 肖炳甲. EAST等离子体边界可见光成像系统研究[J]. 核工业西南物理研究院, 2024, 44(2): 222-228.
[7]	孟凡卫, 张小辉, 胡广海, 孔德峰, 叶扬, 黄艳清, 齐美彬, 赵志豪, 李博, 董期龙. EAST-CTI测试平台上六硼化镧等离子体源的研制[J]. 核工业西南物理研究院, 2024, 44(1): 24-29.
[8]	高泽石, 王亚磊, 李彦龙, 田文喜, 才来中, 吴雪科, 连强, 李昕泽, 王占辉. 径向湍性输运条件下CFETR平行热通量及包层能量沉积模拟[J]. 核工业西南物理研究院, 2024, 44(1): 65-71.
[9]	鲍娜娜, 黄耀, 闫星廷, 魏世闻, 王梓昊. 基于Python的等离子体平衡重建程序设计及实现[J]. 核工业西南物理研究院, 2024, 44(1): 98-104.
[10]	李宜轩, 夏凡, 杨宗谕, 龚新文, 陈程远, 肖国梁, 朱晓博 . 基于深度学习的 HL-2A 装置 ELM 实时控制系统的研制[J]. 核工业西南物理研究院, 2023, 43(4): 380-385.
[11]	谢亚红, 谢俊炜, 韦江龙, 顾玉明, 蒋才超, 潘军军, 许永建, 谢远来, 胡纯栋. 大面积负离子束源磁过滤器的设计及实验研究 [J]. 核工业西南物理研究院, 2023, 43(4): 407-411.
[12]	何敏, 于子泉, 贾忠琪, 孟繁星, 朱超. 基于 PWM 整流器的大功率等离子体点火电源 [J]. 核工业西南物理研究院, 2023, 43(4): 418-423.
[13]	马铭键, 叶孜崇(Yip Chi-shung), 江堤, 毛世峰, 张炜, 靳琛垚, 叶民友 . 激光诱导荧光诊断的波长标定方法的研究 [J]. 核工业西南物理研究院, 2023, 43(4): 469-476.
[14]	袁强华, 刘珊珊, 殷桂琴, 秦彪. 不同频率驱动下容性耦合Ar等离子体随气压变化的放电特性[J]. 核工业西南物理研究院, 2023, 43(4): 482-488.
[15]	何敏, 朱超, 贾忠琪, 孟繁星, 于子泉. 一种大功率等离子体火炬电源系统及控制策略[J]. 核工业西南物理研究院, 2023, 43(3): 312-316.