氧化石墨烯的结构稳定性及硝酸催化作用的第一性原理研究*

林启民 张霞† 芦启超 罗彦彬 崔建功 颜鑫 任晓敏 黄雪

1) (北京邮电大学, 信息光子学与光通信国家重点实验室, 北京 100876)

2) (中北大学, 电子测试技术重点实验室, 太原 030051)

3) (北京市计算中心, 北京 100094)

利用第一性原理方法, 采用超软赝势库系统研究了硝酸熏蒸石墨烯得到的氧化石墨烯结构的稳定性及电子结构.基于石墨烯正交元胞的2 × 2超胞模型建立相应的正交晶系硝酸熏蒸氧化石墨烯模型, 包含15个碳原子和2个氧原子.结果表明熏蒸后包含碳氧双键的氧化石墨烯结构为能量较低的稳定结构, 与实验报道一致.力学稳定性分析表明该结构的处于力学稳定状态.通过分析熏蒸前后的反应物和生成物, 表明硝酸起催化作用; 且硝酸氧化石墨烯为吸热过程, 反应发生需要外界热源.通过分析结构的电子特性, 得出氧化石墨烯为直接带隙本征半导体, 带隙值为1.12 eV, 功函数为5.28 eV.研究结果为硝酸氧化石墨烯的制备及其在光电子器件领域的应用提供了理论依据.

1 引 言

单层石墨烯因其良好的电学、光学、磁学和力学性能受到了广泛的关注[1-9].自2004年Novoselov等[1]成功剥离石墨烯片, 到2010年Bae等[10]在滚筒式的铜基底上生长30英寸石墨烯薄膜并对其掺杂, 以及近期在减少使用酸性物质实现对其改性和其氧化生物还原物在生物降解、环境保护方面的突出表现[11,12], 科研工作者一直致力于实现单层石墨烯制备、增大面积、改性等研究[13,14].石墨烯的制备工艺趋于多元化, 其相关物性也在趋于多元化.Bae等[10]报道了在实验室成功合成大面积的石墨烯, 并就经过硝酸熏蒸后得到的氧化石墨烯结构及其光电性能进行了系统的研究.研究显示, 硝酸处理后的石墨烯有良好的光电性能, 但在检测结果中未发现氮元素和氢元素, 表明氧化石墨烯中碳氧原子成键方式为单键和双键.Yan等[15]也对氧化石墨烯结构进行了系统的研究, 他们主要通过理论计算分析了石墨烯氧化后氧原子可能出现的位置, 但结果中并未包含实验中发现的碳氧双键[10,16].因此, 目前对于氧化石墨烯的形成过程及其结构和性质的研究还不完善, 有必要进行进一步研究.

本文通过第一性原理方法研究了单层石墨烯经过硝酸熏蒸后得到的氧化石墨烯结构.首先, 构建了包含氧吸附和氧替位在内的多种可能存在的氧化石墨烯结构模型, 通过分析键长、能量以及弹性常数找出最稳定结构; 然后, 通过对比硝酸熏蒸石墨烯前后结构的形成能, 分析了硝酸在石墨烯氧化过程中的作用; 最后, 计算并分析了该结构的能带结构、电荷密度和功函数, 讨论了其在光电器件中可能的应用前景.

2 计算方法与模型

本工作基于第一性原理密度泛函理论, 计算采用VASP商用软件包实现[17].计算方法采用密度泛函理论中的局域密度泛函, 势函数采用超软赝势, 平面截断能设为 500 eV, 布里渊区所用K空间采用中心校准的Monkhorst-Pack方法,K点采用 5 × 5 × 1 网格划分, 在计算过程中总能在每个原子上收敛标准为10—6eV, 为防止层间不同周期原子间相互作用, 在z方向设置真空层厚度大于15 Å[15].

图1 替位与吸附位组成的石墨烯结构 (a)石墨烯结构; (b)吸附位结构; (c), (d)在同一模型中一个氧原子替换和一个氧原子吸附近邻碳原子结构; (e), (f)在同一结构中一个氧原子替换和一个氧原子吸附不相邻碳原子结构Fig.1.Graphene structures with substitutional and adsorbed oxygen：(a) Pure graphene; (b) adsorbed oxygen; (c), (d) substitutional and adsorbed oxygen between adjacent carbon; (e), (f) substitutional and adsorbed oxygen between nonadjacent carbon.

