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

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

元丽华1)†    巩纪军1)    王道斌1)    张材荣1)

张梅玲1)3)    苏俊燕1)    康龙2)

1) (兰州理工大学理学院, 兰州　730050)2) (兰州理工大学材料科学与工程学院, 兰州　730050)3) (兰州大学核科学与技术学院, 兰州　730000)(2019 年5 月7日收到; 2020 年1 月15日收到修改稿)

基于第一性原理深入研究了碱金属原子(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体系的束缚.

关键词：多孔石墨烯, 储氢, 第一性原理, 分子动力学PACS：88.30.R–, 81.05.ue, 71.15.Mb, 71.15.Pd

　DOI: 10.7498/aps.69.20190694

H原子饱和悬挂键而得到. 利用密度泛函理论[4,5]

1   引　言

多孔石墨烯(PG)是一种在石墨烯层状结构中存在纳米孔的二维材料, 纳米孔的大小和几何形状可以被[1,2], 不同几何构型的PG展示出金属、半金属及半导体的不同性能. 瑞士科学家Bieri等[3]在2009年首次合成了一种PG结构, 从扫描隧道显微镜图可以清晰地看到纳米孔结构均匀地分配在石墨烯平面内, 或者准确地描述为石墨烯平面内周期性的缺失苯环, 纳米孔间距是7.4 Å. 这种PG结构的单胞由两个C6H3环组成[4,5], 可以通过3 × 3石墨烯超胞去掉一些碳原子后再利用 061801)和兰州理工大学红柳一流学科建设项目资助的课题.†  通信作者. E-mail: yuanlh@lut.edu.cn

