Ion Jump Motion as the Background for Muon Diffusion in Battery Materials Research Using $μ$SR

Probleemstelling

Muon Spin Relaxatie (µSR) wordt al vijftien jaar gebruikt om ion-dynamica (zoals die van Li- en Na-ionen) in batterijmaterialen te onderzoeken. De standaardanalyseprocedure maakt gebruik van de conventionele Kubo-Toyabe (KT) relaxatiefunctie, $G_z^{KT}(t; \Delta_{KT}, \nu_{KT})$. Hierbij wordt de fluctuatiefrequentie $\nu_{KT}$ geïnterpreteerd als de springfrequentie van de ionen ($\nu_i$).

Echter, recente studies hebben strijdige resultaten opgeleverd:

1. Anomalieën: In veel experimenten vertoont $\nu_{KT}$ een onverklaarbare piek bij een specifieke temperatuur $T^*$, gevolgd door een daling. Tegelijkertijd neemt de lijnbreedte $\Delta_{KT}$ boven $T^*$ af. Dit is in tegenspraak met de verwachting dat de springfrequentie van ionen monotoon toeneemt met de temperatuur (Arrhenius-gedrag).
2. Foutieve interpretatie: De afgeleide activeringsenergieën ($E_i$) en prefactoren uit deze analyses wijken vaak sterk af van theoretische verwachtingen voor ionendiffusie.
3. Beperking van het model: De conventionele KT-functie gaat uit van een autocorrelatiefunctie die naar nul gaat bij $t \to \infty$. Dit is echter onjuist voor ionendiffusie waarbij slechts een deel van de ionen rondom de muon springt (een "Edwards-Anderson" type situatie), wat leidt tot een eindige correlatie op lange tijdschalen.

Methodologie

De auteur voert numerieke simulaties uit om de oorsprong van deze anomalieën te achterhalen en een correcter analysemodel te testen.

• Uitgebreid Kubo-Toyabe Model ($G_z^{EA}$): In plaats van de conventionele KT-functie, wordt gebruik gemaakt van een uitgebreide relaxatiefunctie die een Edwards-Anderson parameter $Q$ incorporeert. De autocorrelatiefunctie wordt beschreven als:
  $$C(t) \approx [(1-Q) + Q e^{-\nu_i t}] e^{-\nu_\mu t}$$
  Hierbij staat $Q$ voor het aandeel van de lokale veldfluctuaties veroorzaakt door de springende ionen (bijv. Li/Na), en $1-Q$ voor het statische aandeel van andere ionen. De term $e^{-\nu_\mu t}$ beschrijft de fluctuaties door de beweging van de muon zelf.
• Monte Carlo Simulaties: Er worden tijdspectra gegenereerd voor een temperatuurbereik van 10–400 K, met gebruikmaking van Arrhenius-vergelijkingen voor zowel de ionenspringfrequentie ($\nu_i$) als de muonspringfrequentie ($\nu_\mu$).
  • Parameters: $\nu_{i0} = 5 \times 10^{12} s^{-1}$, $E_i = 0.35$ eV (typisch voor ionendiffusie).
  • Parameters: $\nu_{\mu0} = 2 \times 10^{6} s^{-1}$, $E_\mu = 0.02$ eV (typisch voor kwantumtunneling van muonen).
• Pseudo-Data Analyse: De gegenereerde spectra (die de "waarheid" vertegenwoordigen) worden vervolgens geanalyseerd met de conventionele KT-functie ($G_z^{KT}$), precies zoals in eerdere literatuur is gedaan. Dit simuleert de huidige analysepraktijk om te zien welke fouten hierdoor ontstaan.

Belangrijkste Bijdragen en Resultaten

