8 / 12
М.Б. Логинова, К.В. Марамыгин, А.В. Пономарев, С.В. Русинов, И.О. Сакович
8
рисунок). Таким образом, первым (расчетным) этапом разрабатываемой
методологии является определение уровней энергетических барьеров
Энергетическая диаграмма диссоциации адсорбируемого
элемента на поверхности с учетом переходного (
TS
)
состояния (
IS
,
FS
— начальное и конечное состояния) [10]
бар
i
E
возможных химических реакций в соответствии с иссле-
дуемыми процессами взаимодействия в рассматриваемых гетерогенных
системах. После вычисления значений
бар
i
E
выполняют интер-
вальную оценку необходимых термодинамических условий, обеспе-
чивающих возможность протекания реакций в рассматриваемой
системе, в частности, необходимые и избыточные значения темпе-
ратуры
н и
,
T T
и давления
н и
,
.
P P
Затем формулируют начальные
и граничные условия и задают термодинамические режимы (пиковое
нагружение, адиабатическое сжатие, релаксационный отжиг и др.)
для проведения моделирования.
Второй этап связан непосредственно с моделированием кинематики
взаимодействия элементов гетерогенной системы в рамках сформули-
рованных на первом этапе необходимых термодинамических режимов и
в интервалах установленных значений температуры
н и
...
T Т
и давле-
ния
н и
...
.
P P
Силы межатомарного взаимодействия расчитывают на
основе градиента полной энергии, потенциальную компоненту которой
сначала определяют в МД-приближении, что позволяет решать
кинематические уравнения движения для большого числа взаимо-
действующих элементов на основе алгоритма Верлета [11, 12]