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

Ammorbidimento Raman indotto da vacanze di ossigeno nel catalizzatore Fe2_2(MoO4_4)3_3

Informazioni Fondamentali

  • ID Articolo: 2510.12746
  • Titolo: Oxygen-vacancy-induced Raman softening in the catalyst Fe2_2(MoO4_4)3_3
  • Autori: Young-Joon Song, Roser Valentí (Goethe University Frankfurt)
  • Classificazione: cond-mat.mtrl-sci (Fisica della Materia Condensata - Scienza dei Materiali)
  • Data di Pubblicazione: 15 Ottobre 2025
  • Link Articolo: https://arxiv.org/abs/2510.12746

Riassunto

Il ferromolibdato (Fe2_2(MoO4_4)3_3) è ampiamente utilizzato nei processi catalitici industriali per l'ossidazione deidrogenativa del metanolo a formaldeide, tuttavia il suo meccanismo catalitico dettagliato rimane ancora poco chiaro. Recenti misurazioni Raman in situ e di impedenza suggeriscono che la regione di bulk del Fe2_2(MoO4_4)3_3 potrebbe fungere da serbatoio di atomi di ossigeno, che possono migrare verso la superficie per partecipare alle reazioni catalitiche. Questa conclusione si basa sulla significativa riduzione dell'intensità Raman osservata durante il processo catalitico, che suggerisce la formazione di difetti atomici. Tuttavia, l'origine microscopica di questa riduzione di intensità rimane ancora da chiarire. Il presente studio rivela l'origine dei cambiamenti osservati sperimentalmente nell'intensità Raman attraverso calcoli di teoria del funzionale della densità (DFT). L'analisi fonica mostra che i modi vibrazionali dominati dall'ossigeno (con un piccolo contributo del Mo) compaiono intorno a 782 cm−1^{-1}, esattamente nella regione di frequenza dove è stata misurata sperimentalmente la riduzione dell'intensità Raman. L'utilizzo del metodo del fonone congelato effettivo ha ulteriormente confermato che le vibrazioni dell'ossigeno sono la causa principale della riduzione calcolata dell'intensità Raman.

Contesto di Ricerca e Motivazione

Sfondo del Problema

  1. Importanza dell'Applicazione Catalitica: Il Fe2_2(MoO4_4)3_3 è un importante catalizzatore industriale per l'ossidazione deidrogenativa del metanolo a formaldeide, utilizzato da decenni
  2. Meccanismo Poco Chiaro: Nonostante l'ampio utilizzo, il meccanismo dettagliato del processo catalitico rimane poco compreso, in particolare il ruolo degli atomi di ossigeno
  3. Osservazioni Sperimentali: Gli esperimenti Raman in situ di Schumacher et al. hanno rivelato una significativa riduzione dell'intensità Raman intorno a 782 cm−1^{-1} durante il processo catalitico a 500°C

Motivazione della Ricerca

  1. Necessità di Spiegazione Teorica: La riduzione dell'intensità Raman osservata sperimentalmente manca di una spiegazione teorica microscopica
  2. Chiarimento del Meccanismo Catalitico: È necessario comprendere a livello atomico il ruolo delle vacanze di ossigeno nel processo catalitico
  3. Relazione Struttura-Proprietà: Stabilire una relazione quantitativa tra i difetti di ossigeno e i cambiamenti dello spettro Raman

Contributi Principali

  1. Rivelazione dell'Origine Microscopica dell'Ammorbidimento Raman: Attraverso calcoli DFT si dimostra che le vacanze di ossigeno sono la causa principale della riduzione dell'intensità Raman osservata sperimentalmente
  2. Stabilimento della Relazione Struttura-Spettro: Identificazione del picco Raman intorno a 782 cm−1^{-1} corrispondente al modo di stiramento asimmetrico del MoO4_4 dominato dall'ossigeno
  3. Proposizione del Metodo del Fonone Congelato Effettivo: Sviluppo di un metodo teorico per congelare selettivamente i modi fonici di atomi specifici al fine di analizzare i cambiamenti dell'intensità Raman
  4. Spiegazione dei Fenomeni Sperimentali: I risultati dei calcoli teorici sono altamente coerenti con le osservazioni sperimentali, spiegando perché non sono stati osservati spostamenti di picco o allargamenti durante il processo catalitico

Dettagli Metodologici

Metodi Computazionali

Calcoli di Teoria del Funzionale della Densità (DFT):

  • Utilizzo del Vienna Ab initio Simulation Package (VASP)
  • Metodo delle onde aumentate proiettate (PAW)
  • Funzionale di scambio-correlazione: PBESol (approssimazione del gradiente generalizzato modificato)
  • Correzione di Hubbard U: valore U di 4 eV applicato agli orbitali 3d del Fe
  • Energia di taglio delle onde piane: 500 eV
  • Griglia di punti k: 2×3×2 per la fase monoclina, 3×4×4 per la fase ortorombica

