A comparison of simulation tools for Muon-Induced X-ray Emission (MIXE) in thin films: a study case with lithium batteries

Deze studie vergelijkt de simulatietools SRIM, GEANT4 en PHITS voor muon-geïnduceerde röntgenemissie in dunne lagen en concludeert dat SRIM en PHITS betrouwbare hulpmiddelen zijn voor het voorspellen van muonimplantatieprofielen, terwijl PHITS ondanks een systematische energieverschuiving veelbelovend is voor element-specifieke spectroscopie na correctie.

De kern

Stel je voor dat je een complexe, gelaagde taart wilt onderzoeken zonder hem aan te snijden of te beschadigen. Je wilt precies weten wat er in elke laag zit: de bodem, de vulling, de room, de vruchten. Normaal gesproken zou je de taart moeten doorsnijden, maar dat is niet altijd mogelijk (bijvoorbeeld bij een kostbare antieke vaas of een batterij die nog werkt).

De wetenschappers in dit artikel hebben een magische "röntgenflits" gebruikt om door de lagen heen te kijken. Ze noemen dit MIXE (Muon-Induced X-ray Emission). Maar om deze techniek goed te kunnen gebruiken, moesten ze eerst drie verschillende "rekenmachines" (simulatiesoftware) testen om te zien welke het beste voorspelt waar die flitsen precies stoppen en wat ze zien.

Hier is de uitleg in simpele taal:

1. De Magische Deeltjes: Muonen

Stel je voor dat je een regen van heel zware, snelle deeltjes (muonen) op je taart (of in dit geval, een lithium-ion batterij) laat regenen.

• Hoe het werkt: Deze deeltjes zijn zwaar en sneller dan elektronen. Als ze de batterij binnenkomen, botsen ze tegen atomen aan en verliezen ze snelheid, tot ze ergens tot stilstand komen.
• Het effect: Op het moment dat ze stoppen, worden ze "gevangen" door de atoomkernen. Ze maken een soort van "atoomdans" (een cascade) waarbij ze energie kwijtraken in de vorm van heel krachtige röntgenstraling.
• De boodschap: Elke stof (zoals lithium, nikkel, koper) geeft een heel specifiek geluid (een bepaalde energie) af. Als je luistert naar deze geluiden, weet je precies welke lagen er zijn en hoe diep ze zitten.

2. Het Probleem: De Rekenmachines

Om te weten hoe diep de deeltjes gaan en welke signalen ze moeten geven, gebruiken wetenschappers computersimulaties. In dit artikel hebben ze drie populaire programma's getest:

1. SRIM: Een snelle, simpele rekenmachine.
2. GEANT4: De "gouden standaard", een zeer complexe en nauwkeurige simulator.
3. PHITS: Een krachtige simulator die ook kan voorspellen welke "geluiden" (straling) er worden geproduceerd.

Ze wilden weten: Welke software geeft het meest waarheidsgetrouwe beeld van wat er in de batterij gebeurt?

3. De Vergelijking: Wie doet het goed?

SRIM: De snelle schatting

• Analogie: Stel je voor dat SRIM een snelle schatting is van hoe ver een bal rolt in het gras. Hij is snel en handig, maar als er een groot stuk gras (lucht) voor de taart ligt, wordt hij een beetje slordig.
• Resultaat: Hij is goed voor simpele, dichte materialen, maar faalt als er veel lucht of dunne lagen voor de batterij zitten. Hij kan de "atoomdans" (de straling) niet simuleren.

GEANT4: De nauwkeurige architect

• Analogie: Dit is de gedetailleerde architect die elke steen en elke kieren in de muur berekent.
• Resultaat: Hij is extreem nauwkeurig in het voorspellen van hoe diep de deeltjes gaan, zelfs als de lagen heel verschillend zijn (zoals lucht vs. metaal). Maar hij is zwaar en traag om te gebruiken.

PHITS: De veelzijdige alleskunner

• Analogie: PHITS is als een slimme robot die zowel de diepte berekent als de "muziek" (de straling) voorspelt.
• Resultaat:
  • Diepte: Hij doet het bijna net zo goed als GEANT4.
  • Muziek (Straling): Hij kan de spectra (de "noten") voorspellen, maar er zit een klein foutje in. De "noten" die hij produceert voor zware elementen (zoals nikkel en koper) zijn net iets te hoog van toon dan ze zouden moeten zijn. Het is alsof een pianist de toetsen net iets te ver naar rechts duwt.
  • Toch nuttig: Ondanks die toonhoogte-fout, herkent hij de juiste "melodieën". Hij kan dus wel zeggen: "Ah, hier zit nikkel!" en "Hier zit koper!", en hij kan zelfs zeggen hoeveel van elk er is.

