揭示 Be(BH?)? 卓越的储氢潜力：超高质量储氢容量与常温脱附

《Computational Condensed Matter》：Unveiling the exceptional hydrogen storage potential of Be(BH 4) 2: Ultra-high gravimetric capacity and ambient-temperature desorption

　　

**摘要要点**

• 通过AIMD计算确认了Be(BH4)?在300 K下的热力学与相稳定性
• 固有的高延展性可防止循环吸放氢过程中材料粉化
• HSE06计算揭示了由强离子-共价键稳定的7.60 eV带隙
• 储氢容量达20.83 wt%和110.13 g/L，远超美国能源部2025年目标
• 297.90 K下的脱氢温度实现了常温氢释放

**引言**

能源是推动国民经济发展的根本催化剂。当前，科学界在电力生产和高效储能方面面临重大挑战。尽管过去几十年锂离子电池和超级电容器的发展为储能问题提供了重要解决方案[1]，但不断攀升的全球需求使得寻找替代传统化石燃料的可行方案变得更加紧迫。太阳辐射是一种普遍存在且威力强大的能源，为电力生产提供了可持续途径[2]。因此，材料科学研究已转向发现具有增强功能特性的新型物质。当前学术重点在于开发能够通过经济高效、可靠且高性能的机制将太阳能转化为电能的工程材料[3]。在各种可再生能源中，常被称为"绿色能源"的太阳能已占据主导地位，但其运行效率仍严格取决于光电器件的化学组成[4]。大量研究已探索了各类材料体系中多样的存储模式[[5], [6], [7], [8], [9]]，涵盖从超导到超级电容及热电应用等领域。当前的研究工作主要集中于开发环保且经济的制备技术，以最大化面向多领域工业应用的存储容量[[10], [11], [12], [13], [14], [15], [16], [17], [18], [19]]。

在此阶段，氢已成为现代能源体系中最具前景的"清洁能源载体"，成为科学界关注的焦点，用于太阳能等间歇性可再生能源的长期大规模存储[20,21]。氢具有卓越的能源密度（120–142 MJ/kg），按质量计约为传统化石燃料（汽油、柴油等）的三倍，且燃烧产物仅为水蒸气，这使其在迈向碳中和未来方面无与伦比[[21], [22], [23]]。特别是利用光催化水分解或电解技术与太阳能结合生产的"绿氢"，在工业部门和交通网络的脱碳过程中发挥着关键作用[20,22,24]。

然而，在全球范围内和工业级别上全面实现氢能经济面临的最大科学和技术障碍，是安全、高效且经济地存储和运输氢气的挑战[[25], [26], [27]]。分子氢极轻且高度扩散的特性限制了传统存储方法的实用性。目前广泛使用的高压气瓶（350–700 bar）和低温液态存储系统（?253 °C）不仅存在严重的安全风险，而且由于液化过程的高能耗和不可避免的蒸发损失，也不利于长期使用。为克服这一技术瓶颈，材料科学研究转向了在原子或分子水平上捕获氢气的固态存储技术。固态储氢方法凭借其高体积密度和质量密度、低压运行以及出色的安全性能，提供了比传统压缩或液态系统更可行、更可靠的替代方案[28,29]。在此背景下，能够通过可逆化学键将氢纳入其结构的金属氢化物和复杂合金，因其高存储密度而备受关注[21,22]。

特别是近年来研究的材料，如BeH?(H?)?[30]、Mg(B?H?)?[31]、Li?Al???B?H?[32]、Mg?FeH?[33]、KH?N?[34]、Mg?XH??（X = Ti, Mn, Fe）[35]、LiMH?（M = Sc, Ti且x = 3, 4, 5）[36]、Li?AlH?[37]、X?AlH?（X: Ca, Sr, Ba）[38]，因其灵活的晶格结构、可调的热力学性质和适宜的氢解吸动力学，在设计下一代储能材料方面显示出巨大潜力。优化这些先进材料的结构、力学和光电性质，是实现将太阳能转化为氢并加以存储的高性能多功能能源系统商业化的最关键步骤[39,40]。目前作为预印本发布的新兴工作进一步扩展了这一领域。XCo?H?（X = Ca, Be, Mg）氢化物[41]在此具有特殊相关性，因为该系列包含一个铍基成员：BeCo?H?达到4.16 wt%，但需要609.5 K才能释放氢，这表明仅凭轻阳离子并不能保证在有过渡金属亚晶格束缚氢时获得有利的存储热力学性质。快速发展的KX?H?和XYH?家族提供了进一步的当代基准：除KX?H?（X = Sc, V, Cr）研究[42]外，最近的预印本还研究了KX?H?（X = Fe, Co, Mn）[43]和正交晶系XYH?（X = Ca, K；Y = Ni, Co）氢化物[44]，后者同时建立了形成能、Born和声子稳定性，并报告了约3.91–3.96 wt%的容量。