Calcoli Fonici

  1. Metodo dello Spostamento Finito: Calcolo della matrice delle costanti di forza utilizzando uno spostamento atomico di 0,01 Å
  2. Matrice Dinamica: Costruzione mediante il pacchetto software phonopy per ottenere la dispersione fonica e la densità di stati
  3. Calcolo dell'Intensità Raman: Calcolo dell'intensità Raman attraverso i risultati del tensore dielettrico combinati con gli autovettori fonici

Metodo Tecnico Innovativo

Metodo del Fonone Congelato Effettivo:

  • Congelamento selettivo dei modi di vibrazione fonica di atomi specifici
  • Mantenimento invariato della simmetria cristallina
  • Determinazione del contributo di ciascun atomo attraverso il confronto dei cambiamenti dell'intensità Raman prima e dopo il congelamento
  • Questo metodo non considera gli effetti di carica, ma è efficace nel separare i contributi fonici

Configurazione Sperimentale

Modelli Strutturali

  1. Fase Monoclina a Bassa Temperatura: Gruppo spaziale P2₁/c, Z=8
    • Parametri reticolari: a=15,869 Å, b=9,373 Å, c=15,868 Å, β=108,36°
  2. Fase Ortorombica ad Alta Temperatura: Gruppo spaziale Pbcn, Z=4
    • Parametri reticolari: a=12,867 Å, b=9,287 Å, c=9,373 Å

Parametri Computazionali

  • Polarizzazione di Spin: Utilizzo di un semplice ordinamento ferromagnetico per ottenere un gap a semiconduttore (~2,7 eV)
  • Ottimizzazione Strutturale: Criterio di convergenza delle forze < 1 meV/Å
  • Allargamento Gaussiano: 0,05 eV
  • Allargamento dello Spettro Raman: Allargamento Lorentziano, larghezza a metà altezza di 10 cm−1^{-1}

Analisi Comparativa

Confronto delle differenze strutturali tra la fase monoclina e ortorombica attraverso l'analisi dei modi di simmetria, utilizzando il programma AMPLIMODES per determinare la matrice di trasformazione della transizione di fase.

Risultati Sperimentali

Analisi della Struttura Fonica

  1. Separazione in Frequenza: Lo spettro fonico è diviso in due regioni
    • Regione a Bassa Frequenza (0-403 cm−1^{-1}): Contributi di Fe, Mo e O
    • Regione ad Alta Frequenza (780-1001 cm−1^{-1}): Scomparsa del contributo del Fe, principalmente vibrazioni Mo-O
  2. Modi Vibrazionali Chiave:
    • 789 cm−1^{-1}: Modo di stiramento asimmetrico del MoO4_4, massima intensità Raman
    • 972 cm−1^{-1}: Modo di stiramento simmetrico del MoO4_4, l'atomo di Mo rimane quasi immobile

Risultati del Calcolo dell'Intensità Raman

