Pressure-dependent melting and crystallization of B2-NiAl from neural-network molecular dynamics

Resumo Técnico: Fusão e Cristalização de B2-NiAl Dependente de Pressão

Definição do Problema
O alumineto de níquel (NiAl) é um composto intermetálico ordenado prototípico com aplicações estruturais significativas em altas temperaturas. Embora seu comportamento de fusão à pressão ambiente seja bem caracterizado como congruente, seu comportamento de fusão sob alta pressão permanece insuficientemente determinado. Estudos existentes de dinâmica molecular (MD) basearam-se primariamente em potenciais do método de átomo incorporado (EAM), que podem carecer da transferibilidade necessária para uma termodinâmica de alta pressão precisa. Além disso, os dados experimentais em pressões elevadas são limitados e exibem dispersão devido às dificuldades de medição. Uma questão crítica em aberto é se a fusão do NiAl B2 estequiométrico permanece congruente (onde as fases sólida e líquida compartilham a mesma composição) sob compressão, ou se ocorre partição química. Este estudo aborda essas lacunas empregando potenciais interatômicos de aprendizado de máquina de alta precisão para determinar a curva de fusão dependente da pressão ($T_m(P)$), verificar a consistência termodinâmica e caracterizar a congruência da fusão e a subsequente cristalização.

Metodologia
Os autores utilizaram dinâmica molecular de rede neural (NNMD) com um modelo interatômico de Potencial Profundo (Deep Potential - DeepMD). A metodologia compreendeu quatro componentes principais:

1. Desenvolvimento do Potencial: Um modelo de Potencial Profundo foi treinado em um banco de dados estendido de aproximadamente 37.200 configurações. Este banco de dados foi construído combinando dados de primeiros princípios (DFT) para ligas líquidas de Al–Cu–Ni com estruturas ordenadas de Al–Ni, cristais distorcidos e superfícies. Para garantir a cobertura de estados de alta pressão relevantes para o estudo, um loop de aprendizado ativo foi empregado usando o modelo PET-MAD-S como gerador de configurações. As configurações foram selecionadas via descritores de Sobreposição Suave de Posições Atômicas (SOAP) e Amostragem de Ponto Mais Distante (FPS), sendo então rotuladas com DFT. O modelo final foi treinado para reproduzir energias, forças e tensores de tensão, alcançando alta precisão (RMSE de força < 0,16 eV/Å) na faixa de 0 a 50 GPa.
2. Determinação da Curva de Fusão: As temperaturas de fusão foram determinadas utilizando o método de coexistência de duas fases sólido-líquido. As simulações foram realizadas no ensemble de pressão constante e entalpia constante (NPH) para supercélulas contendo 16.384 átomos. Uma interface sólido-líquido foi criada e o sistema foi evoluído até que um estado de coexistência estável fosse alcançado, permitindo a determinação de $T_m$ em pressões de 0 a 50 GPa.
3. Verificação de Consistência Termodinâmica: Para validar a linha de fusão calculada, simulações independentes de fase única do sólido e do líquido foram conduzidas na $T_m(P)$ determinada. As mudanças de entalpia ($\Delta h$) e volume ($\Delta v$) na fusão foram calculadas. Esses valores foram utilizados para calcular a inclinação de Clapeyron ($dT_m/dP = T_m \Delta v / \Delta h$) e comparados contra a derivada analítica da curva ajustada de $T_m(P)$.
4. Análise de Congruência e Estrutural:
   • Composição: Perfis de composição resolvidos por camada ao longo da direção [001] foram analisados para medir a fração de alumínio na fase líquida em relação ao valor estequiométrico (0,5).
   • Cristalização: Simulações de resfriamento rápido (melt-quench) foram realizadas em várias taxas de resfriamento e pressões. A análise estrutural utilizou o Emparelhamento de Template Poliedral (PTM) para identificar geometrias locais (BCC, FCC, HCP, etc.) e um parâmetro de ordem local específico ($q_{B2}$) para distinguir estruturas B2 quimicamente ordenadas de geometrias do tipo BCC desordenadas.

