2025-11-17T12:46:19.904493

Oxygen-vacancy-induced Raman softening in the catalyst Fe$_2$(MoO$_4$)$_3$

Song, Valentí
Although iron molybdate (Fe$_2$(MoO$_4$)$_3$) has been commercially utilized for the production of formaldehyde from methanol via oxidative dehydrogenation, the detailed mechanism during the catalytic process remains unclear. Recent operando Raman and impedance measurements of the reaction suggested that the bulk region of Fe$_2$(MoO$_4$)$_3$ acts as a reservoir of oxygen atoms that can migrate to the surface to participate in catalysis. This conclusion was drawn from the observed significant reduction in Raman intensity during the catalytic process which implies the formation of atomic defects. However, the microscopic origin of this reduction remains to be clarified. In this work, we performed density functional theory (DFT) calculations to elucidate the origin of the experimentally observed Raman intensity variation. Our phonon analysis reveals that oxygen-dominated vibrational modes, with a small Mo contribution, occur near 782 cm$^{-1}$-the same frequency region where the Raman intensity reduction was measured. Using an effective frozen-phonon approach, we further demonstrate that oxygen vibrations are primarily responsible for the decrease in calculated Raman intensity. Moreover, structural relaxation of Fe$_2$(MoO$_4$)$_3$ including an oxygen vacancy suggests that oxygen diffusion from the bulk to the surface should occur without significant alteration of the local symmetry, consistent with the absence of measurable peak shifts or broadening in the experimental Raman spectra.
academic

Ablandamiento Raman inducido por vacancias de oxígeno en el catalizador Fe2_2(MoO4_4)3_3

Información Básica

  • ID del artículo: 2510.12746
  • Título: Oxygen-vacancy-induced Raman softening in the catalyst Fe2_2(MoO4_4)3_3
  • Autores: Young-Joon Song, Roser Valentí (Universidad Goethe de Frankfurt)
  • Clasificación: cond-mat.mtrl-sci (Física de la materia condensada - Ciencia de materiales)
  • Fecha de publicación: 15 de octubre de 2025
  • Enlace del artículo: https://arxiv.org/abs/2510.12746

Resumen

El ferromolibdato (Fe2_2(MoO4_4)3_3) se ha utilizado ampliamente en procesos catalíticos industriales para la oxidación deshidrogenante del metanol a formaldehído, pero su mecanismo catalítico detallado sigue siendo poco claro. Las mediciones recientes de espectroscopia Raman in situ e impedancia sugieren que la región de volumen de Fe2_2(MoO4_4)3_3 puede actuar como un depósito de átomos de oxígeno, que pueden migrar a la superficie para participar en reacciones catalíticas. Esta conclusión se basa en la reducción significativa de la intensidad Raman observada durante el proceso catalítico, lo que sugiere la formación de defectos atómicos. Sin embargo, el origen microscópico de esta reducción de intensidad aún requiere aclaración. Este estudio revela el origen de los cambios observados experimentalmente en la intensidad Raman mediante cálculos de teoría del funcional de la densidad (DFT). El análisis de fonones muestra que los modos de vibración dominados por oxígeno (con una pequeña contribución de Mo) aparecen cerca de 782 cm−1^{-1}, precisamente en la región de frecuencia donde se observó experimentalmente la reducción de intensidad Raman. El método de fonón congelado efectivo confirma además que la vibración del oxígeno es la causa principal de la reducción calculada de la intensidad Raman.

Antecedentes y Motivación de la Investigación

Contexto del Problema

  1. Importancia de la aplicación catalítica: Fe2_2(MoO4_4)3_3 es un catalizador industrial importante para la oxidación deshidrogenante del metanol a formaldehído, utilizado durante décadas
  2. Mecanismo poco claro: A pesar de su amplia aplicación, el mecanismo detallado del proceso catalítico sigue siendo poco claro, especialmente el mecanismo de acción de los átomos de oxígeno
  3. Observación experimental: Los experimentos de espectroscopia Raman in situ de Schumacher et al. encontraron una reducción significativa de la intensidad Raman cerca de 782 cm−1^{-1} durante el proceso catalítico a 500°C

Motivación de la Investigación

  1. Necesidad de explicación teórica: La reducción de intensidad Raman observada experimentalmente carece de una explicación teórica microscópica
  2. Aclaración del mecanismo catalítico: Se requiere comprender a nivel atómico el papel de las vacancias de oxígeno en el proceso catalítico
  3. Relación estructura-propiedad: Establecer una relación cuantitativa entre defectos de oxígeno y cambios en el espectro Raman

Contribuciones Principales

  1. Revelación del origen microscópico de la reducción de intensidad Raman: Los cálculos DFT demuestran que las vacancias de oxígeno son la causa principal de la reducción de intensidad Raman observada experimentalmente
  2. Establecimiento de la relación estructura-espectro: Se identifica que el pico Raman cerca de 782 cm−1^{-1} corresponde a un modo de vibración de estiramiento asimétrico de MoO4_4 dominado por oxígeno
  3. Propuesta del método de fonón congelado efectivo: Se desarrolla un método teórico para congelar selectivamente modos de fonones de átomos específicos para analizar cambios de intensidad Raman
  4. Explicación de fenómenos experimentales: Los resultados de cálculos teóricos concuerdan altamente con observaciones experimentales, explicando por qué no se observaron desplazamientos de picos o ensanchamiento durante el proceso catalítico

