搜索

文章查询

x

留言板

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

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

碱金属修饰的多孔石墨烯的储氢性能

元丽华 巩纪军 王道斌 张材荣 张梅玲 苏俊燕 康龙

碱金属修饰的多孔石墨烯的储氢性能

元丽华, 巩纪军, 王道斌, 张材荣, 张梅玲, 苏俊燕, 康龙
PDF
HTML
导出引用
导出核心图
  • 基于第一性原理深入研究了碱金属原子(Li, Na, K)修饰的多孔石墨烯(PG)体系的储氢性能, 并且通过从头算分子动力学模拟了温度对Li-PG吸附的H2分子稳定性的影响. 研究结果表明, PG结构的碳环中心是碱金属原子最稳定的吸附位置, PG单胞最多可以吸附4个碱金属原子, Li原子被束缚最强, 金属原子间无团聚的倾向; H2分子通过极化机制吸附在碱金属修饰的PG结构上, 每个金属原子周围最多可以稳定地吸附3个H2分子; Li-PG对H2分子的吸附最强(平均吸附能为–0.246 eV/H2), Na-PG对H2分子的吸附较弱(平均吸附能为–0.129 eV/H2), K-PG对H2分子的吸附最弱(平均吸附能为–0.056 eV/H2), 不适合用做储氢材料; 在不考虑外界压强且温度为300 K的情况下, Li-PG结构可稳定地吸附9个H2分子, 储氢量为9.25 wt.%; 在400 K时, 有7个吸附H2分子脱离Li-PG的束缚, 在600—700 K的范围内, 吸附H2分子全部脱离了Li-PG体系的束缚.
      通信作者: 元丽华, yuanlh@lut.edu.cn
    • 基金项目: 省部级-甘肃省自然科学基金(17JR5RA123)
    [1]

    Jiang L L, Fan Z J 2014 Nanoscale 6 1922

    [2]

    Liu W, Wang Z F, Shi Q W, Yang J L, Liu F 2009 Phys. Rev. B 80 233405

    [3]

    Bieri M, Treier M, Cai J M, Aït-Mansour K, Ruffieux P, Gröning O, Gröning P, Kastler M, Rieger R, Feng X L, Müllen K, Fasel R 2009 Chem. Commun. 45 6919

    [4]

    Du A J, Zhu Z H, Smith S C 2010 J. Am. Chem. Soc. 132 2876

    [5]

    Li Y F, Zhou Z, Shen P W, Chen Z F 2010 Chem. Commun. 46 3672

    [6]

    Hatanaka M 2010 Chem. Phys. Lett. 488 187

    [7]

    Jiang D E, Cooper V R, Dai S 2009 Nano Lett. 9 4019

    [8]

    Merchant C A, Healy K, Wanunu M, Ray V, Petermn N, Battel J, Fischbein M D, Venta K, Luo Z T, Johnson C, Drndić M 2010 Nano Lett. 10 2915

    [9]

    Mukherjee R, Thomas A V, Datta D, Singh E, Li J W, Eksik O, Shenoy V B, Koratkar N 2014 Nature Commun. 5 3710

    [10]

    Bi H, Xie X, Yin K, Zhou Y, Wan S, He L, Xu F, Banhart F, Sun L, Ruoff R S 2012 Adv. Funct. Mater 22 4421

    [11]

    Reunchan P, Jhi S H 2011 Appl. Phys. Lett. 98 093103

    [12]

    Ao Z M, Dou S X, Xu Z M, Jiang Q G, Wang G X 2014 Int. J. Hydrogen energy 39 16244

    [13]

    Yuan L H, Chen Y H, Kang L, Zhang C R, Wang D B, Wang C N, Zhang M L, Wu X J 2017 Appl. Surf. Sci. 399 463

    [14]

    Yuan L H, Kang L, Chen Y H, Wang D B, Gong J J, Wang C N, Zhang M L, Wu X J 2018 Appl. Surf. Sci. 434 843

    [15]

    Yuan L H, Wang D B, Gong J J, Zhang C R, Zhang L P, Zhang M L, Wu X J, Kang L 2019 Chem. Phys. Lett. 726 57

    [16]

    Lu R F, Rao D W, Lu Z L, Qian J C, Li F, Wu H P, Wang Y Q, Xiao C Y, Deng K M, Kan E J, Deng W Q 2012 J. Phys. Chem. C 116 21291

    [17]

    Huang S H, Miao L 2012 J. Appl. Phys. 112 124312

    [18]

    Wang Y, Ji Y, Li M, Yuan P, Sun Q, Jia Y 2011 J. App. Phys. 110 094311

    [19]

    Chen Y D, Yu S, Zhao W H, Li S F, Duan X M 2018 Phys. Chem. Chem. Phys. 20 13473

    [20]

    Wang F D, Zhang T, Hou X Y, Zhang W Q, Tang S W, Sun H, Zhang J P 2017 Int. J. Hydrogen Energy 42 10099

    [21]

    Hashmi A, Farooq M U, Khan I, Son J, Hong J 2017 J. Mater. Chem. A 5 2821

    [22]

    Ning G Q, Xu C G, Mu L, Chen G J, Wang G, Gao J S, Fan Z J, Qian W Z, Wei F 2012 Chem. Commun. 48 6815

    [23]

    Xia K S, Tian X L, Fei S X, You K 2014 Int. J. Hydrogen Energy 39 11047

    [24]

    Elyassi M, Rashidi A, Hantehzadeh M R, Elahi S M 2017 Surf. Interface Anal 49 230

    [25]

    Vanderbilt D 1990 Phys. Rev. B 41 7892

    [26]

    Perdew J P, Burke K, Ernzerhof M 1996 Phys. Rev. Lett. 77 3865

    [27]

    Ao Z M, Peeters F M 2010 Phys. Rev. B 81 205406

    [28]

    Rao D W, Lu R F, Meng Z S, Wang Y, Lu Z L, Liu Y Z, Chen X, Kan E, Xiao C Y, Deng K M, Wu H P 2014 Int.J. Hydrogen Energy 39 18966

    [29]

    基泰尔 C 著 (项金钟, 吴兴惠 译) 2009 固体物理导论 (北京: 化学工业出版社) 第38页

    Kittel C (translated by Xiang J Z, Wu X H) 2009 Introduction to Solid State Physics (Beijing: Chemical Industry Press) p38 (in Chinese)

    [30]

    周晓锋, 方浩宇, 唐春梅 2019 物理学报 68 053601

    Zhou X F, Fang H Y, Tang C M 2019 Acta Phys. Sin. 68 053601

    [31]

    Seentithurai S, Kodi Pandya R, Vinodh Kumar S, Saranya C, mahendran M 2014 Int. J. Hydrogen Energy 39 11016

    [32]

    Sacchi M, Jenkins S J 2014 Phys. Chem. Chem. Phys. 16 6101

  • 图 1  一个碱金属原子修饰PG优化后的几何结构 (a) Li-PG; (b) Na-PG; (c) K-PG

    Fig. 1.  Optimized geometry structure of a alkaline metal atom decorated PG: (a) Li-PG; (b) Na-PG; (c) K-PG.

    图 2  Li-PG吸附氢气分子弛豫后的几何结构(红色字体表示H—H键长)

    Fig. 2.  Optimized geometry structure of the Li-PG with H2 molecules adsorption. Red digits represent the corresponding bond length of H—H.

    图 3  Na-PG吸附H2分子优化后的几何结构图(红色数字表示H—H键长)

    Fig. 3.  Optimized geometry structure of the Na-PG with H2 molecules adsorption. Red digits represent the corresponding bond length of H—H

    图 4  K-PG结构吸附H2分子优化后的几何结构图(红色数字表示H—H键长)

    Fig. 4.  Optimized geometry structure of the K-PG with H2 molecules. Red digits represent the corresponding bond length of H—H

    图 5  (a) Li-PG和(b) Na-PG结构吸附一个H2分子的分波态密度图

    Fig. 5.  Partial density of states (PDOS) of a H2 molecule on (a) Li-PG and (b) Na-PG.

    图 6  Li-PG吸附H2分子的差分电荷密度图(电荷密度等值面是0.01 e / Å3) (a) n = 1; (b) n = 4

    Fig. 6.  Charge density differences of n H2 adsorbed on Li-PG system for (a) n = 1 and (b) n = 4. The isovalue is taken to be 0.01 e/Å3

    图 7  Na-PG (a)及K-PG (b)吸附1个H2的差分电荷密度图(电荷密度等值面是0.01 e/Å3)

    Fig. 7.  Charge density differences of a H2 adsorbed on Na-PG (a) and K-PG (b) system. The isovalue is taken to be 0.01 e/Å3.

    图 8  Li-PG结构吸附12个H2分子优化后的几何结构图

    Fig. 8.  Optimized geometry structure of the Li-decorated PG with 12 H2 adsorption.

    图 9  经过AIMD模拟后12 H2-Li-PG平衡结构(a)及总能随模拟时间变化的函数关系(b)

    Fig. 9.  Equilibrated structure of the 12 H2-Li-PG (a) and fluctuations of total energy as a function of simulation time (b) in AIMD simulations at 300 K.

    表 1  H2吸附能(ΔEad)、平均吸附能($\Delta {\bar E_{{\rm{ad}}}}$)及Li所带的电荷(q)

    Table 1.  Adsorption energies (ΔEad) and average adsorption energies ($\Delta {\bar E_{{\rm{ad}}}}$)of H2 molecules, and the charge of Li (q).

    n(H2)
    1234
    ΔEad /eV–0.308–0.239–0.249–0.044
    $\Delta {\bar E_{{\rm{ad}}}}$/eV·H2–1–0.308–0.274–0.266–0.225
    q /e1.231.401.561.56
    下载: 导出CSV

    表 2  Na-PG吸附H2分子的吸附能

    Table 2.  Adsorption energies of H2 molecules on Na-PG.

    n(H2)
    123
    ΔEh1/eV–0.202–0.262–0.215
    ΔEh2/eV–0.234–0.117–0.166
    下载: 导出CSV
  • [1]

    Jiang L L, Fan Z J 2014 Nanoscale 6 1922

    [2]

    Liu W, Wang Z F, Shi Q W, Yang J L, Liu F 2009 Phys. Rev. B 80 233405

    [3]

    Bieri M, Treier M, Cai J M, Aït-Mansour K, Ruffieux P, Gröning O, Gröning P, Kastler M, Rieger R, Feng X L, Müllen K, Fasel R 2009 Chem. Commun. 45 6919

    [4]

    Du A J, Zhu Z H, Smith S C 2010 J. Am. Chem. Soc. 132 2876

    [5]

    Li Y F, Zhou Z, Shen P W, Chen Z F 2010 Chem. Commun. 46 3672

    [6]

    Hatanaka M 2010 Chem. Phys. Lett. 488 187

    [7]

    Jiang D E, Cooper V R, Dai S 2009 Nano Lett. 9 4019

    [8]

    Merchant C A, Healy K, Wanunu M, Ray V, Petermn N, Battel J, Fischbein M D, Venta K, Luo Z T, Johnson C, Drndić M 2010 Nano Lett. 10 2915

    [9]

    Mukherjee R, Thomas A V, Datta D, Singh E, Li J W, Eksik O, Shenoy V B, Koratkar N 2014 Nature Commun. 5 3710

    [10]

    Bi H, Xie X, Yin K, Zhou Y, Wan S, He L, Xu F, Banhart F, Sun L, Ruoff R S 2012 Adv. Funct. Mater 22 4421

    [11]

    Reunchan P, Jhi S H 2011 Appl. Phys. Lett. 98 093103

    [12]

    Ao Z M, Dou S X, Xu Z M, Jiang Q G, Wang G X 2014 Int. J. Hydrogen energy 39 16244

    [13]

    Yuan L H, Chen Y H, Kang L, Zhang C R, Wang D B, Wang C N, Zhang M L, Wu X J 2017 Appl. Surf. Sci. 399 463

    [14]

    Yuan L H, Kang L, Chen Y H, Wang D B, Gong J J, Wang C N, Zhang M L, Wu X J 2018 Appl. Surf. Sci. 434 843

    [15]

    Yuan L H, Wang D B, Gong J J, Zhang C R, Zhang L P, Zhang M L, Wu X J, Kang L 2019 Chem. Phys. Lett. 726 57

    [16]

    Lu R F, Rao D W, Lu Z L, Qian J C, Li F, Wu H P, Wang Y Q, Xiao C Y, Deng K M, Kan E J, Deng W Q 2012 J. Phys. Chem. C 116 21291

    [17]

    Huang S H, Miao L 2012 J. Appl. Phys. 112 124312

    [18]

    Wang Y, Ji Y, Li M, Yuan P, Sun Q, Jia Y 2011 J. App. Phys. 110 094311

    [19]

    Chen Y D, Yu S, Zhao W H, Li S F, Duan X M 2018 Phys. Chem. Chem. Phys. 20 13473

    [20]

    Wang F D, Zhang T, Hou X Y, Zhang W Q, Tang S W, Sun H, Zhang J P 2017 Int. J. Hydrogen Energy 42 10099

    [21]

    Hashmi A, Farooq M U, Khan I, Son J, Hong J 2017 J. Mater. Chem. A 5 2821

    [22]

    Ning G Q, Xu C G, Mu L, Chen G J, Wang G, Gao J S, Fan Z J, Qian W Z, Wei F 2012 Chem. Commun. 48 6815

    [23]

    Xia K S, Tian X L, Fei S X, You K 2014 Int. J. Hydrogen Energy 39 11047

    [24]

    Elyassi M, Rashidi A, Hantehzadeh M R, Elahi S M 2017 Surf. Interface Anal 49 230

    [25]

    Vanderbilt D 1990 Phys. Rev. B 41 7892

    [26]

    Perdew J P, Burke K, Ernzerhof M 1996 Phys. Rev. Lett. 77 3865

    [27]

    Ao Z M, Peeters F M 2010 Phys. Rev. B 81 205406

    [28]

    Rao D W, Lu R F, Meng Z S, Wang Y, Lu Z L, Liu Y Z, Chen X, Kan E, Xiao C Y, Deng K M, Wu H P 2014 Int.J. Hydrogen Energy 39 18966

    [29]

    基泰尔 C 著 (项金钟, 吴兴惠 译) 2009 固体物理导论 (北京: 化学工业出版社) 第38页

    Kittel C (translated by Xiang J Z, Wu X H) 2009 Introduction to Solid State Physics (Beijing: Chemical Industry Press) p38 (in Chinese)

    [30]

    周晓锋, 方浩宇, 唐春梅 2019 物理学报 68 053601

    Zhou X F, Fang H Y, Tang C M 2019 Acta Phys. Sin. 68 053601

    [31]

    Seentithurai S, Kodi Pandya R, Vinodh Kumar S, Saranya C, mahendran M 2014 Int. J. Hydrogen Energy 39 11016

    [32]

    Sacchi M, Jenkins S J 2014 Phys. Chem. Chem. Phys. 16 6101

  • [1] 张忠强, 于凡顺, 刘珍, 张福建, 程广贵. 氢化多孔石墨烯反渗透特性及机理分析. 物理学报, 2020, (): . doi: 10.7498/aps.69.20191761
    [2] 赵银昌, 戴振宏, 隋鹏飞, 张晓玲. 二维Li+BC3结构高储氢容量的研究. 物理学报, 2013, 62(13): 137301. doi: 10.7498/aps.62.137301
    [3] 颜克凤, 李小森, 孙丽华, 陈朝阳, 夏志明. 储氢笼型水合物生成促进机理的分子动力学模拟研究. 物理学报, 2011, 60(12): 128801. doi: 10.7498/aps.60.128801
    [4] 刘源, 姚洁, 陈驰, 缪灵, 江建军. 氢修饰石墨烯纳米带压电性质的第一性原理研究. 物理学报, 2013, 62(6): 063601. doi: 10.7498/aps.62.063601
    [5] 赵玉娜, 高涛, 吕金钟, 马俊刚. Li-N-H储氢体系热力学性质的第一性原理研究. 物理学报, 2013, 62(14): 143101. doi: 10.7498/aps.62.143101
    [6] 盛喆, 戴显英, 苗东铭, 吴淑静, 赵天龙, 郝跃. 各Li吸附组分下硅烯氢存储性能的第一性原理研究. 物理学报, 2018, 67(10): 107103. doi: 10.7498/aps.67.20172720
    [7] 吴江滨, 钱耀, 郭小杰, 崔先慧, 缪灵, 江建军. 硅纳米团簇与石墨烯复合结构储锂性能的第一性原理研究. 物理学报, 2012, 61(7): 073601. doi: 10.7498/aps.61.073601
    [8] 侯滨朋, 淦作亮, 雷雪玲, 钟淑英, 徐波, 欧阳楚英. 第一性原理对氮掺杂石墨烯作为锂-空电池阴极材料还原氧分子的机理研究. 物理学报, 2019, 68(12): 128801. doi: 10.7498/aps.68.20190181
    [9] 卢金炼, 曹觉先. 单个钛原子储氢能力和储氢机制的第一性原理研究. 物理学报, 2012, 61(14): 148801. doi: 10.7498/aps.61.148801
    [10] 常旭. 多层石墨烯的表面起伏的分子动力学模拟. 物理学报, 2014, 63(8): 086102. doi: 10.7498/aps.63.086102
    [11] 周晶晶, 陈云贵, 吴朝玲, 郑欣, 房玉超, 高涛. 新型轻质储氢材料的第一性原理原子尺度设计. 物理学报, 2009, 58(7): 4853-4861. doi: 10.7498/aps.58.4853
    [12] 杨光敏, 梁志聪, 黄海华. 石墨烯吸附Li团簇的第一性原理计算. 物理学报, 2017, 66(5): 057301. doi: 10.7498/aps.66.057301
    [13] 张加宏, 顾芳, 顾斌, 杨丽娟. 应变石墨烯纳米带谐振特性的分子动力学研究. 物理学报, 2011, 60(5): 056103. doi: 10.7498/aps.60.056103
    [14] 惠治鑫, 贺鹏飞, 戴瑛, 吴艾辉. 硅功能化石墨烯热导率的分子动力学模拟. 物理学报, 2014, 63(7): 074401. doi: 10.7498/aps.63.074401
    [15] 张忠强, 李冲, 刘汉伦, 葛道晗, 程广贵, 丁建宁. 石墨烯碳纳米管复合结构渗透特性的分子动力学研究. 物理学报, 2018, 67(5): 056102. doi: 10.7498/aps.67.20172424
    [16] 李杰杰, 鲁斌斌, 线跃辉, 胡国明, 夏热. 纳米多孔银力学性能表征分子动力学模拟. 物理学报, 2018, 67(5): 056101. doi: 10.7498/aps.67.20172193
    [17] 刘峰斌, 陈文彬, 崔岩, 屈敏, 曹雷刚, 杨越. 活性质吸附氢修饰金刚石表面的第一性原理研究. 物理学报, 2016, 65(23): 236802. doi: 10.7498/aps.65.236802
    [18] 姜平国, 汪正兵, 闫永播. 三氧化钨表面氢吸附机理的第一性原理研究. 物理学报, 2017, 66(8): 086801. doi: 10.7498/aps.66.086801
    [19] 付宝勤, 侯氢, 汪俊, 丘明杰, 崔节超. 钨空位捕获氢及其解离过程的分子动力学. 物理学报, 2019, 68(24): 240201. doi: 10.7498/aps.68.20190701
    [20] 黄艳平, 袁健美, 郭刚, 毛宇亮. 硅烯饱和吸附碱金属原子的第一性原理研究. 物理学报, 2015, 64(1): 013101. doi: 10.7498/aps.64.013101
  • 引用本文:
    Citation:
