Thermal Transport in Ag8TS6 (T= Si, Ge, Sn) Argyrodites: An Integrated Experimental, Quantum-Chemical, and Computational Modelling Study

技術要約：Ag8TS6 (T= Si, Ge, Sn) アルジロディットにおける熱輸送

問題提起
熱電材料は、効率を最大化するために低い格子熱伝導率（$\kappa_L$）を必要とするが、多くの高性能な候補物質は、低い $\kappa_L$ と高いイオン伝導率の間に相関関係があり、それがデバイスの安定性を損なう可能性がある。Agベースのアルジロディット（一般式 $A_{(11-n)/m}^{m+}T_n^+Q_6^{2-}$）は、極めて低い $\kappa_L$ と高いイオン伝導率で知られているが、これらの硫化物変種 $Ag_8TS_6$ ($T = Si, Ge, Sn$) における熱輸送を支配する具体的なメカニズムは、まだ完全には理解されていない。既存の$\kappa_L$ 予測計算モデルはトレードオフに直面している。例えば、従来のモデル（SlackやCahillなど）は全温度範囲において精度に欠けるか、あるいは値を過大評価する傾向がある一方で、厳密な第一原理手法はハイスループット・スクリーニングを行うには計算コストが高すぎる。化学結合、非調和性、およびこれら特定の硫化物における熱伝導とイオン伝導の独立性の相互作用を解明するには、実験データと高度なモデリングを組み合わせた包括的な調査が必要である。

手法
本研究では、実験的な合成と特性評価、およびマルチスケール計算モデリングを組み合わせた統合的なアプローチを採用している。

1. 実験的特性評価:

   • 合成: $Ag_8TS_6$ ($T = Si, Ge, Sn$) の試料を、高温固相反応法により合成した。
   • 構造解析: 相純度の確認および相転移の調査のため、100 Kから400 Kの範囲で粉末X線回折（XRD）を実施した。
   • 輸送特性測定: 熱拡散率はレーザーフラッシュ法（173–500 K）を用いて測定し、イオン伝導率は、イオンの寄与を分離するために電子ブロッキング電極を用いた電気化学インピーダンス分光法（233–303 K）により決定した。電子伝導性は無視できるレベルであることを検証した。
   • 弾性的性質: 超音波パルスエコー法を用いて音速を測定した。

2. 計算モデリング:

   • 量子化学: 電子状態の最適化および電子特性の計算には、密度汎関数理論（DFT/VASP）を用いた。結合解析には、LOBSTERコードを用いた結晶軌道ハミルトニアン人口（COHP）および結晶軌道結合指数（COBI）を用い、結合強度、共有結合性、および多中心相互作用を定量化した。
   • フォノン特性: 調和フォノンバンド構造、群速度、およびデバイ温度を計算した。非調和性を定量化するため、準調和近似（QHA）を用いてモード依存のグリュナイゼンパラメータを導出した。
   • 熱伝導率予測: 2つの異なるモデリング手法を開発し、比較した：
     • グリュナイゼンに基づく2チャネルモデル: Xiaの2チャネルモデル（フォノン＋ディフゾン寄与）と、グリュナイゼンパラメータおよびSlackの理論に基づきフォノン寿命を推定するBjergの手法を組み合わせた修正アプローチ。これにより、高価な非調和力定数の計算を回避している。
     • 機械学習原子間ポテンシャル（MLIP）: 基盤となるMLIP（MACE-MP-03b）を利用して、Wigner輸送方程式を解くフル2チャネル格子動力学アプローチ（Simoncelliら）を通じて熱伝導率を算出した。

主な貢献および結果

• 結合と非調和性: 化学結合解析により、$T-S$ 結合は強固で共有結合的である一方、$Ag-S$相互作用は弱く、フェルミ準位下の占有反結合状態によって特徴付けられることが明らかになった。これは、$Ag^+$ 配位環境の歪みと「ラトラー（rattler）様」の挙動を引き起こす。本研究では、弱い $Ag-Ag$相互作用や、高原子価の多中心結合（$S-Ag-S$）を含む顕著な結合の不均一性を特定しており、これらが強い非調和性を駆動している。
• フォノン特性: すべての化合物は、虚数フォノンモードを持たない動的に安定な構造を示した。低周波のアコースティックモードは $Ag^+$ の振動によって支配されており、その結果、低い音速（約1000–1500 m/s）と低いデバイ温度をもたらしている。低エネルギーモードに対して大きなグリュナイゼンパラメータが観察され、高い非調和性が裏付けられた。
• 熱伝導率モデリング:
  • グリュナイゼンベースのモデルとMLIPベースのアプローチの両方が、これら $Ag_8TS_6$ 化合物の極めて低い格子熱伝導率を予測することに成功し、特に200 K以上の温度域において実験データと良好な一致を示した。
  • 結果は、熱輸送がディフゾンチャネル（非伝搬振動モード）によって支配されていることを裏付けており、これはセレニド系アルジロディットにおける過去の知見と一致している。
  • 実験データにおける低温域（0–50 K）での理想モデルとの乖離は、点欠陥散乱および境界散乱（結晶粒界）に起因すると考えられる。これらの散乱機構を組み込んだ解析モデルにより、低温挙動を再現することに成功した。
• 伝導率のデカップリング: 実験的測定により、これらの材料において熱伝導率とイオン伝導率が独立していることが示された。イオン伝導率は233 Kから303 Kの間で1桁以上増加している一方で、熱伝導率はほぼ一定に保たれている。このことは、熱伝導とイオン伝導を独立して制御できることを示唆しており、熱電応用における極めて重要な知見である。
• 同族置換: テトラル原子（$Si, Ge, Sn$）の置換は、結合特性、音速、グリュナイゼンパラメータ、および熱伝導率に最小限の影響しか与えず、熱輸送特性がこのファミリー全体で堅牢であることを示している。

意義および主張
本論文は、結合解析と効率的な計算モデルを統合することにより、複雑な無機材料の熱伝導率を予測するための堅牢な枠組みを確立したと主張している。主な意義は以下の通りである：

1. 予測因子としての結合解析: 特異な結合環境（弱い $Ag-S$ 結合、反結合状態、および多中心相互作用）は、これらの材料の低い熱伝導率と高い非調和性を予測するのに十分である。
2. 効率的なハイスループット・スクリーニング: 提案されたグリュナイゼンベースの2チャネルモデルおよび基盤となるMLIPアプローチは、完全な第一原理非調和計算に代わる、計算コストの低い選択肢を提供する。これらの手法は、低温データの解釈において散乱効果（欠陥、境界）を考慮すれば、熱伝導率のハイスループット・スクリーニングに適している。
3. 独立したチューニング: 本研究は、Agベースのアルジロディットにおいて熱伝導率とイオン伝導率がデカップリングされているという知見を補強し、イオンの安定性を損なうことなく熱電性能を最適化できる可能性を示している。

著者らは、MLIPおよびグリュナイゼンモデルは実験と強い一致を示すものの、複雑な系に対する基礎的な機械学習ポテンシャルの限界を定義し、特定の組成空間に対する自動訓練手順を洗練させるためには、さらなる検証が必要であると控えめに述べている。