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

Durch die Kombination der dynamischen Mean-Field-Theorie bei endlichen Temperaturen mit DFT+$U$- und Nudged-Elastic-Band-Rechnungen zeigt diese Studie, dass dynamische Korrelationen die Li-Ionen-Migrationsbarrieren in paramagnetischem Li$_2$MnO$_3$ signifikant senken und damit eine einheitliche Erklärung sowohl für die Aktivierungsenergien des Kurz- als auch des Langstreckentransports liefern, ohne extrinsische Unordnung oder geclusterte Leerstellen zu benötigen.

Das Wesentliche

Stellen Sie sich eine Lithium-Ionen-Batterie als eine belebte Stadt vor, in der winzige Lithium-Ionen Pendler sind, die versuchen, von einer Seite der Stadt zur anderen zu gelangen, um Ihr Handy oder Ihr Auto mit Energie zu versorgen. Die „Straßen", auf denen sie reisen, befinden sich in einem Material namens Li₂MnO₃.

Lange Zeit waren Wissenschaftler verwirrt darüber, wie schnell diese Pendler sich bewegen können. Einige Experimente (die sehr kurze Distanzen betrachteten) sagten, die Straßen seien supersmooth und schnell. Andere Experimente (die lange Distanzen betrachteten) sagten, die Straßen seien voller Staus und sehr langsam. Es war, als würde man sagen: „Sie können einen Sprint in 10 Sekunden laufen!" aber auch: „Sie können keinen Marathon laufen, weil die Bahn kaputt ist."

Dieser Artikel löst dieses Rätsel, indem er eine hochmoderne Computersimulation verwendet, um den „Verkehr" auf eine neue Weise zu betrachten.

Die alte Karte vs. die neue Karte

Früher verwendeten Wissenschaftler ein Standard-Computermodell (genannt DFT+U), um die Straßen zu kartieren. Dieses Modell war wie ein einfaches GPS: Es sah, wie die Lithium-Ionen versuchten, über Mauern zu springen, berechnete diese Mauern jedoch als sehr hoch (etwa 0,6 bis 0,9 eV). Dies deutete darauf hin, dass sich die Ionen sehr langsam bewegen würden, was nicht mit den schnellen „Sprint"-Daten aus den Kurzstrecken-Experimenten übereinstimmte.

Die Autoren erkannten, dass das alte Modell eine entscheidende Zutat vermisste: Wärme und Chaos. In der realen Welt sind die Atome in der Batterie nicht an Ort und Stelle gefroren; sie wackeln und vibrieren aufgrund der Wärme (Temperatur). Die Manganatome im Material haben zudem winzige magnetische Spins, die zufällig hin und her kippen. Das alte Modell behandelte diese Spins so, als wären sie in einer perfekten Linie gefroren, was für eine funktionierende Batterie nicht zutrifft.

Die „dynamische" Simulation

Um dies zu beheben, verwendeten die Autoren ein leistungsfähigeres Werkzeug namens DFT+DMFT. Denken Sie daran wie an ein Upgrade von einer statischen 2D-Karte zu einer 3D-Echtzeitsimulation, die die Wärme und das zufällige Kippen der magnetischen Spins berücksichtigt.

Sie simulierten einen einzelnen „leeren Sitz" (eine Leerstelle) in der Lithium-Stadt. Die Lithium-Ionen müssen in diesen leeren Sitz springen, um vorwärts zu kommen.

Die zwei Geschwindigkeiten der Reise

Als sie ihre neue, „heiße und chaotische" Simulation durchführten, entdeckten sie etwas Erstaunliches. Die Energiebarrieren (die Mauern, die die Ionen erklimmen müssen) sanken erheblich, aber nur für bestimmte Arten von Sprüngen.

