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

Technische Zusammenfassung: Vakanzendiffusion über den FeCrAl-Legierungszusammensetzungsraum für unfalltolerante Brennelementumhüllungen

Problemstellung
Eisen-Chrom-Aluminium-Legierungen (FeCrAl) sind führende Kandidaten für unfalltolerante Brennelementumhüllungen (ATF) in leichten Wasserreaktoren (LWR) aufgrund ihrer überlegenen Hochtemperatur-Oxidationsbeständigkeit im Vergleich zu Zircaloy. Die In-Reaktor-Lebensdauer dieser Umhüllungen wird jedoch durch die strahlungsinduzierte mikrostrukturelle Entwicklung bestimmt, wobei der Transport von Punktdefekten – spezifisch die Vakanzendiffusion – der dominante Mechanismus ist. Während die Bildung schützender Oxidschichten gut charakterisiert ist, bleibt die grundlegende Defektphysik von FeCrAl unter Bestrahlung unzureichend verstanden. Bestehende theoretische Behandlungen stützen sich oft auf idealisierte oder zusammensetzungsvereinfachte Umgebungen und versäumen es, die komplexen lokalen chemischen Fluktuationen zu berücksichtigen, die in konzentrierten Mischkristalllegierungen inhärent sind. Diese lokalen Umgebungen verändern signifikant die Vakanzenbindung, die Migrationspfade und das kollektive Verhalten wie Clusterbildung und Trapping. Es besteht ein kritischer Bedarf, die individuellen Beiträge der Bestandteile (Fe, Cr, Al) zum Gesamtvakanzentransport zu entwirren, um die langfristige Materialentwicklung, Phasenstabilität und Kriechverhalten unter kombinierter thermischer und Bestrahlungsbelastung vorherzusagen.

Methodik
Die Autoren verwenden ein speziesenspezifisches kinetisches Monte-Carlo-Framework (KMC), um die Vakanzendiffusion über einen systematischen Bereich von Fe/Cr-Zusammensetzungen (von Fe-reich bis Cr-reich) mit einer festen Konzentration von 5 at.% Al zu modellieren. Die Methodik umfasst die folgenden Schritte:

1. Potenzial- und Datengenerierung: Interatomare Wechselwirkungen werden mittels eines Embedded Atom Method (EAM)-Potenzials ausgewertet, das an eine Dichtefunktionaltheorie-Datenbank (DFT) angepasst wurde. Die Vakanzenmigrationsdaten werden mittels des „Hop-Decorate“-Workflows generiert, bei dem eine BCC-Fe-Referenzstruktur zufällig mit Fe-, Cr- und Al-Atomen entsprechend den Zielzusammensetzungen dekoriert wird. Die Pfade der minimalen Energie werden mittels der Climbing-Image Nudged Elastic Band (CI-NEB)-Methode für die globalen Zusammensetzungen Fe$_{80}$Cr$_{15}$Al$_5$, Fe$_{47,5}$Cr$_{47,5}$Al$_5$ und Fe$_{15}$Cr$_{80}$Al$_5$ berechnet.
2. Surrogatmodellierung: Um die Rechenkosten von On-the-fly-NEB-Berechnungen während der KMC-Simulation zu vermeiden, wurden lineare Surrogatmodelle auf Basis der generierten Barrieren-Datenbank trainiert. Die lokale chemische Umgebung für jeden Sprung wird durch Deskriptoren kodiert, welche die Anzahl der Fe-, Cr- und Al-Atome in festen radialen Schalen um den Sattelpunkt und das migrierende Atom zählen. Es wurden separate lineare Regressionsmodelle für jede migrierende Spezies (Fe, Cr, Al) angepasst, um die spezies-spezifische Abhängigkeit der Barrieren zu erfassen.
3. KMC-Simulationen: Die trainierten Surrogatmodelle steuern Residence-Time-KMC-Simulationen bei Temperaturen von 500 K bis 1500 K. Die Simulationen verfolgen einzelne Vakanzentrajektorien über $10^5$ Schritte, um die mittlere quadratische Verschiebung (MSD) zu berechnen und die Diffusionskoeffizienten ($D$), Aktivierungsenergien ($E_a$) und Präexponentiellen Faktoren ($D_0$) zu extrahieren.

Kernergebnisse

