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

Рамановское смягчение, индуцированное кислородными вакансиями в катализаторе 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í (Университет Гёте, Франкфурт)
  • Классификация: cond-mat.mtrl-sci (физика конденсированного состояния - материаловедение)
  • Дата публикации: 15 октября 2025 г.
  • Ссылка на статью: https://arxiv.org/abs/2510.12746

Аннотация

Железомолибдаты (Fe2_2(MoO4_4)3_3) широко используются в промышленных каталитических процессах окислительного дегидрирования метанола в формальдегид, однако детальный механизм катализа остаётся неясным. Недавние исследования методом операндо рамановской спектроскопии и импедансометрии показали, что объёмная область Fe2_2(MoO4_4)3_3 может функционировать как резервуар кислорода, с атомами, мигрирующими на поверхность для участия в каталитической реакции. Этот вывод основан на значительном снижении интенсивности рамановского сигнала, наблюдаемом во время каталитического процесса, что указывает на образование точечных дефектов. Однако микроскопическое происхождение этого снижения интенсивности остаётся неясным. В данном исследовании с помощью расчётов теории функционала плотности (DFT) выявлено происхождение наблюдаемых экспериментально изменений интенсивности рамановского сигнала. Анализ фононов показывает, что кислородные колебательные моды (с небольшим вкладом Mo) появляются вблизи 782 см1^{-1}, что совпадает с частотной областью, где экспериментально измеряется снижение интенсивности рамановского сигнала. Дальнейшее применение метода эффективно замороженных фононов подтверждает, что кислородные колебания являются основной причиной рассчитанного снижения интенсивности рамановского сигнала.

Научный контекст и мотивация

Проблемная постановка

  1. Важность каталитического применения: Fe2_2(MoO4_4)3_3 является важным промышленным катализатором для окислительного дегидрирования метанола в формальдегид, используется десятилетиями
  2. Неясный механизм: Несмотря на широкое применение, детальный механизм каталитического процесса остаётся неясным, особенно роль кислородных атомов
  3. Экспериментальные наблюдения: Операндо рамановские эксперименты Шумахера и соавторов выявили значительное снижение интенсивности рамановского сигнала вблизи 782 см1^{-1} во время каталитического процесса при 500°C

Научная мотивация

  1. Потребность в теоретическом объяснении: Наблюдаемое экспериментально снижение интенсивности рамановского сигнала требует микроскопического теоретического объяснения
  2. Уточнение механизма катализа: Необходимо понимание на атомном уровне роли кислородных вакансий в каталитическом процессе
  3. Связь структура-свойства: Установление количественной связи между кислородными дефектами и изменениями рамановского спектра

Основные вклады

  1. Выявление микроскопического происхождения снижения интенсивности рамановского сигнала: Расчёты DFT доказывают, что кислородные вакансии являются основной причиной наблюдаемого экспериментально снижения интенсивности рамановского сигнала
  2. Установление связи структура-спектр: Идентификация рамановского пика вблизи 782 см1^{-1} как кислородно-доминируемого асимметричного колебания растяжения MoO4_4
  3. Разработка метода эффективно замороженных фононов: Развитие теоретического метода селективного замораживания фононных мод конкретных атомов для анализа изменений интенсивности рамановского сигнала
  4. Объяснение экспериментальных явлений: Теоретические расчёты хорошо согласуются с экспериментальными наблюдениями, объясняя отсутствие смещения пика или уширения во время каталитического процесса

Методология

Методы расчёта

Расчёты теории функционала плотности (DFT):

  • Использование пакета Vienna Ab initio Simulation Package (VASP)
  • Применение метода проекторов, дополненных волнами (PAW)
  • Функционал обменно-корреляционного взаимодействия: PBESol (модифицированное обобщённое градиентное приближение)
  • Поправка Хаббарда U: U = 4 эВ для 3d орбиталей Fe
  • Энергия отсечки плоских волн: 500 эВ
  • 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 эВ)
  • Оптимизация структуры: Критерий сходимости сил < 1 мэВ/Å
  • Гауссово уширение: 0.05 эВ
  • Уширение рамановского спектра: Лоренцево уширение с полушириной 10 см1^{-1}

Сравнительный анализ

Сравнение структурных различий между моноклинной и ортогональной фазами путём анализа симметричных мод с использованием программы AMPLIMODES для определения матрицы трансформации фазового перехода.

Экспериментальные результаты

Анализ фонон-спектра

  1. Разделение частот: Фонон-спектр разделяется на две области
    • Низкочастотная область (0-403 см1^{-1}): вклад атомов Fe, Mo, O
    • Высокочастотная область (780-1001 см1^{-1}): вклад Fe исчезает, преобладают колебания Mo-O
  2. Ключевые колебательные моды:
    • 789 см1^{-1}: Асимметричное колебание растяжения MoO4_4, максимальная интенсивность рамановского сигнала
    • 972 см1^{-1}: Симметричное колебание растяжения MoO4_4, атом Mo практически неподвижен

Результаты расчётов интенсивности рамановского сигнала

Частота (см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: Снижение интенсивности рамановского сигнала на 8.5% при 789.68 см1^{-1}
  • Замораживание атома O5: Снижение интенсивности рамановского сигнала на 7.6% при 786.18 см1^{-1}
  • Замораживание атома Mo: Снижение интенсивности примерно в два раза меньше, чем для кислородных атомов

Релаксация структуры с кислородной вакансией

Релаксация структуры с кислородной вакансией показывает:

  • Треугольные единицы MoO3_3 стремятся восстановить тетраэдрическую геометрию MoO4_4
  • Пятикоординированные единицы FeO5_5 остаются стабильными
  • Значительные изменения параметров решётки, понижение пространственной группы до P1

Связанные работы

Многофункциональные свойства Fe2_2(MoO4_4)3_3

  1. Магнитные свойства: L-образный ферримагнитный переход ниже 12 K, многосегнетоэлектрическая фаза, индуцированная внешним магнитным полем
  2. Электрохимические применения: Использование в качестве электродного материала для Li/Na-ионных батарей с хорошей циклической стабильностью
  3. Каталитические применения: Промышленный катализатор для производства формальдегида из метанола, концентрация Mo на поверхности влияет на селективность

Операндо спектроскопические исследования

Пионерская работа Шумахера и соавторов впервые предложила концепцию объёмного резервуара кислорода, обеспечив экспериментальную основу для данного исследования.

Методы теоретических расчётов

Применение методов DFT+U при изучении электронной структуры и фонон-спектра переходных металлооксидов.

Выводы и обсуждение

Основные выводы

  1. Кислородные вакансии являются ключевым фактором: Кислородные вакансии являются основной причиной снижения интенсивности рамановского сигнала вблизи 782 см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, а также соответствующие методы теоретических расчётов и экспериментальные методики, обеспечивая прочную теоретическую основу для исследования.


Общая оценка: Это высококачественная теоретическая работа, успешно объясняющая аномальные явления в рамановском спектре, наблюдаемые экспериментально, и обеспечивающая важную теоретическую поддержку для понимания механизма действия железомолибдатного катализатора. Методология отличается инновационностью, расчёты строги, выводы надёжны, работа имеет значительную научную ценность и практическое значение.