Frequenza (cm−1^{-1})Rappresentazione IrriducibileIntensità Normalizzata (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

Risultati dell'Analisi del Fonone Congelato

Attraverso il congelamento selettivo dei modi fonici di diversi atomi:

  • Congelamento dell'Atomo O4: Riduzione dell'intensità Raman dell'8,5% a 789,68 cm−1^{-1}
  • Congelamento dell'Atomo O5: Riduzione dell'intensità Raman del 7,6% a 786,18 cm−1^{-1}
  • Congelamento dell'Atomo Mo: L'entità della riduzione di intensità è circa la metà di quella degli atomi di ossigeno

Rilassamento Strutturale con Vacanza di Ossigeno

Il rilassamento strutturale contenente una vacanza di ossigeno mostra:

  • L'unità triangolare MoO3_3 tende a recuperare la geometria tetraedrica del MoO4_4
  • L'unità pentacoordinata FeO5_5 rimane stabile
  • Cambiamenti significativi nei parametri reticolari, riduzione del gruppo spaziale a P1

Lavori Correlati

Proprietà Multifunzionali del Fe2_2(MoO4_4)3_3

  1. Proprietà Magnetiche: Transizione ferrimagnetica di tipo L al di sotto di 12 K, fase multiferroica indotta da campo magnetico esterno
  2. Applicazioni Elettrochimiche: Utilizzo come materiale di elettrodo per batterie Li/Na ioni, con buone prestazioni di ciclabilità
  3. Applicazioni Catalitiche: Catalizzatore industriale per la sintesi di formaldeide dal metanolo, la concentrazione superficiale di Mo influenza la selettività

Studi Spettroscopici In Situ

Il lavoro pioneristico di Schumacher et al. ha proposto per la prima volta il concetto di serbatoio di ossigeno nel bulk, fornendo la base sperimentale per il presente studio.

Metodi di Calcolo Teorico

Applicazioni del metodo DFT+U nel trattamento della struttura elettronica e delle proprietà fonici degli ossidi di metalli di transizione.

Conclusioni e Discussione

Conclusioni Principali

  1. Le Vacanze di Ossigeno sono Cruciali: Le vacanze di ossigeno sono la causa principale della riduzione dell'intensità Raman osservata intorno a 782 cm−1^{-1}
  2. Identificazione del Modo Vibrazionale: La regione di frequenza corrisponde al modo di stiramento asimmetrico del MoO4_4 dominato dall'ossigeno
  3. Mantenimento della Simmetria: Il processo di diffusione dell'ossigeno non modifica significativamente la simmetria locale, spiegando perché non sono stati osservati spostamenti o allargamenti di picco negli esperimenti
  4. Meccanismo Catalitico: Supporta il meccanismo catalitico del bulk come serbatoio di ossigeno

Limitazioni

  1. Modello Semplificato: Il metodo del fonone congelato non considera gli effetti di carica e le interazioni a lungo raggio
  2. Vacanza Singola: Considera solo una singola vacanza di ossigeno, mentre i processi catalitici reali potrebbero coinvolgere più vacanze o cluster di vacanze
  3. Effetti di Temperatura: I calcoli si basano su 0 K, senza considerare gli effetti delle vibrazioni termiche a temperatura finita
  4. Effetti di Superficie: Si concentra principalmente sulle proprietà del bulk, con descrizione limitata dei siti attivi catalitici sulla superficie

Direzioni Future

  1. Studio di Vacanze Multiple: Ricerca dell'influenza della concentrazione e della distribuzione delle vacanze sullo spettro Raman
  2. Processi Dinamici: Calcolo dei percorsi di diffusione dell'ossigeno e delle barriere energetiche
  3. Modellazione della Superficie: Costruzione di modelli di superficie per lo studio dei siti attivi catalitici
  4. Verifica Sperimentale: Suggerimento di verificare i cambiamenti dello stato di ossidazione del Fe attraverso tecniche come la spettroscopia di Mössbauer

Valutazione Approfondita

Punti di Forza

  1. Integrazione Teoria-Esperimento: Spiega efficacemente i cambiamenti dello spettro Raman osservati sperimentalmente
  2. Innovazione Metodologica: Il metodo del fonone congelato fornisce un nuovo strumento per analizzare i cambiamenti dell'intensità Raman
  3. Rigore Computazionale: Il metodo DFT+U affronta adeguatamente il sistema elettronico fortemente correlato
  4. Immagine Fisica Chiara: Chiarisce a livello atomico il meccanismo d'azione delle vacanze di ossigeno

Insufficienze

  1. Semplificazione del Modello: Il metodo del fonone congelato trascura gli effetti dei cambiamenti della struttura elettronica
  2. Concentrazione di Vacanze: Manca uno studio sistematico dell'influenza di diverse concentrazioni di vacanze
  3. Informazioni Dinamiche Mancanti: Manca l'analisi della dinamica della migrazione degli atomi di ossigeno
  4. Confronto Sperimentale Limitato: La verifica si basa principalmente su un singolo set di dati sperimentali

Impatto

  1. Teoria Catalitica: Fornisce un importante supporto teorico per la comprensione del meccanismo catalitico dell'ossidazione deidrogenativa
  2. Analisi Spettroscopica: Fornisce una base teorica per l'applicazione della spettroscopia Raman nella ricerca catalitica
  3. Progettazione di Materiali: Fornisce idee di progettazione per lo sviluppo di nuovi catalizzatori a base di ossidi
  4. Contributo Metodologico: Il metodo del fonone congelato può essere esteso ad altri sistemi di materiali

Scenari di Applicabilità

  1. Catalizzatori a Base di Ossidi: Applicabile allo studio degli effetti delle vacanze di ossigeno in vari catalizzatori a base di ossidi
  2. Interpretazione di Spettroscopia In Situ: Fornisce metodi per l'interpretazione teorica dei dati di spettroscopia Raman in situ
  3. Ingegneria dei Difetti: Guida il miglioramento delle prestazioni catalitiche attraverso la regolazione delle vacanze di ossigeno
  4. Catalisi Eterogenea: Fornisce un quadro teorico per la comprensione dei processi di trasporto dell'ossigeno nella catalisi eterogenea

Bibliografia

L'articolo cita 34 importanti riferimenti che coprono le proprietà strutturali, magnetiche, elettrochimiche e catalitiche del Fe2_2(MoO4_4)3_3, nonché i metodi di calcolo teorico correlati e le tecniche sperimentali, fornendo una base teorica solida per la ricerca.

Valutazione Complessiva: Questo è un articolo di calcolo teorico di alta qualità che spiega con successo l'anomalia dello spettro Raman osservata sperimentalmente, fornendo un importante supporto teorico per la comprensione del meccanismo di funzionamento del catalizzatore ferromolibdato. Il metodo è innovativo, il calcolo è rigoroso, le conclusioni sono affidabili e possiede un importante valore accademico e significato pratico.