Скорость химической реакции образования 1,2-диметилдисилена из метилсилилена

Д. В. Сидоров, А. А. Грунин, А. А. Шавнев
УДК
544.43
DOI
10.18577/2307-6046-2025-0-5-76-87
версия в формате pdf (929.23 КБ)
Цитирование
Д. В. Сидоров, А. А. Грунин, А. А. Шавнев Скорость химической реакции образования 1,2-диметилдисилена из метилсилилена // Труды ВИАМ. 2025. № 5. DOI: 10.18577/2307-6046-2025-0-5-76-87. URL: https://test.viam.ru/journal/2025/5/6
Ключевые слова
химическая реакция, скорость, 1, 2-диметилдисилен, метилсилилен, квантовая химия, термодинамические параметры, 1
Аннотация

Теоретическим методом квантовой химии выполнено исследование химической реакции образования 1,2-диметилдисилена из метилсилилена, определены геометрические параметры и термодинамические характеристики структурных изомеров 1,2-диметилдисилена. Установлено, что реакция образования 1,2-диметилдисилена из метилсилилена является экзотермической и протекает самопроизвольно при температурах менее 1570 К, при этом скорость исследуемой химической реакции значительно увеличивается с уменьшением температуры.

Введение

Понимание фундаментальных основ протекания физических и химических процессов является базовой составляющей при разработке и изучении материалов, технологий их получения и переработки [1, 2].

Ранее считалось, что образование двойной связи между атомами кремния в кремнийорганических соединениях практически невозможно. Существование короткоживущих кремнийорганических соединений, содержащих двойную связь у атома кремния, доказано в экспериментальных работах, представленных в обзорах [3, 4]. В сообщении [5] описан стабильный при комнатной температуре силаэтен с температурой плавления 95 °С, устойчивый в атмосфере аргона и реагирующий с кислородом воздуха. Позднее в работах [6, 7] синтезированы и изучены методами спектроскопии три новых стабильных дисилена с мезитильными органическими группами (например, 1,3,5-триметилбензена) и установлено довольно легкое взаимное превращение цис- и транс-стереоизомеров полученных соединений.

При помощи метода лазерного фотолиза в работе [8] изучали реакционную способность 1-метилсиленов с метильными, этильными, трет-бутильными, винильными, этинильными, фенильными, триметилсилильными и триметилсилилметильными группами в углеводородных растворителях. Источником метилсиленов в данной работе являлись соответствующие 1-метилсилациклобутаны, которые под действием лазерного воздействия распадались на этилен и силен с оставшейся органической группой. Установлено, что реакционная способность силена возрастает с увеличением электронодонорной и индуктивной способности заместителей у кремния, но при этом существенно зависит от пространственного расположения органических групп. В дальнейшем лазерный фотолиз использован для исследования реакционной способности более чем 30 силенов, синтезированных из кремнийорганических соединений, таких как арил-, винил-, алкинилдисиланы, силациклобутены и R-силилкетены [9, 10]. В работе [11] изучена реакционная способность 1,1-дифенил-2-неопентилсилена в различных растворителях методами лазерного флэш-фотолиза с использованием в качестве предшественника транс-1,1,3,3-тетрафенил-2,4-динопентил-1,3-дисилациклобутана. Полученный силен имел жизненный цикл ~250 мкс в сухом очищенном от примесей кислорода гексане, после чего димеризовался по механизму «голова к хвосту» со скоростью ~5·108 М–1·с–1. В той же работе установлено, что полученный 1,1-дифенил-2-неопентилсилен активно реагирует с кислородом со скоростью ~6,5·105 М–1·с–1. При фотолизе циклотетрасиланов в метилциклогексане обнаружены дисилены, циклотрисиланы и силилены с жизненным циклом до 50 нс [12]. Наряду с фотолизом, силилены можно получить и термическим методом [13].

