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

книги из ГПНТБ / Магалинский, В. Б. Методы статистической теории равновесных состояний

.pdf
Скачиваний:
12
Добавлен:
19.10.2023
Размер:
7.21 Mб
Скачать

Поэтому ниже будут изложены основные результаты, которые относятся к этой области, и будут использованы в дальнейшем для конкретного расчета F.C.. с короткодейстеием..

Радиальная функция как Р.С.,, так и С.Т.С.. имеет, как показывают расчеты методами гл.Ш, а также опыты до рассе­

янию рентгеновских лучей жидкостями (см.гл.Ш, §1)

харак­

тер

затухающих осцилляций,

а их максимумы ( пики) описыва­

ют последовательность

так

называемых к о о р д и н а ­

ц и о н н ы х

с л о е в . .

При расчете по формулам (4.18),

(4 .1 9 )

’•з предыдущего

параграфа основной вклад в

интегра

лы дает

п е р в ы й

к о о р д и н а ц и о н н ы й

 

с л о й

(первый

пик),

поскольку последующие пики, во-пер­

вых, сами по себе убывают

по величине, а во-вторых, силь­

но подавляются быстро убывающей с расстоянием функцией

Ф

('£.') • Поэтому пас в первую очередь будет интересо­

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

см. §8, г л . І )

. Первый координационный слой, таким обра­

зом * обладает

развитыми коллективными свойствами и устой­

чивой структурой,

а описывающая его радиальная функция

с достаточно

хорошо удовлетворяет условиям теоремы 2

кз предыдущего

параграфа*

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

В &8 главы I было показано* что для точно решаемой одномерной С.І.С.- этот формфактор зависит лишь,от одной

переменной (см.

формулы (1.61) „ ( і.б іО

и (1.63)) . По

этой причине, а

также из упомянутых соображений простоты

попытаемся представить этот формфактор

и в общем ’тчае

cj измерений в

виде функции от одной переменной

120

- ß . ( c d - c o 0) .

 

 

 

 

 

 

 

 

 

 

Тогда

относительная

плотность (см .гл .Ill,

§I

„ а

так­

же г л .I, §8 ) для первого координационного

слоя

,ТТ

 

бу­

дет иметь вид

 

 

 

 

 

 

 

 

 

 

П , ( ~ с ,I T ) = A l d j ß A 0 - 6 0 « ) ] , T D ( o ) = I ,

(4.27)

 

 

 

 

 

 

 

 

 

 

 

 

где

^

 

~

( rA

- объем d

-мерного

шара диаметра-^

С0о ~

Ы 1'(ао)>Л(ѵ), B 6 A

J TD C^D-Фуиквдиг ■подлежащие

определению..

 

 

 

 

 

 

 

 

 

 

 

Без

ограничения общности принято

Т ) (О ) = I .

 

 

 

Нормируем

П

на первое координационное

число.

 

которое,

как упоминалось

ваше,, достаточно

хорошо определе­

но для плотного

вещества::

м

 

 

 

 

 

 

 

