氢介导的立方钙钛矿SrLiH?在极端压缩下的结构与电子稳定性研究

《Dalton Transactions》：Hydrogen-mediated structural and electronic stabilization of cubic perovskite SrLiH3 under extreme compressionOpen Access

　　

本工作对立方钙钛矿SrLiH?在静水压力最高达100 GPa条件下的结构、力学、动力学、振动、电子及成键性质进行了全面的第一性原理研究。研究人员通过焓计算，将立方钙钛矿相的相对稳定性与八种竞争性高压候选结构进行了评估。结果表明，在整个研究的压力范围内，立方钙钛矿始终是热力学有利相，且未发生压力诱导的结构相变。广义Born稳定性判据的满足证实了其力学稳定性，而声子计算显示在100 GPa内不存在虚频，确立了该结构的动力学稳定性。红外活性及静默声子模的压力依赖性显示，与LiH?八面体相关的高频振动发生系统性硬化，表明压缩下Li–H相互作用增强。采用Heyd–Scuseria–Ernzerhof杂化泛函的电子结构计算预测，SrLiH?始终保持绝缘体特性，带隙从环境压力下的2.7 eV增大至100 GPa时的3.3 eV。基于量子理论的原子盆分析显示，随压力增加原子盆发生显著收缩，而计算得到的原子电荷仅呈现适度变化，并保持Sr和Li阳离子与H阴离子之间清晰的电荷分离。由于氢亚晶格占据原子体积的最大份额，其对整体可压缩性作出主导贡献。这些发现表明，立方SrLiH?以离子键为主的特性在压缩下得以保持，其压力响应主要由原子盆的集体收缩而非显著的压力诱导电荷再分布所支配。热力学、力学、振动、电子及拓扑分析的联合结果一致表明，立方钙钛矿SrLiH?在极端压缩下表现出显著的结构鲁棒性和增强的电子稳定性，使其成为理解离子型复合氢化物在高压下行为的理想模型体系。

论文解读：立方钙钛矿氢化物SrLiH?在极端压缩下的稳定性与电子响应

一、研究背景与问题

富氢材料因在能源应用、储氢以及极端条件下新型化学成键研究方面的巨大潜力而受到广泛关注。金属氢化物中，氢可表现出从常规氢阴离子（H?）到复杂压力诱导成键构型等多种化学行为。高压研究尤其受到富氢体系超导电性发现的推动，但压力并非普遍导致金属化或超导电性；在许多离子型氢化物中，压缩反而会增强局域化键合相互作用并维持绝缘行为。理解这一多样性的调控机制，对于建立控制氢基材料稳定性与电子响应的普适原理至关重要。

钙钛矿型结构因其结构灵活性、化学多样性以及在扩展骨架中容纳轻元素的能力，为探索上述效应提供了理想平台。氢化物钙钛矿结合了氧化物钙钛矿的结构特征与氢的独特电子性质，氢可占据阴离子位点并参与金属–氢网络。然而，许多氢化物钙钛矿在高压条件下的稳定性尚不清楚，竞争相可能在高压下变得有利。

三元氢化物SrLiH?因含有LiH?八面体单元的钙钛矿衍生结构而成为特别值得研究的对象。电正性的Sr和Li原子与氢共存，为研究压缩下离子键合和氢电子态的演化提供了理想体系。此前对立方SrLiH?的第一性原理研究主要集中于50 GPa以下的力学、弹性和热力学性质，且未涉及动力学稳定性、压力诱导相竞争、杂化泛函电子结构、振动特性及极端压缩下化学键演化等基本问题，压力范围也落后于当前前沿水平。因此，在高达100 GPa的压力下开展系统研究，对揭示氢介导的稳定化机制具有重要意义。

二、研究内容与结论概述

研究人员利用密度泛函理论（DFT），结合声子计算和量子原子分子理论（QTAIM）分析，对立方SrLiH?在100 GPa以内的行为开展了全面研究。通过总能量和焓计算考察了包括立方钙钛矿相和多种可能高压构型在内的九种候选结构，以确定随压力变化的相稳定性；利用弹性常数和声子色散评估了优选相的力学和动力学稳定性；研究了红外活性和静默振动模的压力演化；采用Heyd–Scuseria–Ernzerhof（HSE06）杂化泛函精确确定带隙的压力依赖性；最后通过QTAIM分析研究了压缩下电子密度的再分布和化学键本质。

研究发现，立方SrLiH?在整个压力范围内均为最稳定相，并保持力学和动力学稳定性至100 GPa。与富氢化合物常见的压力诱导金属化行为相反，SrLiH?的带隙从环境压力下的2.7 eV增大至100 GPa时的3.3 eV，表现出增强的电子稳定性。QTAIM分析揭示的氢上电子局域化增强，以及氢主导振动模的硬化，共同表明压缩强化了离子型Li–H骨架。这些发现为理解离子型氢化物的压力诱导稳定化机制提供了新见解。

三、主要技术方法