По причине короткого жизненного цикла нестабильных кремнийорганических соединений с двойными связями, изучение их структуры, свойств и стабильности проводится как экспериментально, так и теоретически – при помощи неэмпирических методов квантовой химии [14]. Проведенное квантово-химическое исследование бизамещенных силенов с помощью метода с базисным набором B3Lyp/cc-pVDZ показало, что эффект обратной поляризации π-электронов может быть использован для повышения кинетической стабильности силенов изогнутой пространственной структуры [15]. Эффект поляризации заряда по направлению к кремнию π-донорными заместителями у атома углерода также описан в работе [16]. В исследовании [17] по стабилизации молекул силабензолов методом с базисным набором B3Lyp/6-31+G(d, p) показано уменьшение частично положительно заряженного кремния за счет π-электронов от орто- и пара-заместителей, таких как аминогруппы. Так, в случае 2,4,6-триаминозамещения в силабензоле преобладают цвиттер-ионные резонансные структуры, которые придают кремнию отрицательный заряд и стабильность такого силабензола как в виде мономера, так и в виде димера. Присоединение воды и кинетику растворения непредельных соединений кремния исследовали методом с базисным набором B3Lyp/6-31G(d) в работе [18]. Установлено также, что полярность силена можно использовать для прогнозирования кинетической стабильности по отношению к присоединению нуклеофильных групп. Исследование механизмов циклоприсоединения ацетилена к силену и гермену проводили теоретическими методами B3Lyp/6-311 + G(d, p) и CCSD/6-311 + G(d, p) в работе [19].

Следует отметить, что теоретическое обсуждение стабильности молекулы метилсилилена по сравнению с силаэтиленом длится на протяжении нескольких десятилетий [20]. Рассматривается также вопрос о наиболее энергетически выгодных синглетных и триплетных состояниях силиленов. Характерным (по сравнению с карбенами) является основное синглетное состояние для силенов. Этот факт объяснен в работе [21], в которой связывают преобразование из σ2-конфигурации нижнего синглета в σ, π-конфигурацию нижнего триплета с увеличением одноэлектронной энергии за счет перехода неспаренного электрона с s-орбитали (для всех нелинейных силиленов) на свободную p-орбиталь. При этом наблюдается уменьшение двухэлектронной энергии за счет разделения неспаренных электронов при переходе от синглета к триплету. Для образования триплетных карбенов и их аналогов (в том числе силиленов) необходимо продвижение электрона из s2p2-конфигурации основного электронного состояния в sp3-конфигурацию, тогда как образование сиглетного состояния этого не требует. В теоретической работе [22] предсказано триплетное основное состояние для соединения (t-Bu3Si)2Si. Позднее в работе [23] уточненным методом функционала плотности (DFT) показано, что три-третбутилсилил(триизопропилсилил)силилен и ди(триизопропил)силилен имеют синглетное основное состояние, тогда как трет-трибутилсилилсилилен – триплетное. Для четырехчленных циклосилиленов установлено синглетное состояние при теоретическом исследовании их химической стабильности [24]. Относительно новое сравнение [25] синглетных и триплетных состояний для силилена и карбена выполнено на уровне теории статической электронной корреляции (CASSCF) с использованием больших гауссовских базисных наборов, включающих функции f-типа. Результаты расчетов показывают, что триплетное состояние для карбена более энергетически выгодно на 42 кДж/моль, по сравнению с синглетным. Для силилена ситуация обратная – энергия триплетного состояния выше синглетного на 84 кДж/моль. Объяснением служит то, что помимо граничных электронов на синглет-триплетный барьер влияют остовные электроны и энергия ядерного отталкивания. В расчете, представленном в работе [26], показано, что синглет-триплетная энергетическая щель и энергетическая щель между верхней занятой и нижней свободной орбиталями уменьшаются при увеличении нуклеофильности функциональной группы. Для оценки влияния заместителей на реакционную способность силенов в работе [27] проведен теоретический анализ реакций димеризации силенов, изомеризации силилена в силен. Установлено, что положение заместителя имеет решающее значение для полярности связи C‒Si, а также реакционной способности силенов. Димеризация силенов протекает по радикальному пути, а наличие π-донорных заместителей обеспечивает лучшую стабилизацию молекулы. В работе [28] исследовали равновесную смесь, состоящую из тетрасилилдисилена и его изомерного бис(силил)силилена ‒ первого изолированного силилена этого типа. Высокая реакционная способность полученного силилена продемонстрирована на реакциях присоединения водорода и аммиака.