4. De Conclusie: Wat betekent dit voor ons?

De onderzoekers hebben ontdekt dat:

• SRIM prima is voor snelle, ruwe schattingen van hoe diep de deeltjes gaan, maar niet voor complexe situaties met lucht.
• GEANT4 en PHITS zijn de beste keuzes voor nauwkeurige diepte-metingen.
• PHITS is de winnaar voor de toekomst, omdat hij ook de chemische samenstelling kan voorspellen. Als de "toonhoogte-fout" (de energie-verschil) wordt gecorrigeerd, kan PHITS een krachtig hulpmiddel worden om batterijen en andere complexe objecten volledig in kaart te brengen zonder ze open te maken.

Kort samengevat:
De wetenschappers hebben getest hoe goed computers kunnen voorspellen waar magische deeltjes stoppen in een batterij. Ze hebben gevonden dat PHITS de beste "twee-in-één" tool is: hij ziet hoe diep de deeltjes gaan én hij kan de chemische samenstelling "luisteren", mits we een klein rekenfoutje in zijn software nog even rechtzetten. Dit helpt wetenschappers om in de toekomst batterijen en kostbare objecten beter te begrijpen zonder ze te beschadigen.

Technische samenvatting

Titel en Context

Onderwerp: Vergelijkende studie van simulatietools (SRIM, GEANT4, PHITS) voor het modelleren van Muon-Induced X-ray Emission (MIXE) in dunne lagen, met als casestudie lithiumbatterijen.
Auteurs: Maxime Lamotte et al. (Paul Scherrer Institute, Zwitserland).
Doel: Het identificeren van betrouwbare en praktische simulatietools voor het voorspellen van muon-stopdieptes en muonische X-stralenspectra in gelaagde materialen, specifiek ter ondersteuning van experimenten aan de SµS-faciliteit (PSI).

---

1. Het Probleem

MIXE is een niet-destructieve techniek die negatieve muonen ($\mu^-$) gebruikt om elementaire samenstelling en isotopen te analyseren op verschillende dieptes in een monster. Muonen worden gestopt in het materiaal, vormen een "muonisch atoom" en vervallen via een cascade die karakteristieke X-stralen uitzendt.

• Uitdaging: Voor de interpretatie van MIXE-data en het ontwerpen van experimenten zijn robuuste simulatietools nodig die de transport, stopping (remmen) en de atomaire cascade van muonen nauwkeurig kunnen modelleren in complexe, gelaagde systemen (zoals batterijen met verschillende dichtheden).
• Gevraagde tools: Er is behoefte aan tools die zowel de stopdiepte-profielen als de spectra (energie en intensiteit) correct voorspellen, rekening houdend met de specifieke geometrie van de experimenten (bijv. de GIANT-opstelling aan het PSI).

2. Methodologie

De auteurs hebben drie Monte Carlo-simulatiekaders vergeleken: SRIM, GEANT4 en PHITS.

• Benchmark: Een lithium-ion batterijcel (multilayer systeem) bestaande uit lagen zoals een pouch, aluminium, NMC (kathode), separators (Whatman/Celgard), grafiet, koper en een pouch-laag.
• Experimentele Referentie: Data afkomstig van de GIANT spectrometer aan het PSI, waarbij muonen met instelbare impulsen (15–60 MeV/c) in de batterij worden geïmplanteerd.
• Simulatie-instellingen:
  • SRIM: Omdat SRIM geen native muon-ondersteuning heeft, werden "pseudo-muonen" gemodelleerd als protonen met een geschaalde massa. SRIM simuleert alleen elektromagnetisch remmen, niet de kernvangst.
  • GEANT4: Gebruikt als referentie-toolkit voor deeltjestransport (fysiek bibliotheek: QGSP_BIC_HP_EMZ).
  • PHITS: Gebruikt met geüpgradede muon-interactiemodules, inclusief kernvangst en de Akylas-Vogel routine voor de muonische cascade.