1. Der kurze Sprung (Der Sprint):
   Für den aller kürzesten Sprung zwischen zwei benachbarten Stellen zeigte die neue Simulation, dass die Mauer nur 0,18 eV hoch war.

   • Das Ergebnis: Dies stimmt perfekt mit den „schnellen Sprint"-Daten aus den Kurzstrecken-Experimenten überein.
   • Die Analogie: Stellen Sie sich einen Pendler vor, der über eine kleine Bordsteinkante steigt. Das ist einfach und schnell. Das alte Modell dachte, die Bordsteinkante sei ein 10-Fuß-Zaun; das neue Modell erkannte, dass es nur ein kleiner Schritt war.

2. Die lange Strecke (Der Marathon):
   Um jedoch eine lange Strecke quer durch die ganze Stadt zurückzulegen, kann der Pendler nicht für immer nur die leichten Schritte nehmen. Irgendwann muss er einen etwas härteren Schritt machen. Die Simulation fand eine zweite, etwas höhere Mauer bei 0,50 eV.

   • Das Ergebnis: Dies stimmt mit den „langsamen Marathon"-Daten aus den Langstrecken-Experimenten überein.
   • Die Analogie: Um quer durch die Stadt zu kommen, müssen Sie viele leichte Schritte machen, aber gelegentlich stoßen Sie auf einen Hügel. Selbst wenn die meisten Schritte leicht sind, wird Ihre Gesamtgeschwindigkeit durch diesen einen Hügel begrenzt.

Warum das wichtig ist

Die große Entdeckung ist, dass man keine komplizierten Erklärungen erfinden muss, um das Geschwindigkeitsproblem zu beheben. Man muss nicht annehmen, dass die Batterie voller „Klumpen" leerer Sitze ist oder dass das Material kaputt ist.

