Файл: Бокштейн Б.С. Термодинамика и кинетика диффузии в твердых телах.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