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

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

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

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

Добавлен: 15.10.2024

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

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

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

Напомним, что неравенство (16.2) позволяет вычислять инте­ ресующие нас величины по теории возмущений. При классиче­ ском описании бинарной функции корреляции разноименных по знаку зарядов возникают трудности, связанные с расходи­ мостью этой функции на малых расстояниях. Эта расходимость обусловлена «термодинамически выгодным» падением разно­ именных зарядов друг на друга. Если не вводить модели с силами отталкивания малого радиуса действия, то положение может быть исправлено лишь с помощью учета квантовых эф­

фектов.

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

К (г) = Ккл(г) + Ккв(г),

(16.3)

где Ккв (г) можно назвать к в а н т о в о й к о р р е л я ц и о н н о й фу нк ц и е й .

Отметим, что квантовая природа частиц проявляется лишь

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

кото­

рая в поставленной задаче много меньше дебаевской

длины

(З^е-С/в). На столь малых расстояниях экранирование несуще­ ственно, и Кип (г) можно сравнительно легко вычислить, ис­ пользуя известное решение чисто кулоновской задачи двух тел.

На больших расстояниях корреляция частиц

описывается в

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

Поэтому можно

написать:

 

 

Ккл (г) — 1 — [ехе2Рexp (— xDr)/r] =*exp [—е ^

exp (—xDr)/r], (16.4)

где е, и е2 — рассматриваемые заряды.

 

препятствуют

Как известно, квантовомеханические запреты

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

Если

рассматривать область

 

то задача. вычисления

К (г)

является двухчастичной и в выражении (16.4) можно по­

ложить хл=0- Тогда в соответствии с

квантовой

статистикой

функция К {г) определяется

диагональным элементом матрицы

плотности

 

 

 

 

К (r12) = const У ехр (— р£п) I

Фп (н, г2) |2,

(16.5)

 

П

 

 

 

где фп(г1, г2) — собственные

волновые

функции

задачи двух

тел с кулоновским взаимодействием. Выделяя движение цент­ ра масс (плоская волна) и интегрируя по суммарному импуль-

* Это обстоятельство является прямым следствием принципа неопреде­ ленностей.

187


су, получаем

К (г) = const 2 exp (— рек) | ip. (г) |2,

(16.6)

к

 

где фк (г) описывает движение частицы с приведенной массой т* те в поле неподвижного кулоновского центра; ек — энер­ гия этого относительного движения. Функции фк (г) для непре­ рывного спектра имеют вид [6]

% (г) = /о '"-Tv

ехР

л/2к) Г 11 + (*^)1 ехР (ikr) х

 

(2я) /!

 

 

 

 

