![](/user_photo/_userpic.png)
книги из ГПНТБ / Эккер, Г. Теория полностью ионизированной плазмы
.pdf§ 4. Микроскопические свойства кулоновской системы |
83 |
Подставляя (4.14) в (4.23), мы найдем выражение для сред ней потенциальной энергии t-й частицы s-конфигурации
W ’= <*»>.=
|
|
P<S+1) (“ Г1, |
d hTj. (4.65) |
= 2 |
Фи+ 2 j Фч |
P(S) (aIi |
|
з—1 |
k |
|
|
Далее, применив оператор Т г к обеим частям выражения (4.65) и подставив полученный результат в (4.64), получим
+ 2 J[' |
ftr]) |
- |
фцй "г}■ (4.66) |
P(s) ( Ч .......ч ) |
] |
Для вычисления второго члена в правой части уравнения надо найти явное выражение для отношения р<8+1>/.р<8), справедливое в рамках приближения парных корреляций.
Рассмотрим выражение для частной молекулярной функции распределения порядка s в виде *)
4 = П 4 |
4 ) 1 1 ( 1 |
+ Ч 7 ) П ( 1 + |
Ч ; * ) П . . . • <4 -67) |
i |
fe(2 ) |
ft(3) |
fc(4) |
В этом выражении величины sgj. . .k суть корреляционные функции системы s-частиц, характеризуемой функцией распределения Ps. При вычислениях, приведенных ниже в настоящей главе, мы используем предположение о том, что корреляционные функции sgt .. .k для любых значений s одинаковы. Последнее является хорошим приближением до тех пор, пока s <^ N.
Порядок корреляционной функции характеризуется числом индексов, совпадающих с индексами коррели рующих частиц. Символ к (q) указывает, что произведение следует брать по всем комбинациям индексов частиц
порядка q.
В приближении парных корреляций мы пренебрегаем всеми корреляциями выше второго порядка. Тогда выраже-1
1) См. примечание на стр. 73.
Ь*
84 Гл. 1, Равновесные состояния кулоновской системы
ние (4.67) можно переписать в виде
|
|
П р2(гь ri) |
/ \ - Ц / М п П К М |
ы - - й(2) — |
|
|
М2) |
П [р1м г * |
Отсюда следует, |
что |
I |
|
||
|
(• • ч rj) |
П ^2 (ГЬ Г;') |
P s +t |
i |
или
Ps+Pl(; - - \rj) = Л (о) пи +ей (*. о)].
г
Согласно соотношению (4.14), мы имеем
(4.68)
(4.69)
(4.70)
|
Р“>(...) |
' k |
Vh) |
P s (...) |
’ |
|
|
Р <2) (*Гу) вГг) = |
р 2 (Гу, гг) TVfclVg |
при |
/с ф q, |
(4.71) |
|||
Р (2) (hry, 9Г ;) = Р2(гу, Гг) |
(7Vft— 1) |
при k = q. |
|
||||
Поэтому, учитывая выражение |
(4.67), получаем, |
что |
|||||
при k^=q |
|
|
|
|
|
|
|
Р ‘2>(ftry, 9гг) = Р “>(йгу) Р '1’ (9гг) [1 + ghq(ftry, 9гг)], |
(4.72) |
||||||
а используя это соотношение, согласно (4.70), находим |
|||||||
Р<»+1>(..., ftry) |
Nh—vh |
Р (1> ( Ч Щ и + ^ Ч 'Л г ) ] . |
|
||||
|
Pft |
|
|
|
|
|
|
|
|
|
|
|
|
|
(4.73) |
Второй член в правой части уравнения (4.66) теперь можно записать в виде
ft
= 2 |
j р а ' |
( ^ Ы Фа d kTj. |
(4.74) |
ft |
|
|
|
При выводе этого уравнения в соответствии с определени ем приближения парных корреляций мы пренебрегли все ми произведениями g-функций. Следует заметить, что
§ 4. Микроскопические свойства кулоновской системы, |
85 |
оператор V; действует лишь на координаты i-й частицы s-конфигурации и, следовательно, не влияет на функции распределения первого порядка, входящие в соотноше ние (4.73), которые относятся к группе остальных частиц.
В предположении однородности и изотропности систе мы корреляционная функция ghq зависит только от абсолютной величины вектора г = йг;- — чг; — расстоя ния между t-й и /-й частицами. Это означает, что
|
|
|
V ighq |
= |
~ |
V} gkq. |
|
(4.75) |
|
Используя полученное равенство, находим |
|
||||||||
У f |
/ г |
|
|
\ |
, j |
ft |
Vi |
Nk —vft |
|
2 j J |
l V i "■ |
p u > (...) |
) h i d |
|
fj = — 2 j |
~Nh |
x |
||
h |
|
|
|
|
|
|
h |
|
|
|
x |
j |
Pa) (hTj) |
|
[ghq(«Г,, 4)1 d V |
(4.76) |
|||
Потенциальная |
энергия |
взаимодействия |
|
согласно ее |
определению, также зависит только от абсолютной вели чины вектора г. Таким образом, сместив систему коорди
нат в точку |
Qrh |
соотношение (4.76) можно представить |
|||
в виде |
|
|
|
|
|
|
|
|
Nh |
X |
|
|
|
|
|
|
|
X |
pay (r)_p(D (_r) |
|
|
|
|
f |
[<Ы(Г) v ghq(r)]dr. |
(4.77) |
|||
Из предположения об однородности системы следует, |
|||||
что Z*'1’ — постоянная величина. |
Поэтому |
|
|||
|
2 |
J |
|
= 0 |
(4.78) |
и, следовательно, |
|
|
|
|
|
|
|
© V ;lnP<S)= — |
|
|
(4.79) |
Просуммировав теперь результат |
(4.79) по всем |
частицам |
|||
s-конфигурации, найдем уравнение |
|
|
|||
|
BV In/><*>.-=-2 V, <*,)(■>, |
(4.80) |
86 Гл. 1. Равновесные состояния кулоновской системы
где
(4.81)
i — 1
есть оператор V в координатном пространстве s-конфигу рации.
Выразим теперь правую часть уравнения (4.80) через среднюю энергию частиц s-конфигурации
W - = т 2 |
' *» + |
2 |
) |
,(N-s) |
tT' + c |
||||
i, |
} |
3, |
Я |
|
(4.82)
Первый член в правой части этого выражения представля ет собой энергию взаимодействия частиц данной s-конфи гурации между собой, второй — среднюю энергию взаимо действия (N — s) оставшихся частиц системы с частицами s-конфигурации, а последний член C<-N~s)описывает среднюю энергию взаимодействия остальных (N — s) частиц друг с другом. Величина C(N~S) не зависит от координат частиц, образующих рассматриваемую s-конфигурацию.
Применим оператор V сначала к выражению (4.82), а затем к (4.65), после чего вычтем полученные уравнения одно из другого. В результате имеем
VW <S> = 2 |
V* ^ > <3,+ 2 j ^ i P^ Z > \\\\)ri)dhгя(4-83) |
|
|
i |
i , g, i |
|
|
ЯФi |
или |
с учетом (4.74) |
|
V(^)(4 |
^ .(ф г)<8,+ |
|
|
i |
|
+ |
2 |
j pa> (Ч-) Ф ^Я т ъ Г ч, S ') d hrj. (4.84) |
i, g , i q=£i
В соответствии с исходным допущением об однородности системы множитель Р<-1) по-прежнему предполагается постоянным по величине.
Здесь следует обратить особое внимание на то обстоя тельство, что во втором члене в правой части уравнения (4.84) в функцию ф(/]- входят координаты с индексами
§ 4. Микроскопические свойства кулоновской системы |
87 |
(q, j), тогда как в функции gmh координаты имеют индексы (i, /). Если бы набор индексов у координат в этих функци ях был одинаков, то в силу соотношений (4.76) и (4.77) второй член можно было бы опустить. Тем не менее рас сматриваемый член все же обращается в нуль, поскольку выполняется следующее равенство:
J <j>q i V i g m k C r i , h r j ) d h r j = - ^ f a q V q g M h { mr q , 4 -)dhi> (4.85)
Чтобы убедиться в этом, вспомним, что функция g зависит только от расстояния между частицами. Следовательно, суммирование по всем значениям i и q во втором члене правой части уравнения (4.84) дает нуль.
Итак, мы показали, что
(4.86)
Учитывая этот результат в (4.80), получаем простое соот ношение между общей молекулярной функцией распре деления s-ro порядка и средней потенциальной энергией частиц рассматриваемой s-конфигурации:
|
©V In РО) = |
- |
V< ф у. |
|
|
(4.87) |
Доказательство второго положения. Для вывода обоб |
||||||
щенного уравнения Больцмана |
— Пуассона применим опе |
|||||
ратор V; к выражению (4.65) и используем уравнение |
||||||
(4.78). В результате найдем |
|
|
|
|
|
|
v , < * ,г - з v ,+ 0 + 2 I |
о |
- |
в |
» |
) |
|
3 = 1 |
3= 5+ 1 |
|
|
|
|
|
Произведя скалярное умножение этого уравнения на оператор V, и принимая затем во внимание хорошо изве стные соотношения теории потенциала, получаем следую щее уравнение:
Ai (Фг)(8) = |
- 4л 2 |
’6 (*г7 - ?Гг) ehei ~ |
|
|
|
|
3=1 |
|
|
|
|
N |
|
|
|
— 4л |
2 |
j - ^ i r - 6 (ftrr - ?ri)eftM 4 - + |
|
|
N |
3 = s + l |
|
|
+ |
2 |
j |
Р<8+1)(..., ъ1})) ( X i<f,ij ) d krj. |
(4.89) |
( Vi Р<*>(...) |
|
?=s+t
Гл. 1. Равновесные состояния кулоновской системы
Обратим внимание на третий член, содержащийся в правой части уравнения (4.89), наличие которого обус ловлено тем обстоятельством, что функция распределения плотности пространственного заряда в данном случае
зависит |
от координаты |
пробной частицы |
hTj. |
Согласно |
|
(4.87), |
можно |
написать |
соотношение |
|
|
|
^ |
L = ^ |
Mp ( g g . y ^ " |
) . |
(4.90) |
Здесь A h— нормировочная постоянная. Исходя из условий нормировки
f Psdti ... drs= 1,
(4.91)
i Pi dvx= 1
v
и выражений (4.7), (4.8), а также (4.68) в пределах приме нимости приближения парных корреляций, можно полу чить условие
\ Р (1) (+;) Р{1) (Ч-) gqh (hrj, qrt) dhrj d?tt = 0. (4.92)
Подставляя соотношение (4.73) в (4.90) и учитывая, что Р (1>— константа, с помощью условия (4.92) находим выра жение
Аъ = |
|
(jVfe —vft) Vs |
|
(4.93) |
j « р [. |
(ф)<«+1)_<ф)<а> - |
|
||
|
d r а |
. d r 8 + 1 |
||
|
|
© |
|
|
Подстановка |
полученного результата |
в уравнение |
||
(4.89) дает |
обобщенное уравнение |
Больцмана — Пуас |
||
сона |
|
|
|
|
§ 4. |
Микроскопические свойства кулоновской |
системы |
89 |
|
|
s |
N |
|
|
Аг (<£г)(5) = |
— |
б.(йг7— 'ггг) ekeq— 4я 2 j |
6 (ftry— 9гг) х |
|
|
j=l |
j=«+l |
|
|
|
х |
j exp [ —(Ф>1а+1=(Ф)|г> j |
rfr, |
rfr,+1 + |
||||
|
* |
j(Vj*u) ( Vi exp [ - <Ф)'т ^ - <Ф>(8)] ) |
|
|||||
+ |
2 |
f |
r |
(ф)<т)_(ф)('8)-| |
|
" |
• |
|
. i = s+ |
i |
- - - - - J- -e- x-—p |
L |
|
J * 1 |
■■■ d r s + i |
||
|
|
|
|
|
|
|
|
(4.94) |
Заметим, |
что (4.94) |
содержит |
уже |
не |
одну |
функцию |
||
(фУ, |
а три различные |
функции |
(<£; )<S), |
(ф )<s>и (ф)<s+1>. |
Обобщенное уравнение Больцмана — Пуассона являет ся довольно сложным, и его решение можно рассмотреть только на основе метода суперпозиции средней потен циальной энергии парного взаимодействия частиц данной s-конфигурации, к изложению которого мы и перейдем.
Доказательство третьего положения. Попытаемся теперь представить среднюю потенциальную энергию частиц s-конфигурации в виде суперпозиции эффективных потен циальных энергий парного взаимодействия для данной конфигурации. Запишем следующее исходное выражение:
<4>>(S,= S ^ $ ( Ч - Л з), |
(4.95) |
4 2 ) |
|
где через к (2) обозначены все возможные парные комбина ции частиц i, j в пределах s-конфигурации.
Можно ожидать, что эффективная потенциальная энер
гия взаимодействующей пары зависит от набора v, обра зуемого частицами s-конфигурации, и, в частности, от числа s. Индекс (q, к) указывает на то, что эффективная потенциальная энергия будет также зависеть от сорта (q, к) двух взаимодействующих частиц s-конфигурации.
После того как мы учли сорт частиц (q, к), влияние системы на эффективную энергию взаимодействия частиц в сущности сводится к действию подсистемы остающихся 0V — s) частиц. Имея это в виду, заметим, что из общего
90 Гл. 1. Равновесные состояния кулоновской системы
приближенного соотношения |
|
N h t t N k- 1 |
(4.96) |
следует |
|
« v^<2). |
(4.97) |
Используя последнее приближенное равенство, а также выражение (4.95) и линеаризуя экспоненциальную функ цию (что, как мы покажем, вполне соответствует рас сматриваемому приближению парных корреляций), с по мощью уравнения (4.94) получаем
Л; |
2 |
V<£ $ |
(5r;> Ч ) = |
— 4 я 2 |
б (hi j — qn ) |
eqeh ■ |
|
|
i—1 |
|
j=i |
|
|
|
|
|
|
iV-e |
|
|
|
|
|
|
4л |
2 [«««ft (N h - ' ’k) Vs- 2 ^ (N h - V h) Vs |
(9ri, »r,)] |
||||
|
|
h |
|
1 = 1 |
|
|
|
|
|
|
j „ P[ |
(ф)<«+1 )_ (ф )(« |
dr.s41 |
|
|
|
|
|
0 |
|
|||
|
|
|
|
|
|
||
|
* |
j |
[(V; |
(,fry, |
?r;))/0] (Nh- v h) Vs d |
. (4.98) |
|
- |
2 |
|
|
(ф)<*+1)_(ф)<«) |
|
||
|
j=s+ l |
j exp { |
0 |
} dt\ ... |
c?rs+i |
|
Заметим, что первый член в правой части этого уравнения содержит знак суммы со штрихом, а в остальных членах он отсутствует. Последний член не содержит суммирова ния по индексу I, поскольку оператор V; исключает все члены, кроме тех, у которых I — i.
Интегрирование по частям в четвертом члене приводит к результату
*r (yVfe~eVft)F5]
j=^+l j exp | ' _ <ф?<*+ 0~ <Ф)<8>] d*i ••• *s+i
4jt |
N |
e h e q { N h — Vh)VS ^ v <pqk («гг, |
h i j ) б(9гг — h t j ) d h r j |
|||
~ ~ в ~ |
2 |
f |
г |
(ф)(«ы)_<ф)(«)-j |
" |
|
|
i-*+l |
3 eXP [ |
--------- 0------- |
J *1 |
dTs+l |
(4.99)
§ 4. Микроскопические свойства кулоновской системы |
91 |
из которого видно, что этот член взаимно сокращается со слагаемыми у = I, и = к, входящими в третий член правой части уравнения (4.98). Таким образом,
S S
А; 2 |
V(l,$ (,гг> hrj)= |
~ /in 2 |
б (йг1 — 5гг) е ^ |
— |
|
|||
j=l |
|
|
|
i=l |
|
|
|
|
|
|
|
4 я |
У] eqeh (Nh—v h) V s |
|
|
||
|
|
|
|
ft |
|
|
|
|
|
|
^ exp |
(ф)(8+1)_(ф)(8) dr, ... drs+l |
|
|
|||
|
|
|
|
в |
|
|
|
|
|
N - s |
s |
eqeh (A/~h — Vft) |
v 4>gu ftr,-, urt)/® |
. (4.100) |
|||
|
|
|
||||||
|
4л S |
2 |
r |
r |
|
rfrj ... |
dis+l |
|
|
k |
i=lu |
\ exp |
|
|
|||
|
|
1фФх1 |
J |
*- |
|
|
|
|
Эффективная энергия парного взаимодействия удовле |
||||||||
творяет тривиальному соотношению |
|
|
||||||
|
|
|
|
fhu — — |
|
?qu- |
|
(4.101) |
|
|
|
|
V 1(2) __ ®fc |
V J.<2> |
|
|
|
В |
нашем |
рассмотрении |
нормировочная |
постоянная |
в выражении (4.93) является произвольной величиной. Потребуем, чтобы
A h= |
(4.102) |
Перегруппируем также слагаемые суммы по I. Тогда, учи тывая соотношения (4.102) и (4.101), мы получим уравне ние (4.100) в виде
S |
|
S |
A, 2 ' 4 |
$ |
(?r«, Ч ) = - 4л 2 ' б (?г, - Ч ) eqeh - |
3 = 1 |
|
3 = 1 |
ft |
|
1=1 ft |
(4.103)
Ясно, что последовательность суммирования в третьем члене правой части уравнения (4.103) не имеет значения,
92 Гл. 1. Равновесные состояния кулоновской системы
Поэтому вместо суммирования по I и по соответствующему сорту частиц и мы сначала проведем суммирование, по индексу / и соответствующему ему сорту частиц к.
Далее учтем, что однородная и изотропная плазма должна быть нейтральной. Это требует выполнения усло вия
л
= 0. |
(4.104) |
h
Используя данное соотношение, второй член правой части уравнения (4.103) можно записать в виде
- 4 я |
+ |
|
(*-105) |
|
ft |
It |
1=1 |
где штрих у знака суммы означает, что в ней должен быть исключен заряд i-й частицы. Здесь также учтено, что вы полняется приближенное равенство v7 — 1 « vg.
Подставив (4.105) в уравнения (4.103), нетрудно свести их к системе уравнений того же типа для отдельных слагаемых
Дг |
(?rn V/) = — 4ne9eh6 (qn — krj) + |
|
N - s |
+ |
4л - ^ - + 4л ( 2 4 - -w " )~v</v ( 4 , ftry). (4.106) |
|
ft |
Ниже мы покажем, что данная система уравнений (для эффективной потенциальной энергии) обладает решениями типа Дебая — Хюккеля. А это означает, что третье поло жение будет доказано.
Напомним его: средняя потенциальная энергия частиц данной s-конфигурации может быть представлена в виде суперпозиции потенциалов отдельных взаимодействующих пар. Потенциальная энергия каждой взаимодействующей пары может быть найдена из обобщенного уравнения Больцмана — Пуассона (4.106). Она зависит от сорта двух
выбранных частиц и набора v групп частиц рассматрщ ваемой s-конфигурации.