-
当今工业生产能源主要来自化石燃料, 而化石燃烧会产生大量有毒有害气体, 其中二氧化氮(NO2)作为燃烧后的主要污染产物, 会对人类的生产生活带来严重危害[1–3]. 一方面, NO2进入臭氧层后, 会加速大气的破坏; 另一方面, 它在阳光的照射下容易转化为光化学烟雾, 在水氧环境下形成硝酸, 危害人类呼吸道健康, 也会对建筑、土壤、植物等造成腐蚀伤害[4–7]. 无论是从工业生产方面考虑, 还是从农业、环境安全考虑, 如何能有效地俘获和检测环境中NO2的问题都亟待解决. 因此寻找一种灵敏度高、选择性强的高性能NO2传感器成为当下研究的热点[8–11].
石墨烯是一种以sp2杂化呈蜂巢晶格排列组成的单层纳米结构, 其独特的结构表现出优异 的光、电、热等性能[12]. 与传统气体传感材料相比, 石墨烯较大的比表面积和较高载流子迁移率会表现出更好的灵敏度和选择性[13–15]. 这些特点使其成为最具潜力的气体传感基底材料, 拥有巨大的应用前景, 受到了广泛的关注和研究[16,17]. Zhang等[18]发现过渡金属掺杂石墨烯, 可以加强石墨烯对碳酰氯的吸附; Choudhuri等[19]研究了CO和CO2气体在非金属掺杂石墨烯上的吸附过程, 发现掺杂石墨烯对气体分子的选择性和灵敏度更好. Shukri等[20]和Shamim等[21]的研究表明, 通过引入适当的掺杂剂或缺陷对石墨烯进行改性, 可以提高石墨烯对CO, NO, SO2和其他分子的灵敏度. Zhang等[22]和Jia等[23]研究表明贵金属元素掺杂可以提高石墨烯对NO2分子的敏感性. 但在NO2分子在石墨烯上吸附的研究中, 石墨烯掺杂的元素多为贵金属. 其制备成本较高, 同时贵金属材料的储备量也相对缺乏. 因此, 要实现基于石墨烯基底的NO2气体传感的大规模应用, 就必须找到一种低成本、吸附效果好的其他元素来替代贵金属. 因此本文选择基于密度泛函理论(DFT)体系下的第一性原理平面波超软赝势方法, 研究了N掺杂、Zn掺杂、N-Zn双掺杂3种掺杂石墨烯表面吸附NO2的原子结构、Mulliken电荷分布、差分电荷密度、态密度以及光学性质. 扩展掺杂石墨烯在气体吸附和传感的应用, 同时也为基于石墨烯材料的传感研究提供理论支持.
-
本文采用了CASTEP软件包[24]计算了石墨烯表面对NO2分子的吸附. 考虑到色散相互作用、长程电子关联效应、范德瓦耳斯力等影响, 为了让计算结果更准确, 采取了广义梯近似下PBE (Perdew-Burke-Ernzerhof)交换关联泛函与DFT-D2色散校正进行几何结构优化[25,26]; 运用平面波超软赝势方法, 离子势用赝势替代描述电子与离子实的相互作用; Kohn-Sham方程和能量泛函采用自洽求解. 构建了4×4×1含有32个原子的超晶胞模型, 如图1所示. 为了避免层间相互作用, 真空层设置为20 Å, 平面波截断能设置为Ecut = 430 eV, 原子力收敛精度设置为0.03 eV/Å, 能量自洽精度收敛为5×10–7 eV/atom, 原子内应力为0.05 GPa, 布里渊区k-point设置为4×4×1. 石墨烯表面吸附NO2的吸附能Eads为
其中Egas为单个NO2气体分子的能量, Egraphene为未吸附NO2分子的石墨烯能量, Egas+graphene为整个吸附体系的总能量[27]. Eads为正时, 吸附过程为吸热过程, 体系能量变高, 吸附不易发生; Eads为负时, 吸附过程为放热过程, 体系能量降低, 结构更稳定, 吸附易发生.
-
优化后石墨烯表面吸附NO2分子的原子结构如图2所示. 石墨烯表面吸附NO2分子前后的吸附距离和吸附能如表1所示. 从表1可以看出, 所有石墨烯表面吸附NO2的吸附能均为负值, 吸附后体系的能量更低更为稳定. 根据文献[28]报道吸附能大于–0.6 eV时为物理吸附, 小于–0.8 eV时为化学吸附. 因此, 未掺杂石墨烯和N掺杂石墨烯吸附NO2为物理吸附, 其相互作用较弱. 而Zn掺杂石墨烯和N-Zn双掺杂石墨烯的吸附能从未掺杂石墨烯表面的–0.2 eV变至– 4 eV, 并且NO2与这两种表面的距离也急剧缩短, 说明Zn掺杂石墨烯和N-Zn双掺杂石墨烯表面与气体分子发生了强烈的相互作用, 吸附方式为化学吸附. 我们推测Zn掺杂石墨烯和N-Zn双掺杂石墨烯表面与NO2分子之间可能形成了化学键.
-
NO2分子的Mulliken电荷分布如表2所示. 未掺杂石墨烯和掺杂石墨烯表面都与NO2分子发生了电荷转移, 从电荷转移数量看来, N-Zn双掺杂表面 (–0.60e) > Zn掺杂表面 (–0.59e) > N掺杂表面 (–0.48e) > 未掺杂表面 (–0.25e). 所有模型都是石墨烯表面上的电子转移到NO2分子上, 这说明NO2分子是作为电子受体的存在, 被石墨烯表面所还原; 并且石墨烯表面吸附的NO2分子键长变化量与电荷转移趋势一致. 从电荷转移量和NO2分子键长变化来看, 掺杂石墨烯表面与NO2分子间的相互作用强于未掺杂石墨烯与NO2, 掺杂有利于提高石墨烯气体传感性能.
图3为石墨烯表面吸附NO2分子的总电子密度等面图 (TCD)、电荷密度差图 (CDD)和电子密度差图(ECD). 从图3(a), (c)可以清楚看出, 未掺杂石墨烯表面和N掺杂石墨烯表面与NO2分子间不存在电子轨道重叠, 表明未掺杂石墨烯和N掺杂石墨烯没有与NO2形成化学键, 它们之间的相互作用较弱. 但Zn掺杂石墨烯表面和N-Zn双掺杂石墨烯表面与NO2存在明显的电子轨道重叠, 说明这两种材料表面与气体分子之间形成了化学键, 它们之间有着强烈的相互作用. 从图3(b)可以看出, 未掺杂石墨烯表面、掺杂石墨烯表面和吸附的NO2分子之间发生了明显的电荷转移. 在石墨烯表面吸附NO2分子的过程中, NO2气体分子上发生电荷积累, 石墨烯表面的吸附位置发生了电荷耗散. 这些结论与吸附能的预测和Mulliken 电荷分析结果相吻合.
-
为了进一步讨论分析未掺杂、N掺杂、Zn掺杂、N-Zn双掺杂石墨烯表面吸附NO2分子后的电子能量情况, 本文计算了不同石墨烯表面吸附NO2分子后的态密度, 如图4所示. 可以看出, 掺杂表面的态密度较未掺杂表面整体向价带偏移; 价带的波峰数量减少, 最大的波峰的数值变大; 导带的波峰数量减少; 费米能级附近杂质能级峰值也发生改变. 这些现象表明掺杂体系的电子性质较未掺杂石墨烯可能发生了改变. 这是由掺杂原子的电子杂化作用, 如图4(d)中的N的2p, Zn的3d电子与C原子2p电子的协同作用, 改变了共掺体系的态密度. 态密度的这些变化, 特别是在费米能级附近区域的变化, 可能提高石墨烯的光电性能, 这对基于石墨烯基底的传感有利.
-
我们计算了石墨烯材料的介电函数、吸收谱和反射谱. 材料的光学性质可以用介电函数来描述, 因为电子跃迁能量远大于声子扰动的能量, 所以电子吸收光子能量由低能级向高能级跃迁时会受到辐射电场的微扰影响可以忽略. 介电函数的实部与虚部是描述材料吸收与释放光子过程的函数, 其复数形式
$ \varepsilon \left(\omega \right)={\varepsilon }_{1}\left(\omega \right)+{\mathrm{i}}{\varepsilon }_{2}\left(\omega \right) $ 描述的是电子跃迁过程, 式中$ {\varepsilon }_{1} $ 和$ {\varepsilon }_{2} $ 分别为$ {\varepsilon }_{1}={n}^{2}-{k}^{2} $ ,$ {\varepsilon }_{2}=2 nk $ [29], 其中n为折射率, k为波数. 介电函数的实部与虚部能够通过Kramers-Kronig色散关系推导, 此外吸收系数$ I\left(\omega \right) $ 和反射率$ R\left(\omega \right) $ 同样可以通过Kramers-Kronig色散关系推导得出,$ \omega $ 为角频率. 介电函数公式如下[30]:可见光范围内 (360—780 nm) 石墨烯表面吸附NO2前后的介电函数虚部如图5所示. 在360—470 nm范围内, 掺杂石墨烯表面的虚部数值略小于在未掺杂石墨烯表面. 在470—780 nm范围内, 掺杂石墨烯表面的虚部数值大于未掺杂石墨烯表面. 其中, 只有N-Zn双掺杂石墨烯虚部数值在670 nm之后改变下降趋势, 变为上升趋势.
吸收率越高说明基态电子吸收光子能量跃迁至激发态的电子数量多, 进而反映出材料对光的响应能力越好. 反射率越高说明电子吸收光子能量跃迁至激发态越多, 释放能量跃迁回低能级的电子也就越多. 石墨烯表面吸附NO2分子的吸收谱和反射谱如图6所示. 从图6(a) 可知, 在可见光范围内, 未掺杂石墨烯表面的吸收系数峰值出现在362 nm处, 约为40300 cm–1, N掺杂石墨烯的吸收系数略高于未掺杂石墨烯. Zn掺杂石墨烯和N-Zn双掺杂石墨烯在420 nm后远大于未掺杂石墨烯的吸收系数, 其中N-Zn双掺杂石墨烯的吸收系数最大, 在360 nm处达到峰值为45500 cm–1. 从图6(b) 可知, 未掺杂石墨烯表面的反射系数峰值出现在408 nm处约为0.12, N掺杂石墨烯和Zn掺杂石墨烯的反射系数略高于未掺杂石墨烯. N-Zn双掺杂石墨烯远大于未掺杂石墨烯的反射系数, N-Zn双掺杂石墨烯在780 nm处达到峰值为0.41. N-Zn双掺杂石墨烯的吸收系数和反射系数较未掺杂石墨烯峰值分别提高了约1.12倍和3.42倍.
通过对石墨烯表面掺杂前后吸附NO2分子的介电函数、吸收谱和反射谱的比较分析可以发现:吸收谱和反射谱的走势符合介电函数的分布; 在可见光范围内, N-Zn双掺杂能大幅度提高石墨烯表面吸附NO2的光学性能, 有利于基于石墨烯基底的NO2气体传感.
-
综上所述, 利用基于DFT 的第一性原理计算方法研究了石墨烯表面吸附NO2分子的微观机制和光学性能.
1) 从吸附能和电荷转移分析来看, 石墨烯表面吸附NO2分子, 吸附能均为负值, 吸附后体系更稳定; NO2分子与石墨烯表面之间发生电荷转移, NO2分子被表面还原; 未掺杂石墨烯和N掺杂石墨烯吸附方式为物理吸附, Zn掺杂和N-Zn双掺杂石墨烯为化学吸附; 掺杂石墨烯表面与NO2分子间的相互作用强于未掺杂石墨烯与NO2.
2) 从介电函数、吸收谱和反射谱来看, 在可见光范围内, 掺杂提高了石墨烯表面的光学性能. 其中, N-Zn双掺杂石墨烯的改善效果最好, 吸收系数和反射系数峰值较未掺杂石墨烯表面的峰值分别提高了约1.12倍和3.42倍.
本文的研究加深了对石墨烯表面吸附NO2过程的理解, 通过掺杂方法提高了石墨烯吸附NO2后的光学性能, 为进一步研究该材料对NO2的传感特性提供了理论依据.
利用掺杂提高石墨烯吸附二氧化氮的敏感性及光学性质的理论计算
Theoretical calculation study on enhancing the sensitivity and optical properties of graphene adsorption of nitrogen dioxide via doping
-
摘要: 为了研究NO2在未掺杂石墨烯和掺杂石墨烯(N掺杂、Zn掺杂、N-Zn双掺杂)上的吸附, 本工作采用密度泛函理论的第一性原理平面波超软赝势对其吸附过程进行模拟. 计算了石墨烯表面吸附NO2分子的吸附能、Mulliken分布、差分电荷密度、态密度和光学性质. 研究结果表明, 与未掺杂石墨烯表面相比, 掺杂石墨烯表面对吸附NO2表现出了更高的敏感性, 吸附能大小顺序为: N-Zn双掺杂表面 > Zn掺杂表面 > N掺杂表面 > 未掺杂表面. 未掺杂石墨烯和N掺杂石墨烯表面与NO2的相互作用较弱, 是物理吸附. Zn掺杂和N-Zn双掺杂石墨烯表面与NO2之间形成了化学键, 是化学吸附. 在可见光范围内, 3种掺杂方式中 N-Zn双掺杂表面对于提高石墨烯光学性能效果最佳, 其吸收系数和反射系数的峰值较未掺杂石墨烯表明分别提高了约1.12倍和3.42倍. N-Zn双掺杂石墨烯不但能增强表面与NO2的相互作用, 同时也能提高材料的光学性能, 这为基于石墨烯基底的NO2气体检测传感提供了理论支撑和实验指导.Abstract: In order to study the adsorption of NO2 on pristine graphene and doped graphene (N-doped, Zn-doped, and N-Zn co-doped), we simulate the adsorption process by applying the first-principles plane-wave ultrasoft pseudopotentials of the density-functional theory in this work. The adsorption energy, Mulliken distribution, differential charge density, density of states, and optical properties of NO2 molecules adsorbed on the graphene surface are calculated. The results show that the doped graphene surface exhibits higher sensitivity to the adsorption of NO2 compared with the pristine graphene surface, and the order of adsorption energy is as follows: N-Zn co-doped surface > Zn-doped surface > N-doped surface > pristine surface. Pristine graphene surface and N-doped graphene surface have weak interactions with and physical adsorption of NO2. Zn-doped graphene surfac and N-Zn co-doped graphene surface form chemical bonds with NO2 and are chemisorbed. In the visible range, among the three doping modes, the N-Zn co-doped surface is the most effective for improving the optical properties of graphene, with the peak absorption and reflection coefficients improved by about 1.12 and 3.42 times, respectively, compared with pristine graphene. The N-Zn co-doped graphene not only enhances the interaction between the surface and NO2, but also improves the optical properties of the material, which provides theoretical support and experimental guidance for NO2 gas detection and sensing based on graphene substrate.
-
Key words:
- NO2 /
- graphene /
- adsorption /
- first principles .
-
-
图 3 石墨烯表面吸附NO2分子 (a) 总电子密度等面图(TCD)、(b) 电荷密度差图(CDD)和(c)电子密度差图(ECD). TCD图等面设为0.02 e/Å3; CDD图等面设为0.01 e/Å3, 蓝色代表电子积累, 黄色代表电子耗尽; ECD图红色代表电荷聚集, 蓝色代表电荷耗尽
Figure 3. (a) Charge density difference (CDD), (b) total charge density (TCD), and (c) electron density difference (EDD) plots of NO2 molecules adsorbed on different graphene surfaces. The isosurfaces of TCD plots are set to 0.02 e/Å3; the isosurfaces of CDD plots are set to 0.01 e/Å3, blue represents electron accumulation and yellow represents electron depletion; in EDD plots, red represents charge accumulation and blue represents charge depletion.
表 1 石墨烯表面吸附NO2的距离和吸附能
Table 1. The distance and adsorption energy of NO2 adsorption on graphene surface.
模型 初始距离/Å 优化后距离/Å 吸附能/eV 未掺杂石墨烯 3.00 2.96 –0.22 N掺杂 3.00 2.73 –0.69 Zn掺杂 3.00 1.94 –4.80 N-Zn双掺杂 3.00 1.93 –4.97 表 2 NO2分子的Mulliken的电荷分布
Table 2. Mulliken charge distribution of NO2.
模型 种类 s电子 p电子 总电子 电荷/e 分子带电荷/e 布居数 键长/Å NO2 N 1.39 3.18 4.57 0.44 0 0.68 1.23 O 1.85 4.36 6.22 –0.22 O 1.85 4.36 6.22 –0.22 未掺杂石墨烯 N 1.46 3.21 4.67 0.33 –0.25 0.63 1.24 O 1.86 4.44 6.29 –0.29 O 1.86 4.44 6.29 –0.29 N掺杂 N 1.51 3.23 4.74 0.26 –0.48 0.60 1.25 O 1.86 4.53 6.39 –0.38 O 1.86 4.51 6.36 –0.36 Zn掺杂 N 1.48 3.40 4.88 0.12 –0.59 0.67 1.26 O 1.86 4.50 6.36 –0.36 O 1.86 4.49 6.35 –0.35 N- Zn
双掺杂N 1.49 3.40 4.88 0.11 –0.60 0.67 1.27 O 1.86 4.51 6.37 –0.37 O 1.86 4.48 6.34 –0.34 -
[1] Cooper M J, Martin R V, Hammer M S, Levelt P F, Veefkind P, Lamsal L N, Krotkov N A, Brook J R, McLinden C A 2022 Nature 601 380 doi: 10.1038/s41586-021-04229-0 [2] Gholami F, Tomas M, Gholami Z, Vakili M 2020 Sci. Total. Environ. 714 136712 doi: 10.1016/j.scitotenv.2020.136712 [3] Wang S, Liu J, Yi H, Tang X L, Yu Q J, Zhao S Z, Gao F Y, Zhou Y S, Zhong T T, Wang Y X 2022 Chemosphere 291 132917 doi: 10.1016/j.chemosphere.2021.132917 [4] Lim H, Kwon H, Kang H, Jang J E, Kwon H J 2023 Nat. Commun. 14 3114 doi: 10.1038/s41467-023-38918-3 [5] Gao Z Y, Li L L, Huang H Y, Xu S P, Yan G, Zhao M L, Ding Z 2020 Appl. Surf. Sci. 527 146939 doi: 10.1016/j.apsusc.2020.146939 [6] Geng X, Li S W, Mawella-Vithanage L, Ma T, Kilani M, Wang B W, Ma L, Hewa-Rahinduwage C C, Shafikova A, Nikolla E, Mao G Z, Brock S L, Zhang L, Luo L 2021 Nat. Commun. 12 4895 doi: 10.1038/s41467-021-25192-4 [7] 熊枫, 彭志敏, 丁艳军, 杜艳君 2022 物理学报 71 203302 doi: 10.7498/aps.71.20220975 Xiong F, Peng Z M, Ding Y J, Du Y J 2022 Acta Phys. Sin. 71 203302 doi: 10.7498/aps.71.20220975 [8] Zhao S K, Shen Y B, Zhou P F, Zhong X X, Han C, Zhao Q, Wei D Z 2019 Sens. Actuat. B-Chem. 282 917 doi: 10.1016/j.snb.2018.11.142 [9] Choi M S, Kim M Y, Mirzaei A, Kim S I, Baek S H, Chun D W, Jin C H, Lee K H 2021 Appl. Surf. Sci. 568 150910 doi: 10.1016/j.apsusc.2021.150910 [10] Brophy R E, Junker B, Fakhri E A, Arnason H Ö, Svavarsson H G, Weimar U, Bârsan N, Manolescu A 2024 Sens. Actuat. B-Chem. 410 135648 doi: 10.1016/j.snb.2024.135648 [11] Rani S, Kumar M, Garg P, Rani S, Kumar M, Garg P, Parmar R, Kumar A, Singh Y, Baloria V, Deshpande U, Singh V N 2022 ACS Appl. Mater. Interfaces 14 15381 doi: 10.1021/acsami.1c24679 [12] Yu W, Sisi L, Haiyan Y, Luo J 2020 Rsc Adv. 10 15328 doi: 10.1039/D0RA01068E [13] Dong Q C, Xiao M, Chu Z Y, Li G C, Zhang Y 2021 Sensors 21 3386 doi: 10.3390/s21103386 [14] Gui Y G, Peng X, Liu K Ding Z Y 2020 Physicas E 119 113959 doi: 10.1016/j.physe.2020.113959 [15] Zhu P C, Tang F, Wang S F, Cao W, Wang Q 2022 Mater. Today Commun. 33 104280 doi: 10.1016/j.mtcomm.2022.104280 [16] Li Q F, Chen W L, Liu W H, Sun M L, Xu M H, Peng H L, Wu H Y, Song S X, Li T H, Tang X H 2022 Appl. Surf. Sci. 586 152689 doi: 10.1016/j.apsusc.2022.152689 [17] Hong H S, Ha N H, Thinh, D D, Nam N H, Huong N T, Hue N T, Hoang T V 2021 Nano Energy 87 106165 doi: 10.1016/j.nanoen.2021.106165 [18] Zhang T, Sun H, Wang F D, Zhang W D, Tang S W, Ma J M, Gong H W, Zhang J P 2017 Appl. Surf. Sci. 425 340 doi: 10.1016/j.apsusc.2017.06.229 [19] Choudhuri I, Patra N, Mahata A, Ahuja R, Pathak B 2015 J. Phys. Chem. C 119 24827 doi: 10.1021/acs.jpcc.5b07359 [20] Shukri M S M, Saimin M N S, Yaakob M K, Yahya M Z A, Taib M F M 2019 Appl. Surf. Sci. 494 817 doi: 10.1016/j.apsusc.2019.07.238 [21] Shamim S U D, Roy D, Alam S, Piya A A, Rahman M S, Hossain M K, Ahmed F 2022 Appl. Surf. Sci. 596 153603 doi: 10.1016/j.apsusc.2022.153603 [22] Zhang X X, Yu L, Gui Y G, Hu W H 2016 Appl. Surf. Sci. 367 259 doi: 10.1016/j.apsusc.2016.01.168 [23] Jia X, An L 2019 Mod. Phys. Lett. B 33 1950044 doi: 10.1142/S0217984919500441 [24] Segall M D, Lindan P J D, Probert M J, Pickard C J, Hasnip P J, Clark S J, Payne M C 2002 J. Phys. Condens. Mat. 14 2717 doi: 10.1088/0953-8984/14/11/301 [25] Yang Z H, Wang Z Y, Su X P 2012 J. Cent. South Univ. 19 1796 doi: 10.1007/s11771-012-1210-3 [26] Basiuk V A, Henao-Holguin L V 2014 J. Comput. Theor. Nanosci. 11 1609 doi: 10.1166/jctn.2014.3539 [27] 朱洪强, 冯庆 2014 物理学报 63 133101 doi: 10.7498/aps.63.133101 Zhu H Q, Feng Q 2014 Acta Phys. Sin. 63 133101 doi: 10.7498/aps.63.133101 [28] Gao X, Zhou Q, Wang J X, Xu L N, Zeng W 2020 Appl. Surf. Sci. 517 146180 doi: 10.1016/j.apsusc.2020.146180 [29] Zheng T, Traian D, Thomas F 2021 Phys. Chem. Chem. Phys. 23 19627 doi: 10.1039/D1CP00682G [30] Ping L K, Mohamed M A, Mondal A K, Taib M F M, Samat M H, Berhanuddin D D, Menon P S, Bahru R 2021 Micromachines 12 348 doi: 10.3390/mi12040348 -