Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Физико-химические основы металлургических процессов

..pdf
Скачиваний:
8
Добавлен:
25.10.2023
Размер:
17.51 Mб
Скачать

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

 

 

 

ТАБЛИЦА

 

ходу

в класс

металлов.

Координа­

 

 

 

 

ционное число может меняться и при

П А Р А М Е Т Р Ы Б Л И Ж Н Е Г О

дальнейшем

нагревании

жидкости,

П О Р Я Д К А В

П Р О С Т Ы Х

 

 

жидкостях [3]

 

однако низкая точность его опреде­

 

 

 

 

о

 

ления не позволяет пока сделать вы­

В е щ е с т в о

T, °K

 

z

вод о том, как зависит

величина г от

ri,

A

 

 

 

 

 

 

температуры

для

различных

объек­

 

 

1253

2,86

10,0

тов.

 

 

 

 

 

 

 

 

A g

 

Не

следует

думать,

что

число

А1

 

973

2,96

10,6

ближайших соседей у каждого атома

Bi

 

558

3,35

7,6

Cd

 

623

3,06

8,3

жидкости в данный момент времени

Fe

 

1803

2,51

8,4

одинаково. Величина

г подвержена

Ge

 

1273

2,70

8,0

флуктуациям

так же, как

и

другие

Sn

 

513

3,21

8,4

локальные

физические

характери­

 

 

 

 

 

 

 

 

 

 

 

 

стики

жидкости. Теоретические рас­

 

 

 

 

 

 

четы флуктуации

координационного

числа

для жидких

металлов

были

проведены

И.

3.

Фишером

и В. К. Прохоренко

[4]; они достигают

в

среднем

1—3 единиц.

Более детальные

сведения о распределении координационных чисел

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

методов

иссле­

дования

жидкостей

(см. ниже).

 

 

 

 

 

 

 

 

 

 

 

 

К О Р Р Е Л Я Ц И О Н Н Ы Е Ф У Н К Ц И И

 

 

 

 

Выше была определена парная функция радиального распреде­

ления

R

(г), описывающая среднее размещение

частиц

жидкости

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

Выберем в жидкости объем dV. Среднее число частиц в этом

объеме

окажется равным dn = R0 dV,

где R0 — плотность

числа

частиц

(число их в

единице объема).

Если

dV достаточно

мало,

то dn будет гораздо

меньше единицы. Поэтому

произведение

R0 dV

можно рассматривать как вероятность обнаружения частицы жидко­ сти в объеме dV:

da),. =

R0 dV.

(ІГ-35)

Теперь выделим в жидкости

два элемента объема dV1

и dV2

->

->

 

вблизи точек с радиус-векторами гх и г2. Определим вероятность dwb 2 того, что в объемах dV\ и dV2 одновременно будут находиться две частицы. Эта вероятность должна быть пропорциональна самим объемам, т. е.

dwlti = Я,( г ь r a ) dVxdV2,

(11-36)

где функция R2 (rlt

г2 ) — так называемая б и н а р н а я

к о р ­

р е л я ц и о н н а я

ф у н к ц и я .

 

Предположим, что объемы dV1 и dV2 расположены далеко друг от друга. Тогда вероятности попадания частиц в эти объемы стано­ вятся независимыми, и по теореме об умножении вероятностей полу­ чаем

dwi$ — dwi dw2 = RodVi dV2.

(11-37)

Из сопоставления этого соотношения с уравнением (П-36) следует,

что при больших величинах г

— г 2 |

 

/?2 fr,,

" г 2 ) = # § .

(П-38)

Можно установить связь между бинарной корреляционной функ­ цией и парной функцией радиального распределения R (г). Вероят­

ность нахождения двух частиц в объемах dV 1

и dV2 можно

предста­

вить в виде произведения двух

вероятностей: во-первых,

того,

что

в объеме dVx

окажется одна частица {dwy), и во-вторых, что на рас-

стоянии 2гг)

от первой частицы окажется вторая

(dw?) (таким

образом, dw'2 — это условная

вероятность

второго

события,

если

уже наступило первое). Величина dw1 равна

по-прежнему

R0

dVx,

a dm>2 в соответствии с определением (П-З) функции радиального рас­

пределения

 

находится из

соотношения

 

 

 

 

Итак:

 

 

dm =

R(\r2—~rl\)dV2.

(11-39)

 

 

 

 

 

 

 

 

 

 

 

 

<toli2 = R0R

(I r 2 гх

I) dVi. dV2.

 

(Ц-40)

Из формулы

(П-38)

и (11-40) видно, что

 

 

 

 

 

ЯлъГг)

 

= яМ\гг-Щ

 

