N, B原子取代调控M-OPE分子器件的量子干涉与自旋输运

彭淑平; 邓淑玲; 刘乾; 董丞骐; 范志强

doi:10.7498/aps.73.20240174

摘要

采用第一性原理计算基础上结合非平衡格林函数方法, 开展了N, B原子取代对间苯乙烯低聚物(M-OPE)分子器件量子干涉与自旋输运的调控研究. 研究结果表明N, B原子在中心苯环不同位置取代对M-OPE分子器件原有的相消量子干涉抑制程度不同. 因此, N, B原子在不同位置取代后的器件电导存在较大差异. 研究还发现B原子取代的器件自旋电流值要明显高于N原子取代的器件, 且B原子在特定位置取代后, 器件在负偏压下的自旋电流值要明显大于正偏压下的自旋电流值, 呈现显著的自旋整流效应. 本文得到的N, B原子取代对分子体系量子干涉和自旋输运调控的物理机制, 可以为杂环芳烃在分子电子学中的进一步应用提供理论指导.

关键词:

Abstract

In this paper, the first-principles method based on density functional theory and non-equilibrium Green’s function is used to investigate the modulation of quantum interference and spin transport in N and B atom substituted meta-phenylene (M-OPE) molecular devices. The zero bias spin transmission spectrum of M-OPE molecular device shows that highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) are located at higher energy positions on both sides of the Fermi level, and there is a clear transmission spectrum valley (anti resonance peak) on the right side of the Fermi level. This indicates that M-OPE molecules are typical destructive quantum interference molecular systems. Research has found that N and B atoms replace carbon atoms at positions 1, 2, and 3 on the central ring of the molecule, which suppress the original destructive quantum interference of M-OPE molecular device to different extents. The substitution of N and B atoms at position 1 has no effect on the original destructive quantum interference of M-OPE molecular device, while the substitution of N and B atoms at positions 2 and 3 significantly suppresses the original destructive quantum interference of M-OPE molecular device. Therefore, there is a significant difference in the electrical conductivity of devices with N and B atoms at different positions, with the order of electrical conductivity values being N2 > N3 > N1 and B2 > B3 > B1. In this study, it is also found that the spin current value of device with B atom substitution is significantly higher than that of device with N atom substitution. After the substitution of B atom at position 2, the spin current value of the device under negative bias is significantly greater than that under positive bias, exhibiting a significant spin rectification effect. Based on the extended curled arrow rule proposed by O’Driscoll et al. to predict the behavior of quantum interference effects, we explain the physical mechanism by which N and B protons at different positions have different effects on the suppression of quantum interference in M-OPE molecular device. The results of the quantum interference and spin transport regulation of molecular systems by the substitution of B and N atoms can provide theoretical guidance for realizing the further application of heterocyclic aromatic hydrocarbons in molecular electronics.

Keywords:

PACS:

05.50.+q	(Lattice theory and statistics)
45.70.Mg	(Granular flow: mixing, segregation and stratification)
05.65.+b	(Self-organized systems)

作者及机构信息

1.
长沙理工大学物理与电子科学学院, 长沙　410114

2.
南京工业大学环境科学与工程学院, 南京　211816

通信作者: 范志强, zqfan@csust.edu.cn

基金项目: 国家自然科学基金(批准号: 12074046)和湖南省研究生科研创新项目(批准号: CXCLY2022141)资助的课题.

Authors and contacts

1.
School of Physics and Electronic Science, Changsha University of Science and Technology, Changsha 410114, China

2.
College of Environmental Science and Engineering, Nanjing University of Technology, Nanjing 211816, China

Corresponding author: Fan Zhi-Qiang, zqfan@csust.edu.cn

Funds: Project supported by the National Natural Science Foundation of China (Grant No. 12074046) and the Hunan Provincial Postgraduate Scientific Research and Innovation Project, China (Grant No. CXCLY2022141).

文章全文

1. 引　言

随着分子电子学研究的快速发展, 对单分子器件本征电输运性质研究受到了广泛的关注. 单分子器件是指采用单分子作为桥梁连接在不同类型电极构成的复合体系, 分子作为一个丰富的电子系统, 由于其分立的轨道能级, 从而具备实现分子级器件和功能电路的潜力^[1–4]. 通过了解影响电子传输特性的因素, 可以根据需要自下而上地定制分子尺度的功能单元^[5–7]. 因此, 控制单分子器件的性质(如电导、光学和磁性)将成为科学研究的重点^[8–10]. 由于单分器件的尺寸进入纳米量级, 因此其输运性质与量子干涉造成的共振输运息息相关^[11–14]. 量子干涉效应是一种分子特有的量子输运特性, 源于电子在分立的分子轨道传输时, 发生相互干涉^[15,16]. 当电子波函数传递的相位路径一致时, 发生相长干涉, 这能有效地增强分子器件电导; 而当电子波函数传递的相位路径相差为π时, 产生相消干涉, 从而降低分子器件电导^[17].

影响分子器件中量子干涉效应的因素很多, 其中分子的几何结构是影响量子干涉效应的重要因素之一, 这导致分子结几何结构发生变化时会表现出不同的输运性能^[18–20]. 不同杂原子(O, N, S, B)取代碳原子所诱导的分子对称性变化可以调节分子中的量子干涉效应, 从而改变母体全碳化合物的性能^[21–23]. O和S杂原子对分子轨道的离域影响很小, 但会缩小分子最高占据态轨道(HOMO)能级与最低未占据态轨道(LUMO)能级的间隙, 从而导致高电导^[24]. 因为N原子的孤对电子与石墨烯纳米带的π系统不共轭, N掺杂的吡啶和嘧啶环形式的石墨烯纳米带价带和导带的能量降低, 而带隙几乎不受影响^[25]. 此外, 选用同一种杂原子取代时, 不同的取代位置对输运行为的影响也不同^[26,27].