П ^ (г, ігсіы _ ((~с) -

А \ТЗ [В•(А - %)]clw=?l(4,2e)

о

 

 

с (1.61-О из

со0

§8, г л . і ) .

 

 

 

 

(сравните

 

 

 

 

 

 

В предела плотной упаковки сфер

Z.J

так

называе­

мое максимальное

контактное число, в, частности

2 =

 

2 ,

=

б

,

"2^ ~

И3 условия

(4.28 ) следует

(Д> — Д .

ZФункция А (гг) находится из уравнения состояния

С.Т.С-, которое по теореме вириала ІОтузиуса для С.Т.С.. имеет такой вид (см, упражнение Ш.2 к гл.Ш,, где для С.Т.С. следует положить U = О) :

 

Ps O

^ V

T

*

O r t =

 

 

=

т г с І 5

( ѵ

Ѵ

^

=

( c 0 o/ 2 ) r T s

( a 0, ^ . C * - 29).

Из (4.28) и

(4.29)

 

 

2 .^С ѵ)/со0

;

Л

А

(-u ')

=

В ( т О

=■

 

 

 

 

 

 

 

 

 

 

(4.30)

Н, — 6 - ( б 0 - 0 } 0 ) =

 

[ ( А У ^ 0 ) “

 

16-896

.I2I

 

Из

(4 .2 7 )

при В

 

А и из

(4.30')

уравнение сос­

тояния С Л'.С,, можно представить

в наглядной форме

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.31)

где

 

 

 

схэ

 

 

 

 

оо

 

 

 

. СО ф г) = г " '

^C O iT dcO , I

- = z “ 1

 

 

,

(4.32)

т.е -

 

‘а

С0 о

 

 

 

о

свободного объема,,

СО

 

имеет

смысл эффективного

ограниченного первым координационным слоем.

 

 

 

Вблизи плотной упаковки строгай статистико-геометри­

ческий анализ,

выполненный

 

З а л ь ц б у р г о м

J

, В у д о м г приводит к уравнению

 

 

 

 

 

 

 

 

 

( ѵ ) ■- сі

( v / ' l Г ,

-

,

 

(4.33)

где

d

- размерность пространства, а ЯУі - минимальный

удельный

объем ; в частности ТГ =Qo ,1Г2 - Q c 'fë /2

/у2.

При

d = I

(4.33)переходит

в

известное

уравнение

Тонкое

(1.6.4) „ справедливое при всеі 15' . Напомним, что

это

уравнение

было пслучено. в г л ,І,

§8 точным расчетом одно­

мерного статистического

интеграла,.

 

 

 

 

 

Сравнивая (4.31) и

(4 .3 3 ),получаем

п

 

‘■со

= Cv/U|)a)o

,

^

= 2d.

 

(4.34)

 

 

 

 

 

 

 

 

 

 

 

 

 

Итак* структура первого координационного слоя может быть описана функцией от одной переменнойТ) (е ) ,0 * ^ < со { 'удовлетворяющей условиям!

Х>(%-)>0, t ) '( V ) S Ö

, 1 3 ( 0 ) = ! ;

(4*35)

оо

•J»

 

■ 2^ ,zj Wt] (^dd^ -2d. (4.36)

122

 

Как.следует уже из элементарного геометрического

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

D

( О

 

и еа интегралов, при выполнении

условий (4 .3 5 )

, (4 ,3 5 )

с

учетом (4 .34)

должно

быть

 

 

 

 

^

=

2 d

^

- 2 ,

/ 2

,

 

 

 

 

(4.34')

в ином случае

нарушается второе

иэ условий

(4 .35),

т.ѳ»

должно быть

t D Y ? )

=

О

>

 

 

 

 

 

 

 

.'

 

 

 

 

 

 

 

О

 

’s=*

 

 

 

(4-37>

 

t)0?)=j''~*

 

 

 

Неравенство

(4Д 4 Оперѳстяет выполняться при

3

* в

частности

при

d = 3

, ^

= I

, оно

становится равенст­

вом. Согласно

(4.37),

это

означает,

что

п р и

d è- 3

(в ч а с т н о с т и , в р е а л ь н о м

т р е х м е р -

н"о

м

с' л у ч а ѳ ) р а с п р е д е л е - н и е

о т н о ­

с и т е л ь н о й

п л о т н о с т и

в п, е р в о м

к о о р д и н а ц и о н н о м с л о е в б л и з и

 

ц л о т н о й

у п а к о в к и

 

с т а н о в и т с я о д-

н о р о д в ы м .

 

 

 

 

 

 

 

 

 

 

 

 

Распределение в случае произвольной размерности мо-

,жет быть описано простой интерполяционный формулой* кото­

рая

удовлетворяет

(4.35)

, (4. 36)

£

^

 

 

 

 

 

Г К О

=■[

I

-

( ^ / о

) ]

*,

 

 

 

V =

(.4о( -

 

) +■ 1

 

 

 

 

 

 

(4 .3 8 ')

 

 

 

 

 

 

 

 

 

 

 

 

 

2z,-Цcl

 

 

 

 

 

 

 

При

d

-

1,

 

-

Я.

 

из

Н.Э8) иолучается точный ре­

зультат X] ( О

= е х р (- ^ /£ ) .

полученный в

§8, гл Д

для

одномерной модели (формула

(1 .65)),

а при d

3

~ предель­

ная форма (4,37).

 

 

 

 

 

 

 

 

 

 

 

 

Распределение (4,38)

наглядно

показывает,

как

"раз­

мывается” первый координационный слой с увеличением раз­

мерности

и координационного числа

1?

 

 

 

 

123

 

Для того

чтобы использовать (4 . 27)

для расчета. Р .С .,

необходимо, преобразовать параметры (.1.27') к переменным

X =

C O / U - = 4 F Ra / a

и

е»

 

в

соответствии с тео­

ремой I об оптимальном моделировании,

т . епредставить

П

 

в виде,, аналогичном

 

(4.17)

:

 

 

 

 

®

____

 

 

_

I

~ Г

 

 

 

 

 

 

 

п

( ъ , ' Ѵ ) = - Ѵ

 

Г_

( Х . б О .

(4.39)

Для этого

S

 

 

 

S

 

 

 

- I

и переменную

выделим

в (4 .2 7 ) множитель V

 

X

в

аргументе

Tg = Д . (сО —(-Оо)

 

после

чего оставшую­

ся зависимость

от удельного

объема

, ЯГ

преобразуем в за­

висимость от конфигурационной энтропии

 

подстановкой

ЯГ

=

ЯУ^ ( б

)

,, где

 

(б")

определяется уравнением,

состояния

С.Т.С*. согласно (4. 29) .

 

 

 

 

 

Тогда о учетом (4.30)

 

 

 

 

 

 

л

" (x,s).= Г

$ 4

 

’ m

- 2IЦ .

( s i '□ f e

a ) 0 ' s

" s

 

(R.sl -

У (4.27')

Я ( « Я X s LÜ.

где

ЯГ, =

K e ®

) ]

,

Например,

для уравнения

 

состояния

(4 .3 3 ) из

(4.29) получается

 

 

 

я г

 

 

г,

 

 

 

- - I n

Г - 6 /d

1

- |

( s ) = г г ( ! І - е х р ( б / с Г ) _

'facile

 

'Теперь,

когда

 

! (п

и все ее

параметры определены,

 

термодинамические функции могут быть вычислены по схеме,

 

изложенной

в предыдущем параграфе, т. е. по формулам

 

(4 .18) „

(4,20) .

(4.21) .

 

вместо X, переменную •

 

-

Если

ввести,

как в

(4.33)„

 

X = Дтгі.3/3 'іГ

, то

основная для расчета

формула (4 ,1 8 )

 

примет вид

 

Сг о

 

 

___

 

 

 

 

U (ѵ,еУ- ( І / 2 ) \ Ф ( Ѵ Х )

( Х . е Ы Х ,

еле)

 

__

 

 

0

 

 

 

 

 

 

 

где

1~ ( Х , б )

определено

формулами

(4.27)

и (4 .38),

 

124

Потенциал' взаимодействия выберем ь общепринптой форме, по­

тенциала Ленарда-Джонса с

И =

6 , ^ = 9

(см. гл Л ,. § 5 ,

( 1 .34')). В новых, обозначениях

_

 

 

 

 

 

Ф (со-] -

8 о [ 2

(0>т /" У

* -

3

 

/<У> ] ,

(4.10)

где

hn

= ^Г/Г7 ^

/ З

 

^ гп -

положение минимума потен-

 

 

 

 

т '

 

 

м1

 

 

 

 

г

\

циальной ямы.

g

-

 

ее глубина. Теперь. подставим

(4.40),

где. со

- Ч . Г, /

j

в

(4Л 0'),;гдѳ

для с] =- Э>,

Г& определено,

в ( 4 . 27' ),

(4.37)

и

(4 .3 3 ').

Используя

(4.20)

и

(4.21),

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

форме;..

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

U' =

6

$ ' ( 2 ^ -

3 t

) ,

 

 

 

 

(4.41)

 

р /

=■

З

б

р

/ г

( t

~ І )

I

 

 

 

 

(4.42)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4 .4 $

где. подлежащий исключению параметр

t

равен.

 

 

 

 

t

 

= ' < p ' l ^ l - e x p f ß

e / â

)

] .

 

 

(4.44)

В формулах

(4.41) -• (4.43)

- Ll\

р ',~ Г /

-

п р и в е д е н ­

н ы е

 

переменные

(см.. гл.1„ §5) :

 

 

 

 

 

и '= и /е0 ,Т =Т/•£„,§'= Я/гг, р'. рѵ /е„<j.«s

'^1= 7

ч

/ѵ 5 - удельный объем для гранецентрированной куби-

°

 

а

 

ческе .л

решетки, с периодом £ m

 

 

 

 

Уравнения состоят- : (4.41)

-

(4.44) получены мето­

дом оптимального однолараметричеокого моделирования. Они:

>

просты по форме и наглядны физически,, так как в этих:

уравнениях явно выделено влияние

 

.:К о я Д е к т и в н ы х

свойств

статистико-геометрического

происхождения ( они. оп­

 

ределяют зависимость,

параметра

t

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

 

энтропет) и влияние:

индивидуальных,

свойств,, связанных

 

с видом потенциала взаимодействия (ѳти свойства, определя­ ют зависимость энергии от плотности ) . Расчеты, выполнен­ ные на основа (4.41 - 4. 44) численными методами,, ока-

125

аадись. в хорошем согласии с

экспериментом в области плот­

ного

газового и конденсированного состояний. Здесь мн рас­

смотрим лишь два результата:.

 

 

I , П а р а м е т р ы

к р и т и ч е с к о й

т о ч к и .

Для

п р и в е д е н н ы х

значений критической плотнос­

ти

Q 1 температуры

7 ^ и

п о л н о г о

давления

р>' -

^ І~Г> V- p '

получается следующие

приближен­

ный.

соотношения:

 

 

_

9' - О/а) -ТД'Д7),Р' - 6(9' ?.

бтс' + 5 Р С' »

loscs'f.

J (4‘46^

Корня.©той системы приближенно равны

 

9 ' = 0,3-Л'е = 1,^ »

0 , 12

Q3l.fc.46')

Соответствующие экспериментальные значения (для благород­ ных газеф ) составляют

9 с' - 0 , 3 2 ГГС' = l , 2 S ; P c' = O ,l2 ,J E > y f t ;-i;'-0 ,2 9 |i.4 6 9 Сравнение ( 4 .4 6 7 и (4<46") говорит о количественном

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

"улучшенное уравнение Ван дер Ьаальса" и результаты наи­ более родственных излагаемому ^методу решеточных и гчрочных

о теорий,- Более подробное сопоставление результатов приведе­

но в таблице.

0

 

 

2,. Л и н и я

Б о й л я

и

п р а в и л о " е д и н и ч-

н о й о ж и м а е м о с т и"

Л а н и я

Б о й л я

определяется равенством нулю

йшфигурздионного давления

р

. При этом полное давление

равно его геэокинетическому значению,

и так называемая

с ж и м а е м о с т ь і ? г г / 1 ,

, ^4-47)

І?ІГ /~Г= I .

126

 

Иожтоыу в литературе, линию Бойля называют

еще

л и ­

н и е й е д и н и ч н о й

с ж и м а е м о с т и ,

 

 

Иэ(4.42) при

 

р '

=

О ,

 

'(:==)

 

, а из (4.43)

на линии Бойля температура^

давление связаны л и н е й ­

н ы м

соотношением

 

 

 

 

 

 

 

 

 

 

 

 

Т

'

о<

Ц

(I

-

д 1) .

 

'

(4.48')

 

Подобная линейная зависимость недавно была обнаруже­

на на опыте ( I

о л л

е

р а щ

 

1967 )

для большей группы

веществ

и в широком интервале, состояний

( от

плотности.

р ,=-0 (т

о ч к а

 

Б о

й л

я ) до

почти удвоенного

значения

критической

плотности). Эта линейная зависимость

получила

название' п р а в и л а

 

е д и н и ч н о й

с ж и м а е ­

м о с т и .

Интерполяция опытных данных к

"Т/ = О

 

приводит

к плотности

 

 

I в

согласии

с (4.48). ». Приф О

линия Бойля закапчивается точкой Бойля.

 

 

 

 

Экспериментальное

значение приведенной теі&ературы

Бойля равно

3,4

тогда как из

(4,48)

при

= О получа­

ется

значение

4,0.

Это расхождение естественно .посколь­

ку метод оптимального

моделирования с помощью М.С. твер­

дых сфер наиболее

эффективен,

 

как указывалось выше,

при

значительных плотностях вещества.. Скорее,, результат (<*.48) следует рассматривать как теоретическое обоснованна пра­ вила единичной ожимаеыости при больших плотностях, и низ­

ких

температурах, т. е. в области конденсированных ‘состоя­

ний,

где пока отсутствуют результаты

прямых измерений,

'

Результаты выполненных расчетов

можно дополнить и

уточнить как за счет использования более полной информа­ ции о радиальной функции С.-Т.С,, так и за счет использо­ вания более реалистической но, очевидно,, более сложной модели с "мягким" отталкиванием»

о

о

127

§,4. Одноиаіжметоическиѳ моделироваішо кулоновской системы

г. .. Квантовая плазма как цростал система.

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

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

Оказывается, что в этих условиях существует область состояний, когда плазма может рассматриваться как систе­ ма,. простая в смысла определений I и 2 , данных в §1 этой гла_)Ы. Тем самым становится возможным ее описание с іюмощьг) теорем об оптимальном моделировании §2 этой главы ••

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

ри а л а .

Вэтом приближении квантовая статистическая сумма равна классическому статистическому интегралу с эффектив­ ный попарно аддитивным взаимодействием,, учитывающим двух­ частичные квантовые, корреляция. Тем самым удовлетворяется требование I в определении простой системы.

Эффективный потенциал взаимодействия зарядов <2 й-

Sg , или п с ѳ в д о п о т е н ц и а . л ф g(x,T),зави­ сит от температуры и удовлетворяет условиям

Ф „ (о,т)< °о, Фр(х;т)->ѳ е , / £ , z - * ..49)

аБ

йб

« ь

128

В дальнейшем будет показано, что можно полностью удовлетворить требованиям I и 2 в определении простой сис­ тема и вместе о тем сохранить существенные свойства ис­ тинного псевдопотенциала, если аппроксимировать нсевдопотѳнциа11 формулой

где

С ( Г )

-

радиус

квантовых корреляций„ Ч

? С ^ О I / С

при

~Cf

Qe?

*- Не ограничивая общности,, положим Ф(Ъ)=І ■

Тогда

 

 

р

 

 

 

 

 

 

I0.. При

Ъ

>~> і

(4,50)

асимптотически справедливо►

2°. При

Z

«

t , I =-■%,

г а е X

= Ѣ (2 iT(U T)~

- _

 

длина

тепловой

де-бройлевской

волны (ji(= m a mg (ri1a-i- т ф і

 

приведенная масса

р а з н о и м е н н ы х

зарядов)

(4.50) переходит в

и н т е р п о л я ц и о н н у ю

 

ф о р м у л у

д е

В и т т а , Справедливость

«той форму­

лы

с хорошей точностью и в широком интервале

температур

 

0 , А

1 о

( І о

=

 

-

потенциал

 

ионизации в вакууме) установлена численным расчетом нсев-

допотенциала

ф -

(О, "Г) в работе Д е в и с а и С т о-

р е р а

(Г968).а _

 

т-

3°.. При

1 «

1

0

симметрия 4rag относительно переста­

новки 0

( Qg

нарушается,, однако 'в этом случае

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

зарядами, когда

2

?_

 

Ф _ ( 0'Т ) £ Ф

 

<4-ш

При этом взаимрдейотвиа одноименных зарядов дает от­ носительно малый вклад в конфигурационный интеграл.

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

17-896

129

 

Соседние файлы в папке книги из ГПНТБ