Thermal Transport in Ag8TS6 (T= Si, Ge, Sn) Argyrodites: An Integrated Experimental, Quantum-Chemical, and Computational Modelling Study

Résumé technique : Transport thermique dans les argyrodites $\text{Ag}_8\text{TS}_6$ ($\text{T} = \text{Si, Ge, Sn}$)

Énoncé du problème
Les matériaux thermoélectriques nécessitent une faible conductivité thermique de réseau ($\kappa_L$) pour maximiser l'efficacité, pourtant de nombreux candidats performants souffrent d'une corrélation entre une faible $\kappa_L$ et une conductivité ionique élevée, ce qui peut compromettre la stabilité des dispositifs. Bien que les argyrodites à base d'argent (formule générale $A_{(11-n)/m}^{m+}\text{T}_n^+\text{Q}_6^{2-}$) soient connues pour leurs $\kappa_L$ exceptionnellement basses et leur conductivité ionique élevée, les mécanismes spécifiques régissant le transport thermique dans les variantes sulfurées $\text{Ag}_8\text{TS}_6$ ($\text{T} = \text{Si, Ge, Sn}$) restent incomplètement compris. Les modèles computationnels existants pour prédire la $\kappa_L$ font face à un compromis : les modèles traditionnels (par exemple, Slack, Cahill) manquent souvent de précision sur toute la plage de température ou surestiment les valeurs, tandis que les méthodes ab initio rigoureuses sont trop coûteuses en calcul pour le criblage à haut débit. De plus, l'interaction entre la liaison chimique, l'anharmonicité et l'indépendance de la conductivité thermique par rapport à la conductivité ionique dans ces sulfures spécifiques nécessite une investigation approfondie combinant des données expérimentales avec une modélisation avancée.

Méthodologie
L'étude emploie une approche intégrée combinant la synthèse et la caractérisation expérimentales avec une modélisation computationnelle multi-échelle :

1. Caractérisation expérimentale :

   • Synthèse : Des échantillons de $\text{Ag}_8\text{TS}_6$ ($\text{T} = \text{Si, Ge, Sn}$) ont été synthétisés via des réactions à l'état solide à haute température.
   • Analyse structurelle : La diffraction des rayons X sur poudre (XRD) a été réalisée de 100 K à 400 K pour confirmer la pureté de la phase et étudier d'éventuelles transitions de phase.
   • Mesures de transport : La diffusivité thermique a été mesurée par la méthode du flash laser (173–500 K), et la conductivité ionique a été déterminée par spectroscopie d'impédance électrochimique (233–303 K) en utilisant des électrodes bloquant les électrons pour isoler les contributions ioniques. La conductivité électronique a été vérifiée comme étant négligeable.
   • Propriétés élastiques : Des méthodes d'impulsion ultrasonore ont été utilisées pour mesurer les vitesses du son.

2. Modélisation computationnelle :

   • Chimie quantique : Des calculs de la théorie de la fonctionnelle de la densité (DFT) (VASP) ont été utilisés pour optimiser les structures et calculer les propriétés électroniques. L'analyse des liaisons a été réalisée via les populations de Hamiltoniens d'orbitales cristallines (COHP) et les indices de liaison d'orbitales cristallines (COBI) via le code LOBSTER pour quantifier la force des liaisons, la covalence et les interactions multicentriques.
   • Propriétés de phonons : Les structures de bandes de phonons harmoniques, les vitesses de groupe et les températures de Debye ont été calculées. Les paramètres de Grüneisen dépendants des modes ont été dérivés de l'approximation quasi-harmonique (QHA) pour quantifier l'anharmonicité.
   • Prédiction de la conductivité thermique : Deux approches de modélisation distinctes ont été développées et comparées :
     • Modèle à deux canaux basé sur Grüneisen : Une approche modifiée combinant le modèle à deux canaux de Xia (contributions phonon + diffuson) avec la méthode de Bjerg pour estimer les temps de vie des phonons basés sur les paramètres de Grüneisen et la théorie de Slack. Cela évite les calculs coûteux de constantes de force anharmoniques.
     • Potentiels interatomiques appris par apprentissage automatique (MLIP) : Un MLIP fondamental (MACE-MP-03b) a été utilisé pour calculer la conductivité thermique via l'approche complète de la dynamique de réseau à deux canaux (Simoncelli et al.), en résolvant l'équation de transport de Wigner.