(Н-41)

Аналогично сказанному выше, можно ввести вероятность dw

2

того, что

в

объемах

dVlt

dV%,

dVn

окажется по одной частице.

Определим

 

dwh 2

п следующим

образом

 

 

 

 

 

dwht

„ = Rn(ru

r! f ..., гя)

dV^Vz,dVn.

(II-42)

Введенная уравнением (11-42) функция координат всех п частиц

Rn

называется

к о р р е л я ц и о н н о й

ф у н к ц и е й

n-г о

п о -

р я д к а .

Как и раньше, при увеличении расстояния

| rt—rk \ между

каждыми двумя объемами вероятности нахождения частиц в каждом

 

 

-> -»

объеме dVt

становятся независимыми, т. е. при \ rtrk

\ » 00

Очевидно,

что Rx = R0.

Частный случай корреляционной функции

п-го порядка — это так

называемая тернарная

корреляционная

71

функция R3 (r 1 (

r2, r3),

описывающая

вероятность данной

конфигу­

рации тройки частиц.

 

 

 

 

 

 

 

Вероятности

dw, введенные

выше,

нормированы следующим

образом:

 

 

 

 

 

 

 

 

 

\dw1=\RbdV1

= R0V

= N,

 

(II-43)

 

 

v

 

 

 

 

 

 

где N — полное

число

частиц

жидкости;

 

 

 

\dw2 = \R2

( г ь r a ) dV1 dV2

= j J/?0

Л У ? (|г, -

PI |) dV2

=

 

 

 

 

V V

 

 

 

 

 

= R0V\

R (I r2

-

rx I)

dV2

= N(N-

1)

(11-44)

 

v

 

 

 

 

 

 

 

и соответственно,

 

 

 

 

 

 

 

 

 

W = W = i o T -

 

 

( , I " 4 5 )

Целесообразно переформулировать корреляционные функции и вероятности так, чтобы первые стали безразмерными, и чтобы нор­ мировочные интегралы равнялись единице. С этой целью введем новые вероятности заданной конфигурации п частиц жидкости:

dWn =~^Fn (rlt r2,...,rn)

dv1 dv2,...,dvn,

(П-46)

где dvt — dxt dyt dzt; xit y{, z{ — координаты /-той частицы. Определенная таким образом корреляционная функция F„ без­

размерна. Далее потребуем, чтобы новая вероятность нормировалась на единицу

\dWn=l,

т. е.

 

 

- ^ - JX l t

га ,.... rn) dv1 dv2, ...,dvn=l.

(II-47)

Из уравнений (11-45) и

(11-47) видно, что

 

d w » = W ^ j T d w » -

( I M 8 >

Уравнение (11-48) подсказывает правильную трактовку вероят­

ностей dwn и dWn.

Вероятность dwn

относится к данной конфигура­

ции

объемов

dVlt

dV2,

когда

несущественно, какие

именно

частицы

занимают

эти

объемы. Наоборот, вероятность dWn

отно­

сится

к

ситуации,

когда в элементе

объема dv1

находится

первая

частица, в элементе dv2

— вторая и т. д. Вероятность dWn не вклю­

чает

в

явном

виде

принципа неразличимости

микрочастиц.

Если

учесть,

что все частицы

одинаковы,

то объем dV1

может быть занят

72

любой

из N частиц,

а объем dV2

— любой

из оставшихся N 1

частиц и т. д.; вследствие этого и появляется

дополнительный мно­

житель

NU(N — я)!.

 

 

 

 

 

 

 

Учитывая теперь соотношения (II-36), (11-46) и (11-48), находим

связь между корреляционными функциями Rn

и Fn:

 

 

 

 

Rn\ru Г,

п

N1

1 г

2.--->

Г

п)

(II-49)

 

 

Г )

(N — n)\ у"-

 

 

Отметим теперь, что корреляционные функции двух

последую­

щих порядков связаны между собой соотношением [5]:

 

 

l/V

JFndvn

= Fn_1

Х ) г2 ) ...,

гп_х),

 

 

(II-50)

вытекающим из определения вероятности dWn (П-46). Кроме того, функции Rn и Fn симметричны относительно перестановки любых

> >

двух одинаковых частиц или радиус-векторов rt и rk.

Введем теперь полную потенциальную энергию системы U (rlt

-*

г 2> • • •>гы)< зависящую от координат всех частиц жидкости. Тогда, согласно статистике Больцмана — Гиббса, вероятность данной кон­ фигурации частиц пропорциональна фактору ехр (—U/kT). Следо­ вательно:

dW N • А ехр

U

