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

基本信息

  • 论文ID: 2510.12746
  • 标题: Oxygen-vacancy-induced Raman softening in the catalyst Fe2_2(MoO4_4)3_3
  • 作者: Young-Joon Song, Roser Valentí (Goethe University Frankfurt)
  • 分类: cond-mat.mtrl-sci (凝聚态物理-材料科学)
  • 发表时间: October 15, 2025
  • 论文链接: https://arxiv.org/abs/2510.12746

摘要

铁钼酸盐(Fe2_2(MoO4_4)3_3)已被广泛用于甲醇氧化脱氢制甲醛的工业催化过程，但其详细的催化机理仍不清楚。最近的原位拉曼和阻抗测量表明，Fe2_2(MoO4_4)3_3的体相区域可能作为氧原子储库，氧原子可以迁移到表面参与催化反应。这一结论基于催化过程中观察到的拉曼强度显著降低，暗示形成了原子缺陷。然而，这种强度降低的微观起源尚待阐明。本研究通过密度泛函理论(DFT)计算揭示了实验观察到的拉曼强度变化的起源。声子分析表明，氧主导的振动模式(伴有少量Mo贡献)出现在782 cm−1^{-1}附近，正是实验测量到拉曼强度降低的频率区域。使用有效冻结声子方法进一步证明了氧振动是计算拉曼强度降低的主要原因。

研究背景与动机

问题背景

  1. 催化应用重要性: Fe2_2(MoO4_4)3_3是甲醇氧化脱氢制甲醛的重要工业催化剂，已使用数十年
  2. 机理不明: 尽管应用广泛，但催化过程的详细机理仍不清楚，特别是氧原子的作用机制
  3. 实验观察: Schumacher等人的原位拉曼实验发现，在500°C催化过程中，782 cm−1^{-1}附近的拉曼强度显著降低

研究动机

  1. 理论解释需求: 实验观察到的拉曼强度降低缺乏微观理论解释
  2. 催化机理阐明: 需要从原子层面理解氧空位在催化过程中的作用
  3. 结构-性能关系: 建立氧缺陷与拉曼光谱变化的定量关系

核心贡献

  1. 揭示拉曼强度降低的微观起源: 通过DFT计算证明氧空位是导致实验观察到的拉曼强度降低的主要原因
  2. 建立结构-光谱关系: 确定了782 cm−1^{-1}附近的拉曼峰对应于氧主导的不对称MoO4_4伸缩振动模式
  3. 提出有效冻结声子方法: 开发了选择性冻结特定原子声子模式的理论方法来分析拉曼强度变化
  4. 解释实验现象: 理论计算结果与实验观察高度一致，解释了为什么催化过程中没有观察到峰位移或展宽

方法详解

计算方法

密度泛函理论(DFT)计算:

  • 使用Vienna Ab initio Simulation Package (VASP)
  • 采用投影缀加波(PAW)方法
  • 交换关联泛函: PBESol (修正广义梯度近似)
  • Hubbard U校正: 对Fe 3d轨道施加4 eV的U值
  • 平面波截断能: 500 eV
  • k点网格: 单斜相2×3×2，正交相3×4×4

声子计算

  1. 有限位移法: 使用0.01 Å的原子位移计算力常数矩阵
  2. 动力学矩阵: 通过phonopy软件包构建并获得声子色散和态密度
  3. 拉曼强度计算: 通过介电张量结果结合声子本征矢量计算拉曼强度

创新技术方法

有效冻结声子方法:

  • 选择性冻结特定原子的声子振动模式
  • 保持晶体对称性不变
  • 通过比较冻结前后的拉曼强度变化来确定各原子的贡献
  • 该方法不考虑电荷效应，但能有效分离声子贡献

实验设置

结构模型

  1. 低温单斜相: P2₁/c空间群，Z=8
    • 晶格参数: a=15.869 Å, b=9.373 Å, c=15.868 Å, β=108.36°
  2. 高温正交相: Pbcn空间群，Z=4
    • 晶格参数: a=12.867 Å, b=9.287 Å, c=9.373 Å

计算参数

  • 自旋极化: 采用简单铁磁排列获得半导体能隙(~2.7 eV)
  • 结构优化: 力收敛标准 < 1 meV/Å
  • 高斯展宽: 0.05 eV
  • 拉曼谱展宽: 洛伦兹展宽，半宽10 cm−1^{-1}

对比分析

通过对称模式分析比较单斜相和正交相的结构差异，使用AMPLIMODES程序确定相变的变换矩阵。

实验结果