相消量子干涉存在于具有间连接性的苯环中, 吡啶基环表现出与苯环相似的行为, 与对位和邻位连接相比, 间连接耦合情况下的相消量子干涉显著地降低分子器件的电导^[28]. 尽管具有苯环结构的分子器件中相消量子干涉无法避免, 但是如何调节相消量子干涉, 提升分子器件电导并使其展现出更丰富的电子输运特性是一个非常具有科学价值的研究方向. B, N原子同为吸电子原子且具有和C原子相近的尺寸, 所以常常在实验中被用于取代C原子从而调控体系的性质. 因此, 本文采用基于密度泛函理论结合非平衡格林函数的第一性原理计算方法, 开展B, N原子在对间苯乙烯低聚物(M-OPE)分子中心苯环不同位置取代对分子器件中量子干涉效应的影响研究. 通过揭示B, N原子取代对量子干涉效应调控的物理机制为杂环芳烃在分子电子学中的进一步应用提供理论指导.

2. 模型与方法

为了探讨N, B原子在M-OPE分子中心苯环不同位置取代后量子干涉效应对自旋输运性质的影响, 选择磁性金属钴作为电极构建分子器件, 如图1所示. 单个改造后的M-OPE分子通过硫原子连接到两个尺寸为3×3的钴电极表面. 分子器件分为左电极、右电极和中心散射区三个区域, 中心散射区包含四层钴电极. N, B原子分别依次取代中心环上1, 2, 3位置的碳原子及对应的H原子从而发生sp²杂化形成π键. 为了简化, N, B原子在1, 2, 3位置取代的器件命名为N1, N2, N3和B1, B2, B3器件.

$图 1 N或B原子取代后M-OPE分子结示意图\r\nFig. 1. Schematic diagram of M-OPE molecular junction after N or B atom substitution.$

图 1 N或B原子取代后M-OPE分子结示意图

Fig. 1. Schematic diagram of M-OPE molecular junction after N or B atom substitution.

下载: 全尺寸图片幻灯片

基于密度泛函理论和非平衡格林函数方法, 将交换关联势选为基于Perdew-Burke-Ernzerhof (PBE)的自旋广义梯度近似(SGGA), 所有原子的轨道基函数为双极化基组. 在自洽计算中, 选择Monkhorst-Pack的K点为3×3×50, 截断能为150 Ry. 器件的自旋极化电流可由朗道公式^[29]求得:

${I}_{\sigma } ({V}_{{\mathrm{b}}}) = \frac{e}{h}\int {T}_{\sigma } (E, {V}_{{\mathrm{b}}} ) \big[{f}_{{\mathrm{L}}}(E, {V}_{{\mathrm{b}}}) -{f}_{{\mathrm{R}}} (E, {V}_{{\mathrm{b}}})\big]{\mathrm{d}}E ,$

其中, ${I}_{\sigma }({V}_{{\mathrm{b}}})$ 为器件的自旋极化电流, $e$ 为电子电量, $h$ 为普朗克常量, ${V}_{{\mathrm{b}}}$ 为左右的电极的电压差, ${f}_{{\mathrm{L}}}(E, {V}_{{\mathrm{b}}})$ 和 ${f}_{{\mathrm{R}}}(E, {V}_{{\mathrm{b}}})$ 分别是两个电极的电子费米-狄拉克分布函数, $\sigma$ 代表自旋向上或自旋向下. ${T}_{\sigma }(E, {V}_{{\mathrm{b}}})$ 是在偏压为 ${V}_{{\mathrm{b}}}$ , 能量为 $E$ 时自旋电子透射系数. ${T}_{\sigma }(E, {V}_{{\mathrm{b}}})={\mathrm{T}}{\mathrm{r}}[{\varGamma }_{{\mathrm{L}}}(E){G}^{{\mathrm{R}}}(E){\varGamma }_{{\mathrm{R}}}(E){G}^{{\mathrm{A}}}(E)]$ . 其中, ${G}^{{\mathrm{R}}}(E)$ 和 ${G}^{{\mathrm{A}}}(E)$ 分别为散射延迟和超前格林函数, ${\varGamma }_{{\mathrm{L}}, {\mathrm{R}}}={\mathrm{i}}[{\Sigma }_{\text{L, R}}^{{\mathrm{R}}}(E)-{\Sigma }_{\text{L, R}}^{{\mathrm{A}}}(E)]$ 为展宽函数, ${\Sigma }_{\text{L, R}}^{{\mathrm{R}}}(E)$ 和 ${\Sigma }_{\text{L, R}}^{{\mathrm{A}}}(E)$ 是左右电极对散射区的自能. 在计算输运性质之前, 先对分子进行结构优化计算, 然后再将分子组成器件进行器件结构优化, 优化过程中使每个原子受力收敛达到0.02 eV/Å. 以上计算均在QuantumATK (2020.09版)软件包中完成^[29,30].

3. 结果与讨论

分子器件输运谱上费米能级处的透射系数反映了分子电导的大小和趋势. 对于具有相消量子干涉效应的分子系统, 在其透射谱在费米能级附近, 可以观测到明显的波谷(反共振峰), 而具有相长量子干涉效应的分子系统, 在费米能级处的透射系数变化较为平缓, 呈现出相对缓和的“U”型. 图2(a)给出了本征M-OPE分子器件的零偏压自旋透射谱. 器件HOMO与LUMO分列于费米能级两侧较高的能量位置. 费米能级右侧出现一个明显的透射谱波谷, 说明M-OPE分子是典型的相消量子干涉效应的分子系统. 当N原子在1位置取代后, N1器件的自旋透射谱相比M-OPE分子器件整体向左小幅度移动, 导致费米能级处的两种自旋态透射系数减小(电导降低), 见图2(b). 此外, 费米能级右侧的透射谱波谷仍然存在, 说明N原子在1位置取代对分子器件量子干涉效应的影响不明显.

$图 2 零偏压下(a) M-OPE和(b) N1的自旋透射谱. 红线和蓝线分别代表自旋向上和自旋向下\r\nFig. 2. Spin transmission spectra of (a) M-OPE and (b) N1 under zero bias. Red and blue lines represent spin up and spin down, respectively.$

图 2 零偏压下(a) M-OPE和(b) N1的自旋透射谱. 红线和蓝线分别代表自旋向上和自旋向下

