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

Sauerstoffleerstellen-induzierte Raman-Erweichung im Katalysator Fe2_2(MoO4_4)3_3

Grundlegende Informationen

  • Papier-ID: 2510.12746
  • Titel: Oxygen-vacancy-induced Raman softening in the catalyst Fe2_2(MoO4_4)3_3
  • Autoren: Young-Joon Song, Roser Valentí (Goethe-Universität Frankfurt)
  • Klassifizierung: cond-mat.mtrl-sci (Kondensierte Materie – Materialwissenschaften)
  • Veröffentlichungsdatum: 15. Oktober 2025
  • Papierlink: https://arxiv.org/abs/2510.12746

Zusammenfassung

Eisenmolybdate (Fe2_2(MoO4_4)3_3) werden häufig in industriellen katalytischen Prozessen zur oxidativen Dehydrierung von Methanol zu Formaldehyd verwendet, doch der detaillierte Katalysemechanismus bleibt unklar. Kürzliche operando-Raman- und Impedanzmessungen deuten darauf hin, dass die Volumenphase von Fe2_2(MoO4_4)3_3 als Sauerstoffspeicher fungieren könnte, wobei Sauerstoffatome zur Oberfläche migrieren können, um an katalytischen Reaktionen teilzunehmen. Diese Schlussfolgerung basiert auf einer signifikanten Abnahme der Raman-Intensität während des katalytischen Prozesses, was auf die Bildung von Atomdefekten hindeutet. Der mikroskopische Ursprung dieser Intensitätsabnahme bleibt jedoch noch zu klären. Diese Studie nutzt Dichtefunktionaltheorie-Berechnungen (DFT), um den Ursprung der experimentell beobachteten Raman-Intensitätsänderungen zu enthüllen. Phononenanalysen zeigen, dass sauerstoffdominierte Schwingungsmoden (mit geringem Mo-Beitrag) in der Nähe von 782 cm−1^{-1} auftreten, genau in dem Frequenzbereich, in dem experimentell eine Abnahme der Raman-Intensität gemessen wurde. Weitere Analysen mit der effektiven gefrorenen Phonon-Methode bestätigen, dass Sauerstoffschwingungen die Hauptursache für die berechnete Raman-Intensitätsabnahme sind.

Forschungshintergrund und Motivation

Problemhintergrund

  1. Katalytische Bedeutung: Fe2_2(MoO4_4)3_3 ist ein wichtiger industrieller Katalysator für die oxidative Dehydrierung von Methanol zu Formaldehyd und wird seit Jahrzehnten verwendet
  2. Ungeklärter Mechanismus: Trotz breiter Anwendung bleibt der detaillierte Katalysemechanismus unklar, insbesondere die Wirkungsweise von Sauerstoffatomen
  3. Experimentelle Beobachtungen: Operando-Raman-Experimente von Schumacher et al. zeigen eine signifikante Abnahme der Raman-Intensität in der Nähe von 782 cm−1^{-1} während des katalytischen Prozesses bei 500°C

Forschungsmotivation

  1. Bedarf an theoretischer Erklärung: Die beobachtete Raman-Intensitätsabnahme erfordert eine mikroskopische theoretische Erklärung
  2. Aufklärung des Katalysemechanismus: Verständnis der Rolle von Sauerstoffleerstellen im katalytischen Prozess auf atomarer Ebene
  3. Struktur-Eigenschafts-Beziehung: Etablierung einer quantitativen Beziehung zwischen Sauerstoffdefekten und Raman-Spektralveränderungen

