Материал: Диффузия в кристаллических телах. Юрьев В.А., Юрьева М.В

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

барьера), следует учитывать потенциальную энергию группы атомов этого ансамбля. В каждый момент есть равновесие между атомами в исходном состоянии 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

Источник: https://studfile.net/preview/16565967/