Fig. 2. Spin transmission spectra of (a) M-OPE and (b) N1 under zero bias. Red and blue lines represent spin up and spin down, respectively.

下载: 全尺寸图片幻灯片

图3显示N2器件的自旋透射谱和HOMO位置的自旋向上和自旋向下透射本征态. 当N原子在2位置取代后, HOMO位置的两种自旋态的透射系数下降了1个数量级, 见图3(a). 从图3(b)可知, 无论是HOMO-up还是HOMO-down, 透射本征态都局域分布于分子的中间区域, 右边的苯分子没有透射本征态分布. 因此, HOMO位置的两种自旋态的透射系数出现了明显的降低. 此外, N2器件的透射谱波谷出现在费米能级左侧–1 eV能量位置. 此时, 费米能级两侧的透射系数变化较为平缓, 费米能级处的两种自旋态透射系数相比N1有很大的增加(电导升高). 这说明, N原子在2位置取代对M-OPE分子器件的相消量子干涉效应产生了抑制.

$图 3 (a) 零偏压下N2的自旋透射谱; (b) HOMO-up和HOMO-down位置的透射本征态. Isovalue的取值固定为0.35\r\nFig. 3. (a) Spin transmission spectrum of N2 under zero bias; (b) transmission eigenstates of HOMO-up and HOMO-down. The isovalue is fixed at 0.35.$

图 3 (a) 零偏压下N2的自旋透射谱; (b) HOMO-up和HOMO-down位置的透射本征态. Isovalue的取值固定为0.35

Fig. 3. (a) Spin transmission spectrum of N2 under zero bias; (b) transmission eigenstates of HOMO-up and HOMO-down. The isovalue is fixed at 0.35.

下载: 全尺寸图片幻灯片

图4给出N3器件的自旋透射谱和LUMO位置的自旋向上和自旋向下透射本征态. 当N原子在3位置取代后, LUMO位置的两种自旋态的透射系数有所减小, 见图4(a). 这是因为LUMO-up和LUMO-down的透射本征态都局域分布于分子的左边和中间区域, 右边的苯分子对透射本征态的贡献很小. 因此, LUMO位置的两种自旋态的透射系数出现了降低. 更为重要的发现是, N3器件的透射谱不再显示透射谱波谷. 此时, HOMO与LUMO之间的透射系数变化非常平缓, 费米能级处的两种自旋态透射系数相比N2又有增加(电导继续升高). 这说明, N原子在3位置取代可以明显地抑制M-OPE分子器件的相消量子干涉效应.

$图 4 (a) 零偏压下N3的自旋透射谱; (b) LUMO-up和LUMO-down位置的透射本征态. Isovalue的取值固定为0.35\r\nFig. 4. (a) Spin transmission spectrum of N3 under zero bias; (b) transmission eigenstates of LUMO-up and LUMO-down. The isovalue is fixed at 0.35.$

图 4 (a) 零偏压下N3的自旋透射谱; (b) LUMO-up和LUMO-down位置的透射本征态. Isovalue的取值固定为0.35

Fig. 4. (a) Spin transmission spectrum of N3 under zero bias; (b) transmission eigenstates of LUMO-up and LUMO-down. The isovalue is fixed at 0.35.

下载: 全尺寸图片幻灯片

接下来, 研究B原子在1, 2, 3位置取代对M-OPE分子器件量子干涉效应的影响. 图5给出了B1和B3的自旋透射谱. B1器件的自旋透射谱相比M-OPE分子器件整体向右小幅度移动, 如图5(a)所示. 此外, 透射谱波谷仍然存在并随着透射谱向右移动. 这说明B原子在1位置取代对分子器件量子干涉效应的影响也不明显. B3器件的自旋透射谱相比B1器件继续向右小幅度移动, 见图5(b). 更重要的是, B3器件的透射谱也不再显示透射谱波谷. 此时, 费米能级两侧的透射系数变化非常平缓, 费米能级处的两种自旋态透射系数相比B1有一定增加(电导升高). 因此, 无论是N原子还是B原子, 在3位置取代都会对M-OPE分子器件的相消量子干涉效应产生了较为明显的抑制.

$图 5 零偏压下(a) B1和(b) B3的自旋透射谱\r\nFig. 5. Spin transmission spectra of (a) B1 and (b) B3 under zero bias.$

图 5 零偏压下(a) B1和(b) B3的自旋透射谱

Fig. 5. Spin transmission spectra of (a) B1 and (b) B3 under zero bias.

下载: 全尺寸图片幻灯片

图6给出了B2的自旋透射谱和LUMO位置的自旋向上和自旋向下透射本征态. 当B原子在2位置取代后, LUMO位置的两种自旋态的透射系数有所减小, 见图6(a). 这是因为LUMO-up和LUMO-down的透射本征态都局域分布于分子的左边和中间区域, 右边的苯分子对透射本征态的贡献很小, 见图6(b). 因此, LUMO位置的两种自旋态的透射系数出现了降低. 更为重要的发现是, B2器件的透射谱不再显示透射谱波谷, 且费米能级处的两种自旋态透射系数在六种器件中最大(电导最高). 这说明B原子在2位置取代可以明显地抑制M-OPE分子器件的相消量子干涉效应.

$图 6 (a) 零偏压下B2的自旋透射谱; (b) LUMO-up和LUMO-down位置的透射本征态. Isovalue的取值固定为0.35\r\nFig. 6. (a) Spin transmission spectrum of B2 under zero bias; (b) transmission eigenstates of LUMO-up and LUMO-down. The isovalue is fixed at 0.35.$

图 6 (a) 零偏压下B2的自旋透射谱; (b) LUMO-up和LUMO-down位置的透射本征态. Isovalue的取值固定为0.35

Fig. 6. (a) Spin transmission spectrum of B2 under zero bias; (b) transmission eigenstates of LUMO-up and LUMO-down. The isovalue is fixed at 0.35.

下载: 全尺寸图片幻灯片

