过渡金属氢化物超导体Mg2TmH6（Tm = Rh, Pd, Ir, Pt）物理性质的第一性原理计算研究

《Advanced Physics Research》：Physical Properties of Transition Metal Hydride Superconductors Mg2TmH6 (Tm = Rh, Pd, Ir, Pt) by First-Principles Calculations

　　

本研究对Mg2TmH6（Tm = Rh, Pd, Ir, Pt）氢化物的结构、储氢潜力、电子、弹性、力学、热物理、超导及光学性质进行了全面的第一性原理研究。这些化合物表现出良好的力学和动力学稳定性，其中Mg2TmH6（Tm = Rh, Ir, Pt）具有热力学稳定性，而Mg2PdH6为亚稳态。其质量储氢容量范围为2.42至3.84 wt.%。所有化合物均表现为韧性、弹性各向异性，并展现出弱金属特性，这反映在弹性常数序列C11 > C12 > C44以及模量趋势Y > B > G上。Mg2PtH6具有最大的体积模量。相比之下，Mg2IrH6表现出最高的杨氏模量、优异的可加工性和良好的干润滑性。电子结构显示费米能级处的态密度较低，预测的超导转变温度为25–44 K。热物理计算揭示了它们的稳定性和高效热输运能力。光学研究表明其在红外和可见光区域具有高反射率，并具有强紫外吸收。光学电导率表明其在紫外光电器件、光电探测器和反射涂层方面具有应用潜力。这些结果表明，Mg2TmH6氢化物兼具良好的储氢性能、机械鲁棒性、超导性和多功能光学性质，使其成为能源存储、超导和先进光电子应用的有前景候选材料。

论文解读文章

一、研究背景与意义

氢能作为一种理想的能源载体，对于减少对进口石油的依赖具有重要意义，尤其对缺乏自然资源的国家而言。氢化物是由氢与金属或类金属结合形成的化合物，根据化学键性质可分为离子型、共价型、金属型和范德华型四类。氢化物在储氢系统、低温制冷器、催化以及核技术工业等领域具有重要应用价值。此外，氢化物还表现出高温超导特性，因其高储氢容量、原材料丰富、成本低廉和潜在的超导性能而成为有前景的材料。

在超导研究领域，二元金属氢化物如H₃S（155 GPa下T_C约203 K）、YH₉（200 GPa下T_C为253–276 K）、YH₆（110 GPa下T_C为251–264 K）和YH₁₀（250 GPa下T_C约326 K）等表现出高温超导特性。大量三元氢化物也被预测在压力下具有优异的超导性能，包括ScCaH₈和ScCaH₁₂（200 GPa下T_C分别约212和182 K）、CaYH₁₂（200 GPa下T_C约258 K）等。值得注意的是，镁基三元氢化物在压力下也表现出超导特性，如MgSiH₆（250 GPa下T_C约63 K）、MgScH₆（100 GPa下T_C约41 K）和MgVH₆（150 GPa下T_C为27.6 K）。

然而，氢化物超导领域面临的主要挑战之一是在低压力或环境压力下实现超导，而非实验中通常所需的极高压力。文献中仅有少数氢化物被提出在环境压力下表现出高温超导性。Sanna等人提出了一族镁基过渡金属氢化物Mg₂TmH₆（Tm = Rh, Ir, Pd, Pt）作为潜在的高温常规超导体，这些材料被预测在结构和热力学上均稳定，且与实验合成的Mg₂RuH₆密切相关。

以往的研究主要集中于这些化合物的结构性质、电子能带结构和超导转变温度T_C以及热电性质。然而，包括质量储氢容量和体积储氢容量、弹性常数、硬度、各向异性、声速、热学及光电行为在内的多项重要物理性质尚未被研究。因此，本研究首次利用密度泛函理论（DFT）全面探索Mg₂TmH₆的结构、电子、力学、储氢、热学和光电性质，并重新审视其超导性质以保持完整性。

二、主要技术方法

研究人员采用基于CASTEP代码的密度泛函理论（DFT）计算方法，使用平面波赝势方法和广义梯度近似（GGA）处理交换关联效应，采用超软赝势描述电子-离子相互作用。结构优化使用Broyden–Fletcher–Goldfarb–Shanno（BFGS）算法，布里渊区采样采用Monkhorst–Pack方案。弹性常数通过CASTEP中的应力-应变方法计算，热物理参数由弹性常数和弹性模量推导。各化合物的能量截断值分别为500、450、600和600 eV，k点网格分别为8×8×8、10×10×10、14×14×14和16×16×16。此外，研究人员还使用ELATE程序生成了弹性各向异性的三维可视化图。

三、研究结果

3.1 晶体结构：Mg₂TmH₆结晶于立方Fm-3m空间群（No. 225），晶格参数a在约6.63–6.68 ?的窄范围内，从Mg₂PdH₆到Mg₂PtH₆略有增加。晶胞体积约为292–298 ?3。Mg₂IrH₆表现出最短的Tm–H键距（约1.72 ?），而Mg₂PdH₆和Mg₂PtH₆的键长略长（约1.79–1.80 ?）。计算结果与参考数据高度一致，验证了计算方法的可靠性。

