барьера), следует учитывать потенциальную энергию группы атомов этого ансамбля. В каждый момент есть равновесие между атомами в исходном состоянии n1 и переходном n* , но нет равновесия между переходным и конечным состоянием, что позволяет записать диффузионную «химическую» реакцию в виде:
n x n* x* n |
2 |
x |
2 |
. |
(2.26) |
|
1 |
1 |
|
|
|
||
В соответствии с этим скорость реакции (W ) или суммарное число скачков-переходов всех атомов из состояния x1 в состояние x2 в единицу времени определяется числом
активированных атомов n* в единице объема на отрезке длиной , соответствующем переходному состоянию. Все эти атомы, достигшие переходного состояния, неизбежно доходят
до конца пути |
со |
средней скоростью x |
на |
отрезке в |
||||||
направлении от x1 |
к x2 . |
|
|
|
|
|
||||
Следовательно, |
x |
|
|
|
|
|||||
|
|
* |
|
|
|
|
||||
|
|
W n |
|
|
. |
|
(2.27) |
|||
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
Частота скачков-переходов одного атома |
|
|||||||||
|
|
W /n1 |
|
|
|
|
||||
или |
|
|
|
|
|
x |
|
|
|
|
|
|
|
n* |
|
|
, |
|
(2.28) |
||
|
|
n1 |
|
|
||||||
|
|
|
|
|
|
|
||||
где n1 - |
число атомов в единице объема |
в исходном |
||||||||
состоянии (в окрестности x1 ).
Поскольку среднее число частиц в некоторой области фазового пространства пропорционально сумме состояний Z в единице объема – функции распределения P Z /V , где V -
71
объем. В соответствии с (2.28), атомы в исходном и переходном состояниях находятся в равновесии и можно, воспользовавшись статистической термодинамикой, записать константу равновесия Kc в виде
|
|
|
|
n* / |
|
P* 2 |
|
|
|
|
|
Kc |
|
|
|
|
, |
|
(2.29) |
|
|
|
n |
P 3 |
|
||||
|
|
|
1 |
1 |
|
|
|
||
где P* 2 и |
P 3 |
- функции |
распределения |
атомного |
|||||
|
1 |
|
|
|
|
|
|
|
|
осциллятора соответственно |
в переходном |
и |
основном |
||||||
состояниях. Цифры в скобках показывают, что число степеней свободы в этих состояниях различно, т.е. в переходном состоянии одна колебательная степень свободы заменена поступательной.
Полагая, что кинетическая энергия атомов в яме и на
вершине |
барьера |
одинакова, оставим в выражении |
для P |
|||
только потенциальную энергию, тогда |
|
|||||
|
|
|
|
|
|
|
|
n |
* |
|
exp U x*, y, z, qi /kT d y d z d qi |
|
|
|
/ |
|
|
|
. (2.30) |
|
|
|
n1 |
|
|
|
|
|
|
exp U x, y, z, qi /kT d x d y d z d qi |
|
|||
|
|
|
|
|
|
|
В числители этого уравнения координата перевальной точки x* фиксирована.
В рассматриваемой модели для основного состояния осциллятора принято гармоническое приближение. Для переходного такое приближение недопустимо, поскольку вблизи вершины барьера потенциал не имеет ничего общего с
гармоническим. Разложив U x, |
y, z, qi |
в ряд в окрестности |
||||
x x1, ограничившись |
квадратичным членом в разложении |
|||||
получим |
|
|
|
|
|
|
U x, y, z, q |
U x , y, z, q |
|
x x 2 . (2.31) |
|||
|
||||||
i |
1 |
i |
2 |
1 |
||
72
Здесь |
|
|
2 |
U |
|
- упругая или силовая постоянная. |
|
|
|
|
|||||
|
|
|
x |
2 |
|
|
|
|
|
|
x x1 |
|
|||
При разложении потенциальной энергии учтено, что
|
U |
|
|
. |
|
|
|
0 |
|||
|
|||||
|
x |
|
|
||
|
x x1 |
|
|
После подстановки (2.31) в (2.30) и интегрирования знаменателя получим
* |
|
|
|
1 |
|
|
|
|
|
2 |
|
|
|||
n |
/ |
|
|
|
|
||
|
n |
2 k T |
|||||
|
|
|
|
|
|
||
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
(2.32) |
|
|
exp U x* |
||||||
|
, y, z, qi /kT d y d z d qi |
||||||
|
|
|
|
|
|
||
exp U x1, y, z, qi /kT d y d z d qi
Теперь и в знаменателе зафиксирована координата исходного узла x1 . Поэтому
|
* |
|
|
1 |
|
* 2 |
|
||
n |
|
2 |
P |
|
|||||
/ |
|
|
|
|
. |
(2.33) |
|||
|
n |
2 k T |
P |
|
|||||
|
|
|
2 |
|
|||||
|
1 |
|
|
|
1 |
|
|
||
В отличие от уравнения (2.29), в уравнении (2.33) функции распределения в исходном и переходном состояниях соответствуют одинаковому числу степеней свободы. Поэтому можно воспользоваться известным соотношением статистической термодинамики, связывающим функцию распределения (P ) с - изобарно-изотермическим потенциалом системы (G )
|
|
|
G kT ln P |
|
|
(2.34) |
|||
и записать |
* 2 |
|
|
|
|
|
|
||
|
P |
|
|
G |
|
|
|
||
|
|
|
exp |
|
|
, |
(2.35) |
||
|
P |
|
kT |
||||||
|
2 |
|
|
|
|
|
|||
1 |
|
|
|
|
|
|
|||
73
где G - разница свободных энергий осциллятора, колеблющегося в плоскости yz в точках с координатами x*
(переходное состояние) и x1 (основное состояние).Другими словами G есть работа обратимого изотермического и изобарического перевода атома, свободно колеблющегося в
плоскости yz из x в x* . |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Подставляя выражения (2.25), (2.33) и (2.35) в уравнение |
|||||||||||||||
(2.21), получим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
1 |
2 |
G |
|
|
||||||
|
|
|
(2.36) |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|||||||||
|
|
|
2 m |
|
|
k T |
|
|
|||||||
|
|
|
|
|
|
|
|||||||||
или, выражая частоту колебаний через силовую |
|||||||||||||||
постоянную: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
1 |
|
|
||||
|
|
|
|
|
2 |
|
(2.37) |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
2 |
|
|
|
||||||||
|
|
|
|
|
|
|
|
m |
|
|
|||||
и подставляя (2.37) в (2.36): |
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
G |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
. |
(2.38) |
|||
|
|
|
exp |
k T |
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Поскольку |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
G H T S U p V T S , |
(2.39) |
||||||||||||||
где H , |
S , |
U и |
|
|
V |
- |
|
изменения энтальпии, |
|||||||
энтропии, внутренней энергии и объема при переводе атома из x1 в x* ; p - давление, и учитывая, что в соответствии с уравнением (2.19)
D 2 ,
получим
74
2 |
|
|
S |
|
|
|
|
H |
|
|
|
||||
|
|
|
|
|
|
|
|
||||||||
D exp |
|
|
exp |
|
|
|
|
|
|||||||
|
|
|
k T |
|
(2.40) |
||||||||||
|
|
|
k |
|
|
|
|
|
|
||||||
2 |
S |
|
|
|
U p |
V |
|
||||||||
|
|
|
|
||||||||||||
exp |
|
|
exp |
|
|
|
|
|
|
|
|||||
|
|
|
kT |
|
|
. |
|||||||||
|
k |
|
|
|
|
|
|
|
|
|
|
||||
Сопоставляя этот результат с уравнением (2.23) можно предэкспоненциальный член записать как
D 2 exp S /k , |
(2.41) |
0 |
|
а энергию активации представить через энтальпию
E H U p V . |
(2.42) |
Из данной модели следует, что:
1.коэффициент диффузии экспоненциально зависит от температуры, причем в показатель экспоненты входит не внутренняя (потенциальная), а свободная энергия активации.
2.Энергия активации совпадает с энтальпией активации, т.е. определяется не только изменением внутренней (потенциальной) энергии при переводе атома из основного в переходное состояние, но и активационным объемом, т.е. разницей объемов в этих двух состояниях. Именно активационный объем определяет зависимость коэффициента диффузии от давления
ln D/ p V /k T . |
(2.43) |
3.Предэкспоненциальный фактор определяется (кроме длины перескока и частоты колебаний) энтропией активации.
Уравнения, полученные с помощью теории переходного состояния, имеют большое значение.
В принципе можно рассчитать параметры диффузииD0, E теоретически, задаваясь теми или иными моделями
переходного состояния. Это позволит получить полезные корреляционные соотношения между энергией активации и
75