close

Вход

Забыли?

вход по аккаунту

?

154234

код для вставкиСкачать
На правах рукописи
Мышлявцева Марта Доржукаевна
Математическое моделирование сложных адсорбционных систем
на поверхности твердых тел: метод трансфер-матрицы
05.13.18 «Математическое моделирование, численные методы
и комплексы программ»
Автореферат
диссертации на соискание ученой степени
доктора физико-математических наук
Пермь – 2013
1
Работа выполнена в Федеральном государственном бюджетном образовательном учреждении высшего профессионального образования «Омский государственный технический университет»
Научный консультант: доктор химических наук
Мышлявцев Александр Владимирович
Официальные оппоненты:
Быков Валерий Иванович,
доктор физико-математических наук, профессор,
ведущий научный сотрудник Института биохимической
физики им. Н.М. Эмануэля РАН, г. Москва
Песков Николай Владимирович,
доктор физико-математических наук, ведущий научный
сотрудник факультета вычислительной математики и кибернетики Московского государственного университета
им. М.В. Ломоносова, г. Москва
Спивак Семен Израилевич
доктор физико-математических наук, профессор,
заведующий кафедрой математического моделирования
Башкирского государственного университета, г. Уфа
Ведущая организация: ФГБОУ ВПО «Санкт-Петербурский государственный
университет»
Защита диссертации состоится «25» марта 2014 года в 14-30 на заседании
Диссертационного совета Д 212.188.08 при ФГБОУ ВПО «Пермский национальный исследовательский политехнический университет» по адресу: 614990, г.
Пермь, Комсомольский пр., д. 29, ауд. 423.
С диссертацией можно ознакомиться в научной библиотеке Пермского национального исследовательского политехнического университета и на сайте
http://www.pstu.ru/
Автореферат разослан ________________
Ученый секретарь
диссертационного совета Д 212.188.08,
кандидат физико-математических наук
А.И. Швейкин
2
СПИСОК ИСПОЛЬЗУЕМЫХ СОКРАЩЕНИЙ
АЦ
БТБК
вн.
г.ц.к.
ЛХ
ММК
МРГ
МТМ
СТМ
–
–
–
–
–
–
–
–
–
активный центр
бензол-трибензойная кислота
внутреннее
гранецентрированная кубическая
Ленгмюр-Хиншельвуд
метод Монте-Карло
модель решёточного газа
метод трансфер-матрицы
сканирующий туннельный микроскоп
СОМ
ст.с.
ТДС
ТКМБ
ТМ
ТМК
ЦГД
ЦК
–
–
–
–
–
–
–
–
сложная органическая молекула
стационарное состояние
термодесорбционный спектр
трикарбокси-метокси-бензол
трансфер-матрица
тримезиновая кислота
циклогексадиен
циануровая кислота
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность работы. Различные явления и процессы, протекающие на поверхности твёрдых тел, лежат в основе многих высокотехнологичных производств. Поверхностные явления играют важную роль в электронике, катализе,
материаловедении, космической технике, атомной и тепловой энергетике и т.д. В
настоящее время большое внимание уделяется так называемым двумерным пористым структурам или самоорганизующимся монослоям. Эти слои состоят из
сложных органических молекул (СОМ), регулярным образом расположенных на
поверхности раздела фаз. Наличие двумерных пористых структур позволяет
сформировать на их основе регулярную трёхмерную наноструктуру, которая используется при изготовлении полевых транзисторов, органических светодиодов,
сенсоров и т.д. В теоретических работах, посвящённых исследованию адсорбции
СОМ на поверхности твёрдого тела, используются чаще всего квантовохимические методы, практически не позволяющие исследовать поведение адсорбционного слоя в целом и ограниченные изучением небольшого количества
частиц. При моделировании всего адсорбционного слоя пользуются методами
статистической механики. Распределение частиц в адсорбционном слое зависит
как от латеральных взаимодействий между адсорбированными частицами, так и
от свойств самой поверхности. Современные методы исследования показали, что
часто одновременно функционируют несколько типов активных центров (АЦ),
при этом необходимо учитывать геометрию адсорбированных молекул. При моделировании таких систем необходимо учитывать также запреты на определённые конфигурации расположения адсорбированных частиц. В частности, при
формировании самоорганизующихся монослоёв СОМ обычно занимают несколько АЦ поверхности и в зависимости от внешних параметров, таких как концентрация, температура и др., могут занимать различное количество мест на поверхности. Теоретический анализ таких сложных систем только начинается и, безусловно, может сыграть значительную роль в понимании процессов самоорганизации на молекулярном уровне. Таким образом, разработка методологических подходов к изучению процессов и явлений, протекающих на поверхности твердых
тел, на основе сложных математических моделей адсорбционных поверхностных
3
слоев является актуальной проблемой, имеющей важное научное и практическое
значение.
Работа выполнена при поддержке Федерального агентства по образованию
РФ в рамках проекта по АВЦП «Развитие научного потенциала высшей школы»
(№ 3.11Ф, 2009-2011гг.); Минобрнауки РФ в рамках госзаказа по госбюджетной
НИР (№ 3.5634.2011, 2012-2013гг.); в рамках ФЦП "Развитие обороннопромышленного комплекса РФ на 2007-2010 годы и на период до 2015 г. ("Ресурс", № 11028, 2013г.); ФЦП "Научные и научно-педагогические кадры инновационной России" на 2009-2013 годы (№ 16.740.11.0762, 2011-2013гг.; №
14.В37.21.0946, 2012-2013гг.).
Объектом исследования является неидеальный адсорбционный слой, его
структура и протекающие в нем процессы.
Цель работы состоит в создании теоретических и методологических основ
численного анализа рассматриваемых в работе явлений, позволяющих описывать
и прогнозировать поведение реальных адсорбционных систем с несколькими типами адсорбционных центров и различными способами адсорбции, на основе
разработанных новых математических моделей, численных методов, алгоритмов
и программных комплексов.
Для достижения указанной цели были поставлены следующие задачи:
 Разработать математические модели для описания формирования сложных
адсорбционных слоёв.
 Разработать и теоретически обосновать математические методы для изучения решёточных моделей.
 Расширить применимость метода трансфер-матрицы и модифицировать алгоритмы для исследования решеточных моделей адсорбционных слоёв, учитывающих наличие нескольких типов активных центров в одной элементарной
ячейке.
 Численно и аналитически исследовать влияние латеральных взаимодействий на критические явления (множественность стационарных состояний (ст.с.),
автоколебания) в реакции, протекающей по параллельному адсорбционному механизму гетерогенного катализа.
 Исследовать модель многоцентровой адсорбции самоорганизующихся монослоёв СОМ на поверхности твёрдого тела, учитывающей возможность различной ориентации молекул (относительно поверхности и относительно друг друга)
и неоднородности поверхности, применить разработанную модель для описания
адсорбции ненасыщенных циклических углеводородов на Si(001)-2×1, монослоёв
тримезиновой кислоты.
 Разработать численные алгоритмы и комплекс программ для реализации
