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 del Artículo: 2311.06135
- Título: Propiedades de transporte termoeléctrico del sulfuro de hierro dopado con electrones (pirita FeS₂)
- Autores: Anustup Mukherjee, Alaska Subedi (CPHT, CNRS, École Polytechnique, Institut Polytechnique de Paris)
- Clasificación: cond-mat.mtrl-sci
- Fecha de Publicación: 13 de noviembre de 2023
- Enlace del Artículo: https://arxiv.org/abs/2311.06135
La pirita FeS₂ ha sido ampliamente estudiada para aplicaciones termoeléctricas debido al gran coeficiente de Seebeck observado a temperatura ambiente. Sin embargo, los valores de coeficiente de Seebeck reportados en la literatura son extremadamente sensibles a las propiedades de la muestra (natural o preparada en laboratorio, cristal masivo o película delgada), existiendo ambigüedad en la magnitud y signo del coeficiente de Seebeck del FeS₂ puro. También hay variación en la magnitud del coeficiente de Seebeck a temperatura ambiente en muestras dopadas con Co. Este trabajo utiliza cálculos de primeros principios para estudiar las propiedades termoeléctricas del FeS₂ dopado con Co, empleando tres esquemas de dopaje diferentes: sustitución explícita de Co, dopaje jellium y aproximación de banda rígida (RBA). Los coeficientes de Seebeck calculados para todas las concentraciones de dopaje con Co son menores que -50 μV/K, lo que sugiere que este sistema puede no ser práctico para aplicaciones termoeléctricas. Interesantemente, la RBA sobreestima significativamente el coeficiente de Seebeck porque no considera los efectos de las modificaciones estructurales inducidas por el dopaje y el apantallamiento de carga en la estructura electrónica.
- Inconsistencia en datos experimentales: Los valores medidos del coeficiente de Seebeck de FeS₂ muestran variaciones enormes, oscilando entre -115 y 524 μV/K, con signos inconsistentes
- Dependencia de la muestra: Las propiedades termoeléctricas dependen fuertemente del tipo de muestra (natural vs. artificial, masiva vs. película delgada)
- Precisión de métodos teóricos: La mayoría de investigaciones teóricas existentes emplean la aproximación de banda rígida (RBA), cuya precisión es cuestionable
- El FeS₂ está compuesto por elementos abundantes en la Tierra y no tóxicos, con potencial en aplicaciones energéticas como fotovoltaica, cátodos de baterías y termoelectricidad
- Aclarar sus propiedades termoeléctricas intrínsecas es crucial para evaluar el valor de aplicaciones prácticas
- Proporciona orientación metodológica para investigaciones teóricas de propiedades termoeléctricas en semiconductores dopados
- La aproximación de banda rígida (RBA) simula efectos de dopaje únicamente mediante el desplazamiento rígido del nivel de Fermi
- Ignora cambios estructurales inducidos por dopaje y efectos de apantallamiento de carga
- Carece de comparación sistemática entre diferentes esquemas de dopaje
- Comparación sistemática de tres esquemas de dopaje: Primera comparación de las diferencias entre sustitución explícita de Co, dopaje jellium y RBA en el cálculo de propiedades termoeléctricas de FeS₂
- Revelación de limitaciones de RBA: Descubrimiento de que RBA sobreestima significativamente el coeficiente de Seebeck, con explicación desde la perspectiva de estructura electrónica
- Aclaración del desempeño termoeléctrico del FeS₂ dopado con Co: Demostración de que el coeficiente de Seebeck es menor que 50 μV/K para todas las concentraciones de dopaje, indicando que este sistema no es adecuado para aplicaciones termoeléctricas
- Orientación metodológica: Proporciona información metodológica importante para investigaciones teóricas de propiedades termoeléctricas en semiconductores dopados
Se empleó teoría funcional de la densidad (DFT) dentro del marco de Vienna Ab initio Simulation Package (VASP):
- Funcional de intercambio-correlación: Aproximación de densidad local (LDA)
- Método de pseudopotencial: Método de ondas aumentadas proyectadas (PAW)
- Corte de energía: ≥420 eV
- Malla de puntos k: 8×8×8 para ciclos autoconsistentes, 24×24×24 para cálculo de densidad de estados
Simula electrones adicionales desplazando el nivel de Fermi en la estructura electrónica de FeS₂, correspondiente a un valor de dopaje de Co equivalente.
Añade explícitamente electrones en la celda unitaria de FeS₂ como carga de fondo, realizando cálculos de estructura electrónica.
Reemplaza progresivamente átomos de Fe con átomos de Co, considerando efectos de sustitución química real.
Se utilizó el código BoltzTraP2 para resolver la ecuación de transporte de Boltzmann linearizada bajo aproximación de tiempo de relajación constante (CRTA):
Coeficiente de Seebeck:
S(μ,T)=σαiναj
Coeficiente de Hall:
RH(μ,T)=σαjσiβναβk
donde los tensores de transporte se definen como:
σαβ(μ,T)=Ω1∫σαβ(ϵ)[−∂ϵfμ(ϵ)]dϵ
- Estructura cristalina: Pirita FeS₂, grupo espacial Pa3
- Composición de celda unitaria: 4 átomos de Fe (posición 4a), 8 átomos de S (posición 8c)
- Concentraciones de dopaje: Se investigaron dopajes de Co del 1%, 5%, 10%, 15%, 25% y 50%
Se consideraron simultáneamente fases no magnética y ferromagnética, ya que estudios experimentales indican que FexCo1-xS₂ es ferromagnético para x≥0.1.
- Las constantes de red de compuestos mixtos siguen la ley de Vegard
- Se relajaron posiciones atómicas para minimizar fuerzas atómicas
- Parámetro interno de FeS₂: u=0.3821
| Esquema de Dopaje | Fase No Magnética (μV/K) | Fase Ferromagnética (μV/K) |
|---|
| RBA | -46.71 | -46.71 |
| Jellium | -14.68 | -14.69 |
| Sustitución de Co | -19.67 | -20.49 |
| Esquema de Dopaje | Fase No Magnética (μV/K) | Fase Ferromagnética (μV/K) |
|---|
| RBA | -39.60 | -39.60 |
| Jellium | -17.52 | -17.52 |
| Sustitución de Co | -25.58 | -25.58 |
- Sobreestimación sistemática de RBA: RBA sobreestima significativamente el valor absoluto del coeficiente de Seebeck en todas las concentraciones de dopaje
- Consistencia en dependencia de temperatura: Los tres métodos muestran que el coeficiente de Seebeck aumenta con la temperatura
- Tipo de portador: El coeficiente de Seebeck es negativo en todos los casos, indicando que los electrones son los portadores principales
- Desempeño termoeléctrico limitado: El coeficiente de Seebeck a temperatura ambiente es menor que 50 μV/K para todas las concentraciones de dopaje
Mediante análisis de densidad de estados (DOS) y estructura de bandas se encontró:
- En RBA: El nivel de Fermi se ubica en una región de pendiente pronunciada de la DOS
- En dopaje jellium: El nivel de Fermi se desplaza a una región de hombro más ancho
- En sustitución de Co: El nivel de Fermi se ubica en una región relativamente plana
Dado que el coeficiente de Seebeck S ∝ d ln N(EF)/dE, la DOS más pronunciada cerca del nivel de Fermi en RBA resulta en valores de S mayores.
| Composición | 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 |
- Datos históricos: Coeficientes de Seebeck reportados en la primera mitad del siglo XX oscilan entre -115 y 524 μV/K
- Investigación moderna:
- Kato et al.: Muestras n tipo naturales con S = -300 μV/K
- Willeke et al.: Monocristales de alta movilidad con S alcanzando -320 μV/K
- Investigación de películas delgadas: Reportan universalmente valores positivos de S (60-75 μV/K)
- Cálculos tempranos: Gudelli et al. y Harran et al. utilizaron RBA prediciendo S optimizado superior a -400 μV/K
- Contribución de este trabajo: Primera comparación sistemática de diferentes esquemas de dopaje, revelando limitaciones de RBA
- Desempeño termoeléctrico limitado: El FeS₂ dopado con Co presenta coeficientes de Seebeck menores que 50 μV/K a temperatura ambiente para todas las concentraciones de dopaje estudiadas, muy por debajo del estándar para excelentes materiales termoeléctricos (>200 μV/K)
- Limitaciones del método RBA: La aproximación de banda rígida sobreestima significativamente el coeficiente de Seebeck porque ignora:
- Modificaciones estructurales inducidas por dopaje
- Efectos de apantallamiento de carga
- Cambios en ensanchamiento de banda y división relativa
- Significado metodológico: El FeS₂ con sustitución de Co es un ejemplo típico del fracaso de RBA, proporcionando una advertencia importante para investigaciones de propiedades termoeléctricas en semiconductores dopados
- Restricciones de métodos computacionales: El uso de LDA puede subestimar la brecha de banda, afectando la precisión de propiedades termoeléctricas
- Rango de concentración de dopaje: El esquema de sustitución explícita de Co está limitado por el tamaño de celda unitaria, imposibilitando estudiar concentraciones muy bajas
- Aproximación de tiempo de relajación: La aproximación de tiempo de relajación constante puede afectar la precisión de propiedades de transporte
- Mejora de métodos computacionales: Utilizar funcionales de intercambio-correlación más precisos (como HSE06)
- Extensión de sistemas de estudio: Investigar efectos de dopaje con otros metales de transición
- Validación experimental: Preparar muestras de alta calidad para verificar predicciones teóricas
- Innovación metodológica: Primera comparación sistemática de tres esquemas de dopaje, proporcionando un referente metodológico importante para el campo
- Perspectivas físicas profundas: Explicación detallada del mecanismo físico de sobreestimación de RBA desde la perspectiva de estructura electrónica
- Rigor computacional: Consideración de efectos magnéticos, uso de parámetros computacionales apropiados y criterios de convergencia
- Valor práctico claro: Indicación explícita de que el FeS₂ dopado con Co no es adecuado para aplicaciones termoeléctricas, evitando inversión experimental ineficaz
- Limitaciones de métodos teóricos: LDA puede no ser suficientemente preciso, particularmente para sistemas de electrones correlacionados
- Rango de temperatura limitado: Enfoque principal en propiedades a temperatura ambiente, carencia de investigación sistemática a altas temperaturas
- Mecanismo de dopaje único: Solo se investigó dopaje con electrones, sin considerar posibilidad de dopaje con huecos
- Contribución académica: Proporciona orientación metodológica importante para investigaciones teóricas de propiedades termoeléctricas en semiconductores dopados
- Valor práctico: Ayuda a investigadores a evitar desperdiciar recursos en sistemas no prometedores
- Reproducibilidad fuerte: Descripción suficiente de detalles computacionales, facilitando reproducción y extensión
- Investigación teórica: Proporciona referencia metodológica para investigaciones de propiedades termoeléctricas en otros semiconductores dopados
- Selección de materiales: Exclusión de candidatos inadecuados en cribado de alto rendimiento de materiales termoeléctricos
- Orientación experimental: Proporciona base teórica para que investigadores experimentales seleccionen estrategias de dopaje apropiadas
El artículo cita 38 referencias importantes que abarcan propiedades fundamentales de FeS₂, investigación experimental, métodos de cálculo teórico y otros aspectos, proporcionando una base bibliográfica sólida para la investigación. Particularmente notable es la citación de trabajos experimentales recientes de Xi et al., que parecen resolver controversias previas sobre el tipo de portador en FeS₂.
Evaluación General: Este es un artículo de investigación teórica de alta calidad que no solo aclara las propiedades termoeléctricas del FeS₂ dopado con Co, sino que más importantemente revela las limitaciones del método RBA ampliamente utilizado, proporcionando una contribución metodológica importante para investigaciones teóricas en este campo. Los resultados de investigación tienen valor importante para guiar la selección teórica y la investigación experimental de materiales termoeléctricos.