В процессе термического распада кремнийорганических соединений, таких как циклические и линейные алкилсиланы и алкилхлорсиланы, образуется метилсилилен, который в дальнейшем может участвовать в реакциях радикального присоединения [29]. В частности, две молекулы метилсилилена, представляющие собой бирадикалы, при присоединении друг к другу могут образовать соединение с двойной связью между атомами кремния (1,2-диметилдисилен). Метилсилилен также обнаружен в результате взаимодействия кремния и метана [30]. При изучении изомеров в системе «метилсилилен‒силаэтилен‒силилкарбен» установлено, что наиболее стабильным является метилсилилен [31, 32]. В связи с отсутствием в научно-технической литературе данных об образовании 1,2-диметилдисилена из метилсилилена представляет интерес исследование данной химической реакции при помощи теоретического метода квантовой химии.

 

Материалы и методы

Оптимизация пространственного расположения атомов и расчет геометрических параметров молекул метилсилилена, цис-1,2-диметилдисилена и транс-1,2-диметилдисилена выполнены неэмпирическим методом квантовой химии по теории возмущений Меллера–Плессета второго порядка (MP2). В качестве расчетного базиса использован функционал Маклина‒Чандлера (6-311G) с добавлением функций диффузионных и поляризационных составляющих. Вычисления для вышеперечисленных химических соединений проводили в синглетном квантовом состоянии до достижения стационарных точек, характеризуемых отсутствием мнимых частот и минимумом электронной энергии. Термодинамические параметры также определены на уровне теории с базисным набором MP2/6-311 + G(d, p).

Поиск структуры переходного состояния проводили путем изменения значений длин связей и углов между атомами в процессе превращения цис-1,2-диметилдисилена в транс-1,2-диметилдисилен. Оптимизация найденной структуры переходного состояния и расчет термодинамических характеристик проведен вышеуказанным теоретическим методом квантовой химии.

Работа выполнена с использованием оборудования ЦКП «Климатические испытания» НИЦ «Курчатовский институт» – ВИАМ.

 

Результаты и обсуждение

Образование молекулы 1,2-диметилдисилена, содержащей двойную связь между атомами кремния, происходит из двух молекул метилсилилена по химической реакции (1). При сближении в пространстве молекул метилсилилена отсутствует переходное состояние, в отличие от химической реакции образования силаэтилена из метилсилана [33]:

CH3SiH + CH3SiH → CH3SiH SiHCH3.                                         (1)

Строение 1,2-диметилдисилена предопределяет возможность существования его структурных изомеров, таких как цис-1,2-диметилдисилен и транс-1,2-диметилдисилен. На рис. 1 представлены структурные формулы этих изомеров.

 

 

Рис. 1. Структурные формулы изомеров 1,2-диметилдисилена

Геометрические параметры цис-1,2-диметилдисилена и транс-1,2-диметилдисилена, вычисленные методом квантовой химии по теории MP2/6-311+G(d, p), представлены в табл. 1.

 

Таблица 1. Геометрические параметры структуры изомеров 1,2-диметилдисилена при температуре 298К

Геометрические параметры

Значения параметров для

цис-1,2-диметилдисилена

транс-1,2-диметилдисилена

Длина связи C1–Si5, нм

0,1883

0,1881

Длина связи Si5–Si7, нм

0,2168

0,2167

Длина связи C1–Н2, нм

0,1093

0,1093

Длина связи Si5–Н6, нм

0,1477

0,1478

Угол H3–C1–Н4, градус

108,4

108,5

Угол C1–Si5–Si7, градус

121,6

121,5

Угол C1–Si5–H6, градус

113,1

113,3

Угол C1–Si5–Si7–С9, градус

33,8

180,0

 