построенных моделей сложных адсорбционных слоёв.
4
Методика исследования. Математическое моделирование проводилось с
использованием аппарата теории обыкновенных дифференциальных уравнений,
статистической физики, теории катастроф, химической кинетики, функционального анализа, линейной алгебры, теории вероятностей. Численные результаты
были получены при помощи метода трансфер-матрицы (МТМ), методов РунгеКутты, Адамса-Башфорта, Розенброка, бисекций, степенного метода, метода
Монте-Карло (ММК). Программирование проводилось в среде Visual Fortran.
Положения, выносимые на защиту:
1. Новый детерминистский математический подход для изучения решёточных моделей без трансляционной инвариантности. В рамках данного подхода
модифицирован метод трансфер-матрицы для неоднородных, трансляционноинвариантных систем, включая использование неквадратных матриц. Предложенные подходы могут быть использованы для изучения любых классических
решёточных моделей, в том числе моделей магнетиков изинговского типа, широко используемых при описании магнитных свойств многих реальных материалов.
2. Усовершенствованный качественный метод исследования математической
модели параллельного адсорбционного механизма гетерогенного катализа, учитывающий её специфику и основанный на общих представлениях теории катастроф.
3. Строгое обоснование применимости вычислительных алгоритмов метода
трансфер-матрицы и разработанного математического метода для изучения решёточных систем без трансляционной инвариантности. Результаты, полученные при
помощи разработанных алгоритмов, верифицированы методом Монте-Карло.
4. Комплексы программ, реализующие: а) разработанный математический
метод для изучения решёточных систем без трансляционной инвариантности; б)
численное построение бифуркационных диаграмм на основе развитого качественного метода исследования параллельного адсорбционного механизма; в) алгоритмы метода трансфер-матрицы для неоднородных систем, обладающих трансляционной инвариантностью; г) алгоритм мультипликативного разложения для
исследования сложных решёточных систем с большим числом состояний.
5. Результаты комплексных исследований адсорбционных систем посредством разработанных математических методов, алгоритмов и комплексов программ, в частности:
– вычислены основные термодинамические характеристики адсорбционных
систем без трансляционной инвариантности на основе разработанного метода и
его программной реализации;
– определена роль латеральных взаимодействий в усложнении кинетического поведения параллельного адсорбционного механизма на основе проведенного
его систематического изучения при помощи разработанных алгоритмов и программного комплекса, их реализующих;
5
– с использованием программного комплекса, реализующего алгоритмы
МТМ, вычислены основные термодинамические характеристики систем с несколькими типами АЦ в одной элементарной ячейке, проведено сравнение результатов вычислительного эксперимента с экспериментальными данными;
– вычислены основные термодинамические характеристики самоорганизующихся монослоёв СОМ, занимающих несколько АЦ поверхности и способных к различной ориентации, исследованы сложные решёточные модели с большим числом состояний при помощи МТМ и его программной реализации;
– обнаружены два новых явления: а) немонотонное изменение степени покрытия от химического потенциала; б) новый тип «чёртовой лестницы» фазовых
переходов.
Полученные результаты диссертационного исследования соответствуют следующим пунктам паспорта специальности 05.13.18 - Математическое моделирование, численные методы и комплексы программ по физико-математическим
наукам: п.1 «Разработка новых математических методов моделирования объектов
и явлений», п.2 «Развитие качественных и приближённых аналитических методов
исследования математических моделей», п.3 «Разработка, обоснование и тестирование эффективных вычислительных методов с применением современных
компьютерных технологий», п.4 «Реализация эффективных численных методов и
алгоритмов в виде комплексов проблемно-ориентированных программ для проведения вычислительного эксперимента», п.5 «Комплексные исследования научных и технических проблем с применением современной технологии математического моделирования и вычислительного эксперимента».
Научная новизна основных результатов работы:
– Разработан новый детерминистский математический метод для изучения
решёточных моделей без трансляционной инвариантности, основанный на распространении основных подходов метода трансфер-матрицы. Ранее такие модели
изучались только стохастическими методами. Модифицирован метод трансферматрицы для неоднородных, трансляционно-инвариантных систем, включая использование неквадратных матриц. Впервые модифицированный МТМ использован для изучения модели адсорбции на ступенчатой поверхности.
– Доказаны теоремы, описывающие условия существования области множественности ст.с. параллельного адсорбционного механизма при произвольном
наборе латеральных взаимодействий. Для трансфер-матрицы, соответствующей
МРГ, доказана единственность, положительность и некратность наибольшего по
модулю собственного значения. Построена система нелинейных алгебраических
уравнений, приближённо определяющая большую статистическую сумму, и вероятности различных конфигураций для решёточных моделей без трансляционной инвариантности. Доказаны теоремы о существовании и в ряде случаев единственности решения этой системы нелинейных алгебраических уравнений.
6
– Создан комплекс программ, реализующих разработанный математический
метод для изучения решёточных систем без трансляционной инвариантности, построение бифуркационных диаграмм, алгоритмы метода трансфер-матрицы для
исследования неоднородных, трансляционно-инвариантных систем, сложных
решёточных систем со сложной элементарной ячейкой.
– В рамках построенной обобщённой решёточной модели адсорбции с использованием метода трансфер-матрицы показано, что учёт латеральных взаимодействий приводит к существенному усложнению области множественности ст.с.
и возникновению автоколебаний для классического механизма ЛенгмюраХиншельвуда; МТМ эффективен при изучении решёточных моделей с несколькими типами АЦ в одной элементарной ячейке (H/Pd(100) и CO/Ni(100)). При
учёте многоцентровости и различных способов адсорбции обнаружены два новых
явления: а) немонотонное изменение степени покрытия от химического потенциала; б) новый тип «чёртовой лестницы» фазовых переходов.
– Разработаны и исследованы решёточные модели, качественно описывающие поведение реальных адсорбционных систем, таких, как 1,4-циклогексадиен
на кремнии (001)-2×1, тримезиновая кислота (ТМК) на поверхности переходных
металлов. Полученные результаты соответствуют данным эксперимента.
Практическая значимость работы: Метод трансфер-матрицы продемонстрировал высокую эффективность для исследования моделей с несколькими типами активных центров с учётом латеральных взаимодействий, систем без трансляционной инвариантности, неоднородных, трансляционно-инвариантных систем, моделей самоорганизующихся монослоёв СОМ на поверхности с учётом нескольких форм адсорбции и направленных латеральных взаимодействий. Результаты работы позволяют глубже понять протекание элементарных процессов на
поверхности твёрдых тел, обобщить имеющиеся экспериментальные данные, качественно предсказать фазовое поведение реального адсорбционного монослоя
СОМ, зная геометрию и химическую структуру молекулы адсорбата и поверхности твёрдого тела. Разработанные модели, вычислительные алгоритмы и комплексы программ и результаты математического моделирования можно использовать при прогнозировании свойств СОМ, которые могут быть использованы
при производстве устройств (сенсоры, органические светодиоды и др.) и функциональных материалов со структурированием в нанометровом диапазоне.
Достоверность и обоснованность результатов, полученных в работе, обеспечивается подтверждением строгих аналитических результатов результатами
численного моделирования, качественным согласованием результатов вычислительного эксперимента с имеющимися экспериментальными данными, подтверждением численных результатов, полученных методом трансфер-матрицы, и
данных расчётов по методу Монте-Карло.
7
Апробация работы. Основные положения работы, результаты теоретических исследований и численного моделирования обсуждались на Международных и Всероссийских конференциях: III International Сonference on unsteady-state
processes in catalysis (St. Petersburg, Russia –1998); Сибирский конгресс по прикладной и индустриальной математике (Новосибирск, 1996; 1998; 2000); VII Всероссийская конференция «Механизмы каталитических реакций» (Санкт – Петербург, 2006); III International Conference “Catalysis: fundamentals and application”
dedicated to the 100th anniversary of Academician Georgii K. Boreskov (Novosibirsk,
2007); II Международный форум по нанотехнологиям «РосНаноТех» (Москва,
2009); The International Conference on nanotechnology for green and sustainable construction (Cairo-Egupt, 2010); Всероссийская научная конференция «Химия под
знаком Сигма» (Омск, 2010); International Conference on Chemical Reactors
(«Сhemreactor-18», Malta, 2008; «Сhemreactor-19»,Vienna, Austria, 2010); Nanotech
Conference and Expo 2010 (Anaheim, USA, 2010); 18th International Vacuum Congress (Beijing, China, 2010); Международный симпозиум «Современная химическая физика» (Туапсе, 2008, 2009, 2010, 2011); International Symposium «Surface
heterogeneity effects in adsorption and catalysis on solids» ( ISSHAC-6, Poland,
Zakopane, 2006; ISSHAC-7, Kazimierz Dolny, Poland, 2009; ISSHAC-8, Poland,
Krakow, 2012); Международная научно-техническая конференция «Динамика
систем, механизмов и машин» (Омск, 2009,2012) и др. Полностью работа доложена и обсуждена на семинарах: кафедры математического моделирования систем и процессов ПНИПУ (Пермь, 2013, руководитель - профессор П.В. Трусов),
кафедры механики композиционных материалов и конструкций ПНИПУ (Пермь,
2013, руководитель - профессор Ю.В. Соколкин), Института механики сплошных
сред УрО РАН (Пермь, 2013, руководитель - академик РАН В.П. Матвеенко), Института проблем переработки углеводородов СО РАН (Омск, 2013, руководитель
- член-корреспондент РАН В.А. Лихолобов), кафедры вычислительной математики СПбГУ (Санкт-Петербург, 2013, руководитель - профессор В.М. Рябов); на
расширенном научном семинаре кафедры «Высшая математика» ОмГТУ (Омск,
2013, руководитель - член-корреспондент РАН А.Ю. Веснин).
Личный вклад автора. Автор принимал непосредственное участие в постановке задач, обосновании и формулировке основных положений, определяющих
научную новизну и практическую значимость, в анализе и обобщении результатов, формулировке выводов. Доказательство всех математических утверждений,
разработка вычислительных алгоритмов, создание комплексов программ и проведение вычислительных экспериментов также принадлежат лично автору.
По теме диссертации опубликована 81 печатная работа, из которых 2 монографии, 1 обзор, 32 статьи в изданиях, рекомендованных ВАК (из них 13 статей в изданиях, входящих в базы Web of Science и Scopus), 40 – в сборниках трудов Международных и Всероссийских конференций; зарегистрированы 15 свиде8
тельств на программы для ЭВМ. Список основных публикаций приведён в конце
автореферата.
Структура и объём диссертации. Диссертация содержит введение, шесть
глав, заключение, список литературы (из 499 источников). Объём диссертационной работы составляет 295 страниц, включая 2 таблицы и 171 рисунков.
СОДЕРЖАНИЕ ДИССЕРТАЦИИ
Введение содержит обоснование актуальности выбора темы исследования,
цель и задачи работы, основные результаты, выносимые на защиту, значение полученных результатов, включая практическое, обсуждение научной новизны.
Первая глава содержит обзор литературных данных по методам исследования решёточных моделей, теоретических и экспериментальных работ, посвящённых анализу адсорбции, десорбции, поверхностной диффузии, фазовых диаграмм, критических явлений, многоцентровой адсорбции, самоорганизующихся
монослоёв СОМ. Обоснован выбор в качестве модели адсорбционного слоя модели решёточного газа (МРГ). Традиционно используемые кластерные методы,
такие, как приближение среднего поля, приближение Бёте-Пайерлса, могут приводить к нефизическим результатам. Ренорм-групповые подходы достаточно
сложны и не позволяют построить сходящиеся последовательности приближений. Метод Монте-Карло, несмотря на свою универсальность, не может гарантировать достижения термодинамического равновесия. Метод трансфер-матрицы
лишён всех этих недостатков. Исходя из этого, в качестве основного метода исследования МРГ был выбран метод трансфер-матрицы.
Во второй главе изложены различные вычислительные алгоритмы МТМ.
Трансфер-матрица представляет собой матрицу, элементы которой являются
больцмановскими множителями, соответствующими всем возможным состояниям пары соседних колец (рис. 1). Применение МТМ основано на следующих точных для полубесконечной системы (рис. 1) соотношениях: 1) большой термодиT
ln 1 , соответствующий одному узлу; 2) вероятM
ность нахождения k соседних узлов в состояниях i1 ,...., ik :
намический потенциал   
Pi1 , ..., ik 
v l , i1 Ti1 i2 Ti2 i3 ....Ti( k 1) ik v r , ik
1k 1
,
где 1 – максимальное собственное значение трансферматрицы (ТМ), T – абсолютная температура, Tij – элементы
ТМ, v l , i1 ( v r , ik ) – i1 (ik ) - компоненты левого (правого) собРис. 1. Квадратная
решётка на цилиндре при М = 4
ственного вектора, соответствующего 1 , и нормированного
на единицу. Известно, что МТМ позволяет в рамках МРГ и
теории переходного состояния вычислять концентрационные
9
зависимости кинетических констант. Для трансфер-матрицы, соответствующей
модели решёточного газа, доказана теорема о единственности, положительности
и некратности наибольшего по модуля собственного значения. (Доказательство
основано на теореме Перрона-Фробениуса и общих свойствах МРГ).
Показана применимость МТМ к изучению решёточных систем без трансляционной инвариантности. Для линейной цепочки из N узлов со свободными
граничными условиями (узлы n типов распределены статистически независимо с
вероятностями p1 ,..., pn , соответственно, и каждый узел может находиться в m состояниях) получены системы нелинейных алгебраических уравнений, определяющих большую статистическую сумму:
pj
pi
n
m
m
m
      ik  ,  ik  xik     jl y jilk /   jl  ,
