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

Adoucissement Raman induit par les lacunes d'oxygène dans le catalyseur Fe2_2(MoO4_4)3_3

Informations de base

  • ID de l'article: 2510.12746
  • Titre: Oxygen-vacancy-induced Raman softening in the catalyst Fe2_2(MoO4_4)3_3
  • Auteurs: Young-Joon Song, Roser Valentí (Université Goethe de Francfort)
  • Classification: cond-mat.mtrl-sci (Physique de la matière condensée - Science des matériaux)
  • Date de publication: 15 octobre 2025
  • Lien de l'article: https://arxiv.org/abs/2510.12746

Résumé

Le ferromolybdate (Fe2_2(MoO4_4)3_3) est largement utilisé dans les procédés catalytiques industriels pour l'oxydéshydrogénation du méthanol en formaldéhyde, mais son mécanisme catalytique détaillé reste peu clair. Les mesures récentes de spectroscopie Raman in situ et d'impédance suggèrent que la région de volume de Fe2_2(MoO4_4)3_3 pourrait servir de réservoir d'atomes d'oxygène, lesquels peuvent migrer vers la surface pour participer aux réactions catalytiques. Cette conclusion est basée sur la diminution significative de l'intensité Raman observée au cours du processus catalytique, suggérant la formation de défauts atomiques. Cependant, l'origine microscopique de cette diminution d'intensité reste à élucider. Cette étude révèle l'origine des variations d'intensité Raman observées expérimentalement par des calculs de théorie de la fonctionnelle de la densité (DFT). L'analyse des phonons montre que les modes de vibration dominés par l'oxygène (avec une contribution mineure de Mo) apparaissent près de 782 cm−1^{-1}, précisément dans la région de fréquence où la diminution d'intensité Raman a été mesurée expérimentalement. L'utilisation de la méthode du phonon gelé effectif confirme en outre que les vibrations d'oxygène sont la cause principale de la diminution calculée de l'intensité Raman.

Contexte et motivation de la recherche

Contexte du problème

  1. Importance de l'application catalytique: Fe2_2(MoO4_4)3_3 est un catalyseur industriel important pour l'oxydéshydrogénation du méthanol en formaldéhyde, utilisé depuis des décennies
  2. Mécanisme peu clair: Malgré son utilisation généralisée, le mécanisme détaillé du processus catalytique reste mal compris, en particulier le rôle des atomes d'oxygène
  3. Observations expérimentales: Les expériences Raman in situ de Schumacher et al. ont révélé une diminution significative de l'intensité Raman près de 782 cm−1^{-1} au cours du processus catalytique à 500°C

Motivation de la recherche

  1. Besoin d'explication théorique: La diminution d'intensité Raman observée expérimentalement manque d'explication théorique microscopique
  2. Élucidation du mécanisme catalytique: Nécessité de comprendre au niveau atomique le rôle des lacunes d'oxygène dans le processus catalytique
  3. Relation structure-propriété: Établir une relation quantitative entre les défauts d'oxygène et les variations du spectre Raman

Contributions principales

  1. Révélation de l'origine microscopique de la diminution d'intensité Raman: Les calculs DFT démontrent que les lacunes d'oxygène sont la cause principale de la diminution d'intensité Raman observée expérimentalement
  2. Établissement de la relation structure-spectre: Identification du pic Raman près de 782 cm−1^{-1} correspondant au mode de vibration d'élongation asymétrique de MoO4_4 dominé par l'oxygène
  3. Proposition de la méthode du phonon gelé effectif: Développement d'une méthode théorique pour geler sélectivement les modes de phonons d'atomes spécifiques afin d'analyser les variations d'intensité Raman
  4. Explication des phénomènes expérimentaux: Les résultats des calculs théoriques sont en excellent accord avec les observations expérimentales, expliquant pourquoi aucun décalage de pic ou élargissement n'a été observé au cours du processus catalytique

Détails des méthodes

Méthodes de calcul

Calculs de théorie de la fonctionnelle de la densité (DFT):

  • Utilisation du Vienna Ab initio Simulation Package (VASP)
  • Méthode des ondes planes augmentées projetées (PAW)
  • Fonctionnelle d'échange-corrélation: PBESol (approximation du gradient généralisé modifiée)
  • Correction Hubbard U: valeur U de 4 eV appliquée aux orbitales 3d du Fe
  • Énergie de coupure des ondes planes: 500 eV
  • Grille de points k: 2×3×2 pour la phase monoclinique, 3×4×4 pour la phase orthorhombique

