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

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

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

независимых атомов конечного размера. Разместим их в гр.уппы по Ni к одинаковых атомов, что можно сделать g способами:

ЛП

(8.6)

П (N1к)!

 

где g — ч ис л о п е р е с т а н о в о к с

п о в т о р е н и я м и .

Пусть Ni *— число атомов, обладающих кинетической энер­ гией Е к и энергией внутреннего движения WK. При этом второй индекс характеризует набор квантовых чисел, относящихся к внутренним степеням свободы атомов. Возьмем ЛГц атомов и вы­ делим из них один атом. Тогда число состояний, которое может

занять этот атом

(относящийся к уровням Е\

и W\) есть

8i =

у p\dpt

= v yi

4яр7dpj

(8.7)

(2nh)3

(2яЛ)з

 

 

 

где V — объем системы. Число состояний, которое может занять

второй атом из совокупности Л'ц, есть (V—Uu)yi,

где «и — так

называемый о б ъ е м з а щ и т н о й сферы,

по

порядку ве­

личины равный

«объему

атома».

 

 

 

Таким образом, число состояний, доступных второму атому,

равно

 

 

 

 

 

 

[У -(Л 7ц-1)«ц]Т1,

 

 

(8-8)

поэтому количество возможных состоянии, которое может за­ нять группа атомов Nи, равно

Pn = V ( V - Ull) ( V - 2 u n)

. [ V - (yVn - l ) « u ] y f . =

 

 

/ F \

 

— (unY i)

Mii

(8.9)

 

При этом неявно предполагается, что V/un — целое число, что достаточно точно, поскольку V/ttn^l. Если это тем не менее не нравится читателю, то всегда можно записать факториал в виде Г-функции:

 

X! = Т ( Х + 1).

Перейдем теперь к задаче размещения N2\ атомов с внутрен­

ней энергией

Wi и ЕШШ= Е2, т. е. группы атомов, относящихся

к состоянию

(W\\ Е2), Для первого из этих атомов остался сво­

бодным объем V—Л'пИц, а число ячеек, доступных ему в фазо­ вом пространстве, равно ( V■Nnun)y2 - Для второго атома из

группы

N2\ имеем [V—(^u + 1) «и]у2 и т. д. Поэтому число со­

стояний

для группы атомов N21 можно записать в виде:

80

Pn = ( V - /Vn Ull) IV -

(Nu + 1) uu] . . .IV -

(Nn + Nn -

1)X

X «11] v^2‘

KV/«11)-Af11ll

(НцТаЛ ■

(8. 10)

[(V/ku)-(/V u + W21)]!

Далее аналогично рассмотрим группы атомов, относящихся к уровням Е3, Д4, ... Число возможных размещений N\ атомов (т. е. атомов, отвечающих уровню W\) вследствие независимо­ сти атомов в рассматриваемой системе можно записать в виде произведения:

Pi=PuPnPa ■ ■ -Рт -

(8-П)

Тогда

 

)' ■ ■ • (un Y i)

” (М цТг)^*1 •

(8. 12)

Теперь перейдем к рассмотрению группы атомов N2, относя­ щихся к уровню W2. Для первого из них соответствующий объем равен (VN\U\2), где «12 — объем защитной сферы при столк­ новении атомов группы N\ с атомами группы Л/2. Соответствую­ щий объем для второй частицы, очевидно, равен (VN\U\ 2 «22)-

Можно теперь повторить все рассуждения, относящиеся к выводу выражения для р\. Тогда

 

— Мц1г\ иN2

 

 

 

Рг =

м '22

___ J

22

П

(8.13)

 

 

■ЛЦ;12

 

 

 

 

\

u22

 

 

 

 

 

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

Р = РхР2 • • -PN ,

(8.14)

которое следует еще домножить на величину N\/ll(NiK)\ и раз-

делить на N\, поскольку атомы тождественны, а два атома не могут находиться в одной ячейке. Тогда

дь

