计算物理 ›› 2022, Vol. 39 ›› Issue (1): 101-108.DOI: 10.19596/j.cnki.1001-246x.8328
高慧1(), 杨在发1, 赵敬芬1, 袁慧敏1, 刘志娥1, 赵显2
收稿日期:
2021-01-04
出版日期:
2022-01-25
发布日期:
2022-09-03
作者简介:
高慧(1979-), 博士, 副教授, 研究方向为金属晶面腐蚀, E-mail: ghbnu0205@163.com
基金资助:
Hui GAO1(), Zaifa YANG1, Jingfen ZHAO1, Huimin YUAN1, Zhie LIU1, Xian ZHAO2
Received:
2021-01-04
Online:
2022-01-25
Published:
2022-09-03
摘要:
采用密度泛函理论研究Nb、Sn、Cu、Fe和Cr 5种合金元素对氧在Zr(0001)晶面吸附能的影响, 发现Nb、Sn和Cu会促进氧在Zr(0001)晶面吸附, Fe和Cr对氧吸附的影响因吸附的位置不同而存在差异。研究5种合金元素对氧化后Zr(0001)晶面化学键的破坏情况, 发现Nb对Zr(0001)晶面破坏程度最小而且可以迅速复原, Sn可以使(0001)晶面相邻的两个化学键都变长, 对晶面的破坏性大。最后讨论5种合金元素在Zr(0001)晶面的偏聚能, 发现Sn、Fe和Cr偏聚能为负值, 容易偏聚到Zr(0001)晶面, 而Nb和Cu偏聚能为正值, 不易在Zr(0001)晶面偏聚。综合以上分析, Nb可以促进氧在Zr(0001)晶面的吸附行为, 氧化后(0001)晶面可以快速复原, 从而阻碍其他氧原子进入, 抑制疖状腐蚀的发生。Sn容易偏聚到Zr(0001)晶面, 可以促进氧在Zr(0001)晶面吸附, 氧化后会造成Zr(0001)晶面较大的破坏, 促进氧进入Zr(0001)晶面, 促进疖状腐蚀的发生。
高慧, 杨在发, 赵敬芬, 袁慧敏, 刘志娥, 赵显. 第一性原理研究Nb、Sn、Cu、Fe和Cr对Zr (0001)晶面抗疖状腐蚀性能的影响[J]. 计算物理, 2022, 39(1): 101-108.
Hui GAO, Zaifa YANG, Jingfen ZHAO, Huimin YUAN, Zhie LIU, Xian ZHAO. Effect of Nb, Sn, Cu, Fe and Cr on Zr (0001) Surface Nodular Corrosion Resistance: First Principles Study[J]. Chinese Journal of Computational Physics, 2022, 39(1): 101-108.
图1 氧原子在Zr(0001)晶面吸附俯视图(a)氧原子吸附在top位;(b)氧原子吸附在hcp位;(c)氧原子吸附在fcc位 (紫色表示表层(第一层)Zr原子,黄色表示第二层Zr原子。红色的表示不同吸附位的氧原子,绿色的表示合金元素。黑色线表示超晶胞边界,圆棍表示原子间的化学键。)
Fig.1 Top view of atmoic structure of Zr (0001) crystal face (a) Oxygen atom is at the top position; (b) Oxygen atom is at the hcp position; (c) Oxygen atom is at the fcc position (The purple and yellow spheres represent Zr atoms in the first layer and the second layer, respectively. The red spheres represent oxygen atoms and the green ones represent alloying elements. The black line represents boundary of supercell and the sticks indicate bonds between atoms.)
position | pure-Zr | Nb | Sn | Cu | Fe | Cr |
top | -7.46 | -8.02 | -3.22 | -4.54 | -5.50 | -6.54 |
hcp | -10.75 | -12.80 | -12.21 | -12.32 | -8.31 | -11.69 |
fcc | -11.13 | -13.18 | -11.94 | -12.27 | -12.07 | -9.14 |
表1 Zr(0001)晶面和含Nb、Sn、Cu、Fe、Cr合金元素的Zr(0001)晶面的氧原子吸附能(Eb/eV)
Table 1 Adsorption energies of oxygen on Zr(0001) crystal face and Zr(0001) crystal face with Nb、Sn、Cu、Fe、Cr(Eb/eV)
position | pure-Zr | Nb | Sn | Cu | Fe | Cr |
top | -7.46 | -8.02 | -3.22 | -4.54 | -5.50 | -6.54 |
hcp | -10.75 | -12.80 | -12.21 | -12.32 | -8.31 | -11.69 |
fcc | -11.13 | -13.18 | -11.94 | -12.27 | -12.07 | -9.14 |
Cu | Sn | Nb | Fe | Cr | |
The first layer | Cu-Zr; Zr-Zr | Sn-Zr; Zr-Zr | Nb-Zr; Zr-Zr | Fe-Zr; Zr-Zr | Cr-Zr; Zr-Zr |
0.311;0.327 | 0.334;0.326 | 0.320;0.322 | 0.392;0.306 | 0.321;0.313 | |
The second layer | Zr-Zr | ||||
0.323;0.324 | 0.325;0.326 | 0.324;0.321 | 0.325;0.322 | 0.331;0.322 | |
The third layer | Zr-Zr | ||||
0.323;0.325 | 0.326;0.325 | 0.324;0.324 | 0.318;0.320 | 0.325;0.328 |
表2 含合金元素氧化后Zr(0001)晶面(三层)化学键键长(单位:nm)
Table 2 Bond lengths of oxidation Zr(0001) crystal face containing alloying elements (in the top three layers) (in nm)
Cu | Sn | Nb | Fe | Cr | |
The first layer | Cu-Zr; Zr-Zr | Sn-Zr; Zr-Zr | Nb-Zr; Zr-Zr | Fe-Zr; Zr-Zr | Cr-Zr; Zr-Zr |
0.311;0.327 | 0.334;0.326 | 0.320;0.322 | 0.392;0.306 | 0.321;0.313 | |
The second layer | Zr-Zr | ||||
0.323;0.324 | 0.325;0.326 | 0.324;0.321 | 0.325;0.322 | 0.331;0.322 | |
The third layer | Zr-Zr | ||||
0.323;0.325 | 0.326;0.325 | 0.324;0.324 | 0.318;0.320 | 0.325;0.328 |
图2 含合金元素氧化后Zr(0001)晶面三层内最近邻两个化学键(X-Zr (X=Nb、Sn、Cu、Fe、Cr), Zr-Zr)与完整的晶体两个锆之间的键长差值(单位为nm)
Fig.2 Difference of bond length between the two nearest chemical bonds (X-Zr (X=Nb, Sn, Cu, Fe, Cr), Zr-Zr) and the two zirconium in a perfect crystal after the oxidation of Zr(0001)crystal surface (unit in nm)
Cu | Sn | Nb | Fe | Cr |
1.69 | -0.87 | 0.16 | -0.78 | -0.32 |
表3 Nb、Sn、Cu、Fe、Cr在Zr(0001)晶面的偏聚能(Es/eV)
Table 3 Partial poly energy of Nb, Sn, Cu, Fe and Cr on the surface of Zr(0001)(Es/eV)
Cu | Sn | Nb | Fe | Cr |
1.69 | -0.87 | 0.16 | -0.78 | -0.32 |
1 |
|
2 |
NIKULINA A V, MARKELOV V A. Zirconium alloy E635 as a material for fuel rod cladding and other components of VVER and RBMK cores[C]. Proceedings of the 11th International Symposium on Zirconium in the Nuclear Industry, Germany: ASTM International, 1995.
|
3 |
|
4 |
|
5 |
GEORGE P, SABOL G P, GERALD R, et al. Development of a cladding alloy for high burnup[C]. Zirconium in the Nuclear Industry: 8th International Symposium, San Diego: ASTM International, 1989.
|
6 |
|
7 |
DOI |
8 |
|
9 |
|
10 |
|
11 |
|
12 |
DOI |
13 |
DOI |
14 |
|
15 |
|
16 |
|
17 |
|
18 |
DOI |
19 |
DOI |
20 |
DOI |
21 |
|
22 |
DOI |
23 |
DOI |
24 |
DOI |
25 |
DOI |
26 |
|
27 |
DOI |
28 |
|
29 |
|
30 |
|
31 |
DOI |
32 |
|
33 |
|
34 |
|
35 |
|
36 |
|
[1] | 詹泓飞, 蔡振宁, 胡光辉. 基于虚时间演化与谱方法的一类基态Wigner函数计算方法[J]. 计算物理, 2022, 39(6): 651-665. |
[2] | 乃皮赛·吾买尔江, 闫好奎, 布玛丽亚·阿布力米提, 王丹琦, 向梅, 安桓. 外电场下CHBr3分子的光谱和解离特性[J]. 计算物理, 2022, 39(5): 624-630. |
[3] | 吉航, 孙仲谋, 周卓彦, 刘玉柱. 外电场对CO分子及离子性质的调制和降解[J]. 计算物理, 2022, 39(3): 327-334. |
[4] | 柴汝宽, 刘月田, 杨莉, 张艺馨, 辛晶, 马晶. 两种不同极性有机小分子在方解石(104)面吸附的密度泛函研究[J]. 计算物理, 2020, 37(2): 221-230. |
[5] | 陈余, 邢永明. 静水压力对Al14Mn2P16磁光性质影响的密度泛函理论研究[J]. 计算物理, 2020, 37(2): 231-239. |
[6] | 尹海峰, 曾春花, 陈文经. 二维二元碳化硅纳米结构的等离激元激发[J]. 计算物理, 2019, 36(5): 603-609. |
[7] | 何志伟, 张秀荣. (BN)25团簇的结构与性能[J]. 计算物理, 2019, 36(2): 219-224. |
[8] | 周康, 冯庆, 田芸, 李科, 周清斌. 过渡金属Cu、Cr掺杂TiO2表面氧化性气体NO2光学气敏传感特性[J]. 计算物理, 2018, 35(6): 702-710. |
[9] | 王啸卿, 刘玉柱, 尹文怡, 李金花. 电场作用下溴甲烷的光谱和解离特性[J]. 计算物理, 2018, 35(5): 619-625. |
[10] | 谢建明, 陈红霞, 庄国策. Mn掺杂(ZnSe)12团簇物性研究[J]. 计算物理, 2018, 35(4): 481-486. |
[11] | 李大伟, 高云亮, 朱芫江, 李进平. Ga掺杂δ-Pu的密度泛函理论计算[J]. 计算物理, 2018, 35(4): 487-493. |
[12] | 林艳, 刘芸, 彭清维, 魏晓楠, 唐延林. 基于密度泛函理论的5~10元瓜环结构参数与前线轨道的计算[J]. 计算物理, 2018, 35(2): 221-229. |
[13] | 陈红霞, 杜思洁, 庄国策. Ni掺杂ZnO团簇的结构和磁性质[J]. 计算物理, 2018, 35(1): 112-118. |
[14] | 李安军, 朱元强, 苏鸿, 杨泽宏. AuCl3催化炔酮分子内环化反应的机理[J]. 计算物理, 2017, 34(6): 679-684. |
[15] | 麻洪吉, 张国英, 方戈亮. 元素替代对MgH2储氢材料释氢影响机理的第一原理研究[J]. 计算物理, 2017, 34(6): 705-712. |
阅读次数 | ||||||
全文 |
|
|||||
摘要 |
|
|||||
版权所有 © 《计算物理》编辑部
地址:北京市海淀区丰豪东路2号 邮编:100094 E-mail:jswl@iapcm.ac.cn
本系统由北京玛格泰克科技发展有限公司设计开发