3 数值结果与分析

3.1 结构及其稳定性

图1(b)为计算得到的石墨烯吸附氧原子和羟基后形成氧化石墨烯的稳定构型.计算结果表明,与同侧吸附相比, 氧原子和羟基的异侧吸附形成包含环氧键的氧化石墨烯结构是最稳定的构型.环氧键的氧原子吸附于石墨烯桥位, 而且氧原子下方的碳原子向上移动, 与相邻碳原子不在同一平面, 导致碳环发生部分扭曲, 形成sp3杂化键, 这一结果与文献[15]中的报道一致.显然, 这一结构中包含的碳氧键均属于碳氧单键, 并未出现实验中测得的碳氧双键[10].

图1(c)—(f)是同时包含氧替位和氧吸附时形成的氧化石墨烯的四种可能结构.计算得到的各结构总能和碳原子与替位以及吸附位氧原子间的键长如表1所列.

表1 结构总能和碳氧原子键长Table 1.Total energy of structures and the bond length of carbon and oxygen.

从表1可以看出, 结构c和d的能量接近且较低, 相较于结构 e 和 f更稳定.其中, 结构 c 的吸附位碳氧键最短 (1.24 Å), 键能较高, 作用较强, 为碳氧双键; 结构d中的碳氧键长为1.38 Å, 吸附的氧原子位于石墨烯碳原子上方.考虑到与吸附氧原子相邻的碳原子在石墨烯中均与周围三个碳原子成键, 故此碳氧键应为单键.而结构e和f中的替位氧原子与周围碳原子处于同一平面, 形成C—O—C 环键, 即形成 sp2杂化键[15], 其地位类似于石墨烯中的碳原子, 未形成碳氧双键.

为更直观地分析成键情况, 做出结构c的电荷密度图如图2所示.从图2可以看出, 替位氧原子并未与吸附氧原子下方的碳原子有电荷交叠, 即并未成键.但替位氧原子与上方(图1中结构c左侧)相邻的两碳原子存在电荷交叠, 分别形成碳-氧单键.同时, 吸附位的氧原子与下方碳原子存在电荷交叠, 形成双键.

图2 电荷密度图 (a)结构c及其电荷密度的俯视图;(b)结构c绕y轴转角50°的电荷密度图(黄色是电荷分布,蓝色是其周期边界剖面)Fig.2.Charge density of structure c：(a) Top view; (b) rotation 50 degrees around the y axis (yellow represents the charge distribution and blue represents the periodic boundary profile).

为进一步分析结构c和d的稳定性, 还计算了它们的力学稳定性.作为稳定的二维正交结构应满足以下力学判据, 即弹性常数关系满足[18]

式中,C11,C22,C12,C66是相应自由度上的弹性常数.

计算得到的结构a, c和d的4个弹性常数C11,C22,C12,C66如表2 所列.

从表2可以看出, 纯净石墨烯结构a的结合性最强.对于氧化石墨烯来说, 结构 d 的C66＜ 0, 不满足二维结构的力学稳定条件; 结构c的C66＞ 0,C11＞ 0, 且C11C22＞C122, 满足力学稳定条件.同时可看出, 该结构C11,C22的值较大, 显现出较强的结合性, 在相应的坐标轴方向不易发生形变.

表2 结构 a, c, d 的弹性常数Table 2.The elastic constants of structures a, c and d.

综合以上分析可以看出, 结构c为氧化石墨烯的最稳定结构, 结构中包含碳氧双键.

3.2 形成能

可通过计算反应前后的形成能变化分析反应过程的难易程度[19].形成能的计算公式如下：

式中Ef为体系的形成能,Ee为末态体系各组分的总能,Eb为初始体系各组分的总能.

通过3.1节中对氧化石墨烯的结构稳定性分析, 可知图1(c)中的结构最为稳定, 结合硝酸处理石墨烯的实验过程, 可写出如下化学方程式：

分别分析反应前后的各组分的能量, 结合(3)和(4)式, 可计算得到上述反应过程的形成能为 7.58 eV, 即这一反应为吸热反应, 故反应的发生需要有外界热源.

