Microstructural Insights into Fast Ion Transport in Solid Electrolytes via Multiscale Modeling

Résumé technique : Aperçus microstructuraux du transport ionique rapide dans les électrolytes solides via la modélisation multi-échelle

Énoncé du problème
Le développement de batteries tout-solide à haute performance est entravé par une compréhension incomplète des caractéristiques microstructurales qui permettent un transport ionique rapide dans les électrolytes solides. Bien que les efforts historiques se soient concentrés sur l'optimisation de la conductivité superionique de masse (bulk), les électrolytes solides pratiques sont polycristallins, contenant des réseaux denses de joints de grains (GB - grain boundaries). Bien que les effets des GB dans les électrolytes d'oxydes soient partiellement compris, le rôle des GB dans les électrolytes de sulfure reste incertain en raison des difficultés expérimentales de caractérisation de la morphologie et de la séparation des contributions des GB par rapport à la conduction de masse. De plus, la modélisation du transport ionique dans les polycristaux est un défi computationnel : elle nécessite de lier les mécanismes de saut à l'échelle atomique à la percolation à l'échelle du continuum, exigeant des potentiels interatomiques capables de décrire avec précision les interactions complexes des GB et des cellules de simulation suffisamment grandes pour capturer le désordre structurel, ce que les méthodes ab initio traditionnelles ne peuvent atteindre aux échelles pertinentes.

Méthodologie
Les auteurs proposent un cadre de modélisation hiérarchique multi-échelle qui intègre la théorie de la fonctionnelle de la densité (DFT), les potentiels d'apprentissage automatique (MLP), la dynamique moléculaire (MD) et les simulations par éléments finis (FE).