计量
  • 文章访问数:  169
  • PDF下载量:  11
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-05-07
  • 修回日期:  2020-01-15
  • 刊出日期:  2020-03-01

碱金属修饰的多孔石墨烯的储氢性能

  • 1. 兰州理工大学理学院, 兰州 730050
  • 2. 兰州理工大学材料科学与工程学院, 兰州 730050
  • 3. 兰州大学核科学与技术学院, 兰州 730000
  • 通信作者: 元丽华, yuanlh@lut.edu.cn
    基金项目: 省部级-甘肃省自然科学基金(17JR5RA123)

摘要: 基于第一性原理深入研究了碱金属原子(Li, Na, K)修饰的多孔石墨烯(PG)体系的储氢性能, 并且通过从头算分子动力学模拟了温度对Li-PG吸附的H2分子稳定性的影响. 研究结果表明, PG结构的碳环中心是碱金属原子最稳定的吸附位置, PG单胞最多可以吸附4个碱金属原子, Li原子被束缚最强, 金属原子间无团聚的倾向; H2分子通过极化机制吸附在碱金属修饰的PG结构上, 每个金属原子周围最多可以稳定地吸附3个H2分子; Li-PG对H2分子的吸附最强(平均吸附能为–0.246 eV/H2), Na-PG对H2分子的吸附较弱(平均吸附能为–0.129 eV/H2), K-PG对H2分子的吸附最弱(平均吸附能为–0.056 eV/H2), 不适合用做储氢材料; 在不考虑外界压强且温度为300 K的情况下, Li-PG结构可稳定地吸附9个H2分子, 储氢量为9.25 wt.%; 在400 K时, 有7个吸附H2分子脱离Li-PG的束缚, 在600—700 K的范围内, 吸附H2分子全部脱离了Li-PG体系的束缚.

