资源描述:
桥位吸附硼原子硅烯的结构与电子性质 ① 台俊杰, 肖 璇 (南昌航空大学 航空制造工程学院,江西 南昌 330063) 摘 要 基于密度泛函理论的第一性原理,研究了新型类石墨烯材料硅烯进行桥位吸附硼原子之后的微观结构和电子性质,并与 原始硅烯结构进行对比分析。 研究发现 3 种覆盖度下,44 桥位吸附硼原子硅烯超晶胞均为放热反应,吸附能最大为-2.56 eV。 吸 附过程中从硅原子到硼原子产生 1.14 个电子转移并形成稳定共价键。 相比于直接带隙半金属结构的原始硅烯,桥位吸附硼原子硅 烯在 K 点处打开 0.15 eV 直接带隙,同时引入 0.67 μB 磁矩。 关键词 密度泛函理论; 硅烯; 类石墨烯; 硼; 吸附;能带;态密度 中图分类号 O611.4; TB383文献标识码 Adoi10.3969/ j.issn.0253-6099.2020.05.031 文章编号 0253-6099(2020)05-0124-04 Structure and Electronic Properties of Silicene with Boron Atoms Adsorbed at Bridge Site TAI Jun⁃jie, XIAO Xuan (School of Aeronautical Manufacturing Engineering, Nanchang Hangkong University, Nanchang 330063, Jiangxi, China) Abstract Based on the first principles of density functional theory (DFT), the microstructure and electronic properties of a kind of new graphene⁃like material, the silicene with boron atoms adsorbed at bridge site, were studied and also compared with the original silicene structure. It is found that under 3 coverage conditions, the silicene supercells with boron atoms adsorbed at 44 bridge sites all experience exothermic reactions with the maximum adsorption energy of -2.56 eV. During the adsorption process, 1.14 electrons are transferred from silicon atom to boron atom and a stable covalent bond is formed. Compared with the original silicene that is the semi⁃metal structure with a direct band gap, the new silicene with boron atoms adsorbed at bridge site has opened a 0.15 eV direct band gap at K point and introduces a 0.67 μB magnetic moment. Key words density functional theory (DFT); silicene; graphene⁃like; boron; adsorption; energy band; density of states 相比于难以调控带隙的石墨烯,较大原子半径的 类石墨烯硅烯因较强的自旋⁃轨道耦合会直接在狄拉 克锥处打开约 0.15 meV 的微小带隙[1-2];同时也相对 更方便采用施加电场、掺杂或吸附等方式对硅烯进行 带隙调控[3-6];更有利的是,硅烯与现今硅基半导体技 术更为兼容,载流子迁移率高达 106m/ s。 在打开较大 硅烯带隙的同时保留极高的载流子迁移率是硅烯面向 半导体元器件大规模应用的关键。 本文计算研究了在 不同覆盖度下桥位吸附硼原子硅烯的结构演变和电子 性质,探索桥位吸附硼原子对结构稳定性和能带打开 情况的影响,研究结果可为后续实验调控硅烯电子结 构提供理论指导。 1 模型与计算方法 采用基于密度泛函理论的第一性原理计算方法, 利用带有自旋极化的平面波赝势法(PAW)描述离子 和电子间的相互作用并采用 Perdew Burke Ernzerhof (PBE)泛函[7-10]。 在布里渊区采样 991 作为 K 点 进行结构弛豫,并采用 24241 进行电子自洽计 算[11],平面波截断能设定为 350 eV。 Si 和 B 原子的 价态分别为 3s23p2和 3s23p3。 计算选取 44 超晶胞,为了消除层间相互作用的 ①收稿日期 2020-04-21 基金项目 航空科学基金(2015ZF56024) 作者简介 台俊杰(1994-),男,安徽六安人,硕士研究生,主要研究方向为新型二维材料理论计算。 通讯作者 肖 璇(1979-),男,湖北随州人,博士,副教授,主要研究方向为新型电子功能材料理论计算、先进高温结构材料的强韧化。 第 40 卷第 5 期 2020 年 10 月 矿矿 冶冶 工工 程程 MINING AND METALLURGICAL ENGINEERING Vol.40 №5 October 2020 影响,在超胞的 C 轴方向真空层取 2 nm[12]。 自洽计 算能量收敛标准设定为 1.010 -5 eV,力的收敛标准设 定为 1.010 -3 eV。 图 1 为硅烯吸附前后的结构图。 弛豫后原始硅烯晶格常数为 0.386 nm,SiSi 键长为 0.228 nm,翘曲高度 0.043 7 nm,与文献报道[12]基本一 致。 与石墨烯的平面 sp2杂化不同,由于硅烯原子翘 曲结构,相邻 Si 原子的 pz轨道只有部分相互重叠,所 以 SiSi 成键为 sp2和 sp3杂化。 在硅烯翘曲面上,共 有 4 个高对称吸附位,分别为顶、穴、桥和隙,如图 1(c) 所示,选取桥位作为计算对比吸附位。 吸附距离设定 过远或太近都不利于结构弛豫,由于 B 原子获电子能 力较强,初始吸附高度可以设定为0.25 nm。 为了确定 不同 Si 原子数量的超胞对 B 原子吸附的影响,B 原子 吸附于硅烯的覆盖度定义为吸附 B 原子个数与超胞 硅烯的原子数量之比,由此采用(22)、(33) 和 (44)共 3 种超胞来分别研究 B 原子在 1/8、1/18 和 图 1 硅烯吸附前后的结构 (a) 吸附前原始硅烯俯视图; (b) 吸附前原始硅烯侧视图; (c) 硅烯吸 附 B 原子的高对称点; (d) 硅烯桥位吸附 B 原子结构弛豫后的侧视图 1/32 这 3 种覆盖度吸附情况下硅原子数量对硼原子 吸附超胞的影响。 原子在硅烯上吸附能 EB [13]定义为 EB = E System - (ESilicene+ EAdatom) (1) 式中 ESystem是吸附后的总能量;ESilicene是吸附前硅烯体 系能量;EAdatom为吸附的单个 B 原子能量。 化学吸附 伴随电荷转移,利用 Bader 电荷分析计算体系在吸附 过程中的电荷转移情况,从而判断 B 原子在不同覆盖 度的电荷转移量。 电荷密度差分 Δρ[14]定义为 Δρ = ρTot - ρ Silicene - ρ B (2) 式中 ρTot为吸附后体系进行电子自洽计算的电荷密度; ρSilicene为吸附基底硅烯的电荷密度; ρB为结构吸附的 B 原子电荷密度。 2 结果与讨论 2.1 结 构 表 1 为桥位吸附 B 原子硅烯经过弛豫后的结构参 数。 从表 1 可以看出,吸附前 SiB 距离为 0.25 nm,3 种覆盖度下吸附硼的硅烯 SiB 键长均为 0.19 nm,相 比于本文设置的原始吸附高度在结构弛豫后 SiB 距 离缩短。 原始硅烯结构经弛豫后的厚度为 0.043 7 nm, 吸附 B 原子后的硅烯层厚度为 0.124 nm,含有吸附原 子的结构总厚度为 0.24 nm(图 1)。 表 1 桥位吸附 B 原子硅烯经弛豫之后的结构参数 位置dSix/ nmEB/ eVρB/ eVμ/ μB 22 桥位0.19-2.911.130.49 33 桥位0.19-2.711.160.96 44 桥位0.19-2.561.140.67 注dSix表示弛豫后的 SiB 键长;EB表示吸附能; ρB表示在吸附弛 豫后 B 原子的获得电子数量;μ 表示弛豫后结构的磁矩。 经结构优化后的 3 种覆盖度下的桥位吸附 B 原 子后的 SiB 键长相同,吸附能均为负值,分别为-2.91、 -2.71 和-2.56 eV,说明表面吸附过程为放热反应,吸 附能随着覆盖度降低而减小,覆盖度为 1/8 和 1/32 的 吸附能仅相差 0.35 eV。 从表 1 还可以看出,桥位吸附 B 原子形成的复合物均呈现出一定的磁性。 值得注意 的是,对谷位、顶位等其他吸附位置的计算结果表明 μ 值为 0,均无磁性。 同时,化学吸附伴随着电荷转移,3 种覆盖度下电荷转移数量相差不大,因此随后对电荷 转移情况描述的差分电荷密度结构图只给出覆盖度为 1/32 的 44 桥位吸附结果。 2.2 能带结构和电子态密度 图 2 为原始硅烯及 44 桥位吸附 B 原子硅烯的 521第 5 期台俊杰等 桥位吸附硼原子硅烯的结构与电子性质 能带和态密度电子结构。 从图 2 可以看出,原始硅烯 的能带结构与石墨烯结构类似,均是带隙大小为 0 的半 金属结构。 其导带底和价带顶在布里渊区中的 K 点相 接触,能带表现为线性色散[15]。 从原始硅烯的总态密 度结构可以看出上下自旋结构对称,因此原始硅烯中没 有产生磁矩。 另外,图 2(c)表明吸附前 SiSi 键之间 除了形成 sp3杂化的 σ 键,在相邻的 pz轨道部分还重叠 形成 π 键,这会使硅烯的平面翘曲结构更加稳定[11]。 5 0 -5 -10 GMKG 能量/eV 2 0 -2G MKG 能量/eV 能量/eV 2 1 0 -1 -2 -5-1005 电子/eV b 能量/eV 15 10 5 0 -5 -10 -15 -10-505 电子/eV e a d 能量/eV 0.4 0.2 0.0 -0.2 -0.4 -10-505 电子/eV c 能量/eV 4 2 0 -0.2 -0.4 -1005-5 电子/eV f 0.15 eV Si - s Si - px Si - py Si - pz Si - s Si - px Si - py Si - pz Tot B Si 图 2 原始硅烯和 44 桥位吸附 B 原子硅烯的能带、态密度电子结构 (a) 原始硅烯的能带结构; (b) 原始硅烯的总态密度; (c) 原始硅烯的 s、px、py和 pz轨道的分态密度; (d) 44 桥位吸附 B 原子结构弛豫后的能带 结构; (e) 44 桥位吸附 B 原子硅烯弛豫后结构总态密度和 B、Si 原子总态密度; (f) 44 桥位吸附 B 原子硅烯结构弛豫后 Si 的 s、px、py和 pz轨道 的分态密度 44 桥位吸附 B 原子结构弛豫后硅烯的能带图计 算结果表明,在 M、K 和 G 高对称点处分别打开了0.37 eV,0.15 eV 和 0.44 eV 带隙,狄拉克锥并未消失,价带 顶接触到费米面,表现为 P 型吸附的特性。 吸附后的 SiB 键长比 SiSi 键长短,SiB 之间有电子转移并 形成稳定的共价键,B 在吸附过程中得 1.14 个电子 (表 1)。 比较图 2(c)和图 2(f)可以看出,原始硅烯的 导带电子主要由 pz轨道构成,价带则由 s 轨道构成。 桥位吸附 B 原子后的 pz轨道电荷向 B 原子转移,原有 的 sp2杂化轨道减弱,能带结构发生变化。 在图 2(e) 和(f)中的总态密度和分态密度的上下自旋均不对称, 结构引入 0.67 μB 磁矩。 在桥位吸附 B 原子结构中, 33 桥位吸附所产生的磁矩最大为 0.96 μB(表 1)。 上述能带和态密度分析结果表明,B 原子吸附对硅烯 体系的电子结构产生显著影响,原始硅烯的半金属能 带结构在吸附 B 原子之后打开了 0.15 eV 的直接带 隙,B 原子吸附硅烯不仅为半导体工业应用调控硅烯 带隙提供了可行性,同时还产生了磁矩。 2.3 差分电荷密度和电子局域密度 为了进一步分析硅烯桥位吸附硼原子的作用机 理,观察硅烯和硼原子之间的电荷转移情况,计算了桥 位吸附硼原子的差分电荷密度以及电子局域密度[16]。 差分电荷结果如图 3(a) ~ (b)所示,其中黄色区域为 得电子区域,蓝色为失电子区域,可以看出吸附 B 原 子后主要影响区域为 B 原子周围最近的 6 个 Si 原子, 电荷主要从 Si 原子流向 B 原子附近局域。 吸附后 SiB 键长为 0.19 nm,原有的 SiSi 键被破坏结构产 生内应力。 图 3(c) ~(d)为原始硅烯和44 桥位吸附 B 原子 硅烯平面电荷密度,对比发现,吸附 B 原子之后 Si 原 子层的结构发生改变,由原厚度 0.437 nm 扭折为 0.124 nm。 图 3(e)为图 3(d)中椭圆区域的电荷密度 切片,可以清晰看出成键情况。 X 区域为 Si 原子 pz轨 道相互部分重叠形成的 sp2轨道杂化和 Y 区域的 sp3 轨道杂化共同形成 SiSi 共价键;Z 部分为 B 原子得 电子 Si 原子失电子形成的强共价键。 621矿 冶 工 程第 40 卷 图 3 44 桥位吸附 B 原子硅烯的差分电荷密度以及电子局域 密度切片结构 (a) 44 桥位吸附 B 原子差分电荷结构俯视图; (b) 44 桥位吸附 B 原子差分电荷结构侧视图; (c) 原始硅烯平面电子局域密度; (d) 44 桥位吸附 B 原子后 Si 原子电子局域密度; (e) 44 桥位吸附 B 原子后 沿 B 原子垂直方向电子局域密度 3 结 论 1) 在 1/8、1/18 和 1/32 这 3 种覆盖度下桥位吸 附 B 原子的硅烯均为放热反应,吸附能为负值。 吸附 过程每个超胞中会形成稳定的 SiB 共价键,产生 1.14 个电子转移,Si 原子为失电子,B 原子为得电子。 2) 44 桥位吸附 B 原子硅烯在 K 点处可以打开 0.15 eV 直接带隙,远大于半金属结构的原始硅烯。 3) 电子态密度的上下自旋不对称使得 3 种覆盖 度下桥位吸附 B 原子硅烯均呈现铁磁性。 参考文献 [1] Lalmi B, Oughaddou H, Enriquez H, et al. Epitaxial growth of a sili⁃ cene sheet[J]. Applied Physics Letters, 2010,97(22)223109. [2] Hikami Shinobu, Larkin Anatoly I, Nagaoka Yosuke. Spin⁃orbit Interac⁃ tion and magnetoresistance in the two dimensional random system[J]. Progress of Theoretical Physics, 1980,63(2)357-742. [3] Song Z, Liu C C, Yang J, et al. Quantum spin hall insulators and quantum valley hall insulators of BiX/ SbX(X = H, F, Cl and Br) monolayers with a record bulk band gap[J]. NPG Asia Materials, 2014,6(12)e147. [4] Farghadan, Rouhollah. Bipolar magnetic semiconductor in silicene na⁃ noribbons[J]. Journal of Magnetism and Magnetic Materials, 2017, 435(8)206-211. [5] Zhang Jian⁃Min, Song Wan⁃Ting, Xu Ke⁃Wei, et al. The study of the P doped silicene nanoribbons with first⁃principles[J]. Computational Materials Science, 2014,95429-434. [6] 付 博,李龙江,张 覃. 白云石电子结构及油酸钠在白云石 (104)面吸附作用的分子模拟研究[J]. 矿冶工程, 2019,39(5) 69-73. [7] Parr R G. Density functional theory of atoms and molecules[J]. Pro⁃ ceedings of the National Academy of Sciences of the United States of America, 1995,191(2)277-278. [8] 张小武,谢 俊,李龙江,等. 氟磷灰石、氯磷灰石和羟基磷灰石电 子结构的第一性原理研究[J]. 矿冶工程, 2017,37(6)84-90. [9] Chl P E B. Projector augmented⁃wave method[J]. Physical Review B, 1994,5017953-17979. [10] Perdew J P, Burke K, Ernzerhof M. Generalized gradient approxi⁃ mation made simple[J]. Physical Review Letters, 1996,77(18) 3865-3868. [11] Sivek J, Şahin Hasan, Partoens B, et al. Adsorption and absorption of boron, nitrogen, aluminium and phosphorus on silicene stability, electronic and phonon properties[J]. Physical Review B, 2013,87 (8)085444. [12] Hernndez Cocoletzi, H Castellanos guila. DFT studies on the Al, B, and P doping of silicene[J]. Superlattices and Microstructures, 2018,114242-250. [13] Sanchezpaisal Y, Sanchezportal D, Ayuela A. Ab initio calculations of zirconium adsorption and diffusion on graphene[J]. Physical Review B, 2009,80(4)045428. (下转第 133 页) 721第 5 期台俊杰等 桥位吸附硼原子硅烯的结构与电子性质 3 结 论 1) 通过热模拟压缩试验获得了 TC18 钛合金在 790~900 ℃下的真应力⁃真应变曲线。 曲线有明显峰 值应力,试样在热压缩过程中发生了动态再结晶,采用 Arrhenius 函数模型建立了 TC18 钛合金的热变形本构 方程。 2) 在 880 ℃、̇ ε=0.01 s -1 变形条件下,有限元模拟 得到双圆锥台试样截面应变呈梯度分布。 3) 双圆锥台试样不同应变区域 α 相晶粒尺寸与 相同应变条件下圆柱形热模拟试样 α 相晶粒尺寸较 为一致,说明双圆锥台试样热变形有限元模拟所得应 变分布结果较为准确,一个双圆锥台试样上可以快速 获取多个不同应变条件的热变形组织。 参考文献 [1] 莱茵斯皮特尔斯. 钛与钛合金[M]. 陈振华,译. 北京化学工 业出版社, 2005. [2] Lin Y C, Huang J, He D G, et al. Phase transformation and dynamic recrystallization behaviors in a Ti55511 titanium alloy during hot com⁃ pression[J]. Journal of Alloys and Compounds, 2019,795471-482. [3] Zherebtsov S V, Murzinova M A, Klimova M V, et al. Microstructure e⁃ volution during warm working of Ti⁃5Al⁃5Mo⁃5V⁃1Cr⁃1Fe at 600 and 800 ℃[J]. Materials Science and Engineering A, 2013,563168-176. [4] Tan L M, Huang Z W, Liu F, et al. Effects of strain amount and strain rate on grain structure of a novel high Co nickel⁃based polycrys⁃ talline superalloy[J]. Materials & Design, 2017,13160-68. [5] 刘岚逸,汪冰峰,张晓泳. TC18 钛合金中剪切局域化及其微观结 构研究[J]. 矿冶工程, 2018,38(4)149-151. [6] Luo J, Wang L F, Liu S F, et al. The correlation between the flow behavior and the microstructure evolution during hot working of TC18 alloy[J]. Materials Science and Engineering A, 2016,654213-220. [7] 尹 勇,谭 真,肖振宇. 微量 Ca、Y 元素对 AZ91 镁合金动态再 结晶及力学性能的影响[J]. 矿冶工程, 2019,39(1)114-117. [8] 郭俊成,肖振宇,杨续跃. 两步锻压提升铸态 AZ80 镁合金延伸率 的研究[J]. 矿冶工程, 2018,38(4)144-148. [9] Sakai T, Belyakov A, Kaibyshev R, et al. Dynamic and post⁃dynamic recrystallization under hot, cold and severe plastic deformation condi⁃ tions[J]. Progress in Materials Science, 2014,60(1)130-207. [10] Jia B H, Song W D, Tang H P, et al. Hot deformation behavior and constitutive model of TC18 alloy during compression[J]. Rare Met⁃ als, 2014,33(4)383-389. [11] Sellars C M, Mctegart W J. On the mechanism of hot deformation[J]. Acta Metallurgica, 1966,14(9)1136-1138. [12] Zener C, Hollomon J H. Effect of strain rate upon plastic flow of steel[J]. Journal of Applied Physics, 1944,15(1)22-32. [13] 彭柳锋,郭道强,魏玉勇,等. 2A12 铝合金锻件成形的有限元模 拟[J]. 矿冶工程, 2017,37(5)99-104. [14] 张 莉,李升军. DEFORM 在金属塑性成形中的应用[M]. 北 京机械工业出版社, 2009. 引用本文 毛 颖,雷 雨,丁 旭,等. TC18 钛合金双圆锥台试样热变 形有限元模拟研究[J]. 矿冶工程, 2020,40(5)128-133. (上接第 123 页) [8] Lv K, Li J, Qing X, et al. Synthesis and photo⁃degradation applica⁃ tion of WO3/ TiO2hollow spheres[J]. Journal of Hazardous Materials, 2011,189(1-2)329-335. [9] 曾爱香,熊惟浩,王采芳,等. 溶胶⁃凝胶法制备空心微珠表面钡铁 氧体包覆层的研究[J]. 材料保护, 2004,37(9)19-20. [10] 张 姣,江学良,余 露,等. 水热法制备二氧化铈纳米空心球及 其吸附性能研究[J]. 材料研究学报, 2016,30(5)365-371. [11] 黄 彬,江 奇,谢德钰,等. 模板法制备聚苯胺空心微球材料及 其电化学性能[J]. 功能材料, 2011,42(1)144-147. [12] 张 珂,汪太生,李自力,等. 细乳液聚合法制备聚合物空心球[J]. 高分子学报, 2018(4)475-481. [13] 叶敏伟. 新型粉末喷涂装备[J]. 电镀与涂饰, 2011,30(10)73-75. [14] 吴子健. 热喷涂技术与应用[M]. 北京 机械工业出版社, 2006. 引用本文 张胜全,何亮亮,王 胜,等. 铜渣尾矿制备空心陶瓷微球及 其表征[J]. 矿冶工程, 2020,40(5)120-123. (上接第 127 页) [14] Fehir J, Richard J. Differential polarization of spin and charge densi⁃ ty in substituted phenoxy radicals[J]. Journal of Physical Chemistry A, 2009,113(32)9249-9260. [15] Vogt P, Padova P D, Quaresima C, et al. Silicene compelling ex⁃ perimental evidence for graphenelike two⁃dimensional silicon[J]. Physical Review Letters, 2012,108(15)155501. [16] Mizoguchi H, Okunaka M, Kitano M, et al. Hydride⁃based electride material, LnH2( Ln = La, Ce, or Y) [ J]. Inorganic Chemistry, 2016,55(17)8833-8838. 引用本文 台俊杰,肖 璇. 桥位吸附硼原子硅烯的结构与电子性质[J]. 矿冶工程, 2020,40(5)124-127. 331第 5 期毛 颖等 TC18 钛合金双圆锥台试样热变形有限元模拟研究
展开阅读全文