储氢材料第一性原理计算的研究进展

储氢材料第一性原理计算的研究进展

李兰兰　程方益　陶占良　陈　军

（南开大学新能源材料化学研究所，先进能源材料化学教育部重点实验室　天津３０００７１）

摘　要　综述了第一性原理计算在储氢材料研究中取得的成果和最新的进展。第一性原理计算在储氢材料

研究中的应用主要有以下４个方面：１）研究纳米结构的储氢性能；２）研究储氢材料中掺杂和缺陷的作用及

对储氢性能的影响；

３）研究储氢机理；４）确定氢化物的几何结构以及预测新型储氢材料。同时展望了第一性原理计算在储氢领域中的应用前景。

关键词　新能源，储氢材料，第一性原理计算

中图分类号：

Ｏ６４１ 文献标识码：Ａ 文章编号：１０００ ０５１８（２０１０）０９ ０９９８ ０６ＤＯＩ：１０．３７２４／ＳＰ．Ｊ．１０９５．２０１０．９０８０７２００９ １２ ０３收稿，２０１０ ０３ １５修回

国家自然科学基金（５０７７１０５６）、中国博士后科学基金（２００８０４３０７１２）和科技部“八六三”（２００９ＡＡ０５Ｚ１０６）资助项目

通讯联系人：陈军，男，教授，博士生导师；Ｅ ｍａｉｌ：ｃｈｅｎａｂｃ＠ｎａｎｋａｉ．ｅｄｕ．ｃｎ；研究方向：能源材料化学

第一性原理计算是指从５个基本物理常数（电子质量、电子电量、普朗克常数、光速和玻耳兹曼常数）出发来预测微观体系的状态和性质的方法。随着计算机科学的迅猛发展，以第一性原理计算为主的

科学计算逐渐成为继理论科学、实验科学之后，人类认识与征服自然的第三种科学方法［１ ２］。在材料科

学研究领域，第一性原理计算占据着越来越重要的地位，利用现代高速的计算机，通过第一性原理计算可以模拟材料的各种物理化学性质，深入理解材料从微观到宏观多个尺度的各类现象与特征，并对材料的物性和结构进行预言，从而为新材料的设计提供指导。在研究新能源材料的储氢材料方面，第一性原理计算可以从微观尺度给出储氢材料的性质，包括电子结构、成键特征、结合能、生成焓和脱附焓以及由

此估算的脱氢温度等物理化学性质［３ ４］。因而，第一性原理计算已经逐渐成为一种研究和探索新型储氢

材料的可靠工具。

第一性原理包括从头算、离散变分方法、超Ｈａｒｔｒｅｅ Ｆｏｃｋ方法等，从准确性和计算速度来看，从头算方法中的密度泛函（ＤＦＴ）理论应用较为广泛。目前，第一性原理计算在储氢领域应用主要体现在：１）研究纳米结构的储氢性能以及纳米结构对储氢性能的影响；２）研究储氢材料中掺杂和缺陷的作用及对储氢性能的影响；３）研究储氢机理，为进一步改善材料储氢性能和开发新的储氢体系提供理论指导；４）确定氢化物的几何结构以及预测新型储氢材料。本文从以上４方面综述了第一性原理计算在储氢材料开发中的研究现状。

１　研究进展

１．１　纳米结构储氢材料的模拟

纳米结构在储氢方面表现出了额外的优势，具有空间维度的材料，当其结构至少一维达到纳米范畴的尺寸时，通常表现出特殊的电子结构和物理化学性质。此外，纳米结构还具有高比表面、高稳定性和

多孔性等诸多的优点，有利于氢在表面和体相的扩散［５］。当纳米结构用作储氢材料时，表现出优良的吸

放氢动力学和热力学性能。第一性原理计算研究纳米结构储氢材料，不仅能预测材料潜在的储氢性能，更重要的是从电子和原子的角度给出纳米尺寸对储氢性能的影响以及氢化过程中材料性能的变化。目前，第一性原理计算模拟的储氢纳米结构主要有纳米颗粒、纳米管、纳米线、纳米笼以及纳米薄膜等。

Ｆｒｏｕｄａｋｉｓ小组［６］模拟了ＳｉＣ纳米管的储氢性能，结果表明，Ｈ２

与ＳｉＣ纳米管的结合能比碳纳米管增加２０％，表明了ＳｉＣ纳米管比碳纳米管更适于做储氢材料。Ｓｕｎ等［７］研究了Ｂ Ｎ纳米笼的储氢性能。在

第２７卷第９期

应用化学Ｖｏｌ．２７Ｎｏ．９２０１０年９月 ＣＨＩＮＥＳＥＪＯＵＲＮＡＬＯＦＡＰＰＬＩＥＤＣＨＥＭＩＳＴＲＹ Ｓｅｐ．２０１０

０Ｋ的温度下，Ｂ３６Ｎ３６纳米笼可以储存１８个Ｈ２

分子，储氢量质量分数达４％；而在室温下，Ｂ Ｎ纳米笼不能储存大量的氢，表明了Ｂ

Ｎ纳米笼不适于用作储氢材料。对Ｍｇ纳米薄膜的计算模拟则表明，在纳米尺度范围内，ＭｇＨ２纳米薄膜每减小２个元胞的厚度，脱附焓减小１０ｋＪ／ｍｏｌＨ２

［８］。２００５年ｄｅＪｏｎｇｈ研究小组［９］的计算结果深化了对纳米结构储氢的理解。他们以ＭｇＨ２