3.2 热力学稳定性和氢化性质：通过凸包分析，Mg₂RhH₆、Mg₂IrH₆和Mg₂PtH₆位于凸包上，因此热力学稳定；而Mg₂PdH₆位于凸包上方56 meV/原子处，被归类为亚稳态。所有化合物均表现出负的结合能和负的形成焓。质量储氢容量顺序为Mg₂RhH₆（3.84 wt.%）> Mg₂PdH₆（3.75 wt.%）> Mg₂IrH₆（2.45 wt.%）> Mg₂PtH₆（2.42 wt.%）。Mg₂IrH₆具有最低的结合能（?3.94 eV/原子）和最放热的形成焓（?2.55 eV），表明其热力学稳定化较强。

3.3 振动性质：声子色散计算显示整个布里渊区不存在虚频，确认了这些化合物的动力学稳定性。低频声学模式主要由较重的Mg和过渡金属原子主导，高频光学模式主要与氢振动相关，反映了组成原子间的质量差异。

3.4 弹性性质：所有化合物的弹性常数遵循C₁₁ > C₁₂ > C₄₄的顺序。正的柯西压力（P_C）表明这些化合物具有以金属键为主的韧性和延展性特征。弹性模量遵循Y > B > G的趋势。Mg₂PtH₆具有最大的体积模量（85.06 GPa），表明其抗体积压缩能力最强；Mg₂IrH₆表现出最高的杨氏模量（97.15 GPa），表明其刚度更高。所有化合物的泊松比均高于临界值，Pugh比值均低于临界阈值，确认了其韧性本质。Mg₂IrH₆具有最高的可加工性指数（3.24）和优异的干润滑性。弹性各向异性分析表明，Mg₂IrH₆具有显著更高的剪切各向异性（A^G = 14.60%）和通用各向异性指数（A^U = 1.71），表明其弹性性质具有明显的方向依赖性。三维可视化显示，Mg₂IrH₆的杨氏模量表面呈强烈裂片状星形，证实了其最大的刚度对比。

3.5 声学性质和材料硬度：Mg₂RhH₆和Mg₂PdH₆的密度为3.572–3.664 g/cm3，而Mg₂IrH₆和Mg₂PtH₆的密度显著更高（5.535–5.576 g/cm3），反映了Ir和Pt较大的原子质量。横波声速以Mg₂PdH₆最高（2.9 km/s），纵波声速从Mg₂RhH₆和Mg₂PdH₆的约5.6 km/s降至Mg₂IrH₆和Mg₂PtH₆的约4.8 km/s。硬度计算表明，Mg₂IrH₆在所有形式体系中均表现出最高的机械硬度，而Mg₂RhH₆最软。所有化合物的硬度值相对较小，表明它们容易被压缩、切割、弯曲或刮擦。

3.6 热物理性质：德拜温度（θ_D）遵循Mg₂PdH₆（476.92 K）> Mg₂RhH₆（464.02 K）> Mg₂IrH₆（429.39 K）> Mg₂PtH₆（411.50 K）的趋势。热膨胀系数从Mg₂RhH₆的5.73×10^?5 K^?1降至Mg₂IrH₆的4.24×10^?5 K^?1。Mg₂IrH₆具有最高的熔点（1577.51 K），表明其热稳定性最强。格律乃森参数（γ）范围为1.70–1.99，其中Mg₂RhH₆最大。最小热导率在Cahill模型中为1.34–1.57 W m^?1 K^?1，在Clarke模型中为0.94–1.10 W m^?1 K^?1。Mg₂IrH₆的Kleinman参数（ζ = 0.31）显著低于其他化合物（0.53–0.55），表明其内部应变响应不同。

3.7 电子性质：态密度（DOS）分析显示费米能级处存在非零态密度，确认了Mg₂TmH₆的金属性质。费米能级处总态密度值分别为Mg₂RhH₆的1.93、Mg₂PdH₆的1.43、Mg₂IrH₆的1.54和Mg₂PtH₆的1.19 states/eV。能带结构显示无带隙，确认了金属性质。Mg₂PtH₆有九条能带穿过费米能级，预期具有最高的电导率。Mulliken布居分析显示Mg原子带正电荷（+1.52至+1.83 e），过渡金属原子带负电荷（?1.70至?2.48 e），H原子带负电荷（?0.20至?0.23 e），电荷溢出参数极小（0.25%–0.40%），表明基组充分。

四、超导性质与结论

电子结构显示费米能级处态密度较低，预测的超导转变温度为25–44 K。热物理计算揭示了这些化合物的稳定性。光学研究表明其在红外和可见光区域具有高反射率，并具有强紫外吸收，光学电导率表明其在紫外光电器件、光电探测器和反射涂层方面具有应用潜力。

综合来看，Mg₂TmH₆氢化物兼具良好的储氢性能、机械鲁棒性、超导性和多功能光学性质，使其成为能源存储、超导和先进光电子应用的有前景候选材料。其中Mg₂IrH₆在机械硬度和热稳定性方面表现突出，而Mg₂RhH₆和Mg₂PdH₆具有更高的储氢容量。这些发现为后续实验合成和实际应用提供了重要的理论指导。