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

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

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

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

Добавлен: 10.04.2024

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

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

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

начальным (d) и конечным (v) положениями. Каждый удален от

вектора

перескока на

расстояние р =

-|-j а (рис. 8).

Для

изотропного дебаевского кристалла при

высоких темпе­

ратурах

(екол = kT)\

 

 

 

 

со —

 

 

 

С й б 2 \

 

(62)

 

 

 

 

кТ ) ’

 

 

 

 

 

 

 

 

где vD— дебаевская

частота;

 

 

 

Q:

м_

атомный

объем,

 

 

 

р

 

 

 

 

 

 

 

 

 

б2 =

<72/s a

 

 

 

 

(63)

и s =

(расстояние

диффундирующего

атома

от перевальной

точки);

С — средняя

упругая постоянная

для перемещения:

Р _

15ц(А + 2ц)

 

 

 

 

(64)

U ~

2(2А+7|Л) ’

 

 

 

 

 

 

 

где

X и

р — коэффициенты Ламэ.

 

 

 

 

 

 

 

 

 

П1

_а_ f2

 

Рис. 8. Геометрия скачка

 

в г. ц. к.

решетке:

Рис. 7. Схема изменения потенциальной энер­

d — атом

в исходном по­

гии диффундирующего атома (модель Флинна)

ложении,

v — вакансия,

[17]

[17]

 

Поскольку даже кубические кристаллы анизотропны, вместо выражения (64) лучше написать среднее значение обратного модуля:

J5_ __ _3__ ,

2

, J _

(65)

2С Сц

С11 — Ci2

 

 

 

Величина G = СХ82 в формуле (62) играет роль гиббсовой сво­ бодной энергии перемещения. Для оценки Ет — энергии активации перемещения нужно знать упругие постоянные. В табл. 2 приведено сопоставление расчетных значений Ет с экспериментальными для г. ц. к. металлов, а в табл. 3 — для о. ц. к. металлов. Для о. ц. к. металлов величину Ет оценивали (если это было возможно) по разности энергии активации самодиффузии [2] и энергии образова­ ния вакансий [18]. Для компактности значения энергии приведены

37


 

 

 

 

 

 

 

Т а б л и ц а

2

Применение теории |17] к г. ц. к. металлам

 

 

 

 

 

Металлы

С„

 

Cl 1

Cj2

С*4

ptn

 

Ет

 

 

 

 

 

 

^расчет *

 

 

 

 

 

 

опыт'

 

 

 

 

1012 дин/см2 *

 

эВ

 

эВ

 

 

 

 

 

 

 

 

 

М е д ь ................

1,762

 

0,256

0,818

0,84

 

0,80—1,10

Серебро . . .

1,315

 

0,171

0,511

0,83

 

0,83—0,88

Золото . . . .

2,016

 

0,159

0,454

0,82

 

0,68—0,87

Никель . . .

2,612

 

0,552

1,317

1,42

 

1,25—1,50

Алюминий . .

1,143

 

0,262

0,316

0,83

 

0,40—0,65

Свинец . . . .

0,555

 

0,050

0,194

0,48

 

-0,56

 

* 1 дин = ю -5

Н.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а

3

Применение теории [17] к о. ц. к. металлам

 

 

 

 

 

Металл

Си

т

(с“ -

с44

Ет

Ет,

Е?

 

 

-

с 12);

расчет’

экспер ^экспер,

 

 

 

эВ

эВ

эВ

эВ

 

 

 

Ш12 дин/см2

 

 

 

 

 

Литий . . . .

0,134

0,011

0,100

0,031

0,29

0,15

0,44

 

Натрий . . .

0,062

0,008

0,059

0,039

 

Калий ..............

0,046

0,0042

0,038

0,042

____

___

 

Х ром ................

3,50

1,41

1,01

0,85

____

____

____

 

Ванадий . . .

2,325

0,565

0,460

0,62

____

____

____

 

Железо . . .

2,43

0,606

1,381

0,70

1,40

1,08

2,48

 

Молибден . .

4,80

1,58

1,087

1,73

1,76

2,24

4,0

 

Тантал . . . .

2,665

0,541

0,873

0,95

1,48

2,91

4,39

 

Ниобий . . . .

2,46

0,56

0,287

0,68

2,52

2,04

4,56

 

Вольфрам . .

5,326

1,616

1,631

2,80

2,09

3,14

5,23

 

в электронвольтах 1;

разброс значений, приведенных в

табл.

2,

соответствует данным

разных

авторов.

 

 

 

 

 

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

Интересно то обстоятельство, что для всех г. ц. к. металлов расчет сделан с одним и тем же значением б2 = 0,104 (для о. ц. к. б2 =

= 0,067).

Таким образом, для всех кристаллов с одинаковой решет*

1 1 эВ =

