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

книги из ГПНТБ / Эккер, Г. Теория полностью ионизированной плазмы

.pdf
Скачиваний:
4
Добавлен:
22.10.2023
Размер:
12.01 Mб
Скачать

§ 3. Статистическая сумма для кулоновской системы

53

Обозначим через vsi

число прямых связей или

цепо­

чек нулевого порядка,

связывающих частицы i и s в диа­

грамме.

Для определенных выше величин выполняются сле­

дующие тривиальные

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

 

 

_ _

_

_

_

|

 

 

m + n=^N,

 

 

 

 

 

 

 

 

г

 

s

 

S

Ns= Ц (ms + ps), m =

 

Цг=5]м*г,

(3.65)

S

 

S

 

 

S

S

 

 

 

 

 

 

 

Р = S Рг-

 

 

 

 

 

 

 

г

 

Вклад

рассматриваемых

диаграмм в величины В^

обозначим

через

B ^ v). Здесь нижний индекс характери­

зует

структуру

прототипа,

а

верхний — набор

частиц

в цепочках диаграммы. Соответствующий вклад опреде­ лится выражением (3.51), в которое следует подставить разложение (3.58), используя введенные нами определе­ ния:

В»- . = _L_

(—l/0)v

 

vsr! ПЧ-х

 

(m’V) (m + р)!

 

X J П [ п (-■£*:)"“]

(3-вв)

is

Вэтом выражении употреблено сокращенное обозначение

xV* = j g {rhi) g (rl2) ... g Ы , - 1)ц.) g (п».) X

X ^ i ■■• r ^ d r hi .. . dr^.. (3.67)

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

надлежат набору т, обозначено в (3.66)

символом I

^г_

 

 

 

J

ТО*

Следуя

способу,

предложенному

Майером

(см.

стр. 49),

попытаемся

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

формулой

(3.52)

провести суммирование по наборам частиц N, выбранным

54 Гл. 1. Равновесные состояния кулоновской системы

надлежащим образом. Чтобы это можно было выполнить, необходимо, разумеется, сначала умножить величину

B't. на множитель nN-

(m,v)

__

__

Рассмотрим теперь

заданные наборы частиц N,

т, р

и рг. Сохраняя их неизменными, с помощью 1) переста­ новок частиц внутри набора рг„, 2) перестановок частиц внутри набора ms и 3) перестановок частиц между набора­ ми ms и ps можно образовать новые диаграммы, при­ надлежащие к группам, дающим одинаковый вклад

в величину В И .

(m,v)

Перестановки в пределах p*s соответствуют изменению порядка частиц в цепочках. Число диаграмм, образующих­ ся таким способом, определяется соотношением

Д р ;! •

<3-68)

Перестановки элементов в пределах ms здесь не рас­ сматриваются, поскольку соответствующие им новые диаграммы учтены при суммировании в выражении (3.51) как новые прототипы диаграмм.

Число новых диаграмм, возникающих от перестановок элементов между группами ms и ps, описываются соотно­ шением

(”*+ р)!

(3.69)

ml

 

Вклад всех диаграмм, образованных перестановками ча­ стиц в случаях 1) и 3), учитывается с помощью умножения

величины nN В -1 , определяемой формулой (3.66), на

(тм,v)

множители (3.68) и (3.69). Таким образом находим

[nNB»~ Л'

пт

(—l/0)v

х

1 (т , v)1

 

II с

X

f П

рг!

(3.70)

 

i

s

 

Чтобы определить вклад, представляемый всеми диаграм­ мами с одинаковыми длинами, но разной структурой

§ 3. Статистическая сумма для кулоновской системы

55

цепочек рг, просуммируем результат (3.70) по всем набо­

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

 

2

^ !П

[ - ( 4 n

/ e ) z | n

s ] M'»i

 

2

 

Psi!

 

 

 

^si-^i=cbnst

5

 

 

 

 

 

 

 

 

 

4 я

 

 

(3.71)

 

 

 

 

“0"

 

 

 

 

 

 

 

 

 

и

в результате получим

 

 

 

 

 

 

пт ( —1/©)VП

 

2 ? * J Ш

- кX^drЬ У-, 1'

 

 

т \

V ,,!

 

 

 

 

 

 

 

 

 

 

 

(3.72)

где хЬ определено формулой

 

 

 

 

 

 

,

■л 4 я

,

(3.73)

 

 

 

*Ь= 2 j i

r

z^ s-

 

 

 

 

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

няется наличием величины vsr! в знаменателе выраже­ ния (3.72). Это связано с тем обстоятельством, что незави­ симое суммирование для каждой цепочки дает излишне много состояний. При этом нам пришлось бы учитывать такие пары диаграмм, все отличие которых состоит лишь в том, что в них переставлены две цепочки, соединяющие одни и те же изображающие точки. Однако такого рода состояния неразличимы. Не входя подробно в детали соответствующих доказательств (см., например, [7]), мы йозволим себе привести лишь результат, сводящийся к тому^ что суммйрование по всем цепочкам от рг = 0 до р* = оо может быть проведено независимо при условии,

если в знаменатель выражения (3.72) вместо v4r! подста-

56

Гл.

1.

Равновесные состояния кулоновской системы

вить величину v?r!

— число прямых связей между части­

цами г и s данного прототипа группы.

 

Учитывая

это, получим выражение, описывающее

вклад

 

группы

диаграмм, соответствующей

данному

прототипу,

т. е.

 

 

[nNB*-

]

 

пт ( —l/0)v

 

 

ml

v",•! ПЧ’, } П х ,* г - 5 <р>.

(3.74)

(m, v)

 

где функция

Xi определяется выражением

 

 

 

.

xi(^) = S(-xb)'1|xi'“)(^)-