• Migrationsbarrieren und lokale Umgebung: Die Studie bestätigt, dass Migrationsbarrieren hochsensibel auf die lokale chemische Umgebung reagieren. Während reines Fe eine Barriere von ~0,63 eV aufweist, liegen Cr bei ~0,95 eV und Al bei ~0,69 eV. In Legierungen verbreitern sich die Verteilungen der Barrieren signifikant, insbesondere für Cr und Al, was darauf hindeutet, dass die lokale chemische Asymmetrie die Migrationspfade moduliert.
• Rolle von Aluminium: Trotz der geringen Konzentration (5 at.%) übt Al einen unverhältnismäßig großen Einfluss auf die Migrationsbarrieren aus. In allen drei Surrogatmodellen erhöht das Vorhandensein eines Al-Atoms in der ersten nächsten Nachbarschaft (1NN) der Sattelschale die Energiebarriere signifikant (+0,048 eV für Fe, +0,073 eV für Cr und +0,151 eV für Al).
• Zusammensetzungsabhängigkeit der Diffusivität:
  • Fe-reich bis äquiatomisch: In Zusammensetzungen im Bereich von Fe$_{80}$Cr$_{15}$Al$_5$ bis Fe$_{32,5}$Cr$_{62,5}$Al$_5$ ist die Vakanzendiffusivität relativ unempfindlich gegenüber der globalen Zusammensetzung. Die Aktivierungsenergien konvergieren in einen Bereich von 0,70–0,84 eV, was darauf hindeutet, dass die Fe-Migration die Transportprozesse in diesen Regimen dominiert.
  • Cr-reicher Bereich: Die Cr-reiche Zusammensetzung (Fe$_{15}$Cr$_{80}$Al$_5$) weist eine deutlich erhöhte Aktivierungsenergie von 1,135 eV auf. Folglich ist ihre makroskopische Diffusivität bei niedrigeren Temperaturen (z. B. ~$10^{-17}$ m$^2$ s$^{-1}$ bei ~570 K im Vergleich zu ~$10^{-13}$ m$^2$ s$^{-1}$ für Fe-reiche Legierungen) um Größenordnungen niedriger.
• Perkolation und „träge“ Diffusion: Die Ergebnisse legen nahe, dass es eine Perkolationsschwelle gibt, bei der das vernetzte Netzwerk niederenergetischer Pfade (dominiert durch Fe-Sprünge) mit zunehmendem Cr-Gehalt zusammenbricht. Im Cr-reichen Regime wird die Vakanzenmigration durch hochbarrierearme Cr-Sprünge dominiert, was zu einer „trägen Diffusion“ führt. Die Autoren beobachten, dass Vakanzentrajektorien in chemisch rauen Landschaften aus Trapping-Taschen entkommen müssen, was zu einem nicht-linearen Verhalten der mittleren quadratischen Verschiebung führt, bevor ein diffuser Zustand erreicht wird.

Bedeutung und Ansprüche
Die Arbeit behauptet, dass die makroskopische Diffusivität von FeCrAl-Legierungen nicht durch eine einfache zusammensetzungsgewichtete Mittelung der Barrieren der reinen Endglieder vorhergesagt werden kann; vielmehr wird sie durch die Heterogenität der lokalen Energielandschaft und die spezifische Konnektivität der niederenergetischen Migrationspfade bestimmt.

• Design-Implikationen: Die Studie identifiziert ein distinktes Zusammensetzungsregime (Cr-reich), in dem die Vakanzenmobilität unterdrückt wird. Diese Unterdrückung impliziert, dass in Cr-reichen $\alpha'$-Phasen (die unter Bestrahlung bekanntlich entstehen) Vakanzen weniger mobil sind, was potenziell die Rate der Vakanz-Interstitials-Rekombination reduziert und die Defektakkumulationskinetik verändert.
• Unfalltoleranz: Unter LOCA-Bedingungen (Loss-of-Coolant Accident), bei denen die Temperaturen über 1200 K steigen, wird die Abhängigkeit der Diffusivität von der Zusammensetzung deutlich. Die Autoren legen nahe, dass die Cr-reiche Zusammensetzung mit ihrer hohen Aktivierungsenergie und unterdrückten Diffusivität eine vorteilhafte Vakanenkinetik für unfalltolerante Brennelementumhüllungen bieten kann, was deren bekannte Oxidationsbeständigkeit ergänzen könnte.
• Methodischer Beitrag: Die Arbeit zeigt, dass reduzierte kinetische Monte-Carlo-Frameworks, die durch spezies-aufgelöste lineare Surrogatmodelle auf Basis von NEB-Datenbanken parametrisiert sind, effektiv den Defekstransport über komplexe Zusammensetzungsräume abbilden können. Dieser Ansatz stellt eine mechanistische Verbindung zwischen atomistischen Prozessen und mesoskaligen Modellen her, die in Brennelement-Leistungs-Codes verwendet werden, und geht damit über empirisches Trial-and-Error beim Legierungsdesign hinaus.

Vom Autor angemerktes Limit
Die Autoren merken explizit an, dass das Modell das Gitter als statisch und starr behandelt und somit lokale Relaxation sowie die zeitliche Entwicklung der strahlungsinduzierten Zusammensetzungsänderung (Segregation) ignoriert. Das lineare Surrogatmodell kann keine Richtungs- oder Symmetriebrechungseffekte oder magnetische Beiträge zur Energielandschaft erfassen. Darüber hinaus konzentriert sich die Studie ausschließlich auf den Vakanzentransport und adresst nicht den Transport von Interstitial-Clustern oder die strahlungsinduzierte Transformation von Versetzungsschleifen, welche ebenfalls die Senkendichte beeinflussen könnten.