Kernbeiträge

  1. Enthüllung des mikroskopischen Ursprungs der Raman-Intensitätsabnahme: DFT-Berechnungen beweisen, dass Sauerstoffleerstellen die Hauptursache für die experimentell beobachtete Raman-Intensitätsabnahme sind
  2. Etablierung der Struktur-Spektral-Beziehung: Identifikation der Raman-Linie in der Nähe von 782 cm−1^{-1} als sauerstoffdominierte asymmetrische MoO4_4-Streckschwingungsmode
  3. Entwicklung der effektiven gefrorenen Phonon-Methode: Theoretische Methode zur selektiven Fixierung spezifischer Atom-Phononenmoden zur Analyse von Raman-Intensitätsänderungen
  4. Erklärung experimenteller Phänomene: Theoretische Berechnungsergebnisse stimmen hochgradig mit experimentellen Beobachtungen überein und erklären, warum während des katalytischen Prozesses keine Peakverschiebung oder Verbreiterung beobachtet wird

Methodische Details

Berechnungsmethoden

Dichtefunktionaltheorie-Berechnungen (DFT):

  • Verwendung des Vienna Ab initio Simulation Package (VASP)
  • Projector Augmented Wave (PAW)-Methode
  • Austausch-Korrelations-Funktional: PBESol (modifizierte verallgemeinerte Gradienten-Näherung)
  • Hubbard-U-Korrektur: U = 4 eV für Fe 3d-Orbitale
  • Ebene-Wellen-Abschneidungsenergie: 500 eV
  • k-Punkt-Gitter: 2×3×2 für monoklin, 3×4×4 für orthorhombisch

Phononenberechnung

  1. Methode der endlichen Verschiebungen: Berechnung der Kraftkonstanten-Matrix mit 0,01 Å Atomverschiebung
  2. Dynamische Matrix: Konstruktion mittels phonopy-Softwarepaket zur Gewinnung von Phononendispersion und Zustandsdichte
  3. Raman-Intensitätsberechnung: Berechnung der Raman-Intensität durch dielektrische Tensoren in Kombination mit Phononen-Eigenvektoren

Innovative technische Methoden

Effektive gefrorene Phonon-Methode:

  • Selektive Fixierung von Phononenschwingungsmoden spezifischer Atome
  • Beibehaltung der Kristallsymmetrie
  • Bestimmung der Atombeiträge durch Vergleich der Raman-Intensität vor und nach dem Fixieren
  • Diese Methode berücksichtigt keine Ladungseffekte, kann aber Phononenbeiträge effektiv trennen

Experimentelle Einrichtung

Strukturmodelle

  1. Tieftemperatur-Monoklin-Phase: Raumgruppe P2₁/c, Z=8
    • Gitterparameter: a=15,869 Å, b=9,373 Å, c=15,868 Å, β=108,36°
  2. Hochtemperatur-Orthorhombische Phase: Raumgruppe Pbcn, Z=4
    • Gitterparameter: a=12,867 Å, b=9,287 Å, c=9,373 Å

Berechnungsparameter

  • Spinpolarisation: Einfache ferromagnetische Anordnung zur Erzielung einer Halbleiter-Bandlücke (~2,7 eV)
  • Strukturoptimierung: Kraftkonvergenzkriterium < 1 meV/Å
  • Gaußsche Verbreiterung: 0,05 eV
  • Raman-Spektral-Verbreiterung: Lorentzsche Verbreiterung mit Halbwertsbreite von 10 cm−1^{-1}

Vergleichende Analyse

Vergleich der Strukturunterschiede zwischen monokliner und orthorhombischer Phase durch Symmetriemodus-Analyse mit dem AMPLIMODES-Programm zur Bestimmung der Transformationsmatrix der Phasenumwandlung.

Experimentelle Ergebnisse

Phononenstrukturanalyse

  1. Frequenztrennung: Phononenspektrum unterteilt sich in zwei Regionen
    • Niederfrequenzbereich (0-403 cm−1^{-1}): Beiträge von Fe, Mo und O-Atomen
    • Hochfrequenzbereich (780-1001 cm−1^{-1}): Fe-Beitrag verschwindet, hauptsächlich Mo-O-Schwingungen
  2. Kritische Schwingungsmoden:
    • 789 cm−1^{-1}: Asymmetrische MoO4_4-Streckschwingungsmode mit maximaler Raman-Intensität
    • 972 cm−1^{-1}: Symmetrische MoO4_4-Streckschwingungsmode, Mo-Atom bewegt sich kaum