为了进一步说明B, N原子在M-OPE分子中心苯环不同位置取代对器件自旋输运性质的影响, 图7给出了6种分子器件的自旋电流-电压特性. 由于N, B原子在1位置取代对M-OPE分子器件相消量子干涉效应的影响不明显, 因此N1和B1器件在整个偏压范围内的自旋电流数值都非常小, 见图7(a)和图7(d). 尽管N原子在2位置取代对M-OPE分子器件相消量子干涉效应抑制作用不明显, 但是将透射谱波谷移动到费米能级左侧–1 eV能量位置. 因此, N2器件自旋电流随着偏压的增加而快速增大, 在整个偏压范围内的数值也明显高于N1器件, 见图7(b). 由于B原子在2位置取代明显抑制M-OPE分子器件相消量子干涉效应, 因此B2器件自旋电流随着偏压的增加也快速增大, 且数值上超过N2器件, 见图7(e). 更重要的是, B2器件在负偏压下的自旋电流值要明显大于正偏压下的自旋电流值, 呈现显著的自旋整流效应, 且自旋向下(down)的整流性质要优于自旋向上(up). 因为B原子在2位置取代后, 中心环呈不对称结构, 使得传输路径不对称, 所以导致B2器件呈现整流效应^[23,31]. 由于N, B原子在3位置取代都可以抑制M-OPE分子器件相消量子干涉效应, 因此N3和B3器件自旋电流随着偏压的增加也快速增大, 见图7(c)和图7(f).

$图 7 器件的自旋电流-电压特性　(a) N1; (b) N2; (c) N3; (d) B1; (e) B2; (f) B3\r\nFig. 7. Spin-resolved current-voltage characteristics of devices: (a) N1; (b) N2; (c) N3; (d) B1; (e) B2; (f) B3.$

图 7 器件的自旋电流-电压特性　(a) N1; (b) N2; (c) N3; (d) B1; (e) B2; (f) B3

Fig. 7. Spin-resolved current-voltage characteristics of devices: (a) N1; (b) N2; (c) N3; (d) B1; (e) B2; (f) B3.

下载: 全尺寸图片幻灯片

O’Driscoll等^[32,33]提出了预测量子干涉效应行为的扩展卷曲箭头规则(ECAR). 当分子线两侧的锚定单元分别是供体基(D)和受体基(A)时, 若D的孤对电子可以离域到A上, 则预测出现相长量子干涉效应. 如果分子线中存在吸电子基团(EWG), 则将两个锚定基团替换为D. 若D的孤对电子都可以独立地离域到同一个EWG, 则预测出现相消量子干涉效应. 图8给出了ECAR规则对B原子取代对M-OPE分子器件量子干涉效应的影响分析^[34]. B原子取代后在sp²杂化轨道保留一个孤立电子, 可作为 $\sigma$ 电子供体, 有助于促进相邻芳基之间的电子传输. 由ECAR规则可知, 当将B原子视为EWG, 两侧碳链可由D取代. B原子在1位置取代会出现相消量子干涉效应, 所以B1器件电导值很小. B原子在2, 3位置取代都会不同程度地抑制相消量子干涉效应, 从而B2和B3器件呈现较高的电导值. ECAR规则可以定性给出原子取代对量子干涉效应的调控结果, 但是无法定量给出不同原子的调控效果. 图7显示B原子取代的器件自旋电流值要明显高于N原子取代的器件. 这是因为N原子取代所引入的电子会部分占据原来的空导带, 改变分子轨道相对于费米能级的位置, 但是吡啶氮的孤立电子不与π系统共轭, 所以分子HOMO与LUMO的间隙几乎不受影响^[25]. 然而, B原子取代后的孤立电子直接掺入π骨架可以增加电子接受特性, 有效地使整个分子内的未配对电子离域, 且B原子空的p_z轨道参与了分子的LUMO^[35], 通常减小LUMO的能量从而缩减分子HOMO与LUMO的间隙. 因此, B原子取代后的器件自旋电流值比N原子取代后的器件更高.

$图 8 B原子在1, 2, 3位置取代的量子干涉效应行为预测\r\nFig. 8. Prediction of quantum interference behavior of B atom substitution at positions 1, 2, and 3.$

图 8 B原子在1, 2, 3位置取代的量子干涉效应行为预测

Fig. 8. Prediction of quantum interference behavior of B atom substitution at positions 1, 2, and 3.

下载: 全尺寸图片幻灯片

4. 总　结

利用第一性原理计算基础上结合非平衡格林函数方法, 本文开展了N, B原子取代对间苯乙烯低聚物(M-OPE)分子器件量子干涉与自旋输运的调控研究. M-OPE分子器件的零偏压自旋透射谱显示, HOMO与LUMO分列于费米能级两侧较高的能量位置, 且费米能级右侧出现一个明显的透射谱波谷(反共振峰). 这说明M-OPE分子是典型的相消量子干涉分子系统. 研究发现N, B原子分别取代分子中心环上1, 2, 3位置的碳原子, 对M-OPE分子器件原有的相消量子干涉抑制程度不同. 其中, N, B原子在1位置取代对M-OPE分子器件原有的相消量子干涉没有影响, 而N, B原子在2, 3位置取代会明显地抑制M-OPE分子器件原有的相消量子干涉. 因此, N, B原子在不同位置取代后的器件电导存在较大差异, 电导值顺序为N2 > N3 > N1和B2 > B3 > B1. 研究还发现B原子取代的器件自旋电流值要明显地高于N原子取代的器件. 另外, B原子2位置取代后, 器件在负偏压下的自旋电流值要明显大于正偏压下的自旋电流值, 呈现显著的自旋整流效应. 根据O’Driscoll等^[32,33]提出的预测量子干涉效应行为的扩展卷曲箭头规则, 解释了N, B原子在不同位置取代对M-OPE分子器件相消量子干涉产生不同影响的物理机制. 本文得到的B, N原子取代对分子体系量子干涉和自旋输运调控的结果, 可以为杂环芳烃在分子电子学中的进一步应用提供理论指导.

参考文献