Principais Resultados

• Curva de Fusão ($T_m(P)$): A temperatura de fusão do B2-NiAl aumenta monotonicamente com a pressão, subindo de $1896 \pm 9$ K a 0 GPa para $3569 \pm 17$ K a 50 GPa. Os resultados concordam estreitamente com os valores experimentais de pressão ambiente ($\sim$1911 K) e com cálculos anteriores baseados em EAM até ~20 GPa. No entanto, em pressões mais altas (30–50 GPa), o modelo DeepMD prevê um aumento significativamente mais íngreme em $T_m$ comparado aos resultados anteriores de EAM, sugerindo uma descrição diferente da estabilidade termodinâmica relativa das fases sólida e líquida comprimidas.
• Consistência Termodinâmica: As inclinações de Clapeyron calculadas, derivadas de medições independentes de $\Delta h$ e $\Delta v$, mostram excelente concordância com a derivada da curva de fusão de Simon–Glatzel ajustada em quase toda a faixa de pressão. Isso confirma a consistência termodinâmica interna das simulações de coexistência. A entalpia de fusão aumenta com a pressão, enquanto a mudança relativa de volume ($\Delta v/v_{sol}$) diminui, indicando um contraste de densidade reduzido entre as fases sob compressão.
• Fusão Congruente: A análise de composição resolvida por camada revela que a fase líquida coexistente permanece essencialmente equiatômica (fração de Al $\approx$ 0,5) em todas as pressões investigadas. Os desvios do valor estequiométrico não excedem $5 \times 10^{-3}$. Isso fornece evidência atomística direta de que a fusão do B2-NiAl permanece congruente dentro do presente modelo, sem partição química significativa entre as fases sólida e líquida.
• Comportamento de Cristalização: Simulações de melt-quench demonstram que o líquido super-resfriado cristaliza em um estado policristalino dominado por uma estrutura geométrica do tipo BCC. O ordenamento químico B2 desenvolve-se progressivamente dentro desta estrutura. A tendência de cristalização aumenta com taxas de resfriamento mais lentas e pressões mais altas. A 30 GPa, a fração BCC é substancialmente maior do que a 0 GPa para as mesmas taxas de resfriamento, indicando que a compressão promove a cristalização em estruturas derivadas de bcc.

Significância e Alegações
O artigo estabelece uma descrição termodinamicamente consistente e dependente da pressão da fusão do B2-NiAl utilizando um modelo de Potencial Profundo de alta fidelidade. As principais contribuições são:

1. Curva de Fusão Refinada: Fornecendo uma relação $T_m(P)$ mais precisa que diverge das previsões anteriores de EAM em altas pressões, sugerindo a necessidade de potenciais mais precisos para capturar a termodinâmica de alta pressão.
2. Validação da Congruência: Oferecendo evidência atomística direta de que a fusão do B2-NiAl permanece congruente sob compressão, resolvendo a incerteza sobre uma potencial partição química.
3. Rigor Metodológico: Demonstrando a utilidade da combinação de simulações de coexistência de duas fases com verificações termodinâmicas independentes de fase única (relação de Clapeyron) para garantir a confiabilidade das previsões do potencial de aprendizado de máquina.
4. Insight Estrutural: Esclarecendo a via de cristalização, mostrando que uma estrutura geométrica do tipo bcc se forma antes do estabelecimento do ordenamento químico B2 completo, e que a pressão aumenta essa tendência de cristalização.

Os autores concluem que estes resultados fornecem uma descrição atomística coerente da fusão induzida por pressão no B2-NiAl, esclarecendo tanto a linha de fusão quanto o caráter da fusão sob compressão dentro do escopo do modelo de potencial desenvolvido.