搜索

x

留言板

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

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

一维线性非共轭石墨烯基(CH2)n分子链的电子输运

贺艳斌 白熙

引用本文:
Citation:

一维线性非共轭石墨烯基(CH2)n分子链的电子输运

贺艳斌, 白熙

Electron transport of one-dimensional non-conjugated (CH2)n molecule chain coupling to graphene electrode

He Yan-Bin, Bai Xi
PDF
HTML
导出引用
  • 一维非共轭烷烃链虽不具富电子或少电子特征, 但常存在于单分子器件或多肽、蛋白质等生物分子中, 对电子传输产生重要影响. 为理解这类物质的电子输运特征, 本研究设计了一维线性非共轭(CH2)n分子结模型, 并利用密度泛函理论结合非平衡态格林函数的方法, 对(CH2)n(n = 1—12)线性分子链与两个石墨烯电极耦合而成的分子结进行了第一性原理计算. 结果表明, CH2分子链随着n值的变化, 其电导值表现出明显的奇偶振荡现象, 并且随着链长的增加呈指数级的衰减, 这一结果与实验研究取得了很好的一致性, 为理解和设计性能更加优良的单分子器件提供了重要理论依据.
    Although the one-dimensional non-conjugated alkane chain, which has an important influence on the electron transport process, does not possess the characteristics of electron-rich and electron-deficient, it often exists in single-molecule devices and biological molecules such as peptides and proteins. In order to understand the electron transport characteristics of alkane chain, a one-dimensional linear non-conjugate (CH2)n molecular junction model is designed in this study. Subsequently, we conduct the systematic study of the electronic transport behavior of (CH2)n (n = 1–12) molecular linear chain coupling to two graphene electrodes, based on the density functional theory and nonequilibrium Green’s function formalism. The results reveal that the structure and conductance of CH2 chain are highly sensitive to the odevity of CH2 unit. When the value of n is odd, the groups of CH2 extend in a zigzag way from the left electrode to the right electrode in the plane of graphene, while the value of n is even, what is different is that the groups of CH2 are arranged above and below the electrode plane. For this reason, the odd-even behavior of conductance occurs in the (CH2)n (n = 1–12) molecular chain. Furthermore, n is also an important factor to affect their transport properties (odd or even behavior of conductance). The longer the (CH2)n chain, the deeper the suppression in transmission spectrum and the lower the equilibrium conductance. What is more, the conductance decreases exponentially with the increase of molecular length, with a decay constant β of 0.67 and 0.60 for odd and even, respectively, which is in good agreement with the experimental research. Additionally, by analyzing their eigenchannels of odd and even (CH2)n molecular chain, we find that the coplanar σ electron with graphene electrode makes a major contribution to the electronic transport channel. The current-voltage curve of (CH2)n molecular chain exhibits nonlinearity, implying their semiconductor characteristics. The interesting mechanical and electronic transport properties are expected to conduce to further experimental synthesis, design and operation of the single molecular nanodevices.
      通信作者: 贺艳斌, heyb@czmc.edu.cn
    • 基金项目: 山西省高等学校科技创新项目(批准号: 2020L0371)、长治医学院科技创新团队项目(批准号: CX201904)、山西省“1331工程”重点创新团队建设计划和山西省“服务产业创新学科群建设计划”项目(201809)资助的课题
      Corresponding author: He Yan-Bin, heyb@czmc.edu.cn
    • Funds: Project supported by the Scientific and Technological Innovation Programs of Higher Education Institutions in Shanxi, China, STIP (Grant No. 2020L0371), the Scientific and Technological Innovation Team Programs of Changzhi Medical College, China (Grant No. CX201904), the Fund for Shanxi “1331Project” Key Innovative Research Team, China, “1331KIRT”, and the Shanxi Province Service Industry Innovation Discipline Group Construction Plan, China (201809)
    [1]

    Xin N, Guan J, Zhou C, Chen X, Gu C, Li Y, Ratner M A, Nitzan A, Stoddart J F, Guo X 2019 Nat. Rev. Phy. 1 211Google Scholar

    [2]

    Wen B, Cao M, Lu M, Cao W, Shi H, Liu J, Wang X, Jin H, Fang X, Wang W, Yuan J 2014 Adv. Mater. 26 3484Google Scholar

    [3]

    Cao M-S, Wang X X, Zhang M, Cao W Q, Fang X Y, Yuan J 2020 Adv. Mater. 32 1907156Google Scholar

    [4]

    Balci O, Polat E O, Kakenov N, Kocabas C 2015 Nat. Commun. 6 6628Google Scholar

    [5]

    Meng L, Xin N, Hu C, Wang J, Gui B, Shi J, Wang C, Shen C, Zhang G, Guo H, Meng S, Guo X 2019 Nat. Commun. 10 1450Google Scholar

    [6]

    Li Y, Zhao L, Yao Y, Guo X 2019 ACS Appl Bio Mater 3 68Google Scholar

    [7]

    Zhou C, Li X, Gong Z, Jia C, Lin Y, Gu C, He G, Zhong Y, Yang J, Guo X 2018 Nat. Commun. 9 807Google Scholar

    [8]

    Li Y, Yang C, Guo X 2020 Acc. Chem. Res. 53 159Google Scholar

    [9]

    Jia C, Migliore A, Xin N, Huang S, Wang J, Yang Q, Wang S, Chen H, Wang D, Feng B, Liu Z, Zhang G, Qu D H, Tian H, Ratner M A, Xu H Q, Nitzan A, Guo X 2016 Science 352 1443Google Scholar

    [10]

    Zotti L A, Cuevas J C 2018 ACS Omega 3 3778Google Scholar

    [11]

    Guo C, Yu X, Refaely-Abramson S, Sepunaru L, Bendikov T, Pecht I, Kronik L, Vilan A, Sheves M, Cahen D 2016 Proc. Natl. Acad. Sci. U. S. A. 113 10785Google Scholar

    [12]

    Schosser W M, Zotti L A, Cuevas J C, Pauly F 2019 J. Chem. Phys. 150 174705Google Scholar

    [13]

    Mu Y, Zhou Y, Zhang T, Zeng Z-Y, Cheng Y 2017 J. Chem. Eng. Data 62 3889Google Scholar

    [14]

    Liu F T, Cheng Y, Yang F-B, Chen X R 2014 Physica E 56 96Google Scholar

    [15]

    贺园园, 程娜, 赵健伟 2017 化学学报 75 893Google Scholar

    He Y Y, Cheng N, Zhao J W 2017 Acta Chim. Sin. 75 893Google Scholar

    [16]

    左敏, 廖文虎, 吴丹, 林丽娥 2019 物理学报 68 237302Google Scholar

    Zuo M, Liao W H, Wu D, Lin L E 2019 Acta Phys. Sin. 68 237302Google Scholar

    [17]

    刘南舒, 周思, 赵纪军 2019 物理化学学报 35 1142Google Scholar

    Liu N S, Zhou S, Zhao J J 2019 Acta Phys.Chem. Sin. 35 1142Google Scholar

    [18]

    Fang X Y, Yu X X, Zheng H M, Jin H B, Wang L, Cao M-S 2015 Phys. Lett. A 379 2245Google Scholar

    [19]

    Malen J A, Doak P, Baheti K, Tilley T D, Segalman R A, Majumdar A 2009 Nano Lett. 9 1164Google Scholar

    [20]

    Soler J M, Artacho E, Gale J D, García A, Junquera J, Ordejón P, Sánchez-Portal D 2002 J. Phys.: Condens. Matter 14 2745Google Scholar

    [21]

    Artacho E, Anglada E, Diéguez O, Gale J D, García A, Junquera J, Martin R M, Ordejón P, Pruneda J M, Sánchez-Portal D, Soler J M 2008 J. Phys.: Condens. Matter 20 064208Google Scholar

    [22]

    García A, Papior N, Akhtar A, Artacho E, Blum V, Bosoni E, Brandimarte P, Brandbyge M, Cerdá J I, Corsetti F, Cuadrado R, Dikan V, Ferrer J, Gale J, García-Fernández P, García-Suárez V M, García S, Huhs G, Illera S, Korytár R, Koval P, Lebedeva I, Lin L, López-Tarifa P, Mayo S G, Mohr S, Ordejón P, Postnikov A, Pouillon Y, Pruneda M, Robles R, Sánchez-Portal D, Soler J M, Ullah R, Yu V W-z, Junquera J 2020 J. Chem. Phys. 152 204108Google Scholar

    [23]

    Perdew J P, Burke K, Ernzerhof M 1996 Phys. Rev. Lett. 77 3865Google Scholar

    [24]

    Troullier N, Martins J L 1991 Phys. Rev. B 43 1993Google Scholar

    [25]

    Brandbyge M, Mozos J L, Ordejón P, Taylor J, Stokbro K 2002 Phys. Rev. B 65 165401Google Scholar

    [26]

    Frederiksen T, Paulsson M, Brandbyge M, Jauho A P 2007 Phys. Rev. B 75 205413Google Scholar

    [27]

    Momma K, Izumi F 2011 J. Appl. Crystallogr. 44 1272Google Scholar

    [28]

    Lang N D, Avouris P 1998 Phys. Rev. Lett. 81 3515Google Scholar

    [29]

    Li Y J, Li S L, Gong P, Li Y L, Fang X Y, Jia Y H, Cao M-S 2018 Physica E 104 247Google Scholar

    [30]

    陈熙, 张胜利 2018 物理化学学报 34 1061Google Scholar

    Chen X, Zhang S L 2018 Acta Phys.Chem. Sin. 34 1061Google Scholar

    [31]

    柳福提, 张淑华, 程艳, 陈向荣, 程晓洪 2016 物理学报 65 106201Google Scholar

    Liu F T, Zhang S H, Cheng Y, Chen X R, Cheng X H 2016 Acta Phys. Sin. 65 106201Google Scholar

  • 图 1  (CH2)n分子结(n = 1−12)计算模型示意图(金色球: 碳原子C, 红色球: 氧原子O, 灰色球: 氮原子N, 淡粉色球: 氢原子H)

    Fig. 1.  Schematic illustration of (CH2)n (n = 1−12) molecule junction (golden ball: carbon atom, red ball, oxygen atom, grey ball: nitrogen atom, light pink ball: hydrogen atom).

    图 2  (CH2)n分子结(n = 1−12)稳定结构图(每个图中下方为沿y轴方向的俯视图, 上方和右方插图分别为沿xz轴方向的侧视图, 坐标方向和小球颜色说明见图1)

    Fig. 2.  The stable structure of (CH2)n (n = 1−12) molecule junction (on each diagram, the bottom figure is a top view along y axis, the upper and right side inset is a side view along x and z axis, respectively. Coordinate direction and color description of the ball are shown in Fig. 1).

    图 3  奇数 (a)和偶数(b) (CH2)n分子结在零偏压下的透射谱图

    Fig. 3.  Transmission coefficient as a function of energy for odd (a) and even (b) (CH2)n molecule junction under zero external bias.

    图 4  奇数 (a)和偶数(b) (CH2)n分子结零偏压电导与分子链长度关系图

    Fig. 4.  The plot of conductance versus chain length for odd (a) and even (b) molecular junction of (CH2)n at the bias of 0 V.

    图 5  奇数 (a)和偶数(b) (CH2)n分子结导电通道图(上方插图为沿x轴的侧视图, 蓝色和黄色区域分别代表得失电荷密度)

    Fig. 5.  The eigenchannels of odd (a) and even (b) molecular junction of (CH2)n (The upper inset is a side view along x axis, blue and yellow areas denote the gain and loss of electron density)

    图 6  奇数 (a)和偶数(b) (CH2)n分子结在不同电压下的电流图

    Fig. 6.  The current of odd (a) and even (b) (CH2)n molecular junction under different external bias.

    表 1  (CH2)n分子结(n = 1—12)的平均键长、键角和结合能(d1, d2, d3, α图2所示)

    Table 1.  The average bond length, bond angle and binding energy in (CH2)n (n = 1–12) molecule junction (the bond length of d1, d2, d3 and the bond angle α are shown in Fig. 2).

    nd1N—Nd2C—Cd3C—Cα/(°)ΔE/eV
    12.44–10.90
    23.671.55–11.04
    35.101.552.542.25–11.12
    46.221.532.5518.81–11.45
    57.711.552.630.91–11.52
    68.911.542.5915.12–11.55
    710.301.552.620.41–11.66
    811.511.542.6014.61–11.57
    912.901.552.620.80–11.69
    1014.111.542.6114.93–11.60
    1115.511.552.621.85–11.72
    1216.711.542.6114.19–11.60
    下载: 导出CSV
  • [1]

    Xin N, Guan J, Zhou C, Chen X, Gu C, Li Y, Ratner M A, Nitzan A, Stoddart J F, Guo X 2019 Nat. Rev. Phy. 1 211Google Scholar

    [2]

    Wen B, Cao M, Lu M, Cao W, Shi H, Liu J, Wang X, Jin H, Fang X, Wang W, Yuan J 2014 Adv. Mater. 26 3484Google Scholar

    [3]

    Cao M-S, Wang X X, Zhang M, Cao W Q, Fang X Y, Yuan J 2020 Adv. Mater. 32 1907156Google Scholar

    [4]

    Balci O, Polat E O, Kakenov N, Kocabas C 2015 Nat. Commun. 6 6628Google Scholar

    [5]

    Meng L, Xin N, Hu C, Wang J, Gui B, Shi J, Wang C, Shen C, Zhang G, Guo H, Meng S, Guo X 2019 Nat. Commun. 10 1450Google Scholar

    [6]

    Li Y, Zhao L, Yao Y, Guo X 2019 ACS Appl Bio Mater 3 68Google Scholar

    [7]

    Zhou C, Li X, Gong Z, Jia C, Lin Y, Gu C, He G, Zhong Y, Yang J, Guo X 2018 Nat. Commun. 9 807Google Scholar

    [8]

    Li Y, Yang C, Guo X 2020 Acc. Chem. Res. 53 159Google Scholar

    [9]

    Jia C, Migliore A, Xin N, Huang S, Wang J, Yang Q, Wang S, Chen H, Wang D, Feng B, Liu Z, Zhang G, Qu D H, Tian H, Ratner M A, Xu H Q, Nitzan A, Guo X 2016 Science 352 1443Google Scholar

    [10]

    Zotti L A, Cuevas J C 2018 ACS Omega 3 3778Google Scholar

    [11]

    Guo C, Yu X, Refaely-Abramson S, Sepunaru L, Bendikov T, Pecht I, Kronik L, Vilan A, Sheves M, Cahen D 2016 Proc. Natl. Acad. Sci. U. S. A. 113 10785Google Scholar

    [12]

    Schosser W M, Zotti L A, Cuevas J C, Pauly F 2019 J. Chem. Phys. 150 174705Google Scholar

    [13]

    Mu Y, Zhou Y, Zhang T, Zeng Z-Y, Cheng Y 2017 J. Chem. Eng. Data 62 3889Google Scholar

    [14]

    Liu F T, Cheng Y, Yang F-B, Chen X R 2014 Physica E 56 96Google Scholar

    [15]

    贺园园, 程娜, 赵健伟 2017 化学学报 75 893Google Scholar

    He Y Y, Cheng N, Zhao J W 2017 Acta Chim. Sin. 75 893Google Scholar

    [16]

    左敏, 廖文虎, 吴丹, 林丽娥 2019 物理学报 68 237302Google Scholar

    Zuo M, Liao W H, Wu D, Lin L E 2019 Acta Phys. Sin. 68 237302Google Scholar

    [17]

    刘南舒, 周思, 赵纪军 2019 物理化学学报 35 1142Google Scholar

    Liu N S, Zhou S, Zhao J J 2019 Acta Phys.Chem. Sin. 35 1142Google Scholar

    [18]

    Fang X Y, Yu X X, Zheng H M, Jin H B, Wang L, Cao M-S 2015 Phys. Lett. A 379 2245Google Scholar

    [19]

    Malen J A, Doak P, Baheti K, Tilley T D, Segalman R A, Majumdar A 2009 Nano Lett. 9 1164Google Scholar

    [20]

    Soler J M, Artacho E, Gale J D, García A, Junquera J, Ordejón P, Sánchez-Portal D 2002 J. Phys.: Condens. Matter 14 2745Google Scholar

    [21]

    Artacho E, Anglada E, Diéguez O, Gale J D, García A, Junquera J, Martin R M, Ordejón P, Pruneda J M, Sánchez-Portal D, Soler J M 2008 J. Phys.: Condens. Matter 20 064208Google Scholar

    [22]

    García A, Papior N, Akhtar A, Artacho E, Blum V, Bosoni E, Brandimarte P, Brandbyge M, Cerdá J I, Corsetti F, Cuadrado R, Dikan V, Ferrer J, Gale J, García-Fernández P, García-Suárez V M, García S, Huhs G, Illera S, Korytár R, Koval P, Lebedeva I, Lin L, López-Tarifa P, Mayo S G, Mohr S, Ordejón P, Postnikov A, Pouillon Y, Pruneda M, Robles R, Sánchez-Portal D, Soler J M, Ullah R, Yu V W-z, Junquera J 2020 J. Chem. Phys. 152 204108Google Scholar

    [23]

    Perdew J P, Burke K, Ernzerhof M 1996 Phys. Rev. Lett. 77 3865Google Scholar

    [24]

    Troullier N, Martins J L 1991 Phys. Rev. B 43 1993Google Scholar

    [25]

    Brandbyge M, Mozos J L, Ordejón P, Taylor J, Stokbro K 2002 Phys. Rev. B 65 165401Google Scholar

    [26]

    Frederiksen T, Paulsson M, Brandbyge M, Jauho A P 2007 Phys. Rev. B 75 205413Google Scholar

    [27]

    Momma K, Izumi F 2011 J. Appl. Crystallogr. 44 1272Google Scholar

    [28]

    Lang N D, Avouris P 1998 Phys. Rev. Lett. 81 3515Google Scholar

    [29]

    Li Y J, Li S L, Gong P, Li Y L, Fang X Y, Jia Y H, Cao M-S 2018 Physica E 104 247Google Scholar

    [30]

    陈熙, 张胜利 2018 物理化学学报 34 1061Google Scholar

    Chen X, Zhang S L 2018 Acta Phys.Chem. Sin. 34 1061Google Scholar

    [31]

    柳福提, 张淑华, 程艳, 陈向荣, 程晓洪 2016 物理学报 65 106201Google Scholar

    Liu F T, Zhang S H, Cheng Y, Chen X R, Cheng X H 2016 Acta Phys. Sin. 65 106201Google Scholar

  • [1] 梁锦涛, 颜晓红, 张影, 肖杨. 硼或氮掺杂的锯齿型石墨烯纳米带的非共线磁序与电子输运性质. 物理学报, 2019, 68(2): 027101. doi: 10.7498/aps.68.20181754
    [2] 吴宇, 蔡绍洪, 邓明森, 孙光宇, 刘文江. 聚噻吩单链量子热输运的第一性原理研究. 物理学报, 2018, 67(2): 026501. doi: 10.7498/aps.67.20171198
    [3] 柳福提, 张淑华, 程艳, 陈向荣, 程晓洪. (GaAs)n(n=1-4)原子链电子输运性质的理论计算. 物理学报, 2016, 65(10): 106201. doi: 10.7498/aps.65.106201
    [4] 陈晓彬, 段文晖. 低维纳米材料量子热输运与自旋热电性质 ——非平衡格林函数方法的应用. 物理学报, 2015, 64(18): 186302. doi: 10.7498/aps.64.186302
    [5] 柳福提, 程艳, 陈向荣, 程晓洪. GaAs纳米结点电子输运性质的第一性原理计算. 物理学报, 2014, 63(13): 137303. doi: 10.7498/aps.63.137303
    [6] 柳福提, 程艳, 陈向荣, 程晓洪, 曾志强. Au-Si60-Au分子结电子输运性质的理论计算. 物理学报, 2014, 63(17): 177304. doi: 10.7498/aps.63.177304
    [7] 柳福提, 程艳, 羊富彬, 程晓洪, 陈向荣. Au-Si-Au结点电子输运性质的第一性原理计算. 物理学报, 2013, 62(10): 107401. doi: 10.7498/aps.62.107401
    [8] 柳福提, 程艳, 羊富彬, 程晓洪, 陈向荣. Si4团簇电子输运性质的第一性原理计算 . 物理学报, 2013, 62(14): 140504. doi: 10.7498/aps.62.140504
    [9] 管东波, 毛健. Magnli相亚氧化钛Ti8O15的电子结构和光学性能的第一性原理研究. 物理学报, 2012, 61(1): 017102. doi: 10.7498/aps.61.017102
    [10] 范志强, 谢芳. 硼氮原子取代掺杂对分子器件负微分电阻效应的影响. 物理学报, 2012, 61(7): 077303. doi: 10.7498/aps.61.077303
    [11] 段玲, 胡飞, 丁建文. 准一维纳米线电子输运的梯度无序效应. 物理学报, 2011, 60(11): 117201. doi: 10.7498/aps.60.117201
    [12] 安义鹏, 杨传路, 王美山, 马晓光, 王德华. C20F20分子电子输运性质的第一性原理研究. 物理学报, 2010, 59(3): 2010-2015. doi: 10.7498/aps.59.2010
    [13] 孙源, 明星, 孟醒, 孙正昊, 向鹏, 兰民, 陈岗. 多铁材料BaCoF4电子结构的第一性原理研究. 物理学报, 2009, 58(8): 5653-5660. doi: 10.7498/aps.58.5653
    [14] 宋久旭, 杨银堂, 刘红霞, 张志勇. 掺氮碳化硅纳米管电子结构的第一性原理研究. 物理学报, 2009, 58(7): 4883-4887. doi: 10.7498/aps.58.4883
    [15] 张计划, 丁建文, 卢章辉. Co掺杂MgF2电子结构和光学特性的第一性原理研究. 物理学报, 2009, 58(3): 1901-1907. doi: 10.7498/aps.58.1901
    [16] 郑小宏, 戴振翔, 王贤龙, 曾雉. B与N掺杂对单层石墨纳米带自旋极化输运的影响. 物理学报, 2009, 58(13): 259-S265. doi: 10.7498/aps.58.259
    [17] 郑新亮, 郑继明, 任兆玉, 郭平, 田进寿, 白晋涛. 钽硅团簇电子输运性质的第一性原理研究. 物理学报, 2009, 58(8): 5709-5715. doi: 10.7498/aps.58.5709
    [18] 倪建刚, 刘 诺, 杨果来, 张 曦. 第一性原理研究BaTiO3(001)表面的电子结构. 物理学报, 2008, 57(7): 4434-4440. doi: 10.7498/aps.57.4434
    [19] 张金奎, 邓胜华, 金 慧, 刘悦林. ZnO电子结构和p型传导特性的第一性原理研究. 物理学报, 2007, 56(9): 5371-5375. doi: 10.7498/aps.56.5371
    [20] 潘志军, 张澜庭, 吴建生. CoSi电子结构第一性原理研究. 物理学报, 2005, 54(1): 328-332. doi: 10.7498/aps.54.328
计量
  • 文章访问数:  2618
  • PDF下载量:  53
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-06-21
  • 修回日期:  2020-09-28
  • 上网日期:  2021-02-01
  • 刊出日期:  2021-02-20

一维线性非共轭石墨烯基(CH2)n分子链的电子输运

  • 长治医学院药学系, 长治 046000
  • 通信作者: 贺艳斌, heyb@czmc.edu.cn
    基金项目: 山西省高等学校科技创新项目(批准号: 2020L0371)、长治医学院科技创新团队项目(批准号: CX201904)、山西省“1331工程”重点创新团队建设计划和山西省“服务产业创新学科群建设计划”项目(201809)资助的课题

摘要: 一维非共轭烷烃链虽不具富电子或少电子特征, 但常存在于单分子器件或多肽、蛋白质等生物分子中, 对电子传输产生重要影响. 为理解这类物质的电子输运特征, 本研究设计了一维线性非共轭(CH2)n分子结模型, 并利用密度泛函理论结合非平衡态格林函数的方法, 对(CH2)n(n = 1—12)线性分子链与两个石墨烯电极耦合而成的分子结进行了第一性原理计算. 结果表明, CH2分子链随着n值的变化, 其电导值表现出明显的奇偶振荡现象, 并且随着链长的增加呈指数级的衰减, 这一结果与实验研究取得了很好的一致性, 为理解和设计性能更加优良的单分子器件提供了重要理论依据.

English Abstract

    • 近二十年来, 各种纳米级电极制备技术的出现, 促进了单分子电子学领域在实验和理论方面的快速发展[1]. 石墨烯作为一种优良的电子导体, 在电磁屏蔽[2]、微型智能传感[3]及雷达吸波材料[4]等方面得到了广泛应用. 特别是Guo等[5-8]利用单层石墨烯作为电极, 构建的“石墨烯-分子-石墨烯(graphene-molecule-graphene, GMG)”异质结, 表现出了结构精确可控、信号可调的特点, 在揭示光开关[5]、化学/生物传感[6]、氢键作用[7]等物质转换规律和生命本征现象方面显示了独特优势, 为将来颠覆性芯片技术和精准分子诊断技术的发展提供了基础[8]. 尤其是Jia等[9]将二芳烯作为功能中心, 在石墨烯电极间构建了一个光、电双模式诱导的单分子全可逆光电子器件. 在这一突破性工作中, 为了避免分子与电极接触界面间的耦合, 在功能分子二芳烯的两边各引入了关键性的3个亚甲基(CH2). 此外, 研究固态分子结中生物分子的电子转移对于理解许多重要的生理过程具有重要意义[10]. 多肽单分子电子输运的实验[11]和理论[12]研究均说明, 多肽分子中处于电子传输方向的主链或分子侧链上的CH2和CH3等烷烃基团对于电子的输运都产生了一定的影响. 然而, 如何理解这些烷烃在分子结中的电子输运, 烷烃基团的数量对于电导有什么影响仍未见系统的研究.

      近来, 有关一维分子链电子输运的研究引起了人们的广泛关注. Mu等[13]通过理论方法研究了一维(SiC)n (n = 1—4)分子链随着n的改变, 其导电性质的变化规律, 阐明了这一分子链的半导体导电特征. Liu等针对(GaAs)n (n = 1—4)原子链的研究, 发现了n值的不同导致的金属或负微分电阻现象; 而针对Si原子链的电子输运研究却发现了明显的奇偶电导振荡现象[14]. 此外, 一维富电子有机共轭分子链[15,16]、点缺陷石墨烯[17]、石墨烯层数[18]等不同分子结的导电性质也被研究报道, 证实石墨烯的输运机制主要来源于单层石墨烯的弹道散射或少层石墨烯的声子散射作用.

      烷烃分子由于不具有金属或分子导体的导电性质, 同时也不属于富电子的共轭体系, 因此在分子器件的研究中较少被探讨. Malen等[19]通过实验研究了(CH2)n分子结的导电特征, 发现其电导值的衰减与链长之间呈指数的关系, 但并没有详细阐释其电子输运机制. 为进一步理解非共轭烷烃在分子器件中的导电本质, 讨论不同链长烷烃分子结的电子输运特性, 本文将基于密度泛函理论(DFT)和非平衡态格林函数(NEGF)相结合的第一性原理方法, 对石墨烯基线性(CH2)n分子结(n = 1—12)的电子输运性质进行了理论研究, 以便为单分子器件的设计提供更多的理论指导.

    • 为研究石墨烯基CH2分子结的电子输运特性, 我们构建了如图1所示的“石墨烯-分子-石墨烯(GMG)”异质结计算模型. 由于该“GMG”异质结在实验过程中均由含氨基的分子与边缘羧基化的石墨烯电极耦合而得, 因此本文计算模型通过将两头含有氨基的CH2分子链通过酰胺键偶联到两侧石墨烯电极间, 设计了1个“三明治”模型的双电极电子输运体系. CH2分子链沿着z轴方向延伸, 且为电子输运的方向, 而xy方向为周期性结构. 左右石墨烯电极间距随着CH2基团的数目增加而递增, 但均处于同一平面, 每个电极区各包含10 × 6阵列的60个碳原子. 散射区所有边缘碳原子均使用氢原子饱和. 模拟过程中, 整个分子模型置于长方形盒子中, 石墨烯平面距离盒子y方向上下边界均为6 Å.

      图  1  (CH2)n分子结(n = 1−12)计算模型示意图(金色球: 碳原子C, 红色球: 氧原子O, 灰色球: 氮原子N, 淡粉色球: 氢原子H)

      Figure 1.  Schematic illustration of (CH2)n (n = 1−12) molecule junction (golden ball: carbon atom, red ball, oxygen atom, grey ball: nitrogen atom, light pink ball: hydrogen atom).

      本研究基于密度泛函理论(DFT)与非平衡态格林函数(NEGF)相结合的第一性原理计算方法, 采用SIESTA 4.1[20,21]软件来完成相关计算. 其中, 密度泛函采用基于数值基组的方法来求解Kohn-Shan方程. SIESTA软件包近来被广泛用于各种材料模拟中, 并取得了与实验结果良好的一致性[22]. 分子结的结构优化过程中, 对电极区原子进行固定, 其余原子自由弛豫到受力小于0.01 eV / Å时达到收敛标准. 计算采用PBE[23]的广义梯度近似(GGA)方法来描述电子的交换关联作用, 核芯势采用模守恒的Troullier-Martins赝势[24]. 实空间积分格点截断设置为250 Ry, 体系的电子温度设置为300 K. 电极区使用了单ζ极化基组(SZP), 而其余原子采用了双ζ基组(DZP)对波函数进行展开.

      在电子输运计算中(TRANSIESTA程序包[25]), 使用Monkhorst-Pack采样方法, 选取1 × 1 × 50的k点来计算电子的透射率. 通过(CH2)n分子结和2个石墨烯电极的电流可以根据以下Landauer-Buttiker公式来计算.

      $ I\left( V \right) = \frac{{2e}}{h}\int T \left( {E,V} \right)\left[ {f(E - {\mu _{\rm{L}}}) - f(E - {\mu _{\rm{R}}})} \right]{\rm{d}}E, $

      其中: $ {f}{(}{E}-\mu _{\rm{L}}{)} $$ {f}{(}{E}-\mu _{\rm{R}}{)} $为左右电极的费米-狄拉克分布函数; $ \mu _{\rm{L}} $$ \mu _{\rm{R}} $为左右电极的化学势; $ {T}\left({E}{, }{}{V}\right) $是电压为V时能量为E的电子从左到右传输的透射系数, 能量积分区间为[$ \mu _{\rm{L}} $, $ \mu _{\rm{R}} $]; 当外电压V = 0 V时, 体系的平衡电导可由零偏压下平衡态的费米能级EF处的透射值T乘以量子电导G0 (G0 = ${{2{e^2}}}/{h}$)得出. 此外, 采用Inelastica软件包[26]计算了不同分子结的导电通道, 采用 VESTA程序[27]对所有模型结构进行了可视化.

    • 为了更好地模拟亚甲基分子链的电子传输特性, 在“石墨烯-(CH2)n-石墨烯”分子结模型中, 我们参考游离烷烃分子中C—C键和C—N键的键长, 首先对两头连接有氨基的(CH2)n分子链进行了结构优化, 然后根据优化的N-N之间的距离设置了石墨烯两侧耦合点的间距, 最后将含有氨基的亚甲基分子链通过酰胺键连接于两侧羧基化的石墨烯电极间, 并再次进行结构优化, 得到如图2所示的n从1到12的一维线性非共轭亚甲基分子结稳定构型. 进一步地, 我们依据公式ΔE = E [GMG] – E [石墨烯电极]– E [分子链], 计算了这些分子结的结合能并列于表1中, 结果发现, 其结合能在数值上非常接近, 均介于–10.90 eV到–11.72 eV之间, 说明优化得到的这些一维非共轭(CH2)n分子结为稳定的平衡结构.

      图  2  (CH2)n分子结(n = 1−12)稳定结构图(每个图中下方为沿y轴方向的俯视图, 上方和右方插图分别为沿xz轴方向的侧视图, 坐标方向和小球颜色说明见图1)

      Figure 2.  The stable structure of (CH2)n (n = 1−12) molecule junction (on each diagram, the bottom figure is a top view along y axis, the upper and right side inset is a side view along x and z axis, respectively. Coordinate direction and color description of the ball are shown in Fig. 1).

      nd1N—Nd2C—Cd3C—Cα/(°)ΔE/eV
      12.44–10.90
      23.671.55–11.04
      35.101.552.542.25–11.12
      46.221.532.5518.81–11.45
      57.711.552.630.91–11.52
      68.911.542.5915.12–11.55
      710.301.552.620.41–11.66
      811.511.542.6014.61–11.57
      912.901.552.620.80–11.69
      1014.111.542.6114.93–11.60
      1115.511.552.621.85–11.72
      1216.711.542.6114.19–11.60

      表 1  (CH2)n分子结(n = 1—12)的平均键长、键角和结合能(d1, d2, d3, α图2所示)

      Table 1.  The average bond length, bond angle and binding energy in (CH2)n (n = 1–12) molecule junction (the bond length of d1, d2, d3 and the bond angle α are shown in Fig. 2).

      分析图2中优化得到的不同n值时分子结的稳定构型, 可以发现, 当n为奇数时, 亚甲基分子链中碳原子基本都处于石墨烯的电极平面内(xy平面), 并以“之”字形沿着分子链方向(z轴)左右排列, 而n为偶数时, 碳原子虽仍为“之”字形排列, 但均在石墨烯电极平面上下波动. 同时, 从图2右侧的插图中, 还可以看出, 电极间存在奇数个CH2基团时, 分布于“之”字形左右的数个CH2几乎重叠, 而为偶数个CH2基团时, “之”字形分布的CH2在电极平面上下沿z轴方向以一定角度扭转延伸, 并随着n值的增加, 这一扭转角α逐渐变小(见表1α值). 此外, 表1中列出了不同n值时, 分子结的键长、键角数值. d1代表了分子结中(CH2)n分子链上两端N原子之间的距离. 从表中这一数值变化中可以发现, n值由奇数增加为偶数时, d1的值增加了约1.40 Å, 而由偶数增加为奇数时, d1值却只增加了约1.20 Å. 这些变化充分说明一维石墨烯基(CH2)n分子结在结构上有明显的奇偶特征. 其余代表烷烃分子中C—C键长的d2, 以及代表烷烃分子中间位C—C键的键长d3并没有随着n的增加而产生明显的差异.

    • 根据优化得到的(CH2)n分子结(n = 1—12)平衡结构, 本文计算了零偏压时, 各(CH2)n分子链的透射谱. 结果如图3所示, 所有分子结的透射谱表现出了类似的峰形特征, 并且由于亚甲基的非共轭结构, 使得透射谱在费米能级附近没有明显的隧穿共振峰, 从而导致了这一分子链较差的导电性质. 同时, 随着分子链n值的增加, 这一分子结的透射系数明显地降低. 尤其值得注意的是, 比较图3(a)奇数分子链和图3(b)偶数分子链的透射谱, 能够发现, 偶数个CH2基团的分子链导电性明显差于奇数分子链, 并且由于结构上的奇偶性, 这一分子结的电导值也表现出明显的振荡衰减(见图3(a)中插图), 这一现象在Si[14], C[28]等原子链的电子输运研究中也被发现, 但Si原子链的电导振荡现象会随着原子数的增加而逐渐减弱, 而在CH2分子链中却并没有发现这一趋势.

      图  3  奇数 (a)和偶数(b) (CH2)n分子结在零偏压下的透射谱图

      Figure 3.  Transmission coefficient as a function of energy for odd (a) and even (b) (CH2)n molecule junction under zero external bias.

    • 电荷输运的机理有两种: 不依赖温度的连续隧穿和依赖于温度的非连续性跳跃. 连续隧穿机理认为, 分子链的长度决定了电子从左电极到右电极传输的距离, 成为影响单分子电子器件导电性质的一个关键因素. 通常情况下, 分子器件的电导将随着长度呈指数性衰减, 其关系可用如(2)式表示[1]:

      $ G = {G_0}{\rm{exp}}( - \beta L), $

      其中: G为分子电导, G0为量子电导常数, L是分子的长度, β为衰减常数. 因此, 本文根据零偏压下(CH2)n分子结在费米能级处的透射系数得到了n从1到12的不同分子链导电值, 并将分子链中N-N键的键长d1作为分子的长度, 总结了(CH2)n分子链的电导和长度之间的关系, 结果如图4所示. 由于该分子链随着长度的增加, 其电导表现出了明显的振荡现象. 因此, 从图4可以看出, 奇数或偶数分子链的电导均随分子长度变化呈指数性衰减的趋势, 并且其线性方程分别为: ln(G) = –0.67L – 7.78 (图4中奇数a, 线性相关系数r = 0.9966)和ln(G) = –0.60L –11.00 (图4中偶数b, 线性相关系数r = 0.9978). 根据线性方程, 可知奇数和偶数(CH2)n分子链的衰减常数基本一致, 分别为0.67和0.60, 相比较于实验报道的衰减常数0.94而言[19], 尽管β值要小一些, 但实验中使用的电极为金电极. 相关研究也证实, β值也将会因耦合基团、溶剂环境等测试条件的改变而改变[1].

      图  4  奇数 (a)和偶数(b) (CH2)n分子结零偏压电导与分子链长度关系图

      Figure 4.  The plot of conductance versus chain length for odd (a) and even (b) molecular junction of (CH2)n at the bias of 0 V.

      为了更好地理解(CH2)n奇偶分子链导电性质的差异和振荡现象, 本文分别以n为11和12的分子结为例, 计算了这个一维非共轭体系的导电本征通道(图5). 从图5可以看出, 不论是奇数还是偶数分子链, 其导电作用主要都是通过C—C键之间的定域σ电子进行电荷传输, 相对于π键的共轭体系, 导电性质要更差一些, 因而, 计算得到的衰减常数β值也要更大(一般来说, β值越小, 共轭体系越明显, 导电性也更好). 此外, 通过比较图5(a)图5(b), 可以发现两者最大的区别是奇数分子链中, 每个CH2基团及对应的σ键与石墨烯电极基本处于同一平面, 因而也更有利于电子传输, 而偶数分子链中, 两者正好为垂直关系. 这充分说明了一维非共轭(CH2)n分子结出现电导奇偶振荡的本质原因. 进一步分析其电子输运的本质, n为奇数时, 应该是以声子散射机制为主; n为偶数时, 应该以定程跳跃散射为主. 因此, 表现出的奇数电导比偶数电导大1个数量级以上[29].

      图  5  奇数 (a)和偶数(b) (CH2)n分子结导电通道图(上方插图为沿x轴的侧视图, 蓝色和黄色区域分别代表得失电荷密度)

      Figure 5.  The eigenchannels of odd (a) and even (b) molecular junction of (CH2)n (The upper inset is a side view along x axis, blue and yellow areas denote the gain and loss of electron density)

    • 外偏压也是影响单分子器件导电性质的一个重要因素. 本文进一步计算了不同链长(CH2)n分子结在外偏压为0—2.0 V范围内不同电压下的电流值, 结果如图6所示. 可以看出, (CH2)n分子结的I-V曲线与文献报道的石墨二炔纳米带的I-V曲线特征相似, 具有半导体特征[30], 尽管奇数分子链相较偶数分子链的电流值要大10倍以上, 但由于烷烃链的导电性较差, 导致该分子链整体电流值均较小, 仅为几到几十nA的数量级. 此外, 从图中还可以看出, 奇数分子链的导通电压约为1.0 V(图6(a)), 而偶数分子链约为1.5 V (图6(b)). 这一偶数分子链较差的导电特性也与之前对导电通道的分析相吻合. 值得注意的是, 当分子链长度为n = 4, 外电压大于1.6 V时, 其电流值超过n = 2的分子链电流值, 根据Landauer-Buttiker公式, 分析其原因为偏压窗口增大后, n = 4的分子链积分计算所得到的透射率增大, 导致了电流的增加[31].

      图  6  奇数 (a)和偶数(b) (CH2)n分子结在不同电压下的电流图

      Figure 6.  The current of odd (a) and even (b) (CH2)n molecular junction under different external bias.

    • 石墨烯基单分子器件中, 线性非共轭 (CH2)n基团的引入可以避免功能分子与电极的耦合, 对电子输运性质有着重要影响. 本文针对这一问题, 通过密度泛函理论与非平衡态格林函数相结合的方法, 对(CH2)n (n = 1—12)线性分子链与石墨烯半无限电极耦合构成的“石墨烯-(CH2)n-石墨烯”分子结的电子输运行为进行了模拟研究. 通过对12个不同链长的分子结稳定结构的优化, 发现奇数个(CH2)n基团的分子链沿着石墨烯平面从左电极向着右电极延伸, 而偶数分子链中CH2基团却在石墨烯平面上下分布且以一定的旋转角向右延伸. 这种结构上的差异也导致了透射谱的相应变化和电导的奇偶振荡. 进一步地, 我们还通过研究不同分子链在零偏压时的电导与链长之间的关系, 分析了此一维线性非共轭分子链的电子传输机理. 结果发现, (CH2)n分子链的电子传输符合连续隧穿的电子输运机制, 并且计算得到其电导随分子链指数性衰减的常数β值约为0.67, 与文献报道的实验结果基本一致. 通过I-V曲线, 也可以发现这个分子链半导体的导电特征和电流的奇偶变化依赖性. 此外, 我们还计算了奇、偶不同分子链的导电本征通道, 从而进一步佐证了这一分子链的导电特性, 并解释了该研究模型中奇偶振荡现象的本质. 希望本文的研究结果能为单分子纳米器件的设计与开发提供重要理论参考.

参考文献 (31)

目录

    /

    返回文章
    返回