声子结构分析

  1. 频率分离: 声子谱分为两个区域
    • 低频区(0-403 cm−1^{-1}): Fe, Mo, O原子均有贡献
    • 高频区(780-1001 cm−1^{-1}): Fe贡献消失，主要为Mo-O振动
  2. 关键振动模式:
    • 789 cm−1^{-1}: 不对称MoO4_4伸缩模式，拉曼强度最大
    • 972 cm−1^{-1}: 对称MoO4_4伸缩模式，Mo原子几乎不动

拉曼强度计算结果

频率 (cm−1^{-1})不可约表示归一化强度 (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

冻结声子分析结果

通过选择性冻结不同原子的声子模式：

  • O4原子冻结: 在789.68 cm−1^{-1}处拉曼强度降低8.5%
  • O5原子冻结: 在786.18 cm−1^{-1}处拉曼强度降低7.6%
  • Mo原子冻结: 强度降低幅度约为氧原子的一半

氧空位结构弛豫

含氧空位的结构弛豫显示：

  • 三角MoO3_3单元趋向恢复四面体MoO4_4几何
  • 五配位FeO5_5单元保持稳定
  • 晶格参数变化显著，空间群降为P1

相关工作

Fe2_2(MoO4_4)3_3的多功能性质

  1. 磁学性质: 12 K以下L型亚铁磁转变，外磁场诱导多铁相
  2. 电化学应用: 作为Li/Na离子电池电极材料，具有良好的循环性能
  3. 催化应用: 工业甲醇制甲醛催化剂，表面Mo浓度影响选择性

原位光谱研究

Schumacher等人的开创性工作首次提出体相氧储库概念，为本研究提供了实验基础。

理论计算方法

相关的DFT+U方法在处理过渡金属氧化物的电子结构和声子性质方面的应用。

结论与讨论

主要结论

  1. 氧空位是关键: 氧空位是导致782 cm−1^{-1}附近拉曼强度降低的主要原因
  2. 振动模式识别: 该频率区域对应氧主导的不对称MoO4_4伸缩振动
  3. 对称性保持: 氧扩散过程不显著改变局域对称性，解释了实验中未观察到峰位移或展宽
  4. 催化机理: 支持体相作为氧储库的催化机理

局限性

  1. 简化模型: 冻结声子方法不考虑电荷效应和长程相互作用
  2. 单一空位: 仅考虑单个氧空位，实际催化过程可能涉及多空位或空位团簇
  3. 温度效应: 计算基于0 K，未考虑有限温度下的热振动效应
  4. 表面效应: 主要关注体相性质，对表面催化活性位点的描述有限

未来方向

  1. 多空位研究: 研究空位浓度和分布对拉曼光谱的影响
  2. 动力学过程: 计算氧原子扩散路径和能垒
  3. 表面建模: 建立表面模型研究催化活性位点
  4. 实验验证: 建议通过穆斯堡尔谱等技术验证Fe氧化态变化

深度评价

优点

  1. 理论与实验结合: 很好地解释了实验观察到的拉曼光谱变化
  2. 方法创新: 冻结声子方法为分析拉曼强度变化提供了新工具
  3. 计算严谨: DFT+U方法合理处理了强关联电子体系
  4. 物理图像清晰: 从原子层面阐明了氧空位的作用机制

不足

  1. 模型简化: 冻结声子方法忽略了电子结构变化的影响
  2. 空位浓度: 未系统研究不同空位浓度的影响
  3. 动力学信息缺失: 缺乏氧原子迁移的动力学分析
  4. 实验对比有限: 主要基于一组实验数据进行验证

影响力

  1. 催化理论: 为理解氧化脱氢催化机理提供了重要理论支撑
  2. 光谱分析: 为拉曼光谱在催化研究中的应用提供了理论基础
  3. 材料设计: 为开发新型氧化物催化剂提供了设计思路
  4. 方法学贡献: 冻结声子方法可推广到其他材料体系

适用场景

  1. 氧化物催化剂: 适用于研究各种氧化物催化剂中的氧空位效应
  2. 原位光谱解释: 为原位拉曼实验数据的理论解释提供方法
  3. 缺陷工程: 指导通过氧空位调控改善催化性能
  4. 多相催化: 为理解多相催化中的氧传输过程提供理论框架

参考文献

本文引用了34篇重要文献，涵盖了Fe2_2(MoO4_4)3_3的结构、磁性、电化学和催化性质，以及相关的理论计算方法和实验技术，为研究提供了坚实的理论基础。

总体评价: 这是一篇高质量的理论计算论文，成功解释了实验观察到的拉曼光谱异常现象，为理解铁钼酸盐催化剂的工作机理提供了重要的理论支撑。方法创新性强，计算严谨，结论可靠，具有重要的学术价值和实用意义。