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

Door het combineren van dynamische veldtheorie bij eindige temperatuur met DFT+U en berekeningen met de nudge-elastische-bandmethode, onthult deze studie dat dynamische correlaties de migratiebarrières voor Li-ionen in paramagnetisch Li$_2$MnO$_3$ aanzienlijk verlagen, waardoor een verenigde verklaring wordt geboden voor zowel de activeringsenergieën voor kort- als langafstandsvervoer zonder dat extrinsieke wanorde of geclusterde vacatures vereist zijn.

De kern

Stel je een lithium-ionbatterij voor als een drukke stad waar kleine lithiumionen de forenzen zijn die proberen van de ene naar de andere kant van de stad te komen om je telefoon of auto van stroom te voorzien. De "wegen" waarover ze reizen, bevinden zich in een materiaal genaamd Li₂MnO₃.

Lange tijd waren wetenschappers in de war over hoe snel deze forenzen konden bewegen. Sommige experimenten (die naar zeer korte afstanden keken) zeiden dat de wegen supersoepel en snel waren. Andere experimenten (die naar lange afstanden keken) zeiden dat de wegen vol stilstonden en zeer traag waren. Het was alsof je zegt: "Je kunt een sprint in 10 seconden lopen!" maar ook: "Je kunt geen marathon lopen omdat het parcours kapot is."

Dit artikel lost dat mysterie op door een supergeavanceerde computersimulatie te gebruiken om het "verkeer" op een nieuwe manier te bekijken.

De oude kaart versus de nieuwe kaart

Vroeger gebruikten wetenschappers een standaardcomputermodel (genaamd DFT+U) om de wegen in kaart te brengen. Dit model was als een basis-GPS: het zag de lithiumionen die probeerden over muren te springen, maar het berekende die muren als zeer hoog (ongeveer 0,6 tot 0,9 eV). Dit suggereerde dat deionen zeer langzaam zouden bewegen, wat niet overeenkwam met de snelle "sprint"-data van de experimenten op korte afstand.

De auteurs realiseerden zich dat het oude model een cruciaal ingrediënt miste: hitte en chaos. In de echte wereld zijn de atomen in de batterij niet op hun plaats bevroren; ze trillen en wiebelen door de hitte (temperatuur). De mangaanatomen in het materiaal hebben ook kleine magnetische spins die willekeurig omklappen. Het oude model behandelde deze spins alsof ze in een perfecte lijn bevroren waren, wat niet waar is voor een werkende batterij.

De "dynamische" simulatie

Om dit op te lossen, gebruikten de auteurs een krachtigere tool genaamd DFT+DMFT. Denk hierbij aan een upgrade van een statische 2D-kaart naar een 3D-simulatie in real-time die rekening houdt met de hitte en het willekeurige omklappen van magnetische spins.

Ze simuleerden een enkele "lege stoel" (een vacature) in de lithiumstad. De lithiumionen moeten in deze lege stoel springen om vooruit te komen.

De twee snelheden van reizen

Toen ze hun nieuwe, "hete en chaotische" simulatie uitvoerden, vonden ze iets verbazends. De energiebarrières (de muren die deionen moeten beklimmen) daalden aanzienlijk, maar alleen voor specifieke soorten sprongen.

1. De korte sprong (de sprint):
   Voor de allerkortste sprong tussen twee aangrenzende plekken toonde de nieuwe simulatie dat de muur slechts 0,18 eV hoog was.

   • Het resultaat: Dit komt perfect overeen met de "snelle sprint"-data van de experimenten op korte afstand.
   • De analogie: Stel je een foren voor die over een kleine stoeprand stapt. Het is makkelijk en snel. Het oude model dacht dat de stoeprand een 3 meter hoge omheining was; het nieuwe model realiseerde zich dat het slechts een kleine stap was.

2. De lange tocht (de marathon):
   Echter, om een lange afstand over de hele stad te reizen, kan de foren niet eeuwig de makkelijke stappen nemen. Uiteindelijk moet hij een iets moeilijkere stap zetten. De simulatie vond een tweede, iets hogere muur op 0,50 eV.

   • Het resultaat: Dit komt overeen met de "trage marathon"-data van de experimenten op lange afstand.
   • De analogie: Om de stad over te steken moet je veel makkelijke stappen nemen, maar af en toe loop je tegen een heuvel aan. Zelfs als de meeste stappen makkelijk zijn, wordt je totale snelheid beperkt door die ene heuvel.

Waarom dit belangrijk is

