Polar vibrational excitations and glass-like response in relaxor PMN

Resumen Técnico: Excitaciones Vibracionales Polares y Respuesta Tipo Vidrio en el Relaxor PMN

Planteamiento del Problema
Los ferroeléctricos relaxores, como el niobato de magnesio y plomo (Pb(Mg$_{1/3}$Nb$_{2/3}$)O$_3$ o PMN), presentan una paradoja física única: conservan un orden cristalino de largo alcance mientras exhiben respuestas dieléctricas y térmicas macroscópicas características de los vidrios estructurales. A diferencia de los ferroeléctricos convencionales (p. ej., PbTiO$_3$ o PTO), donde los modos ópticos polares son discretos y los fonones acústicos no son polares, los relaxores muestran una transición dieléctrica difusa, susceptibilidad dependiente de la frecuencia y un calor específico anómalo a baja temperatura. El desafío central abordado en este trabajo es elucidar los mecanismos microscópicos mediante los cuales el desorden composicional acopla los grados de libertad elásticos, eléctricos-dipolares y composicionales para remodelar el espectro de excitaciones vibracionales, generando así las anomalías tipo vidrio en una red cristalina.

Metodología
Los autores emplean modelos atomísticos precisos basados en interacciones interatómicas aprendidas por aprendizaje automático y polarización, entrenadas con datos de la teoría del funcional de la densidad (DFT) bajo la aproximación SCAN.

• Potenciales Interatómicos: Una red neuronal de Potencial Profundo (DP) modela la superficie de energía potencial, mientras que una red de Wannier Profundo (DW) predice la polarización rastreando el desplazamiento de los centros de Wannier respecto a los núcleos atómicos.
• Modelado Estructural: Para capturar el desorden composicional a escala nanométrica inherente al PMN, las celdas de simulación se construyen como mezclas aleatorias de regiones químicamente ordenadas (COR) y regiones químicamente desordenadas (CDR). Las COR presentan una estructura de doble perovskita con un ordenamiento específico en las subredes de Nb/Mg, mientras que las CDR exhiben una ocupación aleatoria de Mg/Nb.
• Protocolos de Simulación:
  • Dinámica Molecular (MD): Se utilizan simulaciones NVT-MD (50–800 K) para estudiar la susceptibilidad dieléctrica, la ruptura de la ergodicidad y la formación de regiones nanopolares (PNR).
  • Respuesta Lineal: Para analizar los espectros vibracionales, los autores diagonalizan las matrices dinámicas de "estructuras inherentes" (configuraciones minimizadas en energía) derivadas de trayectorias de enfriamiento. Esto permite el cálculo de frecuencias propias, modos propios y fuerzas de oscilador sin las limitaciones de escala temporal de la MD para el análisis espectral.
  • Termodinámica: El calor específico vibracional se calcula integrando la densidad de estados vibracionales (VDOS), contabilizando la anharmonicidad en los modos de baja frecuencia mediante la solución numérica de la ecuación de Schrödinger para potenciales unidimensionales efectivos.

Contribuciones Clave y Resultados