1. Oorzaak van de Anomalieën:
   De simulaties tonen aan dat de waargenomen piek in $\nu_{KT}$ bij $T^*$ en de daaropvolgende daling van $\Delta_{KT}$ niet het gevolg zijn van een verandering in het diffusiemechanisme van de ionen, maar van een kruispunt (crossover) in de dynamica:

   • Bij lage temperaturen is $\nu_i$ laag; de relaxatie wordt gedomineerd door de statische veldverdeling.
   • Rond $T^*$ neemt $\nu_i$ sterk toe en wordt deze vergelijkbaar met de muon-springfrequentie $\nu_\mu$ en de lijnbreedte $\Delta$.
   • De piek in $\nu_{KT}$ treedt op wanneer $\nu_i \approx 3\sqrt{Q}\Delta$.
   • De daling van $\Delta_{KT}$ boven $T^*$ wordt veroorzaakt door "motional narrowing" (bewegingsverenging) van het deel van het veld dat door de ionen wordt veroorzaakt ($Q$-component), terwijl het statische deel ($1-Q$) overblijft. De lijnbreedte convergeert naar $\sqrt{1-Q}\Delta$.

2. Validatie van het Model:
   De simulaties reproduceren semi-kwantitatief de anomalieën die in de literatuur voor materialen zoals $Li_xCoO_2$, $LiFePO_4$ en $NaC_x$ zijn waargenomen (zie Tabel I in het artikel). Dit bevestigt dat deze patronen een handtekening zijn van ionensprongen die binnen het tijdsvenster van de µSR-meting vallen.

3. Herinterpretatie van Bestaande Data:

   • Materialen met anomalieën: Voor materialen waar zowel de piek in $\nu_{KT}$ als de daling in $\Delta_{KT}$ worden waargenomen, is er sprake van detecteerbare ionendiffusie. Echter, de waarden van $\nu_{KT}$ en $\Delta_{KT}$ zoals berekend met de conventionele KT-functie zijn geen directe maat voor de ionenspringfrequentie $\nu_i$.
   • Materialen zonder anomalieën: Als deze kenmerken ontbreken, suggereert het artikel dat de µSR-metingen geen informatie over ionendiffusie bevatten. In plaats daarvan reflecteren deze resultaten waarschijnlijk uitsluitend de zelf-diffusie van de muon ($\nu_\mu$). Dit betekent dat veel eerdere studies die ionendiffusie claimden op basis van een monotoon toenemend $\nu_{KT}$, mogelijk onterecht zijn.

4. Nieuwe Analysestrategie:
   De auteur stelt voor om bestaande data opnieuw te analyseren met de uitgebreide functie $G_z^{EA}(t)$. Dit maakt het mogelijk om $\nu_i$ en $\nu_\mu$ correct te scheiden en te evalueren.

   • Hoewel $\nu_{KT}$ zelf geen goede maat is voor $\nu_i$, kan de temperatuurafhankelijkheid van $\Delta_{KT}$ wel worden gebruikt om een ruwe schatting van de activeringsenergie $E_i$ te maken (via een tweestapsmodel).

Significantie en Conclusie

Dit artikel vormt een fundamentele correctie op de interpretatie van µSR-data in de batterijonderzoeksgemeenschap.

• Paradigmaverschuiving: Het toont aan dat de conventionele analyse (KT-functie) ongeschikt is voor systemen waar ionen en muonen beiden diffunderen, maar met verschillende tijdschalen en mechanismen.
• Kritische Evaluatie: Veel eerder gepubliceerde resultaten over ionendiffusie in batterijmaterialen moeten worden herzien. Het is zeer waarschijnlijk dat veel van deze resultaten in werkelijkheid muon-diffusie beschrijven en geen ionen-diffusie.
• Toekomstperspectief: Hoewel dit betekent dat sommige eerdere conclusies over ionenbeweging mogelijk onjuist zijn, is de studie niet zinloos. De data kan waardevol zijn voor het begrijpen van de rol van "verborgen waterstof" (waarbij muonen fungeren als pseudo-waterstof) in oxide-materialen, wat relevant is voor halfgeleiders en diëlektrica.

Kortom, de paper pleit voor een zorgvuldige heranalyse van vintage µSR-data met het uitgebreide $G_z^{EA}$-model om de werkelijke ionendiffusieparameters ($\nu_i$) en muon-diffusieparameters ($\nu_\mu$) correct te scheiden.