[1]	闫瑞, 吴泽文, 谢稳泽, 李丹, 王音 2018 物理学报 67 097301 Google Scholar Yan R, Wu Z W, Xie W Z, Li D, Wang Y 2018 Acta Phys. Sin. 67 097301 Google Scholar
[2]	Haidar E A, Tawfik S A, Stampfl C, Hirao K, Yoshizawa K, Nakajima T, Nakajima T, Soliman K A, El-Nahas A M 2021 Adv. Theor. Simul. 4 2000203 Google Scholar
[3]	Su T A, Neupane M, Steigerwald M L, Venkataraman L, Nuckolls C 2016 Nat. Rev. Mater. 1 16002 Google Scholar
[4]	Li Y, Zhou Y, Li Y, Hong W, Li H 2022 J. Phys. Chem. C 126 6420 Google Scholar
[5]	李瑞豪, 刘俊扬, 洪文晶 2022 物理学报 71 067303 Google Scholar Li R H, Liu J Y, Hong W J 2022 Acta Phys. Sin. 71 067303 Google Scholar
[6]	Liu J, Huang X, Wang F, Hong W 2019 Acc. Chem. Res. 52 151 Google Scholar
[7]	Fan Z, Chen K 2010 Appl. Phys. Lett. 96 053509 Google Scholar
[8]	Hirai M, Tanaka N, Sakai M, Yamaguchi S 2019 Chem. Rev. 119 8291 Google Scholar
[9]	Liu Q, Li J J, Wu D, Deng X Q, Zhang Z H, Fan Z Q, Chen K Q 2021 Phys. Rev. B 104 045412 Google Scholar
[10]	Tsuji Y, Okazawa K, Kurino K, Yoshizawa K 2022 J. Phys. Chem. C 126 3244 Google Scholar
[11]	Shubin N, Emelianov A, Uspenskii Y, Gorbatsevich A 2021 Phys. Chem. Chem. Phys. 23 20854 Google Scholar
[12]	Ding Z K, Zeng Y J, Pan H, Luo N N, Zeng J, Tang L M, Chen K Q 2022 Phys. Rev. B 106 L121401 Google Scholar
[13]	Pedersen K G L, Strange M, Leijnse M, Hedegård P, Solomon G C, Paaske J 2014 Phys. Rev. B 90 125413 Google Scholar
[14]	Polakovsky A, Showman J, Valdiviezo J, Palma J L 2021 Phys. Chem. Chem. Phys. 23 1550 Google Scholar
[15]	Pan H, Ding Z K, Zeng B W, Luo N N, Zeng J, Tang L M, Chen K Q 2023 Phys. Rev. B 107 104303 Google Scholar
[16]	Qu F Y, Zhao Z H, Ren X R, Zhang S F, Wang L, Wang D 2022 Phys. Chem. Chem. Phys. 24 26795 Google Scholar
[17]	Baer R, Neuhauser D 2002 J. Am. Chem. Soc. 124 4200 Google Scholar
[18]	He R, Wang D, Luo N, Zeng J, Chen K Q, Tang L M 2023 Phys. Rev. Lett. 130 046401 Google Scholar
[19]	彭淑平, 黄旭东, 刘乾, 任鹏, 伍丹, 范志强 2023 物理学报 72 058501 Google Scholar Peng S P, Huang X D, Liu Q, Ren P, Wu D, Fan Z Q 2023 Acta Phys. Sin. 72 058501 Google Scholar
[20]	Zhang W, Zhao Z B, Tan M, Adijiang A, Zhong S, Xu X, Zhao T, Ramya E, Sun L, Zhao X, Fan Z, Xiang D 2023 Chem. Sci. 14 11456 Google Scholar
[21]	Zhang X J, Long M Q, Chen K Q, Shuai Z, Wan Q, Zou B S, Zhang Y 2009 Appl. Phys. Lett. 94 073503 Google Scholar
[22]	Yang Y, Gantenbein M, Alqorashi A, Wei J, Sangtarash S, Hu D 2018 J. Phys. Chem. C 122 14965 Google Scholar
[23]	Fan Z Q, Zhang Z H, Deng X Q, Tang G P, Chen K Q 2013 Appl. Phys. Lett. 102 023508 Google Scholar
[24]	Wang Y H, Huang H, Yu Z, Zheng J F, Shao Y, Zhou X S, Chen J Z, Li J F 2020 J. Mater. Chem. C 8 6826 Google Scholar
[25]	Wang X Y, Yao X, Narita A, Müllen K 2019 Acc. Chem. Res. 52 2491 Google Scholar
[26]	Liu X S, Sangtarash S, Reber D, Zhang D, Sadeghi H, Shi J, Xiao Z Y, Hong W J, Lambert C J, Liu S X 2017 Angew. Chem. Int. Ed. 56 173 Google Scholar
[27]	Chen Z Z, Wu S D, Lin J L, Chen L C, Cao J J, Shao X, Lambert C J, Zhang H L 2023 Adv. Electron. Mater. 9 2201024 Google Scholar
[28]	Manrique D Z, Huang C, Baghernejad M, Zhao X, Al-Owaedi O A, Sadeghi H 2015 Nat. Commun. 6 6389 Google Scholar
[29]	Büttiker M, Imry Y, Landauer R, Pinhas S 1985 Phys. Rev. B 31 6207 Google Scholar
[30]	Smidstrup S, Markussen T, Vancraeyveld P, et al. 2019 J. Phys. Condens. Matter 32 015901 Google Scholar
[31]	Deng X Q, Zhang Z H, Tang G P, Fan Z Q, Qiu M, Guo C 2012 Appl. Phys. Lett. 100 063107 Google Scholar
[32]	O’Driscoll L J, Bryce M R 2021 Nanoscale 13 1103 Google Scholar
[33]	O’Driscoll L J, Sangtarash S, Xu W, Daaoub A, Hong W J, Sadeghi H, Bryce M R 2021 J. Phys. Chem. C 125 17385 Google Scholar
[34]	Markussen T, Stadler R, Thygesen K S 2010 Nano Lett. 10 4260 Google Scholar
[35]	von Grotthuss E, John A, Kaese T, Wagner M 2018 Asian J. Org. Chem. 7 37 Google Scholar

施引文献

期刊类型引用(12)

