高级搜索

铁原子吸附联苯烯单层电子结构的第一性原理

,  ,  ,  ,  , 

downloadPDF

上一篇

下一篇

吴洪芬, 冯盼君, 张烁, 刘大鹏, 高淼, 闫循旺. 铁原子吸附联苯烯单层电子结构的第一性原理[J]. 物理学报, 2022, 71(3): 036801-1. doi: 10.7498/aps.71.20211631
引用本文: 吴洪芬, 冯盼君, 张烁, 刘大鹏, 高淼, 闫循旺. 铁原子吸附联苯烯单层电子结构的第一性原理[J]. 物理学报, 2022, 71(3): 036801-1. doi: 10.7498/aps.71.20211631
shu
Hong-Fen Wu, Pan-Jun Feng, Shuo Zhang, Da-Peng Liu, Miao Gao, Xun-Wang Yan. First-principles study of Fe atom adsorbed biphenylene monolayer[J]. Acta Physica Sinica, 2022, 71(3): 036801-1. doi: 10.7498/aps.71.20211631
Citation: Hong-Fen Wu, Pan-Jun Feng, Shuo Zhang, Da-Peng Liu, Miao Gao, Xun-Wang Yan. First-principles study of Fe atom adsorbed biphenylene monolayer[J]. Acta Physica Sinica, 2022, 71(3): 036801-1. doi: 10.7498/aps.71.20211631
shu

铁原子吸附联苯烯单层电子结构的第一性原理

  • 1.

    曲阜师范大学物理工程学院, 曲阜　273165

  • 2.

    宁波大学物理科学与技术学院物理系, 宁波　 315211

    • 通讯作者: E-mail: yanxunwang@163.com

  • 中图分类号: 68.43.-h, 81.05.ue, 63.20.dk, 75.70.Cn