Contributions clés et résultats

• Liaison et anharmonicité : L'analyse des liaisons chimiques révèle que si les liaisons $\text{T-S}$ sont fortes et covalentes, les interactions $\text{Ag-S}$ sont faibles, caractérisées par des états antibliants occupés sous le niveau de Fermi. Cela conduit à des environnements de coordination $\text{Ag}^+$ distordus et à un comportement de type « agitateur » (rattler-like). L'étude identifie une hétérogénéité de liaison significative, incluant de faibles interactions $\text{Ag-Ag}$ et des liaisons multicentriques hypervalentes ($\text{S-Ag-S}$), qui induisent une forte anharmonicité.
• Propriétés de phonons : Tous les trois composés présentent des structures dynamiquement stables sans modes de phonons imaginaires. Les modes acoustiques de basse fréquence sont dominés par les vibrations de $\text{Ag}^+$, entraînant de faibles vitesses du son (environ 1000–1500 m/s) et de faibles températures de Debye. De grands paramètres de Grüneisen sont observés pour ces modes de basse énergie, confirmant une forte anharmonicité.
• Modélisation de la conductivité thermique :
  • Le modèle basé sur Grüneisen et l'approche par MLIP prédisent tous deux avec succès la conductivité thermique de réseau ultra-basse des composés $\text{Ag}_8\text{TS}_6$, montrant un bon accord avec les données expérimentales, particulièrement au-dessus de 200 K.
  • Les résultats confirment que le transport de chaleur est dominé par le canal des diffusons (modes vibratoires non-propageants), ce qui est cohérent avec les découvertes précédentes sur les argyrodites de séléniure.
  • Les écarts à basse température (0–50 K) dans les données expérimentales par rapport aux modèles idéalisés ont été attribués à la diffusion par défaut ponctuel et à la diffusion aux joints de grains. Un modèle analytique incorporant ces mécanismes de diffusion a réussi à reproduire le comportement à basse température.
• Découplage des conductivités : Les mesures expérimentales démontrent que les conductivités thermique et ionique sont indépendantes dans ces matériaux. Alors que la conductivité ionique augmente de plus d'un ordre de grandeur entre 233 K et 303 K, la conductivité thermique reste presque constante. Cela suggère que les conductivités ionique et thermique peuvent être ajustées indépendamment, une découverte cruciale pour les applications thermoélectriques.
• Substitution isovalente : Le remplacement de l'atome de tétrel ($\text{Si, Ge, Sn}$) entraîne des différences minimales dans le caractère de liaison, les vitesses du son, les paramètres de Grüneisen et la conductivité thermique, indiquant que les propriétés de transport thermique sont robustes à travers cette famille.

Signification et affirmations
L'article affirme établir un cadre robuste pour prédire la conductivité thermique des matériaux inorganiques complexes en intégrant l'analyse des liaisons avec des modèles computationnels efficaces. La principale signification réside dans la démonstration que :

1. L'analyse des liaisons comme prédicteur : L'environnement de liaison spécifique (liaisons $\text{Ag-S}$ faibles, états antibliants, interactions multicentriques) est suffisant pour prédire la faible conductivité thermique et la forte anharmonicité de ces matériaux.
2. Criblage à haut débit efficace : Le modèle à deux canaux basé sur Grüneisen et l'approche MLIP fondamentale offrent des alternatives moins coûteuses en calcul que les calculs anharmoniques complets ab initio. Ces méthodes sont adaptées au criblage à haut débit de la conductivité thermique, à condition que les effets de diffusion (défauts, joints de grains) soient pris en compte dans l'interprétation des données à basse température.
3. Ajustement indépendant : L'étude renforce la conclusion selon laquelle, dans les argyrodites à base d'Ag, les conductivités thermique et ionique sont découplées, permettant l'optimisation potentielle des performances thermoélectriques sans compromettre la stabilité ionique.

Les auteurs notent modestement que bien que le MLIP et les modèles de Grüneisen montrent un fort accord avec les expériences, des vérifications supplémentaires sont nécessaires pour définir les limites des potentiels fondés sur le ML pour les systèmes complexes et pour affiner les procédures d'entraînement automatisées pour des espaces de composition spécifiques.