Файл: Бокштейн Б.С. Термодинамика и кинетика диффузии в твердых телах.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 10.04.2024

Просмотров: 147

Скачиваний: 4

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Более общее выражение связывает коэффициент диффузии с дли­ ной (А) и суммарной частотой скачков (Г):

(227)

Аналогично тому, как это было сделано для кубических кристал­ лов, можно рассмотреть и кристаллы с другой решеткой. Например, для гексагональных кристаллов [21:

° х = 4" IV r

(228)

и

 

D у= Г„с2.

(229)

Здесь индексы _L и || означают, что скачки происходят перпендику­

лярно и параллельно оси с кристалла; с и а — периоды

решетки

вдоль оси с и в плоскости базиса.

 

Сделанные выводы справедливы только для случайных блужда­ ний. Еще раз напомним, что блуждания являются случайными, если перемещение каждой частицы не зависит от движения других частиц и если направления последовательных перескоков не зависят друг от друга. Как правило, это не так. Действительно случайными можно, по-видимому, считать только перемещения вакансии или межузель­ ного атома при самодиффузии. Во всех остальных случаях указанные условия не выполняются и перемещения частиц не являются вполне случайными. В этом случае в выражении для D появляется фактор корреляции (п. 4).

Сделанный расчет справедлив, если на атом не действуют ника-

кие внешние силы. В противном случае (L) ф 0. При учете внеш­ них (движущих) 1 сил выражение (222) для вероятности обнаружить частицу в окрестности точки х в момент времени t (или же выражение (205) для концентрации) приобретает вид уравнения Планка—Фок- кера:

— =

d ( W с)

 

п д%с

(230)

d t

д х

+

д х * '

 

Здесь (и )— средняя скорость, которую частица приобретает под действием движущей силы. По-прежнему D ф D (х). Уравнение

(230)легко обобщить на случай многих переменных.

2.ЗАДАЧА О ДОСТИЖЕНИИ ГРАНИЦ

Во многих диффузионных задачах не требуется точного знания распределения концентрации, а достаточна менее детальная характе­ ристика — время пребывания частицы в некоторой области. Это —• частный случай задачи о вероятности достижения границы области;

1 Как мы покажем ниже (гл. V), движущая сила может и не быть внешней; в этом случае она связана с неидеальностью раствора.

9 Заказ № 737

129



при этом граница может быть отражающей (частица возвращается в область) и поглощающей (частица теряет способность к дальней­ шему блужданию). Среднее время, необходимое для того, чтобы ча­ стица первый раз достигла границы области, называется временем первого перехода (first passage time) — t (0). Очевидно, что время первого перехода растет с увеличением расстояния от границы обла­ сти до исходного положения частицы в области.

В такой формулировке задача нахождения первого времени пере­ хода применима к широкому кругу проблем физико-химической кине­ тики: от исследования диссоциации двухатомных молекул, которая происходит, если молекула набирает благодаря соударениям энергию, большую некоторого критического значения Е, и броуновского дви­ жения до ухода звезды из звездного скопления. В подавляющем большинстве задач рассматриваются марковские системы, не облада­ ющие памятью о прошлых состояниях. В этом случае частице, уже характеризующейся некоторым значением параметра л: может быть координатой, импульсом, энергией и т. д.), безразлично, каким образом это значение достигнуто, дальнейшее течение процесса от

этого не зависит.

расчета

t

(0), следуя обзору Вейса 1.

Как и

Рассмотрим схему

в п. 1, обозначим через W (t,

х'

х) — вероятность того,

что ча^

стица из точки х при t

= 0 перейдет в х' за время t (в общем случае х'

и л; могут быть не только координатами). W удовлетворяет дифферен­

циальному уравнению

 

 

 

 

- = - L xW,

 

 

 

(231)

которое может быть (в случае диффузии, например) уравнением Фика

или Планка—Фоккера.

Соответственно: — L* =

(

D )

или

~

~

~

^ (v)-

Уравнение

(231) справедливо

как

для

непрерывного, так и дискретного х.

останется в

пределах

области

 

Вероятность

того, что частица

(Е) к моменту t:

 

 

 

 

 

