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

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

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

(~peo-l)

, где

ßt(p)~

(00

, a t '- мини­

 

!-(^)\'Ц<тгг о1ъ

мальное значение

ф ( г )

(см..* например,, формулу (1.34')).

 

Бонроо об аналитических свойствах функций \

( ж) и

9

(ъ )

исследован, в

известных работах

t e a

и Ли

(1952

г.,)

по теории фазовых

переходов.

 

 

 

В заключение заметим,

что определение

неприводимых,

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

во й системы,. поскольку' она дает возможность выразить ßj, через А/ - частичные квантовые статистические

суммы

. При этом формула. (2.6 )

имеет

важный

т е о р е т и к о - в е р о я т н о с т н ы й

а н а ­

л о г :

в теории вероятностей она описывает

разложение

характеристической функции* и тогда

=

М ^ имеют

омысл

м о м е н т о в

порядка А/ * а

 

k t K у -

м у л я

н т о в, или

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

Т и л я . .

 

 

 

 

§3. Вириальноа разложение, аналитический метод

На основе разложений Майера (2.11

)и (2.12

) можно

построить разложение давления в ряд по степеням плотнос­

ти', или,, как

принято говорить*

в и р и а л ь н о е

р а з л о ж е н и е .

Для этого надо исключить иѳ

(2Д І)

а. (2„І2)_ активность

2

, т .е . выразить

из

(2.12

) Z

через 9

в виде, ряда

по степеням

9

и затем

под­

ставить

найдениый

ряд в ряд

(2 .1 1 ).

 

 

 

В результате

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

по степеням плотности*

которые

мы

запишем в

виде

 

 

 

 

 

 

 

 

 

и I I

 

 

 

'

у

 

=

f n . / г ) ?

.

(2.13)

Коэффициенты

 

п =о

 

в и р и а л ь н ы м и

^ п

 

называются

к о э ф ф и ц и е н т а м и ;

очевидно*

-

I .

 

Внриальныѳ коэфф'ицис .іты

/| Г)

связаны

с групповыми

интегралами

. Наиболее, простой

а н а л и т и ч е с -

50

к и й способ вычислении ^ (і через $ g основавна при-* менении теории вычетов и принадлежит Д„И-Зубареву; ([958 )

Рассмотрим производную ( 2 .ІУ) *

Ияо

 

 

( Й.1 4 )

 

По теорема Коши,

На плоскости комплексной переменной 9

Ч/ч- /С(Н/*9)сІ9

 

 

9

- 9

.

<*Л 5 )

- ё

----------

 

 

. Здесь контур, интегрирования лежит внутри круга сходимос­

ти ряда (3,13)

и охватывает начало

 

координат,, т.е.,

точ­

ку

9

 

 

 

 

 

 

 

Тѳперь от

комплексной

переменной: плотности р

в

(2,.13) перейдем к комплексной переменной активности

^ *

так

что

согласно ( 2 .12)

^

 

 

 

 

V

 

fz-J

с

(2 .1 2 ' )

Тогда, с

учетом (1 .4 7 ' ) (Ъ %/ è ф ) d ?

 

></>р-екФ[1-&>] т г -(2Л5)

Раскладывая подинтегральное выражение в (2,15)в рад по

степеням

9 , как геометрическую прогрессию* и сравни­

вая с(2.14),. полудим

_ t1

( п

+ | ) { Ң4, = (?fr0

ф ?С *Г ) d ^ / ^ , ( 2 .1 6 )

где контур интегрирования лежит внутри круга сходимооти

(2.12)

и охватывает

начало координат.

 

Положим

_ „

_

п п

Р^і

o ^ ) = £ ( 4 + S ) , S = Z S fc ,

 

 

. ( 2 . 1 7 )

Тогда

.

*

 

 

,

 

^

(I +■s)

 

(2Л6 *

и по теореме о вычетах интеграл (2.16 1 )равен коэффициен­

ту

при ^ п

в разложении функции [ / + S f e ) ] в ряд

по

степеням

^

51

Т е п е р ь в о с п о л ь з у е м с я ф о р м у л а м и б и н о м а Н ь ю т о н а :

+

s f

 

т - і

Гт + п - I

S і-п

(зла )

 

o '

 

 

 

 

 

 

 

І11--0

 

 

 

 

 

 

5

\' р

SI£'/

{.»„у

?

 

 

 

Л 9)

 

 

 

 

m l г.

 

 

 

 

H £i ,Подставим ( fe.-i9 )в

(2.18

выделим в

(2 .^коэф ф и ­

циент при

п

. В

результате

находим вириальный коэф­

фициент \

п + і

::

 

 

 

 

т

 

(п-|)!(пи)Г ! -'J jri' l'mtM-іІ'П

 

 

 

 

 

*

(ш ,|

.

 

t

 

 

 

/■>1

= </Ѵ, /7?^ }

 

 

 

=

^ •

( <1*21 )

.здесь

Л

о _

^

. t ^ І£—^

 

*

 

 

L*

—т

 

 

 

Итак,

формула ( 2.20 ) позволяет

найти коэффициенты

в раэлокении давления по степеням плотностиг если

извест­

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

актив­

ности~

 

 

 

 

 

интегрируя ( 2 Л 5 ) до

t?

с дополнительным услови­

ем ( { ! ? )

I

при g7 -»■

О.

получим связь I ( р ) и

Интегральное, представление <2Л&) позволяет оценить коэффициенты вирияльного разложения ( 2ЛЗ ') и его радиус сходимости,, исходя из оценки для радиуса сходимости ря­

дов

і<іаёера..

 

 

Такие оценки были получены в работах Лебовица и

Пенроуза CLSG4) и

могут быть сформулированы как теоре­

ма:

р а д и у с

с х о д и м о с т и в и р и а л ь н о -

г о

р а . з л о ж е н и я ( 2 . 1 t} я е ы <з н ь ш е, ч е м

52

oi

Р ( > ) н

J f ( z . ) 1 г 2 о( ъ ,

( 2.23 )

 

 

 

 

 

 

 

к-

h -/

 

 

 

П - I

 

 

[■ы

' ( е ^ %

1 ) р ,

( р ) _

( 2.23')

г д е

°< •=

0,28952.

 

 

 

 

 

Остальные

обозначения даны в конце §2* г л .II.

 

 

Достоинство аналитического

определения ( 2 .20) коэф­

фициентов

 

в том, что оно

справедливо

и в

к в а ­

н т о в о м

случае, когда

^

теряют смысл групповых

штегралЬв и либо определяются из разложения большой

статистической

суммы, либо выражаются через

/V-частич­

ные статистические, суммы согласно ( 2.10 ) (см.упражнение 2 к гл.11) .

Недостаток аналитического подхода состоит в том,, что при этом теряется возможность непосредственной г р а ­ ф и ч е с к о й интерпретации коэффициентов вириального разложения,, в отличие от коэффициентов, é . разложений

Майера- §4. Неприводимые интегралы. Метод мультипликативных

 

 

итераций

 

 

 

 

С иомощью определений и правил 3«. 4,

5-

из §2

этой

главы вычисление группового

интеграла.

в конечном,

счета сводится к вычислению вкладов, неприводимых диаг­

рамм ( правило

5),

 

 

 

 

 

 

Выделение

этих диаграмм и их вкладов, осуществляется

следующим образам.

t одинаковых,

 

 

Рассмотрим группу, из

но различи­

мых молекул. Ей соответствует некоторая нумерованная

связная и, вообще говоря, приводимая диаграмма с і

вер­

шинамиГ( имеющими различные индексы

(координаты). Среди

этих вершин, имеются узловые

точки

(узлы). и каждый уаел„

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

из §2, есть

точка сочленения нескольких

связных диаграмм.

Сначала, выберем такой узел

О , ,

от ко-

53

торого

можно отчленить

н е п р и в о д и м у ю диаг­

рамму,

например, о К,

+ I вершиной. Среди

них имеется

и вершина О г ',

 

 

В

интеграла по соответствующим к ( + I

переменных

совместим начало координат о О, , после чего подинтеграль ное( усредняемоа) выражение из произведения функций Май­ ера будет содержать К t индексов (кроме О , К которые не содержатся в оставшейся диаграмме^ группа)* По этим координатам можно интегрировать (усреднять) независимо от координат оставшейся группы* Результат такого интегри­

рования обозначим через

 

(< +j

 

 

 

 

 

 

 

 

 

 

 

 

 

(2*24)

 

 

 

 

к, +і

 

 

 

 

 

 

где два штриха означают неприводимую,

или

д в у с в я з-

a y e

диаграмму

с

k , +1 вершиной..

 

 

 

Теперь остается усреднить оставшуюся связную диаг­

рамму с

( £ ~

к ,)

вершиной,

среди

которых находится и.

бывший узел

О,

(для оставшейся диаграммы

О , может и

не быиь узлом). С атой диаграммой поступим так же,

как

и с исходной

снова выберем узел

О ,

отчленим от

него

связную диаграмму

.DLT"+(c

К г + I вершинами* выделим

ее в к л а д Q £ +, )>

 

1 и так далее,

цока не останется

самая последняя неприводимая диаграмма Т ) ^ / с к^+І

вер­

шинами. Она не имеет узлов

 

и поэтому содержит k |+

I

индекс. Таким образом,

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

равен

 

т

 

,

 

X

-f

 

 

 

 

<D

 

 

 

 

 

 

 

 

 

 

Среди чисел

 

k L

могут

быть одинаковые, т .е .

может

оказаться

несколько

двусвязных; диаграмм ( групп ) одинако­

вого размера.

 

Н ( двусвязных групп из 2

В общем случае

имеется

молекул,,

И к

двусвязных групп из к + І молекулн,

где К =1*2,

I - 1,

причем согласно (2.25

54

 

 

 

1 - 1

 

 

-t-i.

 

 

 

 

 

 

 

k n ,

 

 

 

( 2.26)

 

 

к

=

/

 

 

 

 

 

 

 

 

Тогда

 

I -

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Полный вклад всех таких, разбиений

в групповую сумму 5 ,

 

 

 

 

 

 

 

 

м

 

А-< .

\ п .

 

 

 

 

 

 

 

 

M f ] < s

к+і

Здесь

 

 

 

 

 

 

 

 

е1 к = I

 

О'k+i ~ °Умма различных/ топологических и нумеро—

ваннах)

неприводимых диаграмм. По аналогии о групповой

суммой

S i (см.§2 этой главы)

она называется

н ѳ п р и-

в о д и м о й г р у п п о в о й с у м м о й .

 

По правилу, граф-аналитического соответствия, из­

ложенному в

начале §2,

и в согласий. с( 2.24),

 

 

 

 

5:+, ь г <

 

1% > ,-

 

<»•»>

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

диаграммы с

к +1 вершиной.

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

efef.

г

к

,

 

 

-------/,

 

 

 

sl)*v

 

К’ ІѴ

 

(2.28 )

 

 

 

 

 

 

 

 

 

 

 

 

 

Величина

Рк

называется

н е п р и в о д и м ы м

и н т ѳ г р а д о и . Из оралвення

(2.8)и

(2.28)видно,

что

определение, неприводимого

интеграла. ^ +|аналогично

определению группового интеграла

& ' . При тех жа усло­

виях., что

и для

,

т .е . при малых значениях газового

параметра,

В

так же,

как

и

5» , на. аавися'г от

объема.

 

 

' 1<

 

г

5

в

(2.27)

^

 

 

 

Множитель М р

t

есть число различных

способов,

которыми

различимых молекул можно разбить

на

П (

пар,

П 2 троек,

. . . „

 

групп из

k

+1 молекулы,

и т . д . ,

причем каждая молекула, может входить

как только

в одну,

так

и о д н о

в р е м е н н о

в

н е с к о л ь -

55

к о г р у а в . Последней обстоятельство существенно ус­ ложняет комбинаторную задачу перечисления неприводимых групп по сравнению с аналогичной задачей для просто свя­ занных груди* рассмотренной в §2 гл.П. Эта трудная зада­ ча была решена Дж. Майером и Гаррисоном, которые нашли

 

{nj -

Ü

□ ( £

11КП,

(2.29 )

 

 

Іг

 

 

 

 

Подставляя (2..8 ).(2^28'

) и <2.29 )в

(2. 27),

получим

^

= F ^

l d

( l ^ ) V

n

K l

( 2 . 3 0 )

с ограничительным условием на П к (2

.26 V

 

мТеперь остается выразить коэффициенты

вириального -

разложения (2,13

) , )( п ,

через

неприводимые диаграммы.

Для этого

надо подставить

(2.30)

в ряды Майора и исклю­

чить z?

►Эту задачу можно решить алгебраически или с

помощью теоремы Лагранжа об обращении рядов (вывод Б,Ка~ ва, 1938 г . ) . Окончательно получается следующий замеча­ тельный результата

 

W i ~ ~ И 7 Г f V

, ( г -з п

Таким образом, к о э ф ф и ц и е н т ы в и р и ­

а л ь н о г о р а з л о ж е н и я

п р о п о р ц и о ­

н а л ь н ы

н е п р и в о д и м ы м и н т е г р а ­

л а м .

 

 

Следовательно, формула (2,30 )с (2,31 ) является обращением результата (2.201

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

Мы рассмотрим другой, бодеа простой метод вывода: вириальных коэффициентов (2,31 ) , следуя в бсновном Ван К а м п е н у ( 1961 г . )

Излагаемый ниже метод естественно назвать м е т о ­ д о м м у л ь т и п л и к а т и в н ы х и т е р а ц и й

56

 

Рассмотрим конфигурационный интеграл Q „ опреде­

ленный согласно

(2*1 ),. и обозначим

через Q fn>,V значение..

О ы

в том приближении, когда

т о ч н о

учитывается

одновременное; взаимодействие ( столкновение. ) двух,,

трех,

....

молекул до

П

включительно* Очевидно* Q ^

соот­

ветствует идеальному газу* а Q ^ =

Q

- точное зна­

чение*

 

 

 

 

 

 

 

Поскольку мы ищем разложение термодинамических фун-

■кций в виде ряда

(суммы) по степеням, плотности

и по-

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

Ul Qn„ то естественно о самого

начала искать Q _ в

вида

н е

с у м м ы , а

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

 

е г о при-

блаженных значений*

Q«'

0 =

Q " Q N

N

L Q N

n = 2

 

 

О ^ 2

N

п <">

 

_rn Ww

 

1

^(п-О (2,32)

 

 

/_ (п-І)1

 

N

 

(2 .32')

^

 

N

 

 

 

При вычислении Q^ все перечисляемые конфигурации каждый раз распадаются на группы из П молекул, и эти группы, по смыслу приближения не взаимодействуют друг с другом* Иными словами,, переменные интегрирования рас­

падаются на всевозможные сочетания из

А/ по П

слу­

чайных. величин* и эти

П - сочетания статистически, не-

зависимы* Общее число

таких: сочетаний

/АЛ

Ы>

равно [ п )=г7|7іѵ-п)!

Следовательно* среднее (2*1 )равно произведению средйих

для отдельных П - групп, а

среднее для каждой И - груп­

пы равно точному значению

Q n *

 

Поэтому О Г /с Г - R ^ п )

(2.33)

 

п

8-896

57

где

 

 

R n = Q n / Q ^ °

( 2 . 3 4 )

Теперь подставим (2.33 ) в ( 2.32г) и совершим термо­

динамический предельный переход N -»<*> ,\Л*оо

,/Ѵ/Ѵ= 9

конечно.

 

 

Так как при A/s. П , ' (

П ! ,

то Е .321 )

можно представить в виде

 

 

СО

( 2 .35)

fnQN

н -и

 

Наосматриваемая система термодинамически аддитивна. поэ­

тому R

должно иметь вид

 

 

R 1

(ѵ,т)= |+ V _ I1R

" ( т ) ( 2 . з б )

П +

 

П+І

 

 

 

Тогда для En Q I^ получается разложение по степеням

плоітности.

оо

 

 

N

" п Л ' ( Г

^

^ - 37)

Остается вычислить коэффициенты

р>

для этого paa-

ложения __

 

И И

 

 

 

 

R

=■Him V (R

“ l).

(2. 38)

п +І Ѵ~* а?

и + |

Для нахождения этих коэффициентов существует опре­

деление ( 2.34 ) . в котором знаменатель подсчитывается

по рекуррентной формула

л

(

П

\

(к)_

Г Ѵ

ГЛ

I\ k( 2.33' )

0 .m

о“ КѴі

 

 

 

‘ m

 

 

 

с граничным условием Г у

= I .

 

 

 

Теперь ааймемся последаательным вычислением R п+і результаты вычислений будем излагать на языке диаграмм.

Среднее от диаграммы, или вклад,, вносимый ею по правилу граф-аналитического соответствия, сформулирован­ ному в начале. §2„ будем обозначать графическим символом этой диаграммы, как это сделано в роотношениях (2. 2 ) - (2.5' )из §. 2.

58

При этом часть вклада диаграммы, нѳ зависящую от объема, будем обозначать графическим символом этой диа­ граммы, заключенным в фигурные скобки. Короче,

где t ) n -

{T):}tYXK>v"X^

символ топологической

неприводимой диаграм­

мы с

п

вершинами.

 

 

 

 

 

 

 

 

По правилу 5 из §2

все

в

конечном'ечетѳ сводится

к таким диаграммам. Имея в виду сказанное, приступим к

вычислению R' ,

 

 

 

 

 

 

 

 

 

h =

 

R2=QJQ1° =0

 

+ (

- )

( 2 . ^ )

 

 

 

 

 

 

 

 

 

_

I

 

 

(2.40)

 

 

 

,^(2.)

' 'e

 

I

 

J

 

 

 

R 3 “

^V3/ ^

2.

*

выражение для

\ я а,

уже получено в

формуле

(2.41)

;

Q

находим из

Q ;

 

гдагп =3,

(2.331),

о

:

R 2 = I + з

н

+

з

с— )

+ с— -)

(2.41)

D

 

_

1 + 3 ( ~ ) l A f c L l A . ,

(2.42)

 

з

 

і - < - З М + З С - ) 2+С—

 

 

 

Отсюда видно,, что кроме членов,

общих для числителя и

знаменателя,, числитель содержит только

н е п р и в о д и ­

м у ю

диаграмму,

а

знаменатель

-

малую более высокого

порядка (—)ä ~

V -

 

,

которая

согласно (2.У8)L н ад а ­

ет вклада в

Ң

 

 

 

 

 

 

 

 

 

 

Поэтому

'

R

=

( А

)

 

 

 

 

 

 

 

 

 

(2.43)

 

 

 

 

 

 

 

 

 

 

 

 

 

И = 3

 

R/j

=

Q ц /

 

. Числитель выписан в (2 .5 ; )»

зн ам ен ательн а)

 

находим из

(2.33'),

полагая m ' =4г

k = 3 t

 

4 ^

( 3 )

П

С2)

 

Ч

 

 

(2.44)

 

 

 

 

W /,

=

W i,

п

3

 

 

Ь9

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