First-principles study of Fe atom adsorbed biphenylene monolayer

  • 1.

    School of Physics and Physical Engineering, Qufu Normal University, Qufu 273165, China

  • 2.

    Department of Physics, School of Physical Science and Technology, Ningbo University, Ningbo 315211, China

    • Corresponding author: E-mail. yanxunwang@163.com

  • MSC: 68.43.-h, 81.05.ue, 63.20.dk, 75.70.Cn

  • 摘要: 联苯烯单层由碳原子的四元、六元和八元环组成, 具有与石墨烯相似的单原子层结构. 2021年5月, Science首次报道了该材料的实验合成, 引起了科研工作者的极大关注. 基于第一性原理的密度泛函方法, 研究了铁原子在联苯烯单层的吸附构型并分析了其电子结构. 结构优化、吸附能和分子动力学的计算表明, 联苯烯单层的四元环空位是铁原子最稳定的吸附位点, 吸附能可达1.56 eV. 电子态密度计算表明铁3d电子与碳的2p电子有较强的轨道杂化, 同时电荷转移计算显示铁原子向近邻碳原子转移的电荷约为0.73个电子, 说明联苯烯单层与吸附的铁原子之间形成了稳定的化学键. 另外, 铁原子吸附于联苯烯单层后体系显磁性, 铁原子上局域磁矩大小约为 1.81 μB, 方向指向面外. 因此, 本文确认了联苯烯单层是比石墨烯更好的铁原子吸附载体且体系有磁性, 这为研究吸附材料的电磁、输运、催化等特性提供了新的平台.
    Abstract: Biphenylene monolayer is composed of four-, six- and eight-membered carbon rings and has a monatomic layer structure similar to graphene. It was synthesized in experiment recently and reported in Science in May 2021, which has attracted considerable attention in the research field of two-dimensional materials. By the density functional method of the first principle, we study the adsorption configuration of Fe atoms on biphenylene monolayer and analyze its electronic structure. The calculation of structural optimization, adsorption energy and molecular dynamics show that the biphenylene monolayer is a good matrix of Fe atoms. For Fe atoms, the hollow site in the four-membered ring of the biphenylene monolayer is the most stable adsorption site, and the adsorption energy can reach 1.56 eV. The calculation of charge transfer and density of states show that a stable bond can be formed between biphenylene monolayer and Fe atoms, and 0.73 electron is transferred from Fe atom to the neighbored carbon atom. After Fe atom being absorbed, biphenylene monolayer is magnetic, and the magnetic moment of Fe atom is about 1.81 ${\mu}_{\mathrm{B}}$ and points out of the plane. Compared with graphene, biphenylene monolayer adsorbs Fe atoms more stably, which provides a new platform for studying the electromagnetic, transport and catalytic properties of two-dimensional materials with adatoms.
  • 加载中

  • 图 1  (a)石墨烯和(b)联苯烯单层及其需考虑的吸附位点; (c), (f); (d), (g); (e), (h) 分别是铁原子吸附于联苯烯单层四、六和八元环空位(H位点)的俯视图及侧视图

    Figure 1.  (a) and (b) Graphene and biphenene networks and the adsorption sites of Fe atoms; (c) and (f), (d) and (g), (e) and (h) are top and side views of Fe atoms on the top of 4-, 6-, and 8-membered ring (H sites) of biphenene networks, respectively.

    下载: 全尺寸图片 幻灯片

    图 2  (a), (b)和(c)是在1000 K下, 铁原子分别吸附于联苯烯单层的四、六和八元环后能量随时间的演化曲线; (d) 在1000 K下, 铁原子在八元环空位的运动情况

    Figure 2.  (a), (b) and (c) Time evolution curves of the energy of Fe atoms adsorbed on the 4-, 6- and 8-membered rings of the biphenene network at 1000 K, respectively; (d) motion of iron atom absorbed on hollow site of 8-membered ring at 1000 K.

    下载: 全尺寸图片 幻灯片

    图 3  (a), (b); (c), (d); (e), (f)分别是铁原子吸附于联苯烯单层四、六、八元环空位(H位点)情况下的差分电荷俯视图及侧视图

    Figure 3.  (a), (b); (c), (d); (e), (f) are top and side views of differential charge when Fe atoms are adsorbed on top of 4-, 6- and 8-membered ring (H sites) of biphenene networks, respectively.

    下载: 全尺寸图片 幻灯片

    图 4  (a)石墨烯和(b)联苯烯单层的能带及态密度图(费米能级设为0 eV)

    Figure 4.  Energy band and density of states of (a) graphene and (b) biphenene network. The Fermi energy is set to zero.

    下载: 全尺寸图片 幻灯片

    图 5  铁原子吸附于联苯烯单层四元环空位时的投影态密度图(费米能级为零)　(a), (c) GGA计算; (b), (d) GGA + U计算

    Figure 5.  Projected density of states of Fe atoms adsorbed on the top of 4-membered ring of biphenene network: (a) and (c) GGA calculations; (b) and (d) GGA + U calculations. The Fermi energy is set to zero.

    下载: 全尺寸图片 幻灯片

    图 6  联苯烯单层均匀吸附铁原子后的几种磁序结构 (a)铁磁(FM); (b) 共线反铁磁序一(Coll-I); (c) 共线反铁磁序二(Coll-II); , (d) 奈尔反铁磁序(Nèel)

    Figure 6.  Sketches of several magnetic orders in Fe-adsorbed biphenene monolayer: (a) Ferromagnetic order; (b) collinear anti-ferromagnetic order I; (c) collinear anti-ferromagnetic order II; (d) Nèel antiferromagnetic order.

    下载: 全尺寸图片 幻灯片

    图 7  自洽计算(蓝线)和非自洽计算(红线)铁原子d电子数对Hubbard U的响应.

    Figure 7.  Response of d electron number of Fe atom adsorbed on biphenene network to Hubbard U in self-consistent and non-self-consistent calculations.

    下载: 全尺寸图片 幻灯片

    表 1  铁原子吸附在石墨烯和联苯烯单层各位点的吸附能

    Table 1.  Adsorption energy of Fe atom adsorbed on each point of graphene and biphenene.

    吸附位点吸附能/eV吸附位点吸附能/eV
    石墨烯H0.84联苯烯单层B11.29
    石墨烯B0.28联苯烯单层B21.16
    石墨烯T0.16联苯烯单层B30.88
    联苯烯单层H11.56联苯烯单层B40.88
    联苯烯单层H21.53联苯烯单层T10.94
    联苯烯单层H31.12联苯烯单层T20.96
    下载: 导出CSV

    表 2  铁原子吸附联苯烯单层不同磁序下的能量

    Table 2.  Energies of Fe atom adsorbed biphenene layer in various magnetic orders.

    FM/eVColl-I/eVColl-II/eVNèel/eV
    GGA–230.673–230.309–230.599–230.683
    GGA + U–222.790–222.799–223.818–224.042
    下载: 导出CSV
  • [1] Novoselov K S, Geim A K, Morozov S V, Jiang D, Zhang Y, Dubonos S V, Grigorieva I V, Firsov A A 2004 Science 306 666 doi: 10.1126/science.1102896
    [2] Fan Q T, Yan L H, Tripp M W, Krejci O, Dimosthenous S, Kachel S R, Chen M, Foster A S, Koert U, Liljeroth P, Gottfried J M 2021 Science 372 852 doi: 10.1126/science.abg4509
    [3] Zhang R S, Jiang X W 2019 Front. Phys. 14 13401 doi: 10.1007/s11467-018-0836-5
    [4] Tang C, Kour G, Du A J 2019 Chin. Phys. B 28 107306 doi: 10.1088/1674-1056/ab41ea
    [5] Liu D P, Zhang S, Gao M, Yan X W, Xie Z Y 2021 Appl. Phys. Lett. 118 223104 doi: 10.1063/5.0054730
    [6] Liu D P, Zhang S, Gao M, Yan X W 2021 Phys. Rev. B 103 125407 doi: 10.1103/PhysRevB.103.125407
    [7] Liu D P, Feng P J, Gao M, Yan X W 2021 Phys. Rev. B 103 155411 doi: 10.1103/PhysRevB.103.155411
    [8] Baughman R H, Eckhardt H, Kertesz M 1987 J. Chem. Phys. 87 6687 doi: 10.1063/1.453405
    [9] Narita N, Nagai S, Suzuki S, Nakao K 1998 Phys. Rev. B 58 11009 doi: 10.1103/PhysRevB.58.11009
    [10] Long M Q, Tang L, Wang D, Li Y L, Shuai Z G 2011 ACS Nano. 5 2593 doi: 10.1021/nn102472s
    [11] Li G X, Li Y L, Liu H B, Guo Y B, Li Y J, Zhu D B 2010 Chem. Commun. 46 3256 doi: 10.1039/b922733d
    [12] Song Q, Wang B, Deng K, Feng X L, Wagner M, Gale J D, Mullen K, Zhi L J 2013 J. Mater. Chem. C 1 38
    [13] Liu W, Miao M S, Liu J Y 2015 RSC Adv. 5 70766 doi: 10.1039/C5RA11208G
    [14] Hudspeth M A, Whitman B W, Barone V, Peralta J E 2010 ACS Nano. 4 4565 doi: 10.1021/nn100758h
    [15] Karaush N N, Bondarchuk S V, Baryshnikov G V, Minaeva V A, Sun W H, Minaev B F 2016 Rsc Adv. 6 49505 doi: 10.1039/C6RA06832D
    [16] Konstantinova E, Dantas S O, Barone P M V B 2006 Phys. Rev. B 74 35417 doi: 10.1103/PhysRevB.74.035417
    [17] Zhang S H, Zhou J, Wang Q, Chen X S, Kawazoe Y, Jena P 2015 Proc. Natl. Acad. Sci. U.S.A. 112 2372 doi: 10.1073/pnas.1416591112
    [18] Mandal B, Sarkar S, Pramanik A, Sarkar P 2013 Phys. Chem. Chem. Phys. 15 21001 doi: 10.1039/c3cp53390e
    [19] Deza M, Fowler P W, Shtogrin M, Vietze K 2000 J. Chem. Inf. Comput. Sci. 40 1325 doi: 10.1021/ci000010j
    [20] Terrones H, Terrones M, Hernandez E, Grobert N, Charlier J C, Ajayan P M 2000 Phys. Rev. Lett. 84 1716 doi: 10.1103/PhysRevLett.84.1716
    [21] Bucknum M J, Castro E A 2008 Solid State Sci. 10 1245 doi: 10.1016/j.solidstatesciences.2007.11.035
    [22] Zhu H Y, Balaban A T, Klein D J, Zivkovic T P 1994 J. Chem. Phys. 101 5281 doi: 10.1063/1.467382
    [23] Wang X Q, Li H D, Wang J T 2012 Phys. Chem. Chem. Phys. 14 11107 doi: 10.1039/c2cp41464c
    [24] Liu Y, Wang G, Huang Q S, Guo L W, Chen X L 2012 Phys. Rev. Lett. 108 225505 doi: 10.1103/PhysRevLett.108.225505
    [25] Su C, Jiang H, Feng J 2013 Phys. Rev. B 87 075453 doi: 10.1103/PhysRevB.87.075453
    [26] Tang C P, Xiong S J 2012 AIP Adv. 2 042147 doi: 10.1063/1.4768669
    [27] Nulakani N V R, Kamaraj M, Subramanian V 2015 RSC Adv. 5 78910 doi: 10.1039/C5RA14477A
    [28] Wang X Q, Li H D, Wang J T 2013 Phys. Chem. Chem. Phys. 15 2024 doi: 10.1039/C2CP43070C
    [29] Wang Z H, Zhou X F, Zhang X M, Zhu Q, Dong H F, Zhao M W, Oganov A R 2015 Nano Lett. 15 6182 doi: 10.1021/acs.nanolett.5b02512
    [30] Pereira L F C, Mortazavi B, Makaremi M, Rabczuk T 2016 RSC Adv. 6 57773 doi: 10.1039/C6RA05082D
    [31] Chopra S 2016 RSC Adv. 6 89934 doi: 10.1039/C6RA21155K
    [32] Li Q D, Li Y, Chen Y, Wu L L, Yang C F, Cui X L 2018 Carbon 136 248 doi: 10.1016/j.carbon.2018.04.081
    [33] Tahara K, Yamamoto Y, Gross D E, Kozuma H, Arikuma Y, Ohta K, Koizumi Y, Gao Y, Shimizu Y, Seki S, Kamada K, Moore J S, Tobe Y 2013 Chem. Eur. J. 19 11251 doi: 10.1002/chem.201300838
    [34] Balaban A T, Rentia C C, Ciupitu E 1968 Rev. Roum. Chim. 13 231
    [35] Li X X, Yang J L 2016 Natl. Sci. Rev. 3 365 doi: 10.1093/nsr/nww026
    [36] Yang W J, Gao Z Y, Liu X S, Ma C Z, Ding X L, Yan W P 2019 Fuel 243 262 doi: 10.1016/j.fuel.2019.01.125
    [37] Krasheninnikov A V, Lehtinen P O, Foster A S, Pyykko P, Nieminen R M 2009 Phys. Rev. Lett. 102 126807 doi: 10.1103/PhysRevLett.102.126807
    [38] Peng Y, Lu B Z, Chen S W 2018 Adv. Mater. 30 1801995 doi: 10.1002/adma.201801995
    [39] Luo X, Wei W Q, Wang H J, Gu W L, Kaneko T, Yoshida Y, Zhao X, Zhu C Z 2020 Nano-Micro. Lett. 12 163 doi: 10.1007/s40820-020-00502-5
    [40] Wei X Q, Song S J, Wu N N, Luo X, Zheng L R, Jiao L, Wang H J, Fang Q, Hu L, Y, Gu W L, Song W Y, Zhu C Z 2021 Nano Energy 84 105840 doi: 10.1016/j.nanoen.2021.105840
    [41] Kresse G, Hafner J 1993 Phys. Rev. B 47 558 doi: 10.1103/PhysRevB.47.558
    [42] Kresse G, Furthmüller J 1996 Phys. Rev. B 54 11169 doi: 10.1103/PhysRevB.54.11169
    [43] Blöchl P E 1994 Phys. Rev. B 50 17953 doi: 10.1103/PhysRevB.50.17953
    [44] Perdew J P, Burke K, Ernzerhof M 1996 Phys. Rev. Lett. 77 3865 doi: 10.1103/PhysRevLett.77.3865
    [45] Cococcioni M, Gironcoli S D 2005 Phys. Rev. B 71 035105 doi: 10.1103/PhysRevB.71.035105
    [46] Grimme S 2006 J. Comp. Chem. 27 1787 doi: 10.1002/jcc.20495
    [47] Martyna G J, Klein M L, Tuckerman M 1992 J. Chem. Phys. 97 2635 doi: 10.1063/1.463940
    [48] Zanella I, Fagan S B, Mota R, Fazzio A 2008 J. Phys. Chem. C 112 9163 doi: 10.1021/jp711691r
    [49] Henkelman G, Arnaldsson A, Jónsson H 2006 Comput. Mater. Sci. 36 354 doi: 10.1016/j.commatsci.2005.04.010
    [50] Sun Y J, Zhuo Z W, Wu X J, Yang J L 2017 Nano Lett. 17 2771 doi: 10.1021/acs.nanolett.6b04884
    [51] Zhuang H L, Xie Y, Kent P R C, Ganesh P 2015 Phys. Rev. B 92 035407 doi: 10.1103/PhysRevB.92.035407
    [52] Hu W, Wang C, Tan H, Duan H L, Li G N, Li N, Ji Q Q, Lu Y, Wang Y, Sun Z H, Hu F C, Yan W S 2021 Nat. Commun. 12 1854 doi: 10.1038/s41467-021-22122-2
    [53] Xu K, Ding H, Lü H F, Chen P Z, Lu X L, Han Cheng H, Zhou T P, Liu S, Wu X J, Wu C Z, Xie Y 2016 Adv. Mater. 28 3326 doi: 10.1002/adma.201505732
  • 加载中