De grote ontdekking is dat je geen ingewikkelde verklaringen hoeft te verzinnen om het snelheidsprobleem op te lossen. Je hoeft niet aan te nemen dat de batterij vol zit met "klonten" lege stoelen of dat het materiaal kapot is.

Het artikel toont aan dat Li₂MnO₃ eigenlijk een zeer goed materiaal is (bijna perfect, of "stoichiometrisch"). De reden dat we verschillende snelheden zien in verschillende experimenten, is simpelweg omdat:

• Experimenten op korte afstand alleen de makkelijke, lage heuvels zien (0,18 eV).
• Experimenten op lange afstand de hele reis zien, die wordt vertraagd door de af en toe hogere heuvel (0,50 eV).

De conclusie

Door rekening te houden met de hitte en de magnetische "trilling" van de atomen, creëerden de auteurs één enkel, verenigd verhaal. Ze bewezen dat de lithiumionen lokaal gemakkelijk kunnen razen, maar dat hun totale reis wordt beheerst door een paar iets zwaardere stappen. Dit verklaart waarom de batterij zich anders gedraagt afhankelijk van hoe je het meet, zonder dat je defecten of onzuiverheden in het materiaal hoeft aan te wijzen.

Kortom: de batterij is niet kapot; we hadden gewoon een betere kaart nodig die rekening hield met de hitte en de magnetische dans van de atomen om te begrijpen hoe de lithiumionen echt bewegen.

Technische samenvatting

Technische Samenvatting: Het onthullen van korte- en langeafstands-Li-ionendiffusie in Li2MnO3 vanuit dynamische middelveldtheorie bij eindige temperatuur

Probleemstelling
Li2MnO3 is een kritisch onderdeel van Li-overschot gelaagde kathoden, maar de rol ervan bij het beperken van Li-iontransport blijft omstreden. Experimentele schattingen van Li-ionmigratie in Li2MnO3 vertonen een aanzienlijke discrepantie afhankelijk van de schaal van de sonde. Microscopische $\mu^+$SR-metingen op bijna stoichiometrische poeders rapporteren een lage activeringsenergie voor korte afstand ($E_a \approx 0,156$ eV) in het paramagnetische regime. Daarentegen leveren macroscopische wisselstroom-impedantiemetingen een aanzienlijk hogere schijnbare barrière ($E_a \approx 0,46$ eV) op voor transport over lange afstand.

Eerdere studies op basis van eerste principes, gebaseerd op Dichtheidsfunctionaaltheorie met Hubbard U-correcties (DFT+U), hebben moeite gehad om deze waarden met elkaar in overeenstemming te brengen. Standaard DFT+U-berekeningen voor een enkele Li-vacature voorspellen doorgaans barrières van 0,5–0,8 eV, wat te hoog is om de experimentele gegevens voor korte afstand te verklaren. Om de lage experimentele barrière te matchen, hebben eerdere modellen vaak gegroepeerde vacatureconfiguraties (divacatures of trivacatures) ingeroepen of specifieke magnetische ordeningen aangenomen. Dergelijke vacatureclusters zijn echter onwaarschijnlijk in bijna stoichiometrische monsters, en Li2MnO3 is paramagnetisch onder batterijwerkingsomstandigheden, terwijl veel DFT+U-studies ferromagnetische of niet-magnetische toestanden aannamen die lokale Mn-momenten onderdrukken.

Methodologie
Dit werk hanteert een multi-schaal computationele aanpak om Li+-migratie in paramagnetisch Li2MnO3 met een enkele Li-vacature bij kamertemperatuur (300 K) te bestuderen. De methodologie integreert drie componenten:

1. DFT+U: Gebruikt om migratiepaden en zadel-puntgeometrieën te optimaliseren via de Nudged-Elastic-Band (NEB)-methode. Spin-polariseerde DFT+U-berekeningen werden uitgevoerd met de PBEsol-functie en een 1×2×2 supercel.
2. DFT+DMFT: Dynamische Middelveldtheorie (DMFT) bij eindige temperatuur met een continuous-time quantum Monte Carlo (CTQMC) impureitoplosser werd toegepast om totale energieën te evalueren langs de vaste NEB-geometrieën. Dit kader vangt dynamische elektronische correlaties in de paramagnetische fase op.
3. Behandeling van de elektronische structuur: Om de lage kristallografische symmetrie (C2/m) van Li2MnO3 aan te pakken, die grote niet-diagonale matrixelementen in standaard Wannier-projecties veroorzaakt, diagonaliseerden de auteurs het Mn-d-blok. Dit creëerde een effectief laag-energetisch "d-only"-model waarbij de Wannier-basis impliciet Mn-d–O-p-hybridisatie encodeert, waardoor een nauwkeurige reproductie van de isolerende gap mogelijk wordt.

