Thermoelektrische Transporteigenschaften von elektronendotiertem Pyrit FeS₂

Grundinformationen

• Papier-ID: 2311.06135
• Titel: Thermoelektrische Transporteigenschaften von elektronendotiertem Pyrit FeS₂
• Autoren: Anustup Mukherjee, Alaska Subedi (CPHT, CNRS, École Polytechnique, Institut Polytechnique de Paris)
• Klassifizierung: cond-mat.mtrl-sci
• Veröffentlichungsdatum: 13. November 2023
• Papierlink: https://arxiv.org/abs/2311.06135

Zusammenfassung

Pyrit FeS₂ wird aufgrund des bei Raumtemperatur beobachteten großen Seebeck-Koeffizienten umfassend für thermoelektrische Anwendungen untersucht. Allerdings sind die in der Literatur berichteten Seebeck-Koeffizient-Werte extrem empfindlich gegenüber Probeneigenschaften (natürlich oder im Labor hergestellt, Einkristall oder Dünnfilm), und es bestehen Mehrdeutigkeiten bezüglich Größe und Vorzeichen des Seebeck-Koeffizienten von reinem FeS₂. Auch bei Co-dotierten Proben variieren die Raumtemperatur-Seebeck-Koeffizient-Werte erheblich. Dieses Papier untersucht die thermoelektrischen Eigenschaften von Co-dotiertem FeS₂ mittels Berechnungen aus ersten Prinzipien unter Verwendung von drei verschiedenen Dotierungsschemen: explizite Co-Substitution, Jellium-Dotierung und Rigid-Band-Approximation (RBA). Die berechneten Seebeck-Koeffizienten bei allen Co-Dotierungskonzentrationen sind kleiner als -50 μV/K, was darauf hindeutet, dass dieses System möglicherweise nicht praktisch für thermoelektrische Anwendungen ist. Interessanterweise überschätzt die RBA den Seebeck-Koeffizienten erheblich, da sie die durch Dotierung verursachten Strukturmodifikationen und Ladungsabschirmungseffekte auf die elektronische Struktur nicht berücksichtigt.

Forschungshintergrund und Motivation

Problemdefinition

1. Inkonsistenz experimenteller Daten: Die gemessenen Seebeck-Koeffizient-Werte von FeS₂ zeigen enorme Unterschiede, die von -115 bis 524 μV/K reichen, mit inkonsistenten Vorzeichen
2. Probenabhängigkeit: Thermoelektrische Eigenschaften hängen stark vom Probentyp ab (natürlich vs. künstlich, Einkristall vs. Dünnfilm)
3. Genauigkeit theoretischer Methoden: Bestehende theoretische Studien verwenden häufig die Rigid-Band-Approximation (RBA), deren Genauigkeit fragwürdig ist

Forschungsbedeutung

• FeS₂ besteht aus reichlich vorhandenen und ungiftigen Elementen und hat Potenzial in Energieanwendungen wie Photovoltaik, Batteriekathoden und Thermoelektrik
• Die Klärung seiner intrinsischen thermoelektrischen Eigenschaften ist entscheidend für die Bewertung des praktischen Anwendungswerts
• Bietet methodische Orientierung für theoretische Forschung zu thermoelektrischen Eigenschaften dotierter Halbleiter

Einschränkungen bestehender Methoden

• Die Rigid-Band-Approximation (RBA) simuliert Dotierungseffekte nur durch starre Verschiebung des Fermi-Niveaus
• Ignoriert strukturelle Veränderungen und Ladungsabschirmungseffekte durch Dotierung
• Mangel an systematischem Vergleich verschiedener Dotierungsschemen

Kernbeiträge

1. Systematischer Vergleich von drei Dotierungsschemen: Erstmalige Gegenüberstellung der Unterschiede zwischen expliziter Co-Substitution, Jellium-Dotierung und RBA bei der Berechnung thermoelektrischer Eigenschaften von FeS₂
2. Offenlegung der Einschränkungen der RBA: Nachweis, dass RBA den Seebeck-Koeffizienten erheblich überschätzt, mit Erklärung aus elektronischer Strukturperspektive
3. Klärung der thermoelektrischen Leistung von Co-dotiertem FeS₂: Nachweis, dass der Seebeck-Koeffizient bei allen Dotierungskonzentrationen kleiner als 50 μV/K ist, was darauf hindeutet, dass dieses System nicht für thermoelektrische Anwendungen geeignet ist
4. Methodische Orientierung: Bereitstellung wichtiger methodischer Erkenntnisse für theoretische Forschung zu thermoelektrischen Eigenschaften dotierter Halbleiter

Detaillierte Methodenbeschreibung

Berechnungsmethoden

Dichtefunktionaltheorie (DFT)-Rahmen unter Verwendung des Vienna Ab initio Simulation Package (VASP):

• Austausch-Korrelations-Funktional: Lokale Dichthenäherung (LDA)
• Pseudopotential-Methode: Projector Augmented Wave (PAW)-Methode
• Energiecutoff: ≥420 eV
• k-Punkt-Gitter: 8×8×8 für selbstkonsistente Schleifen, 24×24×24 für Zustandsdichte-Berechnungen