N1

Р/П ( N u )\

(8.15)

р И (МД!

 

N\

 

 

 

 

81

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

*-П

/ V NxiiiK N2U2K . . . NK—t Ц/с—l. к

 

____________________________________________________U K K _______________________________

X

/ V NiUix—• • - — NK_i ик—t,K

 

\kk

Х

N

 

 

П Vi

Ik

(8.16)

i

Таким образом получено полное число возможностей, соответ­ ствующих распределению заданного числа частиц N по груп­

пам NiK.

Нетрудно видеть, что если положить uiK = 0, то

5^3*0=ПК*П(NtKV (8.17)

что соответствует записи ZN в виде (8.2).

Учет конечного радиуса взаимодействия при выводе выра­

жения (8.16) приводит к изменению

термодинамических функ­

ций, описывающих систему. В этом

 

легко

убедиться,

рассмот­

рев энтропию системы, поскольку

 

мерой

энтропии

является

In 53. Пусть

 

 

 

 

Ук — V ^1и1к— ^ 2М2к — •

--М с -1

(8.18)

Вычислим теперь А = ]п$°. Тогда для Ак имеем:

 

 

 

 

N.

 

 

Ак = \п9ь = \п

икк <Ук!“кк)\

 

 

. (VkIUKK У|с)!

+

 

In

Ук 1 -

■N

Икк

(\

икк

к)

 

J .

=NKIn икк +

\икк

=(In 1/ , — l) + ( j V , - ^ - ) l „ ( l - - ^ 5 - ) .

Поскольку

 

 

 

 

i^j / j

Nj <Mkk_ \ _

N kUk k

___J_ / NKUKK\ a

'

V

Ук )

Ук

2 [ VK J

разлагая в ряд In VK, получим приближенное выражение l

Тогда

л = 2 /,* - л , , п У — ^

- 2

2 “ “ ' ^ * " " -

<а20)

к

к

К*

 

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

 

 

1

2

 

 

2V

 

Кроме того,

 

 

. '■

к

 

 

1*- 1

In NK

 

 

 

 

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

 

 

 

 

N„ ~ ехр |" — р г ' _

т г 2 л'" '““ -

Если записать приближенно

 

 

 

=

Уа К

+ Пк2.)3,

где по порядку величины

 

 

 

V„ ^

3

па3 =^5-10-24 см3,

а

 

о

 

(8.21)

(8.22)

(8.23)

(8.24)

то вместо формулы

(8.2) следует написать выражение

 

ZN =

2 2 n 2exp (/р/n2) ехр (— ап6).

(8.25)

 

П

 

Эта запись требует некоторого пояснения. Если N)^>NK («одноквантовость состояния» системы), то дополнительный к обыч­ ному распределению Гиббса экспоненциальный множитель в выражении (8.23) приближенно можно записать в виде

G = ехр [— (VNX/V) (1 + n2)3] — exp (— ап6).

(8.26)

Отметим, что статистическая сумма (8.25) представляет со­ бой сходящееся конечное выражение. Так, если суммировать начиная с некоторого достаточно большого п= Пь то

оо

оо

У 1 п2 exp (— an2) - j1п2 ехр (— ап2) dn =• - у Г (1/2) (1/а)’/г.

Пj

nt

В табл. 1 видно, насколько существенным является ограничение статистической суммы, обусловленное конечным радиусом взаи­ модействия.

8.3

Т а б л и ц а 1

Зависимость фактора G от п

п

1

2

3

4

G

jq—0,00048

jq—0,0075

Ю-0,06

ю - ° .3

5

о1^

о ел

8

Ю-18