l 1
j 1  l 1
i 1  k 1


n
(1)
где xik  exp(  k   ik ) , yijkl  exp(  ijkl ) ,  ijkl – безразмерная энергия латеральных
взаимодействий i-го и j-го типов узлов в состояниях k и l,  k – безразмерный химический потенциал для k-го состояния узлов,  ik – добавка к химическому потенциалу k-го состояния узла, зависящая от типа узла. Система уравнений (1) позволяет вычислять любые термодинамические характеристики системы. Решение
системы (1) осуществлялось методом итераций. Численные результаты подтверждают гипотезу о единственности решения системы (1) в области  ik  0 , и область сходимости итерационного процесса совпадает с областью  ik  0 . Система
уравнений (1), определяющая  ik , принадлежит к классу операторных уравнений
вида
х=Ах.
Параметры
(1):
xik , yijkl  0;  ,
pi  0 ,
n
 pi  1 ,
i 1
i, j  1, n;
k , l  1, m; n  1, m  2 . Числа  ik образуют элемент z векторного пространства Е
размерности nm. Так как все  ik  0 , то z  K , где K – выпуклый конус. Доказаны
теоремы существования и единственности решения системы (1).
Теорема 1. В конусе K существует хотя бы одна неподвижная точка оператора А, определяемого системой (1).
Теорема 2. Если m = 2 и остальные параметры системы (1) произвольны,
то в конусе K система (1) имеет единственное решение.
В третьей главе проведён систематический анализ влияния латеральных
взаимодействий, обратимости адсорбции на множественность ст.с., автоколебания в реакции, протекающей по параллельному адсорбционному механизму
A m  mZ  mAZ ; Bn  nZ  nBZ ; pAZ  qBZ  ( p  q)Z  A p Bq ,
(2)
где AZ, BZ – вещества на поверхности катализатора Z ; A m , Bn , A p Bq – вещества в
газовой фазе. Практически важны случаи: 1) m  n  p  q  1 ; 2) m  n  2;
p  q  1 ; 3) m  2; n  1; p  q  1. Соответствующие механизмы обозначим как
10
(I), (II) и (III). В случае идеального адсорбционного слоя только третий (механизм
Ленгмюра-Хиншельвуда (ЛХ)) демонстрирует множественность ст.с. Для описания механизма (2) используется следующая кинетическая модель:

dx / dt  mk1 PAm pm0  mk 1 pm0 exp m A   k3 p p 0 q 0 exp  p A  q B 
(3)





dy
/
dt

nk
P
p

nk
p
exp
n


k
p
exp
p


q

,

2 Bn n 0
2 n 0
B
3 p0q0
A
B

где x, y, – концентрации веществ AZ, BZ, соответственно; k1 , k 2 ( k 1 , k 2 ) – константы скоростей адсорбции (десорбции) газофазных веществ A m , B n , соответственно; k 3 – константа скорости реакции 3-й стадии в механизме (2); t – время;
PA m , PBn – парциальные давления веществ A m , B n ;  A ,  B – безразмерные хими-
ческие потенциалы адсорбированных веществ AZ и BZ, соответственно;
pn0 , pm0 , p p 0q 0 – вероятности найти пустыми n, m, p + q узлов решётки, имеющих
конфигурацию, необходимую для протекания соответствующего процесса.
Установлено, что механизм (II) качественно не меняет своего поведения независимо от свойств адсорбционного слоя, а механизмы (I) и (III) могут менять
своё поведение. Рассмотрим механизм (III). В качестве модели адсорбционного
слоя выбрана МРГ с двумя типами частиц с учётом латеральных взаимодействий
ближайших соседей. Термодинамический гамильтониан модели имеет вид:
H eff   AA  nA,i nA, j   AB  nA,i nB, j   BB  nB,i nB, j   A  nA,i   B  nB,i ,
nn
nn
nn 
i
i
где  AA ,  AB ,  BB – энергии латеральных взаимодействий ближайших соседей соответствующего типа; nA,i , nB,i – числа заполнения, равные нулю (узел пуст) или
единице (узел заполнен соответствующей частицей).
Построены фазовые диаграммы при Т=0 K на плоскости  A ,  B  для 27 наборов (a;b;с), в которых величины  AA  a,  AB  b,  BB  c принимали значения
10;-10;0, кДж/моль. Показано, что для квадратной и шестиугольной решёток существуют фазы: LG, LLA, LLВ, С(2х2)А, С(2х2)В, С(2х2)АВ, а для треугольной решётки: LG, LLA, LLВ, 3  3 R30 A , 3  3 R30 B , 3  3 R  30 A , 3  3 R  30 B ,


3  3 R30 A2B ,




3  3 R30 AB2 ,








3  3 R30 AB0 . Проведён анализ влияния лате-
ральных взаимодействий и обратимости адсорбции на множественность ст.с. и
автоколебания скорости реакции, протекающей по механизму ЛХ при Т=500 K.
Обозначив u  2k1 PA2 / k3 , v  k 2 PB / k3 , s  k 2 / k3 , w  2k 1 / k3 ,   k3t , из (3) имеем:
dx / d  p00 u  w exp 2 A   exp  A   B 

