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

Resumen Técnico: Difusión de Vacantes a través del Espacio de Composición de Aleaciones de FeCrAl para Revestimientos de Combustible Tolerantes a Accidentes

Planteamiento del Problema
Las aleaciones de hierro-cromo-aluminio (FeCrAl) son candidatas principales para los revestimientos de combustible tolerantes a accidentes (ATF) en reactores de agua ligera (LWR) debido a su superior resistencia a la oxidación a alta temperatura en comparación con el Zircaloy. Sin embargo, la vida útil en reactor de estos revestimientos está gobernada por la evolución microestructural inducida por la radiación, donde el transporte de defectos puntuales —específicamente la difusión de vacantes— es el mecanismo dominante. Aunque la formación de capas de óxido protectoras está bien caracterizada, la física fundamental de los defectos bajo irradiación en el FeCrAl sigue estando insuficientemente comprendida. Los tratamientos teóricos existentes suelen apoyarse en entornos idealizados o composicionalmente simplificados, fallando al no considerar las complejas fluctuaciones químicas locales inherentes a las aleaciones concentradas de solución sólida. Estos entornos locales alteran significativamente la unión de vacantes, las rutas de migración y los comportamientos colectivos como la formación de cúmulos (clustering) y el atrapamiento. Existe una necesidad crítica de desentrañar las contribuciones individuales de los elementos constituyentes (Fe, Cr, Al) al transporte de vacantes global para predecir la evolución del material a largo plazo, la estabilidad de fase y la fluencia (creep) bajo cargas térmicas y de irradiación combinadas.

Metodología
Los autores emplean un marco de Monte Carlo cinético (KMC) resuelto por especie para modelar la difusión de vacantes a través de un rango sistemático de composiciones de Fe/Cr (desde ricas en Fe hasta ricas en Cr) con una concentración fija de 5 at.% de Al. La metodología procede a través de los siguientes pasos:

1. Generación de Potenciales y Datos: Las interacciones interatómicas se evalúan mediante un potencial de Método de Átomo Embebido (EAM) ajustado a una base de datos de la Teoría del Funcional de la Densidad (DFT). Los datos de migración de vacantes se generan utilizando el flujo de trabajo "Hop-Decorate", donde una estructura de referencia de Fe BCC se decora aleatoriamente con átomos de Fe, Cr y Al de acuerdo con las composiciones globales objetivo. Se computan las rutas de mínima energía mediante el método de Imagen Trepadora de Nudged Elastic Band (CI-NEB) para las composiciones globales de $\text{Fe}_{80}\text{Cr}_{15}\text{Al}_5$, $\text{Fe}_{47.5}\text{Cr}_{47.5}\text{Al}_5$ y $\text{Fe}_{15}\text{Cr}_{80}\text{Al}_5$.
2. Modelado Sustituto (Surrogate Modeling): Para evitar el costo computacional de los cálculos de NEB en tiempo real durante el KMC, se entrenan modelos sustitutos lineales en la base de datos de barreras generada. El entorno químico local para cada salto se codifica utilizando descriptores que cuentan átomos de Fe, Cr y Al dentro de capas radiales fijas alrededor del punto de silla y el átomo migrante. Se ajustan modelos de regresión lineal independientes para cada especie migrante (Fe, Cr, Al) para capturar la dependencia de la barrera específica de la especie.
3. Simulaciones KMC: Los modelos sustitutos entrenados impulsan simulaciones de KMC de tiempo de residencia a temperaturas que oscilan entre 500 K y 1500 K. Las simulaciones rastrean trayectorias de vacantes individuales a lo largo de $10^5$ pasos para calcular el desplazamiento cuadrático medio (MSD) y extraer los coeficientes de difusión ($D$), las energías de activación ($E_a$) y los factores pre-exponenciales ($D_0$).

Resultados Clave

