Reduced-Order Modelling of Defect Transport using Surrogate Kinetics: Application to U$_x$Pu$_ {1-x}$N

Technische Samenvatting: Reduced-Order Modellering van Defecttransport met behulp van Surrogate-kinetiek in $U_xPu_{1-x}N$

Probleemstelling
Het modelleren van defecttransport in chemisch gedisordeerde materialen, zoals gemengde actinidenitriden ($U_xPu_{1-x}N$), vormt een aanzienlijke uitdaging. In chemisch uniforme roosters kunnen migratiepaden van defecten worden opgesomd als een eindige set overgangssnelheden. In gedisordeerde systemen varieert de lokale chemische omgeving echter continu, wat een distributie van migratiebarrières creëert in plaats van een gesloten set. Deze chemische complexiteit maakt een volledige enumeratie van omgevingen via Density Functional Theory (DFT) of Moleculaire Dynamica (MD) computationeel onhaalbaar. Bovendien kunnen standaard continuous-time Markov chain (CTMC) modellen de effectief oneindige variëteit aan overgangssnelheden, gegenereerd door lokale wanorde, niet accommoderen. Dit probleem is cruciaal voor geavanceerde nucleaire brandstoffen zoals uraniummononitride (UN), die kan evolueren naar $U_xPu_{1-x}N$ vaste oplossingen door splijting of als gemengde brandstoffen kan worden geproduceerd. Het voorspellen van brandstofprestaties (zwelling, vrijgave van splijtingsgassen, degradatie van thermische geleidbaarheid) vereist begrip van hoe compositionele variatie de defectenergetica en het transport verandert, terwijl directe atomistische simulaties over de enorme configuratieruimte onhaalbaar zijn.

Methodologie
De auteurs stellen een reduced-order workflow voor om atomistische mechanismen te integreren in computationeel efficiënte functionele vormen. De methodologie bestaat uit vier hoofdfasen:

1. Datageneratie (Hop-Decorate): Met behulp van de Hop-Decorate (HopDec) workflow, gebouwd op LAMMPS en ASE, hebben de auteurs een database gegenereerd van ongeveer 3.000 migratiebarrières ($E_b$) en energieverschillen ($\Delta E$). Dit werd bereikt door prototypische nearest-neighbor (1NN) migratiepaden in referentie UN-structuren te identificeren en deze 1.000 keer te "re-decoreren" bij globale composities van $U_{0.5}Pu_{0.5}N$, $U_{0.1}Pu_{0.9}N$ en $U_{0.9}Pu_{0.1}N$. Interatomaire energieën werden berekend met een Finnis–Sinclair EAM-potentiaal gecombineerd met een Buckingham-potentiaal, gevalideerd tegen DFT-data. De Climbing-Image Nudged Elastic Band (CI-NEB) methode werd gebruikt om de energetica voor elke configuratie te evalueren.
2. Constructie van het Surrogate-model: De lokale chemische omgeving werd gedefinieerd door de tellingen van Uraniumatomen binnen specifieke nearest-neighbor (NN) schillen rond het defect in de initiële, finale en zadelpuntconfiguraties. Lineaire surrogate-functies werden aangepast om $E_b$ en $\Delta E$ te voorspellen op basis van deze eenvoudige descriptoren:
   • $\Delta E$ werd gemodelleerd als een lineaire functie van de verandering in U-tellingen tussen de minima-schillen.
   • $E_b$ werd gemodelleerd als een lineaire functie van U-tellingen in de minima- en zadelpunt-schillen, plus een soortspecifieke baseline-correctie ($\beta$).
   • Voor U-gemedieerde en Pu-gemedieerde sprongen op het actiniden-subrooster waren aparte fits vereist vanwege de sterke soortafhankelijkheid.
3. Kinetic Monte Carlo (KMC) Simulatie: De surrogate-modellen werden geïntegreerd in een aangepaste lattice-based KMC (laKMC) simulator. Simulaties werden uitgevoerd op $10 \times 10 \times 10$ nm$^3$ supercellen (64.000 atomen) voor 100.000 stappen. Gemiddelde kwadratische verplaatsingen (MSDs) werden geëxtraheerd om langetermijn-diffusiviteiten ($D$) te berekenen over verschillende temperaturen en composities, met statistische middeling over 50 trajecten per conditie.
4. Trapping-analyse: Om de atomaire oorsprong van macroscopisch transport te verduidelijken, analyseerden de auteurs de kinetisch-thermodynamische ruimte door de disorder-gemiddelde gemiddelde verblijftijden ($\langle \tau \rangle_c$) en de op snelheid gewogen gemiddelde energieveranderingen ($\langle \Delta E \rangle_c$) voor elke compositie te berekenen.

Belangrijkste Resultaten
De studie analyseerde drie defectsoorten: actinide-vacatures ($V_A$), nitride-vacatures ($V_N$) en actinide-nitride divacatures ($V_{A-N}$).

