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.
- ID de l'article: 2311.06135
- Titre: Thermoelectric transport properties of electron doped pyrite FeS₂
- Auteurs: Anustup Mukherjee, Alaska Subedi (CPHT, CNRS, École Polytechnique, Institut Polytechnique de Paris)
- Classification: cond-mat.mtrl-sci
- Date de publication: 13 novembre 2023
- Lien de l'article: https://arxiv.org/abs/2311.06135
La pyrite FeS₂ a été largement étudiée pour des applications thermoélectriques en raison du grand coefficient de Seebeck observé à température ambiante. Cependant, les valeurs de coefficient de Seebeck rapportées dans la littérature sont extrêmement sensibles aux propriétés de l'échantillon (naturel ou préparé en laboratoire, cristal massif ou film mince), et il existe une ambiguïté quant à la magnitude et au signe du coefficient de Seebeck de FeS₂ pur. Des variations sont également observées dans la magnitude du coefficient de Seebeck à température ambiante pour les échantillons dopés au Co. Cet article utilise des calculs ab initio pour étudier les propriétés thermoélectriques de FeS₂ dopé au Co, en employant trois schémas de dopage différents : substitution explicite du Co, dopage jellium et approximation de bande rigide (RBA). Les coefficients de Seebeck calculés pour toutes les concentrations de dopage au Co sont inférieurs à -50 μV/K, suggérant que ce système pourrait ne pas être pratique pour les applications thermoélectriques. Fait intéressant, la RBA surestime considérablement le coefficient de Seebeck, car elle ne tient pas compte des modifications structurales induites par le dopage et des effets d'écrantage des charges sur la structure électronique.
- Incohérence des données expérimentales: Les valeurs mesurées du coefficient de Seebeck de FeS₂ présentent des variations énormes, allant de -115 à 524 μV/K, avec des signes également incohérents
- Dépendance aux propriétés de l'échantillon: Les propriétés thermoélectriques dépendent fortement du type d'échantillon (naturel vs artificiel, massif vs film mince)
- Précision des méthodes théoriques: Les études théoriques existantes utilisent principalement l'approximation de bande rigide (RBA), dont la précision est douteuse
- FeS₂ est composé d'éléments abondants sur terre et non toxiques, avec un potentiel pour les applications énergétiques telles que le photovoltaïque, les cathodes de batterie et la thermoélectricité
- Clarifier ses propriétés thermoélectriques intrinsèques est crucial pour évaluer la valeur d'application pratique
- Fournir des orientations méthodologiques pour la recherche théorique sur les propriétés thermoélectriques des semi-conducteurs dopés
- L'approximation de bande rigide (RBA) simule les effets de dopage uniquement par un déplacement rigide du niveau de Fermi
- Elle ignore les changements structuraux induits par le dopage et les effets d'écrantage des charges
- Manque de comparaison systématique entre différents schémas de dopage
- Comparaison systématique de trois schémas de dopage: Première comparaison des différences entre substitution explicite du Co, dopage jellium et RBA dans le calcul des propriétés thermoélectriques de FeS₂
- Révélation des limitations de la RBA: Découverte que la RBA surestime considérablement le coefficient de Seebeck, avec explication des raisons du point de vue de la structure électronique
- Clarification des performances thermoélectriques de FeS₂ dopé au Co: Démonstration que le coefficient de Seebeck pour toutes les concentrations de dopage est inférieur à 50 μV/K, indiquant que ce système ne convient pas aux applications thermoélectriques
- Fourniture de directives méthodologiques: Apport d'informations méthodologiques importantes pour la recherche théorique sur les propriétés thermoélectriques des semi-conducteurs dopés
Utilisation du cadre de la théorie de la fonctionnelle de la densité (DFT), avec le package Vienna Ab initio Simulation Package (VASP):
- Fonctionnelle d'échange-corrélation: Approximation de la densité locale (LDA)
- Méthode de pseudopotentiel: Méthode des ondes augmentées projetées (PAW)
- Énergie de coupure: ≥420 eV
- Grille de points k: 8×8×8 pour la boucle autocohérente, 24×24×24 pour le calcul de la densité d'états
Simulation des électrons supplémentaires par déplacement du niveau de Fermi dans la structure électronique de FeS₂, correspondant à une valeur de dopage au Co équivalente.
Ajout explicite d'électrons dans la maille unitaire de FeS₂ en tant que charge de fond, suivi d'un calcul de structure électronique.
Remplacement progressif des atomes de Fe par des atomes de Co, tenant compte des effets de substitution chimique réels.
Utilisation du code BoltzTraP2 pour résoudre l'équation de transport de Boltzmann linéarisée sous l'approximation du temps de relaxation constant (CRTA):
Coefficient de Seebeck:
S(μ,T)=σαiναj
Coefficient de Hall:
RH(μ,T)=σαjσiβναβk
où les tenseurs de transport sont définis comme:
σαβ(μ,T)=Ω1∫σαβ(ϵ)[−∂ϵfμ(ϵ)]dϵ
- Structure cristalline: Pyrite FeS₂, groupe d'espace Pa3
- Composition de la maille unitaire: 4 atomes de Fe (position 4a), 8 atomes de S (position 8c)
- Concentrations de dopage: Étude du dopage au Co à 1%, 5%, 10%, 15%, 25% et 50%
Prise en compte simultanée des phases non magnétiques et ferromagnétiques, car les études expérimentales montrent que FexCo₁₋ₓS₂ est ferromagnétique pour x≥0,1.
- Les constantes de réseau des composés mixtes suivent la loi de Vegard
- Relaxation des positions atomiques pour minimiser les forces atomiques
- Paramètre interne de FeS₂: u=0,3821
| Schéma de dopage | Phase non magnétique (μV/K) | Phase ferromagnétique (μV/K) |
|---|
| RBA | -46,71 | -46,71 |
| Jellium | -14,68 | -14,69 |
| Substitution Co | -19,67 | -20,49 |
| Schéma de dopage | Phase non magnétique (μV/K) | Phase ferromagnétique (μV/K) |
|---|
| RBA | -39,60 | -39,60 |
| Jellium | -17,52 | -17,52 |
| Substitution Co | -25,58 | -25,58 |
- Surestimation systématique de la RBA: La RBA surestime considérablement la valeur absolue du coefficient de Seebeck pour toutes les concentrations de dopage
- Cohérence de la dépendance à la température: Les trois méthodes montrent que le coefficient de Seebeck augmente avec la température
- Type de porteur de charge: Le coefficient de Seebeck est négatif dans tous les cas, indiquant que les électrons sont les porteurs de charge dominants
- Performance thermoélectrique limitée: Le coefficient de Seebeck à température ambiante est inférieur à 50 μV/K pour toutes les concentrations de dopage
L'analyse de la densité d'états (DOS) et de la structure de bandes révèle:
- Dans la RBA: Le niveau de Fermi est situé dans une région abrupte de la densité d'états
- Dans le dopage jellium: Le niveau de Fermi se déplace vers une région d'épaule plus large
- Dans la substitution Co: Le niveau de Fermi est situé dans une région relativement plate
Puisque le coefficient de Seebeck S ∝ d ln N(EF)/dE, une densité d'états plus abrupte près du niveau de Fermi dans la RBA conduit à des valeurs de S plus grandes.
| Composition | RBA (μV/K) | Jellium (μ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 |
- Données historiques: Les coefficients de Seebeck rapportés dans la première moitié du 20e siècle varient de -115 à 524 μV/K
- Recherches modernes:
- Kato et al.: Échantillons naturels de type n avec S = -300 μV/K
- Willeke et al.: Monocristaux à haute mobilité avec S atteignant -320 μV/K
- Études de films minces: Rapportent généralement des valeurs positives de S (60-75 μV/K)
- Calculs antérieurs: Gudelli et al. et Harran et al. ont utilisé la RBA pour prédire des S optimisés dépassant -400 μV/K
- Contribution de ce travail: Première comparaison systématique de différents schémas de dopage, révélant les limitations de la RBA
- Performance thermoélectrique limitée: FeS₂ dopé au Co présente un coefficient de Seebeck inférieur à 50 μV/K à température ambiante pour toutes les concentrations de dopage étudiées, bien en dessous des normes des excellents matériaux thermoélectriques (>200 μV/K)
- Limitations de la méthode RBA: L'approximation de bande rigide surestime considérablement le coefficient de Seebeck car elle ignore:
- Les modifications structurales induites par le dopage
- Les effets d'écrantage des charges
- Les changements dans l'élargissement des bandes et la séparation relative
- Signification méthodologique: FeS₂ dopé au Co est un exemple typique d'échec de la RBA, fournissant un avertissement important pour la recherche sur les propriétés thermoélectriques des semi-conducteurs dopés
- Limitations des méthodes de calcul: L'utilisation de la LDA peut sous-estimer la bande interdite, affectant la précision des propriétés thermoélectriques
- Plage de concentrations de dopage: Le schéma de substitution explicite du Co est limité par la taille de la maille unitaire, empêchant l'étude de très faibles concentrations de dopage
- Approximation du temps de relaxation: L'adoption de l'approximation du temps de relaxation constant peut affecter la précision des propriétés de transport
- Amélioration des méthodes de calcul: Utilisation de fonctionnelles d'échange-corrélation plus précises (comme HSE06)
- Extension des systèmes étudiés: Étude des effets du dopage par d'autres métaux de transition
- Vérification expérimentale: Préparation d'échantillons de haute qualité pour valider les prédictions théoriques
- Innovation méthodologique: Première comparaison systématique de trois schémas de dopage, fournissant un repère méthodologique important pour le domaine
- Perspectives physiques profondes: Explication détaillée du mécanisme physique de la surestimation par la RBA du point de vue de la structure électronique
- Calculs rigoureux: Prise en compte des effets magnétiques, utilisation de paramètres de calcul appropriés et de critères de convergence
- Valeur pratique claire: Indication explicite que FeS₂ dopé au Co ne convient pas aux applications thermoélectriques, évitant les investissements expérimentaux inefficaces
- Limitations des méthodes théoriques: La LDA peut ne pas être suffisamment précise, particulièrement pour les systèmes d'électrons corrélés
- Plage de température limitée: Accent principal sur les propriétés à température ambiante, manque d'étude systématique à haute température
- Mécanisme de dopage unique: Étude uniquement du dopage de type n, sans considération du dopage de type p
- Contribution académique: Fourniture de directives méthodologiques importantes pour la recherche théorique sur les propriétés thermoélectriques des semi-conducteurs dopés
- Valeur pratique: Aide les chercheurs à éviter de gaspiller des ressources sur des systèmes non prometteurs
- Forte reproductibilité: Description suffisante des détails de calcul, facilitant la reproduction et l'extension
- Recherche théorique: Fourniture de références méthodologiques pour l'étude des propriétés thermoélectriques d'autres semi-conducteurs dopés
- Sélection de matériaux: Exclusion de candidats inappropriés dans le criblage à haut débit des matériaux thermoélectriques
- Orientation expérimentale: Fourniture de bases théoriques aux chercheurs expérimentaux pour le choix de stratégies de dopage appropriées
L'article cite 38 références importantes, couvrant les propriétés fondamentales de FeS₂, la recherche expérimentale, les méthodes de calcul théorique et d'autres aspects, fournissant une base bibliographique solide pour la recherche. Particulièrement notable est la citation des travaux expérimentaux récents de Xi et al., qui semblent résoudre les désaccords antérieurs concernant le type de porteur de charge de FeS₂.
Évaluation générale: Cet article est une recherche théorique de haute qualité qui non seulement clarifie les propriétés thermoélectriques de FeS₂ dopé au Co, mais plus important encore, révèle les limitations de la méthode RBA largement utilisée, fournissant une contribution méthodologique importante à la recherche théorique du domaine. Les résultats de la recherche ont une valeur importante pour guider le criblage théorique et la recherche expérimentale des matériaux thermoélectriques.