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

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

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

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

Добавлен: 15.10.2024

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

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

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

Вычисление / с учетом условия (IV. 13) приводит к выражению

МI,к 2

72

Vk.i = 8(2я р т )'/г

(IV. 14)

ао

Л“ к.1 ‘

Эта формула демонстрирует существенное отличие квазиклассического выра­ жения для ук,г в кулоновском случае от формулы (11.33), полученной для

случая плоских волн.

2. Рассмотрим случай малых передач энергии при взаимодействии с из­ лучающим ионом, т. е. ы<СР-1. Имеет место следующее соотношение:

____________exp (2nZ1/p1a0)____________

[exp (2nZ1jpiaa) — 1] [exp (2nZ1/p2a0) — 1] *°

_1_

2pi

¥ (1 ) — Re <F

я2 { -

Pi --P2

 

гд е 'F — логарифмическая

производная Г-функции,

d

dxо F (**) I2=

iZi

P1^0

a p2= ( p \ w ) ^ 2 . Про­

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

Vk,; = -

8 (2nm$Ne*)l/' |M / k |2 In (ftcok>;p),

 

h*

который совпадает с соответствующим членом для нейтрального атома без

учета цепочки из петель.

■ Для устранения расходимости при o)k,i— >-0, а также для вычисления по­ правок к логарифмическому члену необходимо учесть цепочку из петель, ко­ торая описывается уравнением (IV.9).

 

СПИСОК ЛИТЕРАТУРЫ

 

 

1.

Зоммерфельд А. Строение

атома и спектры. Т. II. М., Гостехиздат,

1956.

2.

Кудрин Л. П., Тарасов Ю.

А. «Оптика и спектроскопия», 1964, т. 17,

с. 489.

ПРИЛОЖЕНИЕ V

ТЕРМОДИНАМИКА РАСТВОРОВ СИЛЬНЫХ ЭЛЕКТРОЛИТОВ УМЕРЕННОЙ КОНЦЕНТРАЦИИ

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

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

■Пкл = еоеф/(г0е),

(V.1)

где Го — среднее расстояние между ионами; е„ и еь— заряды ионов а и Ь. Существенно, что это выражение содержит в знаменателе диэлектрическую про-

477


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

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

ленных систем, на

концентрированные растворы. При этом не учитывается,

что между ними существует определенное качественное различие.

Действительно,

в условиях применимости теории Дебая — Хюккеля

(т]клО ) радиус

корреляции RK между положением отдельных ионов зна­

чительно больше среднего расстояния между ними, а

средняя энергия

элек­

тростатического

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

мала по сравнению

с кинетической

энер­

гией. Амплитуда

рассеяния при

достаточно высоких

температурах, опреде­

ляемых критерием г)Кл<^1, велика по сравнению с размерами самих ионов. Поэтому в сильно разбавленных системах роль короткодействующих сил, т. е. взаимодействия на малых расстояниях, несущественна. С ростом плотности

ионов п объем дебаевской сферы

uD = (4/3)яг^убывает

как

п~ 3/а, а сред­

ний объем, приходящийся на одну

частицу v — как /г-1.

В

результате, на­

пример, в водных

растворах

при п ^ (0,054-0,l)/Ze моль/л (Z — заряд

иона) v становится

сравнимым

с vn■Одновременно система

испытывает не­

которую перестройку, связанную с тем, что корреляция заряженных частиц существенно усиливается при v > v D.

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

ными частицами близко к радиусу

действия короткодействующих сил Ro-

При этом становится существенным

учет сил отталкивания на малых рас­

стояниях в термодинамических функциях, описывающих систему. При пере­ ходе к очень концентрированным растворам или к расплавам ионных солей радиус корреляции RK снова возрастает и в пределе становится равным нескольким Ro. Структура расплава становится похожей на структуру обыч­ ных жидкостей, которая в основном определяется борцовскими силами от­ талкивания.

Ввиду того что для умеренно разбавленных электролитов v ~ v Ilt т. е. параметр взаимодействия т)кл —1, вириальное разложение термодинамиче­ ских функций вряд ли может привести к успеху. Необходимо использование теоретических методов, заимствованных из теории жидкостей. Здесь можно попытаться добиться некоторого успеха аналитически, поскольку малый без­

размерный параметр типа p = 2ntiR jj <S1.

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

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

твердые

сферы радиусом

Ro. Согласно этой модели, потенциальная энергия

парного

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

частиц сорта а и b

 

 

 

 

 

 

uab=us(r) +иэаль(г),

 

 

(V.2)

где

 

 

 

 

 

+ оо

при 0 < г <

R0

 

еа е Ь

0

при R0< Г<

оо; лаЬ(') =

tr

 

При переходе к концентрированным растворам или к расплавам ионных со­ лен оба эти упрощения теряют силу. Поэтому излагаемый ниже метод не охватывает области высоких плотностей. В этой области скорее приме­

478


ним так называемый метод свободного объема, основанный на определенной аналогии между жидкостью и кристаллом.

Рассмотрим раствор электролита, находящийся в объеме V в состоянии термодинамического равновесия при температуре р~’. Поместим начало коор­ динат в центр иона а с зарядом еа■Взаимодействие выделенного иона с за­ ряженными частицами системы приводит к образованию вокруг иона «де­ баевской атмоаферы» с плотностью

 

па ( т ) =

5 ]

