Revealing Short- and Long-range Li-ion diffusion in Li$_2$MnO$_3$ from finite-temperature dynamical mean field theory

En combinant la théorie de champ moyen dynamique à température finie avec les calculs DFT+$U$ et les calculs de bande élastique poussée, cette étude révèle que les corrélations dynamiques abaissent considérablement les barrières de migration des ions Li dans Li$_2$MnO$_3$ paramagnétique, offrant une explication unifiée des énergies d'activation du transport à courte et à longue portée sans nécessiter de désordre extrinsèque ni de lacunes regroupées.

L'essentiel

Imaginez une batterie lithium-ion comme une ville animée où de minuscules ions lithium sont des navetteurs tentant de traverser la ville d'un bout à l'autre pour alimenter votre téléphone ou votre voiture. Les « routes » qu'ils empruntent se trouvent à l'intérieur d'un matériau appelé Li₂MnO₃.

Pendant longtemps, les scientifiques ont été perplexes quant à la vitesse à laquelle ces navetteurs pouvaient se déplacer. Certaines expériences (observant des distances très courtes) affirmaient que les routes étaient super lisses et rapides. D'autres expériences (observant des distances longues) indiquaient que les routes étaient encombrées de embouteillages et très lentes. C'était comme dire : « Vous pouvez faire un sprint en 10 secondes ! » mais aussi : « Vous ne pouvez pas courir un marathon car la piste est brisée. »

Cet article résout ce mystère en utilisant une simulation informatique ultra-avancée pour observer le « trafic » d'une nouvelle manière.

L'ancienne carte contre la nouvelle carte

Auparavant, les scientifiques utilisaient un modèle informatique standard (appelé DFT+U) pour cartographier les routes. Ce modèle était comme un GPS basique : il voyait les ions lithium tentant de sauter par-dessus des murs, mais il calculait ces murs comme étant très hauts (environ 0,6 à 0,9 eV). Cela suggérait que les ions se déplaceraient très lentement, ce qui ne correspondait pas aux données de « sprint rapide » issues des expériences à courte distance.

Les auteurs ont réalisé que l'ancien modèle manquait un ingrédient crucial : la chaleur et le chaos. Dans le monde réel, les atomes de la batterie ne sont pas figés sur place ; ils tremblotent et vibrent à cause de la chaleur (température). Les atomes de manganèse dans le matériau possèdent également de minuscules spins magnétiques qui basculent de manière aléatoire. L'ancien modèle traitait ces spins comme s'ils étaient figés dans une ligne parfaite, ce qui n'est pas le cas pour une batterie en fonctionnement.

La simulation « dynamique »

Pour corriger cela, les auteurs ont utilisé un outil plus puissant appelé DFT+DMFT. Imaginez cela comme une mise à niveau d'une carte statique en 2D vers une simulation 3D en temps réel qui prend en compte la chaleur et le basculement aléatoire des spins magnétiques.

Ils ont simulé un seul « siège vide » (une lacune) dans la ville lithium. Les ions lithium doivent sauter dans ce siège vide pour avancer.

Les deux vitesses de déplacement

Lorsqu'ils ont lancé leur nouvelle simulation « chaude et chaotique », ils ont découvert quelque chose d'extraordinaire. Les barrières énergétiques (les murs que les ions doivent gravir) ont chuté de manière significative, mais uniquement pour certains types de sauts.

1. Le petit bond (le sprint) :
   Pour le saut le plus court possible entre deux places voisines, la nouvelle simulation a montré que le mur n'avait qu'une hauteur de 0,18 eV.

   • Le résultat : Cela correspond parfaitement aux données de « sprint rapide » des expériences à courte distance.
   • L'analogie : Imaginez un navetteur franchissant un petit trottoir. C'est facile et rapide. L'ancien modèle pensait que le trottoir était une clôture de 3 mètres de haut ; le nouveau modèle a réalisé qu'il ne s'agissait que d'un petit pas.

2. La longue distance (le marathon) :
   Cependant, pour parcourir une longue distance à travers toute la ville, le navetteur ne peut pas se contenter de prendre les pas faciles indéfiniment. Il doit éventuellement effectuer un pas légèrement plus difficile. La simulation a révélé une deuxième barrière, légèrement plus haute, à 0,50 eV.

   • Le résultat : Cela correspond aux données de « marathon lent » des expériences à longue distance.
   • L'analogie : Pour traverser la ville, vous devez prendre de nombreux pas faciles, mais vous rencontrez occasionnellement une colline. Même si la plupart des pas sont faciles, votre vitesse globale est limitée par cette unique colline.

