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

Sintesi Tecnica: Diffusione di Vacanze attraverso lo Spazio Composizionale delle Leghe FeCrAl per il Rivestimento del Combustibile Tollerante agli Accidenti

Problematica
Le leghe ferro-cromo-alluminio (FeCrAl) sono tra i principali candidati per il rivestimento del combustibile tollerante agli accidenti (ATF) nei reattori ad acqua leggera (LWR), grazie alla loro superiore resistenza all'ossidazione ad alta temperatura rispetto allo Zircaloy. Tuttavia, la vita operativa in reattore di questi rivestimenti è governata dall'evoluzione microstrutturale indotta dalle radiazioni, dove il trasporto dei difetti puntiformi — specificamente la diffusione delle vacanze — rappresenta il meccanismo dominante. Mentre la formazione di strati protettivi di ossido è ben caratterizzata, la fisica fondamentale dei difetti nei FeCrAl sotto irraggiamento rimane scarsamente compresa. I trattamenti teorici esistenti si affidano spesso ad ambienti idealizzati o composizionalmente semplificati, fallendo nel rendere conto delle complesse fluttuazioni chimiche locali inerenti alle leghe concentrate a soluzione solida. Questi ambienti locali alterano significativamente il legame delle vacanze, i percorsi di migrazione e i comportamenti collettivi come l'agglomerazione (clustering) e il trapping. Esiste una critica necessità di distinguere i contributi individuali dei costituenti (Fe, Cr, Al) al trasporto complessivo delle vacanze per prevedere l'evoluzione a lungo termine del materiale, la stabilità di fase e il creep sotto carichi termici e di irraggiamento combinati.

Metodologia
Gli autori impiegano un framework di cinetica Monte Carlo (KMC) risolto per specie per modellare la diffusione delle vacanze attraverso un intervallo sistematico di composizioni Fe/Cr (da ricche in Fe a ricche in Cr) con una concentrazione fissa di Al al 5 at.%. La metodologia procede attraverso le seguenti fasi:

1. Generazione di Potenziali e Dati: Le interazioni interatomiche sono valutate utilizzando un potenziale dell'Atomo Embeddato (EAM) adattato a un database di Teoria del Funzionale della Densità (DFT). I dati di migrazione delle vacanze sono generati tramite il workflow "Hop-Decorate", in cui una struttura di riferimento di Fe BCC viene decorata casualmente con atomi di Fe, Cr e Al secondo le composizioni globali target. I percorsi a energia minima sono calcolati utilizzando il metodo Climbing-Image Nched Elastic Band (CI-NEB) per le composizioni globali di $\text{Fe}_{80}\text{Cr}_{15}\text{Al}_5$, $\text{Fe}_{47.5}\text{Cr}_{47.5}\text{Al}_5$ e $\text{Fe}_{15}\text{Cr}_{80}\text{Al}_5$.
2. Modellazione Surrogata: Per evitare il costo computazionale dei calcoli NEB on-the-fly durante la KMC, vengono addestrati modelli surrogati lineari sulla base di dati delle barriere generata. L'ambiente chimico locale per ogni salto (hop) è codificato tramite descrittori che contano gli atomi di Fe, Cr e Al all'interno di gusci radiali fissi attorno al punto di sella e all'atomo migrante. Modelli di regressione lineare separati sono adattati per ogni specie migrante (Fe, Cr, Al) per catturare la dipendenza delle barriere specifica per specie.
3. Simulazioni KMC: I modelli surrogati addestrati guidano simulazioni KMC a tempo di residenza a temperature comprese tra 500 K e 1500 K. Le simulazioni tracciano le traiettorie delle singole vacanze su $10^5$ passi per calcolare lo spostamento quadratico medio (MSD) ed estrarre i coefficienti di diffusione ($D$), le energie di attivazione ($E_a$) e i fattori pre-esponenziali ($D_0$).

Risultati Chiave

