CC BY
85
УДК 681.518.3:541.18
АВТОМАТИЗИРОВАННАЯ СИСТЕМА РАСЧЁТА ТЕРМОДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК АДСОРБЦИИ
© 2008 Л. С. Зеленко1, В. В. Варфоломеева1, В. Ю. Кузнецов1,
А. В. Федькаев1, А. В. Терентьев1, А. К. Буряк2
1 Самарский государственный аэрокосмический университет 2Институт физической химии и электрохимии РАН, г. Москва
Разработана автоматизированная система расчёта термодинамических характеристик адсорбции (ТХА), позволяющая моделировать процесс физической адсорбции для молекул разных классов органических соединений на поверхности графитированной термической сажи (ГТС). Система позволяет проводить расчёты для молекул, обладающих способностью к внутреннему вращению, при всех возможных конформациях. Реализована возможность построения 3Б-модели, позволяющая проверить правильность ввода исходных геометрических параметров молекулы, а также экспорт полученных данных в табличный процессор MS EXCEL для дальнейшей их обработки. Возможности системы продемонстрированы на примере практически важных молекул в сопоставлении с экспериментальными данными.
Термодинамические характеристики, константа генри, теплота адсорбции, адсорбция, графитиро-ванная термическая сажа, атом-атомные потенциалы, эксперимент, молекулярно-статистический метод расчета, моделирование, автоматизированная система
Введение. Автоматизированная система расчёта термодинамических характеристик адсорбции, таких, как константы Генри (^1) и теплоты адсорбции (д), предназначена для исследования адсорбционного и хроматографического поведения соединений разных классов и имеет важное научное и практическое применение в физической и аналитической химии. В качестве объектов исследования выбраны соединения, являющиеся компонентами или микрокомпонентами углеводородных топлив или продуктами их трансформации. Разработка системы является актуальной, так как расширяет область применения молекулярно-статистического метода расчёта ТХА на новые классы органических веществ.
Накопление больших массивов данных по ТХА позволит рассматривать большие ряды гомологов и изомеров. Внутри этих рядов для конкретных молекул можно будет рассматривать адсорбцию всех конформаци-онных изомеров. Поскольку конформацион-ных изомеров и изомеров положения очень много, особенно для углеводородов, то такие расчёты будут важны для исследования свойств углеводородных топлив.
Для расчёта термодинамических характеристик адсорбции необходимо задавать и варьировать различные геометрические параметры молекулы, в первую очередь межатомные расстояния г, валентные углы а и торсионные (двухгранные) углы т (рис. 1). Кроме того, необходимо задавать и в некоторых случаях варьировать и физико-химические параметры, такие, как атомная поляризуемость, диамагнитная восприимчивость и равновесное расстояние.
В настоящее время для расчёта одной конформации (ввода данных, сведения результатов в таблицы) при использовании известной программы расчёта ТХА [1] в среднем затрачивается 2-4 минуты. В случае поворота вокруг одной оси вращения на 360° с шагом 10° необходимо провести расчёты 36 конформаций. Если в молекуле, к примеру, три оси вращения, то необходимо рассчитать 46 6 5 6 (363) конформаций.
На практике часто встречаются молекулы с числом осей вращения больше трёх, для которых провести варьирование углов по всем осям ранее не представлялось возможным из-за временных ограничений. Кроме того, автоматизация расчётов позволяет
2G8
Рис. 1. Основные структурные параметры молекул
уменьшить шаг варьирования. Это, в свою очередь, позволяет оптимизировать геометрическую структуру молекулы и получать корректные атом-атомные потенциалы (ААП), особенно для молекул с неизвестной структурой при количественном сопоставлении их с экспериментальными данными. В случае изомеров только с помощью расчётов ТХА на основе реальных данных о структуре молекул можно провести их надежную идентификацию. Такое исследование выполнено для изомеров метилтриазола, азотсодержащих продуктов трансформации несимметричного диметилгидразина (НДМГ) [2].
Математическое описание молекулярно-статистической теории адсорбции. Главной причиной развития и совершенствования молекулярно-статистических расчётов на ГТС является то, что адсорбция на её поверхности чрезвычайно чувствительна к структуре молекул и позволяет разделять геометрические изомеры [3, 4]. Наибольшее развитие молекулярно-статистические расчёты получили для системы «графитированная термическая сажа - углеводороды и их производные» [5]. В работах [5, 6] разработан полуэмпирический молекулярно-статистический метод расчёта термодинамических характеристик физической адсорбции газов на ГТС при предельно низких степенях заполнения её поверхности, позволяющий проводить исследования поведения молекул, в том числе изомерных, при адсорбции на ГТС. На основании расчётов можно предсказать, взаимодействие какого изомера с поверхностью ГТС будет сильнее и, соответственно, в каком порядке эти соединения будут выходить из колонки с ГТС.
При расчётах сделано допущение, что молекулы являются квазижесткими, так как у большинства молекул амплитуды колебаний ядер весьма малы по сравнению с межъядер-ными расстояниями.
В рамках молекулярно-статистического метода расчёты константы Генри для адсорбции квазижестких молекул приводят по уравнению
к=—гг
4р
/ \1/ 2
ґ2жкТ 4
Ф.
ехр
Ф
кТ
8Іп , (1)
где Ф0 и Ф ”г - значения потенциальной функции Ф взаимодействия между молекулами адсорбата и адсорбента в потенциальном минимуме и её второй производной по расстоянию 2 центра масс молекулы от поверхности адсорбента; © и - углы Эйлера, определяющие ориентацию молекулы относительно поверхности адсорбента.
Величины Ф” Ф0 и равновесное расстояние 20 зависят от © и Ф. Потенциальную функцию Ф представляют в виде суммы ААП взаимодействия каждого атома адсорбируемой молекулы (А) с каждым атомом углерода ГТС (Стс):
Ф
(2)
А С Г
Молекулярно-статистический расчёт базировался на использовании полуэмпири-ческих ААП в форме Бакингема-Корнера [7]:
Фа-а^ =-С'_6 - С8^-8 + Ве^, (3)
где С6, С8 - параметры потенциала диполь-дипольного и диполь-квадрупольного дис-
персионного взаимодеиствии; г - межатомное расстояние; В и q - параметры потенциала отталкивания.
Суммирование р по атомам углерода ГТС над базисной гранью графита проводят с использованием приближения Крауэлла. Потенциальную функцию Ф взаимодействия /-атома молекулы с базисной гранью графита рассчитывают по уравнению
+ -
2Бж$
Ч
ІЧ2і + 1)е
— ,4
а
- С„ 4 с
8 за6
+
(4)
где $ =0,382-10-3 нм-2 - число атомов углерода, приходящееся на единицу площади базисной плоскости графита; ё - межплоскостное расстояние в графите, равное 0,356 нм; £ - дзета-функция Римана; 1. - расстояние от рассматриваемого /-атома до плоскости, проходящей через центры атомов наружного атомного слоя графита.
В выражении (4) учитывается вклад в энергию отталкивания только от наружного атомного слоя графита.
Значения параметров ААП взяты из работ [7, 8] и приведены в таблице 1.
В работе рассматривается нелокализо-ванная адсорбция (свободное движение по математически однородной поверхности).
Таблица 1. Значения параметров ААП при взаимодействии различных атомов с атомом углерода ГТС
Элемент Параметр ААП
^ кДж ■ м Сб, моль п кДж ■ м С8, моль В кДж моль
Водород -0.5672 -0.06403 1618.000
Углерод (£р ) -1.5750 -0.14530 8491.000
Углерод №2) -1.6853 -0.15547 9085.370
Кислород -1.1194 -0.11005 3250.015
Сера -3.8672 -0.45069 26759.806
Межмолекулярное взаимодействие адсорбат-адсорбат не учитывалось, расчёты проводились для области Генри.
Разработанная система позволяет проводить расчёты для всех значений углов внутреннего вращения с любым шагом варьирования в автоматическом режиме.
Описание автоматизированной системы расчёта ТХА. Описание работы системы рассмотрим на примере молекулы 2-фенилэтанола С6Н5-СН2-СН2-ОН, структура которой представлена на рис. 2.
Для решения поставленных задач пользователю в системе предоставляются
следующие функции: возможность задания параметров молекулы и параметров для расчёта ААП сортов атомов (рис. 3), проведение расчётов ТХА (теплота адсорбции, константа Генри) и расчёта положения молекулы в пространстве.
Ввод данных в систему возможен как в режиме диалога, так и через загрузку файла определённого формата. Обмен данными че-
рез файловую систему будет полезен при «стыковке» данной системы с другими системами моделирования процессов физической адсорбции, а также для создания «библиотеки» с данными о молекулах и её изомерах, наиболее часто используемых в экспериментах. Входные данные представлены в виде таблицы, в первом столбце которой находятся порядковые номера атомов молеку-
Рис. 3. Режим ввода параметров молекулы
лы. В рассматриваемом случае (рис. 2) первые шесть номеров - это атомы углерода бензольного кольца. Номера «7» и «8» - Бр3-гиб-ридизованные атомы углерода. Девятый атом
- кислород. Атомы «10» - «14» - атомы водорода, связанные с углеродом бензольного кольца. «15», «16» и «17», «18» - атомы водорода, связанные с С7 и С8, соответственно. «19» - атом водорода спиртовой группы ОН.
Во 2-ом, 3-ом и 4-ом столбцах записываются соответственно номера атомов «I», «К», «Ъ». В пятый столбец заносятся значения длин связей между атомами «1-1». В шестой - значения угла, образованного атомами «I—I—К», в седьмой - значения торсионного угла «I—I—К—Ъ», в 8 столбец - атомные
массы. В последнем девятом столбце можно задать цвета атомов (для отображения 3Б-модели молекулы).
В нижней таблице рисунка 3 приведены характеристики каждого сорта атома. Первые шесть атомов (нумерация в 5-ом столбце) углерода находятся в Бр2-гибридизован-ном состоянии; атомы «7» и «8» - в Бр3-гиб-ридизованном состоянии углерода; «9» - атом кислорода; «10» - «19» - атомы водорода. В другие столбцы (1-4) вводятся параметры ААП для каждого сорта атомов.
На панели ввода данных для расчёта термодинамических характеристик расположены параметры, разъяснение которых представлено в таблице 2.
Таблица 2. Физический смысл параметров
Параметр Физический смысл
NT Число узлов и весов
NP
D Параметр решётки адсорбента
R Длина интервала поиска минимума потенциальной функции
P Число шагов поиска минимума потенциальной функции
ST Числа симметрии
SPS
T1 Начальное значение температуры
T2 Шаг температуры
T3 Конечное значение температуры
После задания входных данных можно проводить расчёты. В системе предусмотрено два вида расчётов: расчёт ТХА заданной молекулы (рис. 4) и её пространственного расположения (рис. 5). Все результаты представляются в виде таблиц на соответствующих вкладках интерфейса. Кроме того, результаты расчётов можно сохранить в файлах формата *.xls или экспортировать их в табличный процессор MS Excel для последующей обработки.
На рис. 4 в верхней таблице представлены результаты расчёта ТХА для каждого значения температуры с заданным шагом. В нижней таблице приведены значения потенциальной функции, т. е. исходные данные для определения ТХА.
Результаты расчётов пространственного расположения молекулы удобно представлять в графическом виде (рис. 6). Поэтому в системе предусмотрена возможность трёхмерного отображения - вывод 3Б-модели
^ Расчет термоаинлмическии хдрдк.герметик дасорбки
фдй/i Реаак.гмрсе4»ме Расевлать моаег*> Справка
Конфигурацию)»** мнтегра/<1-2) КонФ1<гчратмс»*(ый імгеграл(1=3) Константа Гец»< | Константа ГенрЦК.1 ] Изостермми
1 -0.310335874250335 7 54897404406967 7.23368658900928 101626378288729 1502341402:
2 •0.36065321941747 7 3915506823483 165198851509367 6 73509043478456 14 84861888
3 ■0414969222421221 7 30680081920646 1 00024405051917 4 21687736283459 14 64483742
4 ■047361227561306 7 26370423376317 -1 03138408433435 230243226363745 14 41572360
5 -0 536944625066522 7 26558859021561 ■263071703098192 0809059832Є47712 14 16155679
il
+
II 6459675569065
■114489353094358 11.0399744941373 ■10,4038638294885 9 63670072884408 -8,92424894535931 8.43776948697648 •8.2$66$08100079$ 8.38758794801707 8.58703642040422
11 6410327526366 ■11.4220549659558 ■10.9726331954262 -10.2830282024771 9,46074693774868 -8.69805044767369 -8.1 $357872422182 •7.945410734$! 984 8.04090302674535 8.29877984800379
11 632657426854 •11.378712494589 10.8608304717428 10,0879199697117 9.18658012446657 8.35768219512668 7.76134814523114 -7.47678503128182 7.51568307899297 7.79580729413409
116218303638415 • 11.318799444945 •10,7220303999974 9 85706450560712 -8,88187453147778 •8.003447893595
7 36565659169796 *7.03071140472388 •7 00856384429041 7.25142043438919
11.6101954596296 •11.2579114410671 •10.5826552637888 9,64211091715502
8 62669161763812 •7.74367497552737 -7.11483794954328 •$.78211$Э0034357 $.74420643584828 6 9522733421775
11 60010011075* •11.20724509069 •10.47544082644 9.496718463211 8.488214180918 7.652863979677 7.091196228235 -6.823373319203 ■6.820331614671 6.982536974634 .
/Й
Расчет термодннамичес^му характеристик, ддссрбцмн
Файл Редактирование Рассчитать ^модель Справка
ВхОЙЛЫС ленные [Термош*нмінче«* ир-кн | Информация о молек^пе |
№ ІТвмії КмминатаХ КломнпйТіУ км скіпка г
1 0.947331 Э6533&133 -1.34233061743175 0.166763361792329
г 22275Ш5Б1291Є4 -0,902566451235401 □,37і2ЄадВ733041
3 2,45699751495774 0552941930349445 0,106942546144444
4 1.40692597266933 1.36809609273094 ■0,332110163384066
5 0.127157276712334 0527721316542591 -0,506623350325573
6 -0.102339676353264 -0.527736433542751 -0.242133827736332
7 -1,43163415256203 -1,114467526271032 -0,4316^237344041
0 -2,54074333332519 -0,0251594669346953 -0,171561002517932
9 ■2,04372594132969 0.910433940413027 0,79475149041337
10 0.770731459222294 -2,33960948743753 0.4007165084332103
11 3.03763237105775 -1.4313319645212 0.703302730369Э37
12 3.44 41 7СИ 3133763 0.963794255623366 0.241467004641565
13 1 ,вЄЗЄЄ&ЗГ97®2іЄ 2,4131763352904716 -0,5361063Ё6524Е41
14 ■0,6029751330532716 1,45654742092339 ■0,945244011962474
15 -1,6374534213606 -1,93682079722512 0.268765540371762
16 -1.53406586923727 -1.43332554271273 -1.45204476680773
17 -3.45751315135533 -0.435344025555961 0.21575636336236
іе -г,77зоеізй>4га7 0,493375331065615 -1,10369649105906
15 '2,2707433436317' 0,62196961:1351057 1,66969230126
Рис. 5. Декартовые координаты атомов молекулы
молекулы. В режиме отображения молекулы доступны опции: центрирование молекулы -разворот молекулы; отображение нумерации атомов молекулы; отображение осей координат; произвольное вращение молекулы в пространстве; изменение масштаба. Это даёт возможность проверить правильность ввода
исходных геометрических параметров молекулы.
Результаты расчёта и их обсуждение.
По экспериментальным данным работы [9] в газовой фазе молекула этилбензола стабилизируется в антиконформации, с углом ^ = 90°. В силовом поле ГТС она может переходить
из стабильной неплоской конформации в плоскую и обратно. Из работы [10] известно, что этилбензол при адсорбции на ГТС рассматривается как смесь конформацион-ных изомеров. В таблице 3 приведены значения константы Генри для молекулы этилбен-зола при адсорбции на ГТС в интервале температур 300...500 К.
Из таблицы видно, что литературные значения константы Генри попадают в диапазон рассчитанных значений для крайних
торсионных углов молекулы, что подтверждает правильность проведённых расчётов. Кроме того, экспериментальные данные работы [5] и расчётные значения для неплоской конформации молекулы с углом внутреннего вращения т1=30° имеют удовлетворительное совпадение (рис. 7). В условиях хроматографического эксперимента константы Генри, рассчитанные для значения торсионного угла т1=30°, наилучшим образом совпадают с экспериментом [11].
Таблица 3. Значения 1пК1 (мкм) для молекулы этилбензола при адсорбции на ГТС в интервале температур 300.500 К
Этилбензол С6Н5-СН2- О К и)
торсионный угол 0° торсионный угол 90° литературные данные [5]
Т, К 1пК1, мкм 1пК1, мкм 1пК1, мкм
300 8.53 7.01 7.26
350 5.40 4.13 4.52
400 3.10 2.03 2.47
450 1.35 0.44 0.87
500 -0.01 -0.79 -0.41
Рис. 7. Сравнительная диаграмма расчётных (1 - при Т1=0°, 3 - при т1=90°, 4 - при т1=30°) и экспериментальных (2) значений константы Генри для молекулы этилбензола
Молекула этилбензола на поверхности ГТС в плоской (А) и неплоской (Б) конформациях изображена на рис. 8.
В работах[12, 13]установлено, что для молекулы 2-фенилэтанола энергетически наиболее выгоден конформер с внутримолекулярной водородной связью. Сохраняется ли такое внутримолекулярное взаимодействие при адсорбции на ГТС, может быть выяснено путём варьирования геометрии и сопоставления рассчитанных теплот адсорбции и энергии внутримолекулярной водородной связи. Полученные значения ТХА для разных конформаций молекулы в диапазоне темпе-
ратур 300-500 К позволили установить наличие внутримолекулярной водородной связи у 2-фенилэтанола в силовом поле ГТС.
Наиболее наглядно результаты варьирования торсионных углов т и т2 молекулы 2-фенилэтанола можно представить в виде трехмерной диаграммы (рис. 9). Сравнивая такие диаграммы для различных молекул, в том числе изомерных, можно предсказать порядок их выхода из колонны с ГТС и, соответственно, идентифицировать молекулы.
При использовании автоматизированной системы расчёта для построения многомерных диаграмм требуется всего несколько
Рис. S. Молекула этилбензола на поверхности ГТС: плоская (А) и неплоская (Б) (с углом внутреннего вращения 90о) конформации
Рис. 9. Трёхмерная диаграмма зависимости 1пК1 от торсионных углов т и т2 при тз=00 для молекулы 2-фенилэтанола при адсорбции на ГТС (при Т=300 К)
минут. Это позволяет идентифицировать отдельные изомеры в их сложных смесях и прогнозировать порядок выхода веществ из хроматографической колонки на основе величин удерживания.
Заключение. Приведённые примеры не охватывают все имеющиеся и доступные в перспективе возможности автоматизированной системы расчёта. Расширенный вариант программы будет содержать графический модуль, позволяющий визуализировать геометрию молекулы в адсорбированном состоянии после варьирования всех углов внутреннего вращения. Такая визуализация будет доступна для любого конформационного состояния молекулы. В дальнейшем при использовании программы для интерпретации экспериментальных данных и предсказаний величин удерживания изомерных молекул будет поддерживаться автоматизированное построение многомерных диаграмм (углы внутреннего вращения - ТХА) и их сравнение для изомеров разных классов, практически важных в химии ракетных и реактивных топлив.
Библиографический список
1. Афреймович, А. Я. Молекулярно-статистический расчёт адсорбции углеводородов на графите [Текст] : дисс. канд. хим. наук. -М.: МГУ, 1972.
2. Варфоломеева, В. В. Обеспечение экологической безопасности при использовании гидразинных горючих в ракетно-космической отрасли [Текст] / В. В. Варфоломеева, А. В. Терентьев, А. В. Ульянов, А. К. Буряк// Труды XI Всероссийского конгресса «Проблемы выживания человека в техногенной среде современных городов». - Самара, 2006. - С. 57-60.
3. Буряк, А. К. Сочетание масс-спектрометрического и молекулярно-статистического методов для идентификации изомеров [Текст] // Изв. АН СССР. Сер. Химическая. -1990. - № 9. - С. 1995-2000.
4. Буряк, А. К. Применение молекулярно-статистических расчётов для предсказания хроматографического разделения изомерных дифтордифенилов [Текст] / А. К. Буряк,
А. В. Ульянов // Изв. академии наук. Сер. Химическая. - 1996. - №3. - C. 623 - 626.
5. Авгуль, Н. Н. Адсорбция газов и паров на однородных поверхностях [Текст] /
Н. Н. Авгуль, А. В. Киселев, Д. П. Пошкус. -М.: Химия, 1975.
6. Киселев, А. В. Межмолекулярные взаимодействия в адсорбции и хроматографии [Текст] - М.: Высшая школа, 1986.
7. Киселев, А. В. Молекулярные основы адсорбционной хроматографии [Текст] / А. В. Киселев, Д. П. Пошкус, Я. И. Яшин. -М.: Химия, 1986.
8. Буряк, А. К. Применение молекулярно-статистических методов расчёта термодинамических характеристик адсорбции при хромато-масс-спектрометрической идентификации органических соединений [Текст] // Успехи химии. - 2002. - Т. 71. - № 8. -
С. 788.
9. Sun S., Bernstein E. R. Spectroscopy of neurotransmitters and their clusters. 1. Evidence for five molecular conformers of phenylamine in a supersonic jet expansion [Текст] // J. Am. Chem. Soc. 1996. V. 118. P. 5086.
10. Грумадас, А. Ю. Определение потенциальной функции внутреннего вращения молекул из газо-хроматографических данных [Текст] : автореферат дисс. канд. хим. наук. -М.: МГУ, 1983.
11. Kalashnikova E. V, Kiselev A. V., Petrova R. S., Shcherbakova K. D., Poshkus D. P. Chromatographic measurements and molecular statistical calculations of thermodynamic characteristics of adsorption of aromatic and polycyclic hydrocarbons on graphitized thermal carbon black [ Текст] // Chromatographia. 1979. V. 12. N 12. P. 799.
12. Brown R. D., Godfrey P. D. Detection of a higher energy conformer of 2-phenylethanol by millimeter-wave spectroscopy // J. Phys. Chem. A. 2000. V. 104. P. 5742.
13. Dickinson J. A., Hockridge M. R., Kroemer R. T. et al. Conformational choice, hydrogen bonding, and rotation of the S1^!S2 electronic transition moment in 2-phenylethyl alcohol, 2-phenylethylamine, and their water clusters [Текст] // J. Am. Chem. Soc. 1998. V. 120. P. 2622.
References
1. Aphreimovitch, A. Ya. Molecular-statistical calculation of hydrocarbon adsorption on graphite. Diss. cand. chem. science. -Moscow: Moscow State University - 1972.
2. Varfolomeyeva, V. V. Provision of ecological safety when using hydrazine fuels in space rocketry. / V. V. Varfolomeyeva, A. V. Terentyev, A. V. Ulyanov, A. K. Buryak // Transaction of the XI All-Russian Congress “Problems of man’s survival in the technogenous environment of modern cities” - Samara, 2006.
- pp. 57-60.
3. Buryak, A. K. Combination of mass-spectrometric and molecular-statistical methods for isomer identification. // Izvestia (News) of USSR Academy of Sciences. Chemistry series.
- 1990. No. 9 - pp. 1995-2000.
4. Buryak, A. K. Application of molecular
- statistical calculation for the prediction of chromatographic separation of isomeric difluorinediphenyls. / A. K. Buryak, A. V. Ulyanov // Izvestia (News) of Academy of Sciences. Chemistry series. - 1996. - No. 3. - pp. 623626.
5. Avgul, N. N. Adsorption of gas and steam on homogeneous surfaces / N. N. Avgul, A. V. Kiselyov, D. P. Poshkus. - Moscow: Khimiya (Chemistry), 1975.
6. Kiselyov, A. V. Intermolecular interactions in adsorption and chromatography -Moscow: Vysshaya Shkola (Higher School), 1986.
7. Kiselyov, A. V. Molecular foundations of adsorption chromatography / A. V. Kiselyov,
D. P. Poshkus, Ya. I. Yashin. - Moscow: Khimiya
(Chemistry), 1986.
8. Buryak, A. K. Application of molecular-statistical methods of calculating thermodynamic characteristics of adsorption in chromatomass-spectrometric identification of organic compounds // Uspekhi khimii (Advances in chemistry) - 2002. - Vol. 71 - No. 8. p. 788.
9. Sun S., Bernstein E.R. Spectroscopy of neurotransmitters and their clusters. 1. Evidence for five molecular conformers of phenylamine in a supersonic jet expansion // J. Am. Chem. Soc. 1996. V. 118. P. 5086.
10. Grumadas, A. Yu. Determination of the potential function of molecules’ internal rotation on the basis of gasochromatographic data: abstract of diss. cand. chem. science. Moscow: Moscow State University. - 1983.
11. Kalashnikova E.V, Kiselev A.V., Petrova R.S., Shcherbakova K.D., Poshkus D.P Chromatographic measurements and molecular statistical calculations of thermodynamic characteristics of adsorption of aromatic and polycyclic hydrocarbons on graphitized thermal carbon black [ Текст] // Chromatographia. 1979. V. 12. N 12. P. 799.
12. Brown R.D., Godfrey PD. Detection of a higher energy conformer of 2-phenylethanol by millimeter-wave spectroscopy // J. Phys. Chem. A. 2000. V. 104. P. 5742.
13. Dickinson J.A., Hockridge M.R., Kroemer R.T. et al. Conformational choice, hydrogen bonding, and rotation of the S1<S2 electronic transition moment in 2-phenylethyl alcohol, 2-phenylethylamine, and their water clusters [Текст] // J. Am. Chem. Soc. 1998. V. 120. P. 2622.
COMPUTER-AIDED-DESIGN SYSTEM FOR CALCULATING ADSORPTION THERMODYNAMIC CHARACTERISTICS
© 2008 L. S. Zelenko1, V. V. Varfolomeyeva1, V. Yu. Kuznetsov1,
A. V. Fed’kayev1, A. V. Terentyev1, A. K. Buryak2
1Samara State Aerospace University 2Institute of Physical Chemistry and Electrochemistry, Russian Academy of Science, Moscow
A CAD system for calculating adsorption thermodynamic characteristics has been developed which makes it possible to simulate the process of physical adsorption for molecules of various classes of organic compounds on the surface of graphitized thermal black. The system makes it possible to perform calculations for molecules capable of internal rotation, for all possible conformations. A three-dimensional model has been constructed which makes it possible to check whether the input of the molecule’s initial geometric parameters is correct. It also allows exporting the data obtained to the table processor MS EXCEL for their further treatment. The system’s capabilities are shown with practically important molecules as compared to experimental data.
Thermodynamic characteristics, Henri’s constant, adsorption heat, adsorption, graphitized thermal black, atom-atom potentials, experiment, molecular-statistical calculation method, simulation, CAD system
Информация об авторах
Зеленко Лариса Сергеевна, кандидат технических наук, доцент кафедры программных систем. Самарский государственный аэрокосмический университет имени академика С.П. Королева. Область научных интересов: разработка сложных программных систем, в том числе автоматизированных систем обработки информации и управления, автоматизированных системонаучных исследований, технологии программирования и проектирования, информационные технологии.
Варфоломеева Вера Васильевна, кандидат химических наук, доцент кафедры экологии и безопасности жизнедеятельности. Самарский государственный аэрокосмический университет имени академика С .П. Королева. Область научных интересов: теория химического строения и химической связи, водородные связи, поверхностные явления, адсорбция, кинетика и механизмы реакций.
Кузнецов Валадислав Юрьевич, студент. Самарский государственный аэрокосмический университет имени академика С.П. Королева.
Федькаев Алексей Валерьевич, студент. Самарский государственный аэрокосмический университет имени академика С.П. Королева.
Терентьев Алексей Владимирович, ассистент кафедры экологии и безопасности жизнедеятельности. Самарский государственный аэрокосмический университет имени академика С.П. Королева. Область научных интересов: теория химического строения и химической связи, водородные связи, поверхностные явления, адсорбция.
Буряк Алексей Константинович, доктор химических наук, заведующий лабораторией физико-химических основ хроматографии и хромато-масс-спектрометрии. Институт физической химии и электрохимии имени А.Н.Фрумкина РАН. Область научных интересов: хроматография, хромато-масс-спектрометрия, поверхностные явления, адсорбция.
Zelenko, Larisa Sergeyevna, сandidate of technical science, associate professor, department of software systems Samara State Aerospace University named after academician S. P. Korolyov. Area of research: elaboration of complex software systems including CAD systems of processing information and control, CAD system scientific investigations, technologies of programming and design, information technologies.
Varfolomeyeva, Vera Vasiliyevna, сandidate of chemical science, associate professor, department of ecology and vital activity safety Samara State Aerospace University named after academician S. P. Korolyov. Area of research: theory of chemical structure and chemical bonds, hydrogen bonds, surface phenomena, adsorption, kinetics and mechanisms of reactions.
Kuznetsov, Vladislav Yurievitch, studen. Samara State Aerospace University named after academician S. P Korolyov.
Fed’kayev, Alexey Valerievitch., student. Samara State Aerospace University named after academician S. P Korolyov.
Terentiev, Alexey Vladimirovitch, аssistant, department of ecology and vital activity safety Samara State Aerospace University named after academician S. P. Korolyov. Area of research: theory of chemical structure and chemical bonds, hydrogen bonds, surface phenomena, adsorption.
Buruak, Alexey Konstantinovitch, doctor of chemical science, head of laboratory of physicochemical foundations of chromatography and chromato-mass-spectrometry. Institute of physical chemistry and electrochemistry named after A. N. Frumkin, Russian Academy of Science. Area of research: chromatography, chromato-mass-spectrometry, surface phenomena, adsorption.