( r u r2,...,

t

dvxdv2, ...,dvA

(11-51)

 

kT

 

 

 

 

 

 

где А — нормировочный

множитель.

 

 

 

Из уравнения (П-46)

получаем, что

 

 

FN (ri,ra,.~, rN)

=

VNAехр

 

и

(11-52)

 

kT

 

 

 

 

 

Итак, корреляционная функция FN определяется потенциальной энергией системы. Рекуррентные соотношения (II-50) позволяют в принципе найти корреляционные функции низших порядков Fn. Вычисление интегралов (11-50) может быть легко проведено пока лишь для газов, где плотность частиц мала. Прямое определение корреляционных функций жидкости с помощью статистической механики выполнить довольно трудно (см. ниже).

В Ы Ч И С Л Е Н И Е С Р Е Д Н И Х З Н А Ч Е Н И Й

Предположим, что нас интересует среднее значение некоторого

экстенсивного физического свойства жидкости Мп ъ г2, . . ., гп), определяемого взаимным расположением группы из п частиц. В раз­ личных областях жидкости n-конфигурации будут отличаться друг от друга, так что значения свойства Мп будут колебаться вокруг среднего значения. Поскольку вероятность dWn расположения частиц в я-конфигурации описывается корреляционной функцией

Fn (r\y rii

• •

•. rn)y т о

усреднение величины

Mn выполняется

инте­

грированием

с

весом

Fn:

 

 

 

 

 

Мп

= $Mn(rltr2,...,rn)dWn

=

 

=

l/V»

f Mn{ru

г г , r „ ) Fn (ги r2,...,

г„) duj dv2...dvn.

(11-53)

В жидкости, содержащей jV частиц, число различных групп из п. частиц равно Nli(N п)\п\. Поэтому значение свойства Мп для всей жидкости равно

X Fn (rl t r,,...,r„) dv1dv2,...,dvn.

(11-54)

Наиболее часто рассматривают свойства, зависящие либо от координаты одной частицы (аддитивные свойства), либо от взаимных расстояний пар частиц (свойства, определяемые парными взаимодей­ ствиями). В первом случае

Шм) = т \ м ^ й ^ И d v ^

( 1 Ь 5 5 )

>

Из уравнения (11-49) следует, что F1 (rt ) = V/N;Rx = V/N\R0 = 1, поэтому

M i W = 4 - W r ) < f o .

(11-56)

Для свойства, определяемого парными взаимодействиями, можно записать, что

М2 (ЛО =

N

( w l

) i М2(rur2)F2и

 

; а ) А^А,,.

(ІІ-57)

Теперь учтем, что согласно соотношениям

(11-49) и (11-41)

 

 

рг (/ь га ) =

д г ^ І ! )