材料为研究对象，通过基于密度泛函的第一性原理计算，发现在一定的纳米尺度范围内，储氢材料的脱附焓随着材料尺寸的减小而减小，因而可以通过控制材料的尺寸来达到理想的储氢性能。进一步又用计算离子晶

体晶格能的Ｂｏｒｎ Ｈａｒｂｅｒ循环解释了这一纳米尺寸效应：如果将ＭｇＨ２

看成一个理想的离子晶体，就可以用Ｂｏｒｎ Ｈａｒｂｅｒ循环来计算晶体的晶格焓，ＭｇＨ２的晶格焓约是脱附焓的３６倍（晶格焓ΔＨｌａｔＭｇＨ２

＝２７１８ｋＪ／ｍｏｌ，脱附焓ΔＨｌａｔＭｇ＝

１４７ｋＪ／ｍｏｌ）。因此，晶格焓的微小变化将会对脱附焓产生较大的影响。而一般说来，体系中原子数目减小时，由于原子平均配位数的减小，将会导致晶格焓的减小，这种变化对于含原子数目较少的纳米材料更为明显。所以，对于纳米结构的储氢材料，可以通过控制材料的尺寸来控制材料脱氢反应的脱附焓，从而达到理想的储氢性能。他们的研究结果表明，当粒径小于１ ３ｎｍ时，ＭｇＨ２

纳米颗粒的脱附焓将会随粒径的减小而发生显著的变化。例如，纳米晶尺寸为０ ９ｎｍ的纳米颗粒脱附焓为６３ｋＪ／ｍｏｌＨ２

，对应一个比块体脱氢温度（６０３Ｋ）更低的脱氢温度４７２Ｋ。作者研究小组用密度泛函理论研究了Ｍｇ／ＭｇＨ２纳米线体系的储氢性能，发现Ｍｇ／ＭｇＨ２

纳米线体系随着纳米线直径的减小而变得不稳定，特别是当纳米线直径接近一个纳米时，氢原子更加趋向于吸附在

表面，使得体系吸放氢的焓变骤然降至块体体系的一半左右（＜

４０ｋＪ／ｍｏｌ）［１０ １１］。进一步研究表明，可以通过调节Ｈ原子的含量来改变Ｍｇ Ｈ的相互作用。此外，还发现Ｍｇ纳米线的富氢态是热力学可行的。

１．２　储氢材料中掺杂和缺陷的作用

对于大部分储氢材料，氢可以以原子的形式化学吸附（如金属和化学氢化物）储存或以分子的形式物理吸附（如碳结构、金属有机框架等）储存。前者电子从储氢材料转移到Ｈ２分子的反键轨道使Ｈ２分子键断裂，以原子的形式吸附，吸附较强，可以储存大量的Ｈ２

，但结合能过高，在室温下难以放氢；对于后者，Ｈ２与材料之间不发生电荷的转移，形成较弱的吸附，结合能太低以致室温条件下难以储存Ｈ２

。因此，研究者们提出对现有材料进行改性或开发一种新材料，使Ｈ２

的吸附介于范德华吸附和化学吸附之间。这种改性可以通过从结构上引入新的表面、晶界、缺陷和界面，或者从化学上引入一种新的元素或化合物。

１．２．１　碳基储氢材料中掺杂的作用　碳基纳米材料用于储氢领域受到广泛的关注。由于质轻、比表面积大、特殊的管状结构、多孔、热稳定性以及独特的机械性质，碳纳米管（ＣＮＴｓ）曾一度被认为是最具潜

力的储氢材料之一［

１２］。然而，实验和理论研究发现碳纳米管储氢难以满足实际的需求。其原因在于Ｈ２

与纳米管之间物理吸附的结合能过小，而形成化学吸附需要克服３ｅＶ的能垒。这样大的能垒导致了缓

慢的储氢动力学以及Ｈ２

化学吸附到碳纳米管上的过程极其缓慢［１３］。由于其化学惰性，未经修饰的碳纳米结构难以用作室温储氢。

一个增加化学活性的方法是引入修饰原子。修饰原子改变了纳米管的电子结构从而改变了Ｈ２

分子化学吸附的能垒。Ｌ

ｅｅ等［１４］报道了对于（５，５）碳纳米管上２个不同的吸附位，电子注入使化学吸附的能垒分别降低了０ ４９和０ ３０ｅＶ。电子注入对化学吸附能垒的影响对碳纳米管的直径和手性非常敏

感。Ｚｈｏｕ等［１５］的计算结果表明，电子注入对降低扶椅形碳纳米管的活化能非常有效，而空穴注入对降

低锯齿形纳米管的活化能有效。

直接电子注入也能够改变Ｆｅｒｍｉ能级附近的电子结构而有望改变吸附的强度。Ｙｏｏｎ等［１６］的理论

计算表明，电荷注入能够提高Ｈ２

分子与碳富勒烯之间在物理吸附距离的结合能，储氢量质量分数可以达到８ ０％。其原因在于带电的碳富勒烯使Ｈ２分子极化，增强了富勒烯与Ｈ２

分子之间的吸附。随后，他们研究了碱土金属掺杂的碳富勒烯的储氢性能［

１７］。结果表明，Ｃａ原子是非常好的掺杂原子（见表１），可以在碳富勒烯表面均匀分布，由于金属的空ｄ能级与富勒烯之间的电荷转移，在金属修饰的富勒烯周９９９　第９期李兰兰等：储氢材料第一性原理计算的研究进展