dy / d  v  s exp  B (1  x  y )  p00 exp  A   B .
(4)
Уравнения (4) по виду тождественны для трёх типов решётки, однако зависимости p00 ( A ,  B ) , x( A ,  B ) , y( A ,  B ) различны. Для их вычисления применялся МТМ. Показано существование не менее 10 внутренних (вн.) ст.с. для
квадратной и треугольной решёток и не менее 6 для шестиугольной решётки при
11
необратимой адсорбции. Показано, что сложность диаграммы кратности (границы области множественности ст.с.) связана со сложностью фазовой диаграммы.
Для всех типов решётки показано, что если фазовые диаграммы топологически
эквивалентны фазовой диаграмме идеального адсорбционного слоя, то диаграммы кратности топологически эквивалентны диаграмме кратности идеального адсорбционного слоя. Для квадратной решётки верно и обратное утверждение.
Проведён теоретический и численный анализ влияния параметров обратимости
на диаграммы кратности. Для произвольного набора латеральных взаимодействий доказан ряд утверждений об области множественности ст.с.
Теорема 3. Для произвольного набора латеральных взаимодействий и любого типа решётки область множественности ст.с. механизма ЛХ на плоскости
lg u, lg v : 1) при s  0, w  0 ограничена слева и неограниченна справа или не существует; 2) при s  0, w  0 неограниченна.
Показано, что обратимость стадий адсорбции сильно влияет на диаграммы
кратности, качественно меняя их вид, и совместное влияние обратимости адсорбции по обеим стадиям практически аддитивно (рис. 2). Для применения МТМ при
решении (4) переходим от (x,y) к (  A ,  B ). Для МРГ на полубесконечной решётке,
рассматриваемой в МТМ, якобиан перехода   0 ,   A ,  B .
Рис. 2. Диаграммы кратности для набора (10;10;10) при различных значениях параметров s,
w. Числами показано количество вн. ст.с.
При построении фазовых портретов и кинетических кривых использовался
алгоритм (Kaps and Rentrop), реализующий метод Розенброка с автоматическим
выбором шага, который позволяет получить решение с относительной
погрешностью не хуже, чем 10 5 . Проведён параметрический анализ системы (4)
для всех 27 наборов при Т=500 K. Для некоторых из них обнаружены
автоколебания скорости реакции, возникающие как результат бифуркации
Андронова-Хопфа. Установлено, что для квадратной и шестиугольной решётки
возникновение автоколебаний связано с существованием плотной фазы С(2x2)AB,
а для треугольной – фазы 3  3 R30 A2B . На рис. 3 приведены результаты для


набора (10;10;0) при s  0 на квадратной решётке. На рис. 3б видно, что
предельный цикл лежит внутри фазы С(2x2)AB. Как показывают результаты, это –
12
Рис. 3. Результаты для (10;-10;0). а) Бифуркационные диаграммы. б, в) Фазовый портрет, соответствующий точке Z (3,93;4,342 ) , в плоскости  A ,  B  , (x,y), соответственно; г, д) Скорость
реакции (сплошная линия), зависимости степени покрытия x, y по веществу А, В от времени,
соответственно
13
общее явление. На рис. 3г приведена зависимость скорости реакции от времени
для фазовой траектории, выходящей из неустойчивого фокуса, расположенного
внутри предельного цикла, видны автоколебания скорости реакции.
Исследованы диаграммы кратности для механизма (I) при всех возможных
сигнатурах латеральных взаимодействий ближайших соседей на квадратной
решётке. Показано, что область множественности ст.с. (в отличие от случая
идеального адсорбционного слоя) может возникать как результат совместного
воздействия обратимости стадии адсорбции и латеральных взаимодействий в
адсорбционном слое. Ни тот, ни другой фактор по отдельности не приводят к
такому результату. Показано, что в отличие от идеального случая возможно
существование, по крайней мере, трёх вн. ст.с.
В четвёртой главе изложены и исследованы модели, учитывающие наличие
нескольких типов АЦ в одной элементарной ячейке, а также модель
ступенчатой, трансляционно-инвариантной поверхности на квадратной
решётке.
А) Исследована модель системы H/Pd(100) с (f) четырёхкоординированными
и (b) мостиковыми центрами, учитывающая пространственные запреты и
латеральные взаимодействия 
между адсорбированными частицами,
расположенными в ближайших b-центрах. При   0 модель допускает точное
решение. Центры отличаются только разницей в энергии дна потенциальных ям
  E f  Eb . С учётом пространственных запретов элементарная ячейка (рис. 4)
может находиться в 5 состояниях и, соответственно, размер ТМ равен 5 M  5 M .
незаполненные и заполненные b-центры

незаполненные и заполненные f-центры
a
n
e
n

а - элементарная ячейка
Рис. 4. Показаны два типа АЦ, формирующих декорированную решётку.   маркирует 4 возможных запрещённых конфигураций, включающий активированный комплекс диффузии
Ряды решётки делятся на нечётные (n) и чётные (e). Элементарная ячейка
нечётного (чётного) ряда, состоит из 2 (1) АЦ и может находиться в 3 (2)
состояниях. Произведение матрицы «нечётный ряд – чётный ряд» размера
3M  2 M на матрицу «чётный ряд – нечётный ряд» размера 2 M  3M есть
симметричная матрица. С вычислительной точки зрения новизна используемого
подхода заключается в использовании неквадратных матриц. Построены
изотермы и коэффициент поверхностной диффузии как функция полной степени
покрытия при различных значениях , (рис. 5). Показано наличие особенности в
14
области существования упорядоченных структур. Максимум связан с
упорядоченной фазой C (2  2) на (b)-центрах, минимум (рис. 5с) – с областью
сосуществования фаз, а убывание – с запретом ряда конфигураций, включающих
активированный комплекс диффузии.
Рис. 5. Коэффициент поверхностной диффузии как функция полной степени покрытия.
Значения Δ, ε показаны на рисунке
Б) Исследована модель CO/Ni(100), учитывающая (b) мостиковые и (t)
монокоординированные центры, колебательные статистические суммы и
латеральные взаимодействия между адсорбированными частицами (рис. 6).
Показано, что в отличие от кластерного приближения МТМ описывает
формирование упорядоченных структур с большой элементарной ячейкой. На
рис. 7 показана упорядоченная структура при  tot  2 / 3 .
Рис. 6. Модель адсорбции CO/Ni(100). Элементарная ячейка “a” имеет 4 возможные состояния
Рис. 7. Упорядоченная структура С(3х2)
Так как при использовании МТМ фактически вычисляем заселенность для
кластера M   , то автоматически учитываются любые упорядоченные
структуры, соизмеримые с геометрической структурой решётки, расположенной
на бесконечном цилиндре периметра М. Элементы трансфер-матрицы имеют вид:
1
1

Tij  exp vi  v j  vij  , где
2
2

vi  
1
T
M
2
k 1
T
 (n1,i ,k n1,i ,k 1  n2,i ,k n2,i ,k 1  n3,i ,k n3,i ,k 1 ) 
15
M
 (n1,i ,k n2,i ,k 1  n2,i ,k n1,i ,k 1 ) 
k 1
(   2T ln f t ) M
(   T ln f b  E ) M
(5)
 n1,i ,k 
 (n2,i ,k  n3,i ,k ),
T
T
k 1
k 1
 M
 n2,i ,k n2, j ,k  n3,i ,k n3, j ,k )  2  (n1,i ,k n3, j ,k  n3,i ,k n1, j ,k  n2,i ,k n1, j ,k 1  n1, j ,k n3,i ,k 1 ),
T k 1

