УДК 544.723.2.023.221
М. А. Голубчиков*, А. В. Колнооченко, И. В. Лебедев, Н. В. Меньшутина
Российский химико-технологический университет им. Д.И. Менделеева, Москва, Россия 125480, Москва, ул. Героев Панфиловцев, д. 20, корп. 1 * e-mail: golubchikov.mihail@gmail.com
МОДЕЛИРОВАНИЕ ПРОЦЕССА МОНОМОЛЕКУЛЯРНОЙ АДСОРБЦИИ АКТИВНЫХ ВЕЩЕСТВ В ПОРИСТЫХ ТЕЛАХ КЛЕТОЧНЫМ АВТОМАТОМ МАРГОЛУСА
В ходе данной работы исследовался процесс адсорбции активного вещества на поверхности пористого тела. Представлена модель процесса адсорбции в пористых телах на основе вероятностного клеточного автомата с окрестностью Марголуса. Представлены программная реализация и результаты численного моделирования процесса адсорбции по модели Ленгмюра.
Ключевые слова: клеточный автомат, адсорбция, пористое тело, адсорбционное равновесие
В настоящее время в промышленности широко используются различные пористые материалы в качестве матриц-носителей. Примером такой системы может служить катализатор химической реакции, адсорбированный на поверхности пористой структуры или биоразлагаемый пористый материал с активным веществом в качестве носителя лекарственных средств [1, 2, 3].
Актуальной задачей является моделирование процесса адсорбции активных веществ на поверхности внутри пористого тела с целью определения времени заполнения поверхности пористого тела активным веществом при различных условиях.
Процесс адсорбции на поверхности внутри пористых тел хорошо еще не изучен, так как большое влияние на него оказывают молекулярные взаимодействия, а также структура пористого тела и диффузия веществ внутри. В работе использована модель мономолекулярной адсорбции как первая модель, основанная на простой, физически обоснованной модели процесса адсорбции, и клеточный автомат -простой и эффективный способ моделирования, благодаря своей гибкости он способен довольно точно описывать самые различные процессы в любом размере [4, 5].
Процесс адсорбции смоделирован на основе вероятностного клеточного автомата с окрестностью Марголуса. В статье также рассматривается разработанный
кроссплатформенный программный комплекс, который позволяет моделировать процесс адсорбции веществ в зависимости от внешних условий и структуры пористого тела.
Модель мономолекулярной адсорбции Ленгмюра
В основе работы лежит модель физической адсорбции с образованием монослоя. Физическая адсорбция является обратимым процессом, условие равновесия определяется равными скоростями адсорбции молекул адсорбтива A на
вакантных местах поверхности адсорбента Б* и десорбции — освобождения адсорбата из связанного состояния A — S:
A + S* ^ A-S Уравнение адсорбционного равновесия: _ [л - 5] _ е К ~ [5*][Л] " (1-0)С'
где K - константа адсорбционного равновесия, S]=0 и ^*]=(1-0) - доли поверхности адсорбента, занятые и незанятые адсорбатом, а -
концентрация адсорбтива.
Количественно процесс физической мономолекулярной адсорбции описывается уравнением Ленгмюра:
КС
в =-
1 + КС
Принимая допущение о равновесии между фазами, для адсорбции при постоянной температуры можно использовать
фундаментальное уравнение, изотерму Вант-Гоффа, которое дает следующее выражение, которое часто используется для определения энтальпии ДЫ и энтропии ДБ адсорбции из экспериментальных данных:
-ИТ1п(К) =АН- ГД5
Описание клеточного автомата
Для моделирования процессов массопереноса и адсорбции веществ в пористых телах была разработана следующая модель, основанная на вероятностном клеточном автомате с окрестностью Марголуса:
• Система представлена в виде поля с линейным размером N клеток;
• Каждая клетка имеет четыре клетки-соседа (окрестность фон Неймана);
• Расчет идет итеративно, процессы, протекающие на каждой итерации,
принимаются как идущие в один момент времени;
Каждая клетка может находиться в одном из следующих состояний: о Е - растворитель, среда; о М1...МК - растворенное активное вещество;
о Б - неподвижное вещество (пористое тело); Между веществами существует энергия взаимодействия, которая представлена в виде матрицы взаимодействий; Система разбивается на блоки размером 2x2 клеток;
Используется два разбиения системы на блоки от шага к шагу, таким образом организуется частичное перекрывание блоков,
используемых на соседних шагах, которые чередуются (рис. 1): о четное разбиение; о нечетное разбиение; Задаются правила для отдельных блоков, по которым каждый блок с определенной вероятностью может (рис. 2): о повернуться против часовой стрелки (сс^); о повернуться по часовой стрелке (с^); о остаться в исходном состоянии (пт); Вероятности поворота блока зависят от энергии блока после поворота; Клетки с состоянием «неподвижное вещество» при повороте блока не перемещаются, в этом случае в повороте участвуют оставшиеся клетки блока;
Если при повороте блока клетка с состоянием «растворенное активное вещество» будет занимать пространство клетки с состоянием «неподвижное вещество», то такой поворот считается невозможным;
Клетка с состоянием «растворенное активное вещество» считается адсорбированной, если соседняя клетка находится в состоянии «неподвижное вещество» Таким образом, модель описывает случай мономолекулярной адсорбции.
Рис. 1. Перекрывание блоков при четном и нечетном разбиениях решетки автомата
!
1 2
4 3
Рис. 2. Блок может повернуться против (еете), по часовой стрелке (ете), или остаться в исходном состоянии (пт)
Энергии взаимодействия между веществами определяются по справочным данным либо экспериментальными путем.
На каждом шаге разбиения для каждого блока проводится следующий расчет:
Определяются возможности поворота против часовой и по часовой стрелке. Для каждого доступного варианта вычисляется энергия блока как сумма энергий всех возможных связей внутри блока, а также с соседями из смежных блоков. Энергии связи между клетками с различным состоянием постоянные при постоянных условиях и задаются в виде матрицы:
'еЕЕ еЕМ £ММ £МБ
Вычисляется статическая сумма Ъ при всех возможных поворотах и вероятности PJ поворота блока как отношение энергии блока после поворота к статической сумме:
}
е кт
-Е,
12^ Р]=-е*т
Программная реализация
Такая модель легла в основу программного комплекса, который состоит из нескольких частей: модуль подготовки входных данных, расчетный модуль, модуль обработки расчетных данных и модуль визуализации.
Модуль подготовки входных данных представляет собой интерфейс для задания размеров и пористости твердого тела, количества активного вещества в объеме, матрицы взаимодействия веществ и условий проведения расчета.
Расчетный модуль реализует алгоритм расчета адсорбции в системе до наступления адсорбционного равновесия с помощью вероятностного клеточного автомата с окрестностью Марголуса, сохранение
промежуточных данных, проверку выполнения закона сохранения массы в системе.
Модуль обработки расчетных данных строит кинетику мономолекулярной адсорбции активного вещества на поверхности пористого тела, а также графики Вант-Гоффа в координатах логарифма константы равновесия от обратной температуры.
Модуль визуализации позволяет наглядно проследить за системой в разные моменты времени.
Программный комплекс реализован с помощью кроссплатформенных языков программирования Java, C++ и Python, таким образом нет ограничений по выбору платформы для проведения расчетов.
Результаты расчета, обсуждение и развитие
Был проведен численный расчет и сравнение с экспериментальными данными процесса адсорбции нафталина на поверхности кремниевого аэрогеля в среде сверхкритического диоксида углерода плотностью 403 кг/м3 и адсорбции бекламетозола в среде сверхкритического диоксида углерода плотностью 807 кг/м3.
Размеры поля клеточного автомата 200 на 200 клеток, что соответствует 0.2 на 0.2 мкм или 0.04
ухудшается адсорбция. Такой эффект можно объяснить двумя вещами:
1) С ростом температуры увеличивается внутренняя и кинетическая энергии молекул.
2) Энтропия адсорбции отрицательная, а так как система стремится к минимуму энергии, то с ростом температуры переход вещества в адсорбированное состояние становится все менее выгодным.
Рис. 4. Кинетика адсорбции нафталина с 15% начальной объемной долей и 40°, 60° и 80°С
мкм2. Размеры аэрогеля 100 на 100 клеток, что
соответствует 0.1 на 0.1 мкм или 0.01 мкм2. Пористость аэрогеля 95%. Энтальпия адсорбции нафталина -12 кДж/мольК, энтропия -18Дж/моль. Энтальпия адсорбции бекламетозола -85 кДж/мольК, энтропия -4 Дж/моль.
Е 0.Л
0.С6 0.1 0.15 0.2 0.Э 0.3 О.ЭЬ 04 0Ж 0.5
Рис. 5. Изотерма адсорбции нафталинапри T = 40°, 60° и 80°C
Рис. 3. Кинетика адсорбции нафталина при 60°С и 5%, 15% и 30% начальной объемной доли
Расчеты проводились при температурах 40°, 50°, 60°, 70°, 80°С и с начальными объемными долями растворенного вещества от 2.5 до 50% с шагом 2.5%.
На рис. 3 представлена кинетика адсорбции нафталина при 60° C и 5%, 15% и 30% начальной объемной доли. Можно заметить, что с ростом концентрации увеличивается количество адсорбированных частиц. На рис. 4 представлена кинетика адсорбции нафталина с 15% начальной объемной долей и 40°, 60° и 80°С. Заметно, что повышение температуры негативно сказывается на адсорбции нафталина. На рис. 5 представлена изотерма адсорбции при T=40°, 60° и 80°С Также, как и в предыдущем случае, с ростом температуры
На рис. 6 представлена кинетика адсорбции бекламетозола при 60° C и 5%, 15% и 30% начальной объемной доли. Можно заметить, что с ростом начальной концентрации увеличивается количество адсорбированных частиц. На рис. 7 представлена кинетика адсорбции бекламетозола с 15% начальной объемной долей и 40°, 60° и 80° С. Заметно, что изменение температуры не влияет на адсорбцию бекламетозола. На рис. 8 представлена изотерма адсорбции при T=40°, 60° и 80°С Также, как и в предыдущем случае, изменение температуры практически не влияет на адсорбцию. Такой эффект можно объяснить тем, что бекламетозол очень хорошо адсорбируется, энтальпия адсорбции у него очень большая.
н. =
В
- 5 00тЬ
- liOCW заома
Б о id» гаю зав тоо я*» га»
Рис. 6. Кинетика адсорбции бекламетозола при 60°C и 5%, 15% и 30% начальной объемной доли
-4СС: -WC
вое
ЛЕО за» 4С£Ю ИкршишоБ вдеия
Рис. 7. Кинетика адсорбции бекламетозола с 15% начальной объемной долей и 40°, 60° и 80°С
Вант-Гоффа в этом случае не дает хороших результатов.
Рис. 8. Изотерма адсорбции бекламетозола при Т = 40°, 60° и 80°С
После всех расчетов был построен график Вант-Гоффа, зависимость логарифма константы адсорбционного равновесия от обратной температуры (рис. 9, 10). По полученному графику определены энтальпия и энтропия адсорбции, и они были сравнены с экспериментальными данными (табл. 1).
Полученные данные по адсорбции нафталина с высокой точностью совпадают с экспериментальными, что доказывает
работоспособность модели. А данные по адсорбции бекламетозола не совпадают с экспериментальными, значит, для сильно адсорбирующихся активных веществ не стоит ожидать достоверных данных об энергиях адсорбции, так как адсорбция происходит очень сильно, константа равновесия большая и практически не зависит от температуры. График
£.1 1
0.5
ti;к| = 1454.в51дмеви> ■ R*= аэадввонз
о Ж! мои е.®» и.то о.ом! мозг о.омэ
1Л1*
Рис. 9. График Вант-Гоффа для проведенного численного расчета адсорбции нафталина
4 4G чя) = lil IJUliBSi t AEtiaXmik Ri- DifrKMLllJi ■ Я
4,4
1.34
•
■
_
IS 4ia
о.шге о.от7 скепгв о.оогв екаоз ош о ош скоозэ
Рис. 10. График Вант-Гоффа для проведенного численного расчета адсорбции бекламетозола
Таблица 1.
Нафталин Бекламетозол
ДЫ/R, K ДS/R ДЫ/R, K ДS/R
Эксперимент -1443 -2.165 -10223 -0.48
Модель -1455 -2.156 -172 4.872
Предполагается развитие данного программного продукта: оптимизация кода программного модуля отвечающего за процесс расчета с помощью клеточного автомата, усложнение модели адсорбции, создание сервера с базой данных для задания и хранения расчетных данных и результатов обработки расчетов.
Голубчиков Михаил Александрович, аспирант кафедры кибернетики химико-технологических процессов, РХТУ им. Д.И. Менделеева, Россия, Москва
Колнооченко Андрей Викторович, к. ф.-м. н., н. с. Международного учебно-научного центра фармацевтических и биотехнологий, РХТУ им. Д.И. Менделеева, Россия, Москва
Лебедев Игорь Витальевич, инженер ОАО "ВНИИХТ", Россия, Москва
Меньшутина Наталья Васильевна, д. т. н., профессор кафедры кибернетики химико-технологических процессов, РХТУ им. Д.И. Менделеева, Россия, Москва
Литература
1. Anglin E.J., Cheng L., Freeman W.R., Sailor M.J. Porous silicon in drug delivery devices and materials // Adv Drug Deliv Rev. 2008. Vol. 60 №11. P.1266-1277.
2. Trewyn B.G., Slowing I.I., Giri S., Lin V.S.-Y. Mesoporous silica nanoparticle based controlled release, drug delivery, and biosensor systems // Chem. Comm. 2007. P. 3236-3245.
3. Bjorklund R. B., Zangooie S., Arwin H. Adsorption of Surfactants in porous silicon films // Langmuir. 1997. Vol. 13. № 6. P. 1440-1445.
4. Csahyk Z., Vicsek T. Lattice gas model for collective biological motion // Physical Review A. 1995. Vol. 52. №5. P. 5297-5303.
5. Малинецкий Г.Г., Степанцев М.Е. Моделирование движения толпы при помощи клеточных автоматов // Известия ВУЗов. Сер. Прикладная нелинейная динамика. 1997. Т. 5. С. 75-79.
Golubchikov Mikhail Alexandrovich *, Kolnoochenko Andrey Victorovich, Lebedev Igor Vitalievich, Menshutina Natalia Vasilievna
D.I. Mendeleev University of Chemical Technology of Russia, Moscow, Russia. * e-mail: golubchikov.mihail@gmail.com
MODELING OF MONOMOLECULAR ADSORPTION OF ACTIVE SUBSTANCES IN POROUS SOLIDS WITH MARGOLUS CELLULAR AUTOMATA
Abstract
In the present study, the adsorption of the active substance on the surface of the porous body was investigated. The model of the adsorption process in porous solids is based on probabilistic cellular automaton with Margolus neighborhood. Presented software implementation and results of numerical simulation of the adsorption process by the Langmuir model.
Key words: cellular automata, adsorption, porous body, adsorption equilibrium