搜索

x

留言板

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

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

二噻吩硼烷异构体分子结构测定的第一性原理研究

彭淑平 黄旭东 刘乾 任鹏 伍丹 范志强

引用本文:
Citation:

二噻吩硼烷异构体分子结构测定的第一性原理研究

彭淑平, 黄旭东, 刘乾, 任鹏, 伍丹, 范志强

First-principles study of single-molecule-structure determination of dithienoborepin isomers

Peng Shu-Ping, Huang Xu-Dong, Liu Qian, Ren Peng, Wu Dan, Fan Zhi-Qiang
PDF
HTML
导出引用
  • 本文在第一性原理计算基础上结合非平衡格林函数方法, 研究了量子干涉效应对连接镍电极的二噻吩硼烷(dithienoborepin, DTB)分子结自旋输运性质的影响, 并通过氨基和硝基钝化实现了对二噻吩硼烷分子异构体(DTB-A和DTB-B)的区分. 结果表明, 原始的DTB-A和DTB-B分子结在费米能级两侧都有一个自旋向上透射峰和一个自旋向下透射峰, 且两个透射峰的能量位置和高度基本相同. 因此, 原始DTB-A和DTB-B分子结的自旋向上和自旋向下电流曲线基本重合, 不能被明显区分. 然而, 研究发现量子干涉效应能不同程度地增强氨基钝化DTB-A分子结费米能级两侧分子轨道的自旋极化输运能力, 并减弱氨基钝化DTB-B分子结费米能级两侧分子轨道的自旋极化输运能力. 此外, 研究还发现量子干涉效应可以显著提高硝基钝化DTB-B分子结费米能级两侧分子轨道的自旋极化输运能力, 同时减弱硝基钝化DTB-A分子结费米能级两侧分子轨道的自旋极化输运能力. 由于量子干涉效应对氨基和硝基钝化的DTB异构体分子结自旋输运能力有不同的调制作用, 因此可以通过测量氨基和硝基钝化分子结的自旋电流值来区分DTB分子的两种异构体.
    Previous research results show that the conductance difference in molecular junction caused by quantum interference (QI) effect is an important way to identify isomers or improve the recognition sensitivity. Recently, single-molecule conductance of two fully π-conjugated dithienoborepin (DTB) isomers (DTB-A and DTB-B) with tricoordinate boron centers has been measured by using the scanning tunneling microscopy break junction technique. The result shows that QI can enhance chemical responsivity in single-molecule DTB junction. In this work, the first-principles method based on density functional theory and non-equilibrium Green's function is used to study the influence of QI effect on spin-transport property of DTB molecular junction connected to the nickel electrode, and the purpose of distinguishing DTB isomers (DTB-A and DTB-B) is realized by using amino and nitro passivation. The results show that the pristine DTB-A molecule and DTB-B molecule both have a up-spin transmission peak dominated by HOMO and a down-spin transmission peak dominated by LUMO on both sides of the Fermi level, and the energy positions and coefficients of two transmission peaks are basically the same. Therefore, the up-spin and down-spin current curves of the two junctions basically coincide, so that it is impossible to clearly distinguish the two isomers of DTB molecule simply by spin current. The QI can enhance the spin-polarized transport capability of two orbitals of amino-passivated DTB-A molecule to varying degrees but weaken the spin-polarized transport capability of two orbitals of amino-passivated DTB-B molecule. Therefore, the current of DTB-A molecular junction passivated by amino group is always higher than that of DTB-B molecular junction passivated by amino group. However, the QI can greatly enhance the spin-polarized transport capability of two orbitals of nitro-passivated DTB-B molecule but weaken the spin-polarized transport capability of two orbitals of nitro-passivated DTB-A molecule. Therefore, the current of DTB-B molecular junction passivated by nitro is always higher than that of DTB-A molecular junction passivated by nitro. Because the QI has different effects on the spin-transport capability of DTB-A and DTB-B passivated by amino or nitro group, so the two isomers of DTB molecule can be distinguished by measuring the spin current value. The above conclusions provide more theoretical guidance for the practical preparation of spin molecular junctions and the regulation of their spin-transport performance in the future.
      通信作者: 范志强, zqfan@csust.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 12074046, 12204066)、湖南省教育厅科学基金(批准号: 20C0039, 19C0027)、长沙理工大学青年研究人员助推计划 (批准号: 2019QJCZ022)和湖南省研究生科研创新项目(批准号: CXCLY2022141)资助的课题.
      Corresponding author: Fan Zhi-Qiang, zqfan@csust.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant Nos. 12074046, 12204066), the Scientific Research Fund of Hunan Provincial Education Department, China (Grant Nos. 20C0039, 19C0027), the Young Researchers’ Cultivation Programme of Changsha University of Science & Technology, China (Grant No. 2019QJCZ022), and the Hunan Postgraduate Scientific Research and Innovation Project, China (Grant No. CXCLY2022141).
    [1]

    Aviram A, Ratner M A 1974 Chem. Phys. Lett. 29 277Google Scholar

    [2]

    Perrin M L, Frisenda R, Koole M, Seldenthuis J S, Gil J A C, Valkenier H, Hummelen J C, Renaud N, Grozema F C, Thijssen J M, Dulić D, van der Zant H S J 2014 Nat. Nanotechnol. 9 830Google Scholar

    [3]

    Koley S, Chakrabarti S 2018 Chem. Eur. J. 24 5876Google Scholar

    [4]

    Sharma P, Bernard L S, Bazigos A, Magrez A, Ionescu A M 2015 ACS Nano 9 620Google Scholar

    [5]

    Kumar S, Wang Z, Davila N, Kumari N, Norris K J, Huang X, Strachan J P, Vine D, Kilcoyne A L D, Nishi Y, Williams R S 2017 Nat. Commun. 8 658Google Scholar

    [6]

    Li Z L, Sun F, Bi J J, Liu R, Suo Y Q, Fu H Y, Zhang G P, Song Y Z, Wang D, Wang C K 2019 Physica E 106 270Google Scholar

    [7]

    Liu Q, Li J J, Wu D, Deng X Q, Zhang Z H, Fan Z Q, Chen K Q 2021 Phys. Rev. B 104 045412Google Scholar

    [8]

    Komeda J, Takada K, Maeda H, Fukui N, Tsuji T, Nishihara H 2022 Chem. Eur. J. 28 e202201316Google Scholar

    [9]

    Zhao J, Zeng H, Wang D, Yao G 2020 Appl. Surf. Sci. 519 146203Google Scholar

    [10]

    Petersen M Å, Rasmussen B, Andersen N N, Sauer S P A, Nielsen M B, Beeren S R, Pittelkow M 2017 Chem. Eur. J. 23 17010Google Scholar

    [11]

    Fan Z Q, Zhang Z H, Yang S Y 2020 Nanoscale 12 21750Google Scholar

    [12]

    Jiang J, Kula M, Lu W, Luo Y 2005 Nano Lett. 5 1551Google Scholar

    [13]

    Lambert C J 2015 Chem. Soc. Rev. 44 875Google Scholar

    [14]

    Gehring P, Thijssen J M, van der Zant H S J 2019 Nat. Rev. Phys. 1 381Google Scholar

    [15]

    Manrique D Z, Huang C, Baghernejad M, Zhao X, Al-Owaedi O A, Sadeghi H, Kaliginedi V, Hong W, Gulcur M, Wandlowski T, Bryce M R, Lambert C J 2015 Nat. Commun. 6 6389Google Scholar

    [16]

    Fan Z Q, Sun W Y, Zhang Z H, Deng X Q, Tang G P, Xie H Q 2017 Carbon 122 687Google Scholar

    [17]

    Li Z L, Bi J J, Liu R, Yi X H, Fu H Y, Sun F, Wei M Z, Wang C K 2017 Chin. Phys. B 26 098508Google Scholar

    [18]

    Yang Y, Gantenbein M, Alqorashi A, Wei J, Sangtarash S, Hu D, Sadeghi H, Zhang R, Pi J, Chen L, Huang X, Li R, Liu J, Shi J, Hong W, Lambert C J, Bryce M R 2018 J. Phys. Chem. C 122 14965Google Scholar

    [19]

    Liu X, Sangtarash S, Reber D, Zhang D, Sadeghi H, Shi J, Xiao Z Y, Hong W, Lambert C J, Liu S X 2017 Angew. Chem. 129 179Google Scholar

    [20]

    Borges A, Fung E D, Ng F, Venkataraman L, Solomon G C 2016 J. Phys. Chem. Lett. 7 4825Google Scholar

    [21]

    Frisenda R, Janssen V A E C, Grozema F C, van der Zant H S J, Renaud N 2016 Nat. Chem. 8 1099Google Scholar

    [22]

    Liu R, Bi J J, Xie Z, Yin K, Wang D, Zhang G P, Xiang D, Wang C K, Li Z L 2018 Phys. Rev. Appl. 9 054023Google Scholar

    [23]

    Zhang Y P, Chen L C, Zhang Z Q, Cao J J, Tang C, Liu J, Duan L L, Huo Y, Shao X, Hong W, Zhang H L 2018 J. Am. Chem. Soc. 140 6531Google Scholar

    [24]

    Yang G, Sangtarash S, Liu Z, Li X, Sadeghi H, Tan Z, Li R, Zheng J, Dong X, Liu J, Yang Y, Shi J, Xiao Z, Zhang G, Lambert C, Hong W, Zhang D 2017 Chem. Sci. 8 7505Google Scholar

    [25]

    Li Y, Buerkle M, Li G, Rostamian A, Wang H, Wang Z, Bowler D R, Miyazaki T, Xiang L, Asai Y, Zhou G, Tao N 2019 Nat. Mater. 18 357Google Scholar

    [26]

    Huang B, Liu X, Yuan Y, Hong Z W, Zheng J F, Pei L Q, Shao Y, Li J F, Zhou X S, Chen J Z, Jin S, Mao B W 2018 J. Am. Chem. Soc. 140 17685Google Scholar

    [27]

    Bai J, Daaoub A, Sangtarash S, Li X, Tang Y, Zou Q, Sadeghi H, Liu S, Huang X, Tan Z, Liu J, Yang Y, Shi J, Mészáros G, Chen W, Lambert C, Hong W 2019 Nat. Mater. 18 364Google Scholar

    [28]

    Zheng J, Liu J, Zhuo Y, Li R, Jin X, Yang Y, Chen Z B, Shi J, Xiao Z, Hong W, Tian Z Q 2018 Chem. Sci. 9 5033Google Scholar

    [29]

    Liu J, Zhao X, Al-Galiby Q, Huang X, Zheng J, Li R, Huang C, Yang Y, Shi J, Manrique D Z, Lambert C J, Bryce M R, Hong W 2017 Angew. Chem. 129 13241Google Scholar

    [30]

    Cai S, Deng W, Huang F, Chen L, Tang C, He W, Long S, Li R, Tan Z, Liu J, Shi J, Liu Z, Xiao Z, Zhang D, Hong W 2019 Angew. Chem. Int. Ed. 58 3829Google Scholar

    [31]

    Ozawa H, Baghernejad M, Al-Owaedi O A, Kaliginedi V, Nagashima T, Ferrer J, Wandlowski T, García-Suárez V M, Broekmann P, Lambert C J, Haga M A 2016 Chem. Eur. J. 22 12732Google Scholar

    [32]

    Greenwald J E, Cameron J, Findlay N J, Fu T, Gunasekaran S, Skabara P J, Venkataraman L 2021 Nat. Nanotechnol. 16 313Google Scholar

    [33]

    Trasobares J, Vuillaume D, Théron D, Clément N 2016 Nat. Commun. 7 12850Google Scholar

    [34]

    Niu L L, Fu H Y, Suo Y Q, Liu R, Sun F, Wang S S, Zhang G P, Wang C K, Li Z L 2021 Phys. E 128 114542Google Scholar

    [35]

    Wu D, Huang L, Jia P Z, Cao X H, Fan Z Q, Zhou W X, Chen K Q 2021 Appl. Phys. Lett. 119 063503Google Scholar

    [36]

    Fan Z Q, Xie F, Jiang X W, Wei Z, Li S S 2016 Carbon 110 200Google Scholar

    [37]

    Kushmerick J 2009 Nature 462 994Google Scholar

    [38]

    Richter S, Mentovich E, Elnathan R 2018 Adv. Mater. 30 1706941Google Scholar

    [39]

    Liu R, Han Y, Sun F, Khatri G, Kwon J, Nickle C, Wang L, Wang C K, Thompson D, Li Z L, Nijhuis C A, del Barco E 2022 Adv. Mater. 34 2202135Google Scholar

    [40]

    Fan Z Q, Sun W Y, Jiang X W, Zhang Z H, Deng X Q, Tang G P, Xie H Q, Long M Q 2017 Carbon 113 18Google Scholar

    [41]

    Zhou Q, Yang Y, Ni J, Li Z, Zhang Z 2010 Nano Res. 3 423Google Scholar

    [42]

    Baghernejad M, Van Dyck C, Bergfield J, Levine D R, Gubicza A, Tovar J D, Calame M, Broekmann P, Hong W 2019 Chem. Eur. J. 25 15141Google Scholar

    [43]

    Büttiker M, Imry Y, Landauer R, Pinhas S 1985 Phys. Rev. B 31 6207Google Scholar

    [44]

    Smidstrup S, Markussen T, Vancraeyveld P, Wellendorff J, Schneider J, Gunst T, Verstichel B, Stradi D, Khomyakov P A, Vej-Hansen U G, Lee M E, Chill S T, Rasmussen F, Penazzi G, Corsetti F, Ojanperä A, Jensen K, Palsgaard M L N, Martinez U, Blom A, Brandbyge M, Stokbro K 2019 J. Phys. Condens. Matter 32 015901Google Scholar

    [45]

    Pan H, Tang L M, Chen K Q 2022 Phys. Rev. B 105 064401Google Scholar

  • 图 1  氨基或硝基钝化前后DTB-A和DTB-B分子结示意图

    Fig. 1.  Schematic diagrams of DTB-A and DTB-B molecular junctions before and after amino or nitro passivation.

    图 2  零偏压下(a) A1和(b) B1的自旋透射谱, 蓝线和红线分别代表自旋向上和自旋向下, 费米能级在能量尺度上被设定为零

    Fig. 2.  Spin-transmission spectra of (a) A1 and (b) B1 at zero bias, the blue line and the red line represent up-spin and down-spin, respectively. The Fermi level is set at zero in the energy scale.

    图 3  A1的(a) HOMO↑对应的透射本征态和(b) LUMO↓对应的透射本征态; B1的(c) HOMO↑对应的透射本征态和(d) LUMO↓对应的透射本征态

    Fig. 3.  The transmission eigenstates of (a) HOMO↑ and (b) LUMO↓of A1; the transmission eigenstates of (c) HOMO↑and (d) LUMO↓ of B1.

    图 4  零偏压下(a) A2和(b) B2的自旋透射谱, 蓝线和红线分别代表自旋向上和自旋向下, 费米能级在能量尺度上被设定为零

    Fig. 4.  Spin-transport spectra of (a) A2 and (b) B2 at zero bias, the blue line and the red line represent up-spin and down-spin, respectively. The Fermi level is set at zero in the energy scale.

    图 5  A2的(a) HOMO↑对应的透射本征态和(b) LUMO↓对应的透射本征态; B2的(c) HOMO↑对应的透射本征态和(d) LUMO↓对应的透射本征态

    Fig. 5.  The transmission eigenstates of (a) up-spin HOMO and (b) down-spin LUMO of A2; the transmission eigenstates of (c) up-spin HOMO and (d) down-spin LUMO of B2.

    图 6  零偏压下(a) A3和(b) B3的自旋透射谱, 蓝线和红线分别代表自旋向上和自旋向下, 费米能级在能量尺度被设定为零

    Fig. 6.  Spin-transport spectra of (a) A3 and (b) B3 at zero bias, the blue line and the red line represent up-spin and down-spin, respectively. The Fermi level is set at zero in the energy scale.

    图 7  A3的(a) HOMO↑和(b) LUMO↓对应的透射本征态; B3的(c)HOMO↑和(d)LUMO↓对应的透射本征态

    Fig. 7.  The transmission eigenstates of (a) up-spin HOMO and (b) down-spin LUMO of A3; the transmission eigenstates of (c) up-spin HOMO and (d) down-spin LUMO of B3.

    图 8  零偏压下A2, B2, A3和B3在费米能级上的传输路径

    Fig. 8.  The transmission pathways of A2, B2, A3 and B3 at Fermi level under zero bias.

    图 9  (a)—(c) A1, B1, A2, B2, A3和B3分别在±1.8 V范围内的自旋电流-电压特性; (d)在同一自旋状态下, A2与B2和A3与B3随电压的变化的自旋电流之比

    Fig. 9.  (a)–(c) Spin-resolved current-voltage characteristics of A1, B1, A2, B2, A3, and B3 in the range of ±1.8 V, respectively; (d) the variation of spin-current ratios with voltages in the same spin state of A2 to B2 and A3 to B3.

    图 10  (a) A1和(b) B1在0 V, ±0.8 V 和 ±1.6 V偏压下的自旋透射谱, 蓝线和红线分别代表自旋向上和自旋向下, 黑色实线之间的区域为偏压窗口, 费米能级在能量尺度上被设定为零

    Fig. 10.  Spin-transmission spectra of (a) A1 and (b) B1 at 0 V, ±0.8 V 和 ±1.6 V. The blue line and the red line represent up-spin and down-spin, respectively. The region between the black solid lines is the bias window. The Fermi level is set at zero in the energy scale.

    图 11  (a) A2和(b) B2在0 V, ±0.8 V 和 ±1.6 V偏压下的自旋透射谱, 蓝线和红线分别代表自旋向上和自旋向下, 黑色实线之间的区域为偏压窗口, 费米能级在能量尺度上被设定为零

    Fig. 11.  Spin-transmission spectra of (a) A2 and (b) B2 at 0 V, ±0.8 V 和 ±1.6 V. The blue line and the red line represent up-spin and down-spin, respectively. The region between the black solid lines is the bias window. The Fermi level is set at zero in the energy scale.

    图 12  (a) A3和(b) B3在0 V, ±0.4 V 和 ±1.2 V偏压下的自旋透射谱, 蓝线和红线分别代表自旋向上和自旋向下, 黑色实线之间的区域为偏压窗口, 费米能级在能量尺度上被设定为零

    Fig. 12.  Spin-transmission spectra of (a) A3 and (b) B3 at 0 V, ±0.4 V 和 ±1.2 V. The blue line and the red line represent up-spin and down-spin, respectively. The region between the black solid lines is the bias window. The Fermi level is set at zero in the energy scale.

  • [1]

    Aviram A, Ratner M A 1974 Chem. Phys. Lett. 29 277Google Scholar

    [2]

    Perrin M L, Frisenda R, Koole M, Seldenthuis J S, Gil J A C, Valkenier H, Hummelen J C, Renaud N, Grozema F C, Thijssen J M, Dulić D, van der Zant H S J 2014 Nat. Nanotechnol. 9 830Google Scholar

    [3]

    Koley S, Chakrabarti S 2018 Chem. Eur. J. 24 5876Google Scholar

    [4]

    Sharma P, Bernard L S, Bazigos A, Magrez A, Ionescu A M 2015 ACS Nano 9 620Google Scholar

    [5]

    Kumar S, Wang Z, Davila N, Kumari N, Norris K J, Huang X, Strachan J P, Vine D, Kilcoyne A L D, Nishi Y, Williams R S 2017 Nat. Commun. 8 658Google Scholar

    [6]

    Li Z L, Sun F, Bi J J, Liu R, Suo Y Q, Fu H Y, Zhang G P, Song Y Z, Wang D, Wang C K 2019 Physica E 106 270Google Scholar

    [7]

    Liu Q, Li J J, Wu D, Deng X Q, Zhang Z H, Fan Z Q, Chen K Q 2021 Phys. Rev. B 104 045412Google Scholar

    [8]

    Komeda J, Takada K, Maeda H, Fukui N, Tsuji T, Nishihara H 2022 Chem. Eur. J. 28 e202201316Google Scholar

    [9]

    Zhao J, Zeng H, Wang D, Yao G 2020 Appl. Surf. Sci. 519 146203Google Scholar

    [10]

    Petersen M Å, Rasmussen B, Andersen N N, Sauer S P A, Nielsen M B, Beeren S R, Pittelkow M 2017 Chem. Eur. J. 23 17010Google Scholar

    [11]

    Fan Z Q, Zhang Z H, Yang S Y 2020 Nanoscale 12 21750Google Scholar

    [12]

    Jiang J, Kula M, Lu W, Luo Y 2005 Nano Lett. 5 1551Google Scholar

    [13]

    Lambert C J 2015 Chem. Soc. Rev. 44 875Google Scholar

    [14]

    Gehring P, Thijssen J M, van der Zant H S J 2019 Nat. Rev. Phys. 1 381Google Scholar

    [15]

    Manrique D Z, Huang C, Baghernejad M, Zhao X, Al-Owaedi O A, Sadeghi H, Kaliginedi V, Hong W, Gulcur M, Wandlowski T, Bryce M R, Lambert C J 2015 Nat. Commun. 6 6389Google Scholar

    [16]

    Fan Z Q, Sun W Y, Zhang Z H, Deng X Q, Tang G P, Xie H Q 2017 Carbon 122 687Google Scholar

    [17]

    Li Z L, Bi J J, Liu R, Yi X H, Fu H Y, Sun F, Wei M Z, Wang C K 2017 Chin. Phys. B 26 098508Google Scholar

    [18]

    Yang Y, Gantenbein M, Alqorashi A, Wei J, Sangtarash S, Hu D, Sadeghi H, Zhang R, Pi J, Chen L, Huang X, Li R, Liu J, Shi J, Hong W, Lambert C J, Bryce M R 2018 J. Phys. Chem. C 122 14965Google Scholar

    [19]

    Liu X, Sangtarash S, Reber D, Zhang D, Sadeghi H, Shi J, Xiao Z Y, Hong W, Lambert C J, Liu S X 2017 Angew. Chem. 129 179Google Scholar

    [20]

    Borges A, Fung E D, Ng F, Venkataraman L, Solomon G C 2016 J. Phys. Chem. Lett. 7 4825Google Scholar

    [21]

    Frisenda R, Janssen V A E C, Grozema F C, van der Zant H S J, Renaud N 2016 Nat. Chem. 8 1099Google Scholar

    [22]

    Liu R, Bi J J, Xie Z, Yin K, Wang D, Zhang G P, Xiang D, Wang C K, Li Z L 2018 Phys. Rev. Appl. 9 054023Google Scholar

    [23]

    Zhang Y P, Chen L C, Zhang Z Q, Cao J J, Tang C, Liu J, Duan L L, Huo Y, Shao X, Hong W, Zhang H L 2018 J. Am. Chem. Soc. 140 6531Google Scholar

    [24]

    Yang G, Sangtarash S, Liu Z, Li X, Sadeghi H, Tan Z, Li R, Zheng J, Dong X, Liu J, Yang Y, Shi J, Xiao Z, Zhang G, Lambert C, Hong W, Zhang D 2017 Chem. Sci. 8 7505Google Scholar

    [25]

    Li Y, Buerkle M, Li G, Rostamian A, Wang H, Wang Z, Bowler D R, Miyazaki T, Xiang L, Asai Y, Zhou G, Tao N 2019 Nat. Mater. 18 357Google Scholar

    [26]

    Huang B, Liu X, Yuan Y, Hong Z W, Zheng J F, Pei L Q, Shao Y, Li J F, Zhou X S, Chen J Z, Jin S, Mao B W 2018 J. Am. Chem. Soc. 140 17685Google Scholar

    [27]

    Bai J, Daaoub A, Sangtarash S, Li X, Tang Y, Zou Q, Sadeghi H, Liu S, Huang X, Tan Z, Liu J, Yang Y, Shi J, Mészáros G, Chen W, Lambert C, Hong W 2019 Nat. Mater. 18 364Google Scholar

    [28]

    Zheng J, Liu J, Zhuo Y, Li R, Jin X, Yang Y, Chen Z B, Shi J, Xiao Z, Hong W, Tian Z Q 2018 Chem. Sci. 9 5033Google Scholar

    [29]

    Liu J, Zhao X, Al-Galiby Q, Huang X, Zheng J, Li R, Huang C, Yang Y, Shi J, Manrique D Z, Lambert C J, Bryce M R, Hong W 2017 Angew. Chem. 129 13241Google Scholar

    [30]

    Cai S, Deng W, Huang F, Chen L, Tang C, He W, Long S, Li R, Tan Z, Liu J, Shi J, Liu Z, Xiao Z, Zhang D, Hong W 2019 Angew. Chem. Int. Ed. 58 3829Google Scholar

    [31]

    Ozawa H, Baghernejad M, Al-Owaedi O A, Kaliginedi V, Nagashima T, Ferrer J, Wandlowski T, García-Suárez V M, Broekmann P, Lambert C J, Haga M A 2016 Chem. Eur. J. 22 12732Google Scholar

    [32]

    Greenwald J E, Cameron J, Findlay N J, Fu T, Gunasekaran S, Skabara P J, Venkataraman L 2021 Nat. Nanotechnol. 16 313Google Scholar

    [33]

    Trasobares J, Vuillaume D, Théron D, Clément N 2016 Nat. Commun. 7 12850Google Scholar

    [34]

    Niu L L, Fu H Y, Suo Y Q, Liu R, Sun F, Wang S S, Zhang G P, Wang C K, Li Z L 2021 Phys. E 128 114542Google Scholar

    [35]

    Wu D, Huang L, Jia P Z, Cao X H, Fan Z Q, Zhou W X, Chen K Q 2021 Appl. Phys. Lett. 119 063503Google Scholar

    [36]

    Fan Z Q, Xie F, Jiang X W, Wei Z, Li S S 2016 Carbon 110 200Google Scholar

    [37]

    Kushmerick J 2009 Nature 462 994Google Scholar

    [38]

    Richter S, Mentovich E, Elnathan R 2018 Adv. Mater. 30 1706941Google Scholar

    [39]

    Liu R, Han Y, Sun F, Khatri G, Kwon J, Nickle C, Wang L, Wang C K, Thompson D, Li Z L, Nijhuis C A, del Barco E 2022 Adv. Mater. 34 2202135Google Scholar

    [40]

    Fan Z Q, Sun W Y, Jiang X W, Zhang Z H, Deng X Q, Tang G P, Xie H Q, Long M Q 2017 Carbon 113 18Google Scholar

    [41]

    Zhou Q, Yang Y, Ni J, Li Z, Zhang Z 2010 Nano Res. 3 423Google Scholar

    [42]

    Baghernejad M, Van Dyck C, Bergfield J, Levine D R, Gubicza A, Tovar J D, Calame M, Broekmann P, Hong W 2019 Chem. Eur. J. 25 15141Google Scholar

    [43]

    Büttiker M, Imry Y, Landauer R, Pinhas S 1985 Phys. Rev. B 31 6207Google Scholar

    [44]

    Smidstrup S, Markussen T, Vancraeyveld P, Wellendorff J, Schneider J, Gunst T, Verstichel B, Stradi D, Khomyakov P A, Vej-Hansen U G, Lee M E, Chill S T, Rasmussen F, Penazzi G, Corsetti F, Ojanperä A, Jensen K, Palsgaard M L N, Martinez U, Blom A, Brandbyge M, Stokbro K 2019 J. Phys. Condens. Matter 32 015901Google Scholar

    [45]

    Pan H, Tang L M, Chen K Q 2022 Phys. Rev. B 105 064401Google Scholar

  • [1] 李佳锦, 刘乾, 伍丹, 邓小清, 张振华, 范志强. 蒽二噻吩分子连接铁磁锯齿边碳化硅纳米带的巨幅度自旋整流. 物理学报, 2022, 71(7): 078501. doi: 10.7498/aps.71.20212193
    [2] 崔兴倩, 刘乾, 范志强, 张振华. 氧气分子吸附对单蒽分子器件自旋输运性质调控. 物理学报, 2020, 69(24): 248501. doi: 10.7498/aps.69.20201028
    [3] 闫瑞, 吴泽文, 谢稳泽, 李丹, 王音. 导线非共线的分子器件输运性质的第一性原理研究. 物理学报, 2018, 67(9): 097301. doi: 10.7498/aps.67.20172221
    [4] 赵文静, 文灵华. 半无限深势阱中自旋相关玻色-爱因斯坦凝聚体的量子反射与干涉. 物理学报, 2017, 66(23): 230301. doi: 10.7498/aps.66.230301
    [5] 陈伟, 陈润峰, 李永涛, 俞之舟, 徐宁, 卞宝安, 李兴鳌, 汪联辉. 基于石墨烯电极的Co-Salophene分子器件的自旋输运. 物理学报, 2017, 66(19): 198503. doi: 10.7498/aps.66.198503
    [6] 鲁桃, 王瑾, 付旭, 徐彪, 叶飞宏, 冒进斌, 陆云清, 许吉. 采用密度泛函理论与分子动力学对聚甲基丙烯酸甲酯双折射性的理论计算. 物理学报, 2016, 65(21): 210301. doi: 10.7498/aps.65.210301
    [7] 代广珍, 蒋先伟, 徐太龙, 刘琦, 陈军宁, 代月花. 密度泛函理论研究氧空位对HfO2晶格结构和电学特性影响. 物理学报, 2015, 64(3): 033101. doi: 10.7498/aps.64.033101
    [8] 贺泽龙, 白继元, 李鹏, 吕天全. T型双量子点分子Aharonov-Bohm干涉仪的电输运. 物理学报, 2014, 63(22): 227304. doi: 10.7498/aps.63.227304
    [9] 白继元, 贺泽龙, 杨守斌. 平行耦合双量子点分子A-B干涉仪的电荷及其自旋输运. 物理学报, 2014, 63(1): 017303. doi: 10.7498/aps.63.017303
    [10] 胡小颖, 田宏伟, 宋立军, 朱品文, 乔靓. Li-N, Li-2N共掺p型ZnO的第一性原理研究. 物理学报, 2012, 61(4): 047102. doi: 10.7498/aps.61.047102
    [11] 窦俊青, 康雪雅, 吐尔迪·吾买尔, 华宁, 韩英. Mn掺杂LiFePO4的第一性原理研究. 物理学报, 2012, 61(8): 087101. doi: 10.7498/aps.61.087101
    [12] 张易军, 闫金良, 赵刚, 谢万峰. Si掺杂β-Ga2O3的第一性原理计算与实验研究. 物理学报, 2011, 60(3): 037103. doi: 10.7498/aps.60.037103
    [13] 高巍, 巩水利, 朱嘉琦, 马国佳. 掺氮四面体非晶碳的第一性原理研究. 物理学报, 2011, 60(2): 027104. doi: 10.7498/aps.60.027104
    [14] 李琦, 范广涵, 熊伟平, 章勇. ZnO 极性表面及其N原子吸附机理的第一性原理研究. 物理学报, 2010, 59(6): 4170-4177. doi: 10.7498/aps.59.4170
    [15] 安义鹏, 杨传路, 王美山, 马晓光, 王德华. C20F20分子电子输运性质的第一性原理研究. 物理学报, 2010, 59(3): 2010-2015. doi: 10.7498/aps.59.2010
    [16] 刘春旭, 张继森, 刘俊业, 金光. Er3+:YAlO3晶体中Λ型四能级系统的量子相干左手性. 物理学报, 2009, 58(8): 5778-5783. doi: 10.7498/aps.58.5778
    [17] 杨冲, 杨春. Si(001)表面硅氧团簇原子与电子结构的第一性原理研究. 物理学报, 2009, 58(8): 5362-5369. doi: 10.7498/aps.58.5362
    [18] 周晶晶, 陈云贵, 吴朝玲, 郑欣, 房玉超, 高涛. 新型轻质储氢材料的第一性原理原子尺度设计. 物理学报, 2009, 58(7): 4853-4861. doi: 10.7498/aps.58.4853
    [19] 祝国梁, 疏达, 戴永兵, 王俊, 孙宝德. Si在TiAl3中取代行为的第一性原理研究. 物理学报, 2009, 58(13): 210-S215. doi: 10.7498/aps.58.210
    [20] 党宏丽, 王崇愚, 于 涛. γ-TiAl中Nb和Mo合金化效应的第一性原理研究. 物理学报, 2007, 56(5): 2838-2844. doi: 10.7498/aps.56.2838
