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

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

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

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

Добавлен: 15.10.2024

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

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

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

термодинамически неравновесные условия неизотермической плазмы, когда температура ионов превышает температуру элект­ ронов 7’»>Ге. В этом случае, как показал В. И. Коган [4], влияние теплового движения ионов при вычислении W(E) ока­ зывается более существенным, чем влияние корреляции. Интерес представляют, конечно, не только предельные в этом смысле условия импульсного разряда, но и случай термодинамически равновесной достаточно горячей плазмы.

Пусть Е (/)— векторная сумма электрических полей, созда­ ваемых совокупностью электронов и ионов плазмы в фиксиро­ ванный момент времени t в точке наблюдения, в которую поме­ щен атом. Пусть атом находится в начале координат. Тогда Е(^) можно представить в виде:

 

N

Гр/ +

Vjt

 

Е (t) = e

гок ~Ь VKt

(10.30)

к=1 I r0K + vKt I 3

I г0у +

v,t\3

 

 

где индексами / и к нумеруются соответственно ионы и элект­ роны; N — полное число заряженных частиц того и другого знака в объеме V.

Разбивая отрезок времени [0, т] на М одинаковых достаточ­ но близких участков, можно получить нормированное на еди­ ницу статистическое распределение (Ei, Е2...... Еп) значений напряженности поля в последовательные М моментов времени

ta^ax/M, т. е. плотность вероятности

того, что

в момент

ta

поле в месте нахождения атома есть

Еа(п=1,

2,..., М).

Это

распределение нетрудно найти методом А. А. Маркова, обобщая его с трехмерного на ЗМ-мерный случай. Считая, как и при выводе формулы Хольтсмарка, что распределение электронов и ионов является случайным при заданной плотности частиц п (т. е. распределение есть распределение Пуассона), и переходя

к асимптотическому пределу N-+-oo,

V-voo,

п = const, получаем

 

 

, Е М) = —

. .

м

 

 

 

W(Ej,

Е„

| П

Mpaexp(ipaE j] X

 

 

' '

 

 

а=1

 

 

 

 

 

X ехр [— NC (рх,

р2, .

.

рм)] ,

 

(10.31)

где

 

 

 

 

м

 

 

 

 

 

 

 

 

Гр +

vta

X

С (Pi>

Рг>

 

 

 

 

 

 

 

 

I Гр +

via

 

 

 

 

| 3

 

 

X f/г (V) +

fe(v)] dvdr0,

 

 

(10.32)

/i(v). M v) — нормированные на единицу максвелловские рас­ пределения ионов и электронов.

108


Вычисление функции (10.31) в общем виде не представляется возможным. Однако при наличии малого параметра можно по­ лучить не очень сложные аналитические выражения для W(a). Так, если рассматривать адиабатическое приближение при вы­ числении действия микрополя на атом водорода или водородо­ подобный атом, то малым параметром задачи можно считать величину

k~l = n V ( а '/ 4 г 3 , .

(10.33)

где v0i= (2/^пц)112; п,-— плотность ионов; а' — штарковская по­ стоянная:

а' — q (0 — — — ,

выраженная через мгновенный адиабатический сдвиг частоты рассматриваемой штарковской компоненты * q(i).

Тогда обобщенное распределение'Хольтсмарка, учитывающее тепловое движение ионов, можно представить в виде

где

 

 

W (а)

=

W h (ос) +

k ~ ' h S{a),

(10.34)

 

 

 

 

 

 

 

S(; c ) - - U

J L

V'-

L(x) — I

(x)

 

w

V .4

/

dx3

/ x

 

 

 

 

 

 

 

 

(10.35)

 

 

 

 

 

 

 

L(x) и T ( x ) — функции, введенные Чандрасекаром. При хЗ>1:

L ( * ) « ] / 2/п;

7(x)^(2/3)Y2jn .

Функция S(x) табулирована в работе В. И. Когана [5] и при­ ведена на рис. 8. Отметим, что по самому смыслу вывода форму­ лы (10.34) второй член в ней является малой поправкой к ос­ новному (хольтсмарковскому) члену. Легко видеть, что эта по­ правка на тепловое движение возмущающих ионов пропорцио­