Raman-Intensitätsberechnungsergebnisse

Frequenz (cm−1^{-1})Irreduzible DarstellungNormalisierte Intensität (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

Ergebnisse der gefrorenen Phononenanalyse

Durch selektives Fixieren von Phononenmoden verschiedener Atome:

  • O4-Atom fixiert: Raman-Intensitätsabnahme von 8,5% bei 789,68 cm−1^{-1}
  • O5-Atom fixiert: Raman-Intensitätsabnahme von 7,6% bei 786,18 cm−1^{-1}
  • Mo-Atom fixiert: Intensitätsabnahme etwa halb so groß wie bei Sauerstoffatomen

Strukturrelaxation mit Sauerstoffleerstellen

Die Strukturrelaxation mit Sauerstoffleerstellen zeigt:

  • Dreieckige MoO3_3-Einheiten tendieren zur Wiederherstellung tetraedrischer MoO4_4-Geometrie
  • Fünffach koordinierte FeO5_5-Einheiten bleiben stabil
  • Signifikante Änderungen der Gitterparameter, Raumgruppe sinkt auf P1

Verwandte Arbeiten

Multifunktionalität von Fe2_2(MoO4_4)3_3

  1. Magnetische Eigenschaften: Antiferromagnetische Umwandlung vom L-Typ unterhalb von 12 K, magnetfeldinduzierte multiferroische Phase
  2. Elektrochemische Anwendungen: Als Elektrodenmaterial für Li/Na-Ionen-Batterien mit guter Zyklenstabilität
  3. Katalytische Anwendungen: Industrieller Katalysator für Methanol-zu-Formaldehyd-Umwandlung, Oberflächenkonzentration von Mo beeinflusst Selektivität

Operando-Spektroskopie-Studien

Bahnbrechende Arbeiten von Schumacher et al. führten erstmals das Konzept des Volumen-Sauerstoffspeichers ein und bilden die experimentelle Grundlage für diese Studie.

Theoretische Berechnungsmethoden

Anwendung verwandter DFT+U-Methoden bei der Behandlung elektronischer Struktur und Phononeneigenschaften von Übergangsmetalloxiden.

Schlussfolgerungen und Diskussion

Hauptschlussfolgerungen

  1. Sauerstoffleerstellen sind entscheidend: Sauerstoffleerstellen sind die Hauptursache für die Raman-Intensitätsabnahme in der Nähe von 782 cm−1^{-1}
  2. Identifikation der Schwingungsmoden: Der Frequenzbereich entspricht sauerstoffdominierten asymmetrischen MoO4_4-Streckschwingungen
  3. Beibehaltung der Symmetrie: Der Sauerstoffdiffusionsprozess verändert die lokale Symmetrie nicht wesentlich, was erklärt, warum im Experiment keine Peakverschiebung oder Verbreiterung beobachtet wird
  4. Katalysemechanismus: Unterstützt den Katalysemechanismus der Volumenphase als Sauerstoffspeicher

Einschränkungen

  1. Vereinfachtes Modell: Die gefrorene Phonon-Methode berücksichtigt keine Ladungseffekte und Fernwechselwirkungen
  2. Einzelne Leerstelle: Berücksichtigung nur einer einzelnen Sauerstoffleerstelle, während der tatsächliche katalytische Prozess mehrere Leerstellen oder Leerstellen-Cluster beinhalten könnte
  3. Temperatureffekte: Berechnungen basieren auf 0 K, berücksichtigen keine Wärmeschwingungseffekte bei endlicher Temperatur
  4. Oberflächeneffekte: Hauptfokus auf Volumeneigenschaften mit begrenzter Beschreibung von Oberflächenkatalyse-Aktivitätsstellen

Zukünftige Richtungen

  1. Mehrfach-Leerstellen-Forschung: Untersuchung des Einflusses von Leerstellenkonzentration und -verteilung auf das Raman-Spektrum
  2. Dynamische Prozesse: Berechnung von Sauerstoff-Diffusionspfaden und Energiebarrieren
  3. Oberflächenmodellierung: Etablierung von Oberflächenmodellen zur Untersuchung katalytischer Aktivitätsstellen
  4. Experimentelle Validierung: Empfehlung zur Verifizierung von Fe-Oxidationszustandsänderungen durch Mössbauer-Spektroskopie und andere Techniken

Tiefgreifende Bewertung

Stärken

  1. Theorie-Experiment-Kombination: Gute Erklärung der experimentell beobachteten Raman-Spektralveränderungen
  2. Methodische Innovation: Die gefrorene Phonon-Methode bietet ein neues Werkzeug zur Analyse von Raman-Intensitätsänderungen
  3. Rigorose Berechnung: DFT+U-Methode behandelt angemessen stark korrelierte Elektronensysteme
  4. Klares physikalisches Bild: Aufklärung des Wirkmechanismus von Sauerstoffleerstellen auf atomarer Ebene

Mängel

  1. Modellvereinfachung: Die gefrorene Phonon-Methode ignoriert Auswirkungen von Elektronenstrukturänderungen
  2. Leerstellenkonzentration: Keine systematische Untersuchung der Auswirkungen verschiedener Leerstellenkonzentrationen
  3. Fehlende Dynamik-Informationen: Mangel an kinetischer Analyse der Sauerstoffatom-Migration
  4. Begrenzte experimentelle Vergleiche: Validierung hauptsächlich basierend auf einem Satz experimenteller Daten

Auswirkungen

  1. Katalysetheorie: Bietet wichtige theoretische Unterstützung für das Verständnis von Oxidations-Dehydrierungs-Katalysemechanismen
  2. Spektralanalyse: Bietet theoretische Grundlage für die Anwendung der Raman-Spektroskopie in der Katalyseforschung
  3. Materialdesign: Bietet Designideen für die Entwicklung neuer Oxidkatalysa­toren
  4. Methodologischer Beitrag: Die gefrorene Phonon-Methode kann auf andere Materialsysteme verallgemeinert werden

Anwendungsszenarien

  1. Oxidkatalysa­toren: Anwendbar auf die Untersuchung von Sauerstoffleerstellen-Effekten in verschiedenen Oxidkatalysa­toren
  2. Operando-Spektral-Interpretation: Bietet Methoden zur theoretischen Interpretation von operando-Raman-Experimenten
  3. Defekt-Engineering: Leitet die Verbesserung der katalytischen Leistung durch Sauerstoffleerstellen-Regulierung an
  4. Heterogene Katalyse: Bietet theoretischen Rahmen zum Verständnis von Sauerstofftransportprozessen in heterogener Katalyse

Literaturverzeichnis

Dieses Papier zitiert 34 wichtige Literaturquellen, die die Struktur, magnetische Eigenschaften, elektrochemische und katalytische Eigenschaften von Fe2_2(MoO4_4)3_3 sowie verwandte theoretische Berechnungsmethoden und experimentelle Techniken abdecken und eine solide theoretische Grundlage für die Forschung bieten.

Gesamtbewertung: Dies ist ein hochqualitatives theoretisches Berechnungspapier, das erfolgreich das experimentell beobachtete anomale Raman-Spektral-Phänomen erklärt und wichtige theoretische Unterstützung für das Verständnis des Wirkmechanismus von Eisenmolybdat-Katalysatoren bietet. Die Methode ist innovativ, die Berechnung ist rigoros, die Schlussfolgerungen sind zuverlässig und das Papier hat wichtigen akademischen Wert und praktische Bedeutung.