Видно, что геометрически структуры изомеров отличаются незначительно для длины связей: кремний‒углерод, кремний‒кремний и кремний‒водород. Основное отличие заключается в значении двугранного (диэдрального) угла С1–Si5–Si7–C9, который составляет 33,8 градуса для цис-1,2-диметилдисилена и 180 градусов для транс-1,2-диметилдисилена.

Для оценки термодинамических параметров протекания химической реакции образования 1,2-диметилдисилена из метилсилилена при разных температурах определены изменения полной электронной энергии (∆E), энтальпии (∆H), энергии Гиббса (∆G) и энтропии (∆S). Результаты расчетов представлены в табл. 2.

 

Таблица 2. Изменение термодинамических характеристик химической реакции образования 1,2-диметилдисилена при разных температурах

Продукт реакции

∆E

∆H

∆G

∆S, Дж/(моль·К)

кДж/моль

При температуре 298 К

Цис-1,2-диметилдисилен

–237,93

–240,41

–191,66

–163,49

Транс-1,2-диметилдисилен

–237,95

–240,42

–191,58

–163,82

При температуре 1000 К

Цис-1,2-диметилдисилен

–224,80

–233,11

–81,06

–152,04

Транс-1,2-диметилдисилен

–224,81

–233,12

–80,77

–152,35

При температуре 2500 К

Цис-1,2-диметилдисилен

–189,94

–210,73

135,63

–138,54

Транс-1,2-диметилдисилен

–189,97

–210,75

136,39

–138,86

 

Представленные в табл. 2 результаты расчета изменения полной электронной энергии и энтальпии при образовании цис-1,2-диметилдисилена и транс-1,2-диметилдисилена имеют отрицательные значения, что свидетельствует о выделении тепла в ходе данных химических реакций. При этом с повышением температуры тепловой эффект снижается. Сравнение полученных значений изменения термодинамических характеристик для цис- и транс-изомеров 1,2-диметилдисилена показывает незначительную разницу при одной и той же температуре. Анализируя изменение энергии Гиббса, можно сделать вывод, что химические реакции протекают самопроизвольно при температурах 298 и 1000 К, тогда как при температуре 2500 К их протекание затруднительно. Если допустить, что изменение энергии Гиббса в зависимости от температуры изменяется линейно, то графическим методом можно установить температуру 1570 К, ниже которой реакции протекают самопроизвольно.

Для расчета скорости химической реакции образования 1,2-диметилдисилена из метилсилилена использованы математические формулы [34, 35] для константы равновесия K и константы скорости k, выраженные через изменение энергии Гиббса.

Константа равновесия вычислена по уравнению (2) через величину изменения энергии Гиббса:

 

                                                            (2)

где ∆G – изменение энергии Гиббса; R – газовая постоянная; Т ‒ температура.

 

Скорость протекания химической реакции рассчитана по уравнению

 

                                                          (3)

 

где kb – постоянная Больцмана;h – постоянная Планка.

 

Полученные значения констант равновесия и скоростей образования изомеров 1,2-диметилдисилена из метилсилилена представлены в табл. 3.

 

Таблица 3. Значения констант равновесия и скоростей образования изомеров 1,2-диметилдисилена из метилсилилена при разных температурах

Температура, К

Значения параметров для

цис-1,2-диметилдисилена

транс-1,2-диметилдисилена

K

k, см3/(молекула∙с)

K

k, см3/(молекула∙с)

298

3,94·1033

2,44·1046

3,82·1033

2,37·1046

1000

1,72·104

3,57·1017

1,66·104

3,45·1017

2500

1,47·10–3

7,63·1010

1,41·10–3

7,36·1010

 

Полученные значения из табл. 3 можно использовать для определения скорости исследуемой химической реакции по уравнению Аррениуса (4), расчета предэкспоненциального множителя A и энергии активации (Ea):

                                                           (4)

Проведение расчетов в температурных пределах от 1000 до 2500 К позволяет вывести уравнения (5) и (6) скоростей образования цис-1,2-диметилдисилена и транс-1,2-диметилдисилена соответственно:

 

                                                  

 

