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

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

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

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

Добавлен: 15.10.2024

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

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

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

представлен на рис. 21. Для малых к можно написать разло­ жение

Ф (il) = 1 — Л2----

л4 In г] + о (г]4);

 

Г (к)= К? 1 — к 2 ------ — е/с3 ------

е 2к * 1п к -

(17.10)

4яе

 

 

Отсюда нетрудно получить с помощью почленного обращения Фурье асимптотическое представление функции Г (R ):

Г (R) « 12е2Д- 9 — 20s3R~1 + О (e2R~s In R).

(17.11)

Это разложение по обратным степеням R является хорошим приближением, если R^>e и R^> 1. Интересно, что в отличие от

 

\ к °

 

 

 

 

 

\\

 

 

 

 

 

'

\\

 

 

 

 

£Ч __ ^

 

 

 

 

 

 

 

 

p-~=:r:7IZ

 

о

 

1

2

з

R

 

 

Рис. 21. Функция евГ (R) е~'.

 

приближения

Д ебая— Хюккеля,

когда Го= —(e/R) ехр(—R),

T(R)

при конечном

значении

е уменьшается

с увеличением R

не экспоненциально,

а степенным образом (~ R -6).

Рассмотрим

теперь предельный

случай е<§Д и представим

Г(/?)

в виде ряда по малому параметру е. Пусть

 

 

f0(R) = - u { R )

= - e / R ;

 

 

Гх (R) — T(R) — Г„ (R);

h (R) = f ( R ) - /0 (Я) =

= exp (— e/R) — l -f (e/R).

Тогда согласно общему выражению (17.7) имеем

т0(к) = Ш [ \ - п Ш Г 1;

Гх(К) = Г(к) — Го (к) = h (к) [1 - я/ (К)Г1 [1 - /о (к)]-1.

198


Используя тождество

[1 — п/0(к)][1 + лГ0(к)1 =

1,

 

получаем

 

 

 

 

 

Г1 (к) = /х (к) [1 +

«Г0 (/с)]2 [1 -

У (к)Г\-

 

К (к) = ri/i (/с) [1 + яГо («)J-

 

(17.12)

Поскольку f0 = —4лекг2 и /(/с)

определяется выражением

(17.8),

легко видеть, что функция

 

/с2)"1 = егос/4 (1 + к2),

 

Y (к) = [1 — <р Cn)] (1 +

 

если ек< 1 и У(/с) < (1 + /с2)-1

для любых

е и к. Для

е<С1

максимальное значение У (к)

приближенно

равно яе/2

и ряд

Тейлора для Fi (к) можно записать в виде

 

 

 

Гх (к) = h (к) [1 +

СО

У; (к).

(17.13)

«Г0 (к)]2 £

Этот ряд быстро сходится.

Подставляя У (к) в виде (17.12) в формулу (17.13), выпол­

няя

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

Фурье, получаем искомое выражение для

Г] (R) =Г] (1,2) через

бесконечную

последовательность

много­

кратных интегралов:

 

 

 

 

 

 

ОО

Г,°(1,

2)1

 

Г[0) (1, 2 ) - M l, 2) =

 

 

Гх(1, 2) = j

 

 

 

О

 

 

 

 

=

2n Jd3r0 (1 ,3 )/ (3,2) +

n2 J j <Ш4Г0 (1, 3) /х (3, 4) Г0 (4, 2);

 

Г(11) (1, 2) = n Jd3/1(l,

3)/х(2, 3 )+ . . .

(17.14)

Недописанные члены

для Г

все члены для Г)2) и т. д. со­

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

ных переменных.

Любой член

Г ,(1)

содержит

под

интегралом

произведение ( /+ 1) множителей

но R^>e

(не

обязательно

Рассмотрим область, где e<Cl,

/?> 1 ). В этом случае можно сделать приближение

 

 

h (г) «

е2/2г2.

 

 

(17.15)

Подставляя f\(r)

в выражение

(17.14) и интегрируя, получаем

Г(10) (R) ~

h (R) =“

[ ( # -

3) exp ( -

R) Et (R) -

 

О Д

 

 

 

 

 

 

 

- ( Я + 3)ехр(Я)£_(Д)],

 

 

(17.16)

где

 

 

 

 

 

 

 

