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

книги из ГПНТБ / Кудрин, Л. П. Статистическая физика плазмы

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

вычисляются

с помощью квантовомеханических групповых

сумм Si.

 

осуществлена

перестройка

членов теории

Таким образом,

возмущений

и для

ZN получено

выражение,

удобное для вы­

числений.

теперь формулу (15.23) на zl

и

просуммируем

Умножим

полученное выражение по / от 1 до оо. Тогда слева получится статистическая сумма для большого канонического ансамбля,

а справа — разложение функции exp |( ^00 zlVbi}. Отсюда термо­ динамический потенциал получаем в виде:

~

ОО

(15.25)

— PQ = PFP = In Z = ^ zlVblt

где активность z выражается через

число частиц N

по фор­

муле

 

(15.26)

N = z(d lnZN:'dz)V'p

Суммирование диаграмм для кулоновской системы частиц

Для системы заряженных частиц отдельные групповые ин­ тегралы расходятся, поэтому для получения конечных резуль­ татов необходимо, вообще говоря, просуммировать весь ряд (15.25), что в общем случае неосуществимо. Возможно, однако, из каждого группового интеграла 6; ( /= 1, 2, ...) выделить определенные частные интегралы, число которых бесконечно, и вычислить их сумму в явном виде, получив при этом конечное выражение.

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

Рассмотрим графики для групповых сумм Si и выделим в каждой из них диаграммы с наименьшим числом линий взаи­ модействия, что соответствует выделению частных групповых

сумм,

пропорциональных низшим степеням е2. График (15.20)

дает

вклад в

Qo — термодинамический потенциал идеального

газа.

Первый

из графиков (15.21) определяет в групповом ин­

170

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

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

что и любая сложная

диаграмма, содержащая элемент ~ t (т.

е. график, в котором

одна из вертикальных линий имеет лишь один узел), также не дает вклада в Q. По этой причине в групповой сумме Бц к диа­

граммой наинизшего

порядка,

 

дающей

ненулевой вклад в Q,

является кольцевая диаграмма

(см. рис.

16, б).

сумм.

Аналогично

обстоит дело

и

для

других групповых

Следовательно,

рассматривая

лишь

кольцевые диаграммы для

S i , можно написать выражение

 

 

 

 

S[ (ф, р; q'', Р') = ( - 1)'

 

 

 

(Чг, Р; q'', Р'),

(15.27)

где GJ'[)— кольцевой

интеграл

I-го

порядка, соответствующий

диаграмме

 

 

 

 

 

 

 

(15.28)

Множитель (—1)' в формуле (15.27) учитывает знакопеременность ряда теории возмущений, а (/—1) ! — число способов построения кольца (число различных перестановок вертикаль­ ных линий); коэффициент 1/2 возникает вследствие необходи­ мости интегрировать по всем ориентациям кольца, так что в любой точке кольца правое и левое направления эквивалент­ ны. Подставляя формулу (15.27) в выражение (15.26), получаем

=

- Ы М ZliE. f

 

(ф, р; ф, 0) dql.

(15.29)

 

 

/Vl-

v

 

 

 

 

Можно написать

также общее

выражение

для

кольцевого

интеграла:

 

 

 

 

 

 

 

G\l) (ф, р; ф ',

РО =

[ J'1 .

.

. f/_1 G(qb

р; хъ p j X

 

 

 

Ф р-

 

Р'

 

 

X G(xr, Р;

х ;,

р;) G (х,', р,;

qj , р') Ф ( | х, — х2

| ) X

Х С (хг, Р/—1; хг, рг) G (х ;, р,;

q;, P')G(q„ Р; х„ р/_,) X

X Ф ( I х' — xj I )

П

G(q,-, Р,-;

х,-, p,_i) G (х ', Рг- i ; х :, р,) х

XG(x; , р г;

q(' ,

Р ')Ф (|х ; — x i+i\) dxl dx'1d$l. (15.39)

171

Подставляя G\l)в формулу (15.29), легко убедиться, что полу­

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

- T ’ w

i

( [ П « - 1 п ( 1 + П « ) ] ф , (15.31)

К=—со

где

П(*>=2 4я/Л>2 22Я(к>/р2.

1 < а < М

(Р -Р ')Р '

Р2 ехр

2лЫР' у р

2mflp

 

 

Это выражение написано для многокомпонентной плазмы, где М — число компонент; а характеризует сорт частиц.

Вклад диаграмм более высокого порядка в термодинамический потенциал

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

ний (15.18) такой частный ряд, члены которого

содержат кон­

 

 

станту связи е2 в степени на едини­

rv rv

 

цу большую по сравнению с члена­

iru x 1

ми ряда, которые описываются коль­

\svr~

цевыми диаграммами. Очевидно, что

Л/Х?'

'1Л Л '

членам такого ряда должны соот­

 

 

ветствовать диаграммы, содержа­

 

 

щие на одну волнистую линию

 

§

больше, чем кольцевые диаграммы.

 

При этом из групповой суммы для

 

 

Рис. 18. Диаграммы с четырьмя

двух частиц Бц выделяются графи­

линиями взаимодействия, даю­

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

щие ненулевой вклад в груп­

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

повую сумму

S i j K -

 

 

стиц SnK выделяются диаграммы с

 

 

четырьмя линиями

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

и т. д. Ненулевой вклад в SaK с четырьмя волнистыми линиями дают диаграммы, представленные на рис. 18.

График а можно, очевидно, рассматривать как две простей­ шие кольцевые диаграммы, соединенные общей вертикальной линией. В /-частичной групповой сумме следует учитывать диа­

грамму, изображенную на рис. 19,

а, аналогичную диаграмме

а на рис. 18. Здесь имеются два

кольца с одной общей вер­

172

тикальной линией. Правое кольцо содержит i вертикальных ли­ ний, а левое— (/—/) линий, причем г = 2, 3, ..., I—1. Общее число линий взаимодействия равно ( /+ 1).

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

Рис. 19. Два типа диаграмм, описывающих вклад в /-частичную групповую сумму.

всех групповых суммах начиная со второй. Для описания /-ча­ стичной групповой суммы необходимы графики типа представ­ ленных на рис. 19, б. Функция распространения определяется суммированием таких графиков

'

Р;

q'', Р') = 2 (диаграмма на рис. 19, б).

4

' 2<t<l

На этой диаграмме дополнительной линией взаимодействия соединены первая и г-я вертикальные линии. Пользуясь выра­ жениями, соответствующими графикам рис. 19, а и б, можно получить групповые интегралы. Такое вычисление является простым, но достаточно громоздким [4]. В результате вклад графика рис. 19, а в групповую сумму 5; равен

(-!)/+■ (j_i )i

2-2 G}'+V , Р; q'z, Р')>

а соответствующий групповой интеграл есть

Ь, =

-']! 1 G|,+1) (q*, Р; q ' , РW .

(15.32)

Выпишем окончательный результат для вклада системы двойных кольцевых диаграмм, соответствующих рис. 19, а

173

(АИ"), и для модифицированных кольцевых диаграмм, изобра­ женных на рис. 19, б (АО"'):

РД£2"

1

 

СО

1

 

 

 

R*1

Pi

П^*2)

 

 

 

 

 

 

 

 

 

1Pi .Ра

Р2

X

 

4 (2яй)6

 

 

 

мV(l+inf,**))l s P i / V(i+in<K2>)l f e P2 /

 

 

 

 

 

 

 

 

к, ,£к 2~ —оо 0

 

 

 

 

 

 

х dpidpal

 

 

 

 

 

 

(15.33)

 

 

 

1<Sа<М

 

 

 

 

 

 

 

 

 

f c V ’ -

 

 

 

P1P22 2

 

 

 

 

 

,к2) =

pp

 

 

 

pif P;) ФК1 (Рг) фКг (Р„) Ф 2;

 

 

p у J Wa (Pi> p2.

 

 

 

0 0

 

 

 

 

 

 

 

 

 

 

 

 

pAQ"

 

 

CO

1

 

 

 

 

a(Ki,K2)

 

 

 

W-к,,Sк

M

 

 

Лр1 *P2

 

X

 

 

 

 

 

 

 

 

4 (2л/1)

 

 

 

 

 

 

 

 

 

 

 

2= —оо О

 

 

О + ^ Р ? ’ )

( 1+ Е П ^ *))

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X Ф 1Ф 2;

 

 

 

 

 

 

 

л(*1 .К*) _

5

 

< V K!)-

(4лП)геа еГь г аг ь

 

 

Л Р|.Р1

 

 

О ОI

 

 

,2

 

 

Р

 

1<о7ь]<л1

 

 

 

PIP2 IP1 Ра

Г

 

 

 

Р

 

 

 

 

 

 

 

 

 

 

 

. (15.34)

А^а,ь’К‘ = Р ) •

• •

\Wa(PuP2I Pi> Р2) 1^6 (pi*

P2J

Рг> Рз) X

 

о

 

 

 

 

 

 

 

 

 

 

 

 

X 4>Ki (Pi) (p*t 2) ф;, (Pi ) Ф*2(Рз) dp4;

 

 

 

(Фж(Pi) = exp (2jti/cpi/P);

dp2 = dp1dp2;

dp4 =

 

 

 

 

^«(Pl> P2I

Pl> P2)

 

/ ^

V

' - e

x p

l ^

№ -

 

 

 

 

 

 

V2r.ptd

>

 

 

Ptna

 

 

 

 

P‘

Pi (P

Pi)

 

7777”

Рз (P

P2)I*

 

 

2Pma

 

 

 

 

 

2Pma

 

 

 

 

 

 

Эти выражения

записаны

 

для

 

многокомпонентной

плазмы.

Исследуем

вклады

от

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

диаграмм,

кольце­

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

В формулах (15.33) аналитически не удается выполнить суммирование по к и интегрирование по р. Однако если ввести величину v= (рр2/2т)ч=, то с помощью непосредственного

174

вычисления можно убедиться в том,

что для к = О

хорошо

аппроксимируется функцией

Р

/

та \ 3/*

 

т. <°)

 

Ад

1 +

(v2/it) \ 2ярй» )

 

 

 

Тогда из выражений

(15.33)

следует

 

 

Г, , т

й* Ж Ч

4ne»Zap

 

(15.35)

n P(0)=

^ > i

 

 

1+ фрУ2пта)

 

1<а<М

 

 

 

Для простоты будем считать, что массы всех ионов одина­ ковы и равны 1П{. Это эквивалентно рассмотрению плазмы в случае, когда все ионы образованы из атомов одного сорта. Тогда

ПР(0)

й2

 

■+ ft2

1 + ( ^

(15.36)

 

Р7

1 + (%2iPVh)

Р2

р 2/й)

где введены обозначения

 

 

 

 

 

и? =

2

4яе2Z ^ /Хз;

х2 =

4ле22,Р/ХЗ ;

1<а<М

 

 

 

 

 

 

 

К{,е =

 

 

 

,

 

a М — число ионов с различными

зарядами.

Полагая теперь

/с = 0, из выражений

(15.33) получим

вклад в термодинамиче­

ский потенциал от кольцевых диаграмм:

 

 

РА2

_

1

С

с

П 2 (0) dp

 

V

~

2 (2лft)3

\ №

\ -

 

■ т0(0)

 

 

 

1

о

 

 

 

 

 

 

W

П I ( 0 )

 

 

 

 

1

 

dp.

 

 

Й3(2лУ М

 

+

|П „

 

 

 

(0)

 

Параметры %еуц и йех(>характеризуют отношение дебройлевской длины волны электрона к дебаевской длине, которое в рассматриваемых условиях мало. Тогда в выражении, полу­ ченном выше, можно провести разложение по этим парамет­ рам. Если выполнить это с точностью до членов второго поряд­ ка включительно, то

|3AQ'

х3

 

V

12я ( ‘ - т *<-• х

+

— (М*)# 1 —

(15.37)

 

10

 

 

где х2= х 2 + х 2.

Если

полностью

пренебречь квантовыми эф­

фектами, то получим

известное дебаевское выражение:

 

 

— PAQ'/K =

х3/ 12я.

175

Нетрудно показать, что при к Ф О и xefte< l

с той же точностью

получим

 

 

- р Д &/V= (х3 /32я3) хуле2 к -

иj idi -

(х3/64я3) xe%£(5/2),

к=1

О

(15.38)

 

 

где g(5/2) = у]/с~5/2 — функция Римана.

 

К=1

и (15.38)

дают в пределах

В сумме выражения (15.37)

сделанных допущений вклад от кольцевых диаграмм в термо­ динамический потенциал системы. Нетрудно видеть, что послед­ нее из этих выражений менее существенно, поэтому основной вклад в Д£У дает дебаевский член. Оценим порядок членов, стоящих в фигурных скобках формулы (15.37). Если для двух­

компонентной плазмы положить zej t e~ n e= 1020

см~3 и Т—

= 104°К, то К= 1Д • Ю~8 см,

1,45-W cm- \

 

2

\ 2 4 х 3 /

^0,1, а (9/10) (хА ) 2[1—х е/ х ) 2] —0,007. Таким образом, в тер­ модинамике плазмы с плотностью вплоть до Ю20 можно пре­ небречь квантовыми поправками, если ограничиться учетом кольцевых диаграмм.

Отметим, что использование с самого начала больцмановской статистики приводит к неправильной поправке обменного характера. Если корректно учесть эту поправку, то

_

_РД£Г_

_

х 3

1

е

е

х ? ■ У-2

(In 2 — 0,5)

+

V

~

12л

64я

 

i

64я

 

 

 

 

 

 

 

 

 

 

 

+

----- «„X2

 

е

------г;- X3 п2.

(15.39)

 

 

 

64я

е

 

 

2

 

 

Рассматривая аналогично вклад ДО", обусловленный двой­

ными кольцевыми диаграммами, получаем

 

 

 

— РДЙ'/Е = (х2/8) 21

е*п£;

na = Z jW .

(15.40)

 

 

 

 

 

1 < а < М

 

 

 

 

Вклад,

соответствующий

 

модифицированным кольцевым

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

 

 

- РАО"'/v =

- я р

у,

 

 

 

1

 

оо

 

(15.41)

е* е\ nanbР2 f

 

f (dr/r) exp (2rx Y l )

 

 

1 < a ,b < M

 

0

 

0

 

 

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

176

простой способ учета парных столкновений состоит в ограниче­ нии полученного интеграла на расстояниях порядка кулонов­ ской амплитуды рассеяния /0ь ~ е аег>|3. Это не внесет большой ошибки, поскольку зависимость от параметра ограничения ло­ гарифмическая, т. е. достаточно слабая. Тогда

— PQ"7V = —(я/3) V е\ е\ РэпапьIn (xfab). 1<а ,Ь<М

Напомним, что когда вклад квантовых эффектов при рассея­ нии частиц является существенным, т. е. когда ^ ^ \ аь, ограни­ чение рассматриваемого интеграла на малых расстояниях определяется дебройлевской длиной волны.

Модифицированное лестничное приближение

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

m r( I )

V

такую, что

у : г пр (к) dp,

Г Г J 1+ £п<“>

АП' = \МУ (l)dl.

о

Введем также модифицированную двухчастичную

(15.42)

диаграмму

 

 

 

(15.43)

где волнистая линия соответствует

кулоновскому потенциалу,

а пунктирная— потенциалу

некоторого

эффективного взаимо­

действия Фа, 6, вид которого пока неизвестен:

фаЛ I q I . I,

Р=)

( I

q I . I , Р).

177

Используя ранее описанные правила подсчета диаграмм, найдем, что диаграмме (15.43) соответствует величина

РАО' (I)

1

(15.44)

2 (2яЬ)3

V

 

где

U (Р) = I Ф( I Ч I . 6. Р) ехР (ipq) dq.

Определим теперь О (р)

из

условия

равенства подынтеграль­

ных выражений

(15.43)

и (15.44), тогда

 

4я^ V

ngw

(15.45)

0 ( р) =

^ H i + sn<,">

 

Р2

 

Таким образом,

приходим к выводу,

что кольцевая диаграмма,

соответствующая выражению (15.29), и эффективная двухча­ стичная диаграмма (15.43) дают эквивалентные вклады в тер­ модинамический потенциал.

Попытаемся теперь заменить диаграмму, изображенную на

рис. 19, б,

эффективной двухчастичной диаграммой с одной

волнистой

линией и двумя

пунктирными.

Введем обозначение

[см. формулы (15.33)]:

 

-5(к, ,кг) . ,

 

 

РАЙ" (1)

 

 

 

 

Р

apidpi

(15.46)

 

4я (2яЙ)е

2 j

I(1+ЕП<,«‘)) (1+Еп£*>)

 

 

 

 

 

 

Аналогично предыдущему имеем

AQ" = f ДО" © dl.

о

Введем следующую эффективную двухчастичную Грина

\ j \ f \ J

\Ь^>Ч'Лр'> ZZZZ

функцию

(15.47)

Этой диаграмме соответствует величина

 

Р*3(Е)

 

1

Г U (Pi)

U(РА p(Kt’K‘) pfpidp^pt.

V

 

4 (2яй,)в

 

 

 

 

к, ,К,

 

 

Из сравнения этого выражения с

формулой (15.46),

однако,

следует, что

 

 

 

(«1,К2)

 

£/(Pi)tf(P*)

4я-4п

 

.к,)

(15.48)

2 2

X

р

 

Р Г Р г

к, ,к ,

 

 

 

к,,к,

 

 

178

т. е. слева мультипликативная величина по переменным Pi и р2,

а справа нет. Поэтому в

общем случае заменить диаграмму

рис. 19,6 диаграммой (15.47) нельзя.

в правой части

ограни­

Однако если

в суммах

по /С| и к2

читься нулевыми

членами

(ki = /c2 = 0),

то возникает

мульти­

пликативное выражение и

f/(Pi)t/(p2) =

откуда

 

PI

£П<°>

Pi 1 + |Пр

 

 

 

 

 

ц (р) =

. ----1----

 

(15.49)

 

 

 

W

р* 1+П<°>

 

 

Заметим,

что при /с = 0 формула

(15.45) также

сводится к

этому

выражению.

Следовательно,

нулевые члены

рядов

(15.43)

и

(15.48)

можно

описать диаграммами

(15.43) и

(15.47),

в которых

потенциал (пунктир) один и тот

же.

Выра­

жение для последнего при больших и малых импульсах можно получить просто. Так, при малых р

^ Алеа2 1а§

 

4я£

(15.50)

Пр(0) — у?Ip1-, у?

О (р ) —

р- + gx2

1<а<М

 

 

 

т. е. эффективный потенциал О (р)

принимает

дебаевский вид.

В координатном представлении

 

 

 

('

exp (— i pq) dp _

2я2 J

p2 + lx2

-exp (— )/ lyq)

q

Наоборот, при больших импульсах частиц (малые расстояния)

 

ПР(0 ) - 4а4

[Ч 4яe2z $ \* /‘_

 

 

а =

 

(15.51)

 

 

Р4

.1</<М

 

 

 

р%

 

 

 

 

V ( Р) -

 

 

 

 

 

Р4 + 4а%

 

 

В координатном представлении

 

 

Y44|;

------ »_ Г

P4+ 4a4£

dp = — exp ( -

<&•'■)cos

(qal'u) ■

2я2 J

q

 

Для любых расстояний эффективный потенциал можно запи­

сать в явном виде, если воспользоваться

 

аппроксимацией

(15.35):

 

 

 

оо

 

 

 

еаеь Г

Р sjn (pq)

d

(15.52)

Фаь (Я) = ^

 

 

 

q j P ^ '+ ^ p0’)

Ранее было показано, что именно «нулевые» члены рас­ смотренных порядков дают основной вклад в термодинамиче-

179

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