Из полученных математических формул за энергию активации химической реакций образования 1,2-диметилдисилена из метилсилилена можно принять значение ‒212,8 кДж/моль. Отрицательное значение энергии активации связано с тем, что данная химическая реакция является экзотермической. На рис. 2 представлена температурная зависимость логарифма скорости реакции образования 1,2-диметилдисилена. Из-за близких значений для цис- и транс-изомеров 1,2-диметилдисилена кривые скоростей сливаются на данном рисунке, а рассчитанные значения практически совпадают со значениями по выведенным уравнениям, за исключением значения при температуре 298 К, отклонение которого составляет ~6 %. Обозначено также граничное значение для начала самопроизвольного протекания реакции, исходя из отрицательного значения изменения энергии Гиббса. Из полученных графических данных видно, что скорость протекания исследуемой химической реакции увеличивается с уменьшением температуры.

 

Рис. 2. Температурная зависимость логарифма скорости реакции образования 
1,2-диметилдисилена (▲ ‒ рассчитанные значения)

 

В связи с тем, что полученные значения полной электронной энергии для цис- и транс-изомеров 1,2-диметилдисилена имеют разницу не более 0,02 кДж/моль, а скорость их образования сопоставима, можно сделать вывод об их одновременном образовании из метилсилилена в одинаковых количествах. В дальнейшем для упрощения расчетов предлагается принять одинаковые значения полной электронной энергии для цис- и транс-изомера. Следует отметить, что в химии углеводородов энергетически более выгодна транс-форма изомеров. Для количественного определения высоты энергетического барьера, необходимого для перехода цис-1,2-диметилдисилена в транс-1,2-диметилдисилен, проведено сканирование изменения полной электронной энергии при изменении двугранного (диэдрального) угла С1–Si5–Si7–C9. Сканирование проведено начиная со значения 33,8 градуса двугранного угла С1–Si5–Si7–C9, характерного для стационарного состояния транс-1,2-диметилдисилена, до угла 180 градусов для цис-1,2-диметилдисилена. По результатам сканирования обнаружена структура переходного состояния (значение 55,6 градуса двугранного угла С1–Si5–Si7–C9) с характерной единственной мнимой частотой. Используя оптимизацию найденной структуры переходного состояния методом MP2/6-311 + G(d, p), определены термодинамические характеристики при разных температурах. На рис. 3 в графическом виде представлены изменение полной электронной энергии для перехода транс-1,2-диметилдисилена в цис-1,2-диметилдисилен и оптимизированная структура переходного состояния.

Рис. 3. Изменение полной электронной энергии для перехода транс-1,2-диметилдисилена в цис-1,2-диметилдисилен

 

Исходя из полученных данных видно, что высота энергетического барьера, рассчитанная на основе изменения электронной энергии, составляет 100,18 кДж/моль при температуре 298 К и 84,26 кДж/моль при температуре 2500 К. Процесс транс-цис-перехода является эндотермическим, т. е. при повышении температуры требуется большее количество энергии для такого перехода. Изменение энергии Гиббса также увеличивается с ростом температуры и имеет значение >104 кДж/моль, что позволяет сделать вывод о невозможности самопроизвольной цис-транс-изомеризации 1,2-диметилдисилена.

 

Заключения

Методом квантовой химии проведено исследование химической реакции образования 1,2-диметилдисилена из метилсилилена, в ходе которой возможно существование структурных изомеров в виде цис-1,2-диметилдисилена и транс-1,2-диметилдисилена.

На основе рассчитанных термодинамических характеристик установлено, что реакция образования 1,2-диметилдисилена из метилсилилена является экзотермической и протекает самопроизвольно при температурах 1570 К. Скорость исследуемой химической реакции значительно увеличивается с уменьшением температуры, энергия активации в количественном выражении составляет около ‒212,8 кДж/моль.

Выявлено, что образование изомеров 1,2-диметилдисилена происходит в одинаковых количествах с высокой скоростью. Высота энергетического барьера для перехода транс-1,2-диметилдисилена в цис-1,2-диметилдисилен составляет 100,18 кДж/моль при температуре 298 К и 84,26 кДж/моль при температуре 2500 К.

