2025-11-21T07:28:22.047250

Thermoelectric transport properties of electron doped pyrite FeS2

Mukherjee, Subedi
Pyrite FeS$_2$ has been investigated for a wide range of applications, including thermoelectrics due to previous observation of large thermopower at room-temperature. However, the values of thermopower reported in the literature is extremely sensitive to the nature of sample -- whether they are natural or lab grown, bulk crystals or thin films -- and an ambiguity in the magnitude and sign of thermopower of pure FeS$_2$ exists. Variation in the magnitude of room-temperature thermopower has also been observed in Co-doped samples. Therefore, it is of interest to clarify the intrinsic thermopower of this system that could be measured in more pure samples. In this paper, we investigate the thermoelectric properties of Co-doped FeS$_2$ using first principles calculations. We apply three different doping schemes to understand the effect of electron doping in FeS$_2$, namely explicit Co-substitution, jellium doping and electron addition within rigid band approximation (RBA) picture. The calculated thermopower is less than $-50$ $μ$V/K for all values of Co doping that we studied, suggesting that this system may not be useful in thermoelectric applications. Interestingly, we find that RBA substantially overestimates the magnitude of calculated thermopower compared to the explicit Co-substitution and jellium doping schemes. The overestimation occurs because the changes in the electronic structure due to doping-induced structural modification and charge screening is not taken into account by the rigid shift of the Fermi level within RBA. RBA is frequently used in first principles investigations of the thermopower of doped semiconductors, and Co-substituted FeS$_2$ illustrates a case where it fails.
academic

電子ドープされたパイライトFeS₂の熱電輸送特性

基本情報

  • 論文ID: 2311.06135
  • タイトル: Thermoelectric transport properties of electron doped pyrite FeS2
  • 著者: Anustup Mukherjee、Alaska Subedi (CPHT、CNRS、École Polytechnique、Institut Polytechnique de Paris)
  • 分類: cond-mat.mtrl-sci
  • 発表日: 2023年11月13日
  • 論文リンク: https://arxiv.org/abs/2311.06135

要約

パイライトFeS₂は室温で観測される大きなゼーベック係数のため、熱電応用などで広く研究されている。しかし、文献に報告されているゼーベック係数の値は試料の性質（天然または実験室製造、バルク結晶または薄膜）に極めて敏感であり、純粋なFeS₂のゼーベック係数の大きさと符号に曖昧性が存在する。Co掺杂試料の室温ゼーベック係数にも変動が見られる。本論文は第一原理計算を用いてCo掺杂FeS₂の熱電特性を研究し、3つの異なる掺杂スキームを採用している：明示的なCo置換、ジェリウム掺杂、および剛体バンド近似(RBA)。すべてのCo掺杂濃度での計算ゼーベック係数は-50 μV/K未満であり、この体系が熱電応用で実用的でない可能性があることを示唆している。興味深いことに、RBAは掺杂による構造修飾と電荷スクリーニングが電子構造に及ぼす影響を考慮していないため、ゼーベック係数を著しく過大評価している。

研究背景と動機

問題の定義

  1. 実験データの不一致性：FeS₂のゼーベック係数測定値に巨大な差異が存在し、-115から524 μV/Kの範囲にあり、符号も一貫していない
  2. 試料依存性：熱電特性は試料タイプ（天然対人工、バルク対薄膜）に強く依存する
  3. 理論方法の正確性：既存の理論研究の多くは剛体バンド近似(RBA)を採用しているが、その正確性は疑わしい

研究の重要性

  • FeS₂は地球上に豊富で無毒な元素から構成されており、光電池、電池カソード、熱電などのエネルギー応用に潜在力を持つ
  • その本質的な熱電特性を明確にすることは、実際の応用価値を評価するために重要である
  • 掺杂半導体の熱電特性の理論研究に方法論的指針を提供する

既存方法の限界

  • 剛体バンド近似(RBA)は費米エネルギーの剛体移動のみによって掺杂効果をシミュレートする
  • 掺杂による構造変化と電荷スクリーニング効果を無視している
  • 異なる掺杂スキームの系統的比較が不足している