• Migratie-energetica: De surrogate-modellen bereikten lage root-mean-squared errors (RMSE), wat aangeeft dat de migratie-energetica in $U_xPu_{1-x}N$ primair worden beheerst door kort-bereik chemische identiteit in plaats van complexe lang-bereik collectieve effecten.

  • $V_A$ (Actinide-vacature): Migratie wordt sterk door de soort gecontroleerd. De barrièreverdeling splitst zich op op basis van of het migrerende atoom U of Pu is. De lokale omgeving veroorzaakt slechts een lichte verbreding; het primaire compositionele effect is de statistische weging van U- versus Pu-sprongen.
  • $V_N$ (Nitride-vacature): Migratie vertoont genuante milieu-afhankelijke modulatie. De barrièreverdeling evolueert continu met de compositie en vertoont multi-lobed kenmerken bij intermediaire composities, wat overeenkomt met onderscheidende lokale actinide-configuraties. Opvallend genoeg kunnen kleine hoeveelheden U in een Pu-rijke matrix lokaal gunstige, laag-barrière paden creëren.
  • $V_{A-N}$ (Divacature): Het gekoppelde defect vertoont gedrag dat lijkt op de som van geïsoleerde vacatures, waarbij migratie plaatsvindt via een van beide subroosters. De aanwezigheid van het paar induceert een systematische verschuiving in de barrièrehoogte, maar verandert de gevoeligheid voor lokale coördinatie niet fundamenteel.

• Diffusietrends:

  • Niet-lineariteit: Diffusiviteit volgt geen eenvoudige lineaire interpolatie (Wet van Vegard) tussen de UN- en PuN-eindleden.
  • $V_A$ Gedrag: Diffusiviteit vertoont een scherpe, niet-lineaire toename bij kleine toevoegingen van Pu. Dit wordt gedreven door de percolatie van snelle Pu-gemedieerde sprongen. Naarmate de Pu-concentratie toeneemt, gaat het systeem over van een beperking door trage U-sprongen naar een netwerk van snelle Pu-sprongen.
  • $V_N$ Gedrag: Diffusiviteit vertoont "trage diffusie" bij lage Pu-concentraties, waarbij Pu fungeert als een kinetische val die het netto transport vermindert, ondanks de aanwezigheid van enkele laag-barrière gebeurtenissen. Bij hogere Pu-concentraties neemt de diffusiviteit toe naarmate laag-barrière paden prominenter en verbonden worden.
  • $V_{A-N}$ Gedrag: De diffusiviteit van de divacature weerspiegelt een competitie tussen de chemische gevoeligheden van beide subroosters, wat resulteert in een complexe, niet-lineaire compositieafhankelijkheid die verschilt van geïsoleerde vacatures.

• Trapping-analyse: De kinetisch-thermodynamische kaarten ($\langle \tau \rangle_c$ vs. $\langle \Delta E \rangle_c$) onthulden dat:

  • Voor $V_A$ vermindert het verhogen van het Pu-gehalte de verblijftijd zonder de energetische bias significant te veranderen, wat de percolatie van snelle paden bevestigt.
  • Voor $V_N$ verhogen kleine Pu-concentraties de energetische bias ($\langle \Delta E \rangle_c$) terwijl de verblijftijd stabiel blijft, wat aangeeft dat er chemische vallen ontstaan die relatieve energieën moduleren in plaats van alleen barrières te verhogen.

Betekenis en Claims
Het artikel claimt een algemeen, computationeel efficiënt kader te bieden voor het modelleren van defecttransport in chemisch gedisordeerde materialen zonder terug te vallen op "black-box" machine learning-modellen die een gebrek aan fysieke interpreteerbaarheid hebben.

• Reduced-Order Representatie: De auteurs demonstreren dat de complexe migratie-energetica van gedisordeerde actinidenitriden een onverwacht eenvoudige reduced-order representatie toelaten met behulp van lineaire functies van lokale coördinatie-tellingen.
• Mechanistisch Inzicht: Het werk identificeert twee dominante mechanismen die het transport beheersen: soort-gecontroleerde kinetiek (dominant voor actinide-vacatures) en milieu-afhankelijke trapping (dominant voor nitride-vacatures). Het benadrukt dat macroscopische diffusiviteit voortkomt uit de competitie tussen barrièreverdelingen, pad-connectiviteit en verblijftijden, in plaats van simpelweg uit gemiddelde barrièrewaarden.
• Schaalbaarheid: Door atomistische data te distilleren tot compacte surrogate-functies, maakt de aanpak de integratie van atomistische mechanismen in hogere-schaal simulaties mogelijk (bijv. cluster dynamics of continuum fuel performance codes). Dit maakt het efficiënt modelleren van evoluerende composities tijdens brandstofverbruik of transmutatie mogelijk.
• Toepasbaarheid: Hoewel gedemonstreerd op $U_xPu_{1-x}N$, wordt het framework gepresenteerd als toepasbaar op andere compositioneel gedisordeerde materialen, zoals high-entropy legeringen en gedoteerde halfgeleiders, waar brute-force enumeratie van lokale omgevingen onuitvoerbaar is.

De auteurs merken beperkingen op, waaronder de beperking van divacatures tot nearest-neighbor paren (waarbij dissociatie/re-associatie wordt genegeerd) en het gebruik van empirische potentialen die elektronische effecten (bijv. 5f-orbitaal lokalisatie in Pu) niet expliciet behandelen. Ze stellen echter dat de huidige aanpak er succesvol in slaagt de kloof te overbruggen tussen atomistische mechanismen en hanteerbare transportvergelijkingen voor voorspellende brandstofprestatie-modellering.