CC BY
28
УДК 531.19
Я. Г. Грода, кандидат физико-математических наук, доцент (БГТУ)
КОЭФФИЦИЕНТ ДИФФУЗИИ РЕШЕТОЧНОГО ФЛЮИДА С ВЗАИМОДЕЙСТВИЕМ В СЕДЛОВОЙ ТОЧКЕ: СУПЕРПОЗИЦИОННОЕ ПРИБЛИЖЕНИЕ
Рассмотрен процесс диффузии на плоской квадратной решетке решеточного флюида с дополнительным взаимодействием в седловой точке. При исследовании термоактивированных прыжков частиц получено выражение для кинетического коэффициента диффузии решеточного флюида. В рамках суперпозиционного приближения получено аналитическое выражение для оценки кинетического коэффициента диффузии решеточного флюида. Определен кинетический коэффициент диффузии решеточного флюида с притяжением ближайших соседей и равным ему взаимодействием в седловой точке. Результаты аналитических расчетов сопоставлены с данными компьютерного моделирования диффузионного процесса в исследуемой системе по методу Монте-Карло. Установлены границы применимости предложенных аналитических соотношений и пути построения приближений следующих порядков.
The diffusion process of the lattice fluid with additional interaction in the saddle point on a square lattice is considered. The expression for the jump diffusion coefficient in the case of thermally activated particle jumps is derived. A simple analytical expression for the jump diffusion coefficient is deduced within the superposition approximation. The jump diffusion coefficient of the system with nearest neighbor attractive interaction and interaction in the saddle point of the same energy is calculated. The analytical results are compared with the Monte Carlo simulation data. The boundaries of applicability of the jump diffusion coefficient expression are discussed. The ways to construct more sophisticated approximations are discussed.
Введение. Модель решеточного флюида является одной из стандартных моделей физики конденсированного состояния и широко используется для описания физико-химических процессов в объеме и на поверхностях твердых тел [1]. В частности, она оказывается очень полезной при изучении диффузионных процессов в данных системах.
Ранее была предложена общая теория диффузии в решеточных системах и дан строгий статистико-механический вывод выражения для кинетического коэффициента диффузии в пренебрежении влиянием эффектов памяти [2]. В рамках построенной теории были исследованы диффузионные свойства решеточных флюидов с взаимодействием ближайших соседей на решетках различных типов [2, 3]. В дальнейшем отмеченный подход был расширен для описания процессов термодиффузии [4], диффузии в многокомпонентных [5] и многоуровневых [6] системах, диффузии решеточного флюида на неупорядоченных решетках [7, 8].
В то же время можно отметить, что при преодолении частицей межузлового барьера частица попадает в так называемую седловую точку, расположенную в вершине этого барьера. При этом может быть рассмотрено взаимодействие частицы с частицами, являющимися ближайшими соседями по отношению к данной седловой точке. Учет такого дополнительного взаимодействия приводит к изменению эффективного межузлового барьера и, очевидно, будет влиять на диффузионные свойства системы [9].
В настоящей работе сделана попытка учета такого дополнительного взаимодействия между частицами и проводится оценка его влияния на диффузионные свойства решеточного флюида.
Модель. Рассматриваемая в работе модель представляет собой систему из п частиц, расположенных в узлах регулярной плоской квадратной решетки, содержащей N узлов. Каждый узел может быть либо занят частицей, либо быть вакантным. Состояние узла 1 определяется числом заполнения п = 1 или п = 0 в зависимости от того, занят узел частицей или вакантен. Заполнение узла более чем одной частицей запрещено.
С1 С Ь
е 0 1а
f д Н
Рис. 1. Плоская квадратная решетка.
Узлы а, с, е и д - ближайшие соседи узла 0.
Узлы с, Ь, И и д - ближайшие соседи седловой точки I
Находящаяся в узле 0 частица может взаимодействовать с энергией J с частицами, занимающими ближайшие соседние узлы, т. е. узлы а, с, е ид (см. рис. 1). При ее последующем
переходе в узел а при прохождении седловой точки I она также взаимодействует с узлами с, Ь, И и д. Энергия взаимодействия в данном случае принимается равной и в общем случае J Ф
Таким образом, для перехода из узла 0 в ближайший вакантный узел а частице необходимо преодолеть активационный барьер Еа, величина которого может быть определена как
где
Ea = E(!) - E(0),
(1)
где Е(0) - начальная энергия частицы, равная
Е (0) = 3 (пс + пе + п8); (2)
Е(Х) - энергия частицы в седловой точке, вычисляемая по формуле
Е (!) = Ео + 3Е (пъ + Пс + пг + пк), (3)
здесь Е0 - исходная высота межузельного барьера.
Следовательно, активационный барьер может быть записан в следующем виде:
Еа = Е0 - -Ч - А(пс + п& ) + 3^ (пъ + пк X (4)
где А = - - -х.
Кинетический коэффициент диффузии решеточного флюида. В случае произвольной решеточной системы, в которой возможны переходы только между ближайшими узлами, пренебрегая влиянием эффектов памяти и пространственной дисперсии, для коэффициента кинетической диффузии может быть записано следующее соотношение [2]:
Dj =
zwa
2d
W = c Wa 0 П0(1 - Па ^
(5)
где г - число ближайших соседей на решетке рассматриваемого типа; м> - средняя вероятность перехода частицы; а - расстояние между узлами решетки (длина прыжка частицы); й - размерность пространства; с - равновесное значение концентрации частиц; wа0 - частота перескока частицы из узла 0 в узел а, рассчитываемая по формуле
Wa0 = V exp(-ßEa X
(6)
здесь V - частота, имеющая порядок частоты колебаний частицы вблизи узла решетки и определяющая временную шкалу диффузионных процессов; р = 1 / квТ - обратная температура; кв - постоянная Больцмана; Т - температура.
С учетом алгебры чисел заполнения
Па (1 - Па ) = ^ n2 = Па ,
(7)
для частоты перескока ^а0 получаем
0 = (1 + пва)(1 + пс у)(1 + пг у) х
X (1 + п£)(1 + п,Д), (8)
с = exp(ßJ) -1, у = exp(ßA) -1, Е = exp(-ßJЕ) -1.
(9)
Таким образом, для кинетического коэффициента диффузии можем записать следующее выражение:
Dj = D («0 (1 - na )(1 + nec)(1 + ncy)(1 + n у) x c
X (1 + nbE)(1 + nhE)), (10)
D0 = — w-ßE0. 0 2d
(11)
Определяя корреляционные функции для заполнения решеточных узлов соотношением вида
(пп12п3 ■■■ \:>)=спёп^^-•, 1пх (12)
окончательно для кинетического коэффициента диффузии решеточного флюида получаем следующее выражение:
Dj
Dn
= 1 - cg2 (0, a) + cc (g2 (0, e) - cg3 (0, a, e)) +
+ 2yc (g2(0, c) - cg3 (0, a, c)) + + Y2c2 (g3 (0, c, g) - cg4 (0, a, c, g)) + + 2cyc2 (g3 (0, c, e) - cg4 (0, a, c, e)) + + cy2c3 (g4 (0, c, e, g) - cg5 (0, a, c, e, g)) + + 2Ec ((0, b) - cg3 (0, a, b)) + + 2cEc2 (g3 (0, b, e) - cg 4 (0, a, b, e)) + + 2y^c2 (g3(0, b, c) + g3(0, b, g) -- cg4(0, a, b, c) - cg4(0, a, b, g)) + + 2yV (g4 (0, b, c, g) - cg5 (0, a, b, c, g)) + + 2cyEc3 (g4(0, b, e, c) + g4(0, b, e, g) -- cg5 (0, a, b, e, c) - cg5 (0, a, b, e, g)) + + 2cy2Ec4 (g5 (0, b, e, c, g) - cg6 (0, a, b, e, c, g)) + + 2yE2c3 ((0, b, c, h) - cg5 (0, a, b, c, h)) + + cE2c3 (g4(0, b, e, h) - cg5 (0, a, b, e, h)) + + E2c2 (3(0, b, h) -cg4(0, a, b, h)) + + Y 2fc4 (g5 (0, b, c, g, h) - cg6 (0, a, b, c, g, h)) + + 2cyE2c4 (g5 (0, b, c, e, h) - cg6 (0, a, b, c, e, h)) + + cy2E2c5 (g6(0, b, c, e, g, h) -
- cg7 (0, a , b, c, e, g, h)).
(13)
Можно также отметить, что при рассмотрении решеточного флюида без учета взаимодействия в седловой точке (у = а, £ = 0) полученное выражение (13) переходит в известное выражение для кинетического коэффициента диффузии решеточного флюида [10].
Суперпозиционное приближение. Из полученного соотношения видно, что для вычисления кинетического коэффициента диффузии необходим учет корреляций в заполнении до семи решеточных узлов. Наиболее подходящим методом для определения подобного рода корреляционных функций является метод вариации кластеров Кикучи [11]. Однако в силу того, что необходимо учитывать корреляции в заполнении достаточно большого числа узлов, его применение сопряжено со значительными вычислительными трудностями.
Поэтому в качестве первого приближения можно рассмотреть суперпозиционное приближение, в котором непосредственно учитываются лишь парные корреляции в заполнении ближайших соседних узлов, а корреляционные функции более высоких порядков определяются посредством парных корреляций.
В рамках данного подхода корреляционная функция %3(0, а, с) может быть представлена как
&з(0, а, с) = £>(0, а)^(0, с)а, с) = %2, (14)
где % - парная корреляционная функция заполнения двух ближайших узлов частицами:
%2(0, а) = %2(0, с) = %, %2(а, с) = 1, (15)
так как узлы 0-а и 0-С являются ближайшими соседями, а а-с - соседями второго порядка, корреляциями в заполнении которых мы пренебрегаем. Аналогичным образом можно показать, что
£з(0, Ь, к) = с3, £б(0, ь, с, е, £, к) = £5 (16)
В рамках предлагаемого подхода на основе соотношения (13) может быть получено следующее выражение для кинетического коэффициента диффузии решеточного флюида:
О 2
О = (1 -0£ )(1+У0£ )2 (1+ае%)+ О0
+ 2£0(1 -0%2 )(1 + а0% )х х [1 + у 202 %3 +у0% (% +1)] + + £202 (1 - 0%3) (1 + у0%2 )2 (1 + а0%). (17)
Дополнительно можно отметить, что если предположить отсутствие любых корреляций в заполнении решеточных узлов, т. е. положить % = 1, то соотношение (17) принимает вид
О 2 2
О- = (1 -0) (1 + а0) (1 + у0)2 (1 + £0)2. (18)
О0
Данное выражение может быть получено непосредственно из соотношения (10) в рамках мультипликативного приближения, эквивалентного условию отсутствия корреляций в заполнении решеточных узлов.
Входящая в соотношение (17) парная корреляционная функция двух ближайших соседних узлов % может быть найдена, например, в рамках диаграммного приближения [12].
Результаты моделирования. Для верификации предложенных выражений для кинетического коэффициента диффузии может быть выполнено компьютерное моделирование диффузионных процессов по методу Монте-Карло с помощью алгоритма Метрополиса [13], модифицированного с целью учета взаимодействия в седловой точке.
На рис. 2 представлена зависимость от концентрации кинетического коэффициента диффузии решеточного флюида с притяжением ближайших соседей и равным ему взаимодействием в седловой точке и проводится сопоставление результатов моделирования с данными, полученными на основании соотношения (17).
^ 1
са
С-}
8 0,1
/-ч
5°
сГ 0,01
1Е-3
1Е-4
Рис. 2. Зависимость от концентрации кинетического коэффициента диффузии решеточного флюида с притяжением ближайших соседей и равным ему взаимодействием в седловой точке на плоской квадратной решетке: 1 - Т/ Тс = 1,05; 2 - Т/ Тс = 1,20; 3 - Т/ Тс = 1,50; 4 - Т/ Тс = 2,00; 5 - Т/ Тс = 6,00. Точками представлены результаты МКМ, линиями - результаты использования соотношения (17). Пунктирной линией показаны результаты для решеточного газа Ленгмюра (J = Js = 0)
Данное сопоставление ясно показывает, что предлагаемое суперпозиционное приближение
может с успехом использоваться для определения кинетического коэффициента диффузии при температурах T > 1,50^. Ниже указанной температуры имеет место заметное расхождение между результатами аналитических расчетов и данными моделирования.
Также можно отметить, что и при низких температурах аналитические результаты совпадают с результатами моделирования при низкой (с < 0,10) и предельно высокой (с > 0,95) концентрациях частиц.
Такое соответствие результатов двух различных подходов может быть объяснено тем, что соотношение (17) получено путем выражения многочастичных корреляционных функций через парные функции для ближайших соседних узлов, что, безусловно, является достаточно грубым приближением. В то же время вследствие термоактивированости переходов частиц с ростом температуры корреляции в заполнении узлов ослабевают, что и отражается в практически полном совпадении результатов обоих методов. В области низких концентраций влияние корреляционных эффектов мало просто в силу малости частиц, а при высоких концентрациях определяющим корреляционным эффектом является эффект блокировки, который может описываться концентрацией частиц на решетке.
Также можно отметить, что при определении кинетического коэффициента диффузии в соответствии с соотношением (17) при температуре 1,05^ в области низких концентраций имеет место немонотонная зависимость коэффициента диффузии от концентрации. Сопоставление с результатами моделирования на данный момент не подтверждает наличие такого эффекта. Однако следует признать необходимым осуществление более детального компьютерного моделирования диффузионных процессов в указанном диапазоне температур и концентраций.
Сравнение коэффициента диффузии исследуемой системы с аналогичным коэффициентом для системы без взаимодействия в седло-вой точке показало, что, хотя характер зависимости кинетического коэффициента диффузии от концентрации в целом сохраняется, включение в рассмотрение взаимодействия в седловой точке резко ослабляет эту зависимость.
Заключение. Подводя итог, можно сделать вывод о том, что использованное суперпозиционное приближение может служить отправной точкой для построения более точных соотношений для оценки диффузионных свойств системы. Основным подходом к построению приближений следующих порядков является более детальный учет корреляций в заполнении решеточных узлов.
Литература
1. Вихренко, В. С. Равновесные и диффузионные характеристики интеркаляционных систем на основе решеточных моделей / В. С. Вихренко, Я. Г. Грода, Г. С. Бокун. - Минск: БГТУ, 2008. - 326 с.
2. The self-consistent diagram approximation for lattice systems: diffusion properties of interacting lattice gases / G. S. Bokun [et al.] // Physica A. -2000. - Vol. 296, № 1/2. - P. 83-105.
3. SCDA for 3D lattice gases with repulsive interaction / Ya. G. Groda [et al.] // The European Physical Journal B. - 2003. - Vol. 32. -P.527-535.
4. Thermal diffusion and structural relaxation of interacting three-dimensional lattice gases / V. S. Vikhrenko [et al.] // Solid State Ionics. -2003. - Vol. 157. - P. 221-226.
5. Statistical mechanical description and Monte Carlo simulation of diffusion in two-component lattice systems / G. S. Bokun [et al.] // Electrochimica Acta. - 2005. - Vol. 50, № 8. - P. 1725-1731.
6. Groda, Ya. G. Equilibrium and diffusional properties of two-level lattice systems: quasi-chemical and diagram approximation versus Monte Carlo simulation results / Ya. G. Groda, R. N. Las-ovsky, V. S. Vikhrenko // Solid State Ionics. -2005. - Vol. 176. - P. 1675-1680.
7. Diffusion characteristics of particles on energetically disordered lattices / P. Argyrakis [et al.] // Solid State Ionics. - 2008. - Vol. 179. -P.143-147.
8. Грода, Я. Г. Диффузионные свойства ленгмюровского решеточного газа на неупорядоченной решетке с гауссовским распределением высот межузельных барьеров / Я. Г. Грода // Труды БГТУ. Сер. VI, Физ.-мат. науки и информатика. - 2010. - Вып. XVIII. - С. 44-47.
9. Collective surface diffusion on triangular and square interacting lattice gases / A. Danani [et al.] // Surf. Science. - 1998. - Vol. 409. - P. 117-129.
10. Lattice-gas theory of collective diffusion in adsorbed layers / A. Danani [et al.] // Int. J. Mod. Phys. B. - 1997. - Vol. 11, № 19. - P. 2217-2279.
11. Kikuchi, R. A theory of cooperative phenomena / R. Kikuchi // Phys. Rev. - 1951. - Vol. 81, № 6. - P.988-1003.
12. Vikhrenko, V. S. The diagram approximation for lattice systems / V. S. Vikhrenko, Ya. G. Groda, G. S. Bokun // Phys. Let. A. - 2001. - Vol. 286, № 2/3. - P. 127-133.
13. Uebing, C. Monte Carlo study of surface diffusion coefficients in the presence of adsorbate-adsorbate interactions / C. Uebing, R. A. Gomer // J. Chem. Phys. - 1991. - Vol. 95, № 10. -P. 7626-7652.
Поступила 28.02.2011
