Anisotropic Phonon Heat Flow and Thermoelectric Response in Tetragonal GeS$_2$ and GeSe$_2$

Resumen Técnico: Flujo de Calor de Fonones Anisotrópico y Respuesta Termoeléctrica en GeS₂ y GeSe₂ Tetragónicos

Planteamiento del Problema
Los materiales termoeléctricos, que convierten el calor directamente en electricidad, requieren un alto factor de mérito adimensional ($zT$), definido como la relación entre el factor de potencia ($\sigma S^2$) y la conductividad térmica total ($\kappa_e + \kappa_l$). Un desafío central en el diseño termoeléctrico es el acoplamiento mutuo de estos parámetros: aumentar la conductividad eléctrica ($\sigma$) a menudo incrementa la conductividad térmica electrónica ($\kappa_e$). Si bien la reducción de la conductividad térmica de red ($\kappa_l$) mediante la nanoestructuración o la aleación es una estrategia común, existe la necesidad de explorar compuestos no tóxicos y abundantes en la tierra con una $\kappa_l$ intrínsecamente baja y propiedades electrónicas favorables. Especicamente, las propiedades electrónicas, ópticas y de transporte de fonones de los dicalcogenuros de germanio (GeS₂ y GeSe₂) tetragonales en bulto permanecen poco exploradas desde una perspectiva termoeléctrica, a pesar de su anisotropía estructural y su potencial para el desacoplamiento del transporte electrónico y de fonones.

Metodología
Los autores emplearon un marco integral de primeros principios para investigar las propiedades estructurales, electrónicas y térmicas de GeS₂ y GeSe₂ tetragonales (grupo espacial $P42/nmc$).

• Estructura Electrónica: Los cálculos se realizaron utilizando la Teoría de la Funcional de la Densidad (DFT) con el funcional PBE y correcciones de dispersión de Grimme D2 para la relajación estructural y las propiedades elásticas. Para obtener brechas de banda precisas, se realizaron cálculos con el funcional híbrido con pantalla HSE03. Los coeficientes de transporte electrónico se calcularon mediante AMSET (Transporte de Boltzmann ab initio), incorporando mecanismos de dispersión dependientes del estado, incluyendo potenciales de deformación acústica, fonones ópticos polares e impurezas ionizadas.
• Dinámica de Red y Transporte Térmico: La dinámica de red armónica y anharmónica se investigó mediante la Teoría de Perturbación de la Funcional de la Densidad (DFPT) y métodos de desplazamiento finito. La conductividad térmica de red ($\kappa_l$) se calculó utilizando el código ShengBTE, resolviendo la ecuación de transporte de Boltzmann de fonones dentro de la Aproximación del Tiempo de Relajación (RTA).
• Análisis de Enlace: El origen microscópico del comportamiento de transporte se analizó mediante LOBSTER para evaluar los descriptores de la Población de Hamilton de Orbitales de Cristal Integrada (ICOHP) y la Población de Solapamiento de Orbitales de Cristal (ICOOP), caracterizando la naturaleza covalente-iónica de los enlaces Ge–X.
• Evaluación Termoeléctrica: El factor de mérito ($zT$) se calculó combinando los tensores de transporte electrónico de AMSET con los tensores de conductividad térmica de red de ShengBTE RTA, asegurando la consistencia componente a componente (emparejando los coeficientes electrónicos en el plano con la $\kappa_l$ en el plano y los coeficientes fuera del plano con la $\kappa_l$ fuera del plano).

Contribuciones Clave y Resultados

1. Estabilidad Estructural y Mecánica: Se confirmó que tanto el GeS₂ como el GeSe₂ tetragonales son mecánicamente y dinámicamente estables. Las estructuras exhiben una fuerte anisotropía, caracterizada por constantes elásticas en el plano ($C_{11}$) grandes y constantes fuera del plano ($C_{33}$) significativamente menores. Notablemente, ambos materiales presentan canales de deformación de cizalladura blandos (bajo $C_{66}$), lo cual es favorable para reducir la conductividad térmica de red.
2. Propiedades Electrónicas: La interpolación HSE03/Wannier produjo brechas de banda semiconductoras de 2.48 eV para GeS₂ y 1.23 eV para GeSe₂. El análisis LOBSTER reveló un enlace mixto covalente-iónico, donde los enlaces Ge–S muestran una magnitud de ICOHP estabilizadora mayor (−5.27 eV/enlace) en comparación con los enlaces Ge–Se (−4.74 eV/enlace).
3. Conductividad Térmica de Red Anisotrópica: El estudio reveló una pronunciada anisotropía en el transporte de calor. A 300 K, la conductividad térmica de red en el plano ($\kappa_{ab}$) es de 26.86 W m⁻¹ K⁻¹ para GeS₂ y 18.74 W m⁻¹ K⁻¹ para GeSe₂. En contraste, la conductividad fuera del plano ($\kappa_c$) es drásticamente menor, situándose en 1.19 W m⁻¹ K⁻¹ para GeS₂ y 1.52 W m⁻¹ K⁻¹ para GeSe₂. Esta supresión del transporte fuera del plano se atribuye al marco tetraédrico en capas, la reducción de las velocidades de grupo de fonones fuera del plano y la moderada anharmonicidad. La sustitución de S por Se reduce la frecuencia de fonón superior de ~13.6 THz a ~10.3 THz.
4. Desempeño Termoeléctrico: Combinando los datos de transporte electrónico y de red, los valores de $zT$ calculados son moderados. El mayor desempeño se observa para el GeS₂ tipo n fuera del plano a 800 K, produciendo un $zT_c$ de 0.257 con una concentración de portadores de $1 \times 10^{19}$ cm⁻³. Para GeSe₂, la estimación PBE-AMSET da un $zT_c$ de 0.066 para el transporte fuera del plano tipo p a 800 K. Los autores señalan que los resultados de GeSe₂ deben tratarse como cualitativos debido a la brecha de banda casi nula en la descripción PBE, lo que puede requerir una corrección de la brecha para lograr precisión cuantitativa.

Significancia y Reivindicaciones
El artículo afirma que los compuestos GeX₂ tetragonales representan una clase de materiales termoeléctricos fuertemente anisotrópicos donde la respuesta fuera del plano se beneficia significativamente de la supresión del transporte de calor de red. El estudio demuestra que, si bien una $\kappa_l$ baja es una condición necesaria, no es suficiente para un $zT$ alto; el desempeño está controlado finalmente por el equilibrio simultáneo del factor de potencia, la degradación de la movilidad a altas concentraciones de portadores y el transporte de calor electrónico. Los autores concluyen que los GeS₂ y GeSe₂ tetragonales prístinos exhiben valores de $zT$ calculados moderados, siendo el canal de transporte más prometedor a lo largo del eje $c$. Sugieren que futuras mejoras podrían lograrse mediante cálculos de dispersión con corrección de brecha, ingeniería de defectos, aleación y control de deformación para aprovechar aún más la favorable respuesta térmica fuera del plano.