книги из ГПНТБ / Физико-химические основы металлургических процессов
..pdfчение 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) |
||||
Аналогично сказанному выше, можно ввести вероятность dw1§ |
2 |
„ |
||||||||
того, что |
в |
объемах |
dVlt |
dV%, |
dVn |
окажется по одной частице. |
||||
Определим |
|
dwh 2 |
п следующим |
образом |
|
|
|
|||
|
|
dwht |
„ = Rn(ru |
r! f ..., гя) |
dV^Vz,dVn. |
(II-42) |
||||
Введенная уравнением (11-42) функция координат всех п частиц |
Rn |
|||||||||
называется |
к о р р е л я ц и о н н о й |
ф у н к ц и е й |
n-г о |
п о - |
||||||
р я д к а . |
Как и раньше, при увеличении расстояния |
| rt—rk \ между |
каждыми двумя объемами вероятности нахождения частиц в каждом
|
|
-> -» |
|
объеме dVt |
становятся независимыми, т. е. при \ rt—rk |
\ —» 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 dtox\М2 |
(г) 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 |
2л |
оо |
|
|
|
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) |