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

Sintesi Tecnica: Modellazione a Ordine Ridotto del Trasporto di Difetti mediante Cinetica Surrogata in $U_xPu_{1-x}N$

Definizione del Problema
Modellare il trasporto di difetti in materiali chimicamente disordinati, come i nitruri misti di attinidi ($U_xPu_{1-x}N$), presenta una sfida significativa. Nei reticoli chimicamente uniformi, i percorsi di migrazione dei difetti possono essere enumerati come un insieme finito di tassi di transizione. Tuttavia, nei sistemi disordinati, l'ambiente chimico locale varia continuamente, creando una distribuzione di barriere di migrazione piuttosto che un insieme chiuso. Questa complessità chimica rende la piena enumerazione degli ambienti tramite la Teoria del Funzionale della Densità (DFT) o la Dinamica Molecolare (MD) computazionalmente proibitiva. Inoltre, i modelli standard di catena di Markov a tempo continuo (CTMC) non possono accomodare la varietà effettivamente infinita di tassi di transizione generati dal disordine locale. Questo problema è critico per i combustibili nucleari avanzati come il nitruro di uranio (UN), che può evolvere in soluzioni solide di $U_xPu_{1-x}N$ a causa della fissione o può essere fabbricato come combustibile misto. Predire le prestazioni del combustibile (rigonfiamento, rilascio di gas di fissione, degradazione della conducibilità termica) richiede la comprensione di come la variazione composizionale alteri l'energetica e il trasporto dei difetti, eppure la simulazione atomistica diretta attraverso l'immenso spazio configurazionale è impraticabile.

Metodologia
Gli autori propongono un flusso di lavoro a ordine ridotto per incorporare i meccanismi atomici in forme funzionali computazionalmente efficienti. La metodologia consiste in quattro fasi principali:

1. Generazione dei Dati (Hop-Decorate): Utilizzando il workflow Hop-Decorate (HopDec) basato su LAMMPS e ASE, gli autori hanno generato un database di circa 3.000 barriere di migrazione ($E_b$) e differenze di energia ($\Delta E$). Ciò è stato ottenuto identificando percorsi di migrazione prototipici del primo vicino (1NN) in strutture di riferimento di UN e "ri-decorandoli" 1.000 volte a composizioni globali di $U_{0.5}Pu_{0.5}N$, $U_{0.1}Pu_{0.9}N$ e $U_{0.9}Pu_{0.1}N$. Le energie interatomiche sono state calcolate utilizzando un potenziale EAM Finnis–Sinclair combinato con un potenziale di Buckingham, validato rispetto ai dati DFT. Il metodo Climbing-Image Nudged Elastic Band (CI-NEB) è stato utilizzato per valutare l'energetica per ogni configurazione.
2. Costruzione del Modello Surrogato: L'ambiente chimico locale è stato definito dai conteggi di atomi di Uranio all'interno di specifici gusci di vicini più prossimi (NN) attorno al difetto nelle sue configurazioni iniziale, finale e di punto di sella. Funzioni surrogate lineari sono state adattate per predire $E_b$ e $\Delta E$ basandosi su questi semplici descrittori:
   • $\Delta E$ è stato modellato come una funzione lineare del cambiamento nei conteggi di U tra i gusci dei minimi.
   • $E_b$ è stato modellato come una funzione lineare dei conteggi di U nei gusci dei minimi e del punto di sella, più una correzione di base specifica per la specie ($\beta$).
   • Sono stati necessari adattamenti separati per i salti mediati da U e da Pu sul sottoreticolo degli attinidi a causa della forte dipendenza dalla specie.
3. Simulazione Kinetic Monte Carlo (KMC): I modelli surrogati sono stati integrati in un simulatore KMC basato su reticolo (laKMC) personalizzato. Le simulazioni sono state eseguite su supercelle di $10 \times 10 \times 10$ nm$^3$ (64.000 atomi) per 100.000 passi. Sono stati estratti spostamenti medi quadratici (MSD) per calcolare le diffusività a lungo termine ($D$) attraverso varie temperature e composizioni, con media statistica su 50 traiettorie per condizione.
4. Analisi del Trappolamento: Per chiarire le origini su scala atomica del trasporto macroscopico, gli autori hanno analizzato lo spazio cinetico-termodinamico calcolando i tempi di residenza medi mediati dal disordine ($\langle \tau \rangle_c$) e le variazioni di energia medie pesate per il tasso ($\langle \Delta E \rangle_c$) per ogni composizione.

Risultati Chiave
Lo studio ha analizzato tre specie di difetti: vacanze di attinide ($V_A$), vacanze di nitruro ($V_N$) e divacanze attinide-nitruro ($V_{A-N}$).

