本文采用第一性原理计算方法研究了锯齿型延伸对扶手椅边缘三角形石墨烯量子点磁性的物理影响。研究结果表明,当向三角形石墨烯量子点的扶手椅边缘添加两个或三个锯齿型延伸时,大部分修饰结构的基态都符合Lieb定理,但少部分修饰结构违背了Lieb定理,后者可以归因于最高占据分子轨道和最低未占据分子轨道之间较小的能量差。此外,我们还发现添加锯齿型延伸后,一部分修饰结构的基态仍然保持低自旋态,而另一部分结构的基态从低自旋态变成了高自旋态。我们的发现对于调控具有扶手椅边缘石墨烯量子点的磁学特性具有重要的意义。 Using first-principles calculations, we investigate the effect of zigzag extensions on the magnetic properties of triangular graphene quantum dot with armchair edges. The results show that when two or three zigzag extensions are added to the armchair edges of triangular graphene quantum dot, the ground states of most modified structures conform to the Lieb’s theorem, while some structures violate the Lieb’s theorem, which can be attributed to the small energy difference between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO). In addition, we find that the ground states of partial modified structures remain in the low-spin state, while the others change from a low-spin state to a high-spin state. Our findings have important im-plications for manipulating the magnetism of graphene quantum dot with armchair edges.
本文采用第一性原理计算方法研究了锯齿型延伸对扶手椅边缘三角形石墨烯量子点磁性的物理影响。研究结果表明,当向三角形石墨烯量子点的扶手椅边缘添加两个或三个锯齿型延伸时,大部分修饰结构的基态都符合Lieb定理,但少部分修饰结构违背了Lieb定理,后者可以归因于最高占据分子轨道和最低未占据分子轨道之间较小的能量差。此外,我们还发现添加锯齿型延伸后,一部分修饰结构的基态仍然保持低自旋态,而另一部分结构的基态从低自旋态变成了高自旋态。我们的发现对于调控具有扶手椅边缘石墨烯量子点的磁学特性具有重要的意义。
第一性原理,三角形石墨烯量子点,Lieb定理,锯齿型延伸,磁性
Xichao Liu1, Shichao Fang1, Yun Gao2*, Zhongbing Huang1*
1Faculty of Physics and Electronic Technology, Hubei University, Wuhan Hubei
2School of Materials Science and Engineering, Hubei University, Wuhan Hubei
Received: Apr. 25th, 2022; accepted: May 23rd, 2022; published: May 31st, 2022
Using first-principles calculations, we investigate the effect of zigzag extensions on the magnetic properties of triangular graphene quantum dot with armchair edges. The results show that when two or three zigzag extensions are added to the armchair edges of triangular graphene quantum dot, the ground states of most modified structures conform to the Lieb’s theorem, while some structures violate the Lieb’s theorem, which can be attributed to the small energy difference between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO). In addition, we find that the ground states of partial modified structures remain in the low-spin state, while the others change from a low-spin state to a high-spin state. Our findings have important implications for manipulating the magnetism of graphene quantum dot with armchair edges.
Keywords:First Principles, Triangular Graphene Quantum Dot, Lieb’s Theorem, Zigzag Extension, Magnetism
Copyright © 2022 by author(s) and beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
在过去的十几年中,石墨烯基材料由于其新颖的电、磁、力学和热性能,因而在凝聚态物理和材料科学领域得到了广泛的研究 [
石墨烯量子点是人们重点关注的一类碳基纳米结构,它在空间的三个方向上都对电子的运动作了量子限制,由此产生的新奇物性对扩展石墨烯的应用具有重要的意义 [
最近,Rizzo等人通过向扶手椅边缘石墨烯纳米带引入锯齿型延伸来改变其性能,这种不引入除碳以外其它原子的方法激发了科研人员极大的兴趣 [
图1(a)定义了扶手椅边缘有11个碳环的三角形石墨烯量子点,红色球代表A类碳原子,蓝色球代表B类原子,根据量子点中两类碳原子的分布情况,A类碳原子数目记为NA,B类碳原子数目记为NB。对于该类标准的扶手椅边缘三角形石墨烯量子点,根据Lieb定理,晶格的基态总自旋为S=|NA − NB|/2 [
图1. (a) 扶手椅边缘三角形石墨烯量子点,红色球代表A类原子,数目记为NA,蓝色球代表B类原子,数目记为NB。(b) 修饰结构S2MR-S3ML
在本文的计算中,我们采取了基于密度泛函理论的第一性原理计算方法,并且采用广义梯度近似下的PBE交换关联泛函,K点选择1 × 1 × 1 [
我们首先讨论在三角形石墨烯量子点的扶手椅边缘上添加两个锯齿型延伸并且锯齿型延伸指向扶手椅边缘右侧的情况(此时,总自旋S = 1),该结构可称为二聚体(Dimer)构型。表1列出了总自旋S = 1的结构名称(SN)、结构(Structure)、单态与三重态的能量差(ΔE = ES = 0 − ES = 1)、以及相应的基态(GS)和Lieb定理态(LTS)。根据Lieb定理,当S = 0时,修饰结构应处于低自旋态;当S = 1时,修饰结构应处于高自旋态。从表1中我们可以看出D1~D9的基态均为三重态,说明这九种修饰结构均处于符合Lieb定理的高自旋态。
SN | Structure | ΔE = ES = 0 − ES = 1 (Ha) | GS | LTS |
---|---|---|---|---|
D1 | S1TR-S2TR | 0.000095 | Triplet | Triplet |
D2 | S1TR-S2MR | 0.000162 | Triplet | Triplet |
D3 | S1TR-S2HR | 0.000129 | Triplet | Triplet |
D4 | S1MR-S2TR | 0.005871 | Triplet | Triplet |
D5 | S1MR-S2MR | 0.000862 | Triplet | Triplet |
D6 | S1MR-S2HR | 0.000255 | Triplet | Triplet |
D7 | S1HR-S2TR | 0.003695 | Triplet | Triplet |
D8 | S1HR-S2MR | 0.000918 | Triplet | Triplet |
D9 | S1HR-S2HR | 0.000181 | Triplet | Triplet |
表1. 总自旋S = 1时,九种二聚体的名称、结构、单态与三重态的能量差、以及结构的基态和Lieb定理态
我们选取上述九种修饰结构中磁性最稳定的D4和D7,分别计算了它们的自旋密度和态密度,结果如图2所示。图2(a)和图2(b)分别为D4和D7的态密度,上方的黑色实线代表自旋向上的态密度,下方的红色实线代表自旋向下的态密度;图2(c)和图2(d)分别给出D4和D7的自旋密度分布,黄色部分代表自旋向上的电子云,绿色部分代表自旋向下的电子云,褐色球代表碳原子,粉色球代表氢原子。由图2(a)和图2(b)可以看出,D4和D7在费米能级附近自旋向上和自旋向下的态密度呈不对称分布,结合它的三重态基态可以确定D4和D7处于高自旋态。从图2(c)和图2(d)可以看出,自旋密度主要集中在锯齿型延伸附近的碳原子上,而且两个锯齿型延伸上的自旋密度都是同一方向排列。
图2. (a)和(b)分别为D4和D7的态密度,黑色实线代表自旋向上的态密度,红色实线代表自旋向下的态密度;(c)和(d)分别为D4和D7的自旋密度分布,黄色部分代表自旋向上的电子云,绿色部分代表自旋向下的电子云,褐色球代表碳原子,粉色球代表氢原子
接下来,我们讨论添加两个锯齿型延伸且总自旋S = 0的情况。表2列出了总自旋S = 0的结构名称、结构、单态与三重态的能量差、以及相应的基态和Lieb定理态。此时,根据Lieb定理,修饰结构应处于低自旋状态。从表2中可以看到计算结果与Lieb定理预测的结果并非完全一致,其中D10、D11、D12、D15、D18的基态均为三重态。
SN | Structures | ΔE = ES = 0 − ES = 1 (Ha) | GS | LTS |
---|---|---|---|---|
D10 | S1TR-S2TL | 0.000917 | Triplet | Singlet |
D11 | S1TR-S2ML | 0.001302 | Triplet | Singlet |
D12 | S1TR-S2HL | 0.005505 | Triplet | Singlet |
D13 | S1MR-S2TL | −0.008295 | Singlet | Singlet |
D14 | S1MR-S2ML | −0.005824 | Singlet | Singlet |
D15 | S1MR-S2HL | 0.001329 | Triplet | Singlet |
D16 | S1HR-S2TL | −0.021635 | Singlet | Singlet |
D17 | S1HR-S2ML | −0.008326 | Singlet | Singlet |
D18 | S1HR-S2HR | 0.000866 | Triplet | Singlet |
表2. 总自旋S = 0时,九种二聚体的名称、结构、单态与三重态的能量差、以及相应的基态和Lieb定理态
作为对比,我们选取D12和D17并分别计算它们的态密度和自旋密度,结果如图3所示。图3(a)和图3(b)分别为D12和D17的态密度,上方的黑色实线代表自旋向上的态密度,下方的红色实线代表自旋向下的态密度;图3(c)和图3(d)分别显示D12和D17的自旋密度分布,黄色部分代表自旋向上的电子云,绿色部分代表自旋向下的电子云,褐色球代表碳原子,粉色球代表氢原子。由图3(a)可以看出,D12在费米能级附近自旋向上和自旋向下的态密度具有非对称性,再结合它的三重态基态说明它处于高自旋状态;从图3(b)中可以看出,D17在费米能级附近两个态密度呈对称分布,并且它的基态是单态,因此D17处于低自旋态。从图3(c)中可以看出D12的自旋密度主要集中在所添加的锯齿型延伸附近的碳原子上,且两个锯齿型延伸上的自旋密度为同方向排列;从图3(d)中可以看到D17的自旋密度同样集中在锯齿型延伸附近的碳原子上,但是两个锯齿型延伸处的自旋密度反向排列。
对于S = 0的所有修饰结构,我们统计了它们的HOMO和LUMO的能量差,其中违背Lieb定理的D10、D11、D12、D15以及D18对应的能量差分别为−0.163 eV、−0.131 eV、−0.05 eV、−0.131 eV和−0.164 eV,而符合Lieb定理的D13、D14、D16以及D17相应的能量差分别为−0.387 eV、−0.303 eV、−0.754 eV和−0.389 eV。正如之前对融合的氮杂苯低聚物的研究所指出的 [
图3. (a)和(b)分别为D12和D17的态密度,黑色实线代表自旋向上的态密度,红色实线代表自旋向下的态密度;(c)和(d)分别为D12和D17的自旋密度分布,黄色部分代表自旋向上的电子云,绿色部分代表自旋向下的电子云,褐色球代表碳原子,粉色球代表氢原子
在三角形石墨烯量子点的扶手椅边缘上添加三个锯齿型延伸时,我们简称为三聚体(Trimer),结构的总自旋为1/2或3/2。首先,我们讨论总自旋为3/2的情况。表3列出了总自旋S = 3/2的结构名称、结构、二重态与四重态的能量差、以及相应的基态和Lieb定理态。由表3可以看出,所有结构的基态均为四重态,即完全符合Lieb定理。
SN | Structures | ΔE=ES=1/2 − ES=3/2 (Ha) | GS | LTS |
---|---|---|---|---|
T1 | S1TR-S2TR-S3TR | 0.000191 | Quartet | Quartet |
T2 | S1TR-S2TR-S3MR | 0.006070 | Quartet | Quartet |
T3 | S1TR-S2TR-S3HR | 0.003841 | Quartet | Quartet |
T4 | S1TR-S2MR-S3MR | 0.001065 | Quartet | Quartet |
T5 | S1TR-S2MR-S3HR | 0.003843 | Quartet | Quartet |
T6 | S1TR-S2HR-S3MR | 0.006226 | Quartet | Quartet |
T7 | S1TR-S2HR-S3HR | 0.003763 | Quartet | Quartet |
T8 | S1MR-S2MR-S3MR | 0.005009 | Quartet | Quartet |
T9 | S1MR-S2MR-S3HR | 0.006342 | Quartet | Quartet |
T10 | S1MR-S2HR-S3HR | 0.001092 | Quartet | Quartet |
T11 | S1HR-S2HR-S3HR | 0.000360 | Quartet | Quartet |
表3. 总自旋S = 3/2时,十一种三聚体的名称、结构、二重态与四重态的能量差、以及相应的基态和Lieb定理态
我们对磁性较为稳定的T2和T6分别计算了它们的态密度和自旋密度,结果如图4所示。图4(a)和图4(b)分别为T2和T6的态密度,上方的黑色实线代表自旋向上的态密度,下方的红色实线代表自旋向下的态密度;图4(c)和图4(d)分别是T2和T6的自旋密度分布,黄色部分代表自旋向上的电子云,绿色部分代表自旋向下的电子云,褐色球代表碳原子,粉色球代表氢原子。从图4(a)和图4(b)中我们可以看到T2和T6在费米能级附近自旋向上和自旋向下的态密度都呈现出不对称分布。从图4(c)和图4(d)可以看出,自旋密度主要集中在锯齿型延伸处的碳原子附近,且三个锯齿型延伸处的自旋密度都是同方向排列。再结合T2和T6的四重态基态从而可以确定T2和T6确实处于高自旋态,这也意味着引入三个锯齿型延伸同样会使得扶手椅边缘三角形石墨烯量子点从原本的低自旋态转变为高自旋态。
接下来我们分析总自旋S = 1/2的情况,表4中列出了S = 1/2的结构名称、结构、二重态与四重态的能量差、以及相应的基态和Lieb定理态。当总自旋S = 1/2时,根据Lieb定理,修饰后的三角形石墨烯量子点应处于低自旋态,但由表4可以看出,大部分修饰结构符合Lieb定理,而少数修饰结构不符合Lieb定理,对此我们进行了进一步的研究。
图4. (a)和(b)分别为T2和T6的态密度,黑色实线代表自旋向上的态密度,红色实线代表自旋向下的态密度;(c)和(d)分别为T2和T6的自旋密度分布,黄色部分代表自旋向上的电子云,绿色部分代表自旋向下的电子云,褐色球代表碳原子,粉色球代表氢原子
SN | Structures | ΔE=ES=1/2 − ES=3/2 (Ha) | GS | LTS |
---|---|---|---|---|
T12 | S1TR-S2TR-S3TL | −0.001497 | Doublet | Doublet |
T13 | S1TR-S2TR-S3HL | −0.040528 | Doublet | Doublet |
T14 | S1TR-S2MR-S3ML | −0.007699 | Doublet | Doublet |
T15 | S1TR-S2HR-S3HL | −0.041927 | Doublet | Doublet |
T16 | S1MR-S2TR-S3ML | −0.014156 | Doublet | Doublet |
T17 | S1MR-S2MR-S3ML | −0.015901 | Doublet | Doublet |
T18 | S1MR-S2HR-S3TL | −0.021769 | Doublet | Doublet |
T19 | S1MR-S2HR-S3HL | −0.017032 | Doublet | Doublet |
T20 | S1HR-S2TR-S3HL | 0.005336 | Quartet | Doublet |
T21 | S1HR-S2MR-S3HL | 0.003999 | Quartet | Doublet |
T22 | S1HR-S2HR-S3HL | 0.003656 | Quartet | Doublet |
表4. 总自旋S = 1/2时,十一种三聚体的名称、结构、二重态与四重态的能量差、以及相应的基态和Lieb定理态
我们计算了T15和T20的态密度和自旋密度,结果如图5所示。图5(a)和图5(b)分别为T15和T20的态密度,上方的黑色实线代表自旋向上的态密度,下方的红色实线代表自旋向下的态密度;图5(c)和图5(d)分别是T15和T20的自旋密度分布,黄色部分代表自旋向上的电子云,绿色部分代表自旋向下的电子云,褐色球代表碳原子,粉色球代表氢原子。由图5(a)和图5(b)可以看出,T15和T20在费米能级附近自旋向上和自旋向下的态密度呈不对称分布。通过观察相应的自旋密度分布,可以发现T15的自旋密度主要分布在远离另外两处锯齿型延伸的另外一个锯齿型延伸上,结合它的二重态基态可知T15处于低自旋态;T20的自旋密度则集中在三个锯齿型延伸处而且都沿着同一方向排列,这清晰地表明T20处于高自旋状态。
对于S = 1/2的所有修饰结构,我们统计了它们的HOMO与LUMO的能量差。符合Lieb定理的T12-T19的能量差分别为−0.359 eV、−0.38 eV、−0.25 eV、−0.344 eV、−0.372 eV、−0.246 eV、−0.256 eV、−0.241 eV;违背Lieb定理的T20、T21、T22的能量差分别为−0.09 eV、−0.03 eV、−0.064 eV。与两个锯齿型延伸修饰的情况类似,正是HOMO和LUMO的近简并使得后面三种修饰结构从低自旋态转变为高自旋态。
综上所述,在对扶手椅边缘三角形石墨烯量子点引入三个锯齿型延伸时,绝大部分的修饰结构都符合Lieb定理,但有少数结构违背了Lieb定理,并使原本的低自旋态转变成高自旋态。
图5. (a)和(b)分别为T15和T20的态密度,黑色实线代表自旋向上的态密度,红色实线代表自旋向下的态密度;(c)和(d)分别为T15和T20的自旋密度分布,黄色部分代表自旋向上的电子云,绿色部分代表自旋向下的电子云,褐色球代表碳原子,粉色球代表氢原子
我们通过对扶手椅边缘三角形石墨烯量子点添加两个或者三个锯齿型延伸来调控其磁性。研究结果表明,大部分修饰结构的基态符合Lieb定理,但少部分结构的基态违背了Lieb定理。进一步分析表明违背Lieb定理的结构有一个共同点,即它们的最高占据分子轨道和最低未占据分子轨道之间具有较小的能量差。在三角形石墨烯量子点的扶手椅边缘上添加锯齿型延伸后,一部分修饰结构的基态仍然保持低自旋态,而另一部分结构的基态从低自旋态变成了高自旋态。本文的研究结果不仅对调控三角形石墨烯量子点的磁性具有重要的理论意义,而且对调控矩形和六边形等石墨烯量子点的磁性也有借鉴作用。
感谢湖北大学杨辉博士的讨论与交流。
本研究得到了国家自然科学基金(批准号:11674087)的资助。
刘希超,方世超,高 云,黄忠兵. 锯齿型延伸修饰扶手椅边缘三角形石墨烯量子点磁性的第一性原理研究First-Principles Study on the Magnetism of Triangular Graphene Quantum Dot with Armchair Edges Decorated by Zigzag Extensions[J]. 凝聚态物理学进展, 2022, 11(02): 38-48. https://doi.org/10.12677/CMP.2022.112005