1,6-Ю-1» Дж.

38


кой критическое смещение оказывается постоянной долей длины перескока (около 30% пути до перевальной точки).

По-видимому, можно ожидать дальнейшего прогресса, если на­ учиться учитывать вклад ангармоничности при высоких темпера­ турах, количественно описывать динамику колебаний кристаллов с дефектами и — в идеале — рассчитывать силы сцепления.

Как было отмечено, именно эти трудные вопросы физики кри­ сталлов ограничивают развитие теории диффузии.

4.КОРРЕЛЯЦИЯ МЕЖДУ ПАРАМЕТРАМИ САМОДИФФУЗИИ

ИТЕРМОДИНАМИЧЕСКИМИ СВОЙСТВАМИ МЕТАЛЛА

Поскольку термодинамические свойства 1 системы отражают степень ее устойчивости, а высота барьера (Е) для диффузии опре­ деляется в основном теми же факторами, то представляется вполне естественным существование корреляции между термодинамическими

и диффузионными

характеристиками

if

 

 

кристаллов.

 

 

 

 

 

 

 

 

Следует отметить, что эти связи

 

 

 

носят полуэмпирический характер и,

/, Хг

Х}

как правило,

не

являются

однознач­

*

 

 

ными,

а допускают

исключение,

что и

 

-/и у

 

 

видно из рис.

9. На рисунке сплошной

\ А

 

линией

схематически

изображена за­

\

висимость потенциальной энергии кри­

 

сталла

от

координат

(xlt

x s — ста­

 

бильные положения атомов;

х.2— пере­

 

 

 

вальная точка). Кривую можно охарак­

Рис. 9. Схема,

 

иллюстрирующая

связь между диффузионными и тер­

теризовать

тремя

параметрами: глу­

модинамическими

характеристи­

биной ямы (энергия связи, теплота

ками

 

 

сублимации, Hs), положением миниму-

в решетке,

координаты

мов

(стабильные

положения

атомов

узлов;

период

решетки а = х г х 3)

и кривизной (д12ц>/дх2) в

окрестности минимума (упругие модули). Учитывая

сказанное

выше,

 

максимумы на кривой соответствуют перевальным точкам 2.

Как видно из рис. 9, энергия активации диффузии представляет собой разницу между энергиями основного и переходного состояний:

Е = Н* Hs.

(66)

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

Отметим также, что наилучшим образом объемное состояние кристалла отражают его свойства при высоких температурах (вблизи

1 Свободная энергия, теплосодержание, теплоты плавления, испарения, сило­

вые постоянные — модули

и т. д.

2 Разумеется, на схеме

не отражена релаксация системы при переходе атома

в перевальную точку.

 

39



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

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

щения вакансий Е = #£ + Н™«=; Efv + Е%, справедливая для вакансионного механизма. Это важно отметить потому, что в основе всех оценок лежит результат, впервые полученный Бертом и Зинером: изменение свободной энергии перемещения равно работе обратимого (изотермического и изобарического) перемещения атома, совершаемого против сил сцепления (чаще всего предполагается, что это — чисто упругие силы). Таким образом, речь идет только об энергии перемещения. Для твердых растворов внедрения Е = Ет, но для твердых растворов замещения образование дефекта (вакансии) не учитывается и это делает оценки еще более приближенными.

Энергия активации самодиффузии

В главе первой уже упоминалось о связи между энергией акти­ вации самодиффузии и теплотой испарения; для большинства ме­ таллов с г. ц. к. решеткой первая составляет 0,6—0,7 от второй.

Однако наибольшее распространение получили корреляционные соотношения между Е и параметрами плавления.

Самым популярным из них является правило, обнаруженное Бугаковым и Ван-Лимптом:

Е =

32ТПЛ,

(67)

где

[Е] = кал/моль;

[Гпл] = °К.

Это соотношение очень похоже по форме на известное правило

Трутона

в

термохимии.

правило, связывающее Е

с Ьт

По содержанию близко к (67)

теплотой

плавления:

 

 

Е = 16,5Ьпл.

 

(68)

Поскольку

энтропия плавления

для большинства металлов

5 ПЛ =

= Епл/Тпл ^

2,3 кал/моль -град,

коэффициенты в формулах

(67) и

(68) противоречат друг другу. Статистический анализ данных для

большого числа элементов (табл. 4, см. также рис. 1)

позволил

привести их в соответствие; оказалось, что точнее писать

 

Е =

35Тпл

. (69)

и

 

 

Е =

15,2БПЛ.

(70)

Формулу (69) иногда пишут в виде

 

Е =

Ж Т ПЛ,

(71)

что позволяет не думать о размерности.

Расчетные значения совпадают с экспериментальными с точностью ±20% , почти не выходящей за пределы воспроизводимости экспе­ риментов.

40