X F[—\/k, 1, i (kr — kr)\,

(16.7)

где F — вырожденная

гипергеометрическая

функция;

Г(х) —

гамма-функция. Нормировка волновых функций обычна:

J ^k (r) 'V

(r) dr = 6 (k — k')-

(16.8)

Если выполнено условие

борцовского

приближения

(16.2),

то связанными состояниями можно пренебречь (для одноимен­ ных зарядов они вообще отсутствуют). Поскольку k велико, выражение (16.7) можно разложить в ряд и ограничиться пер­ вым приближением по к~1, считая kr произвольным. Тогда

(16.9)

Так как с необходимой точностью

где С

— постоянная Эйлера, то для бинарной функции корре­

ляции

получим

(16.10)

188

где

 

 

 

 

а = те е\ е2$/2Л2.

 

 

 

После

простых

преобразований

интеграла

получим

 

 

 

K(r) =

1

+

 

2*1^2

Г ^

ехр (— <а)

(16.11)

 

 

К <и> J

 

 

 

 

 

 

 

 

П

ПК

 

 

 

где < о >

-

 

 

й -

 

 

 

 

 

 

 

me <v>

 

 

 

 

 

 

 

 

 

 

 

Постоянная

в выражении (16.10)

выбрана

так, чтобы при

г~Э> ^ эта формула

переходила

 

в классический

результат (16.4).

Сравнивая выражения (16.3), (16.4) и (16.11), получаем

 

 

 

КаьМ

Тг <v>

~

ехр (

Р) —

 

 

 

 

 

 

 

 

 

 

 

= ^ £ iЬ - п'/’ф ( ±

 

 

 

 

 

(16.12)

 

 

н <v>

 

\ к

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

где Ф(л:) =

—Д= Г ехр (— t2)lt2dt

интеграл ошибок.

 

 

 

V х J

 

 

 

 

 

 

 

 

 

 

о

 

что при г—>-0 квантовая поправка в форму­

Нетрудно видеть,

ле (16.11) компенсирует расходимость, характерную для клас­

сического выражения бинарной функции корреляции.

Так, при

г= 0

 

К(0) = 1 — 2яе1е2/Й < у> .

(16.13)

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

дует учесть

антисимметричный

триплет

и симметричный син-

глет. Следовательно, вместо

|фк(г)|2 в формуле

(16.6) необхо­

димо записать следующее выражение:

 

J1-j- +

3 I Фк (О — Фк( - г) I2 + I Фк (0 + Фк (— г) |2 -

 

+ -

Г

Sin

 

-

кг) X ] -

ехр (2ik-^

X

 

k

J

х

 

 

 

2

 

 

 

о

 

 

 

 

 

 

X

14------- 1----- (

J

— exp(— lkrx) sin (krx)

(16.14)

 

L

k

k

x

 

 

 

 

 

 

о

 

 

 

189



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

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

К (г) =

I —

Г exp ( -

г) f

Ккв (г) + Коби (г), (16.15)

где

 

 

 

 

 

 

 

/Собм (0

 

Д- ^

^ U j* dk exp (—ak2 + 2ikr)

Jl

 

 

2

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

i

exp (— ikrx) sin (krx)

= - ^ e x p ( - r 2A2) -

 

k

J

x

 

 

 

 

 

 

 

ea

f ^ exp (— ta) — exp [— (г/д)2]

 

(16.16)

 

 

h < v >

J

t2-(r/%)2

 

 

 

 

 

r/%

Первый член описывает обменную корреляцию двух одинако­ вых частиц в идеальном газе, а второй — обменные поправки, учитывающие взаимодействие. При r^> &эти члены экспонен­ циально малы, а при г-^О оба они конечны и вносят существен­ ный вклад в функцию корреляции.

Зная корреляционную функцию, можно вычислить любую термодинамическую величину и получить уравнение состояния системы (см. § 4). Представим термодинамический потенциал П слабо неидеальной плазмы в виде

 

 

е*

 

 

 

 

+ -JV- Jге2

ft

 

 

 

de2

(16.17)

 

 

< л вз> ,

где

 

0

 

 

 

3/

 

 

 

 

(16.18)

 

 

= П0 + пл

 

 

Р/2

 

— потенциал

идеального электронно-ионного газа

с учетом по­

правочного к По

члена, учитывающего слабое

вырождение

электронное

[5];

Нвз — гамильтониан

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

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

жении

(16.17), можно легко

выразить

через К (г). Действи­

тельно,

поскольку

 

 

 

 

 

(16.19)

 

i<i'

е<е'

ie

190


то нетрудно видеть, что

у < £ в з > = - f - J у (K,i (г) + к ее (г) - 2к е1(Г)}, (16.20)

причем выражение в фигурных скобках, согласно сказанному выше, равно

Кц +

Кее -

2Kle = -

(Z + ' )2g2p exp ( -

xDr) +

 

 

+

\кее(Г) + /CSfM(Г) - 2К1е(г)].

 

(16.21)

Первый член приводит в термодинамических функциях

плазмы к обычной

дебаевской

поправке

(см.

вторую

главу).

Поэтому вместо выражения (16.17) можно написать

 

 

Q = Q0 + «Р-1

3j

 

 

 

 

п%3-|- AQd +

 

 

 

 

 

21/2

 

 

 

+

J de> f ^ [Кее (г) + К°е!м (г) -

2К* (г)].

(16.22)

 

о

 

 

 

 

 

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

функций

(16.22)

и

(16.16)

и учитывая,

что при электрон-элек-

тронном взаимодействии

(ее) следует брать приведенную мас­

су т* — те12, а при электрон-ионном взаимодействии

(ei) мож­

но считать т*жтс, получаем

 

 

 

 

 

 

оо

оо

dt exp (— t2)

 

 

 

У -7 К„ =

0I

rdr j

Jt3^2

Ze2Ti$

(16.23)

r j % ee

i2

2,/z

m’/2

 

 

 

 

 

 

 

dr „

0„

Ze2

j" rdr

dt

exp (— t2)

n /s

Zemy

----Aei

%Vel

t2

 

21/2

 

Г

 

о

r f t el

 

 

 

 

 

 

 

 

 

Далее

со

exp (— r2)

_|___

| ^ ^ м = _ 4 л t f ^ r d r

2

Tiv„

 

(16.24)

X

X i dtexp (~ /2)~ exp

= — Л

1-

- n u In 2] .

'

t2 — r2

 

m

 

(16.25)

191