考虑到空气中氧气与水蒸气共存, 可以与生成物中的二氧化氮发生反应形成硝酸, 化学反应方程式如下：

可计算得到硝酸的形成能为—7.05 eV, 为放热反应且形成能较低, 故此过程极易发生.即生成物中的二氧化氮和水可与空气中的氧气发生反应产生硝酸.

因此, 通过对硝酸熏蒸石墨烯的整体过程进行分析后, 可得到如下总化学反应方程式：

可计算得到氧化石墨烯的形成能为—6.52 eV.此反应过程能量降低趋于稳定, 容易发生.

由(6)式可知, 硝酸在此过程中起到催化作用.而硝酸参与的第一步反应有生成物能量增加, 即内能增加为吸热反应, 再考虑到硝酸挥发要吸热, 因此对反应过程予以辅热, 不仅可以提高硝酸挥发,而且可以促进反应进行从而提高反应速率.

3.3 能带结构与功函数

包含碳氧双键的氧化石墨烯能带结构图如图3所示.从图3可以看出, 与零带隙的石墨烯相比,氧化石墨烯变为直接带隙的本征半导体材料, 其导带底与价带顶均出现在倒易空间Γ点和Χ点之间的同一点, 带隙值为 1.12 eV.其能带结构显示价带顶下移和导带底上移导致带隙增大.

功函数反映表面电子逃逸的难易程度, 可用来分析材料的电荷转移机制[20].功函数计算公式为[7,21]

其中Evac和EF分别表示真空能量和费米能.

图4(a)和图4(b)分别为石墨烯和氧化石墨烯沿z轴方向的势能分布图.从图4可以看出, 石墨烯氧化前后的功函数分别4.43 eV和5.28 eV.其中,计算得到的氧化前石墨烯功函数的值与文献[7]中通过广义梯度近似求得的功函数值4.69 eV相近.可见, 氧化后的石墨烯功函数增大 0.85 eV, 这是由于氧元素进入石墨烯结构后使得轨道杂化类型改变导致价带顶下移和导带底上移, 使电子逸出功增大, 即功函数增加, 这与能带结构图中显示的氧化石墨烯带隙增加的结论相符.实验研究表明,硝酸熏蒸 5 min 后, 其功函数净增加 130 meV[10].X射线光电子能谱结果显示纯石墨烯结构中存在碳氧单键与双键, 这是由于石墨烯在陈放或制备后转移过程中已经被氧化, 因此测试所得功函数为石墨烯氧化过程中一段过程, 值相对偏小.其次理论计算虽然本身具有局限性值会偏小, 但是由于计算硬件限制建立模型较小相应的点缺陷浓度较高, 计算所得功函数相对较高.综上所述, 本文的计算结果与实验报道定性符合.氧化石墨烯的直接带隙性质和功函数增大的特性, 有利于其在光电器件中的应用.

图3 包含碳氧双键的氧化石墨烯能带结构图(图中红色虚线代表费米能级)Fig.3.Band structures of the graphene oxide with carbon oxygen double bond (red dashed line represents the Fermi level).

图4 石墨烯氧化前后沿z轴方向的势能分布图 (a)石墨烯势能分布图; (b)氧化石墨烯势能分布图Fig.4.Potential energy along the z-axis：(a) Graphene; (b) graphene oxide.

4 结 论

基于第一性原理方法研究了利用硝酸处理单层石墨烯后得到的包含碳氧双键的氧化石墨烯结构.结果表明：硝酸在处理石墨烯发生氧化的过程中起到催化作用; 硝酸催化作用的第一步反应为吸热反应, 因此在实验过程中进行加热可有效提升反应速率; 且整体反应过程的形成能较低, 反应可以稳定进行.反应得到的包含碳氧双键的氧化石墨烯结构模型满足能量稳定性条件和力学稳定性条件,是稳定结构.氧化石墨烯结构中碳氧单键与碳氧双键同时出现, 且碳氧单键以环氧的形式(C—O—C)存在于石墨烯中, 与实验报道相符.氧化石墨烯结构的带隙打开, 为直接带隙, 带隙值为 1.12 eV; 氧化石墨烯的功函数相对于石墨烯的功函数增大0.85 eV, 与实验报道定性符合.以上研究可以为理解氧化石墨烯的结构及其在光电器件中的进一步应用奠定基础.