Detalles de la Metodología

Métodos Computacionales

Cálculos de teoría del funcional de la densidad (DFT):

  • Utilización del paquete Vienna Ab initio Simulation Package (VASP)
  • Método de ondas aumentadas proyectadas (PAW)
  • Funcional de intercambio-correlación: PBESol (aproximación de gradiente generalizado modificada)
  • Corrección Hubbard U: valor U de 4 eV aplicado a orbitales 3d de Fe
  • Energía de corte de ondas planas: 500 eV
  • Malla de puntos k: 2×3×2 para fase monoclínica, 3×4×4 para fase ortorrómbica

Cálculos de Fonones

  1. Método de desplazamiento finito: Cálculo de la matriz de constantes de fuerza utilizando desplazamientos atómicos de 0.01 Å
  2. Matriz dinámica: Construcción mediante el paquete de software phonopy para obtener dispersión de fonones y densidad de estados
  3. Cálculo de intensidad Raman: Cálculo de intensidad Raman mediante resultados de tensor dieléctrico combinados con vectores propios de fonones

Método Técnico Innovador

Método de fonón congelado efectivo:

  • Congelación selectiva de modos de vibración de fonones de átomos específicos
  • Mantenimiento de la simetría cristalina sin cambios
  • Determinación de contribuciones de cada átomo mediante comparación de cambios de intensidad Raman antes y después de la congelación
  • Este método no considera efectos de carga, pero separa efectivamente las contribuciones de fonones

Configuración Experimental

Modelos Estructurales

  1. Fase monoclínica de baja temperatura: Grupo espacial P2₁/c, Z=8
    • Parámetros de red: a=15.869 Å, b=9.373 Å, c=15.868 Å, β=108.36°
  2. Fase ortorrómbica de alta temperatura: Grupo espacial Pbcn, Z=4
    • Parámetros de red: a=12.867 Å, b=9.287 Å, c=9.373 Å

Parámetros de Cálculo

  • Polarización de espín: Disposición ferromagnética simple para obtener brecha de semiconductor (~2.7 eV)
  • Optimización estructural: Criterio de convergencia de fuerzas < 1 meV/Å
  • Ensanchamiento gaussiano: 0.05 eV
  • Ensanchamiento de espectro Raman: Ensanchamiento lorentziano, ancho medio de 10 cm−1^{-1}

Análisis Comparativo

Comparación de diferencias estructurales entre fase monoclínica y ortorrómbica mediante análisis de modos de simetría, utilizando el programa AMPLIMODES para determinar la matriz de transformación de la transición de fase.

Resultados Experimentales

Análisis de Estructura de Fonones

  1. Separación de frecuencias: El espectro de fonones se divide en dos regiones
    • Región de baja frecuencia (0-403 cm−1^{-1}): Contribuciones de átomos Fe, Mo y O
    • Región de alta frecuencia (780-1001 cm−1^{-1}): Desaparición de contribución de Fe, principalmente vibración Mo-O
  2. Modos de vibración clave:
    • 789 cm−1^{-1}: Modo de estiramiento asimétrico de MoO4_4, intensidad Raman máxima
    • 972 cm−1^{-1}: Modo de estiramiento simétrico de MoO4_4, átomo Mo prácticamente inmóvil

Resultados de Cálculo de Intensidad Raman

Frecuencia (cm−1^{-1})Representación irreducibleIntensidad normalizada (I/Imax_{max})
789.68Ag_g1.000
790.28B1g_{1g}0.628
786.18B3g_{3g}0.964
786.87Ag_g0.515
972.05Ag_g0.852

Resultados del Análisis de Fonón Congelado

Mediante congelación selectiva de modos de fonones de diferentes átomos:

  • Congelación de átomo O4: Reducción de intensidad Raman de 8.5% en 789.68 cm−1^{-1}
  • Congelación de átomo O5: Reducción de intensidad Raman de 7.6% en 786.18 cm−1^{-1}
  • Congelación de átomo Mo: Magnitud de reducción de intensidad aproximadamente la mitad de la de átomos de oxígeno

Relajación Estructural con Vacancia de Oxígeno

La relajación estructural con vacancia de oxígeno muestra:

  • Unidades triangulares MoO3_3 tienden a recuperar geometría tetraédrica MoO4_4
  • Unidades pentacoordinadas FeO5_5 mantienen estabilidad
  • Cambios significativos en parámetros de red, grupo espacial reduce a P1

Trabajos Relacionados