1.	张奇峰，秦国君，崔宁，张立立. 基于漩涡特性的行人流自动扶梯拥堵形成分析. 现代电子技术. 2023(04): 166-170 . 百度学术
2.	魏娟，游磊，郭阳勇，唐志海. 基于小波神经网络的多楼层疏散模型. 系统仿真学报. 2022(02): 269-277 . 百度学术
3.	张琦，渠静. 基于前摄效应的不耐烦行为建模与双向行人流动态. 物理学报. 2022(07): 76-88 . 百度学术
4.	武鑫森，岳昊，刘秋梅，张旭，邵春福. 步行通道内行人流拉链现象的生成机理与仿真研究. 物理学报. 2021(06): 337-350 . 百度学术
5.	魏娟，郭阳勇，李忠玉，唐志海，胡周义. 考虑多群体恐惧传播的行人流模型及应急扩散研究. 安全与环境学报. 2021(02): 733-738 . 百度学术
6.	魏娟，李忠玉，游磊，郭阳勇，唐志海，胡周义. 考虑精英混沌搜索策略的室内行人疏散模型. 系统仿真学报. 2021(07): 1609-1616 . 百度学术
7.	古沐松，高朝邦，汪海鹰，李立，游磊. 后疫情视域下基于深度学习的市民休闲游憩中心应急疏散模型研究. 数字技术与应用. 2021(08): 90-92 . 百度学术
8.	陈永，张薇. 高速跟驰交通流动力学模型研究. 物理学报. 2020(06): 136-148 . 百度学术
9.	霍非舟，刘蒙，吕伟，杜煜. 考虑儿童运动特性的小学教学楼楼梯间人员疏散研究. 中国安全生产科学技术. 2020(04): 149-155 . 百度学术
10.	吴庆东. 航站楼内行人交通行为模型的理论研究. 建筑技术开发. 2020(10): 9-10 . 百度学术
11.	王邦渝. 台阶约束和缓冲平台影响的楼梯疏散仿真研究. 科技资讯. 2020(20): 1-3+6 . 百度学术
12.	孙思雨，尚华艳. 元胞传输模型在行人交通领域的应用研究. 科技与创新. 2019(14): 144-145+147 . 百度学术

其他类型引用(11)

图 1 N或B原子取代后M-OPE分子结示意图

Fig. 1. Schematic diagram of M-OPE molecular junction after N or B atom substitution.

下载: 全尺寸图片幻灯片

图 2 零偏压下(a) M-OPE和(b) N1的自旋透射谱. 红线和蓝线分别代表自旋向上和自旋向下

Fig. 2. Spin transmission spectra of (a) M-OPE and (b) N1 under zero bias. Red and blue lines represent spin up and spin down, respectively.

下载: 全尺寸图片幻灯片

图 3 (a) 零偏压下N2的自旋透射谱; (b) HOMO-up和HOMO-down位置的透射本征态. Isovalue的取值固定为0.35

Fig. 3. (a) Spin transmission spectrum of N2 under zero bias; (b) transmission eigenstates of HOMO-up and HOMO-down. The isovalue is fixed at 0.35.

下载: 全尺寸图片幻灯片

图 4 (a) 零偏压下N3的自旋透射谱; (b) LUMO-up和LUMO-down位置的透射本征态. Isovalue的取值固定为0.35

Fig. 4. (a) Spin transmission spectrum of N3 under zero bias; (b) transmission eigenstates of LUMO-up and LUMO-down. The isovalue is fixed at 0.35.

下载: 全尺寸图片幻灯片

图 5 零偏压下(a) B1和(b) B3的自旋透射谱

Fig. 5. Spin transmission spectra of (a) B1 and (b) B3 under zero bias.

下载: 全尺寸图片幻灯片

图 6 (a) 零偏压下B2的自旋透射谱; (b) LUMO-up和LUMO-down位置的透射本征态. Isovalue的取值固定为0.35

Fig. 6. (a) Spin transmission spectrum of B2 under zero bias; (b) transmission eigenstates of LUMO-up and LUMO-down. The isovalue is fixed at 0.35.

下载: 全尺寸图片幻灯片

图 7 器件的自旋电流-电压特性　(a) N1; (b) N2; (c) N3; (d) B1; (e) B2; (f) B3

Fig. 7. Spin-resolved current-voltage characteristics of devices: (a) N1; (b) N2; (c) N3; (d) B1; (e) B2; (f) B3.

下载: 全尺寸图片幻灯片

图 8 B原子在1, 2, 3位置取代的量子干涉效应行为预测

Fig. 8. Prediction of quantum interference behavior of B atom substitution at positions 1, 2, and 3.

下载: 全尺寸图片幻灯片