Calculs de phonons

  1. Méthode des déplacements finis: Calcul de la matrice des constantes de force avec un déplacement atomique de 0,01 Å
  2. Matrice dynamique: Construction via le package phonopy et obtention de la dispersion des phonons et de la densité d'états
  3. Calcul de l'intensité Raman: Calcul de l'intensité Raman par combinaison des résultats du tenseur diélectrique et des vecteurs propres des phonons

Méthode technique innovante

Méthode du phonon gelé effectif:

  • Gel sélectif des modes de vibration de phonons d'atomes spécifiques
  • Préservation de la symétrie cristalline
  • Détermination de la contribution de chaque atome par comparaison des variations d'intensité Raman avant et après le gel
  • Cette méthode ne considère pas les effets de charge, mais permet une séparation efficace des contributions de phonons

Configuration expérimentale

Modèles de structure

  1. Phase monoclinique basse température: Groupe d'espace P2₁/c, Z=8
    • Paramètres de maille: a=15,869 Å, b=9,373 Å, c=15,868 Å, β=108,36°
  2. Phase orthorhombique haute température: Groupe d'espace Pbcn, Z=4
    • Paramètres de maille: a=12,867 Å, b=9,287 Å, c=9,373 Å

Paramètres de calcul

  • Polarisation de spin: Arrangement ferromagnétique simple pour obtenir une bande interdite semi-conductrice (~2,7 eV)
  • Optimisation de structure: Critère de convergence des forces < 1 meV/Å
  • Élargissement gaussien: 0,05 eV
  • Élargissement du spectre Raman: Élargissement lorentzien, largeur à mi-hauteur de 10 cm−1^{-1}

Analyse comparative

Comparaison des différences structurales entre les phases monoclinique et orthorhombique par analyse des modes de symétrie, utilisant le programme AMPLIMODES pour déterminer la matrice de transformation de la transition de phase.

Résultats expérimentaux

Analyse de la structure des phonons

  1. Séparation en fréquence: Le spectre des phonons se divise en deux régions
    • Région basse fréquence (0-403 cm−1^{-1}): Contribution des atomes Fe, Mo et O
    • Région haute fréquence (780-1001 cm−1^{-1}): Disparition de la contribution du Fe, vibrations principalement Mo-O
  2. Modes de vibration clés:
    • 789 cm−1^{-1}: Mode d'élongation asymétrique de MoO4_4, intensité Raman maximale
    • 972 cm−1^{-1}: Mode d'élongation symétrique de MoO4_4, l'atome Mo est pratiquement immobile

Résultats du calcul de l'intensité Raman

