通过第一性原理模拟解析MgBH3（B=Mo和In）的储氢与功能性质

《RSC Advances》：Decoding the hydrogen storage and functional properties of MgBH3 (B = Mo and In) via first-principles simulations

　　

本研究对MgBH₃（B=Mo，In）的结构、储氢、光电、机械和热力学性能进行了完整的第一性原理（first principles）研究。研究人员通过形成能（formation energies）、Goldschmidt容差因子（Goldschmidt tolerance factor）和八面体因子（octahedral factor）确定了立方相的结构稳定性。MgMoH₃和MgInH₃的计算储氢容量分别为2.45 wt%和2.13 wt%，表明其在现实世界能量存储中具有合理的脱附性能。电子结构计算（GGA-PBE）确定了由于价带和导带重叠而产生的金属导电性（metallic conductivity）。光学分析表明，MgBH₃的高折射率（refractive index）和强吸收（absorption）结合其金属性质，确保了高效的电荷传输和晶格稳定性。弹性常数（elastic constants）满足Born稳定性判据（Born's criteria），证实了机械稳定性。这种稳定性与它们独特的延展性（ductile nature）相结合，确保了牢固的结构完整性并防止微裂纹产生，从而使重复的氢循环高度稳定。强弹性各向异性（elastic anisotropy）表明了氢化物钙钛矿（hydride perovskites）的方向依赖性。热力学评估：德拜温度（Debye temperature）、晶格热导率和最小热导率（lattice and minimum thermal conductivities）以及Grüneisen参数（Grüneisen parameter）提供了对声子输运和热容的深入理解，证实了其在高温下用作热障涂层（thermal barrier coatings）的效用。简言之，MgBH₃（B=Mo，In）氢化物是多功能材料，表现出适中的储氢能力、稳健的机械稳定性、良好的热管理特性以及独特的介电和金属光学响应。

论文解读文章

**研究背景与问题**
随着全球人口增长和科技进步，能源需求持续上升，化石燃料的过度消耗导致CO₂排放增加，氢能因其丰度高、无毒、高效和可持续性而成为替代能源。然而，氢能存储仍是关键挑战：高压压缩存在安全隐患，液化则造成蒸发损失。固体储氢材料，尤其是金属氢化物钙钛矿（metal hydride perovskites），因其良好的体积和重量密度而备受关注。镁基氢化物（如MgBH₃）具有轻质和可逆吸收/脱附特性，但受限于高形成焓、缓慢动力学和空气氧化倾向。此前对CaXH₃（X=Mn，Fe，Co）和AeVH₃（Ae=Be，Mg，Ca，Sr）等体系的研究已证实其储氢潜力，但镁基钙钛矿的探索仍需拓展。为此，研究人员利用第一性原理密度泛函理论（first-principles density functional theory，DFT）系统研究了MgBH₃（B=Mo，In）的结构、弹性、热力学、光电和储氢性质，旨在评估其在实际能量存储和电子器件中的应用前景。论文发表在《RSC Advances》。

**主要技术方法**
研究人员采用基于平面波赝势（plane-wave pseudopotential）的CASTEP程序包进行计算。使用广义梯度近似（generalized gradient approximation，GGA）下的Perdew-Burke-Ernzerhof（PBE）泛函及其修正形式PBEsol进行结构优化，并采用Broyden-Fletcher-Goldfarb-Shanno（BFGS）最小化算法确保能量稳定。截断能设为500 eV，布里渊区采样采用6×6×6的Monkhorst-Pack k点网格。电子自洽场（self-consistent field，SCF）迭代的收敛容差为2.0×10^-5 eV/原子。光学性质通过计算复介电函数（complex dielectric function）获得，并利用Kramers-Kronig变换提取实部和虚部。此外，采用Slack经验模型估算晶格热导率（lattice thermal conductivity，κ_ph），并利用Clarke模型计算最小热导率（minimum thermal conductivity，k_min）。

**研究结果**

**3.1 结构性质与能量学**
通过能量-体积优化确定了晶格参数：MgMoH₃的晶格常数为3.55 ?，体积44.72 ?³；MgInH₃的晶格常数为4.11 ?，体积69.41 ?³，后者因In离子半径更大导致晶格膨胀。形成能计算显示MgMoH₃为?0.769 eV/原子，MgInH₃为?0.475 eV/原子，均为负值，表明二者热力学稳定。总能量分别为?217.08 Ry和?189.30 Ry，进一步证实强内聚相互作用。

