CC BY
16
ТРУДЫ БГТУ. 2014. № 6. Физико-математические науки и информатика. С. 29-31
29
УДК 531.19
Я. Г. Грода, кандидат физико-математических наук, доцент, заведующий кафедрой (БГТУ)
ОПРЕДЕЛЕНИЕ КОЭФФИЦИЕНТА ДИФФУЗИИ РЕШЕТОЧНОГО ФЛЮИДА С ВЗАИМОДЕЙСТВИЕМ В СЕДЛОВОЙ ТОЧКЕ НА ОБЪЕМНОЦЕНТРИРОВАННОЙ КУБИЧЕСКОЙ РЕШЕТКЕ В РАМКАХ СУПЕРПОЗИЦИОННОГО ПРИБЛИЖЕНИЯ
Рассмотрен процесс диффузии решеточного флюида с дополнительным взаимодействием в седловой точке на объемноцентрированной кубической решетке. Предложена методика получения аналитического выражения для оценки кинетического коэффициента диффузии решеточного флюида. Определен кинетический коэффициент диффузии решеточного флюида с притяжением (отталкиванием) ближайших соседей и равным ему взаимодействием в седловой точке. Результаты аналитических расчетов сопоставлены с данными компьютерного моделирования диффузионного процесса в исследуемой системе по методу Монте-Карло.
The diffusion process of the lattice fluid with additional interaction in the saddle point on a body-centered cubic lattice is considered. The technique of obtaining an analytic expression for the jump diffusion coefficient of the lattice fluid is proposed. The jump diffusion coefficient of the system with nearest neighbor attractive (repulsive) interaction and interaction in the saddle point of the same energy is calculated. The analytical results are compared with Monte Carlo simulation data.
Введение. Как было показано ранее, для решеточных флюидов на плоских квадратной [1] и треугольной [2] и простой кубической [3] решетках при исследовании процесса диффузии решеточного флюида может быть рассмотрена модель, в которой наряду с взаимодействием частиц, занимающих ближайшие соседние узлы, может быть учтено взаимодействие частиц, находящихся в так называемой седловой точке, которая совпадает с вершиной межузлового барьера, с их ближайшими соседями. Учет этого дополнительного взаимодействия вызывает изменение эффективного межузлового барьера и, очевидно, будет влиять на диффузионные свойства системы [4].
Для оценки коэффициентов диффузии решеточных флюидов на всех указанных выше решетках в работе [1] было предложено так называемое суперпозиционное приближение, суть которого состоит в непосредственном учете лишь парных (бинарных) корреляций в заполнении ближайших соседних узлов. Все корреляционные функции более высоких порядков определяются посредством произведения соответствующих парных.
В настоящей работе приводятся результаты применения указанного подхода к случаю диффузии решеточного флюида на объемноцентрированной кубической (ОЦК) решетке.
Модель. Модель представляет собой систему из п частиц, расположенных в узлах регулярной ОЦК-решетки, содержащей N узлов. Каждый узел может либо быть занят частицей, либо быть вакантным. Состояние узла 1 определяется числом заполнения п = 1 или п = 0 в зависимости от того, занят узел частицей или вакантен. Заполнение узла более чем одной частицей запрещено.
Находящаяся в узле О частица может взаимодействовать с энергией J с частицами, занимающими ближайшие соседние узлы, т. е. узлы
а, Ь, с, d, е, /, g и к (рис. 1). При ее последующем переходе в узел а при прохождении седловой точки Е которая расположена в середине связи Оа она также взаимодействует с энергией JЕ с узлами Ь, d, е, а, в и у.
gO
Рис. 1. Простая кубическая решетка:
узлы а, Ь, с, d, е, /, g и к -ближайшие соседи узла О; узлы Ь, d, е, а, в и у -ближайшие соседи седловой точки
Таким образом, преодолеваемый частицей активационный барьер
Еа = Е (Е) - Е (0), (1)
где Е(Е) - энергия частицы в седловой точке,
Е(Е) = Е0 + Е (пЬ + пй + Пе + Па + Пр + П1X
здесь Е0 - исходная высота межузельного барьера; Е(0) - начальная энергия частицы, вычисляемая по формуле
Е (0) = J (пЬ + пс + па + пе + пг + п8 + пк).
Это позволяет представить активационный барьер в следующем виде:
Ea = En
-J (nc + nf + ng + nh) +
с
30
Я Г. Грода
+J2 (ив+ Пв+ «) -А(« + + «), (2)
где А = J - Jz.
Кинетический коэффициент диффузии решеточного флюида в суперпозиционном приближении. В дальнейшем, как и в случае решеточного флюида на простой кубической решетке [3], в работе будет рассматриваться система с J = Js и А = 0. Для рассматриваемой модели в рамках общей теории диффузионных процессов в решеточных системах [5] и с учетом предложенного ранее суперпозиционного приближения [1] может быть предложено следующее выражение для оценки кинетического коэффициента диффузии
М = (1 - ^ )(1 + ocg )4 + М)
+3^(1 - cg2 )(1 + ^ )3 (1 + ^2) + +3с2^2 (1 - cg3 )(1 + ccg2 )2 (1 + ccg )2 + +с3^3 (1 - cg4 )(1 + ^2 )3 (1 + ccg), (3)
где
МО = %
a = exp(pJ) -1; ехр(-^) -1;
с - равновесное значение концентрации частиц; g - парная корреляционная функция двух ближайших соседних узлов; г - число ближайших соседей на решетке рассматриваемого типа; а - расстояние между узлами решетки (длина прыжка частицы); ё - размерность пространства; V - частота, имеющая порядок частоты колебаний частицы вблизи узла решетки и определяющая временную шкалу диффузионных процессов; в = 1 / квТ - обратная температура; кв - постоянная Больцмана; Т - температура.
Диаграммное приближение. Входящая в соотношение (3) парная корреляционная функция двух ближайших соседних узлов g может быть найдена, например, в рамках диаграммного приближения [6, 7], согласно которому
_ (1 - с) - 2тдЫтс g _ 7(7 - с) ° гс2(27 - 1)(с-1) ,
где
(г - 2)^ - г;
4дтйт ;
7 _ 0,5(1 + 71 + 4с(1 -с)(Ж-1));
ё -1. ё _ с(1 - с)(Ж -1),
с №+1' 72 '
Ж _ ехр(-в/); _ exp(-PcJ);
т и д - число вершин в простейшей кольцевой диаграмме и ее вес соответственно. Коэффициент X может быть определен из условия равенства критической температуры системы ее точному значению либо его наилучшей оценке. В рассматриваемом примере решеточного флюида на плоской квадратной решетке указанные параметры принимают следующие значения т = 4, д = 9, г = 8 и = 0,630, X = 2,571.
Алгоритм моделирования. С целью верификации предложенных выражений для кинетического коэффициента диффузии может быть выполнено компьютерное моделирование диффузионных процессов по методу Монте-Карло с помощью алгоритма Метрополиса [8], модифицированного с целью учета взаимодействия в седловой точке.
В рамках этого алгоритма случайным образом выбирается узел г, занятый частицей. После этого также случайно определяется направление возможного прыжка частицы в один из ближайших узлов у. Если второй выбранный узел занят частицей, то переход частицы в него, очевидно, невозможен. Тем не менее попытка такого перехода учитывается. Если же он свободен, то переход частицы в него осуществляется с вероятностью
Ра0 _р-1ехр{-в[Jъsъ- Js]}, (5)
где Р0 - нормировочный коэффициент, равный exp(-3PJ) для системы с притяжением ближайших соседей и exp(4PJ) для системы с отталкиванием, его физический смысл состоит в том, чтобы наиболее энергетически выгодный переход частицы осуществлялся с вероятностью, равной 1; ^ -число ближайших соседей частицы, находящейся в седловой точке; 5 - число ближайших соседей частицы, находящейся в исходном узле О. Например, при перескоке частицы из узла О в узел а
• _ ПЬ + Пс + Пё + Пе + ПГ + ng + « ;
(6)
_ ПЬ + Пё + Пе + Па+ Пв+ «.
Если Рг > Ра0, где Рг - случайное число из диапазона [0; 1], то переход частицы между узлами не осуществляется, в противном случае он считается произошедшим. Повторение данной процедуры п раз, где п - число частиц на решетке, формирует один шаг алгоритма Монте-Карло (МКШ).
Для моделирования диффузионных процессов использовалась решетка с периодическими граничными условиями, содержащая 123 = 1728 решеточных узлов. Процедура моделирования состояла из 50 000 МКШ. Дополнительно первые 10 000 МКШ отводились на эквилибризацию системы и не учитывались в дальнейшем.
Определение коэффициента диффузии решеточного флюида с взаимодействием в седловой точке.
31
Последующее усреднение по 1000 траекторий позволяло изучить зависимость среднего квадрата смещения центра масс системы частиц и среднего квадрата смещения отдельной частицы от времени, измеренного в МКШ, и определить соответственно коэффициент кинетической диффузии DJ и коэффициент диффузии меченых атомов Dtr.
Для снижения влияния размеров моделируемой системы на получаемые результаты использовались периодические граничные условия.
Результаты моделирования. На рис. 2 представлены зависимости от концентрации кинетического коэффициента диффузии решеточного флюида с взаимодействием ближайших соседей и равным ему взаимодействием в седло-вой точке и проведено сопоставление результатов моделирования с данными, полученными на основании соотношения (3).
1
а
х
(D
о
\
0,1
0,01
1Е -3
0,0 0,2 0,4 0,6 0,8 1,0 с
а
^ 0,1
4
^
а
х
(D
Q 0,01
1Е -3
0,0 0,2
0,4 0,6
б
0,8 1,0 с
Рис. 2. Зависимость от концентрации кинетического
коэффициента диффузии решеточного флюида с притяжением (а) и отталкиванием (б) ближайших соседей и равным ему взаимодействием в седловой точке на ОЦК-решетке: 1 - Т/ Тс = 1,05; 2 - Т/ Тс = 1,20; 3 - Т / Тс = 2,00; 4 - Т / Тс = 6,00. Точками представлены результаты МКМ, линиями -результаты использования соотношения (3)
Данное сопоставление ясно показывает, что в случае решеточного флюида на ОЦК-решетке развиваемый приближенный подход может с успехом использоваться для определения коэффициентов диффузии при не очень низких температкрах (T > 1,50TC).
Также можно отметить, что в области низких концентраций (с < 0,20) соответствие результатов аналитических расчетов и данных моделирования является весьма хорошим даже при низких температурах. Такое совпадение результатов может быть объяснено тем, что при низких концентрациях корреляции в заполнении решеточных узлов проявляются слабо ввиду малого количества примесных частиц.
Литература
1. Грода Я. Г. Коэффициент диффузии решеточного флюида с взаимодействием в седло-вой точке: суперпозиционное приближение // Труды БГТУ. 2011. № 6 (144): Физ.-мат. науки и информатика. С. 27-30.
2. Грода Я. Г. Диффузия решеточного флюида на плоской треугольной решетке с учетом взаимодействия в седловой точке // Труды БГТУ. 2012. № 6 (153): Физ.-мат. науки и информатика. С.51-53.
3. Грода Я. Г. Коэффициент диффузии решеточного флюида с взаимодействием в сед-ловой точке на простой кубической решетке: результаты моделирования и суперпозиционного приближения // Труды БГТУ. 2013. № 6 (162): Физ.-мат. науки и информатика. С. 27-30.
4. Collective surface diffusion on triangular and square interacting lattice gases / A. Danani [et al.] // Surf. Science. 1998. Vol. 409. P. 117-129.
5. The self-consistent diagram approximation for lattice systems: diffusion properties of interacting lattice gases / G. S. Bokun [et al.] // Physica A. 2000. Vol. 296, No. 1/2. P. 83-105.
6. Вихренко В. С., Грода Я. Г., Бокун Г. С. Равновесные и диффузионные характеристики интеркаляционных систем на основе решеточных моделей. Минск: БГТУ, 2008. 326 с.
7. Vikhrenko V. S., Groda Ya. G., Bokun G. S. The diagram approximation for lattice systems // Phys. Let. A. 2001. Vol. 286, No. 2/3. P. 127-133.
8. Uebing C., Gomer R. A. Monte Carlo study of surface diffusion coefficients in the presence of adsorbate-adsorbate interactions // J. Chem. Phys. 1991. Vol. 95, No. 10. P. 7626-7652.
Поступила 20.03.2014
1