计量
  • 文章访问数:  352
  • PDF下载量:  10
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-10-15
  • 修回日期:  2022-12-20
  • 上网日期:  2022-12-26
  • 刊出日期:  2023-03-05

二噻吩硼烷异构体分子结构测定的第一性原理研究

  • 长沙理工大学物理与电子科学学院, 柔性电子材料基因工程湖南省重点实验室, 长沙 410114
  • 通信作者: 范志强, zqfan@csust.edu.cn
    基金项目: 国家自然科学基金(批准号: 12074046, 12204066)、湖南省教育厅科学基金(批准号: 20C0039, 19C0027)、长沙理工大学青年研究人员助推计划 (批准号: 2019QJCZ022)和湖南省研究生科研创新项目(批准号: CXCLY2022141)资助的课题.

摘要: 本文在第一性原理计算基础上结合非平衡格林函数方法, 研究了量子干涉效应对连接镍电极的二噻吩硼烷(dithienoborepin, DTB)分子结自旋输运性质的影响, 并通过氨基和硝基钝化实现了对二噻吩硼烷分子异构体(DTB-A和DTB-B)的区分. 结果表明, 原始的DTB-A和DTB-B分子结在费米能级两侧都有一个自旋向上透射峰和一个自旋向下透射峰, 且两个透射峰的能量位置和高度基本相同. 因此, 原始DTB-A和DTB-B分子结的自旋向上和自旋向下电流曲线基本重合, 不能被明显区分. 然而, 研究发现量子干涉效应能不同程度地增强氨基钝化DTB-A分子结费米能级两侧分子轨道的自旋极化输运能力, 并减弱氨基钝化DTB-B分子结费米能级两侧分子轨道的自旋极化输运能力. 此外, 研究还发现量子干涉效应可以显著提高硝基钝化DTB-B分子结费米能级两侧分子轨道的自旋极化输运能力, 同时减弱硝基钝化DTB-A分子结费米能级两侧分子轨道的自旋极化输运能力. 由于量子干涉效应对氨基和硝基钝化的DTB异构体分子结自旋输运能力有不同的调制作用, 因此可以通过测量氨基和硝基钝化分子结的自旋电流值来区分DTB分子的两种异构体.