Pourquoi cela compte

La grande découverte est que vous n'avez pas besoin d'inventer des explications compliquées pour résoudre le problème de vitesse. Vous n'avez pas besoin de supposer que la batterie est remplie de « grappes » de sièges vides ou que le matériau est brisé.

L'article montre que Li₂MnO₃ est en réalité un matériau très performant (presque parfait, ou « stœchiométrique »). La raison pour laquelle nous observons des vitesses différentes dans différentes expériences est simplement que :

• Les expériences à courte portée ne voient que les petites collines faciles (0,18 eV).
• Les expériences à longue portée voient l'ensemble du trajet, qui est ralenti par la colline occasionnellement plus haute (0,50 eV).

L'essentiel

En tenant compte de la chaleur et du « tremblement » magnétique des atomes, les auteurs ont créé une histoire unique et unifiée. Ils ont prouvé que les ions lithium peuvent se déplacer facilement à l'échelle locale, mais que leur voyage global est contrôlé par quelques étapes légèrement plus difficiles. Cela explique pourquoi la batterie se comporte différemment selon la méthode de mesure, sans qu'il soit nécessaire de blâmer des défauts ou des impuretés dans le matériau.

En résumé : la batterie n'est pas brisée ; nous avions simplement besoin d'une meilleure carte qui prenait en compte la chaleur et la danse magnétique des atomes pour comprendre comment les ions lithium se déplacent réellement.

Résumé technique

Résumé technique : Révélation de la diffusion à courte et longue portée des ions Li dans Li2MnO3 à partir de la théorie de la fonctionnelle de champ moyen dynamique à température finie

Énoncé du problème
Li2MnO3 est un composant critique des cathodes en couches riches en Li, mais son rôle dans la limitation du transport des ions Li reste débattu. Les estimations expérimentales de la migration des ions Li dans Li2MnO3 présentent une divergence significative selon l'échelle de longueur de la sonde. Les mesures microscopiques par résonance de spin muonique ($\mu^+$SR) sur des poudres quasi-stœchiométriques rapportent une faible énergie d'activation à courte portée ($E_a \approx 0,156$ eV) dans le régime paramagnétique. En revanche, les mesures macroscopiques d'impédance AC donnent une barrière apparente nettement plus élevée ($E_a \approx 0,46$ eV) pour le transport à longue portée.

Les études antérieures basées sur la théorie de la fonctionnelle de la densité avec corrections de Hubbard U (DFT+U) ont eu du mal à réconcilier ces valeurs. Les calculs DFT+U standards pour une seule lacune de Li prédisent généralement des barrières de 0,5 à 0,8 eV, trop élevées pour expliquer les données expérimentales à courte portée. Pour correspondre à la faible barrière expérimentale, les modèles précédents invoquaient souvent des configurations de lacunes groupées (divacances ou trivacances) ou supposaient des ordres magnétiques spécifiques. Cependant, de tels amas de lacunes sont peu probables dans des échantillons quasi-stœchiométriques, et Li2MnO3 est paramagnétique dans les conditions de fonctionnement des batteries, alors que de nombreuses études DFT+U supposaient des états ferromagnétiques ou non magnétiques qui suppriment les moments locaux du Mn.

Méthodologie
Ce travail emploie une approche computationnelle multi-échelle pour étudier la migration des Li+ dans un Li2MnO3 paramagnétique avec une seule lacune de Li à température ambiante (300 K). La méthodologie intègre trois composantes :

1. DFT+U : Utilisée pour optimiser les voies de migration et les géométries des points de selle via la méthode de la bande élastique nudgée (NEB). Des calculs DFT+U polarisés en spin ont été effectués en utilisant le fonctionnel PBEsol et une supercellule 1×2×2.
2. DFT+DMFT : La théorie de la fonctionnelle de champ moyen dynamique (DMFT) à température finie, avec un solveur de défaut par Monte Carlo quantique en temps continu (CTQMC), a été appliquée pour évaluer les énergies totales le long des géométries NEB fixes. Ce cadre capture les corrélations électroniques dynamiques dans la phase paramagnétique.
3. Traitement de la structure électronique : Pour faire face à la faible symétrie cristalline (C2/m) de Li2MnO3, qui provoque de grands éléments de matrice hors diagonale dans les projections de Wannier standards, les auteurs ont diagonalisé le bloc d du Mn. Cela a créé un modèle effectif de basse énergie « d-seulement » où la base de Wannier encode implicitement l'hybridation Mn-d–O-p, permettant une reproduction précise du gap isolant.

