搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

中国的表面物理

管丹丹 贾金锋

杨海林, 陈琦丽, 顾星, 林宁. 氧原子在氟化石墨烯上扩散的第一性原理计算. 物理学报, 2023, 72(1): 016801. doi: 10.7498/aps.72.20221630
引用本文: 杨海林, 陈琦丽, 顾星, 林宁. 氧原子在氟化石墨烯上扩散的第一性原理计算. 物理学报, 2023, 72(1): 016801. doi: 10.7498/aps.72.20221630
Yang Hai-Lin, Chen Qi-Li, Gu Xing, Lin Ning. First-principles calculations of O-atom diffusion on fluorinated graphene. Acta Phys. Sin., 2023, 72(1): 016801. doi: 10.7498/aps.72.20221630
Citation: Yang Hai-Lin, Chen Qi-Li, Gu Xing, Lin Ning. First-principles calculations of O-atom diffusion on fluorinated graphene. Acta Phys. Sin., 2023, 72(1): 016801. doi: 10.7498/aps.72.20221630

中国的表面物理

管丹丹, 贾金锋

Surface physics in China

Guan Dan-Dan, Jia Jin-Feng
Article Text (iFLYTEK Translation)
PDF
HTML
导出引用
  • 以中国科学院表面物理国家重点实验室为主线, 回顾了表面物理在中国的发展历程, 致敬为中国表面物理做出贡献的老一辈科学家. 通过回顾历史可以看出, 中国表面物理蓬勃发展, 不仅做出了很多国际先进水平的工作, 也培养了大批青年人才, 他们已经成为国际凝聚态物理研究的重要力量.
    In this paper, the development of surface physics in China is comprehensively reviewed, focusing on the State Key Laboratory of Surface Physics at the Chinese Academy of Sciences. It especially recognizes and honors the invaluable contributions made by the older generation of scientists in this field. By looking back at the history, it can be seen that the surface physics has developed vigorously in China: not only have many research papers with international advanced level been published, but also a large number of young talents have been cultivated, who have become an important force in the research of condensed matter physics internationally.
      PACS:
      68.43.Jk(Diffusion of adsorbates, kinetics of coarsening and aggregation)
      31.15.es(Applications of density-functional theory (e.g., to electronic structure and stability; defect formation; dielectric properties, susceptibilities; viscoelastic coefficients; Rydberg transition frequencies))
      73.20.At(Surface states, band structure, electron density of states)
      68.65.Pq(Graphene films)
      通信作者: 贾金锋, jfjia@sjtu.edu.cn
    • 基金项目: 国家重点研发计划(批准号: 2019YFA0308600)、科技创新重大项目(批准号: XDB28000000)和上海市科技重大专项(批准号: 2019SHZDZX01)资助的课题.
      Corresponding author: Jia Jin-Feng, jfjia@sjtu.edu.cn
    • Funds: Project supported by the National Key Research and Development Program of China (Grant No. 2019YFA0308600), the Strategic Priority Research Program of Chinese Academy of Sciences (Grant No. XDB28000000), and the Science and Technology Commission of Shanghai Municipality, China (Grants No. 2019SHZDZX01).

    金属的腐蚀就是金属和腐蚀性介质(如O2, Cl等)发生氧化还原反应而发生有害的物理化学等变化, 失去原有性能的过程. 金属腐蚀是一种非常常见的现象, 几乎遍及日常生活中所有产业. 2020年中国工程院调查表明, 我国年腐蚀成本高达3万亿元, 占全国GDP的3.34%.

    石墨烯(G)是目前最薄的防腐蚀涂层和填料[1-5], 其几何孔径小, 使得原子、分子都难以通过. 2008年, Bunch等[3]研究了气体分子对G的渗透, 发现G能阻隔所有气体分子. 有研究[6]表明, O2分子、H2O分子穿过G分别需要克服16.34 eV, 42.8 eV的能垒. Topsakal等[7]发现O原子穿过G需要克服5.98 eV的能垒, O原子穿透G的扩散率约为表面扩散率的10–87倍, 可以忽略. G还具有优异的机械性能(杨氏模量约为1.0 TPa, 断裂强度约为130 GPa)[1,2,8]和优异的耐磨性. 单层和多层G薄膜都非常透明(4层G的透射率 > 90%), 因此G涂层不会干扰下层金属的光学特性[9,10].

    但是G也具有优异的导电性(电子迁移率为105—106 cm2/(V·s))[2], 在缺陷和破损处会发生电化学腐蚀, 从而加速金属的腐蚀速度[11-13]. 解决这一问题的方法包括功能化石墨烯[14-18]、多层石墨烯[19]、制备自修复石墨烯涂层[10]和掺杂原子[3]等. 功能化石墨烯是在G表面引入一些特定的官能团(如-F, -H, -O等[20]), 当这些官能团吸附在G表面时, 可以打开G的带隙, 降低其导电性能, 从而抑制电化学腐蚀的发生. 以氟化石墨烯(FG)为例, F原子的吸附使G上C原子的sp2杂化变为sp3杂化, 破坏了G的π键, 达到降低G的导电性, 抑制缺陷处的电化学腐蚀的目的[21]. Shen等[22]采用F2在还原氧化石墨烯上吸附F原子的办法, 制备了氟吸附浓度(氟化浓度)分别为12.56%, 33.98%和43.36%的3种氟化还原氧化石墨烯, 发现随着氟化浓度的增打, 氟化石墨烯的导电率降低. 此外, 还有研究表明, 氟功能化石墨烯材料由于氟的存在, 会导致表面能降低, 在水环境中表现出高度疏水性[23]. 高度疏水性有利于潮湿环境下抑制腐蚀性介质的吸附, 可以用于制备超疏水防腐蚀涂层[21]. 同时, FG继承了G的片状结构, 这种结构被认为产生了“迷宫效应”, 扩展了腐蚀性介质的扩散路径, 有利于涂层的抗腐蚀性能的提高 [18].

    不仅如此, 氟功能化的石墨烯还继承了石墨烯薄膜优异的阻隔性能. Wu等[24]通过实验证实了F元素的引入可以限制腐蚀物种的扩散, 延缓腐蚀产物的形成; 并通过理论计算发现其机制在于: 当G表面吸附F原子后, 石墨烯层对腐蚀性介质, 如Cl原子的吸附更强, 同时Cl原子穿过缺陷的能垒也显著提高. Li等[20]使用密度泛函理论分别计算了O原子在单个-O, -OH, -H和-F吸附的G涂层上的扩散和穿透, 结果表明, 这些吸附介质打破了石墨烯的电子结构, 增强了吸附位置的电子云, 提高了O原子垂直穿透G的能垒, 从而增加了O原子通过石墨烯的阻碍; 但同时O原子在水平方向上的扩散迁移率也有所提高, 这有利于O原子对金属基底的腐蚀. Yao等[25]通过第一性原理计算分析了O原子在G、两个Cl原子吸附的G和单个Cl, F原子共吸附的G表面的扩散, 结果发现Cl原子的引入会降低O原子在G上的扩散能垒, 因此当环境中的Cl原子吸附在G表面时, 会导致O原子的扩散速率增大, 而F原子吸附在Cl原子的对顶位, 增加了O原子的扩散能垒, 能有效抑制O原子的扩散; 该研究说明氟化对含氯环境下的腐蚀防护有利.

    以上研究表明降低腐蚀性物质在水平方向上的扩散迁移率是功能化石墨烯的抗腐蚀研究上需要关注的一个问题. 但以上理论研究都是围绕单个F原子吸附的FG(低于5%)展开的. 对于氟化浓度较高的石墨烯, 如完全氟化石墨烯(CF)和C∶F为4∶1的单面吸附的部分氟化石墨烯(C4F), 它们都是最稳定氟化石墨烯结构[26], 其宽带隙为3 eV左右[21,27-29] , 而且CF的C原子都是sp3C原子, 与纯粹石墨烯相比, 电子结构变化最大; 而C4F上的F原子呈对顶位吸附, 可以将其表面分隔开, 可推测这能抑制O原子的表面扩散. 因此本文采用基于密度泛函理论的第一性原理计算方法, 研究了O原子在CF和C4F上的扩散和穿透, 探讨CF和C4F对腐蚀性物质的阻隔机理. 此外, 与基底界面的结合强度也是涂层抗腐蚀性能的重要指标之一, 本文也对CF和C4F与Cu基底的相互作用进行了研究.

    本工作使用的是基于密度泛函理论的第一性原理计算软件VASP (vienna ab-initio simulation package)[30], 使用投影缀加平面波(PAW)方法[31]描述电子-离子相互作用, 利用广义梯度近似(GGA)PW91泛函交换关联相互作用[32]. 采用500 eV的平面波截断能, 二维结构和表面的k点网格设置为3 × 3 × 1, 真空层厚度为15 Å, 所有几何结构都进行了全优化, 直到作用在每个原子上的能量和力分别小于10–5 eV和10–2 eV/Å. 还使用DFT-D2方案[33]来修正范德瓦耳斯(vdW)力.

    使用VASP软件的NEB过渡态搜索方法[34]来找到O原子在CF和C4F上的扩散和穿透的路径和能垒. NEB方法只需要知道初态和末态的结构, 就能得到扩散的路径和能垒.

    为了评估FG涂层与基底材料的结合能力, 本文计算了FG涂层与Cu(111)表面的界面结合能, 界面结合能(Eab)的计算公式[35]如下:

    Eab=ECu+ElECu/l,
    (1)

    其中ECu/l表示结合后界面总能量, ECu表示Cu(111)表面能量, El表示吸附的G和FG(包括C4F和CF)层能量, 由于各界面的面积不同, 可以用界面粘附功(Wab)来对比不同界面的结合强度以及界面的稳定性, 粘附功越大, 界面结合力越强, 界面结构越稳定. 其计算公式如下[35]:

    Wab=Eab/ACu/l,
    (2)

    其中ACu/l为界面接触面积.

    首先构建并优化了G, C4F和CF的单层模型, 结构如图1所示. CF中的F原子交替吸附在C原子层两侧, C—F键为1.383 Å, C—C键为1.577 Å, 与文献[36]的1.38 Å和1.58 Å一致. C4F的F原子吸附在C原子层同一侧, 且呈对顶位分布, C—F键长为1.46 Å, 接近文献的1.45 Å, C—C键长为1.508 Å和1.401 Å, 与文献[37]的1.51 Å(Csp3-Csp2)和1.40 Å (Csp2-Csp2)一致.

    图 1 优化后的晶格结构 (a) G; (b) C4F; (c) CF\r\nFig. 1. Optimized lattice structure: (a) G; (b) C4F; (c) CF.
    图 1  优化后的晶格结构 (a) G; (b) C4F; (c) CF
    Fig. 1.  Optimized lattice structure: (a) G; (b) C4F; (c) CF.

    FG与Cu(111)界面模型如图2所示. 分别选取了4 × 4的单层G (C32, a = b = 9.86 Å), C4F (C32F8, a = b = 9.93 Å)和CF (C32F32, a = b = 10.39 Å)与2 × 2的3层Cu(111)表面结构(Cu48, a = b = 10.26 Å)结合, 晶格失配率分别为3.97%, 3.27%和1.26%, 晶格失配率公式[38]如下:

    图 2 优化前的界面模型 (a) Cu/G界面; (b) Cu/C4F界面; (c) Cu/CF界面\r\nFig. 2. Interface model before optimization: (a) Cu/G interface; (b) Cu/C4F interface; (c) Cu/CF interface.
    图 2  优化前的界面模型 (a) Cu/G界面; (b) Cu/C4F界面; (c) Cu/CF界面
    Fig. 2.  Interface model before optimization: (a) Cu/G interface; (b) Cu/C4F interface; (c) Cu/CF interface.
    δ=2|ab|a+b.
    (3)

    研究表明, 桥位是O原子在G上的最佳吸附位[20,7]. 首先计算了O原子在G相邻两个桥位上的扩散, 得到它的最佳扩散路径是桥位(B)-顶位(T)-桥位(B), 与文献[20, 7]一致. 如图3所示, 它的扩散能垒为0.72 eV, 接近文献[20]的0.79 eV. 说明计算方法可取.

    图 3 O原子在G相邻桥位上的扩散的能垒图\r\nFig. 3. The diffusion barrier diagram of diffusion of O atom on the adjacent bridge site of G.
    图 3  O原子在G相邻桥位上的扩散的能垒图
    Fig. 3.  The diffusion barrier diagram of diffusion of O atom on the adjacent bridge site of G.

    分别计算O原子在C4F 无F面和有F面上的扩散, 首先, 研究发现sp2碳原子(不吸附F的C原子)间的桥位是O原子的最稳定吸附位(见附表A1). 图4(a)展示了O原子在C4F无F面的四条扩散路径, 分别是路径1: B-T-B; 路径2: B-H1(空心位1)-B; 路径3: B-H2(空心位2)-B; 路径4: B-T-T-B. 图4(b)展示了O原子在C4F有F面的3条扩散路径, 分别是路径5: B-T-B; 路径6: B-H1-B; 路径7: B-H2-B.

    图 4 (a) O原子在C4F无F面的扩散路径; (b) O原子在C4F有F面的扩散路径\r\nFig. 4. (a) the diffusion paths of O atom on the F-free surface of C4F; (b) the diffusion paths of O atom on the F adsorbed surface of C4F.
    图 4  (a) O原子在C4F无F面的扩散路径; (b) O原子在C4F有F面的扩散路径
    Fig. 4.  (a) the diffusion paths of O atom on the F-free surface of C4F; (b) the diffusion paths of O atom on the F adsorbed surface of C4F.

    在C4F无F面上, O原子的最佳扩散路径是B-T-B, 即路径1, 如图5(a)所示, 其能垒为1.44 eV, 相比于在G上的扩散能垒提高了近一倍. O原子沿路径2(通过空心位H1)扩散到sp2碳六边形对桥位的能垒为4.31 eV, 沿路径3(通过空心位H2)扩散到不同sp2碳六边形环上的能垒为3.42 eV, 远高于沿顶位扩散的能垒, O原子难以通过空心位扩散. 在C4F无F面上, O原子还可以通过路径4扩散到不同sp2碳六边形环上, 其能垒如图5(d)所示, 为2.66 eV, 小于沿空心位扩散的能垒, 但相较于路径1而言, 扩散难度大大提升.

    图 5 O原子在C4F无F面上扩散的能垒图 (a) 路径1; (b) 路径2; (c) 路径3; (d) 路径4\r\nFig. 5. Diffusion barrier diagram of O atom diffusion on the F-free surface of C4F: (a) Path 1; (b) path 2; (c) path 3; (d) path 4.
    图 5  O原子在C4F无F面上扩散的能垒图 (a) 路径1; (b) 路径2; (c) 路径3; (d) 路径4
    Fig. 5.  Diffusion barrier diagram of O atom diffusion on the F-free surface of C4F: (a) Path 1; (b) path 2; (c) path 3; (d) path 4.

    O原子在C4F有F面上的最佳扩散路径也是B-T-B, 即路径5, 如图6(a)所示, 能垒为1.42 eV, 与无F面相近. 图6(b)显示了O原子沿路径6的扩散, 相对于沿顶位的扩散, 沿空心位的能垒更高, 达到了3.36 eV, O原子难以沿空心位扩散. 在不同sp2碳六边形间的扩散, 由于F原子的阻隔, O原子只能通过空心位H2扩散, 如图6(c)所示, 沿路径7的能垒高达4.0 eV, 如此高的能垒将O原子困在sp2碳六边形环上, 几乎无法扩散出F原子的包围, 这有利于涂层的腐蚀防护.

    图 6 O原子在C4F有F面上扩散的能垒图 (a) 路径5; (b) 路径6; (c) 路径7\r\nFig. 6. Diffusion barrier diagram of O atom diffusion on the F adsorbed surface of C4F: (a) Path 5; (b) path 6; (c) path 7
    图 6  O原子在C4F有F面上扩散的能垒图 (a) 路径5; (b) 路径6; (c) 路径7
    Fig. 6.  Diffusion barrier diagram of O atom diffusion on the F adsorbed surface of C4F: (a) Path 5; (b) path 6; (c) path 7

    O原子垂直穿透G的最佳路径是从桥位穿透[4,7], 如图7(a)所示, 计算得到该路径的能垒为6.16 eV, 与文献[4]的5.5 eV和文献[7]的5.98 eV相差不大. 对于O原子穿透C4F的路径, 同样计算了从桥位穿透的情况, 得到能垒为5.08 eV(见图7(b)), 略小于穿透G的能垒, 但仍然能有效阻隔O原子的穿透. 虽然F的加入破坏了C4F上部分C原子的sp2结构, 但其保留的sp2C原子环的结构仍然很稳定, 因此能有效阻隔O原子的穿透.

    图 7 O原子穿透 (a) G; (b) C4F的能垒图\r\nFig. 7. Diffusion barrier diagram of O atom penetrating: (a) G; (b) C4F.
    图 7  O原子穿透 (a) G; (b) C4F的能垒图
    Fig. 7.  Diffusion barrier diagram of O atom penetrating: (a) G; (b) C4F.

    以上计算表明C4F上的F原子大大提高了O原子的水平扩散难度, 能有效抑制O原子在其表面的扩散, 有利于阻碍O原子对金属基底的腐蚀. 这与文献[20]中O在单个F原子修饰的G表上扩散能垒降低的情况不同, 而与文献[25]中Cl原子和F原子对顶位吸附在G中时, O原子的扩散能垒增加相似, 说明C4F上F原子的对顶位吸附应该是O原子的扩散能垒提升的主要原因, 有可能是对顶位的F在石墨烯上形成的稳定共轭电子, 从而增大了O原子的扩散能垒.

    C4F中只有单面吸附F原子, 但无论是有F面还是无F面, O原子都倾向于吸附在sp2 C原子环的桥位上, 并沿着桥位-顶位-桥位扩散. F原子的环绕排列, 大大提高了O原子向周边扩散的能垒, F原子将C4F分隔成一个个区块, 每个区块像陷阱一样, 当O原子吸附在上面时, 巨大的能垒使其难以向下或向其他区块扩散, 因此能很好阻止O原子扩散到金属基底.

    在CF上, 所有的C原子都吸附一个F原子, 变为了sp3C原子, C原子间距离变大. O原子需要穿过表面F原子层进入C原子层, 首先讨论了O原子穿透CF的情形. 如图8(a)所示, 首先, 当O原子接近F原子层时, 受到F原子的排斥, 要克服2.82 eV的能垒进入F原子层并替换F原子, 吸附在C原子顶位上, 之后向下运动到C原子间桥位上, 该过程需克服1.84 eV的能垒, 之后O原子再穿过桥位, 需要克服0.81 eV的能垒. 可见, 相比于G, CF的阻隔能力下降, O原子穿透CF的阻力主要来自于F原子层; 但O原子穿透CF的路径更复杂, 对于减缓O原子的穿透有一定作用.

    图 8 (a) O原子穿透CF的F原子层到达C原子层的能垒图; (b) O原子在CF内相邻桥位上的扩散的能垒图; (c) O原子在CF内沿空心位扩散的能垒图\r\nFig. 8. (a) Diffusion barrier diagram of O atom penetrating the F atomic layer of CF to the C atomic layer; (b) diffusion barrier diagram of O atom diffusion on the adjacent bridge site in CF; (c) diffusion barrier diagram of O atom diffusion along the hollow site in CF.
    图 8  (a) O原子穿透CF的F原子层到达C原子层的能垒图; (b) O原子在CF内相邻桥位上的扩散的能垒图; (c) O原子在CF内沿空心位扩散的能垒图
    Fig. 8.  (a) Diffusion barrier diagram of O atom penetrating the F atomic layer of CF to the C atomic layer; (b) diffusion barrier diagram of O atom diffusion on the adjacent bridge site in CF; (c) diffusion barrier diagram of O atom diffusion along the hollow site in CF.

    O原子穿过F原子层后, O原子取代C—C键连接在两个相邻C原子间. 我们计算了O原子在CF层内的水平扩散, 其中沿空心位扩散需要克服7.60 eV的巨大能垒(见图8(c)), 相比之下, 在相邻桥位间扩散更容易. 其最佳扩散路径的能垒如图8(b)所示, O原子扩散到相邻桥位过程中, 先与两个C原子成键连接, 移动到中间时, 与3个C原子成键连接, 这个过程吸热, 能垒为2.69 eV, 之后再移动到相邻桥位上, 与两个C原子成键连接, 如此高的能垒能有效防止O原子向周围扩散.

    为了确定F吸附对O扩散的影响, 对表面扩散的几条最佳扩散路径过渡态的电子结构进行了分析, 图9图10显示了O原子在石墨烯和F化石墨烯上沿最佳路径扩散的初始态及过渡态的态密度. 对比初始态和过渡态的态密度可以发现, 在各个过渡态的费米能级附近都会出现新的杂化峰, 这些杂化峰由C的p轨道和O的p轨道杂化产生. 从初始态到过渡态过程中, 在G和C4F表面, 过渡态的O原子处于顶位, 其中一个C—O键断裂, 体系能量升高, 而在CF内断裂一个C—C键, 生成一个C—O键. O原子吸附在C4F和CF上时, O原子的p轨道与F原子的p轨道在–5 eV附近有重叠峰, 到达过渡态时, 峰值减小, 扩散过程中, 受到F原子的影响, O原子需要更高的能量完成迁移.

    图 9 O原子沿最佳扩散路径扩散的DOS图 (a) G表面初始态和过渡态的DOS图; (b) C4F有F面初始态和过渡态的DOS图; (c) C4F无F面初始态和过渡态的DOS图; (d) CF内初始态和过渡态的DOS图\r\nFig. 9. DOS diagram of O atom diffusion along the optimal diffusion path: (a) DOS diagram of initial state and transition state of G surface; (b) DOS diagram of initial state and transition state on the F adsorbed surface of C4F; (c) DOS diagram of initial state and transition state on the F-free surface of C4F; (d) DOS diagram of initial state and transition state in CF.
    图 9  O原子沿最佳扩散路径扩散的DOS图 (a) G表面初始态和过渡态的DOS图; (b) C4F有F面初始态和过渡态的DOS图; (c) C4F无F面初始态和过渡态的DOS图; (d) CF内初始态和过渡态的DOS图
    Fig. 9.  DOS diagram of O atom diffusion along the optimal diffusion path: (a) DOS diagram of initial state and transition state of G surface; (b) DOS diagram of initial state and transition state on the F adsorbed surface of C4F; (c) DOS diagram of initial state and transition state on the F-free surface of C4F; (d) DOS diagram of initial state and transition state in CF.

    除了涂层的阻隔性能, 还研究了FG与金属Cu(111)面间的相互作用, 并与G与Cu(111)面的相互作用进行比较, 这对于评估涂层材料的使用价值有很大意义. 本文使用的是4 × 4的G, G与Cu(111)面的界面粘附功为2.626 J/m2, 界面距离为2.938 Å (见图10(a)), 根据bader电荷分析, 发现G与Cu表面第一层原子发生电荷交换, 平均单个Cu原子向G层转移0.03个电荷, G与Cu表面通过弱键相连. C4F与Cu(111)面的界面粘附功更高, 粘附功为3.529 J/m2, C4F与Cu基底表面距离为2.045 Å (见图10(b)), 可以看到部分C原子与表面Cu原子间形成较强的键合, bader电荷分析发现, 表面Cu原子与C4F的C原子间的电荷交换增加, 平均单个Cu原子失去0.092个电荷, 界面粘附功的增加与F原子没有直接关系, 主要是F的吸附破坏了石墨烯的电子结构, 使得C原子与表面的Cu原子能形成更强的键合, 这使得C4F与Cu基底的结合更为紧密, 有利于金属的腐蚀防护. CF与Cu(111)面的界面粘附功为3.559 J/m2, 略高于C4F与Cu(111)面的界面粘附功, 底层F原子与表层Cu原子间距离为2.827 Å (见图10(c)). 根据bader电荷分析, 最上层Cu原子平均失去0.042个电荷, 相比Cu/G界面略微增加, 下层F原子与表面Cu原子间形成弱键.

    图11展示了3种界面体系在费米面附近的总态密度以及3种原子(包括Cu原子的3d轨道、C原子的2p轨道和F原子的2p轨道)的分波态密度, 与Cu/G界面一样, Cu/C4F界面的 C原子在费米能级处态密度不为0, 说明C4F与Cu结合后, 复合体系表现出金属性, 对电化学腐蚀 有利.

    图 10 优化后的界面模型 (a) Cu/G界面; (b) Cu/C4F界面; (c) Cu/CF界面\r\nFig. 10. Optimized interface model: (a) Cu/G nterface; (b) Cu/C4F interface; (c) Cu/CF interface.
    图 10  优化后的界面模型 (a) Cu/G界面; (b) Cu/C4F界面; (c) Cu/CF界面
    Fig. 10.  Optimized interface model: (a) Cu/G nterface; (b) Cu/C4F interface; (c) Cu/CF interface.
    图 11 Cu/G、Cu/C4F和Cu/CF界面的DOS图, 虚线为费米能级\r\nFig. 11. DOS diagram of Cu/G, Cu/C4F and Cu/CF interfaces.
    图 11  Cu/G、Cu/C4F和Cu/CF界面的DOS图, 虚线为费米能级
    Fig. 11.  DOS diagram of Cu/G, Cu/C4F and Cu/CF interfaces.

    由Cu/G界面的态密度图可以发现, G中C的p轨道与Cu的d轨道共振峰较少, 只在费米能级附近出现, 而且展宽较小, 说明G与Cu基底间相互作用不强. 对于Cu/C4F界面, Cu原子的d轨道主要在0到5 eV以内与C原子的p轨道发生杂化, 共振峰的波峰较为平缓, 重叠区域面积较大, 电子离域性较强, 说明Cu原子和C原子间成键较强, C4F与Cu表面的相互作用较强. 对于Cu/CF界面, F原子和C原子的p轨道在0到5 eV和–3 eV以下有多处共振峰和重叠区域, 说明CF比C4F更稳定, CF在费米能级处态密度为0, 与Cu的结合对CF的电子结构没有受到太大影响. 值得注意的是, 在–5到5 eV范围内, 相比于Cu/G和Cu/C4F界面, Cu/CF界面上Cu原子的d轨道与F原子和C原子的p轨道有更多的重叠区域, 尖峰数量较多, 而且在–5—–3 eV范围内Cu的d轨道与F的p轨道有明显的杂化峰, 且重叠区域面积较大, 说明Cu/CF界面中Cu基底与F原子间的杂化作用更强, 因此Cu/CF界面的粘附功大于Cu/G和Cu/C4F界面, 所以Cu/CF界面结构更稳定. 界面粘附功的提升可能是因为F原子具有高负电性, 与表面Cu原子间的相互作用更强.

    为了更好比较界面原子间的键合强度, 使用Lobster程序[39]计算了表面Cu原子与界面处最近邻的C原子或F原子间的COHP和ICOHP, 如图12展示了界面处原子对的-COHP, -COHP为正表示成键作用, 为负表示反键作用. Cu/G界面的Cu—C键、Cu/C4F界面的Cu—C键和Cu/CF界面的Cu—F键的-ICOHP值分别为0.195, 0.395和0.032, 表明界面间键强由大到小为Cu/C4F界面的Cu—C键 > Cu/G界面的Cu—C键 > Cu/CF界面的Cu—F键. 由图12可以发现Cu/CF界面的Cu—F键在费米能级以下有很大的反键区域, 这使得Cu/CF界面的电荷转移虽高于Cu/G界面, 但Cu/CF界面的Cu—F键的强度却更小.

    图 12 3种界面的轨道哈密顿分布(-COHP), 费米能级位于0 eV\r\nFig. 12. Crystal orbital Hamilton populations (-COHP) of three interfaces, the Fermi level is at 0 eV
    图 12  3种界面的轨道哈密顿分布(-COHP), 费米能级位于0 eV
    Fig. 12.  Crystal orbital Hamilton populations (-COHP) of three interfaces, the Fermi level is at 0 eV

    本文采用第一性原理计算方法, 研究了O原子在CF和C4F上的扩散和穿透行为. 计算结果表明, 在C4F上, O原子的扩散受到了F原子的阻碍, O原子吸附在相邻sp2C原子间的桥位上, 最佳扩散路径是沿sp2C原子环上的B-T-B, 需要克服1.42 eV的能垒. O原子在CF内的扩散能垒为2.69 eV, 比在G和C4F上的扩散能垒大很多, 扩散难度大.

    研究还发现随着F原子吸附浓度的增大, 石墨烯C原子层的稳定结构被逐渐破坏, C原子层的阻隔能力逐渐减小, CF的C原子层对O原子的阻隔能力大大减小, O原子主要受到F原子层的阻隔, C4F的C原子层的结构没有受到太大的破坏, 依旧有很好的阻隔能力.

    CF和C4F薄膜与Cu基底界面的相互作用研究表明, 石墨烯的F化可以提高薄膜与金属基底的结合强度, 且界面结合强度与稳定性随F化浓度的增大而变强. 3种界面体系中, Cu/CF界面的界面粘附功最大, 说明Cu/CF界面的结构最稳定, Cu原子的电荷转移比Cu/G界面高一些, 相比于Cu/C4F界面则更少, 成键不强, 界面粘附功主要由F原子和Cu原子间的电荷相互作用提供; Cu/C4F界面的界面粘附功略小于Cu/CF界面的界面粘附功, 其界面结构较为稳定, Cu原子的电荷转移最多, 部分Cu原子和C原子间成键较强, C4F上F的吸附导致石墨烯电子结构被破坏, 使得C原子与Cu表面原子间形成强键, 界面粘附功增大; Cu/G界面的界面粘附功最小, 电荷转移也最少, 说明G与Cu基底的结合最不稳定.

    表 A1  O原子在C4F不同位置吸附的能量数据
    Table A1.  Energy data of O atom adsorption at different positions of C4F.
    O原子吸附位置C4F
    有F面/eV无F面/eV
    桥位–325.31–326.67
    顶位–323.89–325.23
    空心位H1–321.96–322.36
    空心位H2–321.92–323.25
    下载: 导出CSV 
    | 显示表格
    表 A2  O原子在CF内不同位置吸附的能量数据
    Table A2.  Energy data of O atom adsorption at different positions in CF.
    O原子吸附位置CF内/eV
    桥位–413.70
    空心位–406.10
    下载: 导出CSV 
    | 显示表格
    [1]

    Yang W S, Jona F, Marcus P M 1983 Phys. Rev. B 27 1394Google Scholar

    [2]

    Yang W S, Jona F, Marcus P M 1983 Phys. Rev. B 28 2049Google Scholar

    [3]

    Jia J F, Zhao R G, and Yang W S 1993 Phys. Rev. B 48 18109Google Scholar

    [4]

    Jia J F, Zhao R G, Yang W S 1993 Phys. Rev. B 48 18101Google Scholar

    [5]

    贾金峰, 赵汝光, 杨威生 1992 物理学报 41 827Google Scholar

    Jia J F, Zhao R G, Yang W S 1992 Acta Phys. Sin. 41 827Google Scholar

    [6]

    Li J L, Jia J F, Liang X J, Liu X, Wang J Z, Xue Q K, Li Z Q, Tse J S, Zhang Z Y, Zhang S B 2002 Phys. Rev. Lett. 88 066101Google Scholar

    [7]

    Jia J F, Wang J Z, Liu X, Xue Q K, Li Z Q, Kawazoe Y, Zhang S B 2002 Appl. Phys. Lett. 80 3186Google Scholar

    [8]

    Jia J F, Liu X, Wang J Z, Li J L, Wang X S, Xue Q K, Li Z Q, Zhang Z Y, Zhang S B 2002 Phys. Rev. B 66 165412Google Scholar

    [9]

    Xu Z, Bai X D, Wang Z L, Wang E G, 2006 J. Am. Chem. Soc. 128 1052Google Scholar

    [10]

    Xu Z, Bai X D, Wang E G 2006 Appl. Phys. Lett. 88 133107Google Scholar

    [11]

    Wang W L, Bai X D, Liu K H, Xu Z, Golberg D, Bando Y, Wang E G 2006 J. Am. Chem. Soc. 128 6530Google Scholar

    [12]

    Guo Y, Zhang Y F, Bao X Y, Han T Z, Tang Z, Zhang L X, Zhu W G, Wang E G, Niu Q, Qiu Z Q, Jia J F, Zhao Z X, Xue Q K 2004 Science 306 1915Google Scholar

    [13]

    Zhang Y F, Jia J F, Han T Z, Tang Z, Shen Q T, Guo Y, Qiu Z Q, Xue Q K 2005 Phys. Rev. Lett. 95 096802Google Scholar

    [14]

    Zhang Y F, Tang Z, Han T Z, Ma X C, Jia J F, Xue Q K, Xun K, Wu S C 2007 Appl. Phys. Lett. 90 093120Google Scholar

    [15]

    Qi Y, Ma X C, Jiang P, Ji S H, Fu Y H, Jia J F, Xue Q K, Zhang S B 2007 Appl. Phys. Lett. 90 013109Google Scholar

    [16]

    Bao X Y, Zhang Y F, Wang Y P, Jia J F, Xue Q K, Xie X C, Zhao Z X 2005 Phys. Rev. Lett. 95 247005Google Scholar

    [17]

    Ma L Y, Tang L, Guan Z L, He K, An K, Ma X C, Jia J F, Xue Q K, Han Y, Huang S, Liu F 2006 Phys. Rev. Lett. 97 266102Google Scholar

    [18]

    Ma X C, Jiang P, Qi Y, Jia J F, Yang Y, Duan W H, Li W X, Bao X H, Zhang S B, Xue Q K 2007 Proc. Nat. Acad. Sci. USA 104 9204Google Scholar

    [19]

    Jia J F, Li S C, Zhang Y F, Xue Q K 2007 J. Phys. Soc. Jpn. 76 082001Google Scholar

    [20]

    Chiang T C 2014 Science 306 5703

    [21]

    Zhang G H, Qin H J, Teng J, Guo J D, Guo Q L, Dai X, Fang Z, Wu K H 2009 Appl. Phys. Lett. 95 053114Google Scholar

    [22]

    Song C L, Wang Y L, Jiang Y P, Zhang Y, Chang C Z, Wang L L, He K, Chen X, Jia J F, Wang Y Y, Fang Z, Dai X, Xie X C, Qi X L, Zhang S C, Xue Q K, Ma X C 2010 Appl. Phys. Lett. 97 143118Google Scholar

    [23]

    Li Y Y, Wang G A, Zhu X G, Liu M H, Ye C, Chen X, Wang Y Y, He K, Wang L L, Ma X C, Zhang H J, Dai X, Fang Z, Xie X C, Liu Y, Qi X L, Jia J F, Zhang S C, Xue Q K 2010 Adv. Mater. 22 4002Google Scholar

    [24]

    Zhang Y, He K, Chang C Z, Song C L, Wang L L, Chen X, Jia J F, Fang Z, Dai X, Shan W Y, Shen S Q, Niu Q, Qi X L, Zhang S C, Ma X C, Xue Q K 2010 Nat. Phys. 6 584Google Scholar

    [25]

    Cheng P, Song C L, Zhang T, Zhang Y Y, Wang Y L, Jia J F, Wang J, Wang Y Y, Zhu B F, Chen X, Ma X C, He K, Wang L L, Dai X, Fang Z, Xie X C, Qi X L, Liu C X, Zhang S C, Xue Q K 2010 Phys. Rev. Lett. 105 076801Google Scholar

    [26]

    Feng B J, Zhang J, Zhong Q, Li W B, Li S, Li H, Cheng P, Meng S, Chen L, Wu K H 2016 Nat. Chem. 8 564

    [27]

    Chen C Y, Lv H F, Zhang P, Zhuo Z W, Wang Y, Ma C, Li W B, Wang X G, Feng B J, Cheng P, Wu X J, Wu K H, Chen L 2022 Nat. Chem. 14 25Google Scholar

    [28]

    Zhu X T, Cao Y W, Zhang S Y, Jia X, Guo Q L, Yang F, Zhu L F, Zhang J D, Plummer E W, Guo J D 2015 Rev. Sci. Instrum. 86 083902Google Scholar

    [29]

    Li J D, Li J X, Tang J L, Tao Z Y, Xue S W, Liu J X, Peng H L, Chen X Q, Guo J D, Zhu X T 2023 Phys. Rev. Lett. 131 116602Google Scholar

    [30]

    Zhang X, Xu J Y, Tu Y B, Sun K, Tao M L, Xiong Z H, Wu K H, Wang J Z, Xue Q K, Meng S 2018 Phys. Rev. Lett. 121 256001Google Scholar

    [31]

    Yang K, Chen H, Pope T, et al. 2019 Nat. Commun. 10 3599Google Scholar

    [32]

    Zhang Q, Gao A, Meng F, et al. 2021 Nat. Commun. 12 1853Google Scholar

  • 图 1  Si(111)表面生长出的全同铟(a)和铝(b)纳米点阵列

    Fig. 1.  Identical nanocluster arrays of In (a) and Al (b) grown on Si (111) surface.

    图 2  原子级平整的铅薄膜及其超导转变温度和电子态密度随厚度的振荡变化

    Fig. 2.  Pb films with atomic-scale uniformity, and the superconducting transition temperature Tc and the density of states N(EF) as a function of Pb film thickness, demonstrating a nonmonotonic oscillatory behavior in both Tc and N(EF).

    图 3  硒化铋薄膜以五原子层(QL)为单元生长, 从1 QL到6 QL薄膜的角分辨光电子能谱图(a)—(e)和对应的能量分布曲线(f)—(h), 可以清楚看到狄拉克表面态的能隙打开和Rashba型的劈裂[24]

    Fig. 3.  ARPES spectra of Bi2Se3 films at room temperature: (a)–(e) ARPES spectra of 1-6 quintuple layer (QL), the Bi2Se3 films grow in QL-by-QL mode; (f)–(h) the energy distribution curves of (c)–(e) respectively, which clearly shows the opening of the energy gap and the Rashba-type splitting for Dirac surface states[24].

  • [1]

    Yang W S, Jona F, Marcus P M 1983 Phys. Rev. B 27 1394Google Scholar

    [2]

    Yang W S, Jona F, Marcus P M 1983 Phys. Rev. B 28 2049Google Scholar

    [3]

    Jia J F, Zhao R G, and Yang W S 1993 Phys. Rev. B 48 18109Google Scholar

    [4]

    Jia J F, Zhao R G, Yang W S 1993 Phys. Rev. B 48 18101Google Scholar

    [5]

    贾金峰, 赵汝光, 杨威生 1992 物理学报 41 827Google Scholar

    Jia J F, Zhao R G, Yang W S 1992 Acta Phys. Sin. 41 827Google Scholar

    [6]

    Li J L, Jia J F, Liang X J, Liu X, Wang J Z, Xue Q K, Li Z Q, Tse J S, Zhang Z Y, Zhang S B 2002 Phys. Rev. Lett. 88 066101Google Scholar

    [7]

    Jia J F, Wang J Z, Liu X, Xue Q K, Li Z Q, Kawazoe Y, Zhang S B 2002 Appl. Phys. Lett. 80 3186Google Scholar

    [8]

    Jia J F, Liu X, Wang J Z, Li J L, Wang X S, Xue Q K, Li Z Q, Zhang Z Y, Zhang S B 2002 Phys. Rev. B 66 165412Google Scholar

    [9]

    Xu Z, Bai X D, Wang Z L, Wang E G, 2006 J. Am. Chem. Soc. 128 1052Google Scholar

    [10]

    Xu Z, Bai X D, Wang E G 2006 Appl. Phys. Lett. 88 133107Google Scholar

    [11]

    Wang W L, Bai X D, Liu K H, Xu Z, Golberg D, Bando Y, Wang E G 2006 J. Am. Chem. Soc. 128 6530Google Scholar

    [12]

    Guo Y, Zhang Y F, Bao X Y, Han T Z, Tang Z, Zhang L X, Zhu W G, Wang E G, Niu Q, Qiu Z Q, Jia J F, Zhao Z X, Xue Q K 2004 Science 306 1915Google Scholar

    [13]

    Zhang Y F, Jia J F, Han T Z, Tang Z, Shen Q T, Guo Y, Qiu Z Q, Xue Q K 2005 Phys. Rev. Lett. 95 096802Google Scholar

    [14]

    Zhang Y F, Tang Z, Han T Z, Ma X C, Jia J F, Xue Q K, Xun K, Wu S C 2007 Appl. Phys. Lett. 90 093120Google Scholar

    [15]

    Qi Y, Ma X C, Jiang P, Ji S H, Fu Y H, Jia J F, Xue Q K, Zhang S B 2007 Appl. Phys. Lett. 90 013109Google Scholar

    [16]

    Bao X Y, Zhang Y F, Wang Y P, Jia J F, Xue Q K, Xie X C, Zhao Z X 2005 Phys. Rev. Lett. 95 247005Google Scholar

    [17]

    Ma L Y, Tang L, Guan Z L, He K, An K, Ma X C, Jia J F, Xue Q K, Han Y, Huang S, Liu F 2006 Phys. Rev. Lett. 97 266102Google Scholar

    [18]

    Ma X C, Jiang P, Qi Y, Jia J F, Yang Y, Duan W H, Li W X, Bao X H, Zhang S B, Xue Q K 2007 Proc. Nat. Acad. Sci. USA 104 9204Google Scholar

    [19]

    Jia J F, Li S C, Zhang Y F, Xue Q K 2007 J. Phys. Soc. Jpn. 76 082001Google Scholar

    [20]

    Chiang T C 2014 Science 306 5703

    [21]

    Zhang G H, Qin H J, Teng J, Guo J D, Guo Q L, Dai X, Fang Z, Wu K H 2009 Appl. Phys. Lett. 95 053114Google Scholar

    [22]

    Song C L, Wang Y L, Jiang Y P, Zhang Y, Chang C Z, Wang L L, He K, Chen X, Jia J F, Wang Y Y, Fang Z, Dai X, Xie X C, Qi X L, Zhang S C, Xue Q K, Ma X C 2010 Appl. Phys. Lett. 97 143118Google Scholar

    [23]

    Li Y Y, Wang G A, Zhu X G, Liu M H, Ye C, Chen X, Wang Y Y, He K, Wang L L, Ma X C, Zhang H J, Dai X, Fang Z, Xie X C, Liu Y, Qi X L, Jia J F, Zhang S C, Xue Q K 2010 Adv. Mater. 22 4002Google Scholar

    [24]

    Zhang Y, He K, Chang C Z, Song C L, Wang L L, Chen X, Jia J F, Fang Z, Dai X, Shan W Y, Shen S Q, Niu Q, Qi X L, Zhang S C, Ma X C, Xue Q K 2010 Nat. Phys. 6 584Google Scholar

    [25]

    Cheng P, Song C L, Zhang T, Zhang Y Y, Wang Y L, Jia J F, Wang J, Wang Y Y, Zhu B F, Chen X, Ma X C, He K, Wang L L, Dai X, Fang Z, Xie X C, Qi X L, Liu C X, Zhang S C, Xue Q K 2010 Phys. Rev. Lett. 105 076801Google Scholar

    [26]

    Feng B J, Zhang J, Zhong Q, Li W B, Li S, Li H, Cheng P, Meng S, Chen L, Wu K H 2016 Nat. Chem. 8 564

    [27]

    Chen C Y, Lv H F, Zhang P, Zhuo Z W, Wang Y, Ma C, Li W B, Wang X G, Feng B J, Cheng P, Wu X J, Wu K H, Chen L 2022 Nat. Chem. 14 25Google Scholar

    [28]

    Zhu X T, Cao Y W, Zhang S Y, Jia X, Guo Q L, Yang F, Zhu L F, Zhang J D, Plummer E W, Guo J D 2015 Rev. Sci. Instrum. 86 083902Google Scholar

    [29]

    Li J D, Li J X, Tang J L, Tao Z Y, Xue S W, Liu J X, Peng H L, Chen X Q, Guo J D, Zhu X T 2023 Phys. Rev. Lett. 131 116602Google Scholar

    [30]

    Zhang X, Xu J Y, Tu Y B, Sun K, Tao M L, Xiong Z H, Wu K H, Wang J Z, Xue Q K, Meng S 2018 Phys. Rev. Lett. 121 256001Google Scholar

    [31]

    Yang K, Chen H, Pope T, et al. 2019 Nat. Commun. 10 3599Google Scholar

    [32]

    Zhang Q, Gao A, Meng F, et al. 2021 Nat. Commun. 12 1853Google Scholar

  • [1] 刘文超, 罗朝波, 谢紫彤, 彭向阳. 二维SiSnF2中非磁性缺陷的影响和量子尺寸效应. 物理学报, 2025, 74(6): 066401. doi: 10.7498/aps.74.20241503
    [2] 江龙兴, 李庆超, 张旭, 李京峰, 张静, 陈祖信, 曾敏, 吴昊. 基于拓扑/二维量子材料的自旋电子器件. 物理学报, 2024, 73(1): 017505. doi: 10.7498/aps.73.20231166
    [3] 黄雪峰, 陈矗, 李嘉欣, 张敏琦, 李盛姬. 基于激光悬浮的单颗微米粒子/纳米团簇的散射强度分布测量. 物理学报, 2023, 72(17): 174201. doi: 10.7498/aps.72.20230499
    [4] 吴泽飞, 黄美珍, 王宁. 二维莫尔超晶格中的非线性霍尔效应. 物理学报, 2023, 72(23): 237301. doi: 10.7498/aps.72.20231324
    [5] 陈晓娟, 徐康, 张秀, 刘海云, 熊启华. 二维材料体光伏效应研究进展. 物理学报, 2023, 72(23): 237201. doi: 10.7498/aps.72.20231786
    [6] 吴燕飞, 朱梦媛, 赵瑞杰, 刘心洁, 赵云驰, 魏红祥, 张静言, 郑新奇, 申见昕, 黄河, 王守国. 二维范德瓦尔斯异质结构的制备与物性研究. 物理学报, 2022, 71(4): 048502. doi: 10.7498/aps.71.20212033
    [7] 宋蕊, 王必利, 冯凯, 王黎, 梁丹丹. 二维VOBr2单层的结构畸变及其磁性和铁电性. 物理学报, 2022, 71(3): 037101. doi: 10.7498/aps.71.20211516
    [8] 刘天瑶, 刘灿, 刘开辉. 表界面调控米级二维单晶原子制造. 物理学报, 2022, 71(10): 108103. doi: 10.7498/aps.71.20212399
    [9] 孙颖慧, 穆丛艳, 蒋文贵, 周亮, 王荣明. 金属纳米颗粒与二维材料异质结构的界面调控和物理性质. 物理学报, 2022, 71(6): 066801. doi: 10.7498/aps.71.20211902
    [10] 黄玉昊, 张贵涛, 王如倩, 陈乾, 王金兰. 二维双金属铁磁半导体CrMoI6的电子结构与稳定性. 物理学报, 2021, 70(20): 207301. doi: 10.7498/aps.70.20210949
    [11] 宋蕊. 二维VOBr2单层的结构畸变及其磁性和铁电性研究. 物理学报, 2021, (): . doi: 10.7498/aps.70.20211516
    [12] 陈旭凡, 杨强, 胡小会. 过渡金属原子掺杂对二维CrBr3电磁学性能的调控. 物理学报, 2021, 70(24): 247401. doi: 10.7498/aps.70.20210936
    [13] 孟雨欣, 赵漪凡, 李绍春. 褶皱状蜂窝结构的单层二维材料研究进展. 物理学报, 2021, 70(14): 148101. doi: 10.7498/aps.70.20210638
    [14] 刘子媛, 潘金波, 张余洋, 杜世萱. 原子尺度构建二维材料的第一性原理计算研究. 物理学报, 2021, 70(2): 027301. doi: 10.7498/aps.70.20201636
    [15] 钟虓䶮, 李卓. 原子尺度材料三维结构、磁性及动态演变的透射电子显微学表征. 物理学报, 2021, 70(6): 066801. doi: 10.7498/aps.70.20202072
    [16] 王慧, 徐萌, 郑仁奎. 二维材料/铁电异质结构的研究进展. 物理学报, 2020, 69(1): 017301. doi: 10.7498/aps.69.20191486
    [17] 王聪, 刘杰, 张晗. 基于二维纳米材料的超快脉冲激光器. 物理学报, 2019, 68(18): 188101. doi: 10.7498/aps.68.20190751
    [18] 史若宇, 王林锋, 高磊, 宋爱生, 刘艳敏, 胡元中, 马天宝. 基于滑动势能面的二维材料原子尺度摩擦行为的量化计算. 物理学报, 2017, 66(19): 196802. doi: 10.7498/aps.66.196802
    [19] 冯黛丽, 冯妍卉, 张欣欣. 小尺寸铝纳米团簇的相变行为. 物理学报, 2013, 62(8): 083602. doi: 10.7498/aps.62.083602
    [20] 王瑾, 张茹. 基于量子尺寸效应的InP纳米内包层光纤研制及放大性. 物理学报, 2009, 58(3): 1857-1862. doi: 10.7498/aps.58.1857
计量
  • 文章访问数:  3517
  • PDF下载量:  166
出版历程
  • 收稿日期:  2023-11-27
  • 修回日期:  2023-12-09
  • 上网日期:  2023-12-12
  • 刊出日期:  2023-12-05

/

返回文章
返回