W bK abir),

(V .3 )

 

 

1<ь<м

 

 

 

где Каь(г) — бинарная

функция

корреляции,

определяющая

вероятность

найти ион сорта Ь на

расстоянии

г от

иона а.

Условие квазинейтральности

можно представить в следующем виде, воспользовавшись известной форму­

лой электростатики,

описывающей

распределение потенциала <ра(г)

вблизи

заряда а в начале координат*:

 

 

 

 

 

 

 

Фа (г)

а_

| па (г')

dr'

 

 

(V.4)

 

 

8Г

е I г — г' |

 

 

 

 

 

 

 

 

 

 

 

 

 

V

 

 

 

 

Преобразуя

входящий в

это выражение

интеграл

с

помощью формулы

Грина и учитывая, что па(°о)=0, получаем

 

 

 

 

 

 

0 0

СО

 

/

СО

\

 

Фа =

■—

I*

I* «а (0 t * d t + - ^ - \ e a + 4n

I* na (r) r*dr\.

(V.5)

 

8

г

i

 

I

Я»

>

 

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

Игл г3 IKab— U = 0.

г-*-оо

со

Например, флуктуации плотности, пропорциональные J [Каь—l]rfr, при па-

о

рушении этого условия неограниченно велики. Следовательно, па должно стремиться к пулю заведомо быстрее, чем r~s, a (fa — быстрее, чем г-1. Но последнее возможно лишь при условии, что выражение в фигурных скобках формулы (V.5) равно нулю, т. е.

еа + 4п | X) ebnbKab (r) f2dr = 0.

(V.6)

Ло1<Р<М

 

Это условие отражает тот факт, что бинарная функция корреляции должна при г-—>-оо убывать столь быстро, чтобы сумма зарядов центрального иона и окружающего его облака равнялась нулю.

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

Kab (г) = V (г) Gab (') .

причем множитель у(г)=ехр[—ця(г)|}] не зависит от зарядов ионов. Поэтому

связь между К а ь и ф0

задается функцией

Gab. Функция корреляции, по оп­

ределению, должна

быть симметричной функцией еа и еь. С другой стороны,

из выражений (V.4)

и

(V.6)

следует, что ф„ пропорциональна еа.

Тогда

Gab может зависеть

лишь от

произведения

ц>аеа. Других комбинаций

этих

* Предоставляем читателю возможность убедиться в том, что функция (V.4) является решением уравнения Пуассона.

479



величин не существует. Поскольку, далее, Gab безразмерна, то она обязана быть функцией безразмерного потенциала

еа^ а ь = ~ - pY (г).

 

ел0

 

 

Таким образом, если Каь определяется

величиной сра, то эта

зависи­

мость должна иметь ©ид

 

 

 

КаЬ У@аЬ

еаеЬ

¥(г)

(V.7)

 

Е

 

 

Введем характерные безразмерные параметры задачи. Отношение объе­ ма, приходящегося на одну частицу, к дебаевскому объему vD можно запи­ сать в виде

v/aD = 3-/х,

(V.8)

где

 

X2 8яя22е!/?0[}/е = Z2e2fi/sRa\ х 2 = 4-/р;

р = 2ттН^.

Параметр х выражает отношение кулоновского взаимодействия на малых расстояниях к кинетической энергии, параметр р — отношение собственного объема частицы к vD. Таким образом, явно введено короткодействие и пред­ положено оно бесконечно сильным при r<R 0. Считая комбинацию, входящую в аргумент функции Gab, малой Величиной, (разложим неизвестную функцию Gab в ряд и ограничимся первыми двумя членами разложения:

 

 

 

Каь (г) = V (г)

еаеЬ

%&Ч(г) ■

 

 

 

 

 

 

 

Z4*

 

 

 

 

Входящую

сюда

неизвестную

константу

В2,

являющуюся,

вообще

говоря,

функцией

безразмерных параметров

задачи

х и х> Дебай

приравнял

еди­

нице, что справедливо при рассмотрении

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

слабонеидеаль­

ной плазмы.

 

 

 

 

 

 

 

 

лишь

В потенциальную энергию кулоновской системы u,v заряды входят

в виде

парных

произведений.

Поэтому

Gah= Gab (еаеь)-

Представим

в аь

в виде

 

Gab = 1+gab (еаеь) = 1+ (1/2) [gab (еаеь) ~ gab (—еаеь)\ +

 

 

 

 

 

 

 

 

 

+ (1/2) [gab (еаеь) +

gab (— еаеь)\,

 

 

 

 

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

что

СаЬ= 1 при г— >-оо. Очевидно, \

что в

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

которого е+ = е_ =

= Ze при переходе от G++ = G— к

G+_ = G_+, выражение,

стоящее

в

пер­

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

tfeb = v(0 { i-'-fg -x M (0 + * (0}.

(V.9)

где множитель % выделен для удобства.

 

 

Поскольку

 

 

 

 

Нш А ( К а ь ~

0 = 0,

(V.10)

ТО и

*—►00

 

 

 

 

 

lim

t*M(t) = lim

tsR(t) = 0.

 

/-►оо

t-*-oo

 

 

480