Der Artikel zeigt, dass Li₂MnO₃ tatsächlich ein sehr gutes Material ist (nahezu perfekt oder „stöchiometrisch"). Der Grund, warum wir in verschiedenen Experimenten unterschiedliche Geschwindigkeiten sehen, liegt einfach daran:

• Kurzstrecken-Experimente sehen nur die leichten, niedrigen Hügel (0,18 eV).
• Langstrecken-Experimente sehen die gesamte Reise, die durch den gelegentlichen höheren Hügel (0,50 eV) verlangsamt wird.

Das Fazit

Indem sie die Wärme und das magnetische „Wackeln" der Atome berücksichtigten, schufen die Autoren eine einzige, vereinheitlichte Geschichte. Sie bewiesen, dass sich die Lithium-Ionen auf lokaler Ebene leicht bewegen können, ihre gesamte Reise jedoch durch ein paar etwas härtere Schritte gesteuert wird. Dies erklärt, warum sich die Batterie je nach Messmethode unterschiedlich verhält, ohne dass man Defekte oder Verunreinigungen im Material beschuldigen muss.

Kurz gesagt: Die Batterie ist nicht kaputt; wir brauchten nur eine bessere Karte, die die Wärme und den magnetischen Tanz der Atome berücksichtigt, um zu verstehen, wie sich die Lithium-Ionen wirklich bewegen.

Technische Zusammenfassung

Technische Zusammenfassung: Aufdeckung der kurz- und langreichweitigen Li-Ionendiffusion in Li2MnO3 mittels dynamischer Mittelwertfeldtheorie bei endlichen Temperaturen

Problemstellung
Li2MnO3 ist ein kritischer Bestandteil von Li-überschüssigen geschichteten Kathoden, doch seine Rolle bei der Begrenzung des Li-Ionentransports bleibt umstritten. Experimentelle Abschätzungen der Li-Ionmigration in Li2MnO3 zeigen je nach Messskala eine erhebliche Diskrepanz. Mikroskopische $\mu^+$SR-Messungen an nahezu stöchiometrischen Pulvern berichten über eine niedrige Aktivierungsenergie für den Kurzstreckentransport ($E_a \approx 0,156$ eV) im paramagnetischen Regime. Im Gegensatz dazu liefern makroskopische Wechselstrom-Impedanzmessungen eine deutlich höhere scheinbare Barriere ($E_a \approx 0,46$ eV) für den Langstreckentransport.

Bisherige Studien auf Basis der Dichtefunktionaltheorie mit Hubbard-U-Korrekturen (DFT+U) hatten Schwierigkeiten, diese Werte in Einklang zu bringen. Standard-DFT+U-Rechnungen für eine einzelne Li-Leerstelle sagen typischerweise Barrieren von 0,5–0,8 eV voraus, die zu hoch sind, um die kurzreichweitigen experimentellen Daten zu erklären. Um die niedrige experimentelle Barriere zu erreichen, griffen frühere Modelle oft auf geklumpte Leerstellenkonfigurationen (Divakanzen oder Trivakanzen) zurück oder nahmen spezifische magnetische Ordnungen an. Solche Leerstellengruppen sind jedoch in nahezu stöchiometrischen Proben unwahrscheinlich, und Li2MnO3 ist unter Batteriebetriebsbedingungen paramagnetisch, wohingegen viele DFT+U-Studien ferromagnetische oder nichtmagnetische Zustände annahmen, die lokale Mn-Momente unterdrücken.

Methodik
Diese Arbeit verwendet einen multiskaligen rechnerischen Ansatz, um die Li+-Migration in paramagnetischem Li2MnO3 mit einer einzelnen Li-Leerstelle bei Raumtemperatur (300 K) zu untersuchen. Die Methodik integriert drei Komponenten:

1. DFT+U: Wird zur Optimierung von Migrationspfaden und Sattelpunktgeometrien mittels der Nudged-Elastic-Band (NEB)-Methode verwendet. Spin-polarisierte DFT+U-Rechnungen wurden unter Verwendung des PBEsol-Funktional und einer 1×2×2-Supercelle durchgeführt.
2. DFT+DMFT: Die dynamische Mittelwertfeldtheorie (DMFT) bei endlichen Temperaturen mit einem Continuous-Time Quantum Monte Carlo (CTQMC)-Störungs Solver wurde angewendet, um die Gesamtenergien entlang der festen NEB-Geometrien zu bewerten. Dieser Rahmen erfasst dynamische elektronische Korrelationen in der paramagnetischen Phase.
3. Umgang mit der elektronischen Struktur: Um die niedrige kristallographische Symmetrie (C2/m) von Li2MnO3 zu adressieren, die in Standard-Wannier-Projektionen zu großen nicht-diagonalen Matrixelementen führt, diagonalisierten die Autoren den Mn-d-Block. Dies schuf ein effektives niederenergetisches „nur-d"-Modell, bei dem die Wannier-Basis die Mn-d–O-p-Hybridisierung implizit kodiert und eine genaue Reproduktion der isolierenden Bandlücke ermöglicht.

Die Studie bewertet sechs symmetrie-inequivalente Migrationspfade (vier innerhalb der Schicht und zwei zwischen den Schichten), indem sie DMFT-Gesamtenergien für die Anfangs-, Sattelpunkt- und Endbilder berechnet und dabei den korrelierten elektronischen Zustand an jedem Punkt relaxieren lässt.

Hauptergebnisse
Die Anwendung dynamischer Korrelationen bei endlichen Temperaturen renormiert die Aktivierungsenergien für spezifische Migrationspfade signifikant und löst die Diskrepanz zwischen mikroskopischen und makroskopischen experimentellen Daten auf, ohne geklumpte Leerstellen heranzuziehen:

• Kurzreichweitige Diffusion: Für den Pfad mit der niedrigsten Barriere innerhalb der Schicht (Pfad ii: $4h \to T_{2b}^{Li} \to 4h$) wird die DFT+DMFT-Aktivierungsenergie auf 0,18 eV reduziert. Dieser Wert reproduziert quantitativ die bei $\mu^+$SR-Messungen beobachtete kurzreichweitige Aktivierungsenergie ($0,156$ eV). Die Reduktion wird primär durch den Eigenwert-Korrekturterm in der DMFT-Gesamtenergie angetrieben, der die Sattelpunkt-Konfiguration stabilisiert.
• Langreichweitige Diffusion: Der Pfad mit der nächstniedrigeren Barriere (Pfad iv: $4h \to T_{2b}^{Li} \to 2c$) weist unter DMFT eine Aktivierungsenergie von 0,50 eV auf. Diese Barriere steuert die langreichweitige Perkolierung von Li-Ionen. Sie ist konsistent mit der aus Wechselstrom-Impedanzmessungen extrahierten makroskopischen Aktivierungsenergie ($\approx 0,46$ eV).
• Pfadabhängigkeit: Die Barrierenreduktion ist stark pfadabhängig. Während Pfad (ii) einen dramatischen Abfall vom DFT+U-Wert (~0,68 eV) auf 0,18 eV zeigt, weisen andere Pfade (z. B. Pfad iii) nur minimale Änderungen auf. Die Autoren korrelieren dies mit der lokalen Registrierung des wandernden Li-Ions relativ zum MnO6-Netzwerk; Pfad (ii) verläuft durch einen offeneren Bereich mit längeren Li-Mn-Abständen, der durch dynamische Korrelationen effektiver stabilisiert wird.
• Elektronische Struktur: Im Gegensatz zu DFT+U, das eine schwach sitzungsabhängige Änderung der Mn-d-Besetzung bei der Erzeugung einer Leerstelle zeigt, offenbart DMFT eine stark sitzungsabhängige Umverteilung des spektralen Gewichts mit einer signifikanten Abnahme des besetzten Mn-d-Gewichts an Mn-Plätzen, die weiter von der Leerstelle entfernt sind.

Bedeutung und Behauptungen
Die Arbeit behauptet, eine vereinheitlichte mikroskopische Interpretation sowohl des lokalen als auch des makroskopischen Li-Ionentransports in Li2MnO3 zu liefern. Durch die Verwendung einer paramagnetischen DMFT-Beschreibung mit einer einzelnen Leerstelle demonstrieren die Autoren, dass die Hierarchie der Aktivierungsenergien (niedrig für Kurzstrecke, höher für Langstrecke) natürlich aus dem Zusammenspiel spezifischer Migrationspfade und dynamischer Korrelationen bei endlichen Temperaturen entsteht.

Die Studie kommt zu folgenden Schlussfolgerungen:

1. Die experimentell beobachtete kurzreichweitige Aktivierungsenergie kann durch die Migration einer einzelnen Leerstelle in einem paramagnetischen Wirt erklärt werden, wodurch die Notwendigkeit entfällt, extrinsische Unordnung oder geklumpte Leerstellenkonfigurationen zur Erklärung der Daten heranzuziehen.
2. Der Langstreckentransport wird durch einen Schritt mit höherer Barriere (0,50 eV) innerhalb eines Netzwerks von Hops gesteuert, was konsistent mit makroskopischen Impedanzdaten ist.
3. Dynamische Korrelationen bei endlichen Temperaturen sind ein wesentlicher Bestandteil für die Vorhersage der energetischen Verhältnisse der Ionenmigration in korrelierten Oxid-Elektroden, da statische DFT+U-Ansätze die notwendige Renormierung der Barrierenhöhen nicht erfassen können.

Die Autoren weisen darauf hin, dass diese Ergebnisse die Notwendigkeit von DFT-übergreifenden Behandlungen für die Iodendiffusion in korrelierten Oxiden unterstreichen und nahelegen, dass zukünftige Arbeiten diesen Rahmen auf höhere Leerstellenkonzentrationen und gemischtkationische Li-reiche geschichtete Kathoden ausweiten sollten, insbesondere dort, wo Sauerstoff-Redox aktiv wird.