核心的貢献

  1. 3つの掺杂スキームの系統的比較：FeS₂の熱電特性計算における明示的なCo置換、ジェリウム掺杂、RBAの3つの方法の差異を初めて比較した
  2. RBAの限界の解明：RBAがゼーベック係数を著しく過大評価すること、および電子構造の観点からその理由を説明した
  3. Co掺杂FeS₂の熱電性能の明確化：すべての掺杂濃度でゼーベック係数が50 μV/K未満であることを証明し、この体系が熱電応用に不適切であることを示した
  4. 方法論的指針の提供：掺杂半導体の熱電特性の理論研究に重要な方法論的洞察を提供した

方法の詳細

計算方法

密度汎関数理論(DFT)フレームワークを採用し、Vienna Ab initio Simulation Package (VASP)を使用して計算を実施：

  • 交換相関汎関数：局所密度近似(LDA)
  • 疑似ポテンシャル法：投影増強波(PAW)法
  • エネルギーカットオフ：≥420 eV
  • k点グリッド：自己無撞着ループ8×8×8、状態密度計算24×24×24

3つの掺杂スキーム

1. 剛体バンド近似(RBA)

FeS₂の電子構造内で費米エネルギーを移動させることにより追加電子をシミュレートし、等価なCo掺杂値に対応させる。

2. ジェリウム掺杂

FeS₂単位胞に背景電荷として電子を明示的に追加し、電子構造計算を実施する。

3. 明示的なCo置換

Fe原子をCo原子で段階的に置換し、実際の化学置換効果を考慮する。

輸送特性の計算

BoltzTraP2コードを使用して、定緩和時間近似(CRTA)下で線形化ボルツマン輸送方程式を解く：

ゼーベック係数： S(μ,T)=ναjσαiS(\mu, T) = \frac{\nu_{\alpha j}}{\sigma_{\alpha i}}

ホール係数： RH(μ,T)=ναβkσαjσiβR_H(\mu, T) = \frac{\nu_{\alpha \beta k}}{\sigma_{\alpha j}\sigma_{i\beta}}

ここで輸送テンソルは以下のように定義される： σαβ(μ,T)=1Ω∫σαβ(ϵ)[−∂ϵfμ(ϵ)]dϵ\sigma_{\alpha\beta}(\mu, T) = \frac{1}{\Omega}\int \sigma_{\alpha\beta}(\epsilon)[-\partial_\epsilon f_\mu(\epsilon)]d\epsilon

実験設定

構造モデル

  • 結晶構造：パイライトFeS₂、空間群Pa3
  • 単位胞組成：4個のFe原子(4a位)、8個のS原子(8c位)
  • 掺杂濃度：1%、5%、10%、15%、25%、50%のCo掺杂を研究

磁性の考慮

実験がx≥0.1でFexCo1-xS₂が強磁性であることを示しているため、非磁性相と強磁性相の両方を考慮した。

構造最適化

  • 混合化合物の格子定数はVegardの法則に従う
  • 原子位置を弛緩させて原子力を最小化する
  • FeS₂の内部パラメータu=0.3821

実験結果

主要な結果

Fe₀.₇₅Co₀.₂₅S₂の熱電係数(室温)

掺杂スキーム非磁性相 (μV/K)強磁性相 (μV/K)
RBA-46.71-46.71
ジェリウム-14.68-14.69
Co置換-19.67-20.49

Fe₀.₅Co₀.₅S₂の熱電係数(室温)

掺杂スキーム非磁性相 (μV/K)強磁性相 (μV/K)
RBA-39.60-39.60
ジェリウム-17.52-17.52
Co置換-25.58-25.58

主要な発見

  1. RBAの系統的過大評価：すべての掺杂濃度でRBAはゼーベック係数の絶対値を著しく過大評価している
  2. 温度依存性の一貫性：3つの方法すべてがゼーベック係数が温度の増加とともに増加することを示している
  3. キャリアタイプ：すべての場合においてゼーベック係数が負の値であり、電子が主要なキャリアであることを示している
  4. 限定的な熱電性能：すべての掺杂濃度で室温ゼーベック係数は50 μV/K未満である

電子構造分析

RBAが過大評価する理由

状態密度(DOS)とバンド構造の分析により以下が明らかになった：

  • RBAでは：費米エネルギーが状態密度の急峻な領域に位置する
  • ジェリウム掺杂では：費米エネルギーがより広い肩部領域に移動する
  • Co置換では：費米エネルギーが比較的緩やかな領域に位置する

ゼーベック係数S ∝ d ln N(EF)/dEであるため、費米エネルギー付近のより急峻な状態密度がRBAでより大きなS値をもたらす。

異なる掺杂濃度での結果

