化石能源对环境的污染破坏使得社会发展不可持续，人们正积极努力寻找新一代的绿色能源。在这一过程中，氢能因丰富的储量、环境友好的特性、高能量密度以及产生的唯一产物——水，而备受瞩目。虽然氢气生产途径多种多样，但氢气的储存仍是巨大的挑战。传统的高压氢气储存方法要求持续维持高压，由此带来的潜在风险不容忽视。此外，液态储氢虽是可行方法，但能源投入较大、效率相对较低。因此，固态氢气储存被认为是理想的解决方案。固态储存能够在不改变材料结构的前提下，吸收和释放大量的氢气。目前，亟待寻找一种储氢性能出色且能在吻合条件下释放氢气的固态储氢材料。中国科学院力学研究所研究员彭庆与广西大学教授欧阳义芳带领的团队，通过第一性原理计算，开发出具有广阔前景的二维固态储氢材料体系Ti-decorated Irida-Graphene（钛修饰的鸢尾花型石墨烯，简称TIG）。这一体系中，Irida-Graphene（鸢尾花型石墨烯，简称IG）是新型的类石墨烯材料，由三原子、六原子和八原子的碳环组成。研究发现，通过在IG上引入修饰性的钛（Ti）原子，其储氢性能可达7.7wt%。 该团队对钛原子在鸢尾花型石墨烯不同吸附位点的吸附能力进行计算研究，发现鸢尾花型石墨烯的六原子上方的中空位（Hollow）是最稳定的吸附位点。进一步，该研究通过在体系中逐个加入氢分子（H2），计算了连续吸氢能力。结果显示每个钛原子周围可吸附5个氢分子。由于鸢尾花型石墨烯是二维材料，具有上下两个表面，其两个表面可同时吸附氢分子。当钛原子填满鸢尾花型石墨烯的所有六原子环中心环位时，其可储存的氢气达7.7 wt%。根据6.5 wt%的储氢密度可以使氢燃料电池车续航500英里估计，7.7 wt%储氢密度的续航约为590英里（950公里） 应用Van't Hoff方程计算得知，TIG的平均放氢温度为524 K（253摄氏度）。电子结构分析结果表明，氢气分子与TIG之间为Kubas型键合。氢气分子在吸附后获得净电荷。这种电荷转移会导致H-H键距延长，但不会使分子与其结构分离。同时，在TIG储存氢气时，钛原子的3d轨道与纯TIG相比，接近费米能级的电子态密度有所减少，表明在吸附后电荷从钛原子的3d轨道转移到氢分子的1s轨道上。随着额外的氢气分子不断集成到系统中，钛原子的3d轨道上电子将继续转移到氢的1s轨道上，形成稳定吸附。 为了考量TIG材料的结构热稳定性，该研究还通过CI-NEB方法探究了修饰的钛原子的迁移能力，发现钛原子的扩散迁移能垒为5.0 eV。在平均解吸温度524 K下，钛原子的热能为0.68 eV，低于扩散能垒值。由于结构构型的稳定性以及高能垒值在反应温度下的限制，钛原子的迁移受到限制。这意味着在吸氢和放氢过程中，钛原子不会发生迁移，避免了金属团聚引发的储氢结构破坏，证明了该体系作为可行的储氢介质。科研人员还借助第一性原理分子动力学模拟确定了该材料在室温300 K和高温600 K下的热力学稳定性。上述理论预测为开发新型高效的储氢材料提供了新的可选方案。 相关研究成果以Stable and 7.7 wt% hydrogen storage capacity of Ti decorated Irida-Graphene from first-principles calculations为题，发表在International Journal of Hydrogen Energy上。研究工作得到国家自然科学基金和力学所“力英计划”的支持。 论文链接 （a）2×2 的IG超胞；（b）TIG的原子结构。钛原子位于H2中空位点的吸附能最低，结构最稳定。灰色和蓝色球体分别表示碳原子和钛原子。黄色、红色和绿色分别表示中空位、顶位和桥位钛原子吸附位点。通过计算结合能发现，钛原子可以吸附在中空位和顶位，结合能分别为T1 (-2.50 eV)、T2 (-2.17 eV)、H1(-2.33 eV)、H2(-2.84 eV)和H3（-2.77 eV)。Ti原子在六角碳环的中空位的结合能最低。在六角碳环上方装饰了钛原子后，六角碳环的C-C键从1.43 增加到1.46 。IG在钛原子修饰之前（a）和之后（b）的总态密度（TDOS）。费米级Ef被设为零。在IG中掺入钛原子后，可能态的密度得到提高，自旋向上和自旋向下的态密度的对称性消失，且系统变得具有磁性。碳原子的2p轨道在（a）IG和（b）TIG的投影态密度（PDOS）；钛原子3d轨道在（c）孤立状态和（d）TIG中的投影状态密度。随着体系中钛原子掺入后，碳的2p轨道在费米能级以下的状态数显著增加，说明电荷增加。同时，钛原子的3d轨道从离散态变为连续态，表明钛原子与IG之间发生结合。氢分子以不同方向放置于钛原子周围，并经过几何优化后的原子结构。（a）水平方向、（b）垂直方向和（c）倾斜方向。这些不同方向放置的氢分子的吸附能分别为-0.40、-0.13和-0.60 eV。由此可见，倾斜放置是最稳定的吸附状态。TIG吸附H2分子后的优化结构和吸附能。（a）-（e） 1~5 H2分子吸附结构。（f）在IG两侧均放置钛原子，每个钛原子吸附5个H2的结构。灰球、蓝球和红球分别代表碳、钛和氢原子。（g）H2连续吸附能。 各种金属修饰碳纳米材料的储氢性能。可观察到TIG的体系具有相对较高的储氢效率和合适的解吸温度。这表明用TIG是一种有价值的储氢潜在材料。氢原子的1s轨道和钛原子3d轨道的投影态密度（PDOS）。（a）孤立H2的氢1s轨道；（b）TIG+H2系统的氢1s轨道；（c）TIG的钛3d轨道；（d）TIG+H2系统的钛3d轨道。系统吸收H2后，氢的1s轨道在费米级附近的密度增加。钛的3d轨道则在费米水平以下下降。这说明吸氢后，钛原子的3d轨道的电荷向氢的1s轨道转移。TIG体系（a）在吸附第一个氢分子之前和（b）在吸附第一个氢分子之后的差分电荷密度。电荷密度图中黄色代表静电势浓度增加，蓝色代表浓度降低。根据Bader电荷分析，TIG中钛原子将1.13e的电荷转移到相邻的碳原子上，这是钛原子和IG之间形成强键的原因。 钛原子从六角中空位向周围的八角环中空位迁移扩散能垒。插图显示了扩散路径。初始状态的能量设定为0 eV，最终状态比初始状态高0.07 eV。在此扩散路径上的扩散能垒为5.0 eV。这保证了钛原子在服役条件下发生迁移而导致的金属聚集和结构破坏。第一性原理分子动力学模拟在NVT系综下TIG的总能量和温度随模拟时间变化的关系曲线。模拟温度分别为（a）300 K和（b）600 K。（a）和（b）中的插图分别是TIG在两个温度下热浴10 ps后的原子快照。热浴过程中，钛原子在其初始位置附近并且不会移开，这表明在高温条件下TIG仍能保持结构的稳定性。