^ 2 (/ь r2) = N ( ^ _ { )

 

R0R ( I г.

Гх I ) . (11-58)

Тогда

 

 

RJ2 \ М2 ( г ъ r 2 ) R ( | г2

 

 

 

 

 

 

М 2

(Л/) =

 

-

^ | ) dt^ do,.

 

(11-59)

 

 

 

->

->

 

 

 

->

-*

 

 

 

Если функция M2(rlt

 

гЛ

зависит, как

и

R (I г 2 г 2 1 ) ,

только

от

 

 

 

 

 

 

 

 

 

->

->

 

 

расстояния между

 

парой

атомов, то,

обозначая

2—гх|

через

г,

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

М2

(N) = RJ2 J dtox2

(г) R (г) dv =

 

 

 

 

 

СО

 

 

 

 

СО

 

 

 

 

 

= R0V/2

J М 2

(г) Я (г) 4лг2 dr = 2nN J М2 (г) R (г) г2

dr,

(11-60)

 

о

 

 

 

 

о

 

 

 

 

 

 

где R (г) — парная функция радиального распределения.

74

Важным примером применения уравнения (П-60) является вычис­ ление полной энергии жидкости в приближении, когда потенциальная энергия системы может быть представлена в виде суммы энергий парного взаимодействия частиц. Если парный потенциал обозначить через ф (г) и учесть, что средняя кинетическая энергия одноатомной частицы равна 3 / 2 kT, то для полной энергии получается

оо

U = -^-NkT

+ 2nN J ф (г) R (г) г2 dr.

(11-61)

 

о

 

Часто вводят безразмерную

функцию

 

Тогда

со

U =

^.NkT+

*gL$<p(r)g(r)r*dr.

(И-63)

У Р А В Н Е Н И Е

СОСТОЯНИЯ ЖИДКОСТИ

 

Для нахождения уравнения состояния жидкости вспомним, что давление р выражается через свободную энергию F системы следую­ щим образом:

(11-64)

В статистической механике Гиббса показано, что

 

 

F

= —kT

In Z,

 

(11-65)

где Z — интеграл

состояний

жидкости:

 

 

 

7 _

( 2nmkT

\ з л г / 2

і

г

Г

U

(7ltrt,...,rN)

X

 

Л2

)

T v T j e x p [

 

kf

X dv1dv2,...,dvN

=

^ —

j

-Щ-QN.

(11-66)

Величину QN называют конфигурационным интегралом. Тогда для давления справедливо уравнение

Р = КГ{У$*)Т.

(И-67)

Для вычисления производной d In QN/dV применим следующий искусственный прием. Умножим координаты всех частиц qi на неко­ торый параметр X. Так как объем системы пропорционален Я,3, то

dV/V = 3 dXIX.

(11-68)

75

Создавая изменение объема путем вариации параметра X, находим

 

p=r-kT

д In QN

 

_

kT

dQN(X)

 

(11-69)

 

 

 

 

 

 

 

 

 

 

 

 

dV

я=і

WQN

дХ

 

 

 

Далее с учетом

уравнения

(11-66) запишем

 

 

 

 

 

 

 

 

 

kT

dvy

dv2,...,dvN.

(11-70)

 

 

 

 

 

 

 

 

 

 

 

Примем, что потенциальная энергия представляет собой сумму

энергий

парных

взаимодействий,

т. е.

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

(II-71)

 

 

 

 

 

 

 

 

 

 

 

 

 

V

[ r i , r 2 , . . . , r N )

=

£ Ф ( к - — r k

 

 

тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U(X)=

2 < p ( b | r , - r * | ) .

 

 

(11-72)

Дифференцируем

QN (X) по X:

 

 

 

 

 

 

 

% ^ = 3 ^ в

л г - , | е х

р

Г

им

dvx

dv2,...,dvN-

 

 

 

 

 

 

 

 

kT

 

 

 

 

 

 

 

 

Kk

 

 

 

 

 

 

 

 

 

 

X

exp

 

U(X)

dv1

dv2,...,dvN

 

(11-73)

 

 

 

 

kT

 

 

 

 

 

 

 

 

 

 

 

 

(штрих

обозначает

дифференцирование

по

аргументу).

Разделив

уравнение (11-73) на QN (X) и приняв X =

1, находим

 

 

 

1

dQN(X)

 

 

N

 

 

 

 

 

 

 

 

 

 

X

 

 

QN

дХ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X ф' ( Р< г *|) " ^ Г е

х Р [ — w ]

d V l d

v * T ~ > d v

 

(II-74)

Множитель

1/Qyv : exp [— U/kT] dv±

dv2,...,dvN

 

 

 

 

 

 

представляет собой гиббсову вероятность данной

конфигурации

частиц. Следовательно, под

интегралом

усредняется

(по

ансамблю

Гиббса) величина \ г ( — r f t | ( p ' ( | r , — r k ) |, т. е. величина типа М 2 в уравнении (11-60). Поэтому вместо вычисления интеграла в выра­

жении (11-74)

можно автоматически

воспользоваться

результа-

 

 

 

 

->

 

 

 

том (11-60). Тогда,

обозначая |г,- — r k

\ =

г , находим

 

1

dQN

(X)

= 3N'

±г2яЫ\

q>' (r)R(r)r3dr.

(Ц-75)

QN

дХ

 

 

 

 

 

 

 

После подстановки уравнения (П-75) в (11-69) получим

NkT

оо

 

 

V (r)R(r)

rsdr,

(11-76)

Р = V

3 V

о

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

с в я з ь МАКРОСКОПИЧЕСКИХ СВОЙСТВ ж и д к о с т и

СФУНКЦИЕЙ РАДИАЛЬНОГО РАСПРЕДЕЛЕНИЯ

ИПОТЕНЦИАЛОМ МЕЖАТОМНОГО ВЗАИМОДЕЙСТВИЯ |

Полученные уравнения (11-61)

и (11-76) характеризуют

связь

между

энергией жидкости, давлением, парным

потенциалом

ф (г)

и функцией радиального распределения. Было бы весьма

заманчиво

выразить

через

функции ф (г) и

R (г) или их производные дру­

гие свойства жидкости, как термодинамические (энтропия,

свободная

энергия,

поверхностное натяжение

и т. п.), так и кинетические

(вяз­

кость, коэффициент самодиффузии

и т. п.). В общем случае, однако,

этого

сделать

нельзя. Наиболее

наглядным

примером

является

энтропия жидкости S. Последняя не является средним значением какой-либо величины, определяемой суммой свойств каждой ча­ стицы или суммой свойств пар атомов. Поэтому величина S не может быть выражена простым уравнением типа (11-61) или (11-76). Строго говоря, она определяется непосредственно через корреляционную функцию FN, заданную уравнением (11-52):

(II-77)

где 5о — энтропия идеального газа [т. е. нашей жидкости при усло-

->

->

->

вии U (rlt

г2 ,

. . ., rN) = 0], а усреднение проводится по ансамблю

Гиббса.

 

 

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

от

температуры,

из уравнения (11-63) можно найти теплоемкость

при

постоянном

объеме Су\

 

 

 

* = {

ж ) у = 4 N k +

2 п ] Т J Ф (г)

(И-78)

 

 

 

о

 

Далее учитываем, что

ТdS = dU + pdV

Т(dS)v = dU.

Отсюда

 

дії

\

c v

 

\ дТ ) v

T \дТ

v

T

 

dT

2n V

\ у

,f Ф (r ) dg(r)dT r*dr.

(11-79)

Однако уравнение (11-79)

требует

знания функции

dg (г)/дТ\

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

С радиальной функцией распределения оказываются непосред­ ственно связанными некоторые другие термодинамические свойства

жидкости, в частности модуль сдвига и сжимаемость.

Подробный

расчет

этих свойств

приведен

в упомянутой

выше

монографии

И. 3. Фишера

[5]. Например, модуль сдвига

жидкости

(высокоча­

стотный)

выражается

уравнением

 

 

 

 

 

 

 

ftp И

8 {г) dr.

(11-80)

 

 

 

 

дг

 

 

 

Этот

модуль

проявляется в

процессах

деформации

жидкости

с большими скоростями, когда

упругие напряжения не успевают

релаксировать. Соответственно, можно определить и адиабатическую (при постоянной энтропии) сжимаемость жидкости

(И-81)

P s ~ ~ Т

через функции g (г) и ф (г):

РТ1

N

 

дТ

 

 

 

 

dV Js

 

dr

 

 

 

 

 

X

1

ду(г

g{r) г* dr.

(11-82)

 

г г

дг

 

 

 

 

Если p s известна, то легко найти и изотермическую сжимаемость жидкости Вг :

Pr = -£-Ps.

(И-83)

Относительно просто вычисляется и поверхностное натяжение а жидкости в приближении, когда структура жидкости считается одинаковой в объеме и у поверхности раздела и температура далека от критической. В теории Борна — Грина для величины о получается следующее выражение:

со

а = т { т ) 1 ^ ( г ) г ^ г .

(ІІ-84)

о

Оно может быть несколько улучшено, если учесть изменение структуры жидкости в поверхностном слое [5].

У Р А В Н Е Н И Я , СВЯЗЫВАЮЩИЕ К О Р Р Е Л Я Ц И О Н Н Ы Е ФУНКЦИИ

Описанные выше результаты показывают, что знание корреля­ ционных функций Flt F2, . . ., Fn необходимо для расчета целого ряда термодинамических свойств жидкости. Ниже будет показано, что и кинетические свойства жидкостей в ряде случаев выражаются через интегралы, содержащие различные корреляционные функции. Поэтому целесообразно поставить задачу о вычислении корреляцион­ ных функций низших порядков через известные потенциалы межатом­

ного

взаимодействия

ф (г).

Эта задача решалась рядом ученых:

М. Борном, Г. Грином, Дж .

Кирквудом, Н. Н. Боголюбовым и др.

В результате были получены

довольно сложные интегро-дифферен-

циальные уравнения, связывающие функции Fn и /*„+ 1 .

 

 

Для

того чтобы

вывести

эти

уравнения,

введем

обозначение:

 

 

DN(rl,r2,...,rN)

=

 

U

и

 

гг,...,гы\

 

(11-85)

 

 

QN-1exp

 

kT

 

 

 

 

 

 

 

 

 

 

 

 

где

U — потенциальная энергия

системы

частиц;

 

 

 

QN

— конфигурационный

интеграл,

введенный

в

уравне­

 

 

нии (11-66).

 

 

 

 

 

 

 

Далее предположим, что потенциальная энергия может быть

представлена в виде суммы энергий парных взаимодействий

[урав­

нение (II-71)].

(N—s)

 

 

 

DN

 

 

 

Проинтегрируем

раз

величину

по

координатам

(s +

1)-й, (s + 2)-й,

. . ., N-й частиц. Тогда

в

силу

соотношения

(11-50) получится, что

 

 

 

 

 

 

 

F,{ri,rt

rs) = VS~N \... IFN[r1,r2,..,,rN)dvs+1dvs+2,...,dvN.

 

(П-86)