图( 7) 表( 2)
计量
  • 文章访问数:  397
  • HTML全文浏览数:  397
  • PDF下载数:  1
  • 施引文献:  0
出版历程
  • 收稿日期:  2021-09-02
  • 刊出日期:  2022-02-05

铁原子吸附联苯烯单层电子结构的第一性原理

    通讯作者: E-mail: yanxunwang@163.com
  • 1. 曲阜师范大学物理工程学院, 曲阜　273165
  • 2. 宁波大学物理科学与技术学院物理系, 宁波　 315211
  • 收稿日期:  2021-09-02

摘要: 联苯烯单层由碳原子的四元、六元和八元环组成, 具有与石墨烯相似的单原子层结构. 2021年5月, Science首次报道了该材料的实验合成, 引起了科研工作者的极大关注. 基于第一性原理的密度泛函方法, 研究了铁原子在联苯烯单层的吸附构型并分析了其电子结构. 结构优化、吸附能和分子动力学的计算表明, 联苯烯单层的四元环空位是铁原子最稳定的吸附位点, 吸附能可达1.56 eV. 电子态密度计算表明铁3d电子与碳的2p电子有较强的轨道杂化, 同时电荷转移计算显示铁原子向近邻碳原子转移的电荷约为0.73个电子, 说明联苯烯单层与吸附的铁原子之间形成了稳定的化学键. 另外, 铁原子吸附于联苯烯单层后体系显磁性, 铁原子上局域磁矩大小约为 1.81 μB, 方向指向面外. 因此, 本文确认了联苯烯单层是比石墨烯更好的铁原子吸附载体且体系有磁性, 这为研究吸附材料的电磁、输运、催化等特性提供了新的平台.