• Barriere di Migrazione e Ambiente Locale: Lo studio conferma che le barriere di migrazione sono altamente sensibili all'ambiente chimico locale. Mentre il Fe puro presenta una barriera di ~0,63 eV, il Cr è di ~0,95 eV e l'Al è di ~0,69 eV. Nelle leghe, le distribuzioni delle barriere si allargano significamente, in particolare per Cr e Al, indicando che l'asimmetria chimica locale modula i percorsi di migrazione.
• Ruolo dell'Alluminio: Nonostante la bassa concentrazione (5 at.%), l'Al esercita un'influenza sproporzionata sulle barriere di migrazione. In tutti e tre i modelli surrogati, la presenza di un atomo di Al nel primo guscio di vicini più prossimi (1NN) della sella aumenta significativamente la barriera energetica (+0,048 eV per Fe, +0,073 eV per Cr e +0,151 eV per Al).
• Dipendenza Composizionale della Diffusività:
  • Da Fe-ricco a Equiatomico: Nelle composizioni che vanno da $\text{Fe}_{80}\text{Cr}_{15}\text{Al}_5$ a $\text{Fe}_{32.5}\text{Cr}_{62.5}\text{Al}_5$, la diffusività delle vacanze è relativamente insensibile alla composizione globale. Le energie di attivazione convergono in un intervio di 0,70–0,84 eV, suggerendo che la migrazione del Fe domina il trasporto in questi regimi.
  • Regime Ricco in Cr: La composizione ricca in Cr ($\text{Fe}_{15}\text{Cr}_{80}\text{Al}_5$) presenta un'energia di attivazione marcatamente elevata di 1,135 eV. Di conseguenza, la sua diffusività macroscopica è di ordini di grandezza inferiore rispetto alle leghe ricche in Fe a temperature inferiori (ad esempio, ~$10^{-17}$ m$^2$ s$^{-1}$ a ~570 K rispetto a ~$10^{-13}$ m$^2$ s$^{-1}$ per le leghe ricche in Fe).
• Percolazione e Diffusione "Sluggish": I risultati suggeriscono l'esistenza di una soglia di percolazione dove la rete connessa di percorsi a bassa barriera (dominata dai salti del Fe) si interrompe all'aumentare del contenuto di Cr. Nel regime ricco in Cr, la migrazione delle vacanze è dominata da salti del Cr ad alta barriera, portando a una "diffusione lenta" (sluggish diffusion). Gli autori osservano che le traiettorie delle vacanze in paesaggi chimicamente rugosi comportano l'uscita da tasche di intrappolamento (trapping pockets), risultando in un comportamento non lineare dello spostamento quadratico medio prima di raggiungere un regime diffusivo.

Significatività e Rivendicazioni
Il documento afferma che la diffusività macroscopica delle leghe FeCrAl non può essere predetta da una semplice media ponderata della composizione delle barriere degli end-member puri; piuttosto, essa è determinata dall'eterogeneità del paesaggio energetico locale e dalla specifica connettività dei percorsi di migrazione a bassa barriera.

• Implicazioni per la Progettazione: Lo studio identifica un regime composizionale distinto (ricco in Cr) in cui la mobilità delle vacanze è soppressa. Questa soppressione implica che nelle fasi $\alpha'$ ricche in Cr (note per formarsi sotto irraggiamento), le vacanze sono meno mobili, il che può ridurre il tasso di ricombinazione vacanza-interstiziale e alterare la cinetica di accumulo dei difetti.
• Tolleranza agli Accidenti: In condizioni di perdita di refrigerante (LOCA), dove le temperature salgono sopra i 1200 K, la dipendenza composizionale della diffusività diventa pronunciata. Gli autori suggeriscono che la composizione ricca in Cr, con la sua alta energia di attivazione e la diffusività soppressa, possa offrire una cinetica delle vacanze favorevole per il rivestimento tollerante agli accidenti, potenzialmente complementando la sua nota resistenza all'ossidazione.
• Contributo Metodologico: Il lavoro dimosta che i framework di cinetica Monte Carlo a ordine ridotto, parametrizzati da modelli surrogati lineari risolti per specie e addestrati su database NEB, possono mappare efficacemente il trasporto di difetti attraverso complessi spazi composizionali. Questo approccio fornisce un collegamento meccanicistico tra i processi atomistici e i modelli mesoscopici utilizzati nei codici di prestazione del combustibile, andando oltre l'approccio empirico per tentativi ed errori nella progettazione delle leghe.

Limitazioni Note dagli Autori
Gli autori notano esplicitamente diverse limitazioni: il modello tratta il reticolo come statico e congelato, ignorando la rilassazione locale e l'evoluzione della composizione indotta dalle radiazioni (segregazione) nel tempo. Il modello surrogato lineare non può catturare effetti direzionali o di rottura della simmetria, né contributi magnetici al paesaggio energetico. Inoltre, lo studio si concentra esclusivamente sul trasporto delle vacanze e non affronta il trasporto dei cluster interstiziali o la trasformazione dei loop di dislocazione indotta dalle radiazioni, che possono anch'essi influenzare le densità dei sink di difetti.