搜索

文章查询

x

留言板

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

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

三明治结构graphene-2Li-graphene的储氢性能

周晓锋 方浩宇 唐春梅

三明治结构graphene-2Li-graphene的储氢性能

周晓锋, 方浩宇, 唐春梅
PDF
HTML
导出引用
导出核心图
  • 本文使用密度泛函理论中的广义梯度近似对扩展三明治结构graphene-2Li-graphene的几何结构、电子性质和储氢性能进行计算研究. 计算得知: 位于单层石墨烯中六元环面心位上方的单个Li原子与基底之间的结合能最大(1.19 eV), 但小于固体Li的实验内聚能(1.63 eV), 然而, 在双层石墨烯之间的单个Li原子与基底的结合能增加到3.41 eV, 远大于固体Li的实验内聚能, 因此位于双层石墨烯之间的多个Li原子不会成簇, 有利于进一步储氢. 扩展三明治结构graphene-2Li-graphene中每个Li原子最多可以吸附3个H2分子, 储氢密度高达10.20 wt.%, 超过美国能源部制定的5.5 wt.%的目标. 该体系对1—3 个H2分子的平均吸附能分别为0.37, 0.17和0.12 eV, 介于物理吸附和化学吸附(0.1—0.8 eV)之间, 因此该体系可以实现常温常压下对H2的可逆吸附. 通过对态密度分析可知, 每个Li原子主要通过电场极化作用吸附多个H2分子. 动力学和巨配分函数计算表明graphene-2Li-graphene结构对H2分子具有良好的可逆吸附性能. 该研究可以为开发良好的储氢材料提供一个好的研究思路, 为实验工作提供理论依据.
      通信作者: 周晓锋, xfzhouphy@263.net ; 唐春梅, tcmnj@163.com
    • 基金项目: 中央高校基本科研业务费(批准号:2016B01914,2018B19414)、水利科技创新项目(批准号:2015087)、江苏省自然科学基金(批准号:BK20161501)和江苏省六大人才高峰项目(批准号:2015-XCL-010)资助的课题.
    [1]

    Schlapbach L, Zuttel A 2001 Nature 414 353

    [2]

    Chandrakumar K R S, Ghosh S K 2008 Nano Lett. 8 13

    [3]

    Rosi N L, Eckert J, Eddaoudi M, Vodak D T, Kim J, O'keeffe M, Yaghi O M 2003 Science 300 1127

    [4]

    Han S S, Goddard W A 2007 J. Am. Chem. Soc. 129 8422

    [5]

    Kealy T J, Pauson P L 1951 Nature 168 1039

    [6]

    Sun Q, Wang Q, Jena P, Kawazoe Y 2005 J. Am. Chem. Soc. 127 14582

    [7]

    Kim D, Lee S, Hwang Y, Yun K H, Chung Y C 2014 Int. J. Hydrogen Energy 39 13189

    [8]

    Xu B, Lei X L, Liu G, Wu M S, Ouyang C Y 2014 Int. J. Hydrogen Energy 39 17104

    [9]

    Seenithurai S, Pandyan R K, Kumar S V, Saranya C, Mahendran M 2014 Int. J. Hydrogen Energy 39 11016

    [10]

    Chen L, Zhang Y, Koratkar N, Jena P, Nayak S K 2008 Phys. Rev. B 77 033405

    [11]

    Mauron P, Gaboardi M, Remhof A, Bliersbach A, Sheptyakov D, Aramini M, Vlahopoulou G, Giglio F, Pontiroli D, Ricco M, Zuttel A 2013 J. Phys. Chem. C 117 22598

    [12]

    Lein M, Frunzke J, Frenking G 2003 Inorg. Chem. 42 2504

    [13]

    Youn I S, Kim D Y, Singh N J, Park S W, Youn J, Kim K S 2011 J. Chem. Theor. Comput. 8 99

    [14]

    Kealy T J, Pauson P L 1951 Nature 168 1039

    [15]

    Wilkinson G, Rosenblum M, Whiting M C, Woodward R B 1952 J. Am. Chem. Soc. 74 2125

    [16]

    Kubas G J 2001 Kluwer Academic (New York: Plenum Publishing)

    [17]

    Lein M, Frunzke J, Frenking G 2003 Inorg. Chem. 42 2504

    [18]

    Youn I S, Kim D Y, Singh N J, Park S W, Youn J, Kim K S 2011 J. Chem. Theory Comput. 8 99

    [19]

    Sun Q, Wang Q, Jena P, Kawazoe Y 2005 J. Am. Chem. Soc. 127 14582

    [20]

    Delley B 1990 J. Chem. Phys. 92 508

    [21]

    Zhang Q Y, Tang C M, Zhu W H, Cheng C 2018 J. Phys. Chem. C 122 22838

    [22]

    Chang L T, Wei C, Xiao H T 2006 Chin. Phys. 15 2718

    [23]

    Zhao J Y, Zhao F Q, Xu S Y, Ju X H 2013 J. Phys. Chem. A 117 2213

    [24]

    Grimme S, Antony J, Ehrlich S, Krieg H 2010 J. Chem. Phys. 132 154104

    [25]

    Ma L, Zhang J M, Xu K W 2014 Appl. Surf. Sci. 292 921

    [26]

    Gao Y, Wu X, Zeng X C 2014 J. Mater. Chem. A 2 5910

    [27]

    Park J, Burova S, Rodgers A S, Lin M C 1999 J. Phys. Chem. A 103 9036

    [28]

    Abad E, Dappe Y J, Martínez J I, Flores F, Ortega J 2011 J. Chem. Phys. 134 044701

    [29]

    Pliva J, Johns J W C, Goodman L 1991 J. Mol. Spectrosc. 148 427

    [30]

    Toyoda K, Nakano Y, Hamada I, Lee K, Yanagisawa S, Morikawa Y 2009 Surf. Sci. 603 2912

    [31]

    Wang X B, Tang C M, Zhu W H, Zhou X F, Zhou Q H, Cheng C 2018 J. Phys. Chem. C 122 9654

    [32]

    Kealy T J, Pauson P L 1951 Nature 168 1039

    [33]

    Wu G, Wang J, Zhang X, Zhu L 2009 J. Phys. Chem. C 113 7052

    [34]

    汪志诚 2013 热力学·统计物理 (第五版) (北京: 高等教育出版社)

    Wang Z C 2013 Thermodynamics·Statistical Physics (5th Ed.) (Beijing: Higher Education Press) (in Chinese)

  • 图 1  (a)苯环中3种不等价位置; (b) Li原子位于苯环面心位上方的优化结构; (c) C6H6-Li-C6H6势能面扫描曲线和最稳定的C6H6-Li-C6H6三明治结构

    Fig. 1.  (a) Three unequal positions in benzene ring; (b) the optimized structure of the benzene ring with the Li atom located above the face site of the hexagonal ring; (c) potential energy surface scanning curve of C6H6-Li-C6H6 and the most stable sandwich structure of C6H6-Li-C6H6.

    图 2  (a) 单层石墨烯的2 × 3晶胞中Li原子的6个位置; 2个Li原子分别位于(b) ①④组合; (c) ①⑤组合; (d) ①③组合; (e) ①⑥组合; 3个Li原子分别位于(f) ①④⑤组合; (g) ①②⑥组合; (h) 2个Li原子掺杂的最稳定双层石墨烯结构graphene-2Li-graphene

    Fig. 2.  (a) 6 positions of the Li atom in the 2 × 3 unit cell of monolayer graphene; two Li atoms are located at (b) ①④ combination; (c) ①⑤ combination, (d) ①③ combination; (e) ①⑥ composition; three Li atoms are located in (f) ①④⑤ combination; (g) ①②⑥ combination; (h) the most stable graphene-2Li-graphene double-layer graphene structure doped by two Li atoms.

    图 3  graphene-2Li-graphene 的2 × 3晶胞中每个Li原子分别吸附1—4个H2分子的结构图

    Fig. 3.  The structural of the 2 × 3 unit cell of graphene-2Li-graphene with each Li atom adsorbed by 1−4 H2 molecules.

    图 4  不同结构中Li原子或H2分子的态密度图 (a)单独的Li原子; (b) graphene-2Li-graphene中Li原子; (c) graphene-2(Li-H2)-graphene中Li原子; (d) graphene-2(Li-3H2)-graphene中Li原子; (e) 单独的H2分子; (f) graphene-2Li-graphene中graphene; (g) graphene-2(Li-H2)-graphene中H2分子; (h) graphene-2(Li-3H2)-graphene中H2分子

    Fig. 4.  The PDOS of Li atom or H2 molecules: (a) Isolated Li atom; (b) the Li atom in graphene-2Li-graphene; (c) the Li atom in graphene-2(Li-H2)-graphene; (d) the Li atom in graphene-2(Li-3H2)-graphene; (e) isolated H2 molecules; (f) the graphene in graphene-2Li-graphene; (g) the H2 molecules in graphene-2(Li-H2)-graphene; (h) the H2 molecules in graphene-2(Li-3H2)-graphene.

    图 5  graphene-2(Li-3H2)-graphene在77和300 K下5 ps之后的动力学结构图

    Fig. 5.  The structures of graphene-2(Li-3H2)-graphene at 77 and 300 K after 5 ps dynamic times.

    表 1  扩展三明治结构graphene-2(Li-nH2)-graphene[G-2(Li-nH2)-G)](n = 1—4)中的H2分子的Ead, Er, Li和H的平均bader电荷(QLiQH), 双层石墨烯的层间距(DG-G)

    Table 1.  The Ead and Er of H2 molecules average bader charge of Li and H (QLi and QH), interlayer distance of double-layer graphene (DG-G) in the expanded sandwich structure graphene-2(Li-nH2)-graphene[G-2(Li-nH2)-G)](n = 1—4).

    G-2Li-GG-2(Li-H2)-GG-2(Li-2H2)-GG-2(Li-3H2)-GG-2(Li-4H2)-G
    Ead/eV0.370.170.120.06
    Er/eV0.190.190.10–0.08
    QLi/e0.990.620.310.020.01
    QH/e0.200.180.160.12
    DG-G3.693.844.404.904.93
    下载: 导出CSV
  • [1]

    Schlapbach L, Zuttel A 2001 Nature 414 353

    [2]

    Chandrakumar K R S, Ghosh S K 2008 Nano Lett. 8 13

    [3]

    Rosi N L, Eckert J, Eddaoudi M, Vodak D T, Kim J, O'keeffe M, Yaghi O M 2003 Science 300 1127

    [4]

    Han S S, Goddard W A 2007 J. Am. Chem. Soc. 129 8422

    [5]

    Kealy T J, Pauson P L 1951 Nature 168 1039

    [6]

    Sun Q, Wang Q, Jena P, Kawazoe Y 2005 J. Am. Chem. Soc. 127 14582

    [7]

    Kim D, Lee S, Hwang Y, Yun K H, Chung Y C 2014 Int. J. Hydrogen Energy 39 13189

    [8]

    Xu B, Lei X L, Liu G, Wu M S, Ouyang C Y 2014 Int. J. Hydrogen Energy 39 17104

    [9]

    Seenithurai S, Pandyan R K, Kumar S V, Saranya C, Mahendran M 2014 Int. J. Hydrogen Energy 39 11016

    [10]

    Chen L, Zhang Y, Koratkar N, Jena P, Nayak S K 2008 Phys. Rev. B 77 033405

    [11]

    Mauron P, Gaboardi M, Remhof A, Bliersbach A, Sheptyakov D, Aramini M, Vlahopoulou G, Giglio F, Pontiroli D, Ricco M, Zuttel A 2013 J. Phys. Chem. C 117 22598

    [12]

    Lein M, Frunzke J, Frenking G 2003 Inorg. Chem. 42 2504

    [13]

    Youn I S, Kim D Y, Singh N J, Park S W, Youn J, Kim K S 2011 J. Chem. Theor. Comput. 8 99

    [14]

    Kealy T J, Pauson P L 1951 Nature 168 1039

    [15]

    Wilkinson G, Rosenblum M, Whiting M C, Woodward R B 1952 J. Am. Chem. Soc. 74 2125

    [16]

    Kubas G J 2001 Kluwer Academic (New York: Plenum Publishing)

    [17]

    Lein M, Frunzke J, Frenking G 2003 Inorg. Chem. 42 2504

    [18]

    Youn I S, Kim D Y, Singh N J, Park S W, Youn J, Kim K S 2011 J. Chem. Theory Comput. 8 99

    [19]

    Sun Q, Wang Q, Jena P, Kawazoe Y 2005 J. Am. Chem. Soc. 127 14582

    [20]

    Delley B 1990 J. Chem. Phys. 92 508

    [21]

    Zhang Q Y, Tang C M, Zhu W H, Cheng C 2018 J. Phys. Chem. C 122 22838

    [22]

    Chang L T, Wei C, Xiao H T 2006 Chin. Phys. 15 2718

    [23]

    Zhao J Y, Zhao F Q, Xu S Y, Ju X H 2013 J. Phys. Chem. A 117 2213

    [24]

    Grimme S, Antony J, Ehrlich S, Krieg H 2010 J. Chem. Phys. 132 154104

    [25]

    Ma L, Zhang J M, Xu K W 2014 Appl. Surf. Sci. 292 921

    [26]

    Gao Y, Wu X, Zeng X C 2014 J. Mater. Chem. A 2 5910

    [27]

    Park J, Burova S, Rodgers A S, Lin M C 1999 J. Phys. Chem. A 103 9036

    [28]

    Abad E, Dappe Y J, Martínez J I, Flores F, Ortega J 2011 J. Chem. Phys. 134 044701

    [29]

    Pliva J, Johns J W C, Goodman L 1991 J. Mol. Spectrosc. 148 427

    [30]

    Toyoda K, Nakano Y, Hamada I, Lee K, Yanagisawa S, Morikawa Y 2009 Surf. Sci. 603 2912

    [31]

    Wang X B, Tang C M, Zhu W H, Zhou X F, Zhou Q H, Cheng C 2018 J. Phys. Chem. C 122 9654

    [32]

    Kealy T J, Pauson P L 1951 Nature 168 1039

    [33]

    Wu G, Wang J, Zhang X, Zhu L 2009 J. Phys. Chem. C 113 7052

    [34]

    汪志诚 2013 热力学·统计物理 (第五版) (北京: 高等教育出版社)

    Wang Z C 2013 Thermodynamics·Statistical Physics (5th Ed.) (Beijing: Higher Education Press) (in Chinese)

  • [1] 张轶杰, 唐春梅, 高凤志, 王成杰. Li修饰的C6分子对H2O的吸附研究. 物理学报, 2014, 63(14): 147401. doi: 10.7498/aps.63.147401
    [2] 祁鹏堂, 陈宏善. Li修饰的C24团簇的储氢性能. 物理学报, 2015, 64(23): 238102. doi: 10.7498/aps.64.238102
    [3] 高虹, 朱卫华, 唐春梅, 耿芳芳, 姚长达, 徐云玲, 邓开明. 内掺氮富勒烯N2@C60的几何结构和电子性质的密度泛函计算研究. 物理学报, 2010, 59(3): 1707-1711. doi: 10.7498/aps.59.1707
    [4] 罗强, 杨恒, 郭平, 赵建飞. N型甲烷水合物结构和电子性质的密度泛函理论计算. 物理学报, 2019, 68(16): 169101. doi: 10.7498/aps.68.20182230
    [5] 王红艳, 李喜波, 罗江山, 吴卫东, 唐永建. 密度泛函理论研究ScnO(n=1—9)团簇的结构、稳定性与电子性质. 物理学报, 2009, 58(9): 6134-6140. doi: 10.7498/aps.58.6134
    [6] 郭云东, 王红艳, 李喜波, 罗江山, 吴卫东, 唐永建. 密度泛函理论研究Scn,Yn和Lan(n=2—10)团簇的稳定性、电子性质和磁性. 物理学报, 2008, 57(8): 4857-4865. doi: 10.7498/aps.57.4857
    [7] 崔树稳, 李璐, 魏连甲, 钱萍. 双层石墨烯层间限域CO氧化反应的密度泛函研究. 物理学报, 2019, 68(21): 218101. doi: 10.7498/aps.68.20190447
    [8] 王建军, 王飞, 原鹏飞, 孙强, 贾瑜. 石墨烯层间纳米摩擦性质的第一性原理研究. 物理学报, 2012, 61(10): 106801. doi: 10.7498/aps.61.106801
    [9] 元丽华, 巩纪军, 王道斌, 张材荣, 张梅玲, 苏俊燕, 康龙. 碱金属修饰的多孔石墨烯的储氢性能. 物理学报, 2020, 69(6): 068802. doi: 10.7498/aps.69.20190694
    [10] 张秀荣, 李扬, 杨星. WnNim(n+m=8)团簇结构与电子性质的理论研究. 物理学报, 2011, 60(10): 103601. doi: 10.7498/aps.60.103601
    [11] 尹跃洪, 徐红萍. 电场诱导(MgO)4储氢的理论研究. 物理学报, 2019, 68(16): 163601. doi: 10.7498/aps.68.20190544
    [12] 栾晓玮, 孙建平, 王凡嵩, 韦慧兰, 胡艺凡. 锑烯吸附金属Li原子的密度泛函研究. 物理学报, 2019, 68(2): 026802. doi: 10.7498/aps.68.20181648
    [13] 阮文, 谢安东, 余晓光, 伍冬兰. NaBn(n=19)团簇的几何结构和电子性质. 物理学报, 2012, 61(4): 043102. doi: 10.7498/aps.61.043102
    [14] 孙卫国, 刘秀英, 王朝阳, 唐永建, 吴卫东, 张厚琼, 刘淼, 袁磊, 徐嘉靖. 单壁BN纳米管和碳纳米管物理吸附储氢性能的理论对比研究. 物理学报, 2009, 58(2): 1126-1131. doi: 10.7498/aps.58.1126
    [15] 刘立仁, 雷雪玲, 陈杭, 祝恒江. Bn(n=2—15)团簇的几何结构和电子性质. 物理学报, 2009, 58(8): 5355-5361. doi: 10.7498/aps.58.5355
    [16] 赵银昌, 戴振宏, 隋鹏飞, 张晓玲. 二维Li+BC3结构高储氢容量的研究. 物理学报, 2013, 62(13): 137301. doi: 10.7498/aps.62.137301
    [17] 康 龙, 罗永春, 陈玉红, 张材荣, 蒲忠胜. (Li3N)n(n=1—5)团簇结构与性质的密度泛函研究. 物理学报, 2008, 57(7): 4174-4181. doi: 10.7498/aps.57.4174
    [18] 杨光敏, 徐强, 李冰, 张汉壮, 贺小光. 不同N掺杂构型石墨烯的量子电容研究. 物理学报, 2015, 64(12): 127301. doi: 10.7498/aps.64.127301
    [19] 崔焱, 夏蔡娟, 苏耀恒, 张博群, 陈爱民, 杨爱云, 张婷婷, 刘洋. 基于石墨烯电极的齐聚苯乙炔分子器件的整流特性. 物理学报, 2018, 67(11): 118501. doi: 10.7498/aps.67.20180088
    [20] 宋庆功, 赵俊普, 顾威风, 甄丹丹, 郭艳蕊, 李泽朋. 基于密度泛函理论的La掺杂-TiAl体系结构延性与电子性质. 物理学报, 2017, 66(6): 066103. doi: 10.7498/aps.66.066103
  • 引用本文:
    Citation:
计量
  • 文章访问数:  206
  • PDF下载量:  4
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-08-07
  • 修回日期:  2018-12-28
  • 上网日期:  2019-03-23
  • 刊出日期:  2019-03-01

三明治结构graphene-2Li-graphene的储氢性能

    基金项目: 中央高校基本科研业务费(批准号:2016B01914,2018B19414)、水利科技创新项目(批准号:2015087)、江苏省自然科学基金(批准号:BK20161501)和江苏省六大人才高峰项目(批准号:2015-XCL-010)资助的课题.

摘要: 本文使用密度泛函理论中的广义梯度近似对扩展三明治结构graphene-2Li-graphene的几何结构、电子性质和储氢性能进行计算研究. 计算得知: 位于单层石墨烯中六元环面心位上方的单个Li原子与基底之间的结合能最大(1.19 eV), 但小于固体Li的实验内聚能(1.63 eV), 然而, 在双层石墨烯之间的单个Li原子与基底的结合能增加到3.41 eV, 远大于固体Li的实验内聚能, 因此位于双层石墨烯之间的多个Li原子不会成簇, 有利于进一步储氢. 扩展三明治结构graphene-2Li-graphene中每个Li原子最多可以吸附3个H2分子, 储氢密度高达10.20 wt.%, 超过美国能源部制定的5.5 wt.%的目标. 该体系对1—3 个H2分子的平均吸附能分别为0.37, 0.17和0.12 eV, 介于物理吸附和化学吸附(0.1—0.8 eV)之间, 因此该体系可以实现常温常压下对H2的可逆吸附. 通过对态密度分析可知, 每个Li原子主要通过电场极化作用吸附多个H2分子. 动力学和巨配分函数计算表明graphene-2Li-graphene结构对H2分子具有良好的可逆吸附性能. 该研究可以为开发良好的储氢材料提供一个好的研究思路, 为实验工作提供理论依据.