Drei Dotierungsschemen

1. Rigid-Band-Approximation (RBA)

Simulation zusätzlicher Elektronen durch Verschiebung des Fermi-Niveaus in der elektronischen Struktur von FeS₂, entsprechend äquivalenten Co-Dotierungswerten.

2. Jellium-Dotierung

Explizites Hinzufügen von Elektronen als Hintergrundladung in der FeS₂-Einheitszelle mit anschließender Berechnung der elektronischen Struktur.

3. Explizite Co-Substitution

Schrittweise Ersetzung von Fe-Atomen durch Co-Atome unter Berücksichtigung tatsächlicher chemischer Substitutionseffekte.

Berechnung von Transporteigenschaften

Verwendung des BoltzTraP2-Codes zur Lösung der linearisierten Boltzmann-Transportgleichung unter der Constant Relaxation Time Approximation (CRTA):

Seebeck-Koeffizient: $S(\mu, T) = \frac{\nu_{\alpha j}}{\sigma_{\alpha i}}$

Hall-Koeffizient: $R_H(\mu, T) = \frac{\nu_{\alpha \beta k}}{\sigma_{\alpha j}\sigma_{i\beta}}$

wobei die Transporttensoren definiert sind als: $\sigma_{\alpha\beta}(\mu, T) = \frac{1}{\Omega}\int \sigma_{\alpha\beta}(\epsilon)[-\partial_\epsilon f_\mu(\epsilon)]d\epsilon$

Experimentelle Einrichtung

Strukturmodell

• Kristallstruktur: Pyrit FeS₂, Raumgruppe Pa3
• Einheitszellenzusammensetzung: 4 Fe-Atome (4a-Position), 8 S-Atome (8c-Position)
• Dotierungskonzentrationen: Untersuchung von 1%, 5%, 10%, 15%, 25%, 50% Co-Dotierung

Magnetische Überlegungen

Berücksichtigung sowohl nicht-magnetischer als auch ferromagnetischer Phasen, da Experimente zeigen, dass FexCo₁₋ₓS₂ für x≥0,1 ferromagnetisch ist.

Strukturoptimierung

• Gitterkonstanten von Mischverbindungen folgen dem Vegard-Gesetz
• Relaxation von Atompositionen zur Minimierung von Atomkräften
• Interner Parameter von FeS₂: u=0,3821

Experimentelle Ergebnisse

Hauptergebnisse

Seebeck-Koeffizient von Fe₀,₇₅Co₀,₂₅S₂ (Raumtemperatur)

Seebeck-Koeffizient von Fe₀,₅Co₀,₅S₂ (Raumtemperatur)

Wichtigste Erkenntnisse

1. Systematische Überschätzung durch RBA: Bei allen Dotierungskonzentrationen überschätzt die RBA den Absolutwert des Seebeck-Koeffizienten erheblich
2. Konsistente Temperaturabhängigkeit: Alle drei Methoden zeigen, dass der Seebeck-Koeffizient mit steigender Temperatur zunimmt
3. Ladungsträgertyp: Der Seebeck-Koeffizient ist in allen Fällen negativ, was darauf hindeutet, dass Elektronen die Hauptladungsträger sind
4. Begrenzte thermoelektrische Leistung: Der Raumtemperatur-Seebeck-Koeffizient ist bei allen Dotierungskonzentrationen kleiner als 50 μV/K

Analyse der elektronischen Struktur

Gründe für die Überschätzung durch RBA

Durch Analyse der Zustandsdichte (DOS) und Bandstruktur wurde festgestellt:

• In RBA: Das Fermi-Niveau befindet sich in einem steilen Bereich der Zustandsdichte
• In Jellium-Dotierung: Das Fermi-Niveau verschiebt sich in einen breiteren Schulterbereich
• In Co-Substitution: Das Fermi-Niveau befindet sich in einem relativ flachen Bereich

Da der Seebeck-Koeffizient S ∝ d ln N(EF)/dE, führt eine steilere Zustandsdichte in der Nähe des Fermi-Niveaus in der RBA zu größeren S-Werten.

Ergebnisse bei verschiedenen Dotierungskonzentrationen

Verwandte Arbeiten

Aktueller Stand der experimentellen Forschung

• Historische Daten: In der ersten Hälfte des 20. Jahrhunderts berichtete Seebeck-Koeffizient-Werte reichen von -115 bis 524 μV/K
• Moderne Forschung:
  • Kato et al.: Natürliche n-Typ-Proben mit S = -300 μV/K
  • Willeke et al.: Hochmobilitäts-Einkristalle mit S bis zu -320 μV/K
  • Dünnfilmforschung: Allgemein berichtete positive S-Werte (60-75 μV/K)

Theoretische Forschung

• Frühe Berechnungen: Gudelli et al. und Harran et al. verwendeten RBA zur Vorhersage optimierter S-Werte über -400 μV/K
• Beitrag dieser Arbeit: Erstmaliger systematischer Vergleich verschiedener Dotierungsschemen mit Offenlegung der Einschränkungen der RBA

Schlussfolgerungen und Diskussion

