Pressure-Driven Evolution of the Physical Properties of Layered Y2GaB and Y2InB MAB Phases: A Comparative First-Principles Study

技术摘要：压力驱动下层状 Y2GaB 与 Y2InB MAB 相物理性质的演化

问题陈述
被称为 MAB 相的层状三元硼化物由于其纳米层状结构，展现出金属电子行为与类陶瓷刚性的独特结合；在这种结构中，强内的金属-硼键与较弱的层间相互作用并存。虽然以往的研究已经表征了特定的 MAB 化合物在常压下的性质，但对于静水压力和 A 位取代（特别是 Ga 与 In）如何影响其结构稳定性、机械各向异性、电子拓扑结构以及热物理性质，仍缺乏统一的理解。此外，之前的关于 Y2GaB 和 Y2InB 的计算工作报告了负频率声子模式，这引发了对其动力学稳定性以及这些评估对计算参数（单元晶胞与超晶胞处理）敏感性的疑问。本研究通过系统地比较 Y2GaB 和 Y2InB 在 0 至 20 GPa 范围内的压力依赖演化，旨在建立机械可行性的压力窗口，并阐明结构-性质关系，从而填补这些空白。

方法论
本研究采用基于 Kohn-Sham 密度泛函理论（DFT）的第一性原理计算，使用 CASTEP 代码及 Perdew-Burke-Ernerthof (PBE) 广义梯度近似 (GGA)。关键方法论特征包括：

• 结构优化： 在静水压力（0, 5, 10, 15, 和 20 GPa）下，使用 BFGS 算法对六方 $P6_3/mmc$ 结构进行优化。
• 弹性稳定性： 通过有限应力-应变法导出弹性常数。使用针对六方对称性的压力修正准则进行稳定性评估，区分了均匀应变稳定性和声子稳定性。
• 性质推导： 计算了弹性稳定状态下的多晶模量（Voigt-Reuss-Hill 平均值）、声速、Debye 温度以及热导率估算值（采用 Clarke 和 Slack 模型）。
• 电子与光学分析： 使用 Quantum ESPRESSO 计算了带结构、态密度 (DOS) 和费米面（包含及不包含自旋-轨道耦合 SOC）。针对 [100] 和 [001] 方向的电场极化解析了光学性质。
• 成键与动力学： 使用电子密度差图、Mulliken/Hirshfeld 电荷分析以及键重叠布居数来表征成键。使用有限位移法对单元晶胞和扩展超晶胞配置进行声子谱计算，以评估动力学稳定性及其数值敏感性。

主要结果

• 结构演化： 两种化合物在压力下均保持 $P6_3/mmc$ 结构。压缩具有各向异性，基面晶格参数（$a$）的收缩比轴向参数（$c$）更为显著。因此，两种化合物的 $c/a$ 比值均有所增加，且单元晶胞体积单调减小（Y2GaB 在 0–20 GPa 内约为 21.8%，Y2InB 约为 20.3%）。
• 弹性稳定性与机械性质：
  • 稳定性窗口： Y2GaB 在 0, 5, 和 10 GPa 下满足压力修正的弹性稳定性准则，但在 15 和 20 GPa 时变得不稳定。Y2InB 仅在常压（0 GPa）下具有弹性稳定性。
  • 机械响应： 在其稳定范围内，Y2GaB 的体积模量、剪切模量和杨氏模量随压力增加而增大。两种化合物均被预测为具有延展性（Pugh 比率 $G/B < 0.57$）。Y2GaB 在 0 GPa 时表现出高弹性各向异性（$A^U = 5.25$），该值随压力增加而降低，而 Y2InB 在常压下的各向异性较小。
  • 热物理性质： 对于稳定的 Y2GaB，压力增加了声速、Debye 温度和估算的熔融温度。Slack 估算的 300 K 下晶格热导率从 0.605 W m⁻¹ K⁻¹ (0 GPa) 上升至 1.543 W m⁻¹ K⁻¹ (10 GPa)。两种化合物均表现出极低的最小热导率，暗示其在低温下具有作为热障应用材料的潜力，但受限于其相对较低的熔点。
• 电子结构： 两种材料均为多带金属。Y2GaB 在常压下有九个费米面片，Y2InB 有八个。自旋-轨道耦合 (SOC) 会修正局部能带色散，但不会打开能隙或改变费米面交点数量。在 0 至 10 GPa 之间未观察到 Lifshitz 相变（费米面拓扑重构），尽管中心口袋显著缩小。费米能级附近的电子态由 Y-4d 轨道主导，B-2p 和 A 位 p-轨道的贡献次之。
• 成键特征： 分析显示 Y 与 B 之间存在混合离子-共价键，其特征为 B 周围电荷聚集而 Y 周围电荷耗尽。A 位元素（Ga/In）表现出更弥散的金属键。
• 光学性质： 两种化合物均表现出金属 Drude 响应，并在 [100] 和 [001] 极化方向之间表现出强烈的光学各向异性。它们在可见光区域是高效的反射体，在紫外区是吸收体。压力诱导了光学特征的非均匀偏移，通常会提高 Y2GaB 的特征能量。
• 声子动力学： 单元晶胞声子计算显示，两种化合物在 0 和 10 GPa 下均具有正频率，表明在此限制性处理下具有动力学稳定性。然而，扩展超晶胞计算揭示了虚频模式，凸显了低频声子稳定性对晶胞大小和力常数处理的敏感性。值得注意的是，尽管 Y2InB 在 10 GPa 时在单元晶胞声子谱中显示出正频率，但其在弹性上是不稳定的，这强调了弹性稳定性与动力学稳定性准则之间的区别，尤其是对超晶胞尺寸的敏感性。

意义与主张
本文建立了一个理解 A 位取代（Ga 对比 In）和静水压力如何调节 MAB 相性质的比较框架。其主要贡献在于：

1. 定义稳定性窗口： 它确定了 Y2GaB 在机械上可行的压力范围（0–10 GPa），并明确了根据所采用的准则，Y2InB 仅在常压下是稳定的。
2. 解耦稳定性类型： 研究明确展示了弹性稳定性（对均匀应变的抵抗力）与动力学稳定性（对原子位移的抵抗力）是不同的，必须分别进行评估，特别是注意到声子计算对超晶胞尺寸的敏感性。
3. 结构-性质关系： 研究详细说明了压力如何系统地增强 Y2GaB 的刚性、热传输和光学响应，同时降低其弹性各向异性，为未来的实验验证提供了基准。
4. 研究范围限定： 作者明确指出，目前尚不存在这些化合物的实验确立的 MABene 衍生物；本研究侧重于母相，这是对未来评估低维衍生物所必需的基础。热导率数值是以模型估算值（Slack/Clarke）呈现，而非精确的输运预测，因为文中并未进行显式的电子-声子及非谐散射计算。