W(x', t) =

J W(t, x' x)dx.

 

 

 

(232)

v

Пусть r) (x', t) — плотность вероятности, соответствующая попа­ данию первого времени перехода в интервал t , t 4- dt, т. е. t < t (0)<; t + dt. Если частица находится в пределах V к моменту t, то она либо переходит в другую область (V) в следующий момент t + dt,

либо остается в пределах области V. Поэтому

W(x', t) = r\(x', t)dt +

W (x',t + dt).

(233)

1 W e i s s J. La Jolla

Summer School on Chemical Physics,

August 1965, Univ.

of California.

 

 

130


Или

г](л:\ t) =

dW (х', t)

(234)

dt

 

 

 

Момент 1 r -го порядка первого времени перехода по определению ра­ вен

СО

03

 

(р (х')) = f trx\ (х', t) dt =

r \ tr~lW (x, t) dt.

(235)

о

о

 

Последнее равенство получено интегрированием по частям, с учетом

выражения (234). В формуле (235) г = 1, 2, ...

Формально решение уравнения (231) независимо от конкретного

вида оператора Lx можно записать как

 

W (#, *' — * ) =

exp {—Lxt) А х'),

(236)

где

 

б х') — для непрерывного х

 

А х') =

 

SXx' — Для дискретного х,

 

 

 

 

т. е. либо 8-функции, либо оператору Кронекера.

 

Подставляя выражение (236) в (232) и (235), получим

 

(Р(х)) =

Г

Jt ^ d t j e x p i — L J ) A ( x - x ’)dx = r! J L7rA(x — x)dx.

 

 

0

v

°

(237)

Если оператор Lx представить в виде матрицы Lmn (для дискретных

значений х\

х ’ =

т), то уравнение (232) примет вид:

 

W (t, т п) = exp (— Lt)nfirm=

exp (— Lt)nm

(238)

и

 

 

 

 

 

(tr im)) =

г !

S

)nfi]m = r l h

{L-r)nm-

(239)

 

 

n

n

 

 

 

 

(V)

(V)

 

 

Чтобы получить решение в аналитической форме, необходимо сде­ лать дополнительные упрощающие предположения относительно опе­

ратора

L.

Чаще всего

полагают,

что можно ограничиться учетом

только

ближайших переходов: соответственно Ьц — 0 для всех

I i / 1 >

1.

 

 

1 Моментом г-го порядка

величины х

называется ее математическое ожидание

(Мхг) или среднее значение (с учетом вероятности, т. е. центр распределения).

Мхг

= ^ xrW (х) =

j

(х) dx,

где

W (х) — плотность

вероятности данного значения х, нормированная на еди-

 

 

оо

 

 

ницу, т. е.

J IV (х) dx = 1.

 

 

о

 

91

131

 


Типичным примером может служить задача о диссоциации двух­ атомных молекул [116].

В первоначальном рассмотрении предполагали, что реагирующие молекулы являются простыми гармоническими осцилляторами и под­ чиняются статистике Больцмана. Слабое взаимодействие этих моле­ кул с молекулами термостата является причиной их переходов из одних энергетических состояний в другие. Ограничившись рассмо­ трением переходов только в ближайшие энергетические состояния, вместо выражения (231) можно было написать (для дискретного энер­

гетического спектра) систему уравнений:

 

Ц7о = Б101 Г 1 -Б 01Го,

 

^ 1 L01W0

(L10 -f- Lu ) W1-f- L2iWg,

 

^ 2 = ^12^1

(A l L ZS) W 2 -j- LgalFg,

(240)

W n ~ L.n -1, n W n _ i (Lm , n i + L n . n + i ) W n . .

Здесь W( означает вероятность попадания в состояние с энергией Et. Предполагается, что переход в (N + 1)-е энергетическое состояние приводит к реакции диссоциации (молекула уходит из области V).' Таким образом, среднее время попадания частицы в это состояние и есть первое время перехода.

Было также сделано предположение, что выполняется принцип детального равновесия, т. е. скорости прямого и обратного перехо­ дов равны по каждому из возможных путей. Из этого предположения

следовало, что

 

 

 

 

A. i+i .

ехР Y,

 

 

(241)

Т -— :=

 

 

ь!+1. 1

 

 

 

 

v

'

где у =

hv/kT\

v — характеристическая частота осциллятора.

 

Поскольку

W (t) = W0 (/) +

W^t) -]------- Ь WN (t),

из уравне­

ния (235)

следует, что

 

 

 

 

со

 

 

 

 

 

(t (0)) =

|

[W0 it) -j- Wx (t) -j- •••-(- WN (/)] dt.

(242)

 

о

 

 

 

 

 

Используя

преобразования

Лапласа W) (s) = j

e s‘Wj(t)dt,

можно

написать

о

 

 

 

 

 

(t (0)) =

S

WUO).

 

(243)

 

/=o

 

 

 

 

Уравнение (243) оказалось наиболее удобным для вычислений.

Оно привело к

результату

 

 

 

(НО))

 

 

1

 

(244)

k [1 — exp (— 7 )] ■ 2

1

 

 

 

 

 

 

 

 

/=i

 

 

 

в котором k — константа скорости реакции.

132