English Abstract

    • 多孔石墨烯(PG)是一种在石墨烯层状结构中存在纳米孔的二维材料, 纳米孔的大小和几何形状可以被调控[1,2], 不同几何构型的PG展示出金属、半金属及半导体的不同性能. 瑞士科学家Bieri等[3]在2009年首次合成了一种PG结构, 从扫描隧道显微镜图可以清晰地看到纳米孔结构均匀地分配在石墨烯平面内, 或者准确地描述为石墨烯平面内周期性的缺失苯环, 纳米孔间距是7.4 Å. 这种PG结构的单胞由两个C6H3环组成[4,5], 可以通过3 × 3石墨烯超胞去掉一些碳原子后再利用H原子饱和悬挂键而得到. 利用密度泛函理论[4,5]与晶体轨道法[6]计算分析PG的能带结构表明其属于宽带隙的半导体材料, 这种材料解决了石墨烯0带隙的缺陷. 由于纳米孔洞的存在, PG展现出与纯石墨烯不同的特性, PG在储氢[4]、气体分离或净化[5,7]、DNA测序[8]、超级电容器[9]、水净化[10]等许多领域具有潜在的应用价值.

      基于Bieri等[3]实验制备的PG结构, 理论研究表明此结构适合用作储氢材料. Du等[4]利用第一性原理采用VASP软件研究了Li修饰PG的储氢性能, Li原子稳定地吸附在碳环中心, 结合能为1.81 eV, 远高于Li吸附在完整石墨烯的结合能(0.86 eV). Li-PG结构吸附12个H2分子的吸附能为0.243 eV/H2, 储氢量高达12 wt.%. Reunchan和Jhi[11]利用VASP软件, 采用广义梯度近似(GGA)泛函研究了金属修饰PG的储氢性能, 文中只给出了PG和金属原子间的结合能、体系吸附H2分子数目及H2分子吸附能等计算结果, 他们指出Ca-PG最适合做储氢材料, 而Al-PG结构吸附H2较弱, 不适合用作储氢材料. 但Ao等[12]利用Dmol3软件采用局域密度近似(LDA)泛函计算表明Al修饰的PG能够提高H2分子的吸附能力, Al稳定吸附于PG结构的CH环中心上方, Al的结合能为1.78 eV. PG单胞两面最多可吸附两个Al原子, 而且每个Al原子周围最多可以吸附6个H2分子, H2的吸附能在0.41—1.11 eV/H2范围内, 储氢量为10.5 wt.%. 本课题组[13-15]近年来研究了较重金属原子修饰PG的储氢性能, PG单胞单侧最多可以吸附1个金属原子, 否则金属原子有团聚的倾向.

      为了提高金属原子和PG间的结合能力, Lu课题组[16]采用B和N掺杂来改性PG, 利用VASP软件选取GGA泛函研究了Li和Ca修饰PG的储氢性能. 研究结果表明N掺杂使得金属和PG之间的作用变弱, B掺杂虽然可以提高金属和PG之间的作用, 但是掺杂B原子的数目不能超过2个. B掺杂可增强Ca-PG吸附H2分子的强度, 可是B掺杂略减弱Li-PG吸附H2分子的强度 (H2分子吸附在Li-B-PG的吸附能为0.08 eV/H2, 略低于吸附在Li-PG上的吸附能0.09 eV/H2). Huang和Miao[17]采用O原子来改性PG (用1个O原子替代PG中的1个H原子), 利用SIESTA软件研究了Li修饰PG体系的储氢性能, Li的结合能高达到3.84 eV, Li原子周围吸附5个H2, H2的吸附能在0.2—0.5 eV范围. Wang等[18]利用N原子改性PG (PG中连接H的C原子和H原子被1个N原子替代, 称之为g-CN), 利用锂和钙原子共同修饰g-CN结构, Li和Ca稳定的吸附位置分别为不完整碳环中心和纳米孔中心, 储氢量为9.17 wt.%. 最近, Chen等[19]研究Li修饰g-CN层的储氢性能发现单个Li原子和CN层的相互作用很强(结合能为4.73 eV), N掺杂可以提高PG的储氢能力. Wang等[20]对一种PG结构(单胞由C—C键连接3个C6H2环组成, 简称CMP)的储氢性能做了研究, 此结构与Bieri等[3]的PG结构(单胞由C—C键连接两个C6H3环组成)很相似, 从头算分子动力学模拟表明在温度为300 K时, 体系储氢量为10.89 wt.%. 此外, 进一步的研究[21]表明N掺杂此结构可增强对Li的吸附作用, 体系储氢量得到提高.

      PG实验储氢的研究较少, Ning等[22]利用化学气相沉积模板法(template CVD)制备出PG, 表面积为2038 m2·g–1, 纳米孔尺寸约为1 nm, 在3.1 MPa和274 K的条件下, 可吸附15.3 mmol/g的H2分子. Xia等[23]研究了不同结构的PG吸附, 具有孔表面积大的PG结构吸附能力强, 但储氢量仅有0.75 wt.% (77 K 和1 bar (1 bar = 105 Pa)), 说明纯的PG由于化学惰性使得其与氢气的结合较弱, 需经过金属原子或其他化学官能团修饰才适合储氢. Elyassi等[24]实验结果表明, 在室温298 K和压强34 bar的条件下, PG、石墨烯及氧化石墨烯的储氢量分别是2.56 wt.%, 1.70 wt.%及0.74 wt.%, 并指出PG的高吸附量是因为其具有高的孔体积、大的比表面积和较小的孔尺寸.

      已有研究结果表明Li修饰PG可以增强PG对H2分子吸附能力, 但温度对Li-PG吸附H2分子稳定性的影响以及释放H2的温度没有被进一步讨论, 对于碱金属原子Na和K修饰PG储氢性能的研究很少. 本文主要利用第一性原理系统深入研究了碱金属原子Li, Na, K修饰PG的储氢性能并且利用从头算分子动力学(ab initio MD, AIMD)模拟Li-PG结构吸附的H2分子数目随温度的变化情况.

    • 本文计算利用CASTEP软件, 考虑自旋极化效应, 选取超软赝势[25], 采用GGA的Perdew-Burke-Ernzerhof (PBE)交换关联泛函[26]. 我们对截断能和K点的选择进行了收敛测试, 截断能选取500 eV, K点为6 × 6 × 1, 与前人对PG相关计算选取的参数一致[4,5]. 几何结构优化无对称性限制, 计算中SCF收敛标准选取1.0 × 10–6 eV/atom, 能量收敛标准选取5.0 × 106 eV, 结构弛豫止于原子间的力小于0.01 eV/Å. 为了计算中考虑范德瓦耳斯作用, 采用了色散校正DFT-D的Grimme方法(我们在文献[13]中曾讨论了此方法对计算吸附能的影响). 计算模型选取的PG单胞由两个C6H3环构成, 真空层为18 Å, 利用周期性边界条件.

      此外, 我们利用分子动力学研究温度对Li-PG吸附H2分子性能影响时选取了CASTEP中的Dynamics, 采用NVT系综及Nośe-Hoover恒温扩展方法, 外界压强不考虑, 步长为1 fs. 因为Dynamics模拟计算非常耗时, Wang[20]和Ao[27]利用分子动力学计算选取了1.5 ps与1 ps的模拟时长, 为了使计算结果有意义且节约时间, 我们选取的模拟时长为1.5 ps.

      碱金属原子和PG间的平均结合能定义为

      $\Delta {\bar E_{\rm{b}}} = \left( {{E_{n{\rm{M + PG}}}} - {E_{{\rm{PG}}}} - n{E_{\rm{M}}}} \right)/n, $

      其中, EnM+PG, EPGEM对应的是n个碱金属原子修饰PG的总能量、PG单胞总能量以及单个碱金属原子总能量; n为碱金属原子个数.

      H2分子的吸附能和平均吸附能分别定义为

      $ \Delta {E_{{\rm{ad}}}} = \left[ {{E_{n{{\rm{H}}_{\rm{2}}} + {\rm{M + PG}}}} - {E_{(n - 1){{\rm{H}}_{\rm{2}}} + {\rm{M + PG}}}} - {E_{{{\rm{H}}_{\rm{2}}}}}} \right], $

      $ \Delta {\bar E_{{\rm{ad}}}} = \left( {{E_{n{{\rm{H}}_{\rm{2}}}{\rm{ + M + PG}}}} - {E_{{\rm{M + PG}}}} - n{E_{{{\rm{H}}_{\rm{2}}}}}} \right)/n, $

      其中, ${E_{n{{\rm{H}}_{\rm{2}}} + {\rm{M + PG}}}}$, ${E_{(n - 1){{\rm{H}}_{\rm{2}}} + {\rm{M + PG}}}}$, ${E_{{\rm{M + PG}}}}$${E_{{{\rm{H}}_{\rm{2}}}}}$分别指n个H2吸附在碱金属原子-PG结构(M-PG)的总能, (n – 1)个H2吸附在M-PG的总能, M-PG结构的总能及单个H2分子的能量.

    • 首先确定单个碱金属原子吸附在PG上的稳定位置, 图S1 (见补充材料(online))描述了金属原子可能吸附的初始位置. 分别放于不同位置的Li, Na, K原子, 结构优化后Li和Na占据h1和h2位置, K占据h1和b2位置, 见表S1 (见补充材料(online)). 通过对碱金属原子的结合能(表S1)做进一步分析, 发现Li原子最稳定的吸附位置是完整C环中心h1, 结合能量值为1.654 eV (Du[4]采用LDA泛函计算结果为1.81 eV, Rao[28]采用PBE泛函计算结果为1.38 eV, Lu等[16]采用PW91泛函计算结果为1.40 eV); Na位于h1和h2位置的两种几何结构对应的基态能量差别很小(结合能–0.877 eV和–0.879 eV, 见表S1); K最稳定的吸附位置是h1 (结合能–1.037 eV, 见表S1). 图1描绘了碱金属原子吸附在PG上优化后的结构, 图中标出的红色数据是金属原子与碳环中心的距离, Li原子最靠近PG层, 表明Li与PG的相互作用最强. 分析Mulliken原子布居, Li带电荷+1.14 e, Na带电荷+0.99 e, K带电荷+1.02 e (Hirshfeld电荷分析, Li带电荷0.42 e, Na带电荷0.35 e, K带电荷0.36 e), Li转移到PG的电荷最多, 进一步表明Li与PG的相互作用最强. 碱金属原子带正电荷, PG层的C原子带负电荷, 因而PG层和金属原子之间形成静电场, 有利于基体吸附H2分子.

      图  1  一个碱金属原子修饰PG优化后的几何结构 (a) Li-PG; (b) Na-PG; (c) K-PG

      Figure 1.  Optimized geometry structure of a alkaline metal atom decorated PG: (a) Li-PG; (b) Na-PG; (c) K-PG.

      PG同侧吸附两个金属原子优化后的几何结构如图S2 (见补充材料(online))所示, Li原子或K原子位于PG层不同碳环中心, Na原子略微偏离碳环中心. Li-PG与K-PG结构未发生明显的形变, Na-PG结构发生了微小形变, Li, Na, K原子的平均结合能分别为–1.670, –1.590, –1.123 eV, 结合能量值均大于内聚能[11,29,30]. 两个吸附Li原子间距及两个吸附K原子间距均为4.329 Å (Du等[4]计算的Li-Li间距为4.23 Å), 由PG的2 × 2超胞测得Na原子与其最近邻Na原子的间距为4.057 Å (大于钠晶体最近邻两个Na原子间距3.716 Å), 因而PG同侧吸附两个碱金属原子而不会发生团聚现象. PG两侧吸附金属原子优化后的结构如图S3 (见补充材料(online))所示, Li, Na, K的平均结合能分别是为–1.670, –1.123, –1.161 eV, Li或K原子修饰的PG结构未发生明显形变, 但Na-PG几何结构略微形变, 4个Na原子略偏离了完整碳环中心.

    • 为了明确Li, Na, K修饰对PG吸附H2分子能力的影响, 研究了未掺杂的PG结构吸附H2分子的情况, 1个H2分子被放置于PG的6个不同的初始吸附位置(如图S1所示)分别进行结构优化, 结果发现放于给定初始位置的H2分子不会自动移动到其他位置, PG吸H2分子后结构未发生形变, 第1个H2分子吸附于h1位置(完整碳环中心)能最低(氢的吸附能为–0.091 eV), H—H键长保持不变(自由H2的H—H键长为0.753 Å, 此结果和文献[11]计算结果一致); 第2个H2分子稳定的吸附在碳氢环中两H原子之间, 如图S4 (见补充材料(online))所示, H2分子平均吸附能为–0.082 eV/H2; 吸附第3个H2分子后平均吸附能为–0.093 eV/H2, 可能是H2分子间的相互作用使得吸附能量值增大, 于是在PG上继续吸附第4个H2分子, 平均吸附能为–0.067 eV/H2. 这些H2分子的平均吸附能的量值低于0.1 eV/H2, 说明PG对H2分子的吸附能力很弱, 因而无法保证能有效地储氢.

      Li-PG结构单侧吸附H2分子的几何优化结构如图2所示, H2的吸附能、平均吸附能及Li所带电荷如表1所示. 通过分析Mulliken电荷布居, 组成第1个H2的两H原子带–0.05 e和–0.11 e电荷, Li带+1.23 e电荷, H原子从Li得到电荷而带负电, 则H2分子的共价键被极化而变弱, 通过静电场作用被吸附在Li原子周围. H2的吸附能随着吸附H2数目增多而逐渐减小, 第4个H2吸附能低于0.1 eV, 表明Li原子周围可较强地吸附3个H2, 与以前报道一致[4]. 第4个H2分子吸附能仅有–0.044 eV, 而4个吸附H2的平均吸附能是–0.225 eV/H2, 如利用平均吸附能判断, Li周围可稳定吸附4个H2. 通过对结果深入分析发现, 第4个H2的键长0.754 Å; 通过Mulliken电荷布居分析发现, 吸附第4个H2后Li所带的电荷保持不变(如表1所列), 第4个H2的1个H原子得到0.01 e的电荷, 另1个H原子没有电荷得失; 通过Hirshfeld电荷分析发现吸附第4个H2前后, Li原子带的电荷保持不变(0.24 e), 而第4个H2的两H原子不带电荷(0 e). 这说明第4个H2属于物理吸附, 不能稳定的吸附在体系上, 因而利用吸附能可以更准确地判断基体最多吸附H2的数目.

      n(H2)
      1234
      ΔEad /eV–0.308–0.239–0.249–0.044
      $\Delta {\bar E_{{\rm{ad}}}}$/eV·H2–1–0.308–0.274–0.266–0.225
      q /e1.231.401.561.56

      表 1  H2吸附能(ΔEad)、平均吸附能($\Delta {\bar E_{{\rm{ad}}}}$)及Li所带的电荷(q)

      Table 1.  Adsorption energies (ΔEad) and average adsorption energies ($\Delta {\bar E_{{\rm{ad}}}}$)of H2 molecules, and the charge of Li (q).

      图  2  Li-PG吸附氢气分子弛豫后的几何结构(红色字体表示H—H键长)

      Figure 2.  Optimized geometry structure of the Li-PG with H2 molecules adsorption. Red digits represent the corresponding bond length of H—H.

      位于碳环中心(h1构型)和碳氢环中心(h2构型)的Na原子吸附能相近, 因此考虑这两种结构吸附H2分子的情况. 图3和图S5 (见补充材料(online))分别描述了h1和h2构型吸附H2弛豫后的几何结构. 在表2中, ΔEh1是h1构型H2分子的吸附能, ΔEh2是h2构型H2分子的吸附能. 通过分析吸附能, 发现随着H2数目增加, h2构型吸附H2的能力明显降低. 对于h1构型, 3个H2可以稳定地吸附在Na周围. 因而, h1构型是吸附H2分子的最佳结构. 分析Mulliken原子布居发现, Na-PG结构吸附H2之前, Na带+0.99 e的电荷, C均带负电荷. 体系吸附1个H2后, Na带+1.32 e的电荷, H2的两H原子带–0.07 e和–0.15 e的电荷, 表明Na阳离子把0.33 e个电荷转移给H2和C原子. 带负电的C和正电的Na之间形成静电场, H2被极化使得H—H键长由0.753 Å变为0.767 Å.

      图  3  Na-PG吸附H2分子优化后的几何结构图(红色数字表示H—H键长)

      Figure 3.  Optimized geometry structure of the Na-PG with H2 molecules adsorption. Red digits represent the corresponding bond length of H—H

      n(H2)
      123
      ΔEh1/eV–0.202–0.262–0.215
      ΔEh2/eV–0.234–0.117–0.166

      表 2  Na-PG吸附H2分子的吸附能

      Table 2.  Adsorption energies of H2 molecules on Na-PG.

      图4是K-PG吸附1到3个H2优化后的几何结构. 体系吸附的第1个H2垂直于PG平面, 键长为0.761 Å, 吸附能仅有–0.120 eV表明H2被较弱吸附. 当体系吸附第2个H2后, K原子返回碳环中心, 第2个H2吸附能为–0.176 eV, 两H2分子的H—H键长为0.763 Å, 均垂直于PG平面. 第3个H2的吸附能仅有–0.110 eV, 3个H2垂直于PG平面且位于H原子之上. 随着H2分子数目的增多, H2分子的吸附能逐渐减小, K周围较弱地吸附3个H2分子.

      图  4  K-PG结构吸附H2分子优化后的几何结构图(红色数字表示H—H键长)

      Figure 4.  Optimized geometry structure of the K-PG with H2 molecules. Red digits represent the corresponding bond length of H—H

    • 图5和图S6 (见补充材料(online))描述了Li-PG, Na-PG, K-PG吸附1个H2的分波态密度图. 氢分子的σ轨道和Li原子2s轨道间没有发生轨道耦合现象(见图5(a)). 从图5(b)看出, H2分子σ轨道出现在费米能级下约–9.1 eV处, Na原子的3s或2p轨道和氢分子的σ轨道均没有发生轨道耦合现象. K原子的s和p轨道与H2分子的σ轨道也无明显的轨道耦合现象(见图S6). 这说明H2分子吸附在碱金属修饰的PG体系上不是通过Kubas相互作用机制, H2分子主要通过极化机制以及弱范德瓦耳斯相互作用吸附在碱金属修饰的PG上.

      图  5  (a) Li-PG和(b) Na-PG结构吸附一个H2分子的分波态密度图

      Figure 5.  Partial density of states (PDOS) of a H2 molecule on (a) Li-PG and (b) Na-PG.

      为了进一步了解氢的吸附机理, 计算分析了H2吸附体系的差分电荷密度. H2-Li-PG的差分电荷密度如图6及图S7 (见补充材料(online))所示, 黄色和蓝色区域代表电荷的减少和聚集. H2分子两侧分别出现了电荷聚集和减少(图6(a)及图S7), 表明H2分子被极化, 带正电的Li和带负电的C之间的静电场使氢分子极化, H2分子通过极化机制被吸附到Li-PG体系上[4,18]. 从图6(b)观察到围绕在Li周围的3个H2两边均有电荷聚集和减少, 但远离Li原子的第4个氢分子两侧无电荷分布, 表明Li原子和第4个H2分子间无电荷转移, 这与Mulliken原子布居(见表1)分析一致, 进一步表明Li原子周围可以稳定吸附3个H2分子. 以上分析与利用氢分子吸附能判断的结论一致, 说明氢分子的吸附能比平均吸附更能准确地判断吸附H2分子的数目. Seenithurai等[31]在研究Li修饰3 × 3石墨烯吸附H2分子时, LDA泛函计算已高估吸附能, 而作者仅通过平均吸附能判断每个Li原子周围可以吸附5个H2分子, 这一结论不准确. 图7描述了Na-PG及K-PG吸附1个H2分子差分电荷密度图, 可以看出H2分子通过极化机制被吸附到体系上, 但Na (或K)与H2分子间聚集的电荷很少, 表明Na (或K)和H2分子的相互作用较弱.

      图  6  Li-PG吸附H2分子的差分电荷密度图(电荷密度等值面是0.01 e / Å3) (a) n = 1; (b) n = 4

      Figure 6.  Charge density differences of n H2 adsorbed on Li-PG system for (a) n = 1 and (b) n = 4. The isovalue is taken to be 0.01 e/Å3

      图  7  Na-PG (a)及K-PG (b)吸附1个H2的差分电荷密度图(电荷密度等值面是0.01 e/Å3)

      Figure 7.  Charge density differences of a H2 adsorbed on Na-PG (a) and K-PG (b) system. The isovalue is taken to be 0.01 e/Å3.

    • PG单胞两侧最多可以吸附4个碱金属原子, 最多可以吸附12个H2分子. Li-PG吸附H2分子优化后的稳定几何结构如图8所示, 在该结构中12个H2分子都倾斜的分布于Li的周围, 平均吸附能为–0.246 eV/H2 (Du等[4]的计算结果为–0.243 eV/H2), 储氢量为12 wt.%. Na-PG吸附H2分子优化后的稳定几何结构如图S9 (见补充材料(online))所示, H2分子的平均吸附能为–0.129 eV/H2, 表明Na-PG体系对H2分子的吸附弱于Li-PG体系. K-PG吸附H2分子优化后的几何结构如图S10 (见补充材料(online))所示, 氢分子的平均吸附能为–0.056 eV/H2. 尽管4个K原子可以稳定地被吸附在PG单胞上, 但H2分子被K-PG结构吸附较弱, 室温下不能稳定被吸附, 因而K-PG不适合用作储氢材料.

      图  8  Li-PG结构吸附12个H2分子优化后的几何结构图

      Figure 8.  Optimized geometry structure of the Li-decorated PG with 12 H2 adsorption.

      以上计算忽略了零点能(ZPE), 为了考察ZPE对能量的贡献, 利用CASTEP软件通过计算声子得到ZPE (针对最适合储氢的Li-PG结构进行了计算). 采用不同的方法计算了PG的ZPE, 其中线性响应(linear response)法计算的ZPE为3.674 eV, 有限位移(finite displacement)法计算的ZPE为3.695 eV, 这两种方法计算的ZPE值相差很小, 由于之前的计算采用超软赝势, 所以采用有限位移法. Li-PG结构的ZPE为3.730 eV, 而单个Li原子的ZPE为0 eV, 所以ZPE对Li的结合能的贡献仅有0.035 eV (即: 如考虑ZPE校正, 1个Li原子与PG间的结合能为–1.619 eV). 利用有限位移法计算的H2分子的ZPE为0.289 eV (Sacchi和Jenkins [32]计算得到的H2分子ZPE为0.273 eV), Li-PG吸附1个H2分子后的ZPE为4.127 eV, ZPE对第1个H2分子吸附能的贡献为0.108 eV (即: 考虑ZPE校正, 吸附在Li-PG上的第1个H2分子的吸附能为–0.2 eV). 计算结果说明计算H2分子的吸附能需要考虑ZPE校正.

      为了解最适合储氢的Li-PG结构在室温下吸附H2分子的能力, 通过AIMD模拟12 H2-Li-PG在不考虑压强且温度300 K下的稳定性. 图9(a)描述了经过AIMD模拟后12 H2-Li-PG平衡结构, Li-PG结构未发生明显形变, 其中的3个H2 (黄色)脱离体系束缚, 这3个H2和最近邻Li的距离均超过4 Å. 因此, 在室温下, Li修饰的PG单胞最多可吸附9个H2, 储氢量为9.25 wt.%. 总能量随模拟时间的函数关系如图9(b)所示, 图中的曲线在很小能量范围内发生振动, 没有发生大幅度涨落, 进一步表明在室温下Li-PG结构是稳定的. 对12 H2-Na-PG做了AIMD模拟, 在温度为300 K时, 有5个H2分子脱体系的束缚(图S11 (见补充材料(online))), 储氢量为5.47 wt.%.

      图  9  经过AIMD模拟后12 H2-Li-PG平衡结构(a)及总能随模拟时间变化的函数关系(b)

      Figure 9.  Equilibrated structure of the 12 H2-Li-PG (a) and fluctuations of total energy as a function of simulation time (b) in AIMD simulations at 300 K.

      此外, 我们通过AIMD模拟研究了不同温度对吸附H2分子稳定性的影响. 在温度较低时(100 K和150 K), H2分子仍较强的吸附在Li原子周围, 随着温度的升高, 1个H2分子开始远离Li原子, 在250 K时, 有两个H2分子脱离Li原子的束缚, 如图S12 (见补充材料(online))所示. 温度继续升高, 更多的H2分子脱离体系的束缚, 在400 K时, 有7个H2分子脱离体系的束缚, 即58%的H2分子被释放; 随着温度的升高, 释放的H2分子增多, 在600—700 K范围内所有的H2分子全部被释放.

    • 本文深入研究了Li, Na, K分别修饰PG单胞的储氢性能. PG单胞最多可吸附4个碱金属原子, 借助H2分子吸附能、分波态密度、Mulliken原子布居及差分电荷密度计算分析, 我们发现最多有3个H2分子可稳定吸附在碱金属原子周围, H2分子被吸附是利用极化机制, 氢的吸附能比平均吸附能可以更能准确地判断体系吸附的H2分子数目. Li对H2分子的吸附最强, K对H2分子的吸附最弱, 因而K-PG不适合储氢. Li-PG体系在理想条件下的储氢量高达12 wt.%, 为了进一步明确温度对Li-PG体系吸附H2分子的性能的影响, 通过AIDM模拟计算了12 H2-Li-PG结构随温度变化的情况. 在温度为300 K且不考虑外界压强, 最多有9个H2分子被吸附在Li修饰的PG单胞上, 储氢量为9.25 wt.%. 随着温度的升高, H2分子相继脱离了Li原子的束缚, 在400 K时, 有58%的H2分子被释放, 在600—700 K范围内所有的H2分子全部被释放.

参考文献 (32)
补充材料:

目录

    /

    返回文章
    返回