CC BY
12
УДК 544.18: 544.43
Г. М. Храпковский, Е. В. Николаева, Д. Л. Егоров,
А. Г. Шамов, Нгуен Ван Бо
БАРЬЕРЫ РЕАКЦИЙ ГАЗОФАЗНОГО МОНОМОЛЕКУЛЯРНОГО РАСПАДА НИТРОТОЛУОЛОВ.
2. ВЛИЯНИЕ МОЛЕКУЛЯРНОЙ СТРУКТУРЫ НА БАРЬЕРЫ РЕАКЦИЙ ВНУТРИМОЛЕКУЛЯРНОГО ПЕРЕНОСА ВОДОРОДА С ОБРАЗОВАНИЕМ АЦИ-НИТРОТОЛУОЛА
Ключевые слова: квантово-химический расчет, нитротолуол, механизмы реакций, термическое разложение.
Представлены результаты квантово-химического исследования механизма реакции внутримолекулярного переноса водорода от метильной группы к нитрогруппе для о-нитротолуола, 2,3-, 2,4-, 2,5-, 2,6-динитротолуолов и 2,4,6-тринитротолуола. Изучение проводилось с использованием методов B3LYP и wB97XD с различными наборами базисных функций. Показано, что основные выводы о влиянии заместителей на изменение в ряду энтальпии активации реакции образования аци-нитросоединений передаются согласованно всеми использованными в работе методами. По оценкам всех методов минимальный барьер активации среди изученных соединений наблюдается для 2,4,6-тринитротолуола (159,0 кДж/моль, B3LYP/6-31+G(2df,p)). Наибольшее значение по данным расчёта методом B3LYP/6-31+G(2df,p) предсказывается для о-нитротолуола (169,3 кДж/моль) или 2,3-динитротолуола (171,4 кДж/моль) или 2,5-динитротолуола (169,8 кДж/моль).
Keywords: quantum-chemical computation, nitrotoluene, mechanisms of reactions, thermal decomposition.
The results of quantum-chemical study for mechanism of reaction of intramolecular transfer of hydrogen from methyl group to nitro-group for o-nitrotoluene, 2,3-, 2,4-, 2,5-, 2,6-dinitrotoluenes and for 2,4,6-trinitrotoluene are presented in this article. The study was performed using B3LYP and wB97XD methods, with various sets of basic functions. It is shown, that basic conclusions about the influence of substituents on the change in row of enthalpy of activation offormation reaction for aci-nitrocompounds are transmitted in concord by all used methods. According to estimation of all methods the minimal barrier of activation for all studied compounds is observed for 2,4,6-trinitrotoluene (159,0 kJ/mole, B3LYP/6-31+G(2df,p)). The maximal value obtained by computation using B3LYP/6-31+G(2df,p) method is predicted for о-nitrotoluene (169,3 kJ/mole),2,3-dinitrotoluene (171,4 kJ/mole) or for 2,5-dinitrotoluene (169,8 kJ/mole).
В предыдущем сообщении [1] мы рассмотрели блюдается для реакции внутримолекулярного перенекоторые общие закономерности различных альтер- носа водорода группы CH3 к группе NO2 с образо-нативных механизмов термического распада нитрото- ванием аци-нитротолуолов. В работе [1] мы привели луолов. Было отмечено, в частности, что для соедине- оценки барьеров этой реакции, полученных с ис-ний, имеющих метильный заместитель в орто- пользованием метода B3LYP/6-31+G(2df,p). положении к нитрогруппе минимальный барьер на-
Таблица 1 - Энтальпии активации (ДН*), энтальпии реакций (ДНреакц) и энтальпии активации обратных реакций (ДН*обр) образования а^и-форм некоторых нитроаренов, имеющих водородсодержащий заместитель в п-положении к нитрогруппе (кДж/моль) по данным метода различных квантово-химических методов
Соединение wB97XD/6-31G(d,p) wB97XD/6-311+G(df,p) wB97XD/6-31+G(2df,p)
AH" АНреакц ДНобр AH" ДНреакц AH" обр AH" AHреакц AH" обр
о-нитротолуол 192,2 177,2 15,0 189,8 168,8 21,0 188,4 170,2 18,2
2,3-динитротолуол 192,0 177,9 14,1 191,5 170,2 21,3 191,0 172,1 18,9
2,4-динитротоулол 186,2 170,1 16,1 185,0 162,2 21,8 183,3 166,1 17,2
2,5-динитротолуол 192,3 170,2 22,1 190,8 162,8 28,0 189,3 164,8 24,6
2,6-динитротолуол 187,6 171,7 15,9 186,6 165,5 21,1 185,6 166,1 19,4
2,4,6-тринитротолуола 178,9 160,6 18,3 176,7 153,2 23,6 178,0 156,6 21,4
Соединение B3LYP/6-31G(d) B3LYP/6-31+G(2df,p) B3LYP/6-311++G(3df,3pd)
AH АНреакц дНобр AH" ДНреакц AH" обр AH" AHреакц AH" обр
о-нитротолуол 184,9 174,6 10,3 169,3 154,6 14,7 169,1 151,5 17,6
2,3-динитротолуол 183,8 171,9 11,9 171,4 154,5 16,9 171,3 151,5 19,8
2,4-динитротоулол 179,0 165,7 13,3 164,3 147,1 17,2 164,1 144,1 20,0
2,5-динитротолуол 184,2 168,0 16,2 169,8 149,1 20,6 169,7 146,3 23,4
2,6-динитротолуол 174,1 159,2 14,9 164,9 146,2 18,8 165,1 143,3 21,8
2,4,6-тринитротолуола 171,1 153,2 17,9 159,0 137,5 21,4 159,0 134,8 24,3
С целью получения дополнительной информации об особенностях механизма реакции образования аци-нитротолуолов, мы провели исследование для наиболее важных в практическом отношении представителей ряда (о-нитротолуола, 2,4-динитротолуола, 2,6-динитотолуола, 2,4,6-триниротолуола), а также 2,3-и 2,5-динитротолуолов, с использованием различных методов и наборов базисных функций. Поскольку целью работы было изучение механизм реакции в ряду соединений достаточно большого размера, при выборе метода мы ориентировались на различные варианты
методов функционала плотности. С учётом литературных данных [2-20] и результатов проведённых предварительных оценок барьеров реакций нитро-толуолов мы выбрали для дальнейших дополнительных исследований метод wB97XD с наборами базисных функций 6-3 Ш^^), 6-311+G(df,p) и 6-31+(2df,p). С целью получения возможно более полных данных был использован также метод B3LYP с существенно различающимися по размеру базисами 6-3Ш^) и 6-311++G(3df,3pd) (табл. 1).
Таблица 2 - Геометрические параметры (длины связи в пм) и энтальпии образования (кДж/моль) исходных молекул (Р), переходных состояний (ПС) и продуктов (ПР) и их изменения (Д) для реакции образования а^и-форм некоторых нитроаренов с, B3LYP/6-31+G(2df,p) (wB97XD/6-31+G(2df,p)):
Р1
Р2 .О!
0^ О2
?! X!
I \\ рЗ0^^5
р1 Н°
О!
6
Н
!
о
2
Р
5
5
5
Р
!
!
Р
Р
Соединение О1О2 О2Ы3 ы3о4 о4н6 н^'1 ^О1
о-нитротолуол Р1=ОИ2, Р2=Р3= Р4=Р5=Н Р 140,6 147,4 122,6 256,8 109,3 150,8 91,8
(140,0) (146,9) (121,7) (258,5) (109,2) (150,5) (138,8)
ПС 144,2 138,7 132,2 113,9 149,7 141,5 261,1
(143,5) (138,0) (130,8) (112,5) (149,9) (140,7) (327,2)
ПР 147,3 134,6 139,4 97,9 210,6 136,7 246,4
(146,8) (133,9) (137,5) (97,4) (211,4) (135,9) (309,3)
2,3-динитротолуол Р1=ОИ2, Р2=Р3= Р4=Н, Р5=Ы02 Р 139,8 148,1 121,8 274,2 109,4 150,8 113,0
(139,2) (147,4) (121,0) (274,0) (109,3) (150,6) (194,6)
ПС 144,3 138,9 131,7 114,8 147,7 141,2 284,4
(143,5) (138,2) (130,2) (113,7) (147,6) (140,3) (385,7)
ПР 147,2 134,2 139,2 98,0 206,5 136,0 267,5
(147,2) (133,2) (137,2) (97,5) (208,3) (135,0) (366,8)
2,4-динитротолуол Р1=ОИ2, Р2=Р3= Р5=Н, Р 140,7 147,8 122,5 257,6 109,4 150,6 80,4
(140,0) (147,3) (121,6) (257,3) (109,2) (150,3) (165,4)
р4=ыо2 ПС 144,3 139,7 131,6 115,2 148,4 141,1 244,8
(143,7) (138,9) (130,2) (113,6) (148,8) (140,3) (348,7)
ПР 147,4 135,5 138,6 98,0 210,6 136,3 227,5
(147,2) (134,3) (136,7) (97,6) (210,1) (135,4) (331,5)
2,5-динитротолуол Р1=ОИ2, Р2=Р4= Р5=Н, Р3=Ы02 Р 140,5 147,8 122,4 259,9 109,3 150,8 81,6
(139,8) (147,3) (121,5) (258,3) (109,2) (150,5) (166,5)
ПС 143,9 139,2 131,8 114,8 148,6 141,5 251,3
(143,2) (138,5) (130,4) (113,3) (149,0) (140,6) (355,8)
ПР 147,2 134,7 139,5 97,7 216,4 136,6 230,7
(146,5) (134,5) (137,7) (97,2) (220,8) (136,6) (331,2)
2,6-динитротолуол ^=ОН2, Р3=Р4= Р5=Н, Р 140,5 147,8 122,4 284,8 109,3 150,8 100,5
(139,8) (147,2) (121,5) (284,5) (109,2) (150,5) (182,3)
р2=ыо2 ПС 145,0 139,7 131,5 114,9 148,0 140,9 265,5
(144,4) (138,9) (130,1) (113,3) (148,6) (140,2) (367,8)
ПР 148,1 134,9 139,1 97,9 214,4 135,6 246,7
(147,8) (133,6) (137,3) (97,4) (216,4) (134,8) (348,5)
2,4,6-тринитротолуол Р1=ОИ2, Р3=Р5=Н, Р2= Р4=Ы02 Р 140,6 148,1 122,2 284,9 109,4 150,6 102,0
(139,9) (147,4) (121,4) (284,5) (109,2) (150,3) (222,5)
ПС 145,1 140,8 130,9 116,2 146,8 140,5 261,0
(144,5) (140,0) (129,5) (114,4) (147,6) (139,8) (400,5)
ПР 148,2 135,8 138,4 98,0 214,9 135,4 239,6
(147,8) (134,5) (136,6) (97,5) (216,2) (134,6) (379,2)
Основные выводы о влиянии заместителей на изменение в ряду энтальпии активации реакции образования аци-нитросоединений передаются согласованно всеми использованными в работе методами, что
говорит о надежности полученных результатов. Так, например, интервал различий изменения в ряду энтальпий активации, по данным использованных нами разных методов и базисов не превышает 12-14
кДж/моль (табл. 1). По оценкам всех методов минимальный барьер активации наблюдается для 2,4,6-тринитротолуола.
Наибольшее значение по данным расчёта предсказывается для о-нитротолуола или 2,3-динитротолуола или 2,5-динитротолуола. Во всех случаях, увеличение размера базиса уменьшает абсолютное значение барьера активации.
Объяснить наблюдаемые изменения энтальпии активации позволяет анализ изменений геометрических параметров переходного состояния по сравнению с исходной молекулой (Д/, табл. 2 и 3) в ряду изученных нитротолуолов. При увеличении количества нит-рогрупп в ряду молекул о-нитротолуол, 2,4-динитротоулол, 2,6-динитротолуол, 2,4,6-
тринитротолуол изменения длин связей Д/"(02-М3),
Д/"(М3-04), Дг^-Н6) реакционного центра уменьшаются. Противоположная тенденция наблюдается
12 11 для изменений длин связей Д/(0 -0 ), Дг(Р -0 ).
Первая тенденция преобладает, поэтому мы и наблюдаем уменьшение энтальпии активации при увеличении в молекулах о-нитротолуола, 2,4- и 2,6-динитротолуолов, 2,4,6-тринитротолуола количества нитрогрупп. Однако наличие противоположной тенденции приводит к тому, что эффект уменьшения барьера активации для процесса образования аци-нитросоединений выражен очень слабо. Кроме того, малая величина изменения барьера активации в ряду изученных соединений связана ещё и с тем, что изменения геометрических параметров в ряду
/ 4 _ прям _ | .тах . .т1пч
(Д/ряд - |Д/(пс-р)| - |Д/(пс-р)| ) не превышают 2,5 пм (табл. 3).
Таблица 3 - Изменения геометрических параметров (длин связей) реакционного центра (Дг, пм) при протекании реакции образования а^и-форм некоторых нитроаренов, В3ЬУР/6-31+ОДО1",р) ^В97ХБ/6-31+0^,р)) (см. табл. 2)
Соединение Изменение 0102 02Ы3 м3о4 04Н6 Н6^1
о-нитротолуол Р1-0Н2, Р2-Р3- Р4-Р5-Н Д/(пс-р) 3,6 -8,7 9,6 -142,9 40,4 -9.3
(3,5) (-8,9) (9,1) (-146,0) (40,7) (-9.8)
Д/(пс-пр) -3,1 4,1 -7,2 16,0 -60,9 4.8
(-3,3) (4,1) (-6,7) (15,1) (-61,5) (4.8)
2,3-динитротолуол Р1-0Н2, Р2-Р3- Р4-Н, р5-ыо2 Д/(пс-р) 4,5 -9,2 9,9 -159,4 38,3 -9.6
(4,3) (-9,2) (9,2) (-160,3) (38,3) (-10.3)
Дг(ПС-ПР) -2,9 4,7 -7,5 16,8 -58,8 5.2
(-3,7) (5,0) (-7,0) (16,2) (-60,7) (5.3)
2,4-динитротолуол Р1-0Н2, Р2-Р3- Р5-Н, Р4-Ы02 Д/(ПС-Р) 3,6 -8,1 9,1 -142,4 39,0 -9.5
(3,7) (-8,4) (8,6) (-143,7) (39,6) (-10.0)
Д/(пс-пр) -3,1 4,2 -7,0 17,2 -62,2 4.8
(-3,5) (4,6) (-6,5) (16,0) (-61,3) (4.9)
2,5-динитротолуол Р1-0Н2, Р2-Р4- Р5-Н, Р3-Ы02 Д/(пс-р) 3,4 -8,6 9,4 -145,1 39,3 -9.3
(3,4) (-8,8) (8,9) (-145) (39,8) (-9.9)
Д/(пс-пр) -3,3 4,5 -7,7 17,1 -67,8 4.9
(-3,3) (4,0) (-7,3) (16,1) (-71,8) (4.0)
2,6-динитротолуол Р1-0Н2, Р3-Р4- Р5-Н, р2-ыо2 Д/(пс-р) 4,5 -8,1 9,1 -169,9 38,7 -9.9
(4,6) (-8,3) (8,6) (-171,2) (39,4) (-10.3)
Дг(ПС-ПР) -3,1 4,8 -7,6 17,0 -66,4 5.3
(-3,4) (5,3) (-7,2) (15,9) (-67,8) (5.4)
2,4,6-тринитротолуол Р1-0Н2, Р3-Р5-Н, Р2- Р4-Ы02 Д/(пс-р) 4,5 -7,3 8,7 -168,7 37,4 -10.1
(4,6) (-7,4) (8,1) (-170,1) (38,4) (-10.5)
Д/(пс-пр) -3,1 5,0 -7,5 18,2 -68,1 5.1
(-3,3) (5,5) (-7,1) (16,9) (-68,6) (5.2)
В ряду рассматриваемых соедине- |Д/(пс-Р)|тах 4,5 9,2 9,9 169,9 40,4 10.1
ний (4,6) (9,2) (9,2) (171,2) (40,7) (10.5)
Кг т1п ДГ(пс-р) 3,4 7,3 8,7 142,4 37,4 9.3
(3,4) (7,4) (8,1) (143,7) (38,3) (9.8)
Кг тах ДГ(пс-пр) 3,3 5,0 7,7 18,2 68,1 5.3
(3,7) (5,5) (7,3) (16,9) (71,8) (5.4))
Кг т1п Дг(ПС-ПР) 2,9 4,1 7,0 16,0 58,8 4.8
(3,3) (4,0) (6,5) (15,1) (60,7) (4.0)
Д/рядПРЯМ-|Д/(ПС-Р)|таХ - |Д/(ПС-Р)|т'п 1.1 (1.2) 1,9 (1,8) 1,2 (1,1) 27,5 (27,5) 3,0 (2,4) 0,8 (0,7)
Дгряд0бр -|Д/(ПС-ПР)|таХ - |Д/(ПС-ПР)|т'П 0.4 (0.4) 0,9 (1,5) 0,7 (0,8) 2,2 (1,8) 9,3 (11,1) 0,5 (1,4)
Здесь надо отметить, что в анализе мы не рассматриваем образующуюся в ходе реакции связь
46
О -Н , т.к. в исходном соединении она отсутствует, а на больших расстояниях взаимодействие атомов существенного вклада не вносит. Однако это рас-
стояние, вероятнее всего сказывается на увеличении энергии активации реакции образования аци-форм в 2,5-динитротолуоле по сравнению с о-нитротолуолом, 2,4-динитротоулолом, 2,6-динитротолуолом, 2,4,6-тринитротолуолом. Хотя
оно фактически меньше в 2,5-динитротолуоле, чем в 2,6-динитротолуоле и 2,4,6-тринитротолуоле, но в нем нет дополнительного стерического взаимодействия близлежащих заместителей - нитрогрупп в положении 2 и 6 и метильной группы в положении 1 в 2,6-динитротолуоле и 2,4,6-тринитротолуоле.
Как уже отмечалось выше, также максимальное значение барьера активации процесса образования аци-формы наблюдается для 2,3-динитротолуола. Связано это, скорее всего с тем, что в этом случае геометрические параметры переходного состояния по сравнению с исходными молекулами изменяются сильнее, чем в случае какого-либо другого нитротолуола. В результате даже для о-нитротолуола барьер активации рассматриваемого процесса ниже, чем для 2,3-динитротолуола (169,3 и 171,4 кДж/моль по данным метода B3LYP/6-3^^,р) (табл. 1)).
Барьеры обратных реакций (процессов изомеризации аци-нитротолуолов в нитротолуолы) очень малы и изменения химического строения молекул нитротолуолов, в частности накопление нит-рогрупп практически не влияет на величину энтальпии активации. Анализ геометрических параметров реакционного центра в переходном состоянии и аци-нитротолуолах находит этому простое объяснение. Дело в том, что и изменения геометрических параметров в ходе реакций для соединений (АГ(пс-пр)) и изменения геометрических параметров в ряду рассматриваемых соединений (Агрядобр = |АГ(ПС-Р)|тах - |АГ(ПС-Р)|т|п), в этом случае значительно меньше, чем для прямых реакций (А/"рядпрям) образования аци-нитротолуолов.
Структура продуктов реакции изомеризации нитротолуолов значительно ближе по координате реакции к структуре переходного состояния, чем к структуре исходных соединений. Например, длина формирующейся связи 04-Н6 в о-ниротолуоле отличается от её величины в переходном состоянии, по данным метода B3LYP/6-31+G(2df,p), на 142,9 пм; аналогичная величина для обратной реакции составляет только 16 пм (табл. 3).
Кроме того, в ряду изученных соединений изменение длины связи Агрядпрям(О2-Ы3) в ходе реакций образования аци-нитротолуолов составляет 1,9 пм. Для обратной реакции эта величина будет уже почти в два раза меньше (А/"рядобр = 0,9 пм). Для других геометрических параметров наблюдается подобная картина. Например, для связи Р'-О1 указанные выше величины составляют 0,8 пм и 0,5 пм; для связи Ы3-04 они равны 1,2 пм и 0,7 пм соответственно (табл. 3).
Поэтому неудивительно, что в отличие от реакций внутримолекулярного переноса водорода с образованием аци-форм нитротолуолов, которые все являются ярко выраженными эндотермическими процессами, обратные реакции изомеризации нитроновых кислот в нитроарены по данным расчёта должны протекать с выделением значительного количества тепла.
Завершая обсуждение механизма реакции внутримолекулярного переноса водорода метильной группы к нитрогруппе нитротолуолов нужно отме-
тить, что этот процесс, очевидно, является первичным актом мономолекулярного распада о-нитротолуола и его нитропроизводных. Этому процессу соответствует минимальный барьер из всех приведённых в нашем предыдущем сообщении ([1], табл. 3) барьеров реакций мономолекулярного распада нитротолуолов. С другой стороны, оценки барьеров активации указанного процесса, полученные с использованием различных методов и наборов базисных функций, значительно превышают экспериментально наблюдаемые энергии активации нитротолуолов, имеющих водородсодержащий заместитель в орто-положении к нитрогруппе [1]. Отсюда следует вывод, что первичный акт термического распада для этих соединений не является лимитирующей стадией механизма процесса разложения.
Литература
1. Г.М. Храпковский, Е.В. Николаева, Д.Л. Егоров, А.Г. Шамов, Вестник Казанского технологического университета, (2014) (в печати).
2. A.G. Shamov, G.M. Khrapkovskii, Mendeleev Communications, 11, 4, 163-164 (2001).
3. И.В. Аристов, Д.Л. Егоров, А.Г. Шамов, Г.М. Храпков-ский, Вестник Казанского технологического университета, 1, 7-10 (2011).
4. G.M. Khrapkovskii, T.F. Shamsutdinov, D.V. Chachkov, A.G. Shamov, Computational and Theoretical Chemistry, 686, 1-3, 185-192 (2004).
5. Т.Ф. Шамсутдинов, Д.В. Чачков, А.Г. Шамов, Г.М. Храпковский, Известия высших учебных заведений. Серия: Химия и химическая технология, 49, 9, 38-40 (2006).
6. A.F. Shamsutdinov, T.F. Shamsutdinov, D.V. Chachkov, A.G. Shamov, G.M. Khrapkovskii, International Journal of Quantum Chemistry, 107, 13, 2343-2352 (2007).
7. G.M. Khrapkovskii, Mendeleev Communications, 7, 5, 169-171 (1997).
8. Д.Д. Шарипов, Д.Л. Егоров, Д.В, Чачков, А.Г. Шамов, Г.М. Храпковский, Вестник Казанского технологического университета, 7, 45-52 (2010).
9. Д.Л. Егоров, Г.М. Храпковский, А.Г. Шамов, Вестник Казанского технологического университета, 10, 18-21 (2010).
10. Е.В. Николаева. Дис. канд. хим. наук, Казанский гос. технол. ун-т, Казань, 2002. 203 с.
11. Д.В. Чачков. Дис. канд. хим. наук, Казанский гос. технол. ун-т, Казань, 2005. 206 с.
12. В.И. Коваленко, В.Ф. Сопин, Г.М. Храпковский, Структурно-кинетические особенности получения и термодеструкции нитратов целлюлозы. Москва, Рос. акад. наук, Казан. науч. центр, Ин-т орган. и физ. химии им. А. Е. Арбузова, 2005, 213 с.
13. А.Г. Шамов, Е.В. Николаева, Г.М. Храпковский, Вестник Казанского технологического университета, 24,12-20 (2011).
14. А.Г. Шамов, Е.В. Николаева, Г.М. Храпковский, Вестник Казанского технологического университета, 16, 8, 7-15 (2013).
15. Г.М. Храпковский, Д.Л. Егоров, А.Г. Шамов, Вестник Казанского технологического университета, 21, 20-23 (2013).
16. Г.М. Храпковский, Д.Л. Егоров, А.Г. Шамов, Вестник Казанского технологического университета, 22, 13-16 (2013).
17. Д.Л. Егоров, И.В. Аристов, А.Г. Шамов, Г.М. Храпков-ский, Вестник Казанского технологического университета, 5, 7-9 (2013).
18. А.Г. Шамов, Г.М. Храпковский, Вестник Казанского 20. А.Г. Шамов, Е.В. Николаева, Г.М. Храпковский, технолун-та, Спец. вып., 37-44 (2009). ЖПрХ, 82, 10, 1587 - 1608 (2009).
19. Г.М. Храпковский, А.Г. Шамов, Е.В. Николаева, Д.В. Чачков, Успехи химии, 78, 10, 980-1021 (2009).
© Г. М. Храпковский - д.х.н., профессор кафедры катализа КНИТУ, [email protected]; Е. В. Николаева - к.х.н., доцент кафедры катализа КНИТУ; Д. Л. Егоров - к.ф.-м.н., н.с. НИОКХ КНИТУ; А. Г. Шамов - начальник отделения информатизации КНИТУ; Нгуен Ван Бо - аспирант кафедры катализа КНИТУ.
© G. M. Khrapkovskii - dr.Sc. in chemistry, professor, department of catalysis KNRTU, [email protected]; E. V. Nikolaeva - PhD in chemistry, associate professor, department of catalysis KNRTU; D. L. Egorov - PhD in physics-mathematical sciences, researcher, research department of computational chemistry KNRTU; A. G. Shamov - head of department of informatization KNRTU; Nguyen Van Bo - graduate student of department of catalysis KNRTU.
