Polar vibrational excitations and glass-like response in relaxor PMN

Resumo Técnico: Excitações Vibracionais Polares e Resposta Semelhante ao Vidro em PMN Relaxor

Definição do Problema
Ferroelétricos relaxores, como o niobato de magnésio e chumbo (Pb(Mg$_{1/3}$Nb$_{2/3}$)O$_3$ ou PMN), apresentam um paradoxo físico único: eles mantam a ordem cristalina de longo alcance enquanto exibem respostas dielétricas e térmicas macroscópicas características de vidros estruturais. Diferente de ferroelétricos convencionais (ex: PbTiO$_3$ ou PTO), onde os modos ópticos polares são discretos e os fônons acústicos não são polares, os relaxores exibem uma transição dielétrica difusa, suscetibilidade dependente da frequência e calor específico anômalo em baixas temperaturas. O desafio central abordado neste trabalho é elucidar os mecanismos microscópicos pelos quais o desordem composicional acopla os graus de liberdade elástico, dipolo-elétrico e composicional para remodelar o espectro de excitação vibracional, gerando assim as anomalias semelhantes ao vidro em uma rede cristalina.

Metodologia
Os autores utilizam modelos atomísticos precisos baseados em interações interatômicas aprendidas por máquina e polarização, treinadas em dados de teoria do funcional da densidade (DFT) de primeira ordem utilizando a aproximação SCAN.

• Potenciais Interatômicos: Uma rede neural de Potencial Profundo (DP) modela a superfície de energia potencial, enquanto uma rede de Wannier Profundo (DW) prevê a polarização rastreando o deslocamento dos centros de Wannier em relação aos núcleos atômicos.
• Modelagem Estrutural: Para capturar a desordem composicional em escala nanométrica inerente ao PMN, as células de simulação são construídas como misturas aleatórias de regiões quimicamente ordenadas (COR) e regiões quimicamente desordenadas (CDR). As CORs apresentam uma estrutura de dupla perovskita com um ordenamento específico de sublatice Nb/Mg, enquanto as CDRs exibem ocupação aleatória de Mg/Nb.
• Protocolos de Simulação:
  • Dinâmica Molecular (MD): Simulações NVT-MD (50–800 K) são usadas para estudar a suscetibilidade dielétrica, a quebra de ergodicidade e a formação de nanorregiões polares (PNR).
  • Resposta Linear: Para analisar os espectros vibracionais, os autores diagonalizam matrizes dinâmicas de "estruturas inerentes" (configurações de energia minimizada) derivadas de trajetórias de resfriamento. Isso permite o cálculo de frequências próprias, modos próprios e forças de oscilador sem as limitações de escala de tempo da MD para análise espectral.
  • Termodinâmica: O calor específico vibracional é calculado integrando a densidade de estados vibracionais (VDOS), contabilizando a anharmonicidade em modos de baixa frequência via solução numérica da equação de Schrödinger para potenciais unidimensionais efetivos.

Principais Contribuições e Resultados

1. Resposta Dielétrica e Quebra de Ergodicidade:
   As simulações reproduzem qualitativamente a dependência de temperatura da suscetibilidade dielétrica ($\chi'$) e a progressiva perda de ergodicidade ao resfriar. O modelo captura a transição para-elétrica-relaxora difusa e o deslocamento da frequência da susceptibilidade, consistente com a lei de Vogel-Tammann-Fulcher (VTF) observada experimentalmente. Embora a magnitude absoluta do pico de suscetibilidade seja menor que os valores experimentais (atribuída a limitações de tamanho e tempo finito que impedem o acesso às escalas de relaxação de microssegundos), o modelo distingue com sucesso entre a previsão do teorema de flutuação-dissipação e a resposta da derivada numérica, sinalizando o início da dinâmica vítrea.

2. Contínuo de Modos Vibracionais Polares:
   Uma descoberta primária é que a desordem composicional no PMN acopla o movimento elástico à polarização elétrica em todo o espectro vibracional. Diferente do PTO, onde apenas um número finito de modos ópticos de centro de zona são ativos no infravermelho (IR), o PMN exibe um contínuo de excitações vibracionais polares. Isso resulta em uma resposta ampla e dependente da frequência nas regiões do infravermelho distante (10–400 cm$^{-1}$) e de micro-ondas. Os espectros de infravermelho distante calculados (suscetibilidade dielétrica, condutividade óptica e refletividade) mostram excelente concordância quantitativa com dados experimentais a 10 K.

3. Pico de Boson e Excitações Quase-Localizadas (QLE):
   O estudo identifica um excesso de modos de baixa frequência no PMN, manifestando-se como um "pico de boson" na VDOS reduzida ($D(\nu)/\nu^2$).

   • Localização: Ao contrário dos modos de fônons estendidos no PTO, os modos de frequência mais baixa no PMN são excitações quase-localizadas (QLE). Esses modos envolvem grandes deslocamentos de um grupo central de átomos (principalmente Pb) com deslocamentos que decaem exponencialmente na rede circundante.
   • Escalonamento e Anarmonicidade: A VDOS dessas QLE segue uma lei de escalonamento universal $D(\nu) \propto \nu^4$ para frequências abaixo de 20 cm$^{-1}$, independente da dimensionalidade, espelhando o comportamento encontrado em vidros estruturais. Além disso, esses modos exibem anarmonicidade pronunciada, com superfícies de energia potencial dominadas por termos quárticos em vez dos termos quadráticos típicos de fônons harmônicos.

4. Anomalias Térmicas:
   Os modos anômalos de baixa frequência explicam diretamente as propriedades térmicas do PMN. O calor específico vibracional ($C_{vib}$) calculado para $T > 5$ K reproduz o pico de calor específico reduzido ($C/T^3$) experimental.

   • No PMN, o pico ocorre em $\approx 10,8$ K com uma magnitude de 1,43 mJ K$^{-4}$ mol$^{-1}$, significativamente maior e deslocado para temperaturas mais baixas comparado ao PTO (pico em $\approx 16,4$ K).
   • Os autores atribuem esse aumento à contribuição do pico de boson e das QLE. Embora as QLE contribuam significativamente, a curva "PMN sem QLE" indica que o corpo mais amplo do pico de boson (modos com localização intermediária) também desempenha um papel crítico.

Significância e Alegações
O artigo afirma estabelecer um elo microscópico direto entre a desordem composicional de um relaxor cristalino e suas respostas macroscópicas semelhantes ao vidro. Ao demonstrar que a desordem induz um contínuo de modos polares e gera QLE com escalonamento quártico e forte anarmonicidade, o trabalho unifica as anomalias dielétricas, ópticas e térmicas do PMN sob um único quadro: o espectro de excitação polar induzido pela desordem.

Os autores posicionam seu trabalho como o primeiro relato de QLE em um modelo realista derivado de primeiros princípios de um relaxor. Eles argumentam que essas QLE não são meros artefatos de vidro estrutural, mas intrínsecas ao estado relaxor, surgindo do acoplamento específico de graus de liberdade elásticos e dipolares em uma perovskita quimicamente desordenada. Embora o estudo não resolva o comportamento divergente do calor específico abaixo de 1 K (atribuído a modos de tunelamento quântico além da análise vibracional harmônica/anharmônica), ele consegue fazer a ponte entre a natureza cristalina do PMN e suas assinaturas termodinâmicas e dielétricas semelhantes ao vidro. Os achados sugerem que o comportamento "vítreo" em relaxores é uma consequência dinâmica do cenário vibracional específico criado pela desordem composicional, e não uma amorfização estrutural estática.