(3-75)

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

между спектральными компонентами x(f l)(l) и 8 Ш :

Х(/ ^ ) = [£ © ]^ +1.

(3.76)

С другой стороны, из выражения (3.75) следует

^ = 2 ( - х Ь ) йгх ^ .

(3-77)

нг

 

После подстановки соотношения (3.76) в (3.77) нетрудно выполнить суммирование, что дает

Хг

. , ~ 2

(3.78)

 

1 +

 

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

g (I) = J r*g (г)

dr ^ (а2 + S V , ■

(3.79)

откуда при а — 0 следует, что

 

X iH x b + 5 2)-1-

(3.80)

Наконец, обратное преобразование Фурье дает

 

Xi (r) = —

e x p ( - x Dr).

(3.81)

§ 3. Статистическая сумма для кулоновской системы

57

Это соотношение доказывает справедливость ранее выска­ занного 3-го утверждения.

Замечание по поводу 2-го утверждения. Доказательство

2-го утверждения нетрудно связать с выводами, содержа­ щимися в доказательстве 3-го утверждения. Рассмотрим диаграмму-петлю, состоящую из групп частиц as. Подоб­ ную диаграмму можно рассматривать как замкнутую цепочку. Числа а и as соответствуют использованным выше обозначениям рг и p,,s. Совершенно аналогично выраже­ нию (3.70) для вклада такой петли в функцию S получаем выражение

5к е) = Д | [-(^/0)И ",]‘Ц } - L х<«- 1) (г = 0). (3.82)

Как и в предыдущем случае, мы сначала поменяем поря­ док всех частиц, сохранив неизменным их набор о. Затем учтем все возможные наборы, сохранив постоянным только число частиц а. После этого просуммируем результат по всем возможным значениям а от 0 до оо. В результате получим

5<С)= 2 ^ - ( - > {Ь)а -5гХв-1(г = 0).

(3.83)

а 5=2

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

дифференцировать величину S iC) по хЬ> можно получить выражение