本研究采用基于DFT的第一性原理计算，结构、力学、动力学和振动性质通过Quantum ESPRESSO（QE）代码计算，电子结构和成键性质则使用Vienna Ab initio Simulation Package（VASP）代码进行研究。计算采用超软赝势和Perdew–Burke–Ernzerhof（PBE）交换关联泛函，平面波截断能为100 Ry（约1360 eV），总能量计算的布里渊区积分采用16×16×16的Monkhorst–Pack k点网格。稳定性评估通过比较0 K下的焓和300 K下包含振动自由能贡献的Gibbs自由能完成。电子能带结构和态密度采用HSE06杂化泛函计算。声子计算在谐波近似下使用密度泛函微扰理论（DFPT）完成，q点网格为4×4×4。QTAIM分析采用全电子电荷密度，由全势线性化加平面波（FP-LAPW）代码WIEN2k提供，确保原子电荷定义的一致性。

四、研究结果详述

（一）压缩下的结构稳定性

SrLiH?是理想的立方钙钛矿氢化物，空间群为Pm?3m，由角共享LiH?八面体三维网络构成。计算得到的晶格参数a=3.803(3) ?，与实验值3.833(3) ?相差不足1%，与以往计算值3.788–3.833 ?一致。体积–压力关系符合三阶Birch–Murnaghan状态方程，给出零压体积V?=55 ?3、零压体模量B?=43 GPa及体模量压力导数B?'=3.93。

热力学稳定性通过比较0 K下竞争相的相对焓和有限温度下包含振动自由能贡献的Gibbs自由能来评估。研究考虑了立方钙钛矿相和八种候选高压结构，包括正交钙钛矿（Amm2，HP1）、正交后钙钛矿（Cmcm，HP2）、六方钙钛矿（P6?/mmc，HP3）、LiNbO?型三角结构（R3c，HP4）、钛铁矿型三角结构（R?3，HP5）、KCuF?型四方结构（I4/mcm，HP6）、四方钙钛矿（P4mm，HP7）和三角钙钛矿（R3m，HP8）。HP7和HP8在结构优化后自发回复为立方相，表明这些畸变钙钛矿构型在研究压力范围内不是稳定的局域极小值，故被排除。焓计算表明，所有高压候选相的焓均高于低压相，其中HP1与低压相能量最接近，但低压相在100 GPa内始终为最低焓结构。振动自由能贡献分析进一步确认立方钙钛矿相在300 K下是热力学最稳定相。

（二）声子色散与态密度

声子谱在0 GPa和100 GPa下均无虚频，证实了立方相在整个压力范围内的动力学稳定性。环境压力下声子谱延伸至约36 THz，可划分为三个频率区域：8 THz以下主要对应重Sr原子的振动，8–12 THz的中间频率光学模主要来自Li振动并伴有少量H贡献，18–36 THz的高频光学模主要源于轻H原子的振动。100 GPa时最大声子频率增至约58 THz，表明晶格显著硬化。H振动模移至39–58 THz区域且峰强显著增加，Sr振动仍局域于8 THz以下低频区域。压力主要影响氢衍生光学模，证实压缩增强了含氢键的刚度。

Γ点光学声子分析给出振动表示Γvib=4T??+T??，其中三个T??模为红外活性模，一个T??模为静默模。红外活性T??模涵盖Slater-Last型（Sr相对LiH?八面体框架运动）、Slater型（Li相对氢八面体振荡）和Axe型（Li–H键伸缩主导）三种模式。所有光学模在0–100 GPa范围内连续硬化，最高频T??模从约1060 cm?1增至近1900 cm?1，最低频T??模仅从约280 cm?1增至480 cm?1，与重Sr原子振动对压缩不敏感一致。整个压力范围内无模式交叉或频率软化，进一步支持结构的力学鲁棒性。

（三）弹性常数与力学性质

弹性常数计算满足广义Born稳定性判据的所有条件：C???C???2P>0、C??+2C??+P>0、C???P>0，证明立方SrLiH?在100 GPa内力学稳定。C??从107.09 GPa增至564.56 GPa（增幅约427%），C??从10.75 GPa增至214.81 GPa（增幅约1897%），C??从45.36 GPa增至178.20 GPa（增幅约293%）。C??的快速增长表明压缩显著增强不同晶向间的耦合。

Hill近似下，体模量B_H从42.86 GPa增至331.39 GPa，与状态方程所得43 GPa吻合；剪切模量G_H从46.46 GPa增至176.86 GPa；杨氏模量E从102.40 GPa增至450.46 GPa。B_H/G_H比值从0.922增至1.873，根据Pugh判据（B/G>1.75为韧性），SrLiH?在低压下呈脆性，高压下向韧性转变。泊松比σ从0.101增至0.273，表明从方向性键合向中心力型相互作用转变。纵波速度从5964.8 m s?1增至9899.5 m s?1，横波速度从3971.4 m s?1增至5527.9 m s?1。

（四）能带结构与电子性质

HSE06杂化泛函计算表明，环境压力下SrLiH?为间接带隙绝缘体，带隙E_g=2.712 eV，价带顶（VBM）位于M点附近，导带底（CBM）位于Γ点附近。价带（?7 eV至0 eV）完全由H s态主导，阳离子贡献可忽略，体现了高度局域的氢化物（H?）电子特征。导带主要由Sr d和Li p态组成，混有少量H s态。该带隙值显著大于GGA/LDA计算结果（1.51–1.93 eV），体现了HSE06对带隙低估问题的改善。