Конечно, приведенное рассмотрение не может претендовать на точное вычисление статистической суммы, поскольку из фи­ зических соображений ясно, что ограничение статистической сум­ мы должно быть связано не только с геометрией, но и, напри­ мер, с температурой системы. Однако такое рассмотрение четко показывает, что даже геометрический фактор (имеющий тот же смысл, что и в уравнении Ван-дер-Ваальса) говорит о несовер­ шенстве записи статистической суммы в виде (8.2) или (8.5). Как это видно из табл. 1, основное состояние атома практиче­ ски не меняется с введением конечного радиуса взаимодейст­ вия. Мало сказывается это и на нескольких первых членах ста­ тистической суммы. Однако для достаточно высоких возбужден­ ных состояний действие фактора G настолько сильное, что вклад этих состояний практически равен нулю. Возможность выделения статистической суммы ZN по внутренним степеням свободы (уровням атома) в идеальной системе связана с тем, что по поступательным степеням свободы можно выполнить простое интегрирование. Статистическая сумма системы, состоя­ щей из N атомов,

ZN =

Sp [exp (— Яр)] = X exp (— £„Р),

(8.27)

 

n

 

где суммирование

проводится по всем энергетическим

состоя­

ниям системы.

 

 

Если система идеальна, то удобно в качестве отсчета энергии

принять о б щи й нуль.

Тогда, как

это хорошо известно из

любого курса статистической физики

(см., например, [4]):

и

(2nh)3 exp (— 8nP) J dV j* dp exp (— p2/2m) =

=

V

(8.28)

 

где сумму по внутренним

степеням

свободы можно записать

в виде

 

 

 

Z/v = 2 £пехр(— 8пр).

П

Вычисление статистических весов gn можно упростить, если пользоваться схемой Рассела — Саундерса, согласно которой уровни, отличающиеся только ориентацией спина и орбиталь-

84

ного момента, объединяются в один мультиплет. Тогда под еп следует понимать среднюю энергию мультиплета, а статистиче­ ский вес

ёп = (2L + 1) (25 + 1),

(8.29)

где 5 и L — полный спин и орбитальный

момент электронных

состояний атома. Спин-орбитальное взаимодействие описывается схемой Рассела — Саундерса с хорошей точностью лишь для легких атомов, когда тонкая структура несущественна (расщеп­ ление мало по сравнению с расстояниями между уровнями). В тяжелых атомах, как это хорошо известно из квантовой меха­ ники [5], необходимо учитывать j—/-связь. Однако, если не до­ биваться большой точности в численных расчетах, можно прак­

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

(8.29). Значения L и S

для основного состояния атома

можно

рассчитать

с помощью

известного

правила Гунда,

если

известна

электронная конфи­

гурация:

наименьшей энергией

обладает

терм с наибольшим

возможным при данной электронной конфигурации значением S

и наибольшим (возможным

при этом 5)

значением

L. Напом­

ним, что полный момент j

выражается через спин

и орбиталь­

ный момент следующим образом:

 

 

 

 

L—S,

если электронная оболочка заполнена меньше чем

I

наполовину,

 

 

 

 

 

L-\-S,

если заполнение оболочки превышает половину.

В качестве примера подсчитаем gn для основного состояния атома азота 7N, обладающего электронной конфигурацией ls22s22p3. Исследовать нужно лишь незаполненную оболочку атома, содержащую три р-электрона. Поскольку в данном слу­ чае 5 макс = 3/2, а возможные проекции орбитальных моментов р-электронов равны 0, —1, 1, то £ макс= 0, и основное состояние атома можно записать в виде 453/2. Тогда

ft = (2S + 1) (2L -]- 1) = 4.

Статистические веса возбужденных состояний атомов невозмож­ но подсчитать столь просто. Значения gn для п> 1 приведены в соответствующих справочниках.

Расходимость статистической суммы по связанным состоя­ ниям есть проявление кулоновского взаимодействия ядро — электрон в атоме. Может быть, можно представить себе следую­ щую картину, если рассматривать частично ионизованную де­ баевскую плазму: электрон, находящийся на достаточно далекой орбите, испытывает экранирующее действие со стороны заря­ женной компоненты плазмы. Не означает ли это, что на него действует дебаевский потенциал ф? Если это так, то поскольку дебаевский потенциал описывает короткодействующие силы, число уровней в статистической сумме ZN должно быть конеч­

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