与晶体轨道法[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分子的吸附

*  甘肃省自然科学基金(批准号: 17JR5RA123)、国家自然科学基金(批准号: 51562022)、兰州理工大学博士基金(批准号:

© 2020 中国物理学会  Chinese Physical Society

068802-1

http://wulixb.iphy.ac.cn

物 理 学 报   Acta  Phys.  Sin.   Vol. 69, No. 6 (2020)    068802

能为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. 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分子数目随温度的变化情况.

2   模型与计算方法

本文计算利用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 Å, 利用周期性边界条件.

068802-2

物 理 学 报   Acta  Phys.  Sin.   Vol. 69, No. 6 (2020)    068802

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

最强. 分析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分子.

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

¯b=(EnM+PG−EPG−nEM)/n,∆E(1) 

(a)

Top view

(b)(c)

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

[]∆Ead=EnH2+M+PG−E(n−1)H2+M+PG−EH2,(2)H2分子的吸附能和平均吸附能分别定义为

Side view1.755 A

2.303 A

2.659 A

¯ad=(EnH+M+PG−EM+PG−nEH)/n,∆E22(3)图 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.

其中,  EnH2+M+PG ,  E(n−1)H2+M+PG ,  EM+PG 及 EH2 分别指n个H2吸附在碱金属原子-PG结构(M-PG)的总能, (n – 1)个H2吸附在M-PG的总能, M-PG结构的总能及单个H2分子的能量.

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原子略偏离了完整碳环中心.

3   结果与讨论

3.1 碱金属原子修饰PG的几何结构和电子结构

首先确定单个碱金属原子吸附在PG上的稳定位置, 图S1 (见补充材料(online))描述了金属原子可能吸附的初始位置. 分别放于不同位置的Li, Na, K原子, 结构优化后Li和Na占据h1和h2位置, K占据h1和b2位置, 见表S1 (见补充材料(online)). 通过对碱金属原子的结合能(表S1)做进一步分析, 发现Li原子最稳定的吸附位置是完整C环中心h1, 结合能量值为用LDA泛函计算结果为泛函计算结果为1.38 eV, Lu

1.6 eV (Du[4]

采

1.81 eV, Rao[28]

等[16]

采用PBE

采用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的相互作用

3.2 碱金属原子修饰PG吸附H2分子的

性能

3.2.1    PG单侧吸附氢分子的结构

为了明确Li, Na, K修饰对PG吸附H2分子

068802-3

物 理 学 报   Acta  Phys.  Sin.   Vol. 69, No. 6 (2020)    068802

能力的影响, 研究了未掺杂的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.7 Å; 通过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的数目.

位于碳环中心(h1构型)和碳氢环中心(h2构型)的Na原子吸附能相近, 因此考虑这两种结构吸附H2分子的情况. 图3和图S5 (见补充材料

(online))分别描述了h1和h2构型吸附H2弛豫后的几何结构. 在表2中, DEh1是h1构型H2分子的吸附能, DEh2是h2构型H2分子的吸附能.通过分析吸附能, 发现随着H2数目增加, h2构型吸附H2的能力明显降低. 对于h1构型, 3个H2可以稳定地吸附在Na周围. 因而, h1构型是吸附H2分子的最佳结构. 分析Mulliken原子布居发现, Na-PG结构吸附H2之前, Na带+0.99 e的电

(a)Top view(b)(c)

0.763 A

0.761 A

0.762 A

0.770 A

Side view

0.766 A

0.765 A

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

Fig. 2. Optimized geometry structure of the Li-PG with H2molecules adsorption. Red digits represent the correspond-ing bond length of H—H.

¯ad )及表 1    H2吸附能(DEad)、平均吸附能( ∆E

Li所带的电荷(q)

Table 1.    Adsorption energies (DEad) and average

¯ad )of H2 molecules, and theadsorption energies ( ∆Echarge of Li (q).

n(H2)

1

DEad /eV

¯ad ∆E

2–0.239–0.2741.40

3–0.249–0.2661.56

(c)

0.760 A

4–0.044–0.2251.56

–0.308–0.3081.23

(b)

/eV·H2–1

q /e

(a)

Top view

0.763 A

0.763 A

Side view

0.767 A

0.765 A0.7 A

图 3    Na-PG吸附H2分子优化后的几何结构图(红色数

字表示H—H键长)

Fig. 3. Optimized geometry structure of the Na-PG with H2molecules adsorption. Red digits represent the correspond-ing bond length of H—H.

068802-4

物 理 学 报   Acta  Phys.  Sin.   Vol. 69, No. 6 (2020)    068802

荷, 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 Å.

子s轨道出现在费米能级下约–9.1 eV处, Na原子的3s或2p轨道和氢分子的s轨道均没有发生轨道耦合现象. K原子的s和p轨道与H2分子的s轨道也无明显的轨道耦合现象(见图S6). 这说明H2分子吸附在碱金属修饰的PG体系上不是通过Kubas相互作用机制, H2分子主要通过极化机制以及弱范德瓦耳斯相互作用吸附在碱金属修饰的PG上.

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

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

n(H2)

1

DEh1/eVDEh2/eV

4

(a)

2–0.262–0.117

3–0.215

PDOSH2 1sLi sC 2sC 2p

–0.202–0.234

–0.166

2

图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分子.

PDOS20-20

-15

-10

-5

0

5

Energy/eV

4

(b)3

Na sNa pH2 1sC 2sC 2p

1

Top view

(a)(b)(c)

0.762 A

0-20

0.762 A

0.763 A

-15-10-5Energy/eV

0510

Side view

0.763 A0.763 A

图 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.

0.761 A

为了进一步了解氢的吸附机理, 计算分析了

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

Fig. 4. Optimized geometry structure of the K-PG with H2molecules. Red digits represent the corresponding bond 

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分子间无电荷转

length of H—H.

3.2.2    氢分子的吸附性能

图5和图S6 (见补充材料(online))描述了Li-PG, Na-PG, K-PG吸附1个H2的分波态密度图. 氢分子的s轨道和Li原子2s轨道间没有发生轨道耦合现象(见图5(a)). 从图5(b)看出, H2分

068802-5

物 理 学 报   Acta  Phys.  Sin.   Vol. 69, No. 6 (2020)    068802

移, 这与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分子的相互作用较弱.

–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不适合用作储氢材料.

Top viewSide view

(a)(b)

图 8    Li-PG结构吸附12个H2分子优化后的几何结构图Fig. 8. Optimized geometry structure of the Li-decorated 

图 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 istaken to be 0.01 e/Å3.

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, 而单

(a)(b)

图 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 

个Li原子的ZPE为0 eV, 所以ZPE对Li的结合能的贡献仅有0.035 eV (即: 如考虑ZPE校正,1个Li原子与PG间的结合能为–1.619 eV). 利用有限位移法计算的H2分子的ZPE为0.2 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 (黄

0.01 e/Å.

3

3.2.3    碱金属修饰的PG单胞两侧吸附氢

分子的情况

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分子的平均吸附能为

068802-6

物 理 学 报   Acta  Phys.  Sin.   Vol. 69, No. 6 (2020)    068802

色)脱离体系束缚, 这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.%.

4   结　论

本文深入研究了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分子相继

(a)

C

-3094

(b)

Total energy/eV-3095

300 K

脱离了Li原子的束缚, 在400 K时, 有58%的H2分子被释放, 在600—700 K范围内所有的H2分子全部被释放.

-3096

-3097

参考文献

[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. 1322876[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, EksikO, 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, WangC 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, WangC 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-3098

00.30.60.91.21.5

Simulation time/ps

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

Fig. 9. Equilibrated structure of the 12 H2-Li-PG (a) andfluctuations 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分子全部被释放.

068802-7

物 理 学 报   Acta  Phys.  Sin.   Vol. 69, No. 6 (2020)    068802

[16]Lu R F, Rao D W, Lu Z L, Qian J C, Li F, Wu H P, WangY 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, SunH, 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, FanZ 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. HydrogenEnergy 39 11047[24]Elyassi M, Rashidi A, Hantehzadeh M R, Elahi S M 2017Surf. Interface Anal 49 230[25]Vanderbilt D 1990 Phys. Rev. B 41 72[26]Perdew J P, Burke K, Ernzerhof M 1996 Phys. Rev. Lett. 773865[27]Ao Z M, Peeters F M 2010 Phys. Rev. B 81 2006[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 166[29]Kittel C (translated by Xiang J Z, Wu X H) 2009Introduction to Solid State Physics (Beijing: Chemical IndustryPress) p38 (in Chinese) [基泰尔 C 著 (项金钟, 吴兴惠 译)2009 固体物理导论 (北京: 化学工业出版社) 第38页][30]Zhou X F, Fang H Y, Tang C M 2019 Acta Phys. Sin. 68053601 (in Chinese) [周晓锋, 方浩宇, 唐春梅 2019 物理学报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 6101Hydrogen storage capacity of alkali metal atoms

decorated porous graphene*

Yuan Li -Hua 1)†    Gong Ji -Jun 1)    Wang Dao -Bin 1)    Zhang Cai -Rong 1)

Zhang Mei -Ling 1)3)    Su Jun -Yan 1)    Kang Long 2)