組成RBA (μV/K)ジェリウム (μV/K)
Fe₀.₉₉Co₀.₀₁S₂-32.45-28.52
Fe₀.₉₅Co₀.₀₅S₂-10.47-10.24
Fe₀.₉Co₀.₁S₂-17.11-10.13
Fe₀.₈₅Co₀.₁₅S₂-25.41-11.75

関連研究

実験研究の現状

  • 歴史的データ：20世紀前半に報告されたゼーベック係数は-115から524 μV/Kの範囲
  • 現代研究：
    • Katoら：天然n型試料S = -300 μV/K
    • Willekeら：高移動度単結晶S = -320 μV/K
    • 薄膜研究：一般的に正のS値(60-75 μV/K)を報告

理論研究

  • 初期計算：GudelliらとHarranらはRBAを使用して最適化されたSが-400 μV/Kを超えることを予測
  • 本研究の貢献：異なる掺杂スキームを初めて系統的に比較し、RBAの限界を解明

結論と考察

主要な結論

  1. 限定的な熱電性能：Co掺杂FeS₂はすべての研究された掺杂濃度で室温ゼーベック係数が50 μV/K未満であり、優れた熱電材料の基準(>200 μV/K)をはるかに下回っている
  2. RBA方法の限界：剛体バンド近似はゼーベック係数を著しく過大評価する。これは以下を無視しているため：
    • 掺杂による構造修飾
    • 電荷スクリーニング効果
    • バンド幅展開と相対的な分裂の変化
  3. 方法論的意義：Co置換FeS₂はRBAが失敗する典型的な例であり、掺杂半導体の熱電特性研究に重要な警告を提供する

限界

  1. 計算方法の制限：LDAの使用はバンドギャップを低く見積もる可能性があり、熱電特性の正確性に影響する
  2. 掺杂濃度範囲：明示的なCo置換スキームは単位胞のサイズに制限され、非常に低い掺杂濃度を研究できない
  3. 緩和時間近似：定緩和時間近似の採用は輸送特性の正確性に影響する可能性がある

将来の方向性

  1. 計算方法の改善：より正確な交換相関汎関数(HSE06など)の使用
  2. 研究体系の拡張：他の遷移金属掺杂の効果の研究
  3. 実験検証：高品質試料の製造による理論予測の検証

深い評価

利点

  1. 方法論的革新：3つの掺杂スキームを初めて系統的に比較し、この分野に重要な方法論的ベンチマークを提供した
  2. 物理的洞察の深さ：電子構造の観点からRBAの過大評価の物理的メカニズムを深く説明した
  3. 計算の厳密性：磁性効果を考慮し、適切な計算パラメータと収束基準を使用した
  4. 実用的価値の明確性：Co掺杂FeS₂が熱電応用に不適切であることを明確に指摘し、無効な実験投資を回避する

不足

  1. 理論方法の限界：LDAは特に相関電子体系に対して十分に正確でない可能性がある
  2. 温度範囲の限定：主に室温特性に焦点を当てており、高温での系統的研究が不足している
  3. 掺杂機構の単一性：電子掺杂のみを研究し、正孔掺杂の可能性を考慮していない

影響力

  1. 学術的貢献：掺杂半導体の熱電特性の理論研究に重要な方法論的指針を提供した
  2. 実用的価値：研究者が有望でない体系にリソースを浪費することを回避するのに役立つ
  3. 再現性の強さ：計算の詳細が十分に記述されており、再現と拡張が容易である

適用シーン

  1. 理論研究：他の掺杂半導体の熱電特性研究に方法論的参考を提供する
  2. 材料スクリーニング：熱電材料のハイスループットスクリーニングで不適切な候補材料を除外する
  3. 実験指導：実験研究者が適切な掺杂戦略を選択するための理論的根拠を提供する

参考文献

論文はFeS₂の基本特性、実験研究、理論計算方法など各方面をカバーする38篇の重要な文献を引用しており、研究に堅実な文献基盤を提供している。特に注目すべきはXiらの最新実験研究の引用であり、この研究はFeS₂のキャリアタイプに関する以前の争点を解決したと思われる。

総合評価：これは高品質の理論研究論文であり、Co掺杂FeS₂の熱電特性を明確にするだけでなく、より重要なことに広く使用されているRBA方法の限界を解明し、この分野の理論研究に重要な方法論的貢献を提供している。研究結果は熱電材料の理論的スクリーニングと実験研究を指導するために重要な価値を持つ。