[1]	闫瑞, 吴泽文, 谢稳泽, 李丹, 王音 2018 物理学报 67 097301 Google Scholar Yan R, Wu Z W, Xie W Z, Li D, Wang Y 2018 Acta Phys. Sin. 67 097301 Google Scholar
[2]	Haidar E A, Tawfik S A, Stampfl C, Hirao K, Yoshizawa K, Nakajima T, Nakajima T, Soliman K A, El-Nahas A M 2021 Adv. Theor. Simul. 4 2000203 Google Scholar
[3]	Su T A, Neupane M, Steigerwald M L, Venkataraman L, Nuckolls C 2016 Nat. Rev. Mater. 1 16002 Google Scholar
[4]	Li Y, Zhou Y, Li Y, Hong W, Li H 2022 J. Phys. Chem. C 126 6420 Google Scholar
[5]	李瑞豪, 刘俊扬, 洪文晶 2022 物理学报 71 067303 Google Scholar Li R H, Liu J Y, Hong W J 2022 Acta Phys. Sin. 71 067303 Google Scholar
[6]	Liu J, Huang X, Wang F, Hong W 2019 Acc. Chem. Res. 52 151 Google Scholar
[7]	Fan Z, Chen K 2010 Appl. Phys. Lett. 96 053509 Google Scholar
[8]	Hirai M, Tanaka N, Sakai M, Yamaguchi S 2019 Chem. Rev. 119 8291 Google Scholar
[9]	Liu Q, Li J J, Wu D, Deng X Q, Zhang Z H, Fan Z Q, Chen K Q 2021 Phys. Rev. B 104 045412 Google Scholar
[10]	Tsuji Y, Okazawa K, Kurino K, Yoshizawa K 2022 J. Phys. Chem. C 126 3244 Google Scholar
[11]	Shubin N, Emelianov A, Uspenskii Y, Gorbatsevich A 2021 Phys. Chem. Chem. Phys. 23 20854 Google Scholar
[12]	Ding Z K, Zeng Y J, Pan H, Luo N N, Zeng J, Tang L M, Chen K Q 2022 Phys. Rev. B 106 L121401 Google Scholar
[13]	Pedersen K G L, Strange M, Leijnse M, Hedegård P, Solomon G C, Paaske J 2014 Phys. Rev. B 90 125413 Google Scholar
[14]	Polakovsky A, Showman J, Valdiviezo J, Palma J L 2021 Phys. Chem. Chem. Phys. 23 1550 Google Scholar
[15]	Pan H, Ding Z K, Zeng B W, Luo N N, Zeng J, Tang L M, Chen K Q 2023 Phys. Rev. B 107 104303 Google Scholar
[16]	Qu F Y, Zhao Z H, Ren X R, Zhang S F, Wang L, Wang D 2022 Phys. Chem. Chem. Phys. 24 26795 Google Scholar
[17]	Baer R, Neuhauser D 2002 J. Am. Chem. Soc. 124 4200 Google Scholar
[18]	He R, Wang D, Luo N, Zeng J, Chen K Q, Tang L M 2023 Phys. Rev. Lett. 130 046401 Google Scholar
[19]	彭淑平, 黄旭东, 刘乾, 任鹏, 伍丹, 范志强 2023 物理学报 72 058501 Google Scholar Peng S P, Huang X D, Liu Q, Ren P, Wu D, Fan Z Q 2023 Acta Phys. Sin. 72 058501 Google Scholar
[20]	Zhang W, Zhao Z B, Tan M, Adijiang A, Zhong S, Xu X, Zhao T, Ramya E, Sun L, Zhao X, Fan Z, Xiang D 2023 Chem. Sci. 14 11456 Google Scholar
[21]	Zhang X J, Long M Q, Chen K Q, Shuai Z, Wan Q, Zou B S, Zhang Y 2009 Appl. Phys. Lett. 94 073503 Google Scholar
[22]	Yang Y, Gantenbein M, Alqorashi A, Wei J, Sangtarash S, Hu D 2018 J. Phys. Chem. C 122 14965 Google Scholar
[23]	Fan Z Q, Zhang Z H, Deng X Q, Tang G P, Chen K Q 2013 Appl. Phys. Lett. 102 023508 Google Scholar
[24]	Wang Y H, Huang H, Yu Z, Zheng J F, Shao Y, Zhou X S, Chen J Z, Li J F 2020 J. Mater. Chem. C 8 6826 Google Scholar
[25]	Wang X Y, Yao X, Narita A, Müllen K 2019 Acc. Chem. Res. 52 2491 Google Scholar
[26]	Liu X S, Sangtarash S, Reber D, Zhang D, Sadeghi H, Shi J, Xiao Z Y, Hong W J, Lambert C J, Liu S X 2017 Angew. Chem. Int. Ed. 56 173 Google Scholar
[27]	Chen Z Z, Wu S D, Lin J L, Chen L C, Cao J J, Shao X, Lambert C J, Zhang H L 2023 Adv. Electron. Mater. 9 2201024 Google Scholar
[28]	Manrique D Z, Huang C, Baghernejad M, Zhao X, Al-Owaedi O A, Sadeghi H 2015 Nat. Commun. 6 6389 Google Scholar
[29]	Büttiker M, Imry Y, Landauer R, Pinhas S 1985 Phys. Rev. B 31 6207 Google Scholar
[30]	Smidstrup S, Markussen T, Vancraeyveld P, et al. 2019 J. Phys. Condens. Matter 32 015901 Google Scholar
[31]	Deng X Q, Zhang Z H, Tang G P, Fan Z Q, Qiu M, Guo C 2012 Appl. Phys. Lett. 100 063107 Google Scholar
[32]	O’Driscoll L J, Bryce M R 2021 Nanoscale 13 1103 Google Scholar
[33]	O’Driscoll L J, Sangtarash S, Xu W, Daaoub A, Hong W J, Sadeghi H, Bryce M R 2021 J. Phys. Chem. C 125 17385 Google Scholar
[34]	Markussen T, Stadler R, Thygesen K S 2010 Nano Lett. 10 4260 Google Scholar
[35]	von Grotthuss E, John A, Kaese T, Wagner M 2018 Asian J. Org. Chem. 7 37 Google Scholar

