Научная статья УДК 621.763
DOI: 10.18577/2307-6046-2023-0-10-105-118
ФОРМИРОВАНИЕ УГЛЕРОДНЫХ ПОКРЫТИЙ В ПРОЦЕССЕ ПИРОЛИЗА МЕТАНА МЕТОДОМ ХИМИЧЕСКОГО ГАЗОФАЗНОГО ОСАЖДЕНИЯ
Д.В. Сидоров1, А.А. Грунт1, А.А. Шавнев1
Федеральное государственное унитарное предприятие «Всероссийский научно-исследовательский институт авиационных материалов» Национального исследовательского центра «Курчатовский институт», Москва, Россия; [email protected]
Аннотация. Представлен обзор научно-технической литературы, посвященной пиролизу метана, в результате которого происходит химическое осаждение углерода из газовой фазы. Обобщены результаты исследований процесса в зависимости от параметров осаждения и используемых подложек. Проанализированы составы продуктов пиролиза, механизмы и возможные химические реакции.
Ключевые слова: химическое газофазное осаждение, пиролиз метана, пиролитиче-ский углерод, химическая реакция, реагент
Для цитирования: Сидоров Д.В., Грунин А.А., Шавнев А.А. Формирование углеродных покрытий в процессе пиролиза метана методом химического газофазного осаждения // Труды ВИАМ. 2023. № 10 (128). Ст. 10. URL: http://www.viam-works.ru. DOI: 10.18577/2307-6046-2023-0-10-105-118.
Scientific article
FORMATION OF CARBON COATINGS IN THE PROCESS OF METHANE PYROLYSIS BY THE CHEMICAL VAPOR DEPOSITION
D.V. Sidorov1, A.A. Grunin1, A.A. Schavnev1
federal State Unitary Enterprise «All-Russian Scientific-Research Institute of Aviation Materials» of National Research Center «Kurchatov Institute», Moscow, Russia; [email protected]
Abstract. The article presents a review of the scientific and technical literature in the field of methane pyrolysis in the process of chemical vapor deposition ofpyrolytic carbon, summarizes the results of studies of the process depending on the parameters of deposition and the substrates used, analyzes the compositions of pyrolysis products, mechanisms and pathways of chemical reactions.
Keywords: ^emical vapor deposition, methane pyrolysis, pyrolytic carbon, chemical reaction, reagent
For citation: Sidorov D.V., Grunin A.A., Schavnev A.A. Formation of carbon coatings in the process of methane pyrolysis by the chemical vapor deposition. Trudy VIAM, 2023, no. 10 (128), paper no. 10. Available at: http://www.viam-works.ru. DOI: 10.18577/2307-6046-2023-0-10-105-118.
Введение
Постоянное развитие техники, технологий и совершенствование методов для изготовления новых материалов с уникальными характеристиками требуют глубокого понимания их структурных особенностей и морфологии как при получении, так и при эксплуатации и переработке, а также свойств веществ, из которых их получают, физических и химических процессов [1-4].
Получение пиролитического углерода методом химического газофазного осаждения (ХГО) является одной из важнейших технологий в производстве углеродных материалов, особенно композиционных. Такие материалы широко используются в различных отраслях промышленности, обладают низкой плотностью и стабильными механическими характеристиками при повышенных температурах, а также высокой теплопроводностью и удовлетворительной электропроводностью [5, 6]. Метод ХГО совместим со стандартной полупроводниковой технологией, что облегчает его коммерческое применение в интегральных схемах и солнечных батареях [7].
Химическое газофазное осаждение - достаточно сложный процесс, при котором компоненты реагируют вблизи горячей подложки или на ней. Подложку активируют путем нагревания, воздействия излучения или плазмы, в результате чего образуется твердый осадок. Процесс осаждения состоит либо из гомогенных реакций, происходящих в газовой фазе, либо из гетерогенных - на горячей поверхности или в непосредственной близости от нее, что приводит к образованию покрытий.
Метод ХГО позволяет получать углеродные материалы с необходимой структурой и однородностью путем изменения условий осаждения, таких как: температура, состав подложки и реакционной газовой смеси, скорость потока газовой смеси, продолжительность выдержки и т. д. Обычно в качестве источников углерода используют углеводороды, такие как метан, этилен и пропан. При нагревании до температуры пиролиза образуются разнообразные соединения и радикалы - от углеводородов с небольшим количеством атомов углерода до полициклических ароматических углеводородов, а также углерод, получаемый в результате гетерогенных реакций, при которых любая молекула или радикал в газовой фазе может осаждаться на твердой поверхности [7]. Поскольку при нагревании происходит множество сложных конкурентных реакций и их комбинаций, химический механизм пиролиза углеводородов до сих пор недостаточно исследован. Пиролиз метана по сравнению с пиролизом других углеводородов наиболее изучен благодаря высокой химической стабильности, простому механизму и быстрой интерпретации полученных результатов. Это связано с тем, что в процессе пиролиза метана задействовано меньше продуктов и промежуточных соединений.
В последнее время проведено много исследовательских работ по осаждению углерода методом ХГО, результаты которых являются противоречивыми. Такое разночтение объясняется применением различных технологических подходов, аппаратурными и аналитическими возможностями, функциональными требованиями к углеродному покрытию. Цель данной работы - обобщение информации по осаждению углерода методом ХГО. Основное внимание будет уделено химическим реакциям при пиролизе метана и составу получаемых продуктов.
Работа выполнена с использованием оборудования ЦКП «Климатические испытания» НИЦ «Курчатовский институт» - ВИАМ.
Морфология углеродных покрытий
В работе [8] исследовано влияние концентрации водорода на микроструктуру, морфологию и фазовый состав покрытий из углерода, осажденного в результате пиролиза метана методом ХГО. В качестве реактора использовали стандартный микроволновый источник плазмы мощностью 5 кВт, в качестве подложки - полированные монокристаллические кремниевые пластины. В работе применяли методы сканирующей электронной микроскопии, рентгеновской дифрактометрии и спектроскопии комбинационного рассеяния. Температуру подложки измеряли оптическим датчиком через кварцевый колпак. Продолжительность пиролиза составила 5 ч. Концентрацию водорода варьировали от 5 до 20 %, а содержание метана поддерживали на постоянном
уровне. На рис. 1 представлены морфология и поперечное сечение покрытий из пиро-литического углерода, осажденного на кремниевую подложку. При увеличении концентрации водорода морфология покрытий изменялась - от зернистой структуры до игольчатых скоплений. Размер зерна увеличивался с повышением концентрации водорода. Толщина покрытий составляла 1,75; 1,8; 1,65 и 2,9 мм при концентрации водорода 5, 10, 15 и 20 % соответственно. Признаки столбчатого роста кристаллов отсутствовали. Скорость роста кристаллов составила 0,35; 0,36; 0,32 и 0,56 мм/ч при концентрации водорода 5, 10, 15 и 20 % соответственно. При концентрации водорода от 5 до 15 % скорость роста оставалась практически неизменной, от 15 до 20 % - быстро увеличивалась, чему соответствовали большой размер зерна и высокая кристалличность.
Морфология покрытий
Поперечное сечение покрытий
Рис. 1. Морфология и поперечное сечение покрытий из осажденного углерода при концентрации водорода 5 (а), 10 (б), 15 (в) и 20 % (г) [8]
В работах [9-11] исследована температурная зависимость процесса осаждения углерода при пиролизе метана методом ХГО. В работе [9] использовали вертикальный кварцевый проточный реактор, разогреваемый до температуры 600 °С и охлаждаемый газообразным азотом. В качестве катализатора применяли цеолит марки ZSM-5, модифицированный никелем (№/2БМ-5). Метан подавали в течение 15 мин. Исследования проводили методом газовой хроматографии, сканирующей электронной микроскопии, рентгеновской дифрактометрии и спектроскопии комбинационного рассеяния. На рис. 2 представлены углеродные покрытия, полученные при температурах 400 и 600 °С. При температуре 400 °С в основном формировались углеродные структуры типа «рыбья кость», диаметр нитей не зависел от концентрации никеля и составил 25 нм. При температурах >600 °С получены полые структуры. Внешний диаметр нитей увеличивался с 29 до 60 нм при повышении содержания никеля с 2 до 8 %.
Однако диаметр нитей при использовании №/2БМ-5 в концентрации 2 и 11,4 % не изменялся и составил 29 нм. Таким образом, в работе [9] показана явная зависимость морфологии углеродных покрытий, полученных при каталитическом разложении метана с использованием Ni/ZSM-5, от температуры.
В работе [10] в качестве реактора использовали кварцевую трубку длиной 260 мм и внутренним диаметром 34 мм. Подложку нагревали с помощью прямого излучения вольфрамовой нити. Температура нити, которую измеряли оптическим пирометром, составила 1800-2300 °С. Газ-реагент, состоящий из метана и водорода, подавали в реактор через массовый расходомер. В качестве подложек использовали кремниевые
пластины. Плотность зародышеобразования структуры увеличивалась, достигала максимума при температуре 2100 °С, а затем уменьшалась (рис. 3). Морфология покрытия подложки также зависела от продолжительности обработки и температуры нити.
Осаждение при температуре 400 °С
Осаждение при температуре 600 °С
Рис. 2. Морфология покрытий из углерода, осажденного при температурах 400 и 600 °С в течение 15 мин с использованием в качестве катализатора Ni/ZSM-5 в концентрации 2 (а), 4 (б), 8 (в) и 11,4 % (г) [9]
ю о
о
Л
/\ ' \
о/ /
/ О
ъ
Ць^-
1800
2000 Температура, °С
2200
Рис. 3. Зависимость плотности зародышеобразования структуры от температуры нити [10]
В работе [11] изучено влияние температуры, скорости потока метана и парциального давления на морфологию покрытий из осажденного углерода, которую анализировали с помощью спектроскопии комбинационного рассеяния, просвечивающей и сканирующей электронной микроскопии. В качестве реактора ХГО использовали печь с разделенными трубками с зоной нагрева 150 мм и кварцевой трубкой с внешним диаметром 25 мм, в качестве подложки - медную фольгу толщиной 25 мкм. На рис. 4 представлены покрытия из углерода, осажденного при различных значениях температуры, скорости потока метана и парциального давления. Установлено, что плотность зародышеобразования структуры может уменьшаться при увеличении температуры
4
3
2
1
0
(рис. 4, а, б) или уменьшении скорости потока метана (рис. 4, б, в) и парциального давления (рис. 4, б, г). Высокая температура и низкие скорость метана и парциальное давление приводят к низкой плотности зародышеобразования структуры и, следовательно, к большему размеру домена. Рост структуры прекращается при полном покрытии поверхности (рис. 4, г).
Рис. 4. Морфология углеродных покрытий, полученных в различных условиях (температура (°С)/скорость потока метана (см3/мин)/парциальное давление (мм рт. ст.)): 985/35/460 (а); 1035/35/460 (б); 1035/7/460 (в); 1035/7/160 (г) [11]
В работе [12] исследовано влияние состава газовой смеси на морфологию покрытия из осажденного пиролитического углерода. В исследовании, как и в работе [10], использовали осаждение методом «горячей нити». Для нагрева применяли три линейные вольфрамовые проволоки длиной 66 мм и диаметром 0,25 мм. Газ представлял собой смесь аргона, водорода и метана. В качестве подложки использовали нитрид кремния Si3N4. Расстояние между нитями и подложкой составляло 5 мм, давление газа: ~5 кПа.
Рис. 5. Морфология углеродного покрытия, полученного из газовой смеси в различных условиях (отношение потоков газов Аг и Н2/газов СН4 и Н2/скорость потока газовой смеси (мл/мин)): 2,5/0,04/50 (а); 2,5/0,04/100 (б); 2,5/0,14/100 (в); 1,2/0,04/50 (г); 1,2/0,04/100 (д); 7,2/0,04/100 (е) [11]
Сплошные структуры получены при следующих условиях осаждения: температура нити 2200 °С; температура подложки 650 °С; продолжительность осаждения 1 ч; отношение потоков газов: 1,2 - для Аг и Н2, 0,04 - для СН4 и Н2; скорость потоков газа 50 и 100 мл/мин соответственно (рис. 5). При увеличении соотношения газов СН4/Н2 количество графитовых фаз также увеличивалось. Зародышеобразование происходило при отношении потоков газов Аг и Н2, равном 2,5; при отношении 7,2 процесс осаждения нестабилен.
Химические реакции при пиролизе метана
При пиролизе метана протекает множество химических реакций с образованием различных радикалов и молекулярных структур, идентификация наиболее важных из которых имеет решающее значение для формирования заданной структуры пиролити-ческого углерода.
Пиролитическому разложению метана посвящены многочисленные экспериментальные и теоретические исследования. В ранних работах в качестве реактора использовали проточную камеру. Камера в работе [13] представляла собой кварцевую трубку объемом 478 см3, длиной 20,4 см и диаметром 6,0 см, в работе [14] - фарфоровую трубку с внутренним диаметром 5 мм. В качестве нагревательного элемента использовали автоклавную высокотемпературную цилиндрическую печь. В работе [13] температура пиролиза составила 995-1100 К, в работе [14]: 1323-1523 К. Продукты реакции анализировали методами газовой хроматографии и масс-спектрометрии. Разложение молекулы метана описывается следующими химическими уравнениями:
СН4 ^ СН3 + Н (7)
Н + СН4 ^ СН3 + Н2 (2)
2СН3 ^ С2Н6 (3)
СН3 + С2Н6 ^ СН4 + С2Н5 (4)
С2Н5 ^ С2Н4 + Н (5)
Рассмотренный механизм основан на следующих предпосылках: только радикал СН3 является значимым, а углеводородные продукты быстро вступают в цепную реакцию с метильными радикалами. Атомарный водород всегда участвует в реакции 2 с образованием молекулярного водорода. Результаты исследования устанавливают общий характер изменения скорости разложения метана в проточной камере: реакция ускоряется, затем замедляется, постоянная скорость появляется только на мгновение, как перегиб. При низких температурах наблюдается только ускорение, при высоких - короткий индукционный период. Общая энергия активации находится в пределах от 422 до 430 кДж/моль.
В работе [15] предложена следующая модель пиролиза метана:
СН4 ^ С2Н6 ^ С2Н4 ^ С2Н2 ^ С + Н2 (6)
В качестве реактора использовали импульсную трубку из нержавеющей стали длиной 3,3 м. Продукты реакции анализировали с помощью хроматографии, а также методом масс-спектрометрии. Пиролиз осуществляли в диапазоне температур 12002400 °С. Общие кинетические уравнения позволили сделать выводы относительно выхода стабильных промежуточных продуктов. При пиролизе метана концентрация этана всегда чрезвычайно мала. При высоких температурах среди продуктов реакции этан отсутствует, поскольку константа скорости его разложения на 2-3 порядка больше, чем константа скорости разложения метана. Выводы данной работы согласуются с результатами работ [16, 17]. Аналогичные результаты получены для этилена и ацетилена. Ацетилен образуется при температурах >1500 К. Температура >2000 К способствует максимальной конверсии метана в ацетилен, а в области температур от 1000 до 2500 К в продуктах реакции присутствуют водород и углерод.
В работе [18] приведены результаты моделирования пиролиза метана при воздействии плазмы. Для исследования использовали стандартный микроволновый источник плазмы, работающий при атмосферном давлении с частотой 2,45 ГГц. Разряд пропускали через кварцевую трубку, помещенную в металлическую трубу. Внутренний диаметр разрядной трубки составил 26 мм, диаметр плазменного столба: ~20 мм,
скорость подачи метана: 50 дм /мин, мощность плазмы: 2, 3 и 4 кВт. Для расчетов использовали модели 0D и Ш, при реализации которых применяли бомбовый калориметр и реактор вытеснения соответственно. Температура пиролиза метана составила 15002000 К, температура плазмы: 4000-5700 К. Механизм пиролиза описан теоретическим методом, который не учитывает второстепенных реакций. В результате авторы отметили 48 наиболее важных химических реакций, где М - катализатор, адс - адсорбированные молекулы:
Н2 + СН2 — СНз + Н (7) СН2 + СН3 — С2Н4 + Н (31)
СН4 + С — СН3 + СН (8) СН2 + С2Н3 — С2Н2 + СН3 (32)
СН4 + Н — СН3 + Н2 (9) СН2(адс) + М — СН2 + М (33)
СН4 + СН — С2Н4 + Н (10) 2СН3 + М — С2Н + М (34)
СН4 + СН2 — 2СН3 (11) СН3 + М — СН2 + Н + М (35)
СН4 + СН2(адс) — 2СН3 (12) С2Н + С2Н3 — 2С2Н2 (36)
СН4 + С2Н — СН3 + С2Н2 (13) С4Н3 + М — С4Н2 + Н + М (37)
2С2Н2 — С4Н3 + Н (14) С2Н2 + Н + М — С2Н3 + М (38)
Н2 + С2Н — С2Н2 + Н (15) С2Н2 + СН — С2Н + СН2 (39)
С + СН3 — С2Н2 + Н (16) С2Н2 + СН2 — С3Н4 (40)
С + СН2 — С2Н + Н (17) С2Н2 + СН2 — С3Н3 + Н (41)
2Н + М — Н2 + М (18) С2Н2 + С2Н — С4Н2 + Н (42)
2Н + Н2 — 2Н2 (19) С2Н2 + М — С2Н + Н + М (43)
Н + СН — С + Н2 (20) С2Н4 + Н — С2Н3 + Н2 (44)
Н + СН2 — СН + Н2 (21) С2Н4 + Н + М — С2Н5 + М (45)
Н + СН2(адс) — СН2 + Н (22) С2Н4 + СН — С3Н4 + Н (46)
Н + СН3 — СН4 + М (23) С2Н4 + СН2 —* С3Н4 + Н2 (47)
Н + С2Н3 — С2Н2 + Н2 (24) С2Н4 + СН3 — С2Н3 + СН4 (48)
2СН3 — С2Н5 + Н (25) С2Н4 + М — С2Н2 + Н2 + М (49)
СН + СН2 — С2Н2 + Н (26) С2Н4 + М — С2Н3 + Н + М (50)
СН + СН3 — С2Н3 + Н (27) С2Н6 + Н — С2Н5 + Н2 (51)
СН + С2Н3 — СН2 + С2Н2 (28) С2Н6 + СН — С2Н4 + СН3 (52)
2СН2 — С2Н2 + Н2 (29) С2Н6 + СН2 — С2Н5 + СН3 (53)
2СН2 — С2Н2 + 2Н (30) С2Н6 + СН3 — С2Н5 + СН4 (54)
Образование водородсодержащих радикалов начинается с разложения СН4 и заканчивается при появлении С3Н4 и С4Н3. Несмотря на то что конечная концентрация этих продуктов очень низкая, они имеют важное значение для образования углерода. Большинство атомов водорода в химических реакциях переходит в молекулярный водород. Другие продукты пиролиза являются промежуточными и короткоживущими. Показана аналогичная последовательность образования радикалов для атомов углерода. Радикал С2Н является основным, а радикал С не имеет прямого отношения к образованию углерода. Соединение С2 обнаружено спектроскопическим методом в продуктах пиролиза. Однако образование и вклад С2 в предложенной модели были настолько незначительными, что данное соединение было исключено из механизма пиролиза. Конверсия метана и образование водорода происходят в основном в реакциях 7, 9, 15, 18, 20-22, 24, 44, 49, 51. В реакции 9 атомы водорода атакуют молекулу метана, в результате чего выделяется 74,6 % Н2. В реакциях 44, 49 и 51 образуется 6,8; 7,1 и 10,8 % Н2 соответственно. В реакции 9 разлагается 71,5 % метана. В реакциях 48 и 54 большая часть метана (~31 %) образуется вновь. Ацетилен участвует во многих реакциях, но в основном образуется при дегидрировании С2Н4 в прямой реакции 49 и опосредованно - через радикал С2Н3 в реакции 38. Механизм пиролиза метана, предложенный в работе [18], объясняет экспериментальные результаты и позволяет выявить
основные реакции, ответственные за конверсию метана и образование соединений, обнаруженных в эксперименте. С помощью модели удалось спрогнозировать основные продукты пиролиза метана: Н2, С2Н2 и С2Н4.
В работе [19] проведено моделирование процесса пиролиза метана при химическом осаждении из газовой фазы. Модель включает 909 реакций. Достигнута высокая сходимость результатов моделирования и экспериментальных данных. Для осаждения использовали вертикальный проточный реактор диаметром 22 мм и длиной 20 мм, в качестве газа-разбавителя - аргон. Температура пиролиза метана составила 1373 К, продолжительность выдержки в реакторе при атмосферном давлении: 0,045; 0,495 и 0,972 с. Для анализа газообразных продуктов применяли газовый хроматограф, оснащенный газоразделительной системой.
а)
б)
С}Нйс=>СН,-СН-СН, Е=> нс=сн-сн,
" п
сн,=с=сн2 II
Расход соединения: ф Аллш
эНС=С-СН}
Пропин
I—~Г_-
сн,4> сн,
{} С.н, С7нГ
СН,=С =СН,-<=1НС=С-СН)^-
Аллен
Пропин
50-100% ¡-ос-СМ.С^
П^У 30-50 %
10-30 %
5-10 % 1 -5 %
9
(
С6н,
Бензол
$ *
0<=0
Бензол
в)
сн,
сн,
с н,<Ьс н,<Ьс\н, нс^сн-сн^нс^с-сн^нс^с-сн,
Пропин
^ СН,=СН-СН,ОСН,=С=СН, Г 1
Аллен
Со^Сос,н>
С„Н7
I
Бензол ▲
Инден
Рис. 6. Схемы пиролиза метана при продолжительности выдержки в реакторе 0,045 (а), 0,495 (б) и 0,972 с (в) [19]
При любой продолжительности выдержки в реакторе бензол образуется из соединений, содержащих три атома углерода (рис. 6), что совпадает с результатами пиро-
лиза этилена, пропилена и пропана, описанными в работе [20]. В данной работе использовали реактор с вертикальным потоком, работающий при температуре 900 °С, продолжительность выдержки не превышала 1,6 с. На рис. 6, а представлены начальная стадия пиролиза метана и образование бензола. Метан участвует в реакциях дегидрирования с образованием метила и метилирования с образованием этильного радикала. Этильный радикал подвергается дегидрированию, этилен взаимодействует с метиль-ными радикалами и углеводородами. Второй путь разложения метана - дегидрирование с образованием аллена, а затем пропина. Аллен и пропин являются основными источниками радикала С3Н3, из которого в конечном итоге образуется бензол. При увеличении продолжительности выдержки в реакторе до 0,495 с образуются крупные углеводороды, такие как циклопентадиен и инден (рис. 6, б). Атака ацетилена радикалами С3Н5 или С7Н7 является основным путем образования циклопентадиена и индена. Такие реакции не происходят при пиролизе ацетилена и пропилена. Источником бензола также выступают аллены (5,6 %) и радикалы С3Н3 (2,2 %). Дегидрирование и рекомбинация С3Н3 ускоряются, когда продолжительность выдержки в реакторе увеличивается с 0,045 до 0,495 с. Это свидетельствует о более активном протекании элементарных реакций.
Ацетилен является ключевым промежуточным продуктом пиролиза метана при продолжительности выдержки в реакторе 0,972 с. Источниками образования ацетилена являются этилен, винил и винилацетилен, при этом из 14,9 % ацетилена образуется радикал С3Н5. При атаках ацетилена радикалами С5Н5, С3Н5 и С7Н7 образуются С7Н7, циклопентадиен и инден. Показаны также вероятные пути образования бензола при взаимодействии С3Н5 и С3Н3, а также винилацетилена и С2Н2.
В работе [21] для изучения механизма пиролиза метана использовали метод расчета электронной структуры по теории функционала плотности (DFT). Оптимизация геометрии и анализ частоты колебаний всех химических соединений и переходных состояний выполнены с помощью функционала плотности (B3Lyp) с базисным набором б-311+О (ё, р). Константы скорости рассчитаны с помощью программы KiSThelP в соответствии с теорией переходного состояния, для получения поправочных коэффициентов туннелирования использовали метод Вигнера. Константы скорости подобраны для модифицированного выражения Аррениуса в интервале температур 800-2000 К. В результате термодинамических и кинетических расчетов выявлен наиболее вероятный путь пиролиза метана (рис. 7). Определены профили реакций и константы скоростей каждого пути перехода соединений С1 к Сб, предложены наиболее вероятные пути протекания химических реакций. Выдвинуто предположение, что переход от соединений С4 к Сб сопровождается образованием бензола.
сн
-н
II, I!
> СН
,
* с,н<
-н
,
г см
-н
,
С Н -
4 2 _Н
-с,н -н
-с,
+с,н
-н,
+с,н.
1,1
СЛл-
-н
сл.
-н
,
-н,
> с,а——с,
,
-н
С
+с2н3
сл7
+с2н,
I_+ С2НГН-Н, I"
с6н,
+С2Н +н
+с,н, +с,н, -н
+ С,Н3 -н,
Рис. 7. Схема превращения бензола из метана [21]
В работе [22] описана подробная газофазная кинетическая модель пиролиза метана, состоящая из 37 компонентов и 318 реакций. Для разработки подробного механизма использовали программное обеспечение RMG. Температура пиролиза составила 1373 К, давление: 75 кПа. Условия работы реактора, использованные в экспериментах, воспроизведены при моделировании. Полученные концентрации основных продуктов соответствуют экспериментальным результатам. Кинетическое выражение и термодинамические характеристики соединений получены из механизма, генерируемого ЯМО. Анализ чувствительности и диаграммы потока реакционного пути показал, что пиролиз метана протекает без индукционного периода. Соединения С2 являются первичными промежуточными продуктами для образования ацетилена, этилена и бензола. Детальная кинетическая модель систематически сведена к механизму, включающему 29 реакций:
С2н4 + M — H + C2H3 + M (55) C2H6 + Н — C2H5 + Н2 (70)
H + С2Н4 — C2H3 + H2 (56) CH3 + C2H6 — C2H5 + CH4 (71)
2H (+ M) — H2 (+ M) (57) C3H6 + H — CH3 + C2H4 (72)
CH4 + M — CH3 + H + M (55) C3H6 — CH3 + C2H3 (73)
CH3 + C2H4 — CH4 + C2H3 (59) C4H6 + C2H3 — H + H2 + C6H6 (74)
CH3 + H2 — CH4 + H (60) C4H6 + H —> C2H3 + C2H4 (75)
2CH3 + M — C2H6 + M (61) C4H6 — 2C2H3 (76)
CH3 + CH4 — C2H6 + H (62) C3H6 + C2H3 (+M) — CH3 + C4H6 (+M) (77)
C2H5 + M — H+ C2H4 + М (63) C2H4 + M — C2H2 + H2 + M (78)
C2H6 — C2H5 + H (64) H + C2H3 — C2H2 + H2 (79)
2CH3 — C2H5 + Н (65) C3H6 — C2H2 + CH4 (80)
CH3 + CH4 — C2H5 + H2 (66) C2H3 + M — C2H2 + H + M (81)
C2H5 + H — C2H4 + H2 (67) C4H6 + H — C2H2 + C2H3 + H2 (82)
2C2H4 — C2H5 + C2H3 (68) C4H6 — C2H2 + C2H4 (83)
C2H6 + C2H3 — C2H5 + C2H4 (69)
Теоретические результаты, полученные в работе [22], подтверждены экспериментально в работах [23, 24] с использованием реактора ХГО. Реактор представлял собой камеру для осаждения из алюмокерамической трубки диаметром 22 мм и длиной 40 мм, внутри которой установлены трубки из кордиерита. Подложка состоит из 10 участков канальной структуры длиной ~4 мм. Поверхность покрытия из осажденного углерода исследовали с помощью сканирующей электронной микроскопии. Состав газовой фазы анализировали на газовом хроматографе, оснащенном специальной системой вакуумного дозирования. Охлажденные газы разделяли и анализировали с помощью детектора теплопроводности. Продукты конденсации определяли по углеродному балансу. Сравнили данную модель с предложенными в работах [19, 25]. Упрощенная модель (рис. 8), представленная в работе [25], отражает пути пиролиза метана в условиях ХГО при давлении 5,3 103 Па и температуре 1273-1373 К.
Рис. 8. Упрощенная модель пиролиза метана по основному и второстепенному
пути [25]
Модель, предложенная в работе [19], описывает подробный механизм пиролиза метана без регулируемых параметров. Сравнительный анализ работ показывает, что при использовании модели из работы [25] содержание почти всех соединений было занижено. Модель, представленная в работе [19], обеспечивает совпадение рассчитанных данных с экспериментальными, при этом занижено содержание таких соединений, как C2H2 и C2H4. Отклонения в содержании метана возрастали с увеличением продолжительности выдержки в реакторе. Модель расчета в работе [22] наиболее полно фиксирует все значимые продукты реакций и позволяет определить их содержание. Как показано в работе [22], при пиролизе метана происходит только реакция 80, при которой C3H6 разлагается на C2H2 и CH4 и тем самым ингибируется превращение CH4. Этот вывод согласуется с научно-техническими литературными данными о том, что определенное количество метана образуется при пиролизе высших углеводородов [20, 26, 27]. К основным промежуточным соединениям при пиролизе метана относятся C2H2, C2H4 и C6H6. Образование C2H2 с наибольшей вероятностью происходит по реакции 81. Концентрация промежуточного соединения C2H5 определяет относительное количество C2H2 и C2H4 в газовой смеси. Отмечены девять наиболее вероятных реакций синтеза C6H6, в большинстве из которых участвуют соединения C2. Таким образом, соединения C2 являются необходимыми промежуточными продуктами для синтеза C6H6. Альтернативный путь образования C6H6 включает разложение метана с участием соединений C3 и C4. Образование C2H2, C2H4 и C6H6 в основном происходит с участием соединения C2, а промежуточные соединения, такие как C3 и C4, влияют на образование C6H6 в меньшей степени. Отмечено, что соединения C3 участвуют в образовании радикала CH3 и не вносят вклад в образование C6H6.
Заключения
Показано, что метод ХГО позволяет получить углеродные покрытия с заданной структурой и однородностью. Синтез углеродных покрытий методом ХГО является достаточно сложным процессом, а их качество зависит от различных условий. Основное внимание уделено механизму химических реакций, составу получаемых продуктов, их влиянию на ход процесса, роли технологических факторов (давление, температура), конструкции камер для осаждения и используемым подложкам. Для получения углеродных покрытий применяют в основном предельные углеводороды начального ряда, при нагревании которых до температуры пиролиза образуются разнообразные соединения и радикалы. Среди наиболее вероятных механизмов пиролиза метана следует отметить образование метильного радикала на начальной стадии разложения и химические реакции, связанные с последующим образованием ацетилена, этилена и бензола.
Список источников
1. Каблов Е.Н., Ечин А.Б., Бондаренко Ю.А. История развития технологии направленной кристаллизации и оборудования для литья лопаток газотурбинных двигателей // Труды ВИАМ. 2020. № 3 (87). Ст. 01. URL: http://www.viam-works.ru (дата обращения: 11.05.2023). DOI: 10.18577/2307-6046-2020-0-3-3-12.
2. Каблов Е.Н., Семенова С.Н., Сулейманов Р.Р., Чайкун А.М. Перспективы применения эти-ленпропилендиенового каучука в составе морозостойкой резины // Труды ВИАМ. 2019. № 12 (84). Ст. 04. URL: http://www.viam-works.ru (дата обращения: 11.05.2023). DOI: 10.18577/2307-6046-2019-0-12-29-36.
3. Каблов Е.Н., Валуева М.И., Зеленина И.В., Хмельницкий В.В., Алексашин В.М. Углепластики на основе бензоксазиновых олигомеров - перспективные материалы // Труды ВИАМ. 2020. № 1 (85). Ст. 07. URL: http://www.viam-works.ru (дата обращения: 11.05.2023). DOI: 10.18577/2307-6046-2020-0-1-68-77.
4. Сидоров Д.В., Шавнев А.А., Мелентьев А.А. Формирование карбидокремниевых покрытий методом химического газофазного осаждения (обзор). Часть 1 // Труды ВИАМ. 2021. № 6 (100). Ст. 10. URL: http://www.viam-works.ru (дата обращения: 11.05.2023). DOI: 10.18577/2307-6046-2021-0-6-100-111.
5. Фиалков А.С. Углерод, межслоевые соединения и композиты на его основе. М.: Аспект Пресс, 1997. 718 с.
6. Николаев А.И. Получение углерод-углеродных композиционных материалов на основе углеродных нановолокон газофазным способом // Тонкие химические технологии. 2015. Т. 10. С.61-66.
7. Jia K., Zhang J., Zhu Y. et al. Toward the commercialization of chemical vapor deposition gra-phene films // Applied Physics Reviews. 2021. Vol. 8. P. 400-401.
8. Liu J., Hei L., Chen G. et al. Growth of ultrananocrystalline diamond films in an Ar-rich CH4/H2/Ar atmosphere with varying H2 concentrations // New Carbon Materials. 2013. Vol. 28. P.134-138.
9. Majewska J., Michalkiewicz B. Carbon nanomaterials produced by the catalytic decomposition of methane over Ni/ZSM-5: Significance of Ni content and temperature // New Carbon Materials. 2014. Vol. 29. P. 102-108.
10. Park S.-S., Lee J.-Y. Nucleation behavior of diamond particles on silicon substrates in a hot-filament chemical vapour deposition // Journal of Materials Science. 1993. Vol. 28. P. 1799-1804.
11. Li X., Magnuson C.W., Venugopal A. et al. Graphene Films with Large Domain Size by a Two-Step Chemical Vapor Deposition Process // Nano Letters. 2010. Vol. 10. P. 4328-4334.
12. Amaral М., Oliveira F.J., Belmonte M. et al. Hot-filament chemical vapour deposition of nanodi-amond on silicon nitride substrates // Diamond and Related Materials. 2004. Vol. 13. P. 643-647.
13. Chen C.-J., Back M.H., Back R.A. The thermal decomposition of methane. II. Secondary reactions, autocatalysis and carbon formation; non-Arrhenius behavior in the reaction of CH3 with ethane // Canadian Journal of Chemistry. 1976. Vol. 54. P. 3175-3184.
14. Palmer H.B., Lahaye J., Hou K.C. Kinetics and mechanism of the thermal decomposition of methane in a flow system // The Journal of Physical Chemistry. 1968. Vol. 72. P. 348-353.
15. Knorre V.G., Kozlov G.I. Single-pulse shock tube studies on the kinetics of the thermal decomposition of methane // Combustion and Flame. 1962. Vol. 6. P. 253-263.
16. Kevorkian B.V., Heath C.E., Boudart M. The decomposition of methane in shock waves // The journal of physical chemistry. 1960. Vol. 64. P. 964.
17. Hartig R., Troe J., Wagner H.G.G. Thermal decomposition of methane behind reflected shock waves // Symposium (International) on Combustion. 1971. Vol. 13. P. 147-154.
18. Dors M., Nowakowska H., Jasin'sk J., Mizeraczyk J. Chemical kinetics of methane pyrolysis in microwave plasma at atmospheric pressure // Plasma Chemistry and Plasma Processing. 2013. Vol. 34. P. 313-326.
19. Hu C., Li H., Zhang S., Li W. A molecular-level analysis of gas-phase reactions in chemical vapor deposition of carbon from methane using a detailed kinetic model // Journal of Materials Science. 2016. Vol. 51. P. 3897-3906.
20. Norinaga K., Deutschmann O. Detailed Kinetic Modeling of Gas-Phase Reactions in the Chemical Vapor Deposition of Carbon from Light Hydrocarbons // Industrial & Engineering Chemistry Research. 2007. Vol. 46. P. 3547-3557.
21. Li K., Li H., Yan N. et al. Conversion of methane to benzene in CVI by density functional theory study // Scientific Reports. 2019. Vol. 6. P. 19496.
22. Shinde V.M., Pradeep P. Detailed Gas-Phase Kinetics and Reduced Reaction Mechanism for Methane Pyrolysis involved in CVD/CVI Processes // Journal of Analytical and Applied Pyrolysis. 2021. Vol. 154. P. 104998.
23. Becker A., Huttinger K.J. Chemistry and kinetics of chemical vapour deposition of pyrocarbon: IV. Pyrocarbon deposition from methane in the low-temperature regime // Carbon. 1998. Vol. 36. P.213-224.
24. Bruggert M., Hu Z., Huttinger K.J. Chemistry and kinetics of chemical vapour deposition of pyrocarbon: VI. Influence of temperature using methane as a carbon source // Carbon. 1999. Vol. 37. P.2021-2030.
25. Birakayala N., Evans A. A reduced reaction model for carbon CVD/CVI processes // Carbon. 2002. Vol. 40. P. 675-683.
26. Li H., Li A., Bai R., Li K. Numerical simulation of chemical vapour infiltration of propylene into C/C composites with reduced multi-step kinetic models // Carbon. 2005. Vol. 43. P. 2937-2950.
27. Xu W., Zhang Z.-W., Bai R.-C. et al. Kinetic model of gas-phase reactions in the chemical vapor deposition of propane // Carbon. 2014. Vol. 71. P. 345.
References
1. Kablov E.N., Echin A.B., Bondarenko Yu.A. History of development of directional crystallization technology and equipment for casting blades of gas turbine engines. Trudy VIAM, 2020, no. 3 (87), paper no. 01. Available at: http://www.viam-works.ru (accessed: May 11, 2023). DOI: 10.18577/2307-6046-2020-0-3-3-12.
2. Kablov E.N., Semenova S.N., Suleymanov R.R., Chaykun A.M. Prospects for the use of ethylene-propylenediene rubber as part of cold resistant rubber. Trudy VIAM, 2019, no. 12 (84), paper no. 4. Available at: http://www.viam-works.ru (accessed: accessed: May 11, 2023). DOI: 10.18577/2307-6046-2019-0-12-29-36.
3. Kablov E.N., Valueva M.I., I.V. Zelenina, Khmelnitskiy V.V., Aleksashin V.M. Carbon plastics based on benzoxazine oligomers - perspective materials. Trudy VIAM, 2020, no. 1, paper no. 07. Available at: http://www.viam-works.ru (accessed: May 11, 2023). DOI: 10.18577/2307-60462020-0-1-68-77.
4. Sidorov D.V., Schavnev A.A., Melentev A.A. Formation of silicon carbide coatings by chemical vapor deposition (review). Part 1. Trudy VIAM, 2021, no. 6 (100), paper no. 10. Available at: http://www.viam-works.ru (accessed: May 11, 2023). DOI: 10.18577/2307-6046-2021-0-6-100-111.
5. Fialkov A.S. Carbon, interlayer compounds and composites based on it. Moscow: Aspect Press, 1997, 718 p.
6. Nikolaev A.I. Preparation of carbon-carbon composite materials based on carbon nanofibers by gas-phase method. Tonkie khimicheskie tekhnologii, 2015, vol. 10, pp. 61-66.
7. Jia K., Zhang J., Zhu Y. et al. Toward the commercialization of chemical vapor deposition gra-phene films. Applied Physics Reviews, 2021, vol. 8, pp. 400-401.
8. Liu J., Hei L., Chen G. et al. Growth of ultrananocrystalline diamond films in an Ar-rich CH4/H2/Ar atmosphere with varying H2 concentrations. New Carbon Materials, 2013, vol. 28, pp.134-138.
9. Majewska J., Michalkiewicz B. Carbon nanomaterials produced by the catalytic decomposition of methane over Ni/ZSM-5: Significance of Ni content and temperature. New Carbon Materials, 2014, vol. 29, pp. 102-108.
10. Park S.-S., Lee J.-Y. Nucleation behavior of diamond particles on silicon substrates in a hot-filament chemical vapour deposition. Journal of Materials Science, 1993, vol. 28, pp. 1799-1804.
11. Li X., Magnuson C.W., Venugopal A. et al. Graphene Films with Large Domain Size by a Two-Step Chemical Vapor Deposition Process. Nano Letters, 2010, vol. 10, pp. 4328-4334.
12. Amaral M., Oliveira F.J., Belmonte M. et al. Hot-filament chemical vapour deposition of nanodi-amond on silicon nitride substrates. Diamond and Related Materials, 2004, vol. 13, pp. 643-647.
13. Chen C.-J., Back M.H., Back R.A. The thermal decomposition of methane. II. Secondary reactions, autocatalysis and carbon formation; non-Arrhenius behavior in the reaction of CH3 with ethane. Canadian Journal of Chemistry, 1976, vol. 54, pp. 3175-3184.
14. Palmer H.B., Lahaye J., Hou K.C. Kinetics and mechanism of the thermal decomposition of methane in a flow system. The Journal of Physical Chemistry, 1968, vol. 72, pp. 348-353.
15. Knorre V.G., Kozlov G.I. Single-pulse shock tube studies on the kinetics of the thermal decomposition of methane. Combustion and Flame, 1962, vol. 6, pp. 253-263.
16. Kevorkian B.V., Heath C.E., Boudart M. The decomposition of methane in shock waves. The Journal of physical chemistry, 1960, vol. 64, p. 964.
17. Hartig R., Troe J., Wagner H.G.G. Thermal decomposition of methane behind reflected shock waves. Symposium (International) on Combustion, 1971, vol. 13, pp. 147-154.
18. Dors M., Nowakowska H., Jasin'sk J., Mizeraczyk J. Chemical kinetics of methane pyrolysis in microwave plasma at atmospheric pressure. Plasma Chemistry and Plasma Processing, 2013, vol. 34, pp. 313-326.
19. Hu C., Li H., Zhang S., Li W. A molecular-level analysis of gas-phase reactions in chemical vapor deposition of carbon from methane using a detailed kinetic model. Journal of Materials Science, 2016, vol. 51, pp. 3897-3906.
20. Norinaga K., Deutschmann O. Detailed Kinetic Modeling of Gas-Phase Reactions in the Chemical Vapor Deposition of Carbon from Light Hydrocarbons. Industrial & Engineering Chemistry Research, 2007, vol. 46, pp. 3547-3557.
21. Li K., Li H., Yan N. et al. Conversion of methane to benzene in CVI by density functional theory study. Scientific Reports, 2019, vol. 6, p. 19496.
22. Shinde V.M., Pradeep P. Detailed Gas-Phase Kinetics and Reduced Reaction Mechanism for Methane Pyrolysis involved in CVD/CVI Processes. Journal of Analytical and Applied Pyrolysis, 2021, vol. 154, p. 104998.
23. Becker A., Huttinger K.J. Chemistry and kinetics of chemical vapour deposition of pyrocarbon: IV. Pyrocarbon deposition from methane in the low-temperature regime. Carbon, 1998, vol. 36, pp. 213-224.
24. Bruggert M., Hu Z., Huttinger K.J. Chemistry and kinetics of chemical vapour deposition of pyrocarbon: VI. Influence of temperature using methane as a carbon source. Carbon, 1999, vol. 37, pp. 2021-2030.
25. Birakayala N., Evans A. A reduced reaction model for carbon CVD/CVI processes. Carbon, 2002, vol. 40, pp. 675-683.
26. Li H., Li A., Bai R., Li K. Numerical simulation of chemical vapour infiltration of propylene into C/C composites with reduced multi-step kinetic models. Carbon, 2005, vol. 43, pp. 2937-2950.
27. Xu W., Zhang Z.-W., Bai R.-C. et al. Kinetic model of gas-phase reactions in the chemical vapor deposition of propane. Carbon, 2014, vol. 71, p. 345.
Информация об авторах
Сидоров Денис Викторович, ведущий научный сотрудник, к.т.н., НИЦ «Курчатовский институт» - ВИАМ, [email protected] Грунин Алексей Александрович, инженер, НИЦ «Курчатовский институт» - ВИАМ, [email protected]
Шавнев Андрей Александрович, заместитель начальника НИО, к.т.н., НИЦ «Курчатовский институт» - ВИАМ, [email protected]
Information about the authors
Denis V. Sidorov, Leading Researcher, Candidate of Sciences (Tech.), NRC «Kurchatov Institute» -VIAM, [email protected]
Aleksey A. Grunin, Engineer, NRC «Kurchatov Institute» - VIAM, [email protected]
Andrew A. Schavnev, Deputy Head of Scientific-Research Bureau, Candidate of Sciences (Tech.), NRC «Kurchatov Institute» - VIAM, [email protected]
Статья поступила в редакцию 25.07.2023; получена после доработки 17.08.2023;
одобрена и принята к публикации после рецензирования 17.08.2023. The article was submitted 25.07.2023; received in revised form 17.08.2023; approved and accepted for publication after reviewing 17.08.2023.