1) (School of Sciences, Lanzhou University of Technology, Lanzhou 730050, China)

2) (School of Material Science and Engineering, Lanzhou University of Technology, Lanzhou 730050, China)

3) (School of Nuclear Science and Technology, Lanzhou University, Lanzhou 730000, China)

( Received 7 May 2019; revised manuscript received 15 January 2020 )

Abstract

Porous graphene (PG), a kind of graphene-related material with nanopores in the graphene plane, exhibits

novel properties different from those of pristine graphene, leading to its potential applications in many fields.Owing to periodic nanopores existing naturally in the two-dimensional layer, PG can be used as an idealcandidate for hydrogen storage material. High hydrogen storage capacity of Li-decorated PG has beeninvestigated theoretically, but the effect of temperature on the stability of the H2 adsorbed on Li-PG has beennot discussed yet. In this paper, by using the first-principles method, the hydrogen storage capacity on alkalinemetal atoms (Li, Na, K) decorated porous graphene is investigated in depth with generalized gradientapproximation, and the effect of the temperature on the stability of the hydrogen adsorption system iselucidated by the ab initio molecular-dynamics simulation. The results show that the most favorable adsorptionsites of Li, Na and K are the hollow center sites of the C hexagon, and four alkaline metal atoms can beadsorbed stably on both sides of PG unit cell without clustering. Alkaline metal adatoms adsorbed on PGbecome positively charged by transferring charge to PG and adsorbed H2 molecules, and three H2 molecules can

 *  Project suppotted by the Natural Science Foundation of Gansu Province, China (Grant No. 17JR5RA123), the NationalNatural Science of Foundation of China (Grant No. 51562022), the Doctoral Development Fund of Lanzhou University ofTechnology, China (Grant No.061801), and the Hong Liu First-class Disciplines Development Program of LanzhouUniversity of Technology, China.

†  Corresponding author. E-mail:  yuanlh@lut.edu.cn

068802-8

物 理 学 报   Acta  Phys.  Sin.   Vol. 69, No. 6 (2020)    068802

be adsorbed around each alkaline metal atom. By analyzing the Mulliken atomic populations, charge densitydifferences and density of states of H2 adsorbed on Li-PG system, we find that the H2 molecules are adsorbed onalkaline metal atoms decorated graphene complex by attractive interaction between positively charged alkalinemetal adatoms and negatively charged H and weak van der Waals interaction. Twelve H2 molecules areadsorbed on both sides of PG decorated with alkaline metal atoms. The average adsorption energy of H2adsorbed on Li-PG, Na-PG and K-PG are –0.246, –0.129 and –0.056 eV/H2, respectively. It is obvious that thehydrogen adsorption capacity of Li-PG system is strongest, and the hydrogen adsorption capacity of K-PG isweakest, thus K-PG structure is not suitable for hydrogen storage. Furthermore, by the ab initio molecular-dynamic simulation, in which the NVT ensemble is selected but the external pressure is not adopted, the effectof temperature on the stability of H2 molecules adsorbed on Li-PG system is elucidated. The result shows thatthe configuration of Li-PG is very stable, H2 molecules are stably adsorbed around the Li atoms at lowtemperature, and some H2 molecules start to be desorbed from the Li atoms with the increase of temperature.At 200 K, H2 molecules begin to move away from Li atoms, and two H2 molecules escape from the binding ofthe Li atoms at 250 K. At 300 K, nine H2 molecules can be stably absorbed on both sides of Li-PG, and thegravimetric hydrogen storage capacity can reach up to 9.25 wt.%, which is much higher than the the USDepartment of Energy target value of 5.5 wt.% for the year 2017. With the increase of temperature, moreadsorbed H2 molecules are desorbed, seven H2 molecules can be desorbed at 400 K, and all H2 molecules arecompletely desorbed in a temperature range of 600–700 K.

Keywords: porous graphene, hydrogen storage, first-principles, molecular-dynamic

PACS: 88.30.R–, 81.05.ue, 71.15.Mb, 71.15.Pd                          DOI: 10.7498/aps.69.20190694

068802-9