[1]	严岩, 孙峰, 羊志, 孔程昱, 葛云龙, 陈登辉, 邱帅, 李宗良. 金电极对偶氮苯分子结的结构及其电输运性质的力学调控作用. 物理学报, 2024, 73(8): 088502. doi: 10.7498/aps.73.20231999
[2]	张明媚, 郭亚涛, 付旭日, 李梦蕾, 任宝藏, 郑军, 袁瑞玚. 铁磁电极单层二硫化钼纳米带量子结构中的自旋开关效应和巨磁阻. 物理学报, 2023, 72(15): 157202. doi: 10.7498/aps.72.20230483
[3]	秦志杰, 张惠晴, 张广平, 任俊峰, 王传奎, 胡贵超, 邱帅. 通过边缘修饰在非磁性石墨烯基单分子结中引入自旋的理论研究. 物理学报, 2023, 72(13): 138504. doi: 10.7498/aps.72.20230267
[4]	彭淑平, 黄旭东, 刘乾, 任鹏, 伍丹, 范志强. 二噻吩硼烷异构体分子结构测定的第一性原理研究. 物理学报, 2023, 72(5): 058501. doi: 10.7498/aps.72.20221973
[5]	李佳锦, 刘乾, 伍丹, 邓小清, 张振华, 范志强. 蒽二噻吩分子连接铁磁锯齿边碳化硅纳米带的巨幅度自旋整流. 物理学报, 2022, 71(7): 078501. doi: 10.7498/aps.71.20212193
[6]	崔兴倩, 刘乾, 范志强, 张振华. 氧气分子吸附对单蒽分子器件自旋输运性质调控. 物理学报, 2020, 69(24): 248501. doi: 10.7498/aps.69.20201028
[7]	闫瑞, 吴泽文, 谢稳泽, 李丹, 王音. 导线非共线的分子器件输运性质的第一性原理研究. 物理学报, 2018, 67(9): 097301. doi: 10.7498/aps.67.20172221
[8]	赵文静, 文灵华. 半无限深势阱中自旋相关玻色-爱因斯坦凝聚体的量子反射与干涉. 物理学报, 2017, 66(23): 230301. doi: 10.7498/aps.66.230301
[9]	陈伟, 陈润峰, 李永涛, 俞之舟, 徐宁, 卞宝安, 李兴鳌, 汪联辉. 基于石墨烯电极的Co-Salophene分子器件的自旋输运. 物理学报, 2017, 66(19): 198503. doi: 10.7498/aps.66.198503
[10]	鲁桃, 王瑾, 付旭, 徐彪, 叶飞宏, 冒进斌, 陆云清, 许吉. 采用密度泛函理论与分子动力学对聚甲基丙烯酸甲酯双折射性的理论计算. 物理学报, 2016, 65(21): 210301. doi: 10.7498/aps.65.210301
[11]	代广珍, 蒋先伟, 徐太龙, 刘琦, 陈军宁, 代月花. 密度泛函理论研究氧空位对HfO2晶格结构和电学特性影响. 物理学报, 2015, 64(3): 033101. doi: 10.7498/aps.64.033101
[12]	贺泽龙, 白继元, 李鹏, 吕天全. T型双量子点分子Aharonov-Bohm干涉仪的电输运. 物理学报, 2014, 63(22): 227304. doi: 10.7498/aps.63.227304
[13]	白继元, 贺泽龙, 杨守斌. 平行耦合双量子点分子A-B干涉仪的电荷及其自旋输运. 物理学报, 2014, 63(1): 017303. doi: 10.7498/aps.63.017303
[14]	窦俊青, 康雪雅, 吐尔迪·吾买尔, 华宁, 韩英. Mn掺杂LiFePO4的第一性原理研究. 物理学报, 2012, 61(8): 087101. doi: 10.7498/aps.61.087101
[15]	高巍, 巩水利, 朱嘉琦, 马国佳. 掺氮四面体非晶碳的第一性原理研究. 物理学报, 2011, 60(2): 027104. doi: 10.7498/aps.60.027104
[16]	张易军, 闫金良, 赵刚, 谢万峰. Si掺杂β-Ga2O3的第一性原理计算与实验研究. 物理学报, 2011, 60(3): 037103. doi: 10.7498/aps.60.037103
[17]	李琦, 范广涵, 熊伟平, 章勇. ZnO 极性表面及其N原子吸附机理的第一性原理研究. 物理学报, 2010, 59(6): 4170-4177. doi: 10.7498/aps.59.4170
[18]	安义鹏, 杨传路, 王美山, 马晓光, 王德华. C20F20分子电子输运性质的第一性原理研究. 物理学报, 2010, 59(3): 2010-2015. doi: 10.7498/aps.59.2010
[19]	周晶晶, 陈云贵, 吴朝玲, 郑欣, 房玉超, 高涛. 新型轻质储氢材料的第一性原理原子尺度设计. 物理学报, 2009, 58(7): 4853-4861. doi: 10.7498/aps.58.4853
[20]	祝国梁, 疏达, 戴永兵, 王俊, 孙宝德. Si在TiAl3中取代行为的第一性原理研究. 物理学报, 2009, 58(13): 210-S215. doi: 10.7498/aps.58.210

期刊类型引用(12)

1.	张奇峰，秦国君，崔宁，张立立. 基于漩涡特性的行人流自动扶梯拥堵形成分析. 现代电子技术. 2023(04): 166-170 . 百度学术
2.	魏娟，游磊，郭阳勇，唐志海. 基于小波神经网络的多楼层疏散模型. 系统仿真学报. 2022(02): 269-277 . 百度学术
3.	张琦，渠静. 基于前摄效应的不耐烦行为建模与双向行人流动态. 物理学报. 2022(07): 76-88 . 百度学术
4.	武鑫森，岳昊，刘秋梅，张旭，邵春福. 步行通道内行人流拉链现象的生成机理与仿真研究. 物理学报. 2021(06): 337-350 . 百度学术
5.	魏娟，郭阳勇，李忠玉，唐志海，胡周义. 考虑多群体恐惧传播的行人流模型及应急扩散研究. 安全与环境学报. 2021(02): 733-738 . 百度学术
6.	魏娟，李忠玉，游磊，郭阳勇，唐志海，胡周义. 考虑精英混沌搜索策略的室内行人疏散模型. 系统仿真学报. 2021(07): 1609-1616 . 百度学术
7.	古沐松，高朝邦，汪海鹰，李立，游磊. 后疫情视域下基于深度学习的市民休闲游憩中心应急疏散模型研究. 数字技术与应用. 2021(08): 90-92 . 百度学术
8.	陈永，张薇. 高速跟驰交通流动力学模型研究. 物理学报. 2020(06): 136-148 . 百度学术
9.	霍非舟，刘蒙，吕伟，杜煜. 考虑儿童运动特性的小学教学楼楼梯间人员疏散研究. 中国安全生产科学技术. 2020(04): 149-155 . 百度学术
10.	吴庆东. 航站楼内行人交通行为模型的理论研究. 建筑技术开发. 2020(10): 9-10 . 百度学术
11.	王邦渝. 台阶约束和缓冲平台影响的楼梯疏散仿真研究. 科技资讯. 2020(20): 1-3+6 . 百度学术
12.	孙思雨，尚华艳. 元胞传输模型在行人交通领域的应用研究. 科技与创新. 2019(14): 144-145+147 . 百度学术

其他类型引用(11)

计量

文章访问数: 2657
PDF下载量: 86
被引次数: 23

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

搜索

留言板

N, B原子取代调控M-OPE分子器件的量子干涉与自旋输运