快速增长的第一性原理工作目前通过替换碱土或碱金属阳离子，在固定结构原型中筛选过渡金属基氢化物。代表性的近期例子包括四方晶系X?FeH?（X = Ca, Mg, Ba）氢化物[45]和立方钙钛矿型Y?CoH?（Y = Ca, Ba, Mg, Sr）氢化物[46]，以及NaXH?/KXH?钙钛矿[47]和K?Ba???Ga?H?[48]和X?CaZnH?[49]双钙钛矿，所有这些研究都在单一晶格类型中结合了对结构、力学、电子、热力学、声子和储氢的分析。这些研究有效地界定了给定原型的性质窗口，但这些化合物的氢含量受限于其化学组成而非阳离子：最优成员达到4.60 wt%（Mg?FeH?）、5.32 wt%（Mg?CoH?）和5.12 wt%（K?LiAlH?），最重成员降至1.50–1.77 wt%，且氢通过金属d-H杂化或通过[AlH?]/[GaH?]络合物结合在金属或半导体背景上。

本研究处于这一设计空间之外。Be(BH4)?不含硼以外的更重元素，氢以共价[BH?]?四面体形式容纳于闭壳层Be2?中心，没有可用的d态进行金属-氢杂化。因此，我们探测了结晶复杂氢化物所能达到的上限，并确定了在这种独特键合机制下允许常温释放的机理。所提供的额外科学洞见有四方面。第一，该化合物占据的容量区间与过渡金属复杂氢化物可达到的区间在性质上有本质区别，结合了20.83 wt%与110.13 g/L。第二，我们表明其脱氢热力学由产物相的形成能而非金属-氢键强度控制，并量化这一贡献约为BeB?路径每个分子单位0.30 eV；在[MH?]、[MH?]或[AlH?]框架中不存在类似的产物驱动路径。第三，我们通过300 K下的从头算分子动力学验证，近常温的解吸温度并不意味着在运行条件下丧失结构完整性，从而明确解决了容量-稳定性权衡问题而非默认假设。第四，我们解析了三条竞争性分解路径而非单一假设反应，区分了理论氢含量与沿各路径可回收的比例，并结合释放20.83 wt%氢所隐含的更大体积变化来解释力学响应。总体而言，本研究是与寻找可直接部署的非毒性氢化物（如文献[[45], [46], [47], 50]和[49]所示例）相辅相成而非竞争的关系：它确立了性能上限及此类配方可能趋近的底层设计规则。

必须承认，铍基化合物固有的毒性对大规模商业处理及广泛移动部署构成了严峻挑战。尽管如此，确立Be(BH4)?的绝对理论极限提供了至关重要的热力学和力学基准。理解这些极端性能指标对于指导未来安全封装型固态系统的设计或激发合成低毒性结构类似物至关重要。

受高容量、低温储氢介质迫切需求的驱动，本研究对铍硼氢化物Be(BH4)?进行了全面的第一性原理研究。虽然Be(BH4)?是一种传统上通过BeCl?与LiBH?复分解反应高纯度合成的已知无机结晶固体，但其微观力学行为和固态存储的热力学适宜性尚未得到充分阐明。在本工作中，我们系统地绘制了其结构完整性、电子带行为及力学韧性，特别关注其固有延展性和相稳定性。此外，我们严格评估了其相对于美国能源部目标的实际储能指标，详细说明了其体积和质量容量，以及沿多条分解路径的模拟解吸温度。通过这种多角度的计算方法，我们旨在证明Be(BH4)?作为下一代氢能经济中高效常温候选材料的卓越潜力。

**计算方法**

Be(BH4)?的基本电子、结构和力学特性通过第一性原理密度泛函理论（DFT）计算进行了严格研究，计算在Materials Studio套件中的CASTEP程序包[51]中实施。理论框架采用Kohn-Sham形式[[52], [53], [54]]结合平面波赝势方法。为准确捕捉电子转移动力学，采用PBE参数化方案。

**结构性质**

结构表征确认Be(BH4)?系列结晶于高对称四方I4?cd空间群（No. 110），其中原子排列由非对称单元的Wyckoff坐标严格定义。在此框架中，Be阳离子占据16b一般位置(0.5896, 0.7992, 0.9992)，硼原子分布在两个不同的16b位：B1位于(0.6937, 0.8600, 0.8781)，B2位于(0.4548, 0.8353, 0.9914)。氢原子则

**结论**

向碳中和能源未来转变需要发现能够弥合氢存储理论前景与实际可行性之间差距的材料。本研究严格确认了铍硼氢化物（Be(BH4)?）是一种高性能候选材料，兼具极高的质量密度和室温动力学特性。结构与热力学分析证实，四方I4?cd相表现出深度稳定性。

**CRediT作者贡献声明**

Saadi Berri：撰写——原稿、软件、资源、方法论、概念化。?a?atay Yam????er：撰写——审阅与编辑、撰写——原稿、软件、方法论、调查。

**关于生成式AI和AI辅助技术在写作过程中使用的声明**

在本工作准备过程中，作者使用了人工智能来改善论文的语言和可读性。在使用这些AI辅助应用程序后，作者已根据需要修订和编辑了内容，并对发表内容的完整性负全部责任。

**利益冲突声明**

作者声明不存在可能影响本论文所述工作的已知竞争性财务利益或个人关系。

**致谢**

本论文中报告的数值计算完全/部分在土耳其科学技术研究所（TUBITAK）ULAKBIM高性能与网格计算中心（TRUBA资源）完成。

Saadi Berri | ?a?atay Yam????er