Propiedades Multifuncionales de Fe2_2(MoO4_4)3_3

  1. Propiedades magnéticas: Transición ferrimagnética tipo L por debajo de 12 K, fase multiferroica inducida por campo magnético externo
  2. Aplicaciones electroquímicas: Material de electrodo para baterías de iones Li/Na, con buen rendimiento de ciclado
  3. Aplicaciones catalíticas: Catalizador industrial para síntesis de formaldehído a partir de metanol, la concentración de Mo en superficie afecta la selectividad

Estudios de Espectroscopia In Situ

El trabajo pionero de Schumacher et al. propone por primera vez el concepto de depósito de oxígeno de volumen, proporcionando la base experimental para esta investigación.

Métodos de Cálculo Teórico

Aplicaciones del método DFT+U en el tratamiento de estructura electrónica y propiedades de fonones de óxidos de metales de transición.

Conclusiones y Discusión

Conclusiones Principales

  1. Las vacancias de oxígeno son clave: Las vacancias de oxígeno son la causa principal de la reducción de intensidad Raman observada cerca de 782 cm−1^{-1}
  2. Identificación de modo de vibración: La región de frecuencia corresponde a vibración de estiramiento asimétrico de MoO4_4 dominada por oxígeno
  3. Preservación de simetría: El proceso de difusión de oxígeno no cambia significativamente la simetría local, explicando por qué no se observaron desplazamientos de picos o ensanchamiento en experimentos
  4. Mecanismo catalítico: Apoya el mecanismo catalítico del volumen como depósito de oxígeno

Limitaciones

  1. Modelo simplificado: El método de fonón congelado no considera efectos de carga e interacciones de largo alcance
  2. Vacancia única: Solo se considera una vacancia de oxígeno individual, el proceso catalítico real puede implicar múltiples vacancias o aglomerados de vacancias
  3. Efectos de temperatura: Los cálculos se basan en 0 K, sin considerar efectos de vibración térmica a temperatura finita
  4. Efectos de superficie: Se enfoca principalmente en propiedades de volumen, con descripción limitada de sitios activos catalíticos en la superficie

Direcciones Futuras

  1. Investigación de múltiples vacancias: Estudiar el efecto de concentración y distribución de vacancias en el espectro Raman
  2. Procesos dinámicos: Calcular rutas de difusión de átomos de oxígeno y barreras energéticas
  3. Modelado de superficie: Establecer modelos de superficie para investigar sitios activos catalíticos
  4. Verificación experimental: Se sugiere verificar cambios en el estado de oxidación de Fe mediante técnicas como espectroscopia Mössbauer

Evaluación Profunda

Fortalezas

  1. Combinación de teoría y experimento: Explica bien los cambios en el espectro Raman observados experimentalmente
  2. Innovación metodológica: El método de fonón congelado proporciona una nueva herramienta para analizar cambios de intensidad Raman
  3. Cálculos rigurosos: El método DFT+U maneja apropiadamente sistemas electrónicos fuertemente correlacionados
  4. Imagen física clara: Aclara el mecanismo de acción de vacancias de oxígeno a nivel atómico

Insuficiencias

  1. Simplificación del modelo: El método de fonón congelado ignora los efectos de cambios en la estructura electrónica
  2. Concentración de vacancias: No se estudia sistemáticamente el efecto de diferentes concentraciones de vacancias
  3. Información de dinámica faltante: Falta análisis de dinámica de migración de átomos de oxígeno
  4. Comparación experimental limitada: La verificación se basa principalmente en un conjunto de datos experimentales

Impacto

  1. Teoría catalítica: Proporciona apoyo teórico importante para comprender mecanismos de catálisis de oxidación deshidrogenante
  2. Análisis espectral: Proporciona base teórica para la aplicación de espectroscopia Raman en investigación catalítica
  3. Diseño de materiales: Proporciona ideas de diseño para el desarrollo de nuevos catalizadores de óxidos
  4. Contribución metodológica: El método de fonón congelado puede generalizarse a otros sistemas de materiales

Escenarios de Aplicación

  1. Catalizadores de óxidos: Aplicable al estudio de efectos de vacancias de oxígeno en varios catalizadores de óxidos
  2. Interpretación de espectroscopia in situ: Proporciona método para interpretación teórica de datos de espectroscopia Raman in situ
  3. Ingeniería de defectos: Guía la mejora del rendimiento catalítico mediante regulación de vacancias de oxígeno
  4. Catálisis heterogénea: Proporciona marco teórico para comprender procesos de transporte de oxígeno en catálisis heterogénea

Referencias Bibliográficas

Este artículo cita 34 referencias importantes que abarcan propiedades estructurales, magnéticas, electroquímicas y catalíticas de Fe2_2(MoO4_4)3_3, así como métodos de cálculo teórico relevantes y técnicas experimentales, proporcionando una base teórica sólida para la investigación.

Evaluación General: Este es un artículo de cálculo teórico de alta calidad que explica exitosamente el fenómeno anómalo de espectroscopia Raman observado experimentalmente, proporcionando apoyo teórico importante para comprender el mecanismo de funcionamiento del catalizador de ferromolibdato. El método es altamente innovador, los cálculos son rigurosos, las conclusiones son confiables y posee importante valor académico y significado práctico.