vij  
1
T
M
 (n1,i ,k n1, j ,k
k 1
где i,j соответствуют состояниям колец, n1,i ,k , n2,i ,k , n3,i ,k – числа заполнения для
элементарной ячейки, расположенной на k-м месте в кольце, находящимся в i-м
состоянии, k B  1 , f t  1/[1  exp( t / T )] , f b  1/[1  exp( b / T )] – колебательные
статистические суммы. Полученная в эксперименте относительная заселённость
мостиковых центров в области   0,67 (максимально измеренная степень покрытия) воспроизводится именно МТМ, который даёт более быстрый рост доли заполненных (b)-центров при 0,5    0,7 (рис. 8б) по сравнению с результатами,
полученными кластерным методом (рис. 8a). Наблюдаемые аррениусовские параметры десорбции представляются как
  ln k d ( )  ,
 E ( )  , k ( )  k (0)S   /(S (0) ) exp  / T  ,
Ed ( )  T 2 
  d ( )  k d ( ) exp  d
 d
d
T
 T 


(6)
где k d (0)  vd (0) exp  Ed (0) / T  – аррениусовский множитель; S ( ) – коэффициент
прилипания. С учётом S ( )  (1   )S (0) и (6) имеем k d ( )  k d (0)1    /   exp  / T  .
Рис. 8. Доля молекул СО, адсорбированных на мостиковых центрах (ΔЕ=1, ε1=4, ε2 = 1,
ккал/моль и hωt = hωb = 5, meV). (а) – кластерное приближение, (б) метод трансфер-матрицы
Из анализа рис. 9 видно удовлетворительное соответствие результатов моделирования в сравнении с экспериментальными данными для CO/Ni(100). Предэкспоненциальный фактор слабо меняется при   0,5 , затем быстро возрастает на несколько порядков, что связано с сильным отталкиванием ближайших соседей.
Построена и исследована модель трансляционно-инвариантной, ступенчатой поверхности на квадратной решётке. На рис. 10 схематически показана поверхность, соответствующая такой модели при L  3 , где L – ширина террасы
( L  2 ). Считаем, что узлы 1-го, (L+1)-го, (2L+1)-го и т.д. вдоль Y направления
являются выделенными. Рассмотрим мономолекулярную адсорбцию. Узлы каж16
дого 1-го ряда обладают потенциальной ямой глубины  (сильные АЦ). Энергия
частиц, адсорбированных на невыделенные узлы, равняется нулю (слабые АЦ).
Взаимодействуют только ближайшие частицы и энергия
взаимодействия  одинакова
для всех возможных расположений
взаимодействующей
пары. Построенная модель неоднородна, но трансляционноинвариантна. Большая статистическая сумма для решётки
из N рядов с циклическими
граничными условиями имеет
вид  N  Tr B1k , где
B1  T12 T22L2 T21
Рис. 9. Аррениусовские параметры десорбции для
CO/Ni(100). Заштрихованные области – экспериментальные результаты1, сплошные линии – расчётные
T12 , T22 , T21
–
(7)
трансфер-
матрицы,
соответствующие
переходу от выделенного ряда
к невыделенному ряду, между двумя невыделенными рядами и от невыделенного
ряда к выделенному, соответственно. Матрица B1 соответствует переходу между
двумя ближайшими выделенными рядами, T 22 , B1 – симметричны, T21  T12 T ,
i 2
B i  T22Li T21T12 T22
получаются
циклической
перестановкой (7). Собственный вектор матрицы
Bi ,
соответствующий
её
наибольшему
- адсорбированная частица
Рис. 10. Схематический вид
ступенчатой поверхности.
1 –выделенный ряд, 2,3 – невыделенные ряды
собственному значению, описывает распределение
вероятностей расположения адсорбированных
частиц в i-м ряду. Локальная степень покрытия i-го
 v li ,k v ri ,k sk
ряда вычисляется по формуле  i  k
, где
M
v li ,k ( v ri ,k )
– k-я компонента нормированного
левого (правого) собственного вектора, соответствующего наибольшему
собственному значению матрицы B i , s k – число занятых узлов в кольце,
находящемся в k-м состоянии; v li ,k , v ri ,k   1 .
При построении фазовых диаграмм использовался подход, основанный на
идеях конечно-размерного фазовых диаграмм использовался подход, основанный
1
Zhdanov V.P. Arrhenius parameters for rate processes on solid surfaces // Surface Sci. Rep., 1991. V.12, № 5. P.184 242.
17
на идеях конечно-размерного скейлинга, согласно которому точка фазового

 ( M )  ( M )
перехода определяется из равенств:
,  1 ( M )  ln max , где  (M ) –

M
M
2
безразмерная корреляционная длина, 2 – второе по модулю собственное
значение матрицы B1 . Построенные фазовые диаграммы изображены на рис. 11.
В целях верификации полученных результатов в нескольких точках
использовался ММК. Видно, что результаты, полученные МТМ и ММК,
практически совпадают. Показано, что при 2   /   3 и L  3 заполнение
выделенного ряда зависит от химпотенциала немонотонно, что может оказать
силь-ное влияние на кинетику химических реакций, протекающих на
ступенчатых поверхностях. Установлено отсутствие упорядоченных структур
при  /   4 .
Рис. 11. Фазовые диаграммы исследуемой модели
Исследована модель поверхности с одним видом адсорбированных частиц
на квадратной решётке со случайно распределёнными узлами двух типов.
Поверхность рассматриваем как квадратную решётку с 2 типами узлов,
распределённых случайным образом. Узлы 1-го и 2-го типа появляются с
вероятностью р и 1-р. Рассмотрим мономолекулярную адсорбцию частиц одного
вида. Для применения разработанного алгоритма так же, как и ранее, рассмотрим
цилиндр из M узлов в периметре. Следуя МТМ, данное построение позволяет
свести двумерную задачу к одномерной. В результате получаем одномерную
цепочку, состоящую из случайно распределённых узлов n  2 M типов. Каждый
узел может находиться в m  2 M состояниях. Построены термодесорбционные
спектры (ТДС) при М=4 и различных наборах латеральных взаимодействий.
Считая, что латеральные взаимодействия активированного комплекса десорбции
отсутствуют, константа скорости мономолекулярной десорбции вычислялась по
18
формуле (6). Уравнение мономолекулярной десорбции
d
1

dT


   
k d ( )     
 T   

  
 
  T
решалось методом Адамса-Башфорта. Показано значительное разнообразие форм
спектров, многие из которых характерны для реальных систем. Отметим наличие
высокотемпературных и низкотемпературных «плеч» (рис. 12), характерных для
экспериментальных TДС. Для исследования таких моделей ранее использовались
только стохастические методы.
Рис. 12. ТДС для модели
с центрами 2 типов, случайно распределёнными
по квадратной решётке.
Еа0 = 35 ккал/моль; β =
20 K/c; v0 =1013с-1; ε11 =
ε12 = ε22, ΔЕ=8 ккал/моль.
(а) р = 0,5, латеральные
взаимодействия
(ккал/моль) показаны на
рисунке, (б) ε11 = 4, значения вероятности р показаны на рисунке
В пятой главе описаны и исследованы модели многоцентровой адсорбции
молекул с возможностью различной ориентации в адсорбционном слое. Построена и исследована модель адсорбции гетероядерных двухатомных молекул (димеров) на квадратной решётке. Димер может адсорбироваться 2 способами, занимая: 1) один АЦ (трансверсальная адсорбция, рис. 13а,б); 2) два АЦ (планарная
адсорбция, рис. 13в). Возможны 2 различные ориентации молекул (рис. 13а,б). Независимо от расa)
положения адсорбированной молекулы теплота
адсорбции, приходящаяся на данный тип атома,
б)
в)
постоянна; энергия латеральных взаимодействий
между атомами различных молекул, находящихся
Рис. 13. Типы расположения
в ближайших узлах, не зависит от ориентации ададсорбированных молекул на
сорбированных молекул; атомы одного сорта отповерхности
талкиваются, атомы различных сортов притягиваются, энергия этих взаимодействий по модулю равны. Для трансверсально адсорбированных молекул учитывается взаимодействие как атомов, находящихся на
поверхности, так и атомов, находящихся над ней. Параметры модели: T ,  ,  –
разность теплот адсорбции, приходящихся на атомы различных сортов,  – энергия латеральных взаимодействий атомов одного сорта. Каждый узел решётки
может находиться в 11 состояниях (рис.14), 11-е состояние соответствует пусто19
му узлу. Построены изотермы и степени покрытия поверхности (рис.15). Отличительной особенностью является наличие двух максимумов в зависимости степени
покрытия для частиц, адсорбированных параллельно поверхности (рис.15в,г), что связано с формированием
структур С(2х2).
Построена и исследована модель адсорбции гомоядерных димеров. Учитывается бесконечно сильное отталкивание между ближайшими соседями. Параметры
Рис. 14. Возможные состояния узла решётки
модели: T ,  , h  h2  h1 , где h1 и h2 – теплота адсорбции,
соответственно, на 1 и 2 АЦ. При анализе основного состояния (Т=0 K) в случае
шестиугольной решётки при h > 0 было обнаружено явление, называемое в статистической физике «чёртовой лестницей»: переход от структуры, состоящей из
планарно адсорбированных
димеров, к структуре, состоящей из трансверсально
адсорбированных димеров,
происходит через бесконечную
последовательность
упорядоченных
структур,
образованных
адсорбированными димерами с различной ориентацией. Следует отметить, что в нашей
модели «чёртова лестница»
фазовых переходов есть
следствие наличия двух типов адсорбции, в то время
как в ранее изученных моделях «чёртова лестница»
Рис. 15. 1 - изотермы, 2 (3) - степени покрытия поверхнобыла следствием наличия
сти планарно (трансверсально) адсорбированными молекулами; 4 - раздельные степени покрытия трансверсально
нескольких типов конкуриадсорбированными молекулами для двух различных орирующих латеральных взаиентаций.   10 , кДж/моль
модействий. В табл.1 приведены соответствующие структуры, их элементарные ячейки, значения степени
покрытия и термодинамического потенциала в основном состоянии. АЦ, занятые
димерами, адсорбированными на 1 (на 2) АЦ показаны красными (черными и синими), свободные АЦ – белыми кругами. Упорядоченные фазы пронумерованы в
порядке их появления с увеличением химпотенциала. Для одного АЦ при переходе от одной структуры к другой имеем
20
n 
1 
 (6n  4)  4n1n  ,
(
6
n

4
)
h



 i 

( 4  2n ) 2 
i 1


где n – номер упорядоченной фазы. Аналогичное поведение адсорбционного слоя
можно наблюдать в реальных системах СОМ (адсорбция тримезиновой кислоты
на Au(111))2. Эта модель была исследована при T  0 . В случае шестиугольной
решётки рассматривался размер решётки, кратный размеру элементарной ячейки
упорядоченных фаз при n  1,4 и n   . Построены изотермы и степени покрытия
для ширин решёток M  6;8;10;12 (рис. 16). Видно, что образуются все упорядоченные фазы с элементарными ячейками, соизмеримыми размеру решётки. На
изотермах и функциях степени покрытия имеются ступеньки и плато. Отметим,
что степень покрытия у фазы №1 и №4 одинакова, а плотности структур различны. Итак, фаза № 4 может образовываться при ненулевых температурах, возможно образование и других промежуточных структур между фазами №1 и №∞.
Таблица 1
Упорядоченные структуры димеров на шестиугольной решётке
№
n
Упорядоченная
фаза
Данные о
структуре
№
n
 6  6
