Файл: Кудрин, Л. П. Статистическая физика плазмы.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 15.10.2024

Просмотров: 196

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Сравним

L\ с членами, квадратичными по плотности,

L2'~ (e2p )|2.

Отношение L2/L i~Y a. где а = е2хср<^1, у =

= р£о^>1. Первое из этих неравенств более сильное, поэтому

существуют условия, когда отношение L2/Lj мало.

Так, если в

плазме «атомами»

с зарядом

(Z—1)

являются

однократно

ионизованные атомы гелия (Нец), то при Р=^0,06

атм nz~\—-

(l/2)ne^ 9 - 1014 см~3 L2/Li^ 10-2, а

отношение

логариф­

мического члена,

входящего в

выражение (14.10),

к

Lj при­

мерно равно 10-1. Таким образом, когда число одноэлектрон­ ных ионов сравнимо с числом электронов, член в AQ, обуслов­ ленный сдвигом уровня, более существен, чем остальные по­ правки к дебаевскому члену. Поэтому вместо выражения (14.10) можно написать:

= -

S

5' " 4 -

+ f

W)*2ju ' (z‘z,>’(ln5isri^ -

 

/

 

 

ч

 

 

- 2С + - £ - ) - у

(Ре2)3 (2 Zt е%) (2 Z? h) -

( ^ р )

X

X Ъе^

5/ 2

[6ХРфЕт) ~

1 ” P£ml ~ Т (

~

i

ш

 

 

 

 

(14.12)

Это выражение следует дополнить условием химического рав­ новесия

 

 

 

Мд “t” М'е

Мд_1

 

 

(14.13)

и уравнениями

 

 

ай

 

ай

 

 

 

ай

=

пZ — 1

 

= — п.

(14.14)

 

1

 

 

apz

d\ie

 

 

где

tij — плотности соответствующих компонент

плазмы. Выра­

зив

химические

потенциалы

щ

через

реальную плотность ча­

стиц, согласно уравнениям (14.14),

получим

из уравнения

(14.13)

 

 

 

 

 

 

 

 

^ г ( ^

)

" еХр(_Р' ) [ 1

+ , <"г— " г. » .)1 = 0 .

04.15)

причем / < 1; / — потенциал ионизации. Функция f достаточно громоздка и не будем выписывать ее здесь [11]. В случае водородной плазмы (Z=l)

/ ~ ^ H D+ y ( e 2PxD)2,

где

х2, = 8ле2Рде.

159



§ 15. МЕТОД КВАНТОВЫХ ГРУППОВЫХ ИНТЕГРАЛОВ МОНТРОЛЛА И УОРДА И УРАВНЕНИЕ СОСТОЯНИЯ ПЛАЗМЫ

Интересный, хотя и несколько сложный, подход в термоди­ намике системы кулоновских частиц был предложен Монтроллом и Уордом, которые обобщили метод групповых интегра­ лов Майера на квантовый случай [12]. Метод не позволяет по­ лучить термодинамические функции сильно неидеальной систе­ мы, поскольку он базируется на теории возмущений. Однако с его помощью можно, вообще говоря, продвинуться при вы­ числении термодинамического потенциала в область больших плотностей плазмы. Квантовомеханическое рассмотрение позволяет преодолеть трудности классической теории, связан­ ные с малыми расстояниями между частицами. Кроме того, ес­ ли плотность частиц не слишком мала, имеет смысл вводить в

термодинамические функции поправки, связанные с вырожде­ нием системы, т. е. с отклонением от статистики Больцмана.

До'

плотностей п ~ 1020 смгг

и температур р- 1~

1 эв и

выше

дебройлевская длина волны

электронов Ае остается еще мень­

шей

среднего расстояния между частицами г0.

Поэтому,

если

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

большим поправкам.