• Validatie: De simulaties werden vergeleken op basis van:
  1. Gemiddelde stopdiepte (Mean Stopping Depth - MSD) en spreiding (range straggling).
  2. Voorspelde energieën van muonische X-stralen (K- en L-lijnen).
  3. Relatieve intensiteiten van de spectra.

3. Belangrijkste Bijdragen

• Vergelijking van Transport: Eerste uitgebreide vergelijking van deze drie codes voor MIXE-toepassingen in multilayer-systemen met scherpe dichtheidscontrasten (bijv. lucht naar vaste stof).
• Identificatie van Systematische Fouten: Het opsporen van een systematische energieverschuiving in de PHITS-spectra voor zware elementen, veroorzaakt door een schalingsprobleem in de cascade-routine (aama.f).
• Validatie van Intensiteiten: Aantonen dat hoewel de absolute energieën in PHITS afwijken, de relatieve intensiteiten en spectrale vormen zeer goed overeenkomen met theorie en experiment.
• Praktische Richtlijnen: Het bieden van een leidraad voor onderzoekers bij het kiezen van de juiste tool voor hun specifieke MIXE-experiment (snelle schatting vs. volledige spectrale simulatie).

4. Resultaten

A. Muon Stopdieptes (Transport)

• GEANT4 vs. PHITS: Beide tools leveren bijna identieke resultaten voor de gemiddelde stopdiepte en de spreiding, zelfs in complexe geometrieën met grote luchtspaties en dunne lagen. De afwijkingen zijn verwaarloosbaar (< 1%).
• SRIM:
  • Voor geïsoleerde lagen komt SRIM goed overeen met PHITS (afwijking < 5%).
  • Beperking: SRIM faalt bij systemen met grote lage-dichtheidsgebieden (zoals de 10 cm luchtspatie voor de sample) omdat het de gemiddelde diepte beïnvloedt door trajecten in de lucht en door numerieke precisieproblemen (vaste eenheid van $\mu$m). Het kan geen realistische optica of kernvangst simuleren.
  • Conclusie: SRIM is nuttig voor snelle, ruwe schattingen van bereik in compacte geometrieën, maar niet voor nauwkeurige transportmodellen in complexe opstellingen.

B. X-stralenspectra (Cascade)

• Energieverschuiving: PHITS vertoont een systematische energieoverschatting voor de K-lijnen van medium- en zware elementen (bijv. Ni, Cu) ten opzichte van theoretische waarden (MuDirac) en experimentele data.
  • Bijvoorbeeld: Voor $^{58}$Ni K$\alpha_1$ voorspelt PHITS ~1673 keV, terwijl MuDirac en experiment ~1432 keV aangeven (een afwijking van ~17%).
  • Deze fout neemt toe met het atoomnummer (Z) en is minder groot voor L-lijnen.
• Relatieve Intensiteiten: Ondanks de energieverschuiving, reproduceert PHITS de relatieve intensiteiten tussen verschillende lijnen en elementen zeer goed. De verhouding tussen K- en L-lijnen en de spectrale vorm zijn betrouwbaar.
• Isotopische verschuivingen: PHITS onderscheidt isotopische verschuivingen (bijv. tussen $^{58}$Ni en $^{60}$Ni), maar onder- of overschat de grootte ervan in vergelijking met MuDirac.

5. Betekenis en Conclusie

De studie concludeert dat:

1. GEANT4 en PHITS de meest geschikte tools zijn voor kwantitatieve MIXE-analyse, omdat ze zowel transport als cascade accurate modelleren.
2. SRIM een waardevol hulpmiddel blijft voor snelle, voorlopige schattingen van stopdieptes in individuele lagen, mits de beperkingen bij lage dichtheden worden begrepen.
3. PHITS potentieel heeft voor voorspellende MIXE-spectroscopie, maar de huidige versie vereist correctie van de energie-bias in de cascade-routine voor zware elementen. Een praktische oplossing zou zijn om PHITS te koppelen aan tabellen met correcte overgangsenergieën (bijv. gegenereerd door MuDirac) terwijl de PHITS-logica voor relatieve intensiteiten behouden blijft.

Impact: Deze resultaten vormen de basis voor de ontwikkeling van een webgebaseerd simulatieplatform voor toekomstige MIXE-gebruikers en helpen bij het optimaliseren van experimenten voor dieptegestructureerde analyse van complexe objecten zoals batterijen, archeologische artefacten en geologische monsters.