1

5
5
h 
18
18
2
1
5
   h 
4
16
4
3

11
17
h 
50
50

7
13
h 
36
36
  0,55(5)
…
 20  20
8
  0,5625
10 10 
Данные о
структуре
12 12 
  0,55(5)
8  8
Упорядоченная
фаза

13
41
h

100
100
  0,54
…
 2  2

  0,56
1
 
2
  0,5
В случае квадратной решётки были получены изотермы (рис. 17), степени
покрытия (рис. 18), скрытые теплоты адсорбции, внутренняя энергия, энтропия.
2
Ye Y.C., Sun W., Wang Y.F., Shao X., Xu X.G., Cheng F., Li J.L., Wu K. A. Unified Model: Self-Assembly of
Trimesic Acid on Gold // J. Phys. Chem. C. 2007. V.111. Р.10138-10141.
21
Рис. 16. Парциальные изотермы и функции степени покрытия от химпотенциала при
Т = 200 K для разных ширин решёток М в случае шестиугольной решётки
Рис. 17. Изотермы модели димеров в случае
квадратной решётки
Рис. 18. Степень покрытия поверхности
в случае квадратной решётки
При T<600 K на изотермах наблюдаются плато и ступеньки. При T>600 K первое
плато исчезает и в системе наблюдается один фазовый переход II рода в плотную
упорядоченную структуру. Легко видеть, что степень покрытия поверхности и
изотерма ведут себя по-разному. Упорядоченная фаза, соответствующая 1-му (2му) горизонтальному участку, состоит из димеров, адсорбированных на 2 (1) АЦ,
22
и имеет структуру с(4х2) (с(2х2)) (рис. 19а). В данном случае происходит фазовый переход, заключающийся в изменении типа адсорбции димеров и проявляющийся на изотермах в виде вертикального участка в окрестности μ=h. Фаза, состоящая из молекул, ориентированных двумя способами одновременно, не образуется в данной модели.
c ( 4  2)
c ( 2  2)
c (5  3 )
(а)
c(4  3)
p( 3  3)
(б)
Рис. 19. Поверхностные структуры димеров при T=0 K в случае: (а) квадратной, (б) треугольной решётки. Рамками выделены элементарные ячейки
В случае треугольной решётки были получены изотермы (рис. 20), степени
покрытия (рис. 21), внутренняя энергия и парциальные изотермы. С понижением
температуры на  (  ) и  (  ) появляются ступеньки и плато при ρ = 0,2; 0,25;
0,33(3) и   0,4 ; 0,375; 0,33(3), соответственно. Плато возникают в области существования упорядоченных фаз c(5  3) , c(4  3) и p( 3  3) (рис. 19б).
Структура c(5  3) образуется из решёточного газа в результате фазового перехода II рода и переходит в фазу c(4  3) посредством фазового перехода I рода.
Рис. 20. Изотермы модели димеров в случае
треугольной решётки
Рис. 21. Степень покрытия поверхности в
случае треугольной решётки
Показано, что зависимость  (  ) может быть немонотонной, т.е. чем больше давление, тем больше пустых мест на поверхности. Это неожиданное явление объясняется совместным воздействием двух факторов: два различных способа адсорбции (один из которых многоцентровый) и бесконечно сильное отталкивание ближайших соседей. Показано, что изотерма адсорбции и  (  ) непропорциональны.
23
В шестой главе рассматривается обобщённая решёточная модель адсорбции
СОМ на поверхность твёрдого тела. Полагается, что частица способна адсорбироваться k способами, отличающихся количеством занятых узлов, их конфигурацией и теплотами адсорбции q1 ,, qk . Решётка может быть как однородной, так
и неоднородной. Учитываются различные типы латеральных взаимодействий
между адсорбированными молекулами. Исследован случай k  2, m1  1, m2  4 на
квадратной решётке. Молекулы адсорбируются двумя способами, занимая 1 и 4
АЦ. Запрещается ближайшее соседство между двумя адсорбированными молекулами. Параметры модели: Δ – разница между теплотами адсорбции молекул, адсорбированных на 4 и 1 АЦ; μ; T. Каждая ячейка может находиться в 6 состояниях (обозначены числами 0,1,2,3,4,5 на рис. 22). Возможны фазы c4 44 ,
c3  341 , c2 2 (рис. 23). Вычисления проводились при М=12 (наименьшее зна-
чение ширины полосы, совместимое со всеми упорядоченными структурами).
Рис. 22. Допустимые конфигура-
Рис. 23. Возможные упорядоченные фазы
ции адсорбированных частиц
Из анализа рис. 24
следует, что на изотермах
и степенях покрытия
имеются три плато, соответствующие
упорядоченным структурам. По
мере роста  вначале из
решёточного газа возникает фаза c4 44 с плот-
Рис. 24. θ1 (θ4) – парциальные покрытия поверхности молекулами, адсорбированными на 1 (4) АЦ; θ – полная степень
покрытия поверхности; ρ – изотерма адсорбции
24
ностью   0,125 и   0,5 .
Здесь имеет место непрерывный фазовый переход. Далее фаза c4 44
превращается в фазу
c3  341 с   0,22(2) и
  0,55(5) , посредством фазового перехода I рода. Аналогично, фаза c3  341