Рассмотрим систему кулоновских частиц в объеме V при температуре (Н в состоянии термодинамического равновесия. Термодинамические свойства системы N частиц можно вычис­ лить, зная ее статистическую сумму (для канонического ан­ самбля) :

Z„ = 2 f^ (q * O e x p {- p tf}¥ J q ^ q " = Spexp(-p//) (15.1)

{"} V

или статистическую сумму системы с переменным числом частиц, соответствующую большому каноническому ансамблю:

Arf

00

z" Ц J

(-/0 exp( -

хч

 

2 =

2

т Vn(q*) dq» ^

 

 

if=o

|n j

 

 

 

 

 

- S

p £ zN exp (-

№■

(15.2)

 

 

 

P i

 

 

Здесь dqN = dqidq2 ■■■ dqN\ dqi^dXidyidZi. Суммирование про­ водится по всем наборам квантовых чисел {n} = ni, п2, ...

Функции 'Кп(ч^) образуют полную ортонормированную систе­ му, симметризованную соответствующим образом. Гамиль­ тониан

й = - V ^ r - A , ( + tf(qw).

(15.3)

лшЛ

KKN

160


a z = ePu,

где и- — химический

потенциал. Если

 

— соб­

ственные

функции оператора

(15.3),

то

(15.1)

и (15.2)

можно

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

 

~

 

со

 

 

 

 

 

 

Z.v = 2 ехр{— P£n)>

 

2

/

ехр ( - Р а д ,

(15.4)

Za, =

V

 

{n}

 

 

 

JV=0

{п}

 

 

 

 

 

где Еп — собственные значения гамильтониана

(15.3).

 

 

Функцию Грина N частиц удобно записать в виде, завися­

щем от двух систем координат

qv

и q "'

и

двух

обратных

температур р и р':

 

 

 

 

 

 

 

 

 

 

Слг(Ч" ,Р ; qA,',P') = S

< ( ч П ы р { - Ф - Р ) Щ У аЮ -

 

(15.5)

 

{"}

 

 

 

 

 

 

 

 

 

 

Предположим, что GN определено так, что

G ^q^, Р;

qN', р') =

= 0 при

р < Р '. Если

Р = р/, то

из

условия

ортогональности

функций Ч'п следует

 

 

 

 

 

 

 

 

 

 

 

Gtffo",

Р; qJV',P') = n6 (Qi — q') .

 

(15.6)

 

 

 

 

 

1<»<JV

 

 

 

 

 

Из сравнения выражений

(15.1)

и (15.5)

получаем

 

 

 

Z jv =

f Gn(q w ,

P; q^ , 01 dqN,

 

 

(15.7)

 

 

v

 

 

 

 

 

 

 

 

 

t . e. знание функции Грина определяет все термодинамиче­ ские свойства системы.

Величина GN(qNP; qN, 0) представляет собой амплитуду вероятности для системы частиц, находившихся в данных точ­ ках пространства qN при бесконечной температуре, возвра­ титься в те же точки пространства после медленного охлажде­ ния системы до температуры р-1. Таким образом, статистиче­

ская сумма есть среднее от амплитуды

вероятности такого

процесса по всем возможным начальным точкам.

 

 

 

Известно,

что функция

GN является решением в виде функ­

ции источника уравнения Блоха:

 

 

 

 

( 4 г +

**)

(4W- Р;

Ч"\ Р') =6 (Р -

Р') п 6 (Qi. _

q;),

 

 

\ op

J

 

 

1 < 1< л Л

'

 

где оператор в скобках действует только

на переменные

qA' и

р. Это уравнение в интегральной форме имеет вид:

 

 

 

G*(qw, Р; q " \ P0=G°w (qv . Р; q"', Р ') -

 

 

_

f n (qw , Р;

х", n V { x " ) G N (х" , Р";

q"', p')dx" dp",

(15.8)

 

Ь'

 

 

 

 

 

 

где

G°w — функция Грина,

описывающая систему

невзаимо­

действующих частиц.

 

 

 

 

6 З ак . 635


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

Точно уравнение (15.8) решается лишь для идеальной си­ стемы. Предположим для простоты, что система невзаимодей­ ствующих частиц подчинена статистике Больцмана и удовлет­ воряет граничному условию для GK

Gn (qw, Р; q^', Р')

0

при qw-> оо.

 

 

Тогда решение уравнения (15.8) с

учетом равенства (15.6)

имеет вид

п

 

 

 

G°(q", Р; q " \ Р') = 1/!

2п (Р— |У) h2

■ г

х

 

1<;<лг

 

 

mt (qt -

q t' ) 2 \

 

(15.9)

 

 

 

 

2 А * ( Р - р ') 1

где фактор ЛП учитывает число возможных перестановок из N одинаковых частиц. Полагая qAr = qftr' и Р' = 0, легко получить из формул (15.7) и (15.9) статистическую сумму

VN

Г \ (

пЧ YG

(15.10)

2.N ЛП

* Ц

2яйД2 ) '

\ < i < N

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

Спиновые состояния системы невзаимодействующих элек­ тронов можно учесть, записав функцию Грина, антисиммет­ ричную по координатам qiV' и qX т. е.

G°(q\ 5; q!'P')

1 (q2, f5; qr P )

N1

G° (q^p; q>'p')

=

n (qN . Р; q*', Р') =

(q1, P; q2', P') .

(q1, P; q^'.P')

(q2, P; q2\ P') .

.

G° (q2, P; q^,P'}

G°(qXP;q2',P') .

.

G° (qM, P; q"',p')

del (qw , р; q*', р').

Это выражение эквивалентно детерминанту Слэтера, по­ строенному из одночастичных волновых функций системы.

Рассматриваемый формализм по существу ничем не отли­ чается от формализма функций Грина, который обсуждался в шестой главе. Однако математический язык изложения яв­

162