Hauptschlussfolgerungen

1. Begrenzte thermoelektrische Leistung: Co-dotiertes FeS₂ zeigt bei allen untersuchten Dotierungskonzentrationen Raumtemperatur-Seebeck-Koeffizienten kleiner als 50 μV/K, weit unter dem Standard für ausgezeichnete thermoelektrische Materialien (>200 μV/K)
2. Einschränkungen der RBA-Methode: Die Rigid-Band-Approximation überschätzt den Seebeck-Koeffizienten erheblich, da sie folgende Faktoren ignoriert:
   • Strukturmodifikationen durch Dotierung
   • Ladungsabschirmungseffekte
   • Veränderungen in Bandverbreiterung und relativer Aufspaltung
3. Methodische Bedeutung: Co-substituiertes FeS₂ ist ein typisches Beispiel für das Versagen der RBA und bietet wichtige Warnungen für die Forschung zu thermoelektrischen Eigenschaften dotierter Halbleiter

Einschränkungen

1. Einschränkungen der Berechnungsmethode: Die Verwendung von LDA kann die Bandlücke unterschätzen und die Genauigkeit thermoelektrischer Eigenschaften beeinflussen
2. Bereich der Dotierungskonzentrationen: Das Schema der expliziten Co-Substitution ist durch die Größe der Einheitszelle begrenzt und kann sehr niedrige Dotierungskonzentrationen nicht untersuchen
3. Relaxationszeitnäherung: Die Verwendung der Constant Relaxation Time Approximation kann die Genauigkeit von Transporteigenschaften beeinflussen

Zukünftige Richtungen

1. Verbesserte Berechnungsmethoden: Verwendung genauerer Austausch-Korrelations-Funktionale (wie HSE06)
2. Erweiterung des Forschungssystems: Untersuchung der Effekte anderer Übergangsmetall-Dotierungen
3. Experimentelle Validierung: Herstellung hochqualitativer Proben zur Validierung theoretischer Vorhersagen

Tiefgreifende Bewertung

Stärken

1. Methodische Innovation: Erstmaliger systematischer Vergleich von drei Dotierungsschemen bietet wichtige methodische Benchmarks für das Feld
2. Tiefe physikalischer Einsichten: Tiefgreifende Erklärung des physikalischen Mechanismus der RBA-Überschätzung aus elektronischer Strukturperspektive
3. Rigorose Berechnungen: Berücksichtigung magnetischer Effekte mit angemessenen Berechnungsparametern und Konvergenzstandards
4. Klarer praktischer Wert: Klare Angabe, dass Co-dotiertes FeS₂ nicht für thermoelektrische Anwendungen geeignet ist, vermeidet ineffektive experimentelle Investitionen

Mängel

1. Einschränkungen theoretischer Methoden: LDA kann möglicherweise nicht ausreichend genau sein, besonders für stark korrelierte Elektronensysteme
2. Begrenzte Temperaturbereiche: Hauptfokus auf Raumtemperatureigenschaften mit mangelnder systematischer Hochtemperaturforschung
3. Einzelne Dotierungsmechanismen: Nur Elektronendotierung untersucht, Möglichkeit der Lochdotierung nicht berücksichtigt

Auswirkungen

1. Akademischer Beitrag: Bietet wichtige methodische Orientierung für theoretische Forschung zu thermoelektrischen Eigenschaften dotierter Halbleiter
2. Praktischer Wert: Hilft Forschern, Ressourcen bei nicht vielversprechenden Systemen zu sparen
3. Starke Reproduzierbarkeit: Ausreichend detaillierte Berechnungsbeschreibungen ermöglichen einfache Reproduktion und Erweiterung

Anwendungsszenarien

1. Theoretische Forschung: Bietet methodische Referenzen für thermoelektrische Eigenschaftsforschung anderer dotierter Halbleiter
2. Materialauswahl: Ausschluss ungeeigneter Kandidaten bei Hochdurchsatz-Screening thermoelektrischer Materialien
3. Experimentelle Orientierung: Bietet theoretische Grundlagen für Experimentatoren bei der Auswahl geeigneter Dotierungsstrategien

Referenzen

Das Papier zitiert 38 wichtige Literaturquellen, die alle Aspekte von FeS₂-Grundeigenschaften, experimenteller Forschung und theoretischen Berechnungsmethoden abdecken und eine solide Literaturgrundlage für die Forschung bieten. Besonders hervorzuheben ist die Zitierung der neuesten experimentellen Arbeiten von Xi et al., die scheinbar frühere Kontroversen über den Ladungsträgertyp in FeS₂ gelöst haben.

Gesamtbewertung: Dies ist ein hochqualitatives theoretisches Forschungspapier, das nicht nur die thermoelektrischen Eigenschaften von Co-dotiertem FeS₂ klärt, sondern noch wichtiger die Einschränkungen der weit verbreiteten RBA-Methode offenlegt und wichtige methodische Beiträge für theoretische Forschung in diesem Bereich leistet. Die Forschungsergebnisse haben wichtigen Wert für die Orientierung theoretischer Screening- und experimenteller Forschung zu thermoelektrischen Materialien.