L'étude évalue six voies de migration non équivalentes par symétrie (quatre intracouches et deux intercouches) en calculant les énergies totales DMFT pour les images initiale, de point de selle et finale, permettant à l'état électronique corrélé de se relaxer à chaque point.

Résultats clés
L'application des corrélations dynamiques à température finie renormalise significativement les énergies d'activation pour des voies de migration spécifiques, résolvant la divergence entre les données expérimentales microscopiques et macroscopiques sans invoquer de lacunes groupées :

• Diffusion à courte portée : Pour le saut intracouche à la barrière la plus basse (Voie ii : $4h \to T_{2b}^{Li} \to 4h$), l'énergie d'activation DFT+DMFT est réduite à 0,18 eV. Cette valeur reproduit quantitativement l'énergie d'activation à courte portée (0,156 eV) observée dans les mesures $\mu^+$SR. La réduction est principalement due au terme de correction par valeur propre dans l'énergie totale DMFT, qui stabilise la configuration du point de selle.
• Diffusion à longue portée : La voie suivante à la barrière la plus basse (Voie iv : $4h \to T_{2b}^{Li} \to 2c$) présente une énergie d'activation de 0,50 eV sous DMFT. Cette barrière contrôle la percolation à longue portée des ions Li. Elle est cohérente avec l'énergie d'activation macroscopique ($\approx 0,46$ eV) extraite des mesures d'impédance AC.
• Dépendance à la voie : La réduction de la barrière est fortement dépendante de la voie. Alors que la voie (ii) voit une chute dramatique de la valeur DFT+U (~0,68 eV) à 0,18 eV, d'autres voies (par exemple, la voie iii) montrent un changement minimal. Les auteurs corrèlent cela avec l'alignement local de l'ion Li migrant par rapport au réseau MnO6 ; la voie (ii) passe par une région plus ouverte avec des distances Li-Mn plus longues, qui est plus efficacement stabilisée par les corrélations dynamiques.
• Structure électronique : Contrairement à la DFT+U, qui montre un changement faiblement dépendant du site dans l'occupation des d du Mn lors de la création d'une lacune, la DMFT révèle une redistribution forte et dépendante du site du poids spectral, avec un appauvrissement significatif du poids occupé des d du Mn sur les sites Mn plus éloignés de la lacune.

Signification et affirmations
L'article prétend fournir une interprétation microscopique unifiée du transport des ions Li à la fois local et macroscopique dans Li2MnO3. En utilisant une description DMFT paramagnétique à une seule lacune, les auteurs démontrent que la hiérarchie des énergies d'activation (faible pour la courte portée, plus élevée pour la longue portée) découle naturellement de l'interaction entre des voies de migration spécifiques et des corrélations dynamiques à température finie.

L'étude conclut que :

1. L'énergie d'activation à courte portée observée expérimentalement peut être rationalisée par la migration d'une seule lacune dans un hôte paramagnétique, éliminant la nécessité d'invoquer un désordre extrinsèque ou des configurations de lacunes groupées pour expliquer les données.
2. Le transport à longue portée est contrôlé par une étape à barrière plus élevée (0,50 eV) au sein d'un réseau de sauts, cohérent avec les données d'impédance macroscopique.
3. Les corrélations dynamiques à température finie sont un ingrédient essentiel pour prédire l'énergétique de la migration ionique dans les électrodes d'oxydes corrélés, car les approches statiques DFT+U échouent à capturer la renormalisation nécessaire des hauteurs de barrière.

Les auteurs notent que ces résultats soulignent la nécessité de traitements au-delà de la DFT pour la diffusion ionique dans les oxydes corrélés et suggèrent que les travaux futurs devraient étendre ce cadre à des concentrations de lacunes plus élevées et à des cathodes en couches riches en Li à cations mixtes, en particulier là où le redox de l'oxygène devient actif.