1. Respuesta Dieléctrica y Ruptura de la Ergodicidad:
   Las simulaciones reproducen cualitativamente la dependencia de la temperatura de la susceptibilidad dieléctrica ($\chi'$) y la pérdida progresiva de la ergodicidad al enfriar. El modelo captura la transición difusa paraeléctrico-relaxor y el desplazamiento de la frecuencia del pico de susceptibilidad, consistente con la ley de Vogel-Tammann-Fulcher (VTF) observada en experimentos. Aunque la magnitud absoluta del pico de susceptibilidad es menor que los valores experimentales (atribuido a las limitaciones de tamaño finito y tiempo finito que impiden acceder a escalas de relajación de microsegundos), el modelo distingue con éxito entre la predicción del teorema de fluctuación-disipación y la respuesta de la derivada numérica, señalando el inicio de la dinámica vítrea.

2. Continuo de Modos Vibracionales Polares:
   Un hallazgo principal es que el desorden composicional en el PMN acopla el movimiento elástico a la polarización eléctrica a través de todo el espectro vibracional. A diferencia del PTO, donde solo un número finito de modos ópticos en el centro de zona son activos en el infrarrojo (IR), el PMN exhibe un continuo de excitaciones vibracionales polares. Esto resulta en una respuesta amplia y dependiente de la frecuencia en las regiones del infrarrojo lejano (10–400 cm$^{-1}$) y de microondas. Los espectros de IR lejano calculados (susceptibilidad dieléctrica, conductividad óptica y reflectividad) muestran una excelente concordancia cuantitativa con los datos experimentales a 10 K.

3. Pico de Bosón y Excitaciones Cuasi-Localizadas (QLE):
   El estudio identifica un exceso de modos de baja frecuencia en el PMN, manifestándose como un "pico de bosón" en la VDOS reducida ($D(\nu)/\nu^2$).

   • Localización: A diferencia de los modos fonónicos extendidos en el PTO, los modos de menor frecuencia en el PMN son excitaciones cuasi-localizadas (QLE). Estos modos involucran grandes desplazamientos de un grupo central de átomos (principalmente Pb) con desplazamientos que decaen exponencialmente en la red circundante.
   • Escalamiento y Anarmonicidad: La VDOS de estas QLE sigue una ley de escalamiento universal $D(\nu) \propto \nu^4$ para frecuencias por debajo de 20 cm$^{-1}$, independiente de la dimensionalidad, reflejando el comportamiento encontrado en vidrios estructurales. Además, estos modos exhiben una pronunciada anarmonicidad, con superficies de energía potencial dominadas por términos cúrticos en lugar de los términos cuadráticos típicos de los fonones armónicos.

4. Anomalías Térmicas:
   Los modos anómalos de baja frecuencia explican directamente las propiedades térmicas del PMN. El calor específico vibracional ($C_{vib}$) calculado para $T > 5$ K reproduce el pico del calor específico reducido ($C/T^3$) experimental.

   • En el PMN, el pico ocurre a $\approx 10.8$ K con una magnitud de 1.43 mJ K$^{-4}$ mol$^{-1}$, significativamente mayor y desplazado a temperaturas más bajas en comparación con el PTO (pico en $\approx 16.4$ K).
   • Los autores atribuyen este incremento a la contribución del pico de bosón y las QLE. Si bien las QLE contribuyen significativamente, la curva "PMN sin QLE" indica que el cuerpo más amplio del pico de bosón (modos con localización intermedia) también desempeña un papel crítico.

Significancia y Reivindicaciones
El artículo afirma establecer un vínculo microscópico directo entre el desorden composicional de un relaxor cristalino y sus respuestas macroscópicas tipo vidrio. Al demostrar que el desorden induce un continuo de modos polares y genera QLE con escalamiento cúrtico y fuerte anarmonicidad, el trabajo unifica las anomalías dieléctricas, ópticas y térmicas del PMN bajo un solo marco: el espectro de excitación polar inducido por el desorden.

Los autores posicionan su trabajo como el primer reporte de QLE en un modelo realista derivado de primeros principios de un relaxor. Argumentan que estas QLE no son meros artefactos de vidrios estructurales, sino que son intrínsecas al estado relaxor, surgiendo del acoplamiento específico de los grados de libertad elásticos y dipolares en una perovskita químicamente desordenada. Si bien el estudio no resuelve el comportamiento divergente del calor específico por debajo de 1 K (atribuido a modos de tunelamiento cuántico más allá del análisis vibracional armónico/anharmónico), logra cerrar la brecha entre la naturaleza cristalina del PMN y sus firmas termodinámicas y dieléctricas de tipo vidrio. Los hallazgos sugieren que el comportamiento "vítreo" en los relaxores es una consecuencia dinámica del paisaje vibracional específico creado por el desorden composicional, más que una amorfización estructural estática.