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

Resumo Técnico: Difusão de Lacunas através do Espaço de Composição da Liga FeCrAl para Revestimento de Combustível Tolerante a Acidentes

Definição do Problema
As ligas de ferro-cromo-alumínio (FeCrAl) são candidatas proeminentes para revestimentos de combustível tolerantes a acidentes (ATF) em reatores de água leve (LWR) devido à sua resistência superior à oxidação em altas temperaturas em comparação ao Zircaloy. No entanto, a vida útil em reator desses revestimentos é governada pela evolução microestrutural induzida por radiação, onde o transporte de defeitos pontuais — especificamente a difusão de lacunas — é o mecanismo dominante. Embora a formação de camadas de óxido protetoras seja bem caracterizada, a física fundamental dos defeitos em FeCrAl sob irradiação permanece pouco compreendida. Tratamentos teóricos existentes frequentemente dependem de ambientes idealizados ou composicionalmente simplificados, falhando em considerar as complexas flutuações químicas locais inerentes às ligas de solução sólida concentrada. Esses ambientes locais alteram significativamente a ligação de lacunas, as vias de migração e os comportamentos coletivos, tais como o agrupamento (clustering) e o aprisionamento (trapping). Há uma necessidade crítica de distinguir as contribuições individuais dos elementos constituintes (Fe, Cr, Al) para o transporte de lacunas global, a fim de prever a evolução de longo prazo do material, a estabilidade de fase e a fluência (creep) sob carga térmica e de irradiação combinadas.

Metodologia
Os autores empregam uma estrutura de Monte Carlo cinético (KMC) resolvida por espécie para modelar a difusão de lacunas através de uma gama sistemática de composições de Fe/Cr (de rico em Fe a rico em Cr) com uma concentração fixa de 5 at.% de Al. A metodologia procede através das seguintes etapas:

1. Potencial e Geração de Dados: As interações interatômicas são avaliadas usando um potencial de Método de Átomo Embutido (EAM) ajustado a um banco de dados de Teoria do Funcional da Densidade (DFT). Os dados de migração de lacunas são gerados usando o fluxo de trabalho "Hop-Decorate", onde uma estrutura de referência de Fe BCC é decorada aleatoriamente com átomos de Fe, Cr e Al de acordo com as composições globais alvo. Os caminhos de energia mínima são computados usando o método de Imagem Escalante Nudged Elastic Band (CI-NEB) para as composições globais de Fe$_{80}$Cr$_{15}$Al$_5$, Fe$_{47,5}$Cr$_{47,5}$Al$_5$ e Fe$_{15}$Cr$_{80}$Al$_5$.
2. Modelagem Substituta (Surrogate Modeling): Para evitar o custo computacional de cálculos NEB em tempo real durante o KMC, modelos substitutos lineares são treinados no banco de dados de barreiras gerado. O ambiente químico local para cada salto é codificado usando descritores que contam átomos de Fe, Cr e Al em camadas radiais fixas ao redor do ponto de sela e do átomo migrante. Modelos de regressão linear separados são ajustados para cada espécie migrante (Fe, Cr, Al) para capturar as dependências de barreira específicas da espécie.
3. Simulações KMC: Os modelos substitutos treinados direcionam simulações de KMC de tempo de residência em temperaturas variando de 500 K a 1500 K. As simulações rastreiam trajetórias de lacunas únicas ao longo de $10^5$ passos para calcular o deslocamento quadrático médio (MSD) e extrair coeficientes de difusão ($D$), energias de ativação ($E_a$) e fatores pré-exponenciais ($D_0$).

Resultos Principais