Fréquence (cm−1^{-1})Représentation irréductibleIntensité normalisée (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

Résultats de l'analyse du phonon gelé

Par gel sélectif des modes de phonons de différents atomes:

  • Gel de l'atome O4: Diminution de l'intensité Raman de 8,5% à 789,68 cm−1^{-1}
  • Gel de l'atome O5: Diminution de l'intensité Raman de 7,6% à 786,18 cm−1^{-1}
  • Gel de l'atome Mo: Diminution d'intensité environ moitié moins importante que celle des atomes d'oxygène

Relaxation de structure avec lacune d'oxygène

La relaxation de structure contenant une lacune d'oxygène montre:

  • Les unités MoO3_3 triangulaires tendent à retrouver la géométrie tétraédrique MoO4_4
  • Les unités FeO5_5 pentacoordinées restent stables
  • Changements significatifs des paramètres de maille, groupe d'espace réduit à P1

Travaux connexes

Propriétés multifonctionnelles de Fe2_2(MoO4_4)3_3

  1. Propriétés magnétiques: Transition ferrimagnétique de type L en dessous de 12 K, phase multiferroïque induite par champ magnétique externe
  2. Applications électrochimiques: Matériau d'électrode pour batteries Li/Na-ion, avec excellente performance cyclique
  3. Applications catalytiques: Catalyseur industriel pour la synthèse du formaldéhyde à partir du méthanol, la concentration superficielle de Mo affecte la sélectivité

Études spectroscopiques in situ

Le travail fondateur de Schumacher et al. a proposé pour la première fois le concept de réservoir d'oxygène en volume, fournissant la base expérimentale pour cette étude.

Méthodes de calcul théorique

Applications pertinentes de la méthode DFT+U dans le traitement de la structure électronique et des propriétés de phonons des oxydes de métaux de transition.

Conclusions et discussion

Conclusions principales

  1. Les lacunes d'oxygène sont essentielles: Les lacunes d'oxygène sont la cause principale de la diminution d'intensité Raman observée près de 782 cm−1^{-1}
  2. Identification des modes de vibration: La région de fréquence correspond aux vibrations d'élongation asymétrique de MoO4_4 dominées par l'oxygène
  3. Préservation de la symétrie: Le processus de diffusion d'oxygène ne modifie pas significativement la symétrie locale, expliquant l'absence de décalage de pic ou d'élargissement observée expérimentalement
  4. Mécanisme catalytique: Soutient le mécanisme catalytique du volume comme réservoir d'oxygène

Limitations

  1. Modèle simplifié: La méthode du phonon gelé ne considère pas les effets de charge et les interactions à longue portée
  2. Lacune unique: Seule une lacune d'oxygène unique est considérée, le processus catalytique réel pourrait impliquer plusieurs lacunes ou des amas de lacunes
  3. Effets de température: Les calculs sont basés sur 0 K, sans considération des effets de vibration thermique à température finie
  4. Effets de surface: Principalement axé sur les propriétés en volume, description limitée des sites actifs catalytiques en surface

Directions futures

  1. Étude de lacunes multiples: Étudier l'effet de la concentration et de la distribution des lacunes sur le spectre Raman
  2. Processus dynamiques: Calculer les chemins de diffusion des atomes d'oxygène et les barrières énergétiques
  3. Modélisation de surface: Établir des modèles de surface pour étudier les sites actifs catalytiques
  4. Vérification expérimentale: Recommandation de vérifier les changements d'état d'oxydation du Fe par spectroscopie Mössbauer et autres techniques

Évaluation approfondie

Points forts

  1. Combinaison théorie-expérience: Explique bien les variations du spectre Raman observées expérimentalement
  2. Innovation méthodologique: La méthode du phonon gelé fournit un nouvel outil pour analyser les variations d'intensité Raman
  3. Calculs rigoureux: La méthode DFT+U traite raisonnablement le système électronique fortement corrélé
  4. Image physique claire: Élucide le mécanisme d'action des lacunes d'oxygène au niveau atomique

Insuffisances

  1. Simplification du modèle: La méthode du phonon gelé ignore les changements de structure électronique
  2. Concentration de lacunes: Pas d'étude systématique de l'effet de différentes concentrations de lacunes
  3. Informations cinétiques manquantes: Manque d'analyse cinétique de la migration des atomes d'oxygène
  4. Comparaison expérimentale limitée: Principalement basée sur la vérification d'un ensemble de données expérimentales

Impact

  1. Théorie catalytique: Fournit un soutien théorique important pour comprendre le mécanisme catalytique de l'oxydéshydrogénation
  2. Analyse spectroscopique: Fournit une base théorique pour l'application de la spectroscopie Raman dans la recherche catalytique
  3. Conception de matériaux: Fournit des idées de conception pour développer de nouveaux catalyseurs à base d'oxydes
  4. Contribution méthodologique: La méthode du phonon gelé peut être généralisée à d'autres systèmes de matériaux

Scénarios d'application

  1. Catalyseurs à base d'oxydes: Applicable à l'étude des effets des lacunes d'oxygène dans divers catalyseurs à base d'oxydes
  2. Interprétation de spectroscopie in situ: Fournit des méthodes pour l'interprétation théorique des données de spectroscopie Raman in situ
  3. Ingénierie des défauts: Guide l'amélioration des performances catalytiques par contrôle des lacunes d'oxygène
  4. Catalyse hétérogène: Fournit un cadre théorique pour comprendre les processus de transport d'oxygène en catalyse hétérogène

Références

Cet article cite 34 références importantes couvrant les propriétés structurales, magnétiques, électrochimiques et catalytiques de Fe2_2(MoO4_4)3_3, ainsi que les méthodes de calcul théorique pertinentes et les techniques expérimentales, fournissant une base théorique solide pour la recherche.

Évaluation globale: Cet article est un travail de calcul théorique de haute qualité qui explique avec succès l'anomalie du spectre Raman observée expérimentalement et fournit un soutien théorique important pour comprendre le mécanisme de fonctionnement du catalyseur ferromolybdate. La méthode est innovante, les calculs rigoureux, les conclusions fiables, avec une valeur académique et une signification pratique importantes.