1. Potentiels d'apprentissage automatique (MLP) : Pour équilibrer la précision et le coût computationnel, les auteurs ont développé des potentiels de tenseur de moment (MTP) légers en utilisant un schéma d'apprentissage actif en boucle fermée. Ce processus a consisté à extraire des environnements atomiques locaux à partir de microstructures complètes contenant des GB, à effectuer des calculs DFT sur ces structures locales, et à entraîner de manière itérative les MTP avec une complexité croissante jusqu'à atteindre une précision proche de l'ab initio.
2. Dynamique moléculaire (MD) : Les MLP entraînés ont été employés dans des simulations MD à grande échelle (jusqu'à la nanoseconde) de modèles de GB substituts, incluant des bicristaux, des cristaux colonnaires et des polycristaux avec des réseaux de GB aléatoires. Ces simulations ont capturé le désordre des ions Li, les corrélations et la cinétique de diffusion à des températures allant de 300 K à 600 K. Les métriques clés incluaient les auto-diffusivités, les barrières de diffusion et l'entropie excédentaire pour quantifier le désordre structurel.
3. Simulations par éléments finis (FE) : Les auto-diffusivités et les largeurs de GB dérivées de la MD ont été utilisées comme entrées pour des simulations FE à l'échelle du continuum. Pour gérer la séparation d'échelle, deux approches ont été utilisées :
   • GB résolus : Pour les nanograins (~25 nm), les GB ont été résolus par volume en utilisant des solveurs nodaux accélérés par Fourier (FANS).
   • GB effondrés : Pour les micrograins (~1 µm), les GB ont été représentés comme des éléments finis d'épaisseur effondrée.
     Cette approche a permis de prédire les diffusivités et conductivités macroscopiques pour des systèmes polycristallins avec des tailles de grains pertinentes pour les dispositifs.

Contributions clés et résultats
L'étude a été appliquée aux électrolytes solides argyrodites Li6PS5X (où X ∈ {Cl, Br, I}), en investiguant systématiquement les effets de l'ordonnancement des anions, de la substitution des halogénures et de la morphologie des GB.

• Énergétique et structure des GB : Pour le Li6PS5Cl avec des anions ordonnés, les énergies de formation des GB se sont révélées faibles (12–22 meV Å⁻² pour les GB généraux, ~5 meV Å⁻² pour Σ3), suggérant l'abondance de GB aléatoires dans les échantillons synthétisés. Les GB induisent un désordre structurel, comme en témoigne l'augmentation de l'entropie excédentaire et la présence de vides sub-nanométriques, bien que ces vides aient un impact minimal sur la diffusivité.
• Effets de l'ordonnancement des anions et de la substitution :
  • Désordre des anions : Dans le Li6PS5Cl à anions désordonnés, les GB agissent comme des barrières, réduisant la auto-diffusivité des ions Li d'environ 60 % par rapport à la masse. Cependant, les fréquences de saut caractéristiques restent élevées (~10⁸ Hz), ce qui est cohérent avec la faible résistance des GB observée expérimentalement.
  • Substitution des halogénures : Une relation d'échelle linéaire a été établie entre la barrière de diffusion des ions Li et le rayon de l'anion dans les phases ordonnées. L'augmentation du rayon de l'anion (Cl → Br → I) augmente la barrière de diffusion d'environ 80 meV par augmentation de 0,1 Å.
  • Effets de GB opposés : L'impact des GB dépend de la masse (bulk). Dans les phases à diffusion rapide (Li6PS5Cl à anions désordonnés), les GB suppriment la diffusion. Dans les phases à diffusion lente (Li6PS5I à anions ordonnés), les GB augmentent la diffusion jusqu'à 60 fois, agissant comme des canaux de connexion.
• Transport non-Arrhenius : Les simulations ont reproduit avec succès le comportement de transport non-Arrhenius observé dans le Li6PS5I à anions ordonnés, identifiant une température de transition (~250 K) où le mécanisme de transport passe d'un régime dominé par la masse (haute température) à un régime dominé par les GB (basse température). Ce comportement est attribué à l'énergie d'activation plus faible des GB par rapport à la masse dans ce système spécifique.
• Effets de taille de grain : Les simulations de continuum ont révélé que la réduction de la taille des grains à l'échelle nanométrique atténue le comportement non-Arrhenius. Pour le Li6PS5Cl, l'affinement des grains améliore le contact intergranulaire sans compromettre la conductivité superionique, tandis que pour le Li6PS5I, la nanotalisation peut activer le transport dans les systèmes dépourvus de comportement superionique intrinsèque.

Signification et affirmations
L'article affirme que son cadre de modélisation multi-échelle fournit un lien rigoureux et quantitatif entre les mécanismes à l'échelle atomique et les performances macroscopiques des dispositifs, résolvant ainsi le défi de longue date de la modélisation du transport ionique dans les électrolytes solides polycristallins.

• Aperçu mécanistique : L'étude clarifie que les GB n'entravent pas universellement le transport ionique ; plutôt, leur effet (amélioration vs suppression) est dicté par la diffusivité intrinsèque de la phase de masse.
• Directives de conception : Les résultats offrent des orientations a priori pour la conception microstructurale. Pour les sulfures à diffusion rapide comme le Li6PS5Cl, l'affinement des grains est bénéfique pour le contact sans sacrifier la conductivité. Pour les systèmes à diffusion lente comme le Li6PS5I, la nanotalisation peut être une stratégie pour activer le transport.
• Paramètres compositionnels : Le travail identifie le rayon de l'anion comme un paramètre compositionnel clé régissant les barrières de diffusion dans les argyrodites ordonnées, remettant en question l'idée que la polarisabilité du réseau est le seul facteur dominant.
• Validation expérimentale : Le comportement non-Arrhenius prédit et les températures de transition dans le Li6PS5I concordent avec les observations expérimentales, suggérant que les GB sont à l'origine de ce phénomène dans ces matériaux.

Les auteurs concluent que cette approche permet de concevoir des électrolytes solides avancés en tenant compte explicitement de la complexité microstructurale, allant au-delà des modèles de masse idéalisés vers des représentations polycristallines réalistes.