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

Oxygen-vacancy-induced Raman softening in the catalyst Fe2_2(MoO4_4)3_3

Basic Information

  • Paper ID: 2510.12746
  • Title: Oxygen-vacancy-induced Raman softening in the catalyst Fe2_2(MoO4_4)3_3
  • Authors: Young-Joon Song, Roser Valentí (Goethe University Frankfurt)
  • Classification: cond-mat.mtrl-sci (Condensed Matter Physics - Materials Science)
  • Publication Date: October 15, 2025
  • Paper Link: https://arxiv.org/abs/2510.12746

Abstract

Iron molybdate (Fe2_2(MoO4_4)3_3) has been widely employed in the industrial catalytic process of methanol oxidative dehydrogenation to formaldehyde, yet its detailed catalytic mechanism remains unclear. Recent in-situ Raman and impedance measurements suggest that the bulk region of Fe2_2(MoO4_4)3_3 may function as an oxygen reservoir, with oxygen atoms capable of migrating to the surface to participate in catalytic reactions. This conclusion is based on the significant reduction in Raman intensity observed during the catalytic process, implying the formation of atomic defects. However, the microscopic origin of this intensity reduction awaits clarification. This study employs density functional theory (DFT) calculations to elucidate the origin of the experimentally observed Raman intensity changes. Phonon analysis reveals that oxygen-dominated vibrational modes (with minor Mo contributions) appear near 782 cm−1^{-1}, precisely in the frequency region where experimental measurements detected reduced Raman intensity. The effective frozen phonon method further confirms that oxygen vibrations are the primary cause of the calculated Raman intensity reduction.

Research Background and Motivation

Problem Background

  1. Catalytic Application Significance: Fe2_2(MoO4_4)3_3 is an important industrial catalyst for methanol oxidative dehydrogenation to formaldehyde, employed for decades
  2. Unclear Mechanism: Despite widespread application, the detailed mechanism of the catalytic process remains obscure, particularly the role of oxygen atoms
  3. Experimental Observations: In-situ Raman experiments by Schumacher et al. revealed a significant reduction in Raman intensity near 782 cm−1^{-1} during catalytic processes at 500°C

Research Motivation

  1. Need for Theoretical Explanation: The experimentally observed Raman intensity reduction lacks microscopic theoretical interpretation
  2. Elucidation of Catalytic Mechanism: Understanding the role of oxygen vacancies in catalysis at the atomic level is essential
  3. Structure-Property Relationship: Establishing quantitative relationships between oxygen defects and Raman spectral changes

Core Contributions

  1. Revealing the Microscopic Origin of Raman Intensity Reduction: DFT calculations demonstrate that oxygen vacancies are the primary cause of the experimentally observed Raman intensity reduction
  2. Establishing Structure-Spectral Relationships: Identifying that the Raman peak near 782 cm−1^{-1} corresponds to oxygen-dominated asymmetric MoO4_4 stretching vibrational modes
  3. Proposing the Effective Frozen Phonon Method: Developing a theoretical approach to selectively freeze phonon modes of specific atoms for analyzing Raman intensity changes
  4. Explaining Experimental Phenomena: Theoretical calculations are highly consistent with experimental observations, explaining why peak shifts or broadening were not observed during catalysis

Methodology Details

Computational Methods

Density Functional Theory (DFT) Calculations:

  • Vienna Ab initio Simulation Package (VASP)
  • Projector augmented wave (PAW) method
  • Exchange-correlation functional: PBESol (revised generalized gradient approximation)
  • Hubbard U correction: 4 eV applied to Fe 3d orbitals
  • Plane wave cutoff energy: 500 eV
  • k-point mesh: 2×3×2 for monoclinic phase, 3×4×4 for orthorhombic phase

Phonon Calculations

  1. Finite Displacement Method: Force constant matrix calculated using 0.01 Å atomic displacements
  2. Dynamical Matrix: Constructed via phonopy software package to obtain phonon dispersions and density of states
  3. Raman Intensity Calculation: Computed from dielectric tensor results combined with phonon eigenvectors

Innovative Technical Approach

Effective Frozen Phonon Method:

  • Selective freezing of phonon vibrational modes for specific atoms
  • Preservation of crystal symmetry
  • Determination of atomic contributions by comparing Raman intensities before and after freezing
  • Method neglects charge effects but effectively isolates phonon contributions

Experimental Setup

Structural Models

  1. Low-Temperature Monoclinic Phase: P2₁/c space group, Z=8
    • Lattice parameters: a=15.869 Å, b=9.373 Å, c=15.868 Å, β=108.36°
  2. High-Temperature Orthorhombic Phase: Pbcn space group, Z=4
    • Lattice parameters: a=12.867 Å, b=9.287 Å, c=9.373 Å

Computational Parameters

  • Spin Polarization: Simple ferromagnetic arrangement yielding semiconductor band gap (~2.7 eV)
  • Structure Optimization: Force convergence criterion < 1 meV/Å
  • Gaussian Broadening: 0.05 eV
  • Raman Spectral Broadening: Lorentzian broadening with half-width of 10 cm−1^{-1}

Comparative Analysis

Symmetry mode analysis comparing monoclinic and orthorhombic phases using AMPLIMODES to determine transformation matrices for phase transitions.

Experimental Results