English Abstract

    • 在当前能源危机越演越烈、环境污染越来越严重的情况下, 高效清洁可再生能源氢能由于具有无污染和可循环使用等优点引起了广泛关注. 氢能在实际使用中面临的一个关键问题就是氢存储. 发展高性能储氢材料是解决氢存储的主要途径. 然而, 寻找在室温条件下稳定且有较高储氢密度的材料是一个巨大的挑战. 美国能源部提出2017年实现的储氢材料的储氢密度目标是5.5 wt.%[1,2]. 一般来说, 理想的储氢材料对H2分子的平均吸附能(average adsorption energy, Ead)应介于化学吸附和物理吸附之间(0.1—0.8 eV)[1,2].

      目前, 碳纳米结构由于具有较大的比表面积被认为是较有前途的储氢材料. 然而, 纯碳纳米材料由于与H2分子之间存在弱相互作用并不能有效储氢. 因此, 研究者提出了多种途径来改善它们的储氢性能. 其中, 采用金属原子对碳纳米材料进行修饰可以有效增强材料对氢的吸附强度[3-5], 但金属原子在纳米结构表面容易形成团簇, 致使不能有效储氢. 因此, 金属原子与基底之间的结合能(binding energy, Eb)应大于固体金属的实验内聚能[6]. 据悉, 碱金属Li是体积最小、质量最轻的金属, 采用Li原子对纳米材料进行掺杂后可以在更大空间内吸附较多的H2分子进而增加储氢密度. 据报道, Li原子与含有Stone-Wals缺陷的石墨烯之间的结合能是2.16 eV[7], 高于固体Li的实验内聚能(1.63 eV). 另外, Li掺杂的石墨炔[8]、石墨烯[9]和碳纳米管[10]等都表现出良好的储氢性能. 实验上也已经证明使用Li原子掺杂的纳米材料可以有效储氢, 例如Li12C60[11]和Li掺杂的多孔材料[12,13]等都是良好的储氢材料.

      C6H6-metal-C6H6和C5H5-metal-C5H5等三明治结构由于具有较好的稳定性已经引起了人们的特别关注, 例如: 三明治结构二茂铁(C5H5-Fe-C5H5)[14,15]已被实验合成, 结构中每个C5H5环有5个${\text{π}}$电子, Fe原子的电子结构为[Ar]3d64s2, 最外层有8个未配对电子, 因此, 二茂铁结构共有18个电子, 满足18电子规则[16], 形成类似惰性气体的稳定电子组态, 因此能被实验合成. 此外, Lein等[17]对三明治结构[E5-Ti-E5]2–(其中E = CH, N, P, As, Sb)进行了理论研究, 得知满足18电子规则的过渡金属掺杂三明治结构最稳定. 但是, 如何让其他种类金属掺杂的三明治结构也稳定下来呢?本文提出C6H6-Li-C6H6结构[18]可以通过吸附H2分子进一步稳定自身的电子结构, 成为良好的储氢体系.

      众所周知, 石墨烯是以六元环为基本单位扩展形成的二维平面结构, 目前研究者已经对金属修饰的石墨烯的储氢性能开展了广泛研究, 认为多个金属原子在石墨烯表面上吸附时容易发生聚合现象, 使得体系的储氢性能急剧下降[19]. 为了进一步提高金属原子在石墨烯表面上的结合能, 研究者提出了多种方案, 比如在石墨烯中除去一个或者多个C原子形成空位, 使得金属原子在空位处的结合强度大大增强[20]. 如何进一步增强Li原子与石墨烯之间的结合强度是一个需要迫切解决的问题, 对储氢至关重要. 本文提出由C6H6-Li-C6H6扩展形成的graphene-2Li-graphene三明治结构, 并对它的几何结构、电子性质和储氢性能进行计算研究, 第2部分给出计算方法, 结果和讨论在第3部分给出, 最后在第4部分给出结论.

    • 本文中的计算主要采用Dmol3软件包完成[21]. 计算过程使用密度泛函理论(density functional theory, DFT)[22]中的广义梯度近似方法 (generalized gradient approximation, GGA). 选取PBE (Perdew-Burke-Ernzerhof)交换关联梯度修正和极化函数扩展的双数值原子轨道DNP基组. DNP基组就是函数中包含高于自由原子中的最高占据轨道的一级角动量, 与G6-31(d, p)基组的精度一致[23]. 本文对C, H, Li采用全电子的核处理方式. 众所周知, GGA容易低估体系的结合能, 而局域密度近似 (local density approximation, LDA)又会高估体系的结合能, 因此, 计算中采用了范德瓦耳斯 (van der Waals, vdW)修正(标准的DFT计算基础上进行势能修正, 即DFT-D方法). 本文主要采用交换关联势中的Grimme扩散修正 (PBE-Grimme[24])方法. PBE-Grimme已经被用于计算多种材料的储氢性能, 比如Pd掺杂的石墨烯[25,26]、C48B12和C60Ca32[27]等. 几何结构优化以能量、梯度和位移是否收敛为依据, 自洽场过程以体系的能量和电荷密度是否收敛为依据. 文中对所有的稳定结构都进行了频率计算, 确定无虚频.

      为了验证本文计算方法的正确性, 首先对苯环C6H6进行优化, 得知C—H键长为1.09 Å, C—C键长为1.40 Å, 最高占据分子轨道(highest occupied molecular orbital, HOMO)和最低未占分子轨道(lowest unoccupied molecular orbital, LUMO)之间的能隙(energy gap, Eg)为5.20 eV, 与Park等[27]计算得到的C—H键长(1.08 Å)、Abad等[28]计算得到的C—C键长(1.40 Å)、C—C实验键长(1.39 Å[29])和Toyoda等[30]计算得到的能隙值(4.70 eV)非常接近, 以上结果说明本文使用的计算方法可靠.

    • 首先对单个Li原子在C6H6上方的最稳定位置进行研究. 分析可知, 苯环上方有3种不等价位置, 分别为: 1) C原子顶位上方; 2) C—C键桥位上方; 3) 六元环面心位上方, 如图1(a)所示. 计算得知, 三种同分异构体的优化结构中, Li原子都稳定于面心位, 与苯环中心距离为1.78 Å, 如图1(b)所示. Li原子与基底的结合能binding energy (Eb)定义如下[14,31]:

      图  1  (a)苯环中3种不等价位置; (b) Li原子位于苯环面心位上方的优化结构; (c) C6H6-Li-C6H6势能面扫描曲线和最稳定的C6H6-Li-C6H6三明治结构

      Figure 1.  (a) Three unequal positions in benzene ring; (b) the optimized structure of the benzene ring with the Li atom located above the face site of the hexagonal ring; (c) potential energy surface scanning curve of C6H6-Li-C6H6 and the most stable sandwich structure of C6H6-Li-C6H6.

      ${E_{\rm{b}}} = {E_{{{\rm{C}}_6}{{\rm{H}}_6}}} + {E_{{\rm{Li}}}} - {E_{{{\rm{C}}_6}{{\rm{H}}_6} {\text{-}} {\rm{Li}}}},$

      其中$E_{{\rm C_6}{\rm H_6}} $, ELi以及$ E_{{\rm C_6}{\rm H_6}{\text{-}}{\rm Li}}$分别代表C6H6、单个Li原子以及单个Li原子吸附的C6H6的总能量. 结合能反映金属原子与体系的结合强度, 结合能越大表示与体系的结合强度越大, 体系的热力学稳定性越好. 计算可知: 在苯环面心位置处的Li原子与苯环之间的结合能为1.19 eV, 小于固体Li的实验内聚能(1.63 eV)[7], 可以推测, 当六元环扩展成石墨烯平面时, 多个Li原子在石墨烯表面进行吸附时容易成簇.

      Chen等[10]研究发现三明治结构C6H6-Sc-C6H6比C6H6-Sc结构更稳定, 因此本文接下来进一步将C6H6-Li扩展成三明治结构C6H6-Li-C6H6, 并研究其稳定性. 为了得到Li原子在双层石墨烯之间掺杂时的最佳层间距, 首先对C6H6-Li-C6H6结构的总能量进行扫描, 扫描曲线如图1(c)所示, 以Li原子到苯环平面的距离为变量, 扫描步长为0.1 Å. 由图可知, Li与下层苯环相距1.30 Å时, 体系总能量最低. 计算得知, C6H6-Li-C6H6中Li原子与基底的结合能为3.43 eV, 明显大于固体Li的实验内聚能(1.63 eV)[7], 是C6H6-Li结构中Li与C6H6之间结合能的3倍多, 因此可以进一步推测, 当上下两个C6H6分别扩展成两层石墨烯平面时, 多个Li原子在双层石墨烯之间进行掺杂时不会发生成簇.

      图2(a)中, 单层石墨烯的2 × 3晶胞中一共有6个六元环, 6个面心位分别定义为 ①—⑥, 本文接下来先研究2个Li原子位于平面上的最稳定结构. 根据2个Li原子之间的相对位置, 一共有4种不同组合来放置Li原子, 分别为①②组合、①④组合、①⑤组合和①⑥组合, 它们的优化结构如图2(b)图2(e)所示. 在这四种结构中, Li原子与石墨烯之间的平均结合能(average binding energy, Eavb)定义为石墨烯和Li原子各自的能量之和减去Li掺杂体系的总能量得到的差再除以Li原子数目的值. 计算得知, Li原子与基底的平均结合能分别为1.38, 1.37, 1.36和1.35 eV, 均小于固体Li的实验内聚能(1.63 eV)[7], 表明仍不能有效避免多个Li原子的成簇.

      图  2  (a) 单层石墨烯的2 × 3晶胞中Li原子的6个位置; 2个Li原子分别位于(b) ①④组合; (c) ①⑤组合; (d) ①③组合; (e) ①⑥组合; 3个Li原子分别位于(f) ①④⑤组合; (g) ①②⑥组合; (h) 2个Li原子掺杂的最稳定双层石墨烯结构graphene-2Li-graphene

      Figure 2.  (a) 6 positions of the Li atom in the 2 × 3 unit cell of monolayer graphene; two Li atoms are located at (b) ①④ combination; (c) ①⑤ combination, (d) ①③ combination; (e) ①⑥ composition; three Li atoms are located in (f) ①④⑤ combination; (g) ①②⑥ combination; (h) the most stable graphene-2Li-graphene double-layer graphene structure doped by two Li atoms.

      同时, 当2个Li原子位于石墨烯的同一侧时, 石墨烯会向Li原子所在的一侧发生平面弯曲, 这在周期性结构中是一种较不稳定的现象. 为了改善弯曲的石墨烯平面, 在石墨烯另一侧添加第3个Li原子, 分别位于④和②位置, 如图2(f)图2(g)所示, 图2(f)图2(g)结构中Li与基底的平均结合能分别为1.42 和1.44 eV, 都小于固体Li的实验内聚能(1.63 eV)[7], 表明多个Li原子吸附在单层石墨烯表面上时会发生成簇.

      接下来在双层石墨烯之间放置2个Li原子. 根据2个Li原子之间的相对位置, 仍然考虑4种不同组合的吸附结构, 分别为①②组合、①④组合、①⑤组合和①⑥组合. 在优化结构中Li与基底之间的结合能分别为3.30, 3.35, 3.41和3.43 eV, 均大于固体Li的实验内聚能(1.63 V)[7], 所以当两个Li原子位于双层石墨烯之间时不会成簇. 其中, 两个Li原子为①⑥组合时与基底的平均结合能最大, 因此该结构最稳定, 如图2(h)所示. 我们注意到, 优化结构中Li原子稍稍向外偏移, 说明它们之间有排斥作用, 因此, 两个Li原子的距离相对较大(①⑥组合)时最稳定, 和平均结合能计算结果一致. 同时计算发现, 图2(h)结构(graphene-2Li-graphene)中Li与基底的结合能和C6H6-Li-C6H6中苯环与Li的结合能相同, 说明当苯环扩展成石墨烯平面, 双层石墨之间如①⑥组合放置Li原子时, 上下两层石墨烯对Li原子的结合能基本不会受到影响. 但是, 在图2(h)结构(graphene-2Li-graphene)中两层石墨烯之间的距离增加到3.69 Å (如表1所列), 比C6H6-Li-C6H6结构中两层苯环的距离(2.60 Å)大出了0.89 Å, 说明当苯环扩展成石墨烯平面时, 层间放置的Li原子会排斥双层石墨烯在垂直于平面的方向上向外移动. graphene-2Li-graphene结构中Li原子之间的距离为5.71 Å, 远大于Li—Li二聚体的键长(2.61 Å), 结合双层石墨烯对Li原子的结合能远大于相邻Li原子之间结合能的事实, 证明多个Li原子位于双层石墨烯之间时不会发生成簇, 有利于进一步储氢.

      G-2Li-GG-2(Li-H2)-GG-2(Li-2H2)-GG-2(Li-3H2)-GG-2(Li-4H2)-G
      Ead/eV0.370.170.120.06
      Er/eV0.190.190.10–0.08
      QLi/e0.990.620.310.020.01
      QH/e0.200.180.160.12
      DG-G3.693.844.404.904.93

      表 1  扩展三明治结构graphene-2(Li-nH2)-graphene[G-2(Li-nH2)-G)](n = 1—4)中的H2分子的Ead, Er, Li和H的平均bader电荷(QLiQH), 双层石墨烯的层间距(DG-G)

      Table 1.  The Ead and Er of H2 molecules average bader charge of Li and H (QLi and QH), interlayer distance of double-layer graphene (DG-G) in the expanded sandwich structure graphene-2(Li-nH2)-graphene[G-2(Li-nH2)-G)](n = 1—4).

      因为储氢材料必须具有一定的动力学稳定性, 我们进一步对图2(h)结构(graphene-2Li-graphene)进行了300 K温度下的动力学模拟. 本文使用Nose-Hoover Chain thermostats进行有限的温度调节, 研究体系的动力学稳定性. 在300 K, 1 fs步长, 5 ps后的动力学结构与图2(h)结构基本上没有改变, 因此, graphene-2(Li-3H2)-graphene结构具有较高的热力学和动力学稳定性.

      对于已实验合成的二茂铁(C5H5-Fe-C5H5)[14,15]结构, Fe原子的外层电子排布为[Ar]3d64s2, 每个C5H5单元具有5个${\text{π}}$电子, Fe原子的最外层3d轨道有8个电子, 因此, 该体系一共具有18个电子, 形成类似惰性气体的稳定电子态, 因此以二茂铁结构为代表的过渡金属掺杂的三明治结构的稳定性可以使用18电子规则[32]进行解释. 那么, 如C6H6-Li-C6H6等碱金属掺杂的三明治结构如何能稳定下来呢?我们认为在该结构中吸附H2分子是一个很好的途径, 从该思路出发, 碱金属掺杂的三明治结构显然可以成为良好的储氢体系. 接下来就对扩展三明治结构graphene-2Li-graphene的(2 × 3)单元的储氢性能进行深入研究. 图3(a)图3(d)分别给出了graphene-2Li-graphene体系的(2 × 3)单元中每个Li原子吸附1—4个H2分子后的优化结构. 表1给出了体系对H2的平均吸附能Ead、连续吸附能(consecutive adsorption energy, Ec)、和Li原子的Bader电荷(QLiQH)、双层石墨烯的层间距DG-G.

      图  3  graphene-2Li-graphene 的2 × 3晶胞中每个Li原子分别吸附1—4个H2分子的结构图

      Figure 3.  The structural of the 2 × 3 unit cell of graphene-2Li-graphene with each Li atom adsorbed by 1−4 H2 molecules.

      体系对H2分子的Ead定义为[33]

      $\begin{split} {E_{{\rm{ad}}}} =& \left( {{E_{{\rm{G}}{\text{-}}2{\rm{Li}}{\text{-}} {\rm{G}}}} + n{E_{{{\rm{H}}_2}}} }\right.\\ &-\left.{ {E_{{\rm{G}}{\text{-}}2\left( {{\rm{Li}}{\text{-}}n{{\rm{H}}_2}} \right) {\text{-}} {\rm{G}}}}} \right)/(2n),\end{split}$

      式中EG-2Li-G, ${E_{{{\rm{H}}_2}}}$, $E_{{\rm{G}} {\text{-}} 2\left( {{\rm{Li}}{\text{-}} n{{\rm{H}}_2}} \right){\text{-}} {\rm{G}}} $分别是graphene-2Li-graphene、单个H2分子、和graphene-(Li-nH2)-graphene的总能量. 从表1可知, graphene-2Li-graphene中每个Li原子对1—3个H2分子的平均吸附能位于0.10—0.80 eV之间, 在能够实现可逆吸附能量值范围内[1,2], 因此, graphene-2Li-graphene中的每个Li原子最多吸附3个H2分子, 储氢密度可达10.20 wt.%, 超过美国能源部对储氢材料制定的2017年储氢密度5.5 wt.%的目标[3].

      体系对H2分子的连续吸附能(Er)定义为[33]

      $\begin{split}{E_{\rm r}} =& \left( {{E_{{\rm{G}}{\text{-}}2\left[ {{\rm{Li}}{\text{-}} \left( {n - 1} \right){{\rm{H}}_2}} \right] {\text{-}} {\rm{G}}}} + 2{E_{{{\rm{H}}_2}}} }\right.\\ &- \left.{{E_{{\rm{G}}{\text{-}}\left( {{\rm{Li}}{\text{-}}n{{\rm{H}}_2}} \right) {\text{-}} {\rm{G}}}}} \right)/2,\end{split}$

      式中$E_{{\rm{G}} {\text{-}} 2\left[ {{\rm{Li}}{\text{-}} {n-1}{{\rm{H}}_2}} \right]{\text{-}} {\rm{G}}} $, $E_{{\rm H}_2} $$E_{{\rm{G}} {\text{-}} \left( {{\rm{Li}}{\text{-}} n{{\rm{H}}_2}} \right){\text{-}} {\rm{G}}} $分别是graphene-2[Li-(n–1)H2]-graphene、单个H2分子、和graphene-2(Li-nH2)-graphene的总能量. 体系对H2的连续吸附能越低, 就越难吸附H2分子, 很小的正值或者负值表示体系不能继续吸附H2分子. 从表1可知, 扩展三明治结构graphene-2Li-graphene中的每个Li原子对第4个H2分子的连续吸附能为–0.08 eV, 所以, 该体系中每个Li原子最多可以吸附3个H2分子. 计算可知, H2分子中的H—H键长在0.75—0.76 Å之间, 和自由H2分子的键长(0.75 Å)比较接近, 说明氢是以H2分子形式被吸附. 从图3可知, 当每个Li原子吸附3个H2分子时, 距离第3个H2分子为2.30 Å, 稍大于与另外两个H2分子之间的距离(2.03和1.95 Å). 当每个Li原子吸附4个H2分子时, 与4个H2分子之间的距离分别为2.00, 2.00, 2.03和2.52 Å. 以上结果证明当每个Li原子吸附第3个和第4个H2分子时, 排斥作用越来越明显. 从表1得知, 当Li原子没有吸附H2分子时带有0.99e的电荷, 因此, 最外层s轨道上的1个电子都转移给了C原子, Li与C之间形成离子键, 能在双层石墨烯之间较强结合, 有利于进一步储氢. 但是随着吸附H2分子数目的增加, Li原子所带的电荷逐渐变少, 说明最外层的空轨道逐渐被来自H原子的电子占据. 进一步计算发现, 当双层石墨烯之间的每个Li原子吸附1—3个H2分子时, 双层石墨烯之间的距离会从3.69 Å进一步拉大到3.84, 4.40和4.90 Å, 因此, 在双层石墨烯之间占据的原子越多, 对石墨烯产生的垂直方向上的排斥力越大. 但是当每个Li原子吸附4个H2分子时的距离和吸附3个H2分子时的距离几乎相等, 因为第4个H2分子跑到了很远的位置, 几乎不与Li成键, 因此对层间距影响不大.

      表1给出的Bader电荷可知Li原子带正电荷, 说明Li原子与石墨烯结合时, 会将电荷转移给石墨烯而处于缺电子状态, 所以可以有效吸附H2分子, 而Li原子吸附1—4个H2时每个H原子分别带有平均为0.20e, 0.18e, 0.16e和0.12e的正电荷, 所以, 电荷应该从H2转移到Li原子最外层的空轨道上, 使得Li原子的最外层轨道逐步被填满从而达到饱和状态稳定下来.

      为了进一步理解体系的储氢机理, 本文计算了体系的分态密度(partial density of states, PDOS), 结果如图4所示. 态密度图是通过对体系的离散能级进行洛伦兹展开, 展开系数为0.15 eV, 由水平轨道给出的权重求和而得到. (HOMO+LUMO)/2被定义为费米能级(Fermi energy, EF), 在图中取在0 eV处, 用黑色虚线表示. 由图4可知, 孤立Li原子的PDOS主要分布在–2—6 eV之间, 在费米能级处有两个较宽的峰. Li原子的d轨道态密度分布位于较高能级处, 因此在–15—5 eV的范围内没有被占据. 位于双层石墨烯层间的Li原子的PDOS的分布扩展到–15—–2 eV的范围内, 这表明 Li原子的p轨道与石墨烯的C原子的P轨道之间发生轨道杂化, 如图4(b)图4(f)所示. 每个Li原子吸附1个H2分子后, G-2(Li-H2)-G中H2的PDOS的一个峰分布在–13—5 eV之间, 对应于H2${\text{σ}}$成键态, 进一步说明H—H键未断裂, 如图4(g)所示. H2分子的${\text{σ}}$成键轨道与Li原子的s和p轨道的PDOS在–13— –2 eV之间有明显重叠; 在0—5 eV范围内, H2分子的PDOS有几个低峰, 对应于H2分子的${{\text{σ}}^*}$反键轨道, 它与Li原子的s和p轨道之间也存在杂化, 如图4(c)图4(g)所示. 结合Bader电荷可知, Li原子最外层的1个电子转移给了C原子, 因此带有0.99e的电荷在其周围激发电场, 通过电场极化作用吸附H2分子. H2分子的H—H键因此变长, 吸附作用得到加强. G-2Li-G中Li原子的p轨道主要分布在–15—5 eV之间, 较大的峰位于–10— –5 eV之间, 而吸附H2分子之后移到了费米能级之上的0—5 eV之间, Li原子所带的电荷越来越少, 说明Li原子靠近费米能级处的空轨道逐渐被来自H原子的电子占据, 形成较大的电子占据态, 如图4(d)图4(h)所示. 综上所述, 三明治结构graphene-2Li-graphene中每个Li原子主要通过电场极化作用吸附1—3个H2分子.

      图  4  不同结构中Li原子或H2分子的态密度图 (a)单独的Li原子; (b) graphene-2Li-graphene中Li原子; (c) graphene-2(Li-H2)-graphene中Li原子; (d) graphene-2(Li-3H2)-graphene中Li原子; (e) 单独的H2分子; (f) graphene-2Li-graphene中graphene; (g) graphene-2(Li-H2)-graphene中H2分子; (h) graphene-2(Li-3H2)-graphene中H2分子

      Figure 4.  The PDOS of Li atom or H2 molecules: (a) Isolated Li atom; (b) the Li atom in graphene-2Li-graphene; (c) the Li atom in graphene-2(Li-H2)-graphene; (d) the Li atom in graphene-2(Li-3H2)-graphene; (e) isolated H2 molecules; (f) the graphene in graphene-2Li-graphene; (g) the H2 molecules in graphene-2(Li-H2)-graphene; (h) the H2 molecules in graphene-2(Li-3H2)-graphene.

      为了研究储氢体系是否具有良好的氢释放性能, 接下来对graphene-2(Li-3H2)-graphene结构在温度77 和300 K下进行动力学模拟. 本文使用Nose-Hoover Chain thermostats进行有限的温度调节, 研究储氢结构的氢释放性能. 图5给出了graphene-2(Li-3H2)-graphene结构在77 和300 K下5 ps之后的动力学结构图. 由图5可知, 结构在77 K和1 fs步长下5 ps后释放2个H2分子. 而在300 K和1 fs步长下5 ps后释放5个H2分子. 因此, 在较高的温度下, 更多的H2会被释放出来. 上述结果表明graphene-2(Li-3H2)-graphene结构对H2具有良好的可逆吸附性能.

      图  5  graphene-2(Li-3H2)-graphene在77和300 K下5 ps之后的动力学结构图

      Figure 5.  The structures of graphene-2(Li-3H2)-graphene at 77 and 300 K after 5 ps dynamic times.

      进一步根据巨配分函数[34]求得在一定温度下体系能够吸附的平均H2分子数, 在温度T时, graphene-2Li-graphene能够吸附的平均氢分子数由下面的公式计算得到:

      $\overline N = \displaystyle\frac{\displaystyle{\mathop \sum \nolimits_{N = 1}^{{\rm{max}}} N{{\rm{e}}^{N\left( {{\text{μ}} + {E_{{\rm{a}}N}}} \right)/\left( {kT} \right)}}}}{{1 + \displaystyle\mathop \sum \nolimits_{N = 1}^{{\rm{max}}} {{\rm{e}}^{N\left( {{\text{μ}} + {E_{{\rm{a}}N}}} \right)/\left( {kT} \right)}}}},$

      式中k为玻尔兹曼常数, T是开尔文温度, EaN为体系吸附N个H2时的平均吸附能, $\mu $为一定压强和温度下H2的化学势. 由(4)式计算可知, 在1 atm (1 atm = 1.01 × 105 Pa), 400 K的条件下, 每个Li原子能够吸附的平均H2分子数$\overline N $为0.08, 接近于0, 因此, H2几乎完全释放. 因此, Li原子掺杂的三明治结构graphene-2Li-graphene可以实现在常温常压下对H2分子的可逆吸附.

    • 本文主要得到以下几点结论:

      1)单层石墨烯中六元环面心位上方是单个Li原子的最稳定位置, 但是Li原子与石墨烯的结合能(1.19 eV)小于固体Li的实验内聚能(1.63 eV), 因此当多个Li原子吸附时会发生成簇现象;

      2) Li原子在双层石墨烯之间掺杂形成的扩展三明治结构graphene-2Li-graphene中Li原子与双层石墨烯的结合能高达3.43 eV, 远大于固体Li的实验内聚能(1.63 eV), 因此当多个Li原子在双层石墨烯之间掺杂时不会成簇;

      3)周期三明治结构graphene-2Li-graphene中每个Li原子最多可以吸附3个H2分子, 储氢密度高达10.20 wt.%, 超过美国能源部制定的5.5 wt.%的目标;

      4)通过对态密度分析可知, Li原子掺杂的三明治graphene-2Li-graphene主要通过电场极化的作用吸附H2分子;

      5)动力学和巨配分函数计算表明, Li原子内掺杂的双层石墨烯结构对H2分子具有良好的可逆吸附性能.

参考文献 (34)

目录

    /

    返回文章
    返回