превращается в фазу c2 2 посредством фазового перехода I рода. Заметим, что
при возрастании  меняется тип адсорбции. Установлено, что немонотонное изменение степени покрытия (рис. 24) связано с существованием фазы c3  341 .
Разработана и изучена модель адсорбции СОМ на ступенчатую поверхность. За основу взята модель ступенчатой поверхности, введённая выше (рис.
10). Молекула может адсорбироваться на поверхность 2 способами: на 1 и на 4
L=1
L=2
L=3
L=4
L=5
L=6
Рис. 25. Изотермы (кривые 1, 2 – для случая 1, 2, соответственно) и степени покрытия (кривые
3, 4 – для случая 1, 2)
25
АЦ (рис. 22), все АЦ энергетически эквивалентны. Рассмотрены два случая, когда молекулы, адсорбированные на соседние АЦ, расположенные на разных террасах: а) не взаимодействуют между собой; б) могут взаимодействовать между
собой в зависимости от способа адсорбции. На рис. 25 приведены изотермы и
степени покрытия для различных величин ширины террас L (от 1 до 6 АЦ) и
Δ/RT=9, при котором образуются упорядоченные структуры. При L  1 степень
покрытия  , как и в случае однородной поверхности, является немонотонной.
Величина θ (разность между максимальным и конечным значением степени покрытия) зависит от морфологии поверхности, а в данной модели – от ширины
террасы L и имеет максимальное значение 0,166(6) при L  2 . Такое поведение 
определяется последовательностью образования различных структур (рис. 25).
Показано, что, когда размер адсорбированной частицы соизмерим с шириной
террасы, возможно образование структур A 2 , A 3 , B 4 , A 5 , B5 , не характерных
для случая однородной поверхности. Вычисленные значения внутренней энергии
U и других термодинамических характеристик как функции плотности  для
значений L  1,..., 6 подтверждают существование указанных выше структур. Показано, что при L  1 реализуется один режим адсорбции, при L  2 – два режима
адсорбции и при L  3 – три режима адсорбции. Заметим, что результаты, полученные для двух моделей, качественно не отличаются друг от друга.
Разработанная обобщённая модель адсорбции СОМ применена при моделировании реальной адсорбционной системы. Адсорбционный монослой 1,4циклогексадиена (ЦГД) на поверхности Si(001) - 2x1 был выбран как модельная
система для исследования качественных особенностей поведения ненасыщенных
циклических углеводородов на поверхности Si(001). В данной модели учитывается экспериментально установленный факт3, что молекула 1,4-ЦГД адсорбируется
3-мя различными способами, занимая 1 АЦ (πкомплекс), 2 АЦ (ди-σкомплекс) и 4 АЦ (двойной-ди-σ-комплекс)
(рис. 26). Образование
каждого из комплексов
сопровождается определённым тепловым эфа)
б)
в)
фектом: 1 ,  2 – для диРис. 26. Решёточная модель адсорбции 1,4-ЦГД на Si(001) - 2x1
σ-комплекса, двойногоди-σ-комплекса, соответственно. Теплоту адсорбции 1,4-ЦГД с формированием
3
Kato H.S., Wakatsuchi M., Kawai M., Yoshinobu J. Different adsorbed states of 1,4-cyclohexadiene on Si(001) controlled by substrate temperature // J. Phys. Chem. 2007. V.111. P. 2557-2564.
26
π-комплекса учитывали неявно в химпотенциале. Согласно литературным данным, полагали 1   ,  2  2 . Вычисленные характеристики приведены на рис.
27, 28. Показано, что разработанная модель качественно воспроизводит все особенности поведения реальной системы 1,4-ЦГД на Si(001) - 2x1 . Так же, как и выше, существуют области θ, в которых с ростом  уменьшается степень покрытия.
Возникают 2 упорядоченные фазы с симметрией типа «шахматная доска», одна
из них состоит из π-комплексов, а вторая – из двойных-ди-σ-комплексов (рис. 29).
В промежуточной области значений  существует фаза, состоящая из ди-σкомплексов и двойных-ди-σ-комплексов.
Рис. 27. 2 - изотерма  адсорбции, 1,3,4,5 - полная  , парциальные 1 , 2 , 4 степени покрытия
поверхности
а)
Рис. 28. Энтропия адсорбционного слоя
1,4-ЦГД на Si(001) - 2x1
б)
в)
г)
д)
Рис. 29. Структура упорядоченных фаз 1,4-ЦГД на Si(001) - 2x1 . а) фаза из двойных-ди-σкомплексов и фаза из ди-σ-комплексов, не обнаруженная в рамках построенной модели; б)
смешанная фаза; г) СТМ-изображение4 1,4–ЦГД на Si(001) - 2x1 .
Фаза c(2×2) из π-комплексов: в) в рамках исследуемой модели, д) изображение 1,4–ЦГД на
Si(001) - 2x1 , полученное методом дифракции медленных электронов3
4
Akagi K., Tsuneyuki S., Yamashita Y., Hamaguchi K., Yoshinobu J. Structural and chemical property of unsaturated
cyclic-hydrocarbon molecules regularly chemisorbed on Si(001) surface // Appl. Surf. Sci. 2004. V. 234. P.162.
27
Разработана и исследована решёточная модель, учитывающая направленные
взаимодействия в адсорбционном монослое СОМ. Для построения модели направленных межмолекулярных взаимодействий в качестве основы был выбран
адсорбционный слой тримезиновой кислоты (ТМК) на грани (111) монокристалла
с г.ц.к. решёткой. Молекула ТМК – это стандартная молекула, используемая для
исследования свойств и условий формирования двумерных пористых структур5,6,7. Предполагается, что молекула ТМК занимает один узел треугольной решётки. Учитываются 2 возможные ориентации (рис. 30а,б) молекулы, определяющиеся положением карбоксильных групп, и 3 типа направленных взаимодействий (рис. 30в), отличающиеся взаимной ориентацией молекул. Каждому типу
а)
в)
б)
Рис. 30. а, б) Возможные
ориентации молекулы ТМК
в рамках МРГ; в) 3 типа направленных взаимодействий с энергиями w1 , w2 , w3
соответствует определённая энергия w1 , w2 , w3 . Модель может описывать адсорбцию любых молекул, симметрия которых совпадает с симметрией ТМК4,5,7. Вычислены изотермы (рис. 31), парциальные изотермы, внутренняя энергия и дифференциальная теплота адсорбции при различных значениях величин взаимодействий w1 , w2 , w3 . Построенная модель воспроизводит качественные особенности
а)
б)
Рис. 31. Изотермы модели адсорбции ТМК при w1 / RT  7,5 : а) w3 / RT  0 , для различных
значений w2 / w1 ; б) w2 / RT  3 , для различных значений w3 / w1
5
Liang H., He Y., Ye Y., Xu X., Cheng F., Sun W. et al. Two-dimensional molecular porous networks constructed by surface assembling // Coord. Chem. Rev. 2009. V.253. P.2959.
6
Cicoira F., Santato C., Rosei F. Two-Dimensional Nanotemplates as surface cues for the controlled assembly of organic
molecules // Top. Curr. Chem. 2008. V.285. P.203.
7
Kannappan K., Werblowsky T.L., Rim K.T., Berne B.J., Flynn G.W. An experimental and theoretical study of the formation of nanostructures of self-assembled cyanuric acid through hydrogen bond networks on graphite // J. Phys. Chem.
B. 2007. V. 111. P. 6634.
28
адсорбции ТМК и других молекул, симметрия которых совпадает с симметрией
ТМК. Для этих систем наиболее общими и характерными структурами являются
упорядоченная фаза сотового типа и плотноупакованная фаза (рис. 32), которые
возникают вследствие фазовых переходов I (рис. 31). Изотерма Ленгмюра является хорошей аппроксимацией, описывающей образование системы «гость-хозяин»
после образования сотовой структуры (рис. 32). В реальных системах w3 / w1  0,5 .
Сотовая структура
ТМК на Аg(111)
Плотноупакованная структура (модель)
БТБК на Ag(111)
ТКМБ на графите
БТБК на Ag(111)
Структура типа «гость-хозяин» (модель)
ЦК на графите
ЦК на графите
ТМК на Ag(111)
Рис. 32. Упорядоченные структуры, образующиеся при адсорбции ТМК и молекул с такой же
симметрией. Слева – изображения, полученные в рамках разработанной модели, остальные получены разными авторами с помощью СТМ4,5,7,8
ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ
Предложен новый комплексный подход к изучению решеточных моделей адсорбционных слоёв, а также моделей магнетиков изинговского типа, основанный
на развитии метода трансфер-матрицы. В рамках разработанного подхода:
1. Построены математические модели, позволяющие описывать и прогнозировать особенности поведения реальных адсорбционных систем с несколькими
типами адсорбционных центров и различными способами адсорбции.
2. Разработан новый математический метод для изучения решёточных моделей без трансляционной инвариантности, основанный на развитии основных под8
Su G.-J., Zhang H.-M., Wan L.-J., Bai C.-L., Wandlowski T. Potential-Induced Phase Transition of Trimesic Acid
Adlayer on Au(111) // J. Phys. Chem. B. 2004. V.108. P.1931-1937.
29
ходов метода трансфер-матрицы. Доказаны теоремы, обосновывающие вычислительный алгоритм. Ранее детерминистские методы для изучения таких систем не
использовались.
3. Модифицирован метод трансфер-матрицы для неоднородных, трансляционно-инвариантных систем, включая использование неквадратных матриц. Для
изучения модели адсорбции на ступенчатой поверхности впервые был использован метод трансфер-матрицы.
4. Показано, что латеральные взаимодействия в адсорбционном слое сильно
влияют на наблюдаемую кинетику гетерогенно-каталитических процессов, в том
числе приводят к её качественным изменениям. В частности, для механизма Ленгмюра-Хиншельвуда область множественности стационарных состояний претерпевает существенные изменения, и сложность диаграммы кратности определяется сложностью фазовой диаграммы адсорбционного слоя. Для некоторых наборов
энергий латеральных взаимодействий обнаружены автоколебания скорости реакции, возникающие как результат бифуркации Андронова-Хопфа.
5. Разработана и исследована модель многоцентровой адсорбции, учитывающая различные способы адсорбции. Обнаружены два новых явления: а) немонотонное изменение степени покрытия поверхности от химического потенциала
(давления), б) новый тип «чертовой лестницы» фазовых переходов.
6. Построены и с помощью метода трансфер-матрицы исследованы модели,
воспроизводящие качественные особенности поведения реальных адсорбционных систем: CO/Ni(100), 1,4-циклогексадиена на поверхности Si(001)-2×1, тримезиновой кислоты на Au(111).
ОСНОВНЫЕ ПУБЛИКАЦИИ АВТОРА ПО ТЕМЕ ДИССЕРТАЦИИ
Монографии
1. Мышлявцев А.В., Мышлявцева М.Д. Вычислительные аспекты метода трансферматрицы. - Кызыл: ТувИКОПР СО РАН, 2000. - 101с.
2. Мышлявцев А.В., Мышлявцева М.Д. Механизм Ленгмюра-Хиншельвуда: неидеальность
адсорбционного слоя и автоколебания. - Омск: Изд-во ОмГТУ, 2012. - 120с.
Обзоры
3. Fefelov V.F., Gorbunov V.A., Myshlyavtsev A.V., Myshlyavtseva M.D. Chapter 15. Statistical Thermodynamics of Lattice Gas Models of Multisite Adsorption. // Thermodynamics - Fundamentals and Its Application in Science, Ricardo Morales-Rodriguez (Ed.), ISBN: 978-953-51-0779-8,
InTech, 2012, Р. 389-410.
Статьи, опубликованные в изданиях, входящих в базы Web of Science и Scopus
4. Myshlyavtsev A.V., Dongak M.D. (Myshlyavtseva). Statistics of adsorption on top and bridge
sites of a square lattice: transfer matrix approach // J. Stat. Phys. 1997. V.87, № 3/4. P. 593-607.
5. Myshlyavtsev A.V., Myshlyavtseva M.D. Apparent Arrhenius Parameters for Desorption from
a Square Lattice with Top and Bridge Sites: Transfer-Matrix Approach // Phys. Low-Dim. Struct.
1998. V.9/10. P.55-64.
6. Myshlyavtsev A.V., Myshlyavtseva M.D. Thermodynamics for Lattice Models with Randomly
Distributed Quenched Impurities: Exact Solution // Phys. Low-Dim.Struct. 2000. V. 9/10. Р. 127-135.
7. Myshlyavtsev A.V., Myshlyavtseva M.D. Modeling of adsorption and phase diagrams for
stepped surfaces: Transfer matrix approach // Applied surf. science. 2007.V.253, № 13. P. 5591-5595.
30
8. Мышлявцев А.В., Мышлявцева М.Д. Латеральные взаимодействия в адсорбционном
слое и критические явления в реакции, протекающей по механизму Ленгмюра-Хиншельвуда //
Кинетика и катализ, 2007. Т. 48, № 4. С. 576-585.
9. Fefelov V.F., Gorbunov V.A., Myshlyavtsev A.V., Myshlyavtseva M.D. The simplest selfassembled monolayer model with different orientations of complex organic molecules. Monte Carlo
and transfer-matrix techiques // Chemical Engineering Journal, 2009. V.154. P. 107-114.
10. Fefelov V.F., Gorbunov V.A., Myshlyavtsev A.V., Myshlyavtseva M.D., Evsejeva S. The
simplest model of adsorption of molecules with different orientations in adlayer on the stepped surface
// Applied surface science. 2010. V. 256, Issue 17. P. 5298-5304.
11. Fefelov V.F., Gorbunov V.A., Myshlyavtsev A.V., Myshlyavtseva M.D. Model of
homonuclear dimer adsorption in term of two possible molecule orientations with respect to surface:
Square lattice // Phys. Rev. E. 2010. V.82. P. 041602-1 – 041602-5.
12. Горбунов В.А., Мышлявцев А.В., Мышлявцева М.Д., Фефелов В.Ф. Моделирование адсорбции ненасыщенных циклических углеводоров на реконструированной поверхности кремния (001)-2х1 методами Монте-Карло и трансфер-матрицы // Журнал физической химии, 2011.
Т. 85, № 1. С. 99-106.
13. Fefelov V.F., Gorbunov V.A., Myshlyavtsev A.V., Myshlyavtseva M.D., Akimenko S.S.
Devil’s staircase behavior of a dimer adsorption model. // Adsorption. V. 19. Numbers: 2-4. 2013.
DOI 10.1007/s10450-013-9471-1. P. 495-499.
14. Gorbunov V.A., Akimenko S.S., Myshlyavtsev A.V., Fefelov V.F., Myshlyavtseva M.D. Adsorption of triangular-shaped molecules with directional nearest-neighbor interactions on a triangular
lattice. // Adsorption. V. 19. Numbers: 2-4. 2013. DOI 10.007/s10450-013-9480-0. P. 571-580.
15. Горбунов В.А., Мышлявцев А.В., Мышлявцева М.Д., Фефелов В.Ф. О немонотонности
функции степени покрытия в моделях многоцентровой адсорбции с различной ориентацией по
отношению к поверхности твердого тела. // Физикохимия поверхности и защита материалов. –
2013. Т. 49, № 4. С. 352-360.
Статьи, опубликованные в изданиях, входящих в базу Scopus
16. Мышлявцев А.В., Мышлявцева М.Д., Романовский Р.К. Применение метода трансферматрицы к решеточным системам без трансляционной инвариантности // Сибирский журнал
индустриальной математики, 2006. Т. IX, № 1(25). С. 106-115.
Статьи, опубликованные в изданиях перечня ВАК, кроме Web of Science и Scopus
17. Мышлявцев А.В., Мышлявцева М.Д. Неидеальность адсорбционного слоя и автоколебания в механизме Ленгмюра-Хиншельвуда. Необратимая адсорбция // Омский научный вестник, 2006. № 1(34). С. 57-60.
18. Мышлявцев А.В., Мышлявцева М.Д. Изотермы адсорбции гетероядерных двухатомных
полярных молекул на однородных поверхностях: метод трансфер-матрицы // Омский научный
вестник, 2006. № 2(35). С. 64-68.
19. Мышлявцев А.В., Мышлявцева М.Д. Множественность стационарных состояний и автоколебания в механизме Ленгмюра-Хиншельвуда в случае шестиугольной решетки // Омский
научный вестник, 2006. № 7(42). С. 62–68.
20. Мышлявцев А.В., Мышлявцева М.Д. Множественность стационарных состояний и автоколебания в механизме Ленгмюра-Хиншельвуда в случае треугольной решетки. Необратимая адсорбция // Омский научный вестник, 2006. № 6(41). С.60-67.
21. Мышлявцев А.В., Мышлявцева М.Д. Влияние латеральных взаимодействий на область
множественности параллельного адсорбционного механизма для случая мономолекулярной
адсорбции по обоим веществам // Изв. вузов. Химия и химтехнология,2007. Т.50, №5. С. 27-31.
22. Мышлявцев А.В., Мышлявцева М.Д. Сравнительный анализ влияния типа решетки на
область множественности в механизме Ленгмюра-Хиншельвуда в условиях неидеальности адсорбционного слоя // Известия вузов. Химия и химтехнология, 2007. Т. 50, № 11. С. 104-109.
23. Мышлявцев А.В., Мышлявцева М.Д. Влияние температуры и ширины полосы, используемой в методе трансфер-матрицы, на критические явления в реакции, протекающей по механизму Ленгмюра-Хиншельвуда // Омский научный вестник, 2012. № 2 (110). С. 25-28.
31
24. Мышлявцев А.В., Мышлявцева М.Д. Теоретический анализ влияния обратимости мономолекулярной и бимолекулярной стадий адсорбции на диаграммы кратности механизма Ленгмюра-Хиншельвуда // Доклады Академии наук высшей школы РФ, 2012. № 1(18). С. 6-18.
25. Мышлявцев А.В., Мышлявцева М.Д. Анализ влияния типа решетки на автоколебания
скорости реакции, протекающей по механизму Ленгмюра-Хиншельвуда. // Вестник ИжГТУ.
Серия «Математика», 2012. Вып. 3. С. 150-154.
Статьи, опубликованные в других рецензируемых изданиях
26. Myshlyavtsev A.V., Dongak M.D. (Myshlyavtseva). Isotherms and chemical diffusion coefficient for the simplest model with two types of adsorption centres into one elementary cell // Phys.
Low-Dim. Struct. 1996. V.4/5. P.65-73.
27. Мышлявцев А.В., Мышлявцева М.Д. Диаграммы кратности для параллельного адсорбционного механизма в условиях неидеальности адсорбционного слоя. Необратимая адсорбция
// Доклады Академии Наук Высшей школы РФ, 2005. № 2(5). С. 56-64.
Свидетельства о государственной регистрации программ для ЭВМ
1. Мышлявцев А.В., Мышлявцева М.Д. Построение фазовых диаграмм и степеней покрытия для модели ступенчатой поверхности, рег. № 2013613124 (26.03.2013) // Программы для
ЭВМ, базы данных, топологии интегральных микросхем. RU ОБПБТ. № 2, 20.06.2013.
2. Мышлявцев А.В., Мышлявцева М.Д. Наблюдаемые аррениусовские параметры десорбции для модели с несколькими типами активных центров (модель грани (100) г.ц.к.), рег.
№ 2013613405 (04.04.2013) // Программы для ЭВМ, базы данных, топологии интегральных
микросхем. RU ОБПБТ. № 2, 20.06.2013.
3. Мышлявцев А.В., Мышлявцева М.Д. Термодесорбционные спектры систем без трансляционной инвариантности, рег. № 2013613506 (09.04.2013) // Программы для ЭВМ, базы
данных, топологии интегральных микросхем. RU ОБПБТ. № 2, 20.06.2013.
4. Мышлявцев А.В., Мышлявцева М.Д. Вычисление термодинамических характеристик
модели пирамидальных молекул, адсорбирующихся на основание или на вершину, с использованием метода трансфер-матрицы при учете бесконечных отталкивательных взаимодействий ближайших соседей, рег. № 2013613573 (10.04.2013) // Программы для ЭВМ, базы данных, топологии интегральных микросхем. RU ОБПБТ. № 2, 20.06.2013.
5. Мышлявцев А.В., Мышлявцева М.Д. Построение границ областей устойчивости стационарных состояний для механизма Ленгмюра-Хиншельвуда, рег.№2013613671 (12.04.2013)
// Прогр. для ЭВМ, базы данных, топол. интегр. микросхем. RU ОБПБТ. № 2, 20.06.2013.
Печатается в авторской редакции
Подписано в печать _______. Формат 60х841/16. Бумага офсетная.
Отпечатано на дупликаторе. Усл.печ.л. 2. Уч.-изд.л. 2.
Тираж 100 экз. Заказ ____
______________________________________________________
Издательство ОмГТУ. 644050, г. Омск, пр. Мира, 11; т. 23-02-12
Типография ОмГТУ
32
Документ
Категория
Без категории
Просмотров
7
Размер файла
3 595 Кб
Теги
154234
1/--страниц
Пожаловаться на содержимое документа