Temperature-driven structural phase transitions in SmNiO$_3$: insights from deep potential molecular dynamics simulations

Resumo Técnico: Transições de Fase Estrutural Impulsionadas pela Temperatura em SmNiO3

Definição do Problema
Os nickelatos de terras raras (RNiO3) exibem uma transição metal-isolante (MIT) abrupta, acompanhada por mudanças estruturais e eletrônicas, especificamente a desproporção de ligação (distorção de modo de respiração dos octaedros NiO6) e a desproporção de carga. Um desafio central no campo é o desmembramento das forças motrizes desta transição: se ela é impulsionada primordialmente por efeitos eletrônicos correlacionados (instabilidade de carga) ou pela dinâmica estrutural da rede. A Teoria do Funcional da Densidade tradicional com correção de Hubbard (DFT+U) trata a estrutura e os estados eletrônicos de maneira enviesada e emaranhada, tornando difícil isolar o papel específico da dinâmica da rede ao desencadear a transição de fase.

Metodologia
Para desacoplar as respostas estruturais e eletrônicas, os autores empregaram uma abordagem de potencial interatômico baseado em aprendizado de máquina combinada com simulações de dinâmica molecular (MD) multiescala, utilizando SmNiO3 como um sistema prototípico.

• Geração de Dados: Conjuntos de dados de treinamento foram gerados usando cálculos de DFT+U com polarização de spin (VASP, funcional PBEsol, Ueff = 2,0 eV) em supercélulas representando tanto a fase monoclínica de baixa temperatura ($P2_1/n$) quanto a fase ortorrômbica de alta temperatura ($Pbnm$).
• Treinamento do Potencial: Um modelo de Potencial Profundo (DP) foi treinado utilizando o framework de aprendizado ativo DP-GEN. O modelo foi treinado em 1.429 configurações, alcançando erros de raiz quadrada média (RMSE) de 2,03 meV/átomo para energia e aproximadamente 80–85 meV/Å para forças, validando sua precisão contra o DFT+U.
• Simulação: Simulações de MD de larga escala foram realizadas em uma supercélula de 7.680 átomos utilizando o método de Dinâmica Molecular de Potencial Profundo (DPMD) dentro do ensemble isotérmico-isobárico (NPT). Crucialmente, este framework de MD clássica utiliza a superfície de energia potencial derivada do DFT, mas não inclui explicitamente os graus de liberdade eletrônicos durante a simulação, isolando assim a resposta estrutural às flutuações térmicas.

Resultos Principais

• Evolução Estrutural: As simulações reproduziram com sucesso a desproporção de ligação característica da fase $P2_1/n$ em baixas temperaturas (50 K), mostrando volumes distintos para os octaedros NiO6 expandidos ("L") e comprimidos ("S").
• Transição Impulsionada pela Temperatura: À medida que a temperatura aumentava, as flutuações térmicas alargavam a distribuição de volume dos octaedros. A 320 K, a distinção entre as unidades "L" e "S" era significativamente suprimida. A 340 K, as distribuições de volume fundiram-se completamente, indicando o desaparecimento da desproporção de ligação e uma transição espontânea para a fase centrossimétrica $Pbnm$. O processo inverso (resfriamento) confirmou o surgimento da desproporção a 340 K.
• Temperatura e Pressão Críticas: A temperatura crítica ($T_c$) simulada para a transição estrutural é de 340 K. O estudo demonstrou adicionalmente que a pressão hidrostática modula esta transição; aumentar a pressão de 0 para 10 GPa suprime monotonicamente $T_c$ de 340 K para 225 K, consistente com as tendências experimentais na família RNiO3.
• Comparação com Experimento: Os autores observam que o $T_c$ simulado (340 K) é inferior ao valor experimental (~400 K) e que a transição parece menos abrupta sem a histerese típica de transições de primeira ordem.

Significância e Alegações
O artigo afirma que estas descobertas fornecem insights atomísticos sobre a evolução estrutural que desencadeia a transição de fase em RNiO3. Ao isolar a resposta estrutural dos efeitos eletrônicos, os autores demonstram que a transição de fase estrutural é intrinsecamente impulsionada pela temperatura e ocorre espontaneamente via distorções coletivas da rede.

Os autores concluem que a dinâmica da rede desempenha um papel significativo no direcionamento da transição de fase e da subsequente MIT. Eles postulam que, embora os efeitos eletrônicos sejam essenciais para a descrição completa da MIT, a transição estrutural em si pode ser entendida como um fenômeno impulsionado pela temperatura. O trabalho sugere que uma imagem física completa da MIT em nickelatos requer o comportamento cooperativo de mecanismos de rede e eletrônicos, oferecendo uma perspectiva mais clara do que estudos onde esses efeitos permanecem emaranhados. Os autores não pretendem ter resolvido o mecanismo eletrônico completo, mas sim ter isolado e caracterizado com sucesso o componente estrutural da transição.