• Barreiras de Migração e Ambiente Local: O estudo confirma que as barreiras de migração são altamente sensíveis ao ambiente químico local. Enquanto o Fe puro exibe uma barreira de ~0,63 eV, o Cr é ~0,95 eV e o Al é ~0,69 eV. Nas ligas, as distribuições de barreiras alargam-se significativamente, particularmente para Cr e Al, indicando que a assimetria química local modula as vias de migração.
• Papel do Alumínio: Apesar da sua baixa concentração (5 at.%), o Al exerce uma influência desproporcional nas barreiras de migração. Em todos os três modelos substitutos, a presença de um átomo de Al na primeira camada de vizinhos mais próximos (1NN) eleva significativamente a barreira de energia (+0,048 eV para Fe, +0,073 eV para Cr e +0,151 eV para Al).
• Dependência Composicional da Difusividade:
  • Rico em Fe a Equiatômico: Em composições variando de Fe$_{80}$Cr$_{15}$Al$_5$ a Fe$_{32,5}$Cr$_{62,5}$Al$_5$, a difusividade de lacunas é relativamente insensível à composição global. As energias de ativação convergem para uma faixa de 0,70–0,84 eV, sugerindo que a migração de Fe domina o transporte nestes regimes.
  • Regime Rico em Cr: A composição rica em Cr (Fe$_{15}$Cr$_{80}$Al$_5$) exibe uma energia de ativação marcadamente elevada de 1,135 eV. Consequentemente, sua difusividade macroscópica é ordens de magnitude inferior à das ligas ricas em Fe em temperaturas mais baixas (ex: ~$10^{-17}$ m$^2$ s$^{-1}$ a ~570 K comparado a ~$10^{-13}$ m$^2$ s$^{-1}$ para ligas ricas em Fe).
• Percolação e Difusão "Lenta" (Sluggish): Os resultados sugerem a existência de um limiar de percolação onde a rede conectada de caminhos de baixa barreira (dominada por saltos de Fe) se rompe conforme o conteúdo de Cr aumenta. No regime rico em Cr, a migração de lacunas é dominada por saltos de Cr de alta barreira, levando a uma "difusão lenta". Os autores observam que as trajetórias de lacunas em paisagens quimicamente rugosas envolvem a fuga de bolsões de aprisionamento, resultando em um comportamento não linear do deslocamento quadrático médio antes de atingir um regime difusivo.

Significância e Alegações
O artigo afirma que a difusividade macroscópica das ligas FeCrAl não pode ser prevista por uma simples média ponderada pela composição das barreiras dos componentes puros; em vez disso, ela é determinada pela heterogeneidade do cenário de energia local e pela conectividade específica das vias de migração de baixa barreira.

• Implicações de Design: O estudo identifica um regime composicional distinto (rico em Cr) onde a mobilidade de lacunas é suprimida. Esta supressão implica que, nas fases $\alpha'$ ricas em Cr (conhecidas por se formar sob irradiação), as lacunas são menos móveis, o que pode reduzir a taxa de recombinação lacuna-intersticial e alterar a cinética de acumulação de defeitos.
• Tolerância a Acidentes: Sob condições de Acidente de Perda de Fluido de Resfriamento (LOCA), onde as temperaturas sobem acima de 1200 K, a dependência composicional da difusividade torna-se pronunciada. Os autores sugerem que a composição rica em Cr, com sua alta energia de ativação e difusividade suprimida, pode oferecer uma cinética de lacunas favorável para o revestimento tolerante a acidentes, potencialmente complementando sua conhecida resistência à oxidação.
• Contribuição Metodológica: O trabalho demonstra que estruturas de Monte Carlo cinético de ordem reduzida, parametrizadas por modelos substitutos lineares resolvidos por espécie e treinados em bases de dados NEB, podem mapear efetivamente o transporte de defeitos através de complexos espaços composicionais. Esta abordagem fornece um elo mecanístico entre processos atômicos e modelos de mesoescala usados em códigos de desempenho de combustível, indo além do método empírico de tentativa e erro no design de ligas.

Limitações Notadas pelos Autores
Os autores explicitamente notam várias limitações: o modelo trata a rede como estática e congelada, ignorando a relaxação local e a evolução da composição induzida por radiação (segregação) ao longo do tempo. O modelo substituto linear não consegue capturar efeitos direcionais ou de quebra de simetria, nem contribuições magnéticas ao cenário de energia. Além disso, o estudo foca exclusivamente no transporte de lacunas e não aborda o transporte de clusters de intersticiais ou a transformação de loops de deslocamento induzida por radiação, que também podem influenciar as densidades de sumidouros de defeitos.