d S (с>

8^ 2 ( - * Ь Г х (о>(г = 0) =

а ^ 1

 

=

5 Г 1х(г = 0)_зс0(г = 0)1’ (3-84)

что позволяет непосредственно перейти к функции х (г). Имея в виду, что Х<0>есть функция g, в пределе при г —> О и а -*• 0 находим

dS(с) xD

58 Гл. 1. Равновесные состояния кулоновской системы

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

5(С):

''D

(3.86)

12л

 

 

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

получим хорошо известный закон

Дебая — Хюккеля:

=

п -

(3.87)

 

2 4 л

 

Рассмотренные петлевые

диаграммы

определяют вклад

в вириальное разложение, связанный с взаимодействием по закону Дебая — Хюккеля.

3.5. Разложение по большим группам

Предложенное Майером разложение по группам в g- представлении устраняет расходимость, связанную с даль­ ним взаимодействием. К сожалению, при этом возникает новая расходимость для областей близких взаимодействий, когда г -*■ 0. Это неудивительно, так как степенной ряд от функции ехр(—1/г) при г — 0 расходится.

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

Чтобы избежать расходимости, возникающей на корот­ ких расстояниях, Абё [8] предложил использовать разло­ жение по болыйим группам.

Для понимания этого метода достаточно рассмотреть лишь однокомпонентную систему (г? = е2). Полученное

§ 3. Статистическая сумма для кулоновской системы

59

ранее выражение для S можно тогда написать в виде

S = s{c)+ S(P) =

- s - ’+ s S J a n t

g2X (rik)

1

dri . ..

drm^\,

© J

v,.„!

 

 

 

 

(3.88)

где S <c>есть вклад от петлевых диаграмм, соответствующий дебаевскому взаимодействию, а — вклад от диаграммпрототипов. Первое, суммирование в выражении (3.88) проводится по всем диаграммам-прототипам с различными числами т, а второе распространяется на все диаграммы-

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

Простейшая диаграмма-прототип и ее вклад в S суть

< 0 :5 (m’v)= 5(2' 3) = ^ - J

О

 

оо

 

= —

j "7 ехР ( — Зхг>г) dr.

' (3.89)

 

о

 

Входящий в это выражение интеграл расходится.

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

Проиллюстрируем это на примере двуузловых диаграмм возрастающего порядка, представленных в виде следующе­

го ряда:

 

 

 

 

^

+ о

+

# . . . .

(3.90)

60 Гл. 1. Равновесные состояния кулоновской системы

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

Вклад двуузловых диаграмм (3.90) может быть записан в виде

оо со

2 - т ( - т Г - 4"г!* =

 

О v = 3

00

е%ч(г\

= 2яп2|{ е _ _ 0_ - 1 + ^ | ^ - - - | - [ ^ Ш ] 2} г2 dr (3.91)

о

 

или с учетом обозначений

~ ^ е2ъ ^ = 2 ( —1г У '^ г

(/>°) (3-92)

v=l

 

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

 

оо

 

S (2>= -^- ^ w3 (г)-4nr2dr,

(3.93)

о

 

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

Изучим теперь прототипы диаграмм с тремя узлами.

Рассмотрим основные их типы

 

 

(X ) (2 ),

( J )

(3.94)

От каждого из этих прототипов можно образовать осталь­ ные виды путем замены двойных связей в основных диагграммах на многократные. Проведем суммирование по всем этим дополнительным группам диаграмм и исходным четырем типам (3.94). В результате получим

Sl3) =

j w0(12) w2(23) w2(31) d r ^ +

 

 

+

j w2(12) w2(23) w2(31) dndr2,

(3.95)

§ 3. Статистическая сумма для кулоновской системы

61

где первый член соответствует вкладу первых трех диаграмм.

Обобщение описанной выше процедуры на прототипы диаграмм с числом узлов более трех очевидно.

Снова мы видим, что вклады могут быть пред­ ставлены как интегралы от произведений функций wt. Это наводит на мысль о том, что в качестве упорядочивающей схемы теперь можно ввести диаграммы в ш,-представлении. Однако в противоположность диаграммам в /- и g-пред- ставлении мы должны отличать разные типы связей шг между изображающими точками. Для этого мы используем следующие обозначения:

wQ

; щ

; ш2 —

:з “*■

■> • • • •

 

 

 

 

(3.96)

Применяя

данные

обозначения, можно

записать

вклады

5<2) и <S(3)

в виде

 

 

 

 

с

4яг2 dr,

 

 

g(3> _ п3 Г

]

dti dr2-

(3.97)

 

 

 

 

Искомая функция S снова представляет собой сумму по некоторым неприводимым диаграммам, выражающуюся через групповые интегралы в ^-представлении. В полной аналогии с методом Майера запишем выражение (3.43) для S в форме

s = s (c)+ 2 f c $ r " fc+1’

(3-98)

ft= 1

 

где мы определили в качестве неприводимого группового интеграла в ^-представлении величину

V* = JT J 2

П u’i Я П

(3-99)

В сумму под интегралом (3.99) входят лишь «неприводи­ мые шг диаграммы».

62 Тл. 1. Равновесные состояния кулоновской системц

Соотношение (3.99) при k = 1 и к = 2 сводится к фор­ мулам (3.97). Для к = 3 мы имеем

Г-*.]+О +П +0 +Ш +£ 3 ) *«***■

(3.100)

Абё использовал свою идею при рассмотрении электрон­ ной системы с размазанным ионным фоном. Он также вычислил член S l2>, чтобы оценить вклад членов высшего порядка.

Непосредственное вычисление дает формулу

S<2>

4)е2

(3.101)

4 я в

где константа А равна

 

 

 

 

А = 4 -С

_ 1_

In 3-

11

(3.102)

 

12

 

72

 

Здесь через С обозначена постоянная Эйлера. Порядок величины членов третьего приближения есть

5 '3’- 0 ( ( х сгго)5).

(3.103)

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

K3D

Sm (3.104) 12я

то увидим, что вклад от прототипов диаграмм второго порядка в (х0гш)-1 раз меньше, а от прототипов третьего порядка в (ко^и,)"2 раз меньше по сравнению со значением S iK\ соответствующим приближению Дебая — Хюккеля.

На фиг. 3 приведена величина отношения вкладов прототипов диаграмм и петлевых диаграмм в уравнение состояния:

А2

S™ — n(dS™/dn)

_ 0,92—1п(Хр/гщ)

(3.105)

Ас

S (c>п (dSK,/dn)

A.-olrw

 

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