нальна температуре.

Асимптотика распределения Хольтсмарка, соответствующая случаю распределения ближайшего соседа, есть

WH (*) ~ (15/4 |/2 я ) х~‘/2\

(10.36)

Как видно из рис. 8 и формулы (10.36), поправка, вносимая в распределение ближайшего соседа при x ~ 2-f-3, составляет ощутимую величину. Отметим, что рассмотренное выше прибли­ жение Гофмана и Теймера тоже по сути дела построено на рас­ пределении А. А. Маркова.

* Напомним, что в электрическом поле термы атома водорода расщеп­ ляются на штарковские компоненты, причем величина расщепления пропор­ циональна напряженности внешнего поля Е (линейный штарк-эффект).

109



Наиболее последовательно учет корреляции заряженных ча­ стиц в дебаевской плазме при вычислении распределения элект­ рического микрополя проведен, по-видимому, в работе [1], в ко­ торой показано, что функция W(a) выражается через парные функции корреляции кулоновских частиц и диполя, помещенного в плазму. Таким диполем может быть атом водорода или водо-

Рис. 8. Функция S(x), учитывающая влияние теплового движения ионов на распределение микрополя в плазме.

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

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

W{ Е):

A (iq) == ] exp (iqE) W (E) dE —

 

= Zn Xf exp { - p [H (X) + ip -‘qE (X)]\ dX,

(10.37)

 

v

 

 

где V — объем, занимаемый плазмой; ZN— конфигурационный

интеграл;

Н(X) — гамильтониан взаимодействия

системы; А'—

совокупность пространственных координат частиц;

N — число

частиц в системе.

по

сравнению

Пусть

— некоторая достаточно малая

со средним межчастичным расстоянием величина. Присвоим вы­

деленной частице в плазме

индекс 1 и введем

в гамильтониан

системы дополнительный член, учитывающий

короткодействие:

2

? ( * - '/ / ) .

'(10.38)

1—2

 

 

ПО


где

О

при

7? — г1г < О;

V (R — ru)

оо

при

R —ги > 0;

 

г1г*= | г1—г*|; гг — координата i'-й частицы.

Нетрудно видеть, что соответствующая новому обобщенному гамильтониану системы функция Лн(т1) определена в области действительных значений |R e q K ° o комплексной переменной г). В силу экспоненциального вида подынтегрального выражения в формуле (10.37) нетрудно убедиться в том, что эта функция

также и аналитична

в указанной

области *. Следовательно,

функция Ля('п) при

вещественных

q однозначно определяет

функцию ЛяОц), которая, в свою

очередь, однозначно опреде­

ляет A(iq) как Пт Ля (iq).

 

 

 

R-+0

 

 

Запишем электрическое поле в виде векторной суммы

 

 

N

 

 

E(X) =

^ Z <V - ^ ,

(Ю.39)

 

 

i=2

 

где Zi — заряд

t-й частицы.

Продифференцировав

1пЛн(т)) с

учетом (10.39)

и затем почленно проинтегрировав

полученный

результат от 0 до q, получим искомую связь фурье-компоненты функции распределения с парной функцией корреляции:

 

 

N

q

 

A

(iq) = lim exp

— i V z A

■'

— Г d r ^ X

 

я-»о

i=2

q J

 

 

0

 

 

X K r

ii rii)q v

 

(10.40)

 

 

p

 

 

где

 

 

 

 

 

I exp (Pfr„t

drsdr4 . . . drN

K r

( « , rlt) =

 

 

(10.41)

{ exp~(— P # ;, R) dhdr2 ■ . -drN

ZtZj

H a , R — о

(Zj +

« V )

X

2

ru

 

 

i-

1./■si t-/= 1

 

 

N

N

 

(10.42)

 

rii):

* ~

"£=2

iqp

1=2

 

 

* Это обстоятельство также полезно проверить читателю.