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

技術要約：層状Y2GaBおよびY2InB MAB相の物理的性質における圧力駆動型の進化

問題提起
MAB相として知られる層状三元ホウ化物（borides）は、強固な面内金属-ホウ素結合と、より弱い層間相互作用が共存するナノ積層構造により、金属的な電子挙動とセラミックスのような剛性を併せ持つユニークな特性を示している。これまでの研究では、特定のMAB化合物が常圧下で特性評価されているが、静水圧およびAサイト置換（具体的にはGa対In）が、構造安定性、機械的異方性、電子トポロジー、および熱物理的性質にどのように影響するかについての統一的な理解は不足している。さらに、これまでのY2GaBおよびY2InBに関する計算的研究では負の周波数を持つフォノンモードが報告されており、それらの動力学的安定性と、計算パラメータ（単位格子対スーパーセルの扱い）に対する評価の感性に疑問を投げかけている。本研究は、0から20 GPaまでのY2GaBおよびY2InBの圧力依存的な進化を系統的に比較することで、これらのギャップに対処し、機械的に実行可能な圧力窓を確立し、構造・物性相関関係を明らかにすることを目的とする。

手法
本研究は、Perdew-Burke-Ernzerhof (PBE) 汎関数近似 (GGA) を用いたCASTEPコードによる、Kohn-Sham密度汎関数理論 (DFT) に基づく第一原理計算を用いている。主な手法的特徴は以下の通りである：

• 構造最適化： 六方晶 $P6_3/mmc$ 構造を、BFGSアルゴリズムを用いて静水圧下（0, 5, 10, 15, 20 GPa）で最適化した。
• 弾性安定性： 弾性定数は有限応力-歪み法によって導出された。安定性は、六方晶対称性のための圧力補正基準を用いて、均質歪み安定性とフォノン安定性を区別して評価した。
• 物性の導出： 弾性的に安定な状態において、多結晶弾性率（Voigt-Reuss-Hill平均）、音速、デバイ温度、および熱伝導率の推定値（ClarkeおよびSlackモデル）を算出した。
• 電子および光学解析： バンド構造、状態密度 (DOS)、およびフェルミ面を、スピン軌道相互作用 (SOC) の有無を用いてQuantum ESPRESSOで計算した。光学特性は、[100]および[001]方向の電場偏極に対して解明した。
• 結合および動力学： 電子密度差マップ、Mulliken/Hirshfeld電荷解析、および結合重なり個数を利用して結合を特徴づけた。フォノン分散は、単位格子および拡張されたスーパーセルの両方の構成について、有限変位法を用いて計算し、動力学的安定性と数値的感度を評価した。

主な結果

• 構造の進化： 両化合物とも圧力下で $P6_3/mmc$ 構造を維持している。圧縮は異方的であり、底面格子定数 ($a$) は軸方向のパラメータ ($c$) よりも大幅に収縮する。その結果、両者において $c/a$ 比が増大し、単位格子体積は単調に減少する（0–20 GPaにおいてY2GaBで約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において、圧縮は音速、デバイ温度、および推定融点を上昇させる。300 KにおけるSlack推定格子熱伝導率は、0 GPaでの0.605 W m⁻¹ K⁻¹から10 GPaでの1.543 W m⁻¹ K⁻¹へと上昇する。両化合物は非常に低い最小熱伝導率を示しており、比較的低い融点に制限されるものの、低温域での熱障壁用途への可能性を示唆している。
• 電子構造： 両物質はマルチバンド金属である。常圧において、Y2GaBは9つのフェルミ面シートを持ち、Y2InBは8つを持つ。スピン軌道相互作用 (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）として提示されており、明示的な電子－フォノンおよび非調和散乱計算が行われていないことを認めている。