• Energetica di Migrazione: I modelli surrogati hanno raggiunto bassi errori quadratici medi (RMSE), indicando che l'energetica di migrazione in $U_xPu_{1-x}N$ è governata principalmente dall'identità chimica a corto raggio piuttosto che da complessi effetti collettivi a lungo raggio.

  • $V_A$ (Vacanza di Attinide): La migrazione è fortemente controllata dalla specie. La distribuzione delle barriere biforca in base al fatto che l'atomo migrante sia U o Pu. L'ambiente locale causa solo un lieve allargamento; l'effetto composizionale primario è la ponderazione statistica dei salti U rispetto a quelli Pu.
  • $V_N$ (Vacanza di Nitruro): La migrazione mostra una genuina modulazione dipendente dall'ambiente. La distribuzione delle barriere evolve continuamente con la composizione, mostrando caratteristiche multilobate a composizioni intermedie corrispondenti a distinte configurazioni locali di attinidi. In particolare, la presenza di una piccola quantità di U in una matrice ricca di Pu può creare percorsi favorevoli locali con barriere basse.
  • $V_{A-N}$ (Divacanza): Il difetto accoppiato mostra un comportamento simile alla somma delle vacanze isolate, con la migrazione che procede attraverso uno dei due sottoreticoli. La presenza della coppia induce uno spostamento sistematico nell'altezza della barriera ma non altera fondamentalmente la sensibilità alla coordinazione locale.

• Tendenze della Diffusività:

  • Non-linearità: La diffusività non segue una semplice interpolazione lineare (legge di Vegard) tra i membri finali UN e PuN.
  • Comportamento di $V_A$: La diffusività mostra un aumento rapido e non lineare con piccole aggiunte di Pu. Ciò è guidato dalla percolazione di salti mediati da Pu che sono più veloci. Man mano che la concentrazione di Pu aumenta, il sistema passa da essere limitato dai lenti salti di U a una rete di rapidi salti di Pu.
  • Comportamento di $V_N$: La diffusività mostra una "diffusione lenta" a basse concentrazioni di Pu, dove il Pu agisce come una trappola cinetica, riducendo il trasporto netto nonostante l'esistenza di alcuni eventi a bassa barriera. Ad alte concentrazioni di Pu, la diffusività aumenta man mano che i percorsi a bassa barriera diventano più prevalenti e connessi.
  • Comportamento di $V_{A-N}$: La diffusività della divacanza riflette una competizione tra le sensibilità chimiche di entrambi i sottoreticoli, risultando in una dipendenza composizionale non lineare e complessa, distinta dalle vacanze isolate.

• Analisi del Trappolamento: Le mappe cinetico-termodinamiche ($\langle \tau \rangle_c$ vs. $\langle \Delta E \rangle_c$) hanno rivelato che:

  • Per $V_A$, l'aumento del contenuto di Pu riduce il tempo di residenza senza alterare significativamente il bias energetico, confermando la percolazione di percorsi veloci.
  • Per $V_N$, piccole concentrazioni di Pu aumentano il bias energetico ($\langle \Delta E \rangle_c$) mentre il tempo di residenza rimane stabile, indicando la formazione di trappole chimiche che modulano le energie relative piuttosto che aumentare semplicemente le barriere.

Significato e Rivendicazioni
L'articolo sostiene di fornire un quadro generale e computazionalmente efficiente per modellare il trasporto di difetti in materiali chimicamente disordinati senza ricorrere a modelli di machine learning "black-box" privi di interpretabilità fisica.

• Rappresentazione a Ordine Ridotto: Gli autori dimostrano che la complessa energetica di migrazione dei nitruri di attinidi disordinati ammette una rappresentazione a ordine ridotto inaspettatamente semplice utilizzando funzioni lineari dei conteggi di coordinazione locale.
• Intuizione Meccanicistica: Il lavoro identifica due meccanismi dominanti che governano il trasporto: la cinetica controllata dalla specie (dominante per le vacanze di attinide) e il trappolamento dipendente dall'ambiente (dominante per le vacanze di nitruro). Evidenzia come la diffusività macroscopica emerga dalla competizione tra distribuzioni di barriere, connettività dei percorsi e tempi di residenza, piuttosto che semplicemente dai valori medi delle barriere.
• Scalabilità: Distillando i dati atomistici in funzioni surrogate compatte, l'approccio consente l'integrazione di meccanismi atomici in simulazioni a scala superiore (ad esempio, dinamica dei cluster o codici di prestazione del combustibile a scala continua). Ciò consente la modellazione efficiente delle composizioni che evolvono durante il burnup del combustibile o la trasmutazione.
• Applicabilità: Sebbene dimostrato su $U_xPu_{1-x}N$, il framework è presentato come applicabile ad altri materiali composizionalmente disordinati, come le leghe ad alta entropia e i semiconduttori drogati, dove l'enumerazione brute-force degli ambienti locali è impraticabile.

Gli autori segnalano alcune limitazioni, tra cui il vincolo delle divacanze a coppie di vicini più prossimi (ignorando la potenziale dissociazione/riassociazione) e l'uso di potenziali empirici che non trattano esplicitamente gli effetti elettronici (ad esempio, la localizzazione dell'orbitale 5f in Pu). Tuttavia, affermano che l'attuale approccio riesce a colmare il divario tra i meccanismi atomistici e le equazioni di trasporto trattabili per la modellazione predittiva delle prestazioni del combustibile.