Vacancy Diffusion Across FeCrAl Alloy Composition Space for Accident-Tolerant Fuel Cladding

Technische Samenvatting: Vacante Diffusie over de FeCrAl-Legeringscompositieruimte voor Accident-Tolerant Brandstofbekleding

Probleemstelling
IJzer-chroom-aluminium (FeCrAl) legeringen zijn vooraanstaande kandidaten voor accident-tolerant brandstofbekleding (ATF) in lichte waterreactoren (LWRs) vanwege hun superieure oxidatieweerstand bij hoge temperaturen vergeleken met Zircaloy. De in-reactor levensduur van deze bekledingen wordt echter bepaald door de door straling geïnduceerde microstructurele evolutie, waarbij puntdefecttransport—specifiek vacante diffusie—het dominante mechanisme is. Hoewel de vorming van beschermende oxidelagen goed is gekarakteriseerd, blijft de fundamentele defectfysica van FeCrAl onder bestraling slecht begrepen. Bestaande theoretische behandelingen vertrouwen vaak op geïdealiseerde of compositioneel vereenvoudigde omgevingen, die er niet in slagen de complexe lokale chemische fluctuaties te verklaren die inherent zijn aan geconcentreerde vaste oplossingen. Deze lokale omgevingen veranderen de vacante binding, migratiepaden en collectieve gedragingen zoals clustering en trapping aanzienlijk. Er is een kritieke behoefte om de individuele bijdragen van de constituerende elementen (Fe, Cr, Al) aan het totale vacante transport te ontrafelen om de langetermijnmateriaal-evolutie, fasestabiliteit en kruip onder gecombineerde thermische en bestralingsbelasting te voorspellen.

Methodologie
De auteurs maken gebruik van een soort-specifieke kinetische Monte Carlo (KMC) framework om vacante diffusie te modelleren over een systematisch bereik van Fe/Cr-composities (van Fe-rijk naar Cr-rijk) met een vaste concentratie van 5 at.% Al. De methodologie verloopt via de volgende stappen:

1. Potentiaal en Datageneratie: Interacties tussen atomen worden geëvalueerd met behulp van een Embedded Atom Method (EAM) potentiaal, die is aangepast aan een Density Functional Theory (DFT) database. Vacante migratiedata wordt gegenereerd met behulp van de "Hop-Decorate" workflow, waarbij een BCC Fe-referentiestructuur willekeurig wordt gedecoreerd met Fe, Cr en Al atomen volgens de beoogde globale composities. Minimale energiepaden worden berekend met de Climbing-Image Nudged Elastic Band (CI-NEB) methode voor de globale composities van Fe$_{80}$Cr$_{15}$Al$_5$, Fe$_{47.5}$Cr$_{47.5}$Al$_5$ en Fe$_{15}$Cr$_{80}$Al$_5$.
2. Surrogaatmodellering: Om de computationele kosten van on-the-fly NEB-berekeningen tijdens KMC te vermijden, zijn lineaire surrogaatmodellen getraind op de gegenereerde barrièredatabase. De lokale chemische omgeving voor elke sprong wordt gecodeerd met behulp van descriptoren die het aantal Fe-, Cr- en Al-atomen tellen binnen vaste radiale schillen rond het zadelpunt en het migrerende atoom. Afzonderlijke lineaire regressiemodellen zijn gefit voor elke migrerende soort (Fe, Cr, Al) om de soort-specifieke barrièreafhankelijkheden te vangen.
3. KMC-simulaties: De getrainde surrogaatmodellen drijven residence-time KMC-simulaties aan bij temperaturen variërend van 500 K tot 1500 K. Simulaties volgen enkelvoudige vacante trajecten over $10^5$ stappen om de gemiddelde kwadratische verplaatsing (MSD) te berekenen en diffusiecoëfficiënten ($D$), activeringsenergieën ($E_a$) en pre-exponentiële factoren ($D_0$) te extraheren.

Belangrijkste Resultaten