**3.2 几何稳定性**
Goldschmidt容差因子（t）、八面体因子（μ）和物理容差因子（τ）评估表明：MgMoH₃的t=0.997，μ=0.448，τ=3.756；MgInH₃的t=0.949，μ=0.519，τ=3.693，均在立方钙钛矿稳定范围内（t: 0.813–1.107，μ: 0.37–0.859，τ<4.18），确认结构几何稳定。

**3.3 储氢性能**
通过质量分数公式计算，MgMoH₃和MgInH₃的重量储氢容量分别为2.45 wt%和2.13 wt%。脱附温度约567 K（MgMoH₃）和350 K（MgInH₃），表明适中的热稳定性。体积储氢容量分别为112 kg m^-3和72 kg m^-3，MgMoH₃因晶格更紧凑而容量更高。与多种已报道的钙钛矿氢化物（如RbCaH₃、CsPtH₃等）相比，两种化合物的储氢性能具有竞争力，适用于固定式或中等容量储氢场景。

**3.4 机械性质**
弹性常数满足Born稳定性判据（C₁₁?C₁₂>0, C₁₁+2C₁₂>0, C₄₄>0），证实机械稳定。MgMoH₃的体模量B=117.09 GPa、剪切模量G=57.14 GPa、杨氏模量E=147.44 GPa；MgInH₃的对应值分别为47.54 GPa、16.65 GPa、44.73 GPa，后者更软。泊松比ν分别为0.29和0.34，均高于0.26的脆性-延性分界，显示延性特征；Pugh比B/G分别为2.05和2.85（>1.75），进一步确认延性。Cauchy压力均为正值，表明金属键合。弹性各向异性因子A分别为0.40和0.79，偏离1，表明强各向异性。可加工性指数μ_M分别为3.19和5.16，显示优良的加工性能。维氏硬度分别为6.12 GPa和2.89 GPa，属软至中等硬度范围，有利于防止氢循环中的脆性断裂。

**3.5 电子性质**
能带结构计算显示两种化合物在费米能级（Fermi level）处无带隙，价带和导带强烈重叠，确认金属性质。部分态密度（partial density of states，PDOS）分析揭示MgMoH₃中Mg-p、Mo-p/d和H-p轨道主导；MgInH₃中Mg-s/p、In-p/d和H-p轨道主导，且存在显著的p-d杂化。HSE06杂化泛函进一步确认金属行为的内在性，有利于电子转移和氢吸附/脱附动力学。

**3.6 光学性质**
在0–12 eV光子能量范围内系统研究了光学参数。介电函数实部ε₁(ω)在低能区较高，表明强极化；虚部ε₂(ω)在红外区快速下降，反映吸收减弱。折射率n(ω)在0 eV时超过10，在可见光区降至1.2–2.6。消光系数k(ω)在红外区最大，紫外区趋近零。吸收系数α(ω)在紫外区达到峰值：MgMoH₃为28.85×10⁴ cm^-1（9.50 eV），MgInH₃为18.65×10⁴ cm^-1（4.66 eV）。光学电导率和反射率谱显示典型的金属行为，证实强的光子-电子相互作用。

**3.7 热力学性质**
德拜温度Θ_D：MgMoH₃为565.32 K，MgInH₃为307.91 K，后者较低表明热膨胀更缓和，适合热障涂层（thermal barrier coatings，TBCs）。晶格热导率κ_ph分别为11.14 W m^-1 K^-1和1.68 W m^-1 K^-1，MgInH₃的较低值归因于重原子和晶格畸变导致的强声子散射。Grüneisen参数γ分别为1.71和2.05，反映平均非谐性。最小热导率k_min分别为1.26 W m^-1 K^-1和0.59 W m^-1 K^-1，表明高温下的热传输效率。

**结论与讨论**
研究人员通过DFT计算系统揭示了MgBH₃（B=Mo，In）的多功能性质。结论指出：两种化合物具有热力学稳定性（负形成能）、几何稳定性（容差因子在合理范围）、适中的储氢容量（重量密度2.45 wt%和2.13 wt%）、金属导电性、强光学吸收（～10⁴ cm^-1）、延性机械行为和高弹性各向异性。热力学分析表明其适合高温应用，如热障涂层。该研究拓展了对镁基氢化物钙钛矿的理论认识，为多功能氢化物材料的未来研究提供了指导。