Литература
  1. Каблов Е.Н., Хмелева К.М., Заварзин С.В., Козлов И.А., Лонский С.Л. Влияние тер-мической обработки на характеристики алюмоцинковых покрытий, полученных методом ХГН // Авиационные материалы и технологии. 2022. № 1 (66). Ст. 07. URL: http://www.journal.viam.ru (дата обращения: 01.11.2024). DOI: 10.18577/2713-0193-2022-0-1-78-91.
  2. Каблов Е.Н., Шульдешов Е.М., Петрова А.П., Лаптева М.А., Сорокин А.Е. Зависимость комплекса свойств звукопоглощающего материала типа ВЗМК от концентрации гидрофобизирующего состава на основе кремнийорганического герметика // Авиационные материалы и технологии. 2020. № 2 (59). С. 41–49. DOI: 10.18577/2071-9140-2020-0-2-41-49.
  3. Avakyan V.G., Guselnikov S.L., Guselnikov L.E. Classical planar doubly bonded Si-substituted silenes ‒ a boundary system between olefins and heavier group 14 analogs geometric and electronic structures, rotational barriers about the Si=C double bond, and com-parison with the analogs and isoelectronic phosphenes, further evidence against C-ylides of silicon // Journal of organometallic chemistry. 2003. Vol. 686. P. 257–271.
  4. Gusel’nikov L., Nametkin N. Formation and properties of unstable intermediates
  5. containing multiple p-p bonded group 4B metals // Chemical reviews. 1979. Vol. 79. P. 529–577.
  6. Brok A., Abdesaken F., Gutekunst B., etc. A solid silaethene: isolation and characterization // Journal of the chemical society, chemical communications. 1981. Vol. 110. P. 3877–3923.
  7. West R. The disilenes: chemistry of silicon-silicon double bonds // Pure apply chemistry. 1984. Vol. 56. P. 163–173.
  8. West R. Chemistry of the silicon-silicon double bond // Angewandte chemie. 1987. Vol. 26. P. 1201–1211.
  9. Leigh W., Boukherroub R., Kerst K. Substituent effects on the reactivity of the silicon−carbon double bond. Resonance, inductive, and steric effects of substituents at silicon on the reactivity of simple 1-methylsilenes // Journal of the American Chemical Society. 1998. Vol. 120. P. 9504–9512.
  10. Baines K. Brook silenes: inspiration for a generation // Chemical communication. 2013. Vol. 49. P. 6366–6369.
  11. Morkin T., Leigh W. Substituent effects on the reactivity of the silicon-carbon double bond // Accounts of Chemical Research. 2001. Vol. 34. P. 129–136.
  12. Owens T., Grinyer J., Leigh W. Direct detection of 1,1-diphenyl-2-neopentylsilene and the effects of solvent polarity on its reactivity with nucleophiles // Organometallics. 2005. Vol. 24. P. 2307–2318.
  13. Shizuka H., Murata K., Aray Y. et al. Photolysis of cyclotetrasilanes // Journal of the chemical society, Faraday transactions. 1989. Vol. 85. P. 2369–2379.
  14. Чернышев Е.А., Комаленкова Н.Г. Генерирование и реакции силенов в газовой фазе // Успехи химии. 1990. Т. 59. С. 918–930.
  15. Fischer R., Power P. π-Bonding and the lone pair effect in multiple bonds involving heavier main group elements: developments in the new millennium // Chemical reviews. 2010. Vol. 110. P. 3877–3923.
  16. Cheng M., Chu S. Bonding pattern continuum from covalent to dative carbon-silicon bonds for substituted silenes: a theoretical study // Organometallics. 2005. Vol. 24. P. 3746–3752.
  17. Ottoson H. Zwitterionic Silenes: Interesting goals for synthesis? // Chemistry – A European journal. 2003. Vol. 9. P. 4144–4155.
  18. El-Nahas A., Johansson M., Ottosson H. Reverse Si‒C bond polarization as a means for stabilization of silabenzenes: a computational investigation // Organometallics. 2003. Vol. 22. P. 5556–5566.
  19. Bendikov M., Quadt S., Rabin O. et al. Addition of nucleophiles to silenes. A theoretical study of the effect of substituents on their kinetic stability // Organometallics. 2002. Vol. 21. P. 3930–3939.
  20. Pavelka L., Hanson M., Staroverov V. et al. Mechanism of the addition of alkynes to silenes and germenes: a density functional study // Canadian journal of chemistry. 2014. Vol. 93. P. 132–142.
  21. Gaspar P., West R. Silylenes // PATAI's chemistry of functional groups. John Wiley&Sons, 1998. P. 2463–2568.
  22. Kutzelnigg W. Chemical bonding in higher main group elements // Angewandte chemie. 1984. Vol. 23. P. 272–295.
  23. Holthausen M., Koch W., Apeloig Y. Theory predicts triplet ground-state organic silylenes // Journal of the American chemical society. 1999. Vol. 121. P. 2623–2624.
  24. Yoshida M., Tomaoki N. DFT study on triplet ground state silylenes revisited:  the quest for the triplet silylene must go on // Organometallics. 2002. Vol. 21. P. 2587–2589.
  25. Ayoubi-Chianeh M., Kassaee M. Stable four‐membered cyclosilylenes at theoretical levels // Journal of the Chinese chemical society. 2021. Vol. 68. P. 541–550.
  26. Apeloig Y., Pauncz R., Karni M. et al. Why is methylene a ground state triplet while silylene is a ground state singlet? // Organometallics. 2003. Vol. 22. P. 3250–3256.
  27. Ayoubi-Chianeh M. Theoretical study of new stable triplet silylenes // Silicon. 2022. Vol. 14. P. 8965–8969.
  28. Kizhuvedath A., Mallikasseri J., Mattew J. Electronic effects of substituents on the reactivity of silenes: a computational analysis // Structural chemistry. 2024. Vol. 35. P. 119–133.
  29. Reiter D., Holzner R., Porzelt A. et al. Disilene–silylene interconversion: a synthetically accessible acyclic bis(silyl)silylene // Chemical reviews. 2019. Vol. 141. P. 13536–13546.
  30. Sidorov D.V., Kirilin A.D., Shavnev A.A. Rates of formation of silaethylene and the methyl and silyl radicals by the pyrolysis of methylsilane // Theoretical foundations of chemical engineering. 2024. Vol. 56. P. 554–559.
  31. Maier G., Reisenauer H., Glatthaar J. Reactions of silicon atoms with methane and silane in solid argon: a matrix-spectroscopic study // Chemistry – А European journal. 2002. Vol. 8. P. 4383–4391.
  32. Gordon M. The methylsilylene‒silaethylene‒silylcarbene isomerization // Chemical physics letters. 1978. Vol. 54. P. 9–13.
  33. Сидоров Д.В., Кирилин А.Д., Грунин А.А., Шавнев А.А. Структурная изомеризация системы «метилсилилен‒силаэтилен‒силилкарбен» // Труды ВИАМ. 2023. № 9 (127). Ст. 03. URL: http://www.viam-works.ru (дата обращения: 18.11.2024). DOI: 10.18577/2307-6046-2023-0-9-32-40.
  34. Сидоров Д.В., Кирилин А.Д., Шавнев А.А., Мелентьев А.А., Флотский А.А., Грунин А.А. Переходное состояние в химической реакции образования силаэтилена из метилсилана // Труды ВИАМ. 2022. № 9 (115). Ст. 09. URL: http://www.viam-works.ru (дата обращения: 18.11.2024). DOI: 10.18577/2307-6046-2022-0-9-111-120.
  35. Foresman J.B., Frish A.E. Exploring chemistry with electronic structure methods. 3rd ed. Wallingford, CT: Gausssian Inc., 2015. 335 p.
  36. Чоркендорф И., Наймантсведрайт Х. Современный катализ и химическая кинетика. Долгопрудный: Интеллект, 2010. 504 с.