De studie evalueert zes symmetrie-inequivalente migratiepaden (vier intralaag en twee interlaag) door DMFT totale energieën te berekenen voor de initiële, zadel-punt- en finale beelden, waarbij de gecorreleerde elektronische toestand op elk punt mag relaxeren.

Belangrijkste resultaten
De toepassing van dynamische correlaties bij eindige temperatuur renormaliseert de activeringsenergieën voor specifieke migratiepaden aanzienlijk, waardoor de discrepantie tussen microscopische en macroscopische experimentele gegevens wordt opgelost zonder gegroepeerde vacatures in te roepen:

• Korte-afstandsdiffusie: Voor de hop met de laagste barrière in de intralaag (Pad ii: $4h \to T_{2b}^{Li} \to 4h$) wordt de DFT+DMFT activeringsenergie gereduceerd tot 0,18 eV. Deze waarde reproduceert kwantitatief de activeringsenergie voor korte afstand ($0,156$ eV) die is waargenomen in $\mu^+$SR-metingen. De reductie wordt voornamelijk gedreven door de eigenwaarde-correctieterm in de DMFT totale energie, die de zadel-puntconfiguratie stabiliseert.
• Lange-afstandsdiffusie: Het pad met de op één na laagste barrière (Pad iv: $4h \to T_{2b}^{Li} \to 2c$) heeft een activeringsenergie van 0,50 eV onder DMFT. Deze barrière controleert de percolatie van Li-ionen over lange afstand. Het is consistent met de macroscopische activeringsenergie ($\approx 0,46$ eV) die is ontleend aan wisselstroom-impedantiemetingen.
• Pad-afhankelijkheid: De barrièrereductie is sterk pad-afhankelijk. Waar Pad (ii) een dramatische daling ziet van de DFT+U-waarde (~0,68 eV) naar 0,18 eV, tonen andere paden (bijvoorbeeld Pad iii) minimale verandering. De auteurs correleren dit met de lokale registratie van de migrerende Li-ion ten opzichte van het MnO6-netwerk; Pad (ii) passeert door een meer open gebied met langere Li-Mn-afstanden, wat effectiever wordt gestabiliseerd door dynamische correlaties.
• Elektronische structuur: In tegenstelling tot DFT+U, dat een zwakke, site-afhankelijke verandering in Mn-d-bezetting toont bij het creëren van een vacature, onthult DMFT een sterke site-afhankelijke herverdeling van spectrale gewicht, met aanzienlijke uitputting van bezet Mn-d-gewicht op Mn-sites die verder van de vacature verwijderd zijn.

Betekenis en claims
Het artikel claimt een verenigde microscopische interpretatie te bieden van zowel lokaal als macroscopisch Li-iontransport in Li2MnO3. Door gebruik te maken van een beschrijving met een enkele vacature en paramagnetische DMFT, tonen de auteurs aan dat de hiërarchie van activeringsenergieën (laag voor korte afstand, hoger voor lange afstand) natuurlijk voortvloeit uit het samenspel van specifieke migratiepaden en dynamische correlaties bij eindige temperatuur.

De studie concludeert dat:

1. De experimenteel waargenomen activeringsenergie voor korte afstand kan worden gerechtvaardigd door migratie van een enkele vacature in een paramagnetische gastheer, waardoor de noodzaak om extrinsieke wanorde of gegroepeerde vacatureconfiguraties in te roepen om de gegevens te verklaren, komt te vervallen.
2. Lange-afstandstransport wordt gecontroleerd door een stap met een hogere barrière (0,50 eV) binnen een netwerk van hops, wat consistent is met macroscopische impedantiedata.
3. Dynamische correlaties bij eindige temperatuur een essentieel ingrediënt zijn voor het voorspellen van energetica van ionenmigratie in gecorreleerde oxide-elektroden, aangezien statische DFT+U-benaderingen falen in het vangen van de nodige renormalisatie van barrièrehoogten.

De auteurs merken op dat deze bevindingen de noodzaak benadrukken van behandelingen die verder gaan dan DFT voor ionendiffusie in gecorreleerde oxiden, en suggereren dat toekomstig werk dit kader moet uitbreiden naar hogere vacatureconcentraties en gemengde-kation Li-rijke gelaagde kathoden, met name waar zuurstofredox actief wordt.