El{R) = J

^ t ^ - d r ,

 

£_(/?) = J

txv(~ x)-dx.

199



При Д» 1 выражение (17.16) дает асимптотически Г)0* (R) « « 12е2//?6, что представляет собой первый член более общего

выражения (17.11). Оказывается,

что и при /?<1

можно полу­

чить разумные оценки для Г (R).

Действительно, в этом случае

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

(17.16)

порядка

е2(1 + 1п|/?|)

и по­

грешность, возникающая

от

подстановки f i

в

виде

(17.15),

имеет порядок e/R, тогда

как

Г(1) ~ е 3(1 + Д-1), а Г(‘)~ е г'4'2 для

i>2. Если R~>e, можно заменить Г, на Г{0) и воспользоваться

для Г 50)выражением (17.16). Погрешность, которая

при этом

получается, имеет порядок е при ^ > 1

и e/R при

1.

Для всех членов бесконечного

ряда

(17.14), которые содер­

жат интегрирование по двум или

большему числу переменных,

существенный вклад определяется областью гц вблизи единицы, даже если R мало. Следовательно, подстановка выражения (17.15) в эти интегралы приводит к хорошему приближению для всех значений R (предполагается, конечно, что е<С1). Од­ нако в этом разложении имеются два члена, каждый из кото­ рых содержит однократный интеграл. Для этих двух интегра­

лов существенная область интегрирования есть

г13~

г32~ г12,

так что приближение (17.15) для f i

при их вычислении

в обла­

сти R ^ e недопустимо. Необходимо

использовать

корректное

выражение

 

 

 

к (R) = exp (— e/R) — 1 -f- (e/R).

В

этой области /?<С1,

е<с 1. Используя точное выражение

для

fI

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

получим для первого

интегрального члена в выражении

для

Г |0) [см. формулу (17.14)]:

У(/?) = 2«|г/ЗГ 0(1, 3)М З, 2) ^ е 2

+ ф / — ) +

1

6

Vе /

где С=0,5772 — постоянная Эйлера;

Фх (х) = X

■Х2 + 1 +

In е + 2С ----

-f-

(17.17)

\8 ]•

ОЧт)+

 

+ — (х2— 2х — 1)ехр (— 1/х);

Фа ( х ) = — х

х2 +

+ ^ ) £ - {1/х) +

 

10

( 1 + т + ^

+ —

( 60х2 — 24х — 47------

exp (—1/х).

60

V

х

 

200


Хотя другие члены в выражении для Г (, !) порядка е3, член,

содержащийся в формуле (17.14), в действительности порядка е3/(/? + е) и должен быть учтен при очень малых R. Тогда

п I d3f, (1,3) /, (2,3) = 8*F (R/г),

(17.18)

где

 

 

1 °°

*+«

 

F (х) = — j dmfi (и)

j‘ dyyf1 (у).

 

Ох — и

Наконец, пренебрегая членами порядка R3 и е, найдем при

Л « 1 :

Г, (R) fx (R) = s2 Jin e — 0,0956 + R2J-i- In R + -i- In e—0,120б ]+

+ ф1 ( W т F ( R l e ) + | ф 2( R / г ) J .

(17.19)

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

При х<С 1

Ф1( х ) ^ х — — ;

Ф2 (х) яа ----

— ;

/

3

2V ' 2

10

/= » ж 4 1 п 2 —

--------- х + — х2.

w

 

3 2

4

Поэтому при R<^e выражение (17.19) с точностью до членов R2/e2 сводится к виду

I \ (R) — h(R)~ е2 [In е + 0,4953 + (R/2е)].

(17.20)

При е< 7?< 1

Фх (х) « In х — 0,0772 + ~ ( y In х + 0,7520^ ;

Ф2 (х )« In х — 0,9106 + - j (Inх + 0,9894^;

F (х) « л2/ 16х,

так что в этой области, пренебрегая членами порядка e2/Rz, R3 и членами более высоких порядков, получаем

Га (R) - Д (R) « в2 Jin R - 0,1728 + у [ у In

+ 1,3689] +

+ R2J-y In/?— 0,2724] J .

(17.21)

Итак, бинарная функция корреляции w (1,2), представляю­ щая наибольший физический интерес при вычислении термоди­

201