Phonon Structure Analysis

  1. Frequency Separation: Phonon spectrum divided into two regions
    • Low-frequency region (0-403 cm−1^{-1}): Contributions from Fe, Mo, and O atoms
    • High-frequency region (780-1001 cm−1^{-1}): Fe contribution vanishes, dominated by Mo-O vibrations
  2. Key Vibrational Modes:
    • 789 cm−1^{-1}: Asymmetric MoO4_4 stretching mode with maximum Raman intensity
    • 972 cm−1^{-1}: Symmetric MoO4_4 stretching mode with negligible Mo atomic motion

Raman Intensity Calculation Results

Frequency (cm−1^{-1})Irreducible RepresentationNormalized Intensity (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

Frozen Phonon Analysis Results

Through selective freezing of phonon modes for different atoms:

  • O4 Atom Frozen: Raman intensity reduction of 8.5% at 789.68 cm−1^{-1}
  • O5 Atom Frozen: Raman intensity reduction of 7.6% at 786.18 cm−1^{-1}
  • Mo Atom Frozen: Intensity reduction approximately half that of oxygen atoms

Oxygen Vacancy Structure Relaxation

Structure relaxation with oxygen vacancies reveals:

  • Triangular MoO3_3 units tend to recover tetrahedral MoO4_4 geometry
  • Five-coordinate FeO5_5 units remain stable
  • Significant lattice parameter changes with space group reduction to P1

Related Work

Multifunctional Properties of Fe2_2(MoO4_4)3_3

  1. Magnetic Properties: L-type ferrimagnetic transition below 12 K, magnetic field-induced multiferroic phase
  2. Electrochemical Applications: Electrode material for Li/Na-ion batteries with excellent cycling performance
  3. Catalytic Applications: Industrial methanol-to-formaldehyde catalyst with surface Mo concentration affecting selectivity

In-Situ Spectroscopic Studies

Pioneering work by Schumacher et al. first proposed the bulk oxygen reservoir concept, providing experimental foundation for this study.

Theoretical Computational Methods

Applications of DFT+U methods in treating electronic structure and phonon properties of transition metal oxides.

Conclusions and Discussion

Main Conclusions

  1. Oxygen Vacancies are Key: Oxygen vacancies are the primary cause of Raman intensity reduction near 782 cm−1^{-1}
  2. Vibrational Mode Identification: This frequency region corresponds to oxygen-dominated asymmetric MoO4_4 stretching vibrations
  3. Symmetry Preservation: Oxygen diffusion does not significantly alter local symmetry, explaining the absence of observed peak shifts or broadening in experiments
  4. Catalytic Mechanism: Supports the catalytic mechanism of bulk phase functioning as oxygen reservoir

Limitations

  1. Simplified Model: Frozen phonon method neglects charge effects and long-range interactions
  2. Single Vacancy: Only considers isolated oxygen vacancies; actual catalytic processes may involve multiple vacancies or vacancy clusters
  3. Temperature Effects: Calculations based on 0 K without considering thermal vibrations at finite temperatures
  4. Surface Effects: Primarily focuses on bulk properties with limited description of surface catalytic active sites

Future Directions

  1. Multi-Vacancy Studies: Investigating effects of vacancy concentration and distribution on Raman spectra
  2. Kinetic Processes: Computing oxygen atom diffusion pathways and energy barriers
  3. Surface Modeling: Establishing surface models to study catalytic active sites
  4. Experimental Verification: Recommending verification of Fe oxidation state changes through Mössbauer spectroscopy

In-Depth Evaluation

Strengths

  1. Theory-Experiment Integration: Effectively explains experimentally observed Raman spectral changes
  2. Methodological Innovation: Frozen phonon method provides new tools for analyzing Raman intensity variations
  3. Computational Rigor: DFT+U method appropriately handles strongly correlated electron systems
  4. Clear Physical Picture: Elucidates oxygen vacancy mechanisms at the atomic level

Weaknesses

  1. Model Simplification: Frozen phonon method neglects electronic structure changes
  2. Vacancy Concentration: Lacks systematic investigation of different vacancy concentration effects
  3. Missing Kinetic Information: Absence of kinetic analysis for oxygen atom migration
  4. Limited Experimental Comparison: Verification primarily based on single experimental dataset

Impact

  1. Catalytic Theory: Provides important theoretical support for understanding oxidative dehydrogenation catalytic mechanisms
  2. Spectroscopic Analysis: Establishes theoretical foundation for Raman spectroscopy applications in catalysis research
  3. Materials Design: Offers design strategies for developing novel oxide catalysts
  4. Methodological Contribution: Frozen phonon method generalizable to other material systems

Applicable Scenarios

  1. Oxide Catalysts: Applicable to studying oxygen vacancy effects in various oxide catalysts
  2. In-Situ Spectral Interpretation: Provides methodology for theoretical interpretation of in-situ Raman experimental data
  3. Defect Engineering: Guides performance improvement through oxygen vacancy engineering
  4. Heterogeneous Catalysis: Provides theoretical framework for understanding oxygen transport in heterogeneous catalysis

References

This paper cites 34 important references covering the structure, magnetic properties, electrochemistry, and catalytic properties of Fe2_2(MoO4_4)3_3, as well as relevant theoretical computational methods and experimental techniques, providing a solid theoretical foundation for the research.

Overall Assessment: This is a high-quality theoretical computational paper that successfully explains anomalous Raman spectral phenomena observed experimentally, providing important theoretical support for understanding the working mechanism of iron molybdate catalysts. The methodology is innovative, calculations rigorous, conclusions reliable, and the work possesses significant academic value and practical significance.