100 GPa下，SrLiH?仍保持间接半导体特性，未发生压力诱导的半导体-金属转变。带隙从2.712 eV增大至3.337 eV，与富氢材料普遍金属化的趋势相反。VBM仍位于M点附近，CBM从Γ点附近移至R点。价带仍以H s态为主，但Li和Sr在价带区域的态发生明显展宽且强度略有增加，反映压缩下轻微的共价杂化增强，但强离子性本质得以保持。基本电子跃迁主要由H-1s→Li衍生态构成，凸显氢在决定电子性质中的关键作用。

（五）低压SrLiH?结构的成键性质

QTAIM拓扑分析在分子图中识别出三类键临界点：Li–H、Sr–H和H?H相互作用。Li–H键临界点处电子密度ρ_BCP=0.0183 e bohr?3，Laplacian为?2ρ=+0.0777 e bohr??；Sr–H相互作用ρ_BCP=0.0157 e bohr?3，?2ρ=+0.0406 e bohr??。两者均具有低电子密度和正Laplacian特征，符合闭壳层相互作用判据，表明Li–H和Sr–H键主要为离子性，源于负电性氢化物离子与电正性阳离子之间的静电吸引。Li–H键临界点电子密度略高，归因于Li较小的离子半径导致的更强局域静电相互作用。H?H键临界点ρ_BCP=0.0120 e bohr?3、?2ρ=+0.0143 e bohr??，排除真实H–H化学键的存在，属于晶体拓扑强加的弱次级相互作用。电子密度等值线图显示电荷强烈局域于氢原子周围，Bader原子盆表明氢盆容纳大部分价电子密度，Sr和Li盆明显缺电子，证实了从Sr和Li向氢亚晶格显著的电荷转移。

（六）立方SrLiH?的整体体量与可压缩性

通过QTAIM分析，环境压力下Sr、Li和三个H盆占单位胞体积55.89 ?3中的55.47 ?3。三个H盆占QTAIM总体积的约59.85%，Sr占33.45%，Li占5.96%。Sr、Li和H盆的局域体模量分别为41.28、39.90和40.24 GPa，局域可压缩性分别为0.02422、0.02506和0.02485 GPa?1，各原子盆的内在可压缩性无显著差异。各盆对宏观可压缩性的贡献主要由体积分数决定：Sr贡献0.00810 GPa?1，Li贡献0.00149 GPa?1，三个H盆合计贡献0.01487 GPa?1，总和给出全局可压缩性κ_global=0.02446 GPa?1，对应体模量约40.9 GPa，与宏观状态方程所得43 GPa一致。氢亚晶格对整体可压缩性的主导贡献源于其占优的体积分数，而非异常高的局域可压缩性。

QTAIM原子电荷分析显示，Sr电荷从0 GPa时的+1.4278e降至100 GPa时的+0.6523e，Li电荷从+0.8369e降至+0.5300e，每个H原子的电荷从?0.7549e增至?0.3941e。电荷中性条件q_Sr+q_Li+3q_H=0在整个压力范围内保持一致。尽管电荷数值变化显著，但各原子电荷符号在整个压力范围内不变：Sr和Li始终为正，H始终为负。Sr电荷减少约0.78e，Li减少约0.31e，每个H增加约0.36e，三者累积变化相互补偿。这表明压缩导致电子密度在原子盆间的渐进再分布，而非简单的氢亚晶格电荷积累。QTAIM电荷的平滑演化无突变或符号反转，表明电子响应是连续的，不存在压力诱导的电荷转移相变。

五、结论总结与讨论

本工作通过热力学、力学、振动、电子和拓扑分析的联合研究，一致证明立方钙钛矿SrLiH?在高达100 GPa的极端压缩下表现出优异的结构鲁棒性和增强的电子稳定性。立方钙钛矿相在整个压力范围内保持热力学、力学和动力学稳定，无压力诱导结构相变。与富氢化合物普遍金属化的趋势相反，带隙随压力增大而增宽，从2.71 eV增至3.34 eV，表明强离子型电子结构得以保持。QTAIM分析证实Sr和Li阳离子盆与H阴离子盆之间的电荷分离持续存在，氢亚晶格因占据最大体积分数而对整体可压缩性作出主导贡献。压缩机制以原子盆的集体收缩为特征，而非压力诱导的显著电荷再分布，这支持了SrLiH?以离子键合为主的本质。

这些结果确立了SrLiH?作为理解离子型复合氢化物压力诱导稳定化机制的优秀模型体系，为理解氢在高压下氢化物钙钛矿行为中的调控作用提供了重要见解。氢不仅通过机械方式稳定晶体结构，还在超高压下保持其半导体电子特征，这一发现突出了氢介导稳定化的双重作用。本工作为未来实验高压研究和拓展至相关富氢钙钛矿材料的研究奠定了坚实的理论基础。