English Abstract

    全文HTML

    1.   引　言

    • 石墨烯由于其奇特的二维蜂窝状结构以及在电、磁、热、力学等方面独特的物理化学性质, 引起了科研工作者极大的兴趣. 自2004年石墨烯被制备以来[1], 探索新型的低维材料便一直是凝聚态物理领域的前沿课题, 其中二维的碳材料是该领域科研工作者的重点研究对象之一[2-7]. 近年来, 虽然理论上已经预测了20余种[8-31]碳的二维网格结构, 如石墨炔(graphyne) [8,9]、石墨双炔 (graphdiyne)[10,11]、石墨n炔(graphenylene)[12,13]、联苯烯单层(biphenylene) [14,15]、轴烯(radialenes)[16]等, 但遗憾的是, 到目前为止, 实验上所合成的二维碳材料非常少, 主要有石墨烯[1]、石墨双炔[11]、γ-石墨炔(γ-graphyne)[32]、tetrakis-DBAs [33]等.

      近日Fan等[2]在实验上成功制备了石墨烯二维同素异形体—联苯烯单层并发表在Science上, 引起相关科研工作者的广泛关注. 该材料最早由Balaban等[34]于1968年理论预测得到. 联苯烯单层有类石墨烯的二维单层结构, 二者均由碳元素组成, 每个原子都是三配位的碳原子, 其稳定性都来源于sp2杂化轨道的成键机制. 石墨烯是最典型的二维材料, 也是人们研究最广泛的二维材料, 具有独特的晶体结构和电子结构, 从而具备优异的力学、热学、电学性质, 在能源、传感器、半导体、生物医学等多个领域有广泛的应用潜力. 鉴于联苯烯单层与石墨烯的相似性, 联苯烯自身的研究具有重要的科学意义. 另外, 石墨烯的化学修饰进一步丰富了石墨烯基材料的功能, 如可作为电极材料、储能材料、催化材料、自旋电子器件材料等. 相比之下, 联苯烯的化学修饰会带来更多更复杂的物理特性. 这是因为联苯烯单层由四元环、六元环和八元碳环组成, 而石墨烯中只有六元环, 故联苯烯可以提供更多样化的化学修饰环境. 因此, 联苯烯及化学修饰的联苯烯材料有什么样的奇异物性很值得期待, 另一方面确定其基本电子结构也是当前的一个紧迫任务.

      铁原子作为一种极具代表性的过渡金属原子, 吸附或嵌入在二维单层中形成复合材料, 在自旋电子器件[35]、催化[36]、能源(电池)等诸多领域有着巨大的研究价值和广泛的应用前景. 以石墨烯为代表的二维碳材料吸附过渡金属原子后常表现出有趣的磁性质[37]; 此外, 碳基铁原子催化剂是氧还原等反应中活性最高的催化剂之一[38], 由于极高的原子利用率, 明确的活性位点和长期稳定性等特点, 其具有优异的催化性能与巨大的应用前景[39,40]. 而联苯烯单层作为一种新兴的二维碳材料, 其吸附铁原子的性质尚未被研究.

      本文采用基于第一性原理的密度泛函理论方法, 研究了铁原子在联苯烯单层的吸附构型并分析了吸附前后的电子结构. 结果表明, 联苯烯单层是铁原子的良好载体. 对于铁原子的吸附, 联苯烯单层有着比石墨烯更大的吸附能, 二者之间存在明显的电荷转移和稳定的化学键. 铁原子的吸附使联苯烯单层变为磁性金属材料. 这对小型化自旋电子器件及单原子催化剂等领域的研究具有参考价值.

    2.   计算方法

    • 本文计算使用基于第一性原理的VASP 软件包(Vienna ab-inito simulation package)[41,42]. 采用投影缀加平面波方法[43]赝势来描述离子实与电子的相互作用, 交换关联泛函采用Perdew-Burke-Ernzerhof形式[44]的广义梯度近似(generalized gradient approximation, GGA). 采用GGA+U方法和VDW-D3方法考虑了d电子态在位库仑能修正和范德瓦耳斯作用[45,46]. 结构优化和自洽计算的平面波基组展开的截断能设置为500 eV, 结构优化的总能量和力的收敛阈值分别为10–5 eV和0.01 eV/Å. 进行自洽计算时采用了18 × 15 × 1的k点网格, 确保k点间隔在倒空间中小于 2π × 0.005 Å–1. 联苯烯单层位于晶胞的xy平面, 采用联苯烯单层的 3 × 3 × 1超胞构建吸附模型, 在z方向上取18 Å的真空层以模拟分离的联苯烯单层. 用从头算分子动力学(ab-initio molecular dynamics, AIMD)模拟来检验热稳定性. 采用粒子数-体积-温度(NVT)系综, 温度保持在1000 K, 持续5 ps, 时间步长为1 fs[47].

    3.   结果与讨论

      3.1.   吸附模型

    • 如图1(a)所示, 根据对称性, 石墨烯有3个可能的原子吸附位点, 分别是顶位(T)、空位(H)和桥位(B). 结果表明, 铁原子可以稳定吸附于空位(H)上, 吸附高度为1.53 Å, 铁-碳键长为 2.11 Å, 这与已有的结果相吻合[48], 说明本文中计算参数的选取是合理的; 联苯烯单层是一种由四元、六元和八元的碳环组成的单原子层二维材料, 图1(b)的黑色虚线框展示了该材料的最简原胞. 由于对称性的降低, 联苯烯单层需要考虑更多的吸附位点, 分别是2个顶位(T1, T2)、 3个空位(H1, H2, H3)和4个桥位(B1, B2, B3, B4), 在图1(b)中分别由绿、红、蓝点标注. 图1(c)—(h) 展示了铁原子在 H1, H2, 和H3位点的吸附模型, 经过完全的结构弛豫计算发现, 铁原子可以稳定地吸附在联苯烯单层的H1, H2和H3位点. 在H3位点吸附时, 铁原子倾向于吸附在八元环空位的一侧, 且会导致联苯烯单层发生明显的晶格畸变.

      • 3.2.   吸附能

    • 为检查吸附的稳定性, 首先计算吸附能, 计算公式如下:

      式中, Etot 表示载体吸附原子后的总能量; Eslab 表示未作吸附时的原始联苯烯单层的能量; EFe 表示铁的自由原子的能量, 越大的吸附能, 意味着铁原子在能量上越利于吸附. 如表1所列, 计算铁原子在石墨烯和联苯烯单层上的吸附能时发现, 对石墨烯而言, 铁原子在其空位(H)的吸附能, 明显大于在顶位(T)和桥位(B)的吸附能, 即从能量角度看, 石墨烯的空位是铁原子的最稳定吸附位点; 对于联苯烯单层来说, 铁原子在其四元环和六元环空位 (H1, H2和H3) 的吸附能明显较大, H1, H2和H3位点都是铁原子的稳定吸附位点, 其中四元环空位的吸附能最大, 是铁原子的最稳定吸附位点; 联苯烯单层的顶位(T1, T2)和桥位 (B1, B2, B3, B4)的吸附能虽都为正值, 但结构优化表明这些吸附构型并不能稳定存在, 而且, 与H1和H2位点相比, 吸附能明显较小, 即从能量角度看铁原子并不倾向于在这些位点形成稳定的吸附.

      石墨烯H位点对铁原子的吸附能为0.84 eV, 由于吸附能较低, 在实际应用中石墨烯难以作为铁原子的载体. 铁原子在联苯单层 H1和H2位点的吸附能大小分别为1.56和1.53 eV, 几乎是铁原子在石墨烯上吸附能的2倍. 从能量角度看, 联苯烯单层对铁原子的吸附比石墨烯更为稳定, 说明联苯烯单层是过渡金属原子的理想吸附载体, 可作为电、磁、热、催化等功能材料的理想平台.

      • 3.3.   分子动力学

    • 通过第一性原理的AIMD, 进一步检查了联苯烯单层空位(H1, H2和H3)吸附铁原子的热力学稳定性. 在1000 K下, 铁原子吸附于联苯烯单层的四元、六元和八元环后能量随时间的演化曲线见图2, 其中还给出了对应1000 K下经过5 ps后的吸附情况的俯视图和侧视图. 如图2(a)和图2(b)所示, 当铁原子吸附于四元环和六元环空位时, 能量随时间的变化曲线在一定的区间内波动, 且不会出现明显的能量下降, 吸附铁原子后的联苯烯单层的二维结构能在1000 K下很好地保持, 铁原子能够稳定地保持在四、六元环的空位(H1, H2位点). 如图2(c)所示, 当铁原子吸附于联苯烯八元环空位时, 能量曲线随时间的变化曲线在一定区间内, 有着相对较大的起伏, 由于八元环有着较大的空间, 在热扰动下, 铁原子会在八元环内作往复运动, 运动过程如图2(d)中的①, ②, ③, ④和⑤所示. 该吸附特点是石墨烯等材料所不具备的.

      • 3.4.   电荷转移

      3.4.1.   Bader电荷

    • 利用Bader电荷分析[49], 定量计算了铁原子与联苯烯单层(H1, H2和H3)间的电荷转移情况. 当吸附在石墨烯空位时, 每个铁原子失去0.80 e, 最近邻和次近邻碳原子得电子, 每个最近邻碳原子得电子约为0.10 e. 与石墨烯的吸附情况类似, 当铁原子吸附在联苯烯单层时, 电荷从铁原子转移至近邻碳原子上. 在四元环空位吸附时, 铁原子失去0.73 e, 每个最近邻的碳原子得电子为0.14 e, 其余电子转移到次近邻碳原子上; 联苯烯单层中的六元环为压扁的六元环, 从铁原子看, 有四个最近邻和两个次近邻碳原子, 铁原子失去0.83 e, 每个最近邻碳原子得0.23 e, 另外两个次近邻碳原子各失去0.06 e; 在八元环空位时, 铁原子上失去0.69 e, 每个最近邻碳原子得到0.15 e. 这些电荷转移说明铁原子与联苯烯单层有着紧密的结合, 与石墨烯情况相比, 最近邻碳原子上获得的电荷数增多, 所以有着更大的吸附能.

      • 3.4.2.   差分电荷密度分析

    • 通过差分电荷计算, 进一步分析了铁原子在联苯烯单层上的电荷转移情况, 计算方法如下:

      式中, ρtot 表示联苯烯单层吸附铁原子后的总电荷密度, ρslab表示未作吸附时的原始载体的电荷密度, ρFe 表示孤立铁原子的电荷密度.

      如图3所示, 黄色部分表示电荷密度增加, 蓝色表示电荷密度减少. 由图3中黄色区域可以看出, 铁原子和联苯烯单层间出现了电荷的聚集, 这与Bader电荷分析中铁原子对联苯烯单层转移电荷的结果相一致, 表明联苯烯单层与铁原子之间形成了稳定的化学键.

      • 3.5.   电子结构

      3.5.1.   石墨烯和联苯烯单层

    • 石墨烯与联苯烯单层在–4.0—3.0 eV能量范围内的能带与电子态密度如图4(a) 和图4(b)所示. 由费米能级附近的电子态分布可见, 联苯烯单层与石墨烯有着截然不同的电子性质, 石墨烯是零能隙的半导体, 费米能级处的电子态密度为零, 在倒空间的K点存在“狄拉克锥”; 而联苯烯单层有两条能带穿过费米能级, 不再具有“狄拉克锥”的能带特征, 在费米能级处态密度有一定的数值, 呈现金属性. 对联苯烯单层来说, 不用掺杂就有一定的载流子浓度, 可作为导电薄层或条带应用于量子器件.

      • 3.5.2.   铁原子吸附于联苯烯单层四元环空位

    • 以铁原子的最稳定吸附位点为代表, 图5(a)展示了当铁吸附于联苯烯单层的四元环空位时, 铁的3d轨道及其最近邻碳原子2p轨道的投影态密度. 在费米能级以下, 铁的3d轨道和碳的2p轨道有较强的杂化, 结合Bader电荷中吸附的Fe原子向近邻碳原子转移的电荷约为 0.73 e, 进一步说明铁原子与联苯烯单层之间能够形成较强的相互作用. 在费米能级附近 (–3.0—2.0 eV), 体系的电子态主要由铁的3d轨道主导. 同时, 铁3d电子态有较强的自旋极化, 从而形成局域磁矩, 使铁原子与联苯烯单层的复合体系表现出明显的磁性. 此外, 在–0.23 eV处, 自旋向下的电子态存在着极强的态密度尖峰. 图5(c)显示了铁3d轨道的分轨道${\mathrm{d}}_{xy}, $$ {\mathrm{d}}_{yz}, {\mathrm{d}}_{xz}, {\mathrm{d}}_{{x}^{2}-{y}^{2}}, {\mathrm{d}}_{{z}^{2}}$, 由图5(c)可以看出, 该尖锐的态密度峰来源于铁的$ {\mathrm{d}}_{{z}^{2}} $轨道. 这是由于垂直于联苯烯单层的$ {\mathrm{d}}_{{z}^{2}} $轨道难以与近邻碳原子的p轨道杂化造成的, 与体系的准二维特性密切相关.

      自旋极化计算表明, 吸附铁原子后的联苯烯单层为磁性金属材料. 当铁原子吸附于联苯烯单层四元环空位时, 磁矩大小约为 1.81 ${\mu }_{\mathrm{B}}$. 进一步检查材料的非共线磁性发现, 当铁原子吸附于联苯烯单层四元环空位时, 易磁化轴指向面外, 铁原子的磁各向异性能大小为0.64 meV, 与常见的二维磁性材料 CrXTe3 和Fe2Si相当[50,51].

      • 3.6.   磁　序

    • 石墨烯在比表面积、机械强度、电子迁移速度、热导率等方面有着明显的优势, 但石墨烯本身是非磁性的, 限制了石墨烯在自旋电子器件方面的应用. 自石墨烯发现以来, 研究人员采用缺陷、吸附、掺杂、边界等方式实现石墨烯基复合材料的磁性[52]. 由于石墨烯的吸附位点(六元环空位)太过密集且吸附能较小, 不利于形成规则的磁原子阵列, 又因为联苯烯单层的四元环空位形成分布稀疏、间距合理的吸附位点网格, 且吸附牢固, 所以, 以联苯烯为基底, 四元环空位吸附上铁原子, 易于形成稳定的磁有序结构. 这是联苯烯相对于石墨烯一个明显的优点. 本研究小组设计了几种磁序, 如图6所示, 分别命名为铁磁(FM)、共线反铁磁序一(Coll-I)、共线反铁磁序二(Coll-II)和奈尔反铁磁序(Nèel), 用标准GGA方法分别计算了相应的能量, 并显示在表2第1行中. 这些数据显示奈尔反铁磁序最稳定, 从而确定均匀吸附铁原子的联苯烯单层具有反铁磁性.

      • 3.7.   用GGA + U方法考虑电子关联修正后的计算结果

    • GGA + U方法能够用直接增加d轨道在位库仑排斥能的方式修正d电子的关联效应, 比杂化泛函方法能得到更清晰的物理图像, 又采用GGA + U方法计算了联苯烯吸附铁原子的吸附位点、电子结构及磁性, 以进一步核实前面的计算结果. 首先, 通过线性响应方法自洽地得到了联苯烯吸附铁原子体系的U值[45]. 如图7所示, 蓝线和红线分别为自洽和非自洽情形下d轨道电子数对微小Hubbard U值的线性响应, 斜率为k1和k2, 则U = 1/k1 –1/k2 = 3.6 eV. 另外, 吸附体系中易于产生电极化和瞬态偶极相互作用, 色散力对吸附构型起着重要的作用, 因而在计算中采用DFT-D2形式来引入范德瓦耳斯修正[46]. 在考虑电子关联和色散关联修正情况下, 对图1(b)中联苯烯的各个吸附位点优化了铁原子吸附的几何结构, 并用(1)式计算了吸附能. 结果表明, 空位吸附比桥位和定位吸附更稳定. 空位吸附时, 四环、六环、八环空位的吸附能分别为1.45, 1.21, 1.19 eV, 故四环空位是最稳定的吸附位置, 与前面计算的结论一致. 图5(b)给出了GGA + U计算的四环空位的铁原子投影态密度, 由图5(b)可以看出, 自旋向上和向下的d电子态发生更大的劈裂, 自旋向上d电子态处于更深的能量范围(–6.0 eV至–4.0 eV). 由图5(d)可以看出, 自旋向下的 $ {\mathrm{d}}_{{z}^{2}} $电子态移到了费米能之上, 这将导致更大磁矩出现, 铁原子的磁矩为 4.87 ${\mu }_{\rm B}$. 另外, 对于图6中4种磁序, 也计算了相应的能量, 列于表2中. 考虑关联效应后, 奈尔反铁磁序的能量仍然最低, 进一步确认了其铁磁奈尔序基态.

    4.   总结与展望

    • 本文通过密度泛函理论框架下的第一性原理计算研究了联苯烯单层对铁原子的吸附情况, 明确了联苯烯单层可能的吸附位点. 通过吸附能和分子动力学的计算, 证明了铁原子可以稳定的吸附于联苯烯单层四元环和六元环的空位上. 铁原子与联苯烯单层存在明显的电荷转移, 铁3d电子态与碳2p电子态有较强的轨道杂化, 二者间能够形成稳定的化学键. 铁原子在四元环空位吸附时, 易磁化轴指向面外, 磁矩大小约为1.81 ${\mu}_{\mathrm{B}}$. 相对于石墨烯而言, 联苯烯是更优异的吸附载体, 体现在两个方面, 一是吸附能较大, 二是吸附位点(四环空位)分散且间距合理, 这导致吸附后能形成稳定的复合材料. 联苯烯作为磁性材料, 而且磁性原子排列规则, 便于形成稳定的磁序; 作为催化材料, 催化位点稳固, 分散合理. 联苯烯的另一个优点是金属性, 优良的导电性有利于作为分子导线、电极材料等电输运器件. 金属性也有利于联苯烯吸附铁原子作为催化材料, 有实验研究表明[53], 金属态的电催化剂具有更为优异的电催化活性, 其原因是电子在金属态的催化剂内部可以更快地传输. 因此, 本文结果确认了联苯烯单层是比石墨烯更好的铁原子吸附载体且体系有磁性, 为研究新型吸附材料的电磁、输运、催化等特性提供了重要的基础信息.

    参考文献 (53)

目录

/

返回文章
返回