85

дебаевском потенциале нетрудно оценить. Действительно, если рассматривать достаточно большие квантовые числа п, то дви­ жение электрона на далеких орбитах можно считать квазиклассическим. Пусть электрон движется в центральносимметричном поле и (г). Тогда число состояний электрона, приходящееся на объем фазового пространства [г, r + dr, 0<р(г) ^ р 0(г)], где Ро(г)— импульс в точке поворота:

dN —

р3 .4пгЧг/(2лК)3.

(8.30)

При фиксированном г частица, движущаяся квазиклассиче­

ски, может обладать импульсами, удовлетворяющими

условию

 

 

 

£

-

{

1

г

-

|

( г ,

| )

< 0

что определяет

финитное

движение

частицы

с

отрицательной

полной

энергией

Е.

Подставляя

 

граничное

значение

р0=

= I' 2т\и(г) \ в выражение (8.30)

и интегрируя по всем г,

по­

лучим полное число состояний

дискретного спектра:

 

 

 

 

 

 

 

 

_

 

 

ОО

 

 

 

 

=

16л»

1

з (r) r 4 r

= 2 V 2 _

 

т_^_ Г 2 |

и

р/, dr

(8 _3 1 )

 

3 (2яй)з

J

 

 

 

Зя

 

II»

Jо

 

W l

v

 

Если и (г) — дебаевский

потенциал:

 

 

 

 

 

 

 

 

 

и (г) =

— е2 ехр (— иг)/г,

 

 

(8.32)

то из выражений (8.32) и (8.31) следует:

 

 

 

 

 

 

N =

(32/9 У л ) (а0к)-3/г

2 (1о/а0)3/г .

(8.33)

Пусть %7’ = р -‘^ 1

эв,

я —■1018 см-3,

тогда /в~5-10~7 см,

а число

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

Наиболее простая мысль, которая при этом возникает, за­ ключается в том, что атом не может занимать большего объема, чем это ему позволено межчастичным расстоянием. Если рас­ смотреть ту же дебаевскую плазму, что и в предыдущем при­ мере, т. е. плотность заряженных частиц в плазме порядка 1018 см-3, то легко видеть, что максимальное квантовое число, на котором следует ограничить суммирование,

Пмакс* {rja0)4‘ =• 14.

(8.34)

Это уже разумное значение, не противоречащее эксперименту.

86

Отметим, кстати, что такой простой подход, который, конеч­ но, не дает точного значения пма1;с, оправдан в связи с тем, что статистическая сумма слабо чувствительна к точному значению Пмакс, так как каждый член экспоненциально мал в случае 1/р<С/ и взять 12 или 14 членов не столь уже важно. Важно лишь, чтобы это было разумное число, а число членов суммы было бы не слишком большим. В разное время и разными ав­ торами предлагались различные модификации описанного выше способа ограничения статистической суммы. Все они по сути дела похожи на результат (8.34), и авторы спорили лишь о чис­ ленном коэффициенте. Например, Унзольд предложил выра­ жение

2

=

Z ' U

i i

Пмакс ~~

6

а0 '

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

Разумный подход к ограничению статистической суммы ато­ мов был предложен Инглисом и Теллером в 1939 г. [2]. Если рассматривать водородоподобный атом в поле иона, то вследст­ вие взаимодействия собственного дипольного момента атома с квазистатическим зарядом происходит штарковское расщепле ние уровней, которое по порядку величины равно

ДЕ (dE) ~

еа0пге/г20 — Ае2а0п2пг/*,

(8.35)

где п-—главное квантовое

число; п — плотность ионов.

Ввиду

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

ным дипольным моментом.

 

между уровнями

 

Если ДЕ сравнимо с расстоянием

 

е*_I J ___1

 

^ _2_ ег

(8.36)

во I п2

(п + 1)2

пз а0

 

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

87

линейчатый спектр исчезает и переходит в континуум. Прирав­ нивая выражения (8.35) и (8.36), получаем

П _ п 2^1Ь

‘‘макс '»

Более аккуратный расчет, выполненный Инглисом и Теллером, позволяет написать формулу [4]:

 

Пмакс= 1,04 - 103-тг-а/*5.

(8.37)

Экспериментаторы до сих пор часто пользуются этой

формулой

для оценки

плотности зарядов п в плазме. Однако утверждение,

что именно

это пма1;с определяет наивысшее возбужденное со­

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

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

Е = У F; (г,) =

е V г;гг3 (1 -г кг.) exp (— nrt),

(8.38)

£

i

 

где х= (г®/Зр)-1/2— обратная дебаевская длина. При этом счи­

таем, что возмущающее поле создается системой заэкранирован­ ных по Дебаю зарядов (ионов). Для поля Е можно написать функцию распределения, определяющую вероятность воздейст­ вия на атом данной величины < Е > в зависимости от плотности заряженной компоненты и температуры плазмы. Вопрос о рас­ пределении микрополей в плазме кратко обсудим в следующей главе. Пока примем на веру, что в дебаевской плазме не слиш­ ком низкой плотности п вероятность найти возмущающую ча­ стицу (ион) на расстоянии г с хорошей точностью описывается выражением [3]

dP (г) = 4лп ехр | — и (г) р — 4лп ( (г')гехр [— и (г') Р] dr'j r2dr,

(8.39)

где и(г) — потенциальная энергия взаимодействия иона

с ато­

мом.

При расчете эффекта Штарка в сильном электрическом поле определяется напряженность критического поля Еа, при которой наступает ионизация атома, находящегося в состоянии п. В ре­ зультате расчета порога ионизации получается следующая фор­ мула [2]:

Е0 Ze,/(16n4aj;) .

88

Сопоставление с квантовомеханическим расчетом и с экспери­ ментом показывает, что с достаточной степенью точности можно написать значение Е0, вдвое отличающееся от классического:

E0^Zel(8n‘ a*).

(8.40)

С учетом экранирования возмущающего иона поле, действую­ щее на атом *:

Е (г) — er~2 (1 -fxr)exp(— хг),

(8.41)

где г — расстояние между ядром атома и возмущающим ионом. Приравнивая, выражения (8.40) и (8.41), получаем квантовое число для ограничения статистической суммы:

пмакс = 2~3/‘ (r/a0)'u (1 — хгГ 7*ехр (хг/4).

(8.42)

Для оценки пма,;с можно подставить вместо г в формуле (8.42) среднее расстояние между ионами г0.

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

 

Z(r) = 2

V е gn exp [— en (r) PJ,

(8.43)

 

 

n = I

 

 

где en — энергия

атомных

уровней во

внешнем

электрическом

поле.

параграфа приведем

простой

способ оценки

В заключение

статистической суммы атома ZN. Выражение (8.5) удобно раз­ бить на несколько частей. Если 1/(5<С/ (плазма является низко­ температурной), то достаточно учесть несколько первых членов суммы, например до «точки перевала» [см. рассуждения после

формулы (8.5)],

а остальную часть суммы

учесть квазикласси­

чески, заменив

сумму интегралом. Когда

1/р —/, т. е. возбуж­

денные состояния заметно заселены, большой вклад в статисти­ ческую сумму дают состояния Еп ^1. Тогда

Zn 2 + exp (— р/) V gn ~ 2 + (пмакс— 1) ехр (— р/), (8.44)

П

Если нужна несколько более аккуратная оценка, то можно написать следующее приближенное равенство:

Z h ~ 2 + 2 £ eexp(-eeP) + (tfa-

£ *Л ехр (-Р /),

а=2

\

а=2 /

где значение N2 очевидно.

* Рассмотрим для определенности водород (Z—1).

89

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