• Barreras de Migración y Entorno Local: El estudio confirma que las barreras de migración son altamente sensibles al entorno químico local. Mientras que el Fe puro exhibe una barrera de ~0.63 eV, el Cr es de ~0.95 eV y el Al es de ~0.69 eV. En las aleaciones, las distribuciones de las barreras se ensanchan significativamente, particularmente para Cr y Al, indicando que la asimetría química local modula las rutas de migración.
• Papel del Aluminio: A pesar de su baja concentración (5 at.%), el Al ejerce una influencia desproporcionada en las barreras de migración. En los tres modelos sustitutos, la presencia de un átomo de Al en la primera capa de primer vecino (1NN) del punto de silla eleva significativamente la barrera de energía (+0.048 eV para Fe, +0.073 eV para Cr y +0.151 eV para Al).
• Dependencia Composicional de la Difusividad:
  • De Rico en Fe a Equiatómico: En composiciones que van desde $\text{Fe}_{80}\text{Cr}_{15}\text{Al}_5$ hasta $\text{Fe}_{32.5}\text{Cr}_{62.5}\text{Al}_5$, la difusividad de las vacantes es relativamente insensible a la composición global. Las energías de activación convergen a un rango de 0.70–0.84 eV, sugiriendo que la migración de Fe domina el transporte en estos regímenes.
  • Régimen Rico en Cr: La composición rica en Cr ($\text{Fe}_{15}\text{Cr}_{80}\text{Al}_5$) exhibe una energía de activación notablemente elevada de 1.135 eV. En consecuencia, su difusividad macroscópica es órdenes de magnitud inferior a la de las aleaciones ricas en Fe a temperaturas bajas (por ejemplo, ~$10^{-17}$ m$^2$ s$^{-1}$ a ~570 K frente a ~$10^{-13}$ m$^2$ s$^{-1}$ para las aleaciones ricas en Fe).
• Percolación y Difusión "Lenta" (Sluggish): Los resultados sugieren un umbral de percolación donde la red conectada de rutas de baja barrera (dominada por saltos de Fe) se rompe a medida que aumenta el contenido de Cr. En el régimen rico en Cr, la migración de vacantes está dominada por saltos de Cr de alta barrera, lo que conduce a una "difusión lenta". Los autores observan que las trayectorias de las vacantes en paisajes químicamente rugosos implican la salida de bolsas de atrapamiento, resultando en un comportamiento de desplazamiento cuadrático medio no lineal antes de alcanzar un régimen difusivo.

Significación y Reivindicaciones
El artículo sostiene que la difusividad macroscópica de las aleaciones de FeCrAl no puede predecirse mediante un simple promedio ponderado por composición de las barreras de los extremos puros; más bien, está determinada por la heterogeneidad del paisaje de energía local y la conectividad específica de las rutas de migración de baja barrera.

• Implicaciones de Diseño: El estudio identifica un régimen composicional distinto (rico en Cr) donde la movilidad de las vacantes se ve suprimida. Esta supresión implica que en las fases $\alpha'$ ricas en Cr (conocidas por formarse bajo irradiación), las vacantes son menos móviles, lo que potencialmente reduce la tasa de recombinación vacante-intersticial y altera la cinética de acumulación de defectos.
• Tolerancia a Accidentes: Bajo condiciones de pérdida de refrigerante (LOCA), donde las temperaturas aumentan por encima de 1200 K, la dependencia composicional de la difusividad se vuelve pronunciada. Los autores sugieren que la composición rica en Cr, con su alta energía de activación y difusividad suprimida, puede ofrecer una cinética de vacantes favorable para el revestimiento tolerante a accidentes, complementando potencialmente su conocida resistencia a la oxidación.
• Contribución Metodológica: El trabajo demuestra que los marcos de Monte Carlo cinético de orden reducido, parametrizados por modelos sustitutos lineales resueltos por especie y entrenados en bases de datos NEB, pueden mapear eficazmente el transporte de defectos a través de complejos espacios de composición. Este enfoque proporciona un vínculo mecanístico entre los procesos atómicos y los modelos de mesoescala utilizados en los códigos de rendimiento de combustible, yendo más allá del ensayo y error empírico en el diseño de aleaciones.

Limitaciones Señaladas por los Autores
Los autores señalan explícamente varias limitaciones: el modelo trata la red como estática y congelada, ignorando la relajación local y la evolución de la composición inducida por la radiación (segregación) a lo largo del tiempo. El modelo sustituto lineal no puede capturar efectos direccionales o de ruptura de simetría, ni contribuciones magnéticas al paisaje de energía. Además, el estudio se centra únicamente en el transporte de vacantes y no aborda el transporte de cúmulos intersticiales ni la transformación de los bucles de dislocación inducida por la radiación, que también pueden influir en las densidades de sumideros de defectos.