• Migratiebarrières en Lokale Omgeving: De studie bevestigt dat migratiebarrières zeer gevoelig zijn voor de lokale chemische omgeving. Terwijl puur Fe een barrière van ~0,63 eV vertoont, is die van Cr ~0,95 eV en die van Al ~0,69 eV. In legeringen verbreden de distributies van barrières zich aanzienlijk, met name voor Cr en Al, wat aangeeft dat lokale chemische asymmetrie de migratiepaden moduleert.
• Rol van Aluminium: Ondanks de lage concentratie (5 at.%) oefent Al een onevenredig grote invloed uit op de migratiebarrières. In alle drie de surrogaatmodellen verhoogt de aanwezigheid van een Al-atoom in de eerste naaste buur (1NN) de energiebarrière aanzienlijk (+0,048 eV voor Fe, +0,073 eV voor Cr, en +0,151 eV voor Al).
• Compositionele Afhankelijkheid van Diffusiviteit:
  • Fe-rijk naar Equiatomisch: In composities variërend van Fe$_{80}$Cr$_{15}$Al$_5$ tot Fe$_{32,5}$Cr$_{62,5}$Al$_5$ is de vacante diffusiviteit relatief ongevoelig voor de globale compositie. Activeringsenergieën convergeren naar een bereik van 0,70–0,84 eV, wat suggereert dat Fe-migratie het transport domineert in deze regimes.
  • Cr-rijk Regime: De Cr-rijke compositie (Fe$_{15}$Cr$_{80}$Al$_5$) vertoont een merkbaar verhoogde activeringsenergie van 1,135 eV. Bijgevolg is de macroscopische diffusiviteit orders van grootte lager dan die van Fe-rijke legeringen bij lagere temperaturen (bijv. ~$10^{-17}$ m$^2$ s$^{-1}$ bij ~570 K vergeleken met ~$10^{-13}$ m$^2$ s$^{-1}$ voor Fe-rijke legeringen).
• Percolatie en "Sluggish" Diffusie: De resultaten suggereren een percolatiedrempel waarbij het verbonden netwerk van paden met lage barrières (gedomineerd door Fe-sprongen) afbreekt naarmate het Cr-gehalte toeneemt. In het Cr-rijke regime wordt de vacante migratie gedomineerd door hoog-barrière Cr-sprongen, wat leidt tot "sluggish" (traag) diffusie. De auteurs observeren dat vacante trajecten in chemisch ruwe landschappen betrokken zijn bij het ontsnappen uit trapping-pockets, wat resulteert in niet-lineair gemiddelde kwadratische verplaatsingsgedrag voordat een diffuus regime wordt bereikt.

Betekenis en Claims
Het artikel beweert dat de macroscopische diffusiviteit van FeCrAl-legeringen niet voorspeld kan worden door een eenvoudige compositie-gewogen gemiddelde van de zuivere end-member barrières; in plaats daarvan wordt deze bepaald door de heterogeniteit van het lokale energie Landschap en de specifieke connectiviteit van migratiepaden met lage barrières.

• Ontwerpimplicaties: De studie identificeert een onderscheidend compositioneel regime (Cr-rijk) waar de vacante mobiliteit wordt onderdrukt. Deze onderdrukking impliceert dat in Cr-rijke $\alpha'$-fasen (die bekend staan om vorming onder bestraling), vacantes minder mobiel zijn, wat potentieel de snelheid van vacante-interstitieel recombinatie vermindert en de kinetiek van defectaccumulatie verandert.
• Accident Tolerance: Onder LOCA-condities (Loss-of-Coolant Accident) waarbij de temperaturen boven 1200 K stijgen, wordt de compositionele afhankelijkheid van diffusiviteit uitgesproken. De auteurs suggereren dat de Cr-rijke compositie, met zijn hoge activeringsenergie en onderdrukte diffusiviteit, gunstige vacante-kinetica kan bieden voor accident-tolerant bekleding, wat potentieel complementair is aan de bekende oxidatieweerstand.
• Methodologische Bijdrage: Het werk demonstreert dat gereduceerde kinetische Monte Carlo-frameworks, geparametriseerd door soort-specifieke lineaire surrogaatmodellen getraind op NEB-databases, effectief de defecttransport over complexe compositieruimtes kunnen in kaart brengen. Deze aanpak biedt een mechanistische link tussen atomistische processen en mesoschaalmodellen die worden gebruikt in brandstofprestatiecodes, waarmee men verder gaat dan empirische trial-and-error in legeringsontwerp.

Beperkingen Genoemd door Auteurs
De auteurs merken expliciet op dat het model het rooster behandelt als statisch en bevroren, waarbij lokale relaxatie en stralings-geïnduceerde compositionele evolutie (segregatie) over tijd worden genegeerd. Het lineaire surrogaatmodel kan geen directionele of symmetrie-brekende effecten of magnetische bijdragen aan het energie Landschap vastleggen. Bovendien richt de studie zich uitsluitend op vacante transport en behandelt het geen interstitiële cluster-transport of de stralings-geïnduceerde transformatie van dislocatielussen, die ook de defect sink-densiteiten kunnen beïnvloeden.