English Abstract

    • Aviram和Ratner[1]在1974年首次研究分子整流器件, 其目的之一是实现分子尺度的电子器件, 从而推动电子器件进一步的小型化. 经过半个世纪的研究和发展, 分子器件的实验和理论研究都取得了重要发现. 在研究分子结的电子输运行为时, 科学家们发现了负微分电阻[2-6]、整流[7,8]、开关[9,10]、场效应管[11,12]等特性. 与以往的微电子器件不同的是, 电子在单分子结不同能级之间传输时, 由于相互干涉, 其电导会被增强或抑制, 这就是单分子结电输运过程中的量子干涉效应[13,14]. 因此, 量子干涉对分子结的电子输运性质有明显的影响. 以往的研究表明, 锚接基团与有机分子的连接位[15-17], 掺杂原子的取代位置[18,19]以及分子结构[20]都是影响量子干涉效应的因素. 此外, 量子干涉效应还可以通过机械力[21,22]、化学环境[23,24]和电化学[25-27]等外界刺激[28-30]来调节. 因此, 利用量子干涉效应调制有机单分子的电荷输运成为分子电子学研究的热点. 在分子的不同结构位置引入基团进行钝化, 使电子在分子轨道传输过程中发生干涉, 从而改变电子在结中的透射率[31]. 通过改变分子构型来控制分子结的电导率, 还可以在分子水平上制造诸如分子整流器[32-34]、分子开关[35,36]、分子晶体管[37,38]和单分子逻辑门[39]等电子器件.

      自然界中很多分子存在同分异构体. 尽管部分异构体之间的结构差异较小, 但是性质上却有很大区别, 尤其是形成分子结后的输运性质会呈现明显差异. Fan等[40]以锯齿形石墨烯纳米带为电极, 研究氧化还原反应对蒽醌异构体自旋输运性质的影响, 发现氧化还原反应诱导的自旋向上电子的电导开关可以设计成自旋电流开关. Zhang等[23]用机械可控断裂连接技术鉴定了两个含有二酮吡咯异构体的单分子(SDPP和SPPO). SPPO分子的可逆质子化和去质子化, 会导致单分子电导发生1个数量级以上的变化, 从而达到区分两种异构体的目的. Zhou等[41]报道了联苯和一氯联苯异构体的拉曼光谱的实验和理论研究, 得到了它们的拉曼光谱特征, 并以银纳米棒为基底, 对这些化合物进行检测和识别. 以往的研究结果表明, 量子干涉效应引起的分子结电导差异是识别异构体或提高识别灵敏度的重要途径. 利用扫描隧道显微镜裂结技术, 对两种具有三配位硼中心的二噻吩硼烷(dithienoborepin, DTB)异构体(DTB-A和DTB-B)进行单分子电导测量[42]. 结果表明, 量子干涉效应能提高单分子DTB结的化学响应性. 由此启发, 本文研究了氨基和硝基钝化前后的DTB-A和DTB-B两种异构体附着在镍电极上的自旋输运性质. 结果表明, 经氨基或硝基钝化后, 两种DTB异构体的电导率会因量子干涉效应而有较大差异, 从而实现了区分DTB异构体的目的.

    • 为了探讨量子干涉效应对氨基或硝基钝化前后DTB-A和DTB-B分子自旋输运性质的影响, 选择磁性金属镍作为电极构建分子结, 如图1所示. 单个DTB-A分子或DTB-B分子通过硫原子连接到两个尺寸为3×3的镍电极表面. 分子器件分为左电极、右电极和中心散射区3个区域. 中心散射区包含四层镍电极, 从而建立了分子与电极之间的耦合、平衡费米能级和平衡时的电荷中性. 氨基或硝基通过取代硼原子上的氢原子来钝化DTB分子. 为了简化, 用异构体DTB-A和DTB-B连接到镍电极上的分子结分别命名为A1和B1. 用氨基钝化的DTB-A和DTB-B连接到镍电极上的分子结分别命名为A2和B2. 与硝基钝化的分子结DTB-A和DTB-B分别命名为A3和B3.

      图  1  氨基或硝基钝化前后DTB-A和DTB-B分子结示意图

      Figure 1.  Schematic diagrams of DTB-A and DTB-B molecular junctions before and after amino or nitro passivation.

      本文在第一性原理计算基础上结合非平衡格林函数方法研究了量子干涉效应对连接镍电极的二噻吩硼烷(dithienoborepin, DTB)分子结自旋输运性质的影响. 将交换关联势选为基于Perdew-Burke-Ernzerhof (PBE)的自旋广义梯度近似(SGGA), 所有原子的轨道基函数为双极化基组. 在自洽计算中, 选择Monkhorst-Pack的K点为3×3×50, 截断能为150 Ry (1 Ry = 13.6 eV). 器件的自旋极化电流可由朗道公式求得[43],

      $ {I_{\text{σ }}}\left( {{V_{\text{b}}}} \right) = \frac{e}{h}\int {{T_{\text{σ }}}\left( {E,{V_{\text{b}}}} \right)} \left[ {{f_{\text{L}}}\left( {E,{V_{\text{b}}}} \right) - {f_{\text{R}}}\left( {E,{V_{\text{b}}}} \right)} \right]{\text{d}}E , $

      其中, $ {I_{\text{σ }}}\left( {{V_{\text{b}}}} \right) $为器件的自旋极化电流, $ e $为电子电量, $ h $为普朗克常量, $ {V_{\text{b}}} $为左右的电极的电压差, $ {f_{\text{L}}}\left( {E, {V_{\text{b}}}} \right) $$ {f_{\text{R}}}\left( {E, {V_{\text{b}}}} \right) $分别是两个电极的电子费米-狄拉克分布函数, $ \sigma $代表自旋向上或自旋向下. $ {T_{\text{σ }}}\left( {E, {V_{\text{b}}}} \right) $是在偏压为$ {V_{\text{b}}} $, 能量为$ E $时自旋电子透射系数,

      $ {T_{\text{σ }}}\left( {E,{V_{\text{b}}}} \right) = {\text{Tr}}\left[ {{\varGamma _{\text{L}}}\left( E \right){{\rm{G}}^{\text{R}}}\left( E \right){\varGamma _{\text{R}}}\left( E \right){{\rm{G}}^{\text{A}}}\left( E \right)} \right] , $

      其中${{\rm{G}}^{\text{R}}}\left( E \right)$${{\rm{G}}^{\text{A}}}\left( E \right)$分别为散射延迟和超前格林函数, ${\varGamma _{{\text{L, R}}}} = {\text{i}}\left( {\Sigma _{{\text{L, R}}}^{\text{R}}(E) - \Sigma _{{\text{L, R}}}^{\text{A}}(E)} \right)$为展宽函数, $ \Sigma _{{\text{L, R}}}^{\text{R}}(E) $$ \Sigma _{{\text{L, R}}}^{\text{A}}(E) $是左右电极对散射区的自能. 在计算输运性质之前, 先对分子进行结构优化计算, 然后再将分子组成器件进行器件结构优化, 优化过程中使每个原子受力收敛达到 0.02 eV/Å. 以上计算均在QuantumATK (2020.09版)软件包中完成[44].

    • 图2分别显示了零偏压下A1和B1的自旋透射谱. 众所周知, 分子器件在低偏压下的电子输运性质主要取决于费米能级附近透射峰的分布及其高度. 如图2(a)所示, 费米能级两侧存在一个由最高占据分子轨道(highest occupied molecular orbital, HOMO↑)主导的自旋向上透射峰和一个由最低未占据分子轨道(lowest unoccupied molecular orbital, LUMO↓)主导的自旋向下透射峰. 两个透射峰与费米能级之间的能量间隔距离以及两个透射峰对应的高度基本相同. 在图2(b)中, B1在费米能级两侧也有一个HOMO↑主导的透射峰和一个LUMO↓主导的透射峰. 这两个透射峰的能量位置与A1对应的透射峰基本相同. B1和A1中由HOMO↑主导的透射峰高度也基本相同. 但是B1中由HOMO↑主导的透射峰展宽比A1狭窄, 导致透射峰的占据面积有一定的减小. 此外, B1中由LUMO↓主导的透射峰高度比A1对应的透射峰高, 导致相应的传输峰占据面积明显增大. 这表明DTB两种异构体分子结的自旋输运能力存在一定的差异.

      图  2  零偏压下(a) A1和(b) B1的自旋透射谱, 蓝线和红线分别代表自旋向上和自旋向下, 费米能级在能量尺度上被设定为零

      Figure 2.  Spin-transmission spectra of (a) A1 and (b) B1 at zero bias, the blue line and the red line represent up-spin and down-spin, respectively. The Fermi level is set at zero in the energy scale.

      A1和B1的透射峰差异可以用它们对应的透射本征态来解释, 如图3所示. 在图3(a)中, A1的HOMO↑对应的透射本征态在左电极和中间分子上呈离域分布. 然而, 右电极上的透射本征态是局域的. 在图3(b)中, A1的LUMO↓对应的透射本征态与其HOMO↑对应的透射本征态相差不大. 由于两个透射本征态的空间展开近似, 所以HOMO↑对应的透射峰的高度与LUMO↓对应的透射峰的高度非常接近. 在图3(c)中, B1的HOMO↑对应的透射本征态在左电极和分子的左半部分呈离域分布. 但是分子的右半部分和右电极上的透射本征态局域程度增加, 导致B1的HOMO↑对应的透射峰比A1狭窄. 如图3(d)所示, B1的LUMO↓对应的透射本征态在左电极、中间分子和右电极上的离域程度增加, 所以B1的LUMO↓对应的透射峰高度相比A1有明显增大.

      图  3  A1的(a) HOMO↑对应的透射本征态和(b) LUMO↓对应的透射本征态; B1的(c) HOMO↑对应的透射本征态和(d) LUMO↓对应的透射本征态

      Figure 3.  The transmission eigenstates of (a) HOMO↑ and (b) LUMO↓of A1; the transmission eigenstates of (c) HOMO↑and (d) LUMO↓ of B1.

      图4显示了零偏压下A2和B2的自旋透射谱. 如图4(a)所示, 氨基钝化的DTB-A分子结A2在费米能级两侧也有一个HOMO↑对应的透射峰和一个LUMO↓对应的透射峰. 两个透射峰和费米能级之间的能量间隔距离基本相同. A2的LUMO↓对应的透射峰高度略高于A1, A2的HOMO↑对应的透射高度与A1的HOMO↑对应的透射峰相比有显著提高. 在图4(b)中, 氨基钝化的DTB-B分子结B2自旋透射谱与B1相比有很大的区别. 虽然B2在费米能级两侧也有一个HOMO↑对应的透射峰和一个LUMO↓对应的透射峰, 但两个透射峰的高度远低于B1. 与B1相比, B2的LUMO↓对应的透射峰高度降低了56%, 而B2的HOMO↑对应的透射峰高度降低了68%. 这说明通过氨基钝化可明显提高DTB-A和抑制DTB-B分子的自旋输运能力.

      图  4  零偏压下(a) A2和(b) B2的自旋透射谱, 蓝线和红线分别代表自旋向上和自旋向下, 费米能级在能量尺度上被设定为零

      Figure 4.  Spin-transport spectra of (a) A2 and (b) B2 at zero bias, the blue line and the red line represent up-spin and down-spin, respectively. The Fermi level is set at zero in the energy scale.

      图5显示了A2和B2在费米能级两侧的HOMO↑对应的透射本征态和LUMO↓对应的透射本征态. 如图5(a)所示, A2的HOMO↑对应的透射本征态在左电极和中间分子上呈离域分布, 且离域程度高于A1的HOMO↑对应的透射本征态. 更重要的是, 右电极上的透射本征态也表现出良好的离域分布, 具有良好的空间扩展性. 在图5(b)中, A2的LUMO↓对应的透射本征态在左电极、中间分子和右电极处也表现出较好的离域分布. 因此, A2的HOMO↑和LUMO↓对应的透射峰高度明显高于A1中相应的透射峰高度. 在图5(c)中, B2的HOMO↑对应的透射本征态完全局域在左电极和分子左半部分. 在分子右半部分和右电极上几乎没有透射本征态的分布. 在图5(d)中, B2的LUMO↓对应的透射本征态也完全局域在左电极和分子左半部分, 右电极上仅有很少的透射本征态分布. 更重要的是, B2的HOMO↑与LUMO↓对应的透射本征态是由两种具有π相位差的状态叠加而成, 因此会发生相消干涉[44]. 这表明量子干涉效应可以削弱氨基钝化的DTB-B分子的自旋输运能力. 因此, B2的HOMO↑与LUMO↓对应的透射峰高度相较B1均大幅度降低. 量子干涉效应对氨基钝化DTB-A和氨基钝化DTB-B的不同影响效果使A2和B2的自旋输运能力存在显著差异, 从而能够实现区分DTB分子异构体的目的.

      图  5  A2的(a) HOMO↑对应的透射本征态和(b) LUMO↓对应的透射本征态; B2的(c) HOMO↑对应的透射本征态和(d) LUMO↓对应的透射本征态

      Figure 5.  The transmission eigenstates of (a) up-spin HOMO and (b) down-spin LUMO of A2; the transmission eigenstates of (c) up-spin HOMO and (d) down-spin LUMO of B2.

      图6分别显示了零偏压下A3和B3的自旋透射谱. 在图6(a)中, 硝基钝化DTB-A分子结A3的HOMO↑对应的透射峰高度相比A1的HOMO↑对应的透射峰高度有明显降低. 此外, 费米能级的右侧出现了一个由硝基诱导的新分子轨道(蓝色箭头所示). 由于新轨道与HOMO之间存在较大的能量距离, 它们之间的相互作用较弱. 尽管A3的LUMO↓对应的透射峰高度相比A1的LUMO↓对应的透射峰高度增大, 但是费米能级位置的自旋向下透射峰高度却低于A1. 在图6(b)中, 硝基钝化DTB-B分子结B3的HOMO↑对应的透射峰高度相比B1的HOMO↑对应的透射峰高度提升了53%. B3的LUMO↓对应的透射峰高度相比B1的LUMO↓对应的峰高度也提升了28%. 此外, 硝基钝化还导致B3的LUMO↓对应的透射峰向费米能级移动. 更重要的是, B3的LUMO+1↓对应的透射峰同时向费米能级移动, 并与LUMO↓对应的透射峰十分接近. 因此, 费米能级处的自旋向下透射谱高度相比B1有了很大的提高.

      图  6  零偏压下(a) A3和(b) B3的自旋透射谱, 蓝线和红线分别代表自旋向上和自旋向下, 费米能级在能量尺度被设定为零

      Figure 6.  Spin-transport spectra of (a) A3 and (b) B3 at zero bias, the blue line and the red line represent up-spin and down-spin, respectively. The Fermi level is set at zero in the energy scale.

      图7所示为A3和B3在费米能级两侧的HOMO↑和LUMO↓对应的透射本征态. 在图7(a)中, A3的HOMO↑对应的透射本征态在中间分子上呈离域分布, 但其空间扩展要比A1的HOMO↑对应的透射本征态差. 此外, 透射本征态在右电极和右锚定硫原子上几乎没有分布. 因此, A3的HOMO↑对应的透射峰高度低于A1. 然而, A3的LUMO↓对应的透射本征态在在整个分子结和硝基上呈现离域分布, 导致相应的透射高度比A1有较大的增大. 在图7(c), (d)中, B3的HOMO↑和LUMO↓对应的透射本征态在整个分子结和硝基上也呈现离域分布. 更重要的是, B3的HOMO↑与LUMO↓对应的透射本征态是由两种具有2π相位差的状态叠加而成, 从而导致干涉加强[44]. 因此, 与B1相比, 相应的透射高度增大. 这表明量子干涉效应可以大大提高硝基钝化的DTB-B分子的自旋输运能力, 从而能够实现区分DTB分子异构体的目的.

      图  7  A3的(a) HOMO↑和(b) LUMO↓对应的透射本征态; B3的(c)HOMO↑和(d)LUMO↓对应的透射本征态

      Figure 7.  The transmission eigenstates of (a) up-spin HOMO and (b) down-spin LUMO of A3; the transmission eigenstates of (c) up-spin HOMO and (d) down-spin LUMO of B3.

      为了进一步说明量子干涉效应对氨基或硝基钝化的DTB-A和DTB-B的自旋输运能力有截然不同的影响, 本文给出了零偏压下A2, B2, A3和B3在费米能级的传输路径图, 如图8所示. 传输路径(transmission pathways, TP)是一种分析分子键对输运系数贡献的工具. 通常, 该路径的主要用途是计算出电子在特定能量下的传播位置(以及方向)[45]. 每个箭头的粗细表示每对原子之间局域传输的大小, 箭头的指向表示自旋电子流的方向. 在图8(a)中, 氨基钝化只阻止了DTB-A分子上硼原子周围的透射路径. 虽然反射在左右的五环结构上, 但自旋电子可以顺利地从左电极通过分子的下半部分原子和锚定硫原子传播到右电极表面. 然而, 如图8(b)所示, 传输路径在氨基钝化的DTB-B分子中呈现不连续分布, 不能到达右侧锚定的硫原子. 这表明氨基钝化能引起相消量子干涉效应, 降低DTB-B分子的自旋输运能力. 因此, 氨基钝化的DTB-A分子的自旋输运能力大于氨基钝化的DTB-B分子. 在图8(c)中, 硝基钝化不仅阻断了DTB-A分子上硼原子周围的传输路径, 而且传输路径在DTB-A分子的其他位置也呈不连续分布. 硝基钝化显著提高了传输路径在DTB-B分子上的分布范围和迁移能力. 尽管在分子中发生反射, 但自旋电子可以从左电极平稳地穿过整个分子, 并从被锚定的硫原子向右电极的更深的原子传播. 这表明硝基钝化能诱导出相长量子干涉效应, 从而增强了DTB-B分子的自旋输运能力. 因此, 硝基钝化的DTB-B分子的自旋输运能力远远大于硝基钝化的DTB-A分子.

      图  8  零偏压下A2, B2, A3和B3在费米能级上的传输路径

      Figure 8.  The transmission pathways of A2, B2, A3 and B3 at Fermi level under zero bias.

      为了进一步研究氨基和硝基钝化对DTB-A/DTB-B分子结自旋输运性质的影响, 图9中计算了6个分子器件的自旋电流-电压特性以及在相同自旋状态下被氨基和硝基钝化的DTB-A/B分子结的自旋电流比. 如图9(a)所示, A1和B1的自旋向上和自旋向下电流随着电压的升高逐渐增大, 正负电压下的自旋向上和自旋向下电流基本相同, A1和B1的自旋-电流曲线基本一致, 不能明显区分. 因此, 通过测量自旋电流无法区分DTB分子的两种异构体. 从图9(b)可以发现, A2的自旋向上和自旋向下电流随电压的增大而迅速增大. 而B2的自旋向上和自旋向下电流随电压的增大而缓慢增大. 因此, A2的自旋向上和自旋向下电流明显大于B2, 表现出更好的导电性. 即用氨基钝化法测量DTB分子的自旋电流可以区分DTB分子的两种异构体. 从图9(c)可以看出, B3的两个自旋电流随电压的增大而迅速增大. 而A3的自旋向上和自旋向下电流随电压的增大而缓慢增大. 因此, B3的两个自旋电流明显大于A3, 表现出更好的导电性. 同样, 利用硝基钝化法测量DTB分子的自旋电流, 也可以区分DTB分子的两种异构体. 通过对比图9(b), (c)中DTB-B分子的自旋电流值可以发现, 氨基钝化能显著降低DTB-B分子的自旋电流值, 而硝基钝化能显著提高DTB-B分子的自旋电流值. 该结论可用于调控DTB-B分子的自旋输运性质.

      图  9  (a)—(c) A1, B1, A2, B2, A3和B3分别在±1.8 V范围内的自旋电流-电压特性; (d)在同一自旋状态下, A2与B2和A3与B3随电压的变化的自旋电流之比

      Figure 9.  (a)–(c) Spin-resolved current-voltage characteristics of A1, B1, A2, B2, A3, and B3 in the range of ±1.8 V, respectively; (d) the variation of spin-current ratios with voltages in the same spin state of A2 to B2 and A3 to B3.

      为了说明氨基或硝基钝化后DTB-A和DTB-B自旋电流的分辨率, 图9(d)显示了电流比R2R3随电压的变化, 其中相同自旋态下R2 = IA2/IB2, R3 = IB3/IA3. 在低偏压区, 随着电压的升高, R2-up比值先降低, 在±0.5 V后基本稳定在4. 而随着电压的升高, R2-down比值缓慢增大, 在±0.6 V后始终大于R2-up. 这说明A2和B2的自旋向下电流检测可以很好地区分氨基钝化后DTB分子的两种异构体. R3-up和R3-down比值均随电压的增大而减小. 而R3-down在±0.2 V和±0.4 V的比值分别为12.5和11.2, 超过1个数量级. 因此, ±0.4 V偏压内, 通过检测硝基钝化结构的自旋向下电流, 可以更好地区分DTB分子的两种异构体.

      为了详细解释6个分子器件自旋电流随电压变化的特性, 图10图11给出了A1, B1, A2和B2在0 V, ±0.8 V 和 ±1.6 V偏压下的自旋透射谱以及图12给出了A3和B3在0 V, ±0.4 V 和 ±1.2 V偏压下的自旋透射谱. 图10(a)显示A1器件HOMO↑对应的透射峰和LUMO↓对应的透射峰会随着正负偏压的增大而向右同幅度移动. 此外, 两个透射峰的高度也随着正负偏压的增大而逐渐降低. 在±0.8 V偏压下, 偏压窗内自旋透射谱的占据面积基本相等. 尽管在±1.6 V的偏压窗扩大了, 但是由于透射峰高度的进一步降低导致偏压窗内的自旋透射谱占据面积增加不明显. 因此, A1的自旋电流在±0.8 V 之后增幅变缓. 图10(b)显示B1器件HOMO↑对应的透射峰和LUMO↓对应的透射峰随着正负偏压增大的变化情况与A1器件一致. B1在±0.8 V和 ±1.6 V偏压窗内自旋透射谱的占据面积与A1基本相等. 因此A1和B1的自旋-电流曲线基本一致, 不能明显区分.

      图  10  (a) A1和(b) B1在0 V, ±0.8 V 和 ±1.6 V偏压下的自旋透射谱, 蓝线和红线分别代表自旋向上和自旋向下, 黑色实线之间的区域为偏压窗口, 费米能级在能量尺度上被设定为零

      Figure 10.  Spin-transmission spectra of (a) A1 and (b) B1 at 0 V, ±0.8 V 和 ±1.6 V. The blue line and the red line represent up-spin and down-spin, respectively. The region between the black solid lines is the bias window. The Fermi level is set at zero in the energy scale.

      图  11  (a) A2和(b) B2在0 V, ±0.8 V 和 ±1.6 V偏压下的自旋透射谱, 蓝线和红线分别代表自旋向上和自旋向下, 黑色实线之间的区域为偏压窗口, 费米能级在能量尺度上被设定为零

      Figure 11.  Spin-transmission spectra of (a) A2 and (b) B2 at 0 V, ±0.8 V 和 ±1.6 V. The blue line and the red line represent up-spin and down-spin, respectively. The region between the black solid lines is the bias window. The Fermi level is set at zero in the energy scale.

      图  12  (a) A3和(b) B3在0 V, ±0.4 V 和 ±1.2 V偏压下的自旋透射谱, 蓝线和红线分别代表自旋向上和自旋向下, 黑色实线之间的区域为偏压窗口, 费米能级在能量尺度上被设定为零

      Figure 12.  Spin-transmission spectra of (a) A3 and (b) B3 at 0 V, ±0.4 V 和 ±1.2 V. The blue line and the red line represent up-spin and down-spin, respectively. The region between the black solid lines is the bias window. The Fermi level is set at zero in the energy scale.

      图11(a)显示A2器件HOMO↑对应的透射峰和LUMO↓对应的透射峰会随着正负偏压的增加而向右同幅度移动. 此外, 两个透射峰的高度也随着正负偏压的增加而逐渐降低. 因此, A2的自旋电流随着正负电压增加的变化情况于A1一致. 图11(b)显示B2器件HOMO↑对应的透射峰和LUMO↓对应的透射峰随着正负偏压增加同样向右移动并逐渐降低. 由于B2在零偏压下HOMO↑对应的透射峰和LUMO↓对应的透射峰相比A2大幅降低, B2在±0.8 V 和 ±1.6 V偏压窗内的自旋透射谱占据面积与A2相比减小明显. 因此, B2的自旋向上和自旋向下电流明显小于A2. 这说明通过检测A2和B2的自旋电流可以区分氨基钝化后DTB分子的两种异构体.

      图12(a)显示A3器件HOMO↑对应的透射峰和LUMO↓对应的透射峰同样会随着正负偏压的增大而向右同幅度移动. 尽管HOMO↑对应的透射峰高度随着正负偏压的增大而逐渐降低, 但是会逐渐进入偏压窗. 在±0.4 V和±1.2 V偏压窗内, 自旋向上透射谱的占据面积大于自旋向下透射谱的占据面积. 因此, A3在±1.2 V之前的自旋向下电流始终小于自旋向上电流. 图12(b)显示B3器件HOMO↑对应的透射峰和LUMO↓对应的透射峰在±0.4 V电压之前没有随着正负偏压增大发生明显移动. 因此, 随着偏压的增大, 偏压窗内的自旋透射谱占据面积迅速增长. 此外, 在±0.4 V偏压窗内的自旋向下透射谱的占据面积大于自旋向上透射谱的占据面积. 因此, B3在±0.4 V之前的自旋向下电流始终大于自旋向上电流. 当偏压增至±1.2 V, 偏压窗内两种自旋透射谱的占据面积接近相等. B3在±0.4 V之后的两种自旋态电流曲线重合在一起. 综上所述, ±0.2 V和±0.4 V的自旋电流比值R3-down要显著大于R3-up. 因此, ±0.4 V偏压内, 通过检测硝基钝化结构的自旋向下电流, 可以更好区分DTB分子的两种异构体.

    • 本文在第一性原理计算基础上结合非平衡格林函数方法研究了量子干涉效应对连接在镍电极上的DTB分子结自旋输运性质的影响, 并通过氨基和硝基钝化实现了区分DTB异构体(DTB-A和DTB-B)的目的. 结果表明, 原始DTB-A和DTB-B分子在费米能级两侧均有一个HOMO↑对应的透射峰和一个LUMO↓对应的透射峰, 两个透射峰的能量位置和高度基本相同. 因此, 两个结构的自旋向上和自旋向下电流曲线基本重合, 从而无法简单地通过原始分子的自旋电流来明确区分DTB分子的两种异构体. 而量子干涉效应能不同程度地增强氨基钝化DTB-A分子的双轨道自旋极化输运能力, 并减弱氨基钝化DTB-B分子的双轨道自旋极化输运能力. 因此, 氨基钝化的DTB-A分子结的电流总是高于氨基钝化的DTB-B分子结的电流. 量子干涉效应能显著提高硝基钝化DTB-B分子的双轨道自旋极化输运能力, 但削弱了硝基钝化DTB-A分子的双轨道自旋极化输运能力. 因此, 硝基钝化的DTB-B分子结的电流总是高于硝基钝化的DTB-A分子结的电流. 由于量子干涉效应对DTB-A和被氨基或硝基钝化的DTB-B的自旋输运能力有不同的影响, 因此可以通过测量钝化结构的自旋电流值来区分DTB分子的两种异构体. 上述结论为今后实际制备自旋分子结及其自旋输运性能的调控提供了更多的理论指导.

参考文献 (45)

目录

    /

    返回文章
    返回