О выборе верхнего предела в интеграле столкновений для кулоновских частиц

УДК 533.9

О ВЫБОРЕ ВЕРХНЕГО ПРЕДЕЛА В ИНТЕГРАЛЕ СТОЛКНОВЕНИЙ ДЛЯ КУЛОНОВСКИХ ЧАСТИЦ

С. А. Майоров

В работе рассмотрена проблема выбора параметра обрезания кулоновского взаимодействия для различных физических систем - классической кулоновской и ионной плазмы, пылевой плазмы, компенсированных полупроводников, звездных систем, плазмы многозарядных ионов и кластеров в сильном лазерном поле.

Представленные в работе [1] результаты компьютерного моделирования столкновений в плазме вызвали оживленную дискуссию, в результате которой появилось еще некоторое количество мнений по этому вопросу. Поскольку в результате не появилось единой позиции, а актуальность задачи о выборе верхнего предела интегрирования только увеличилась в связи с новыми приложениями, то имеется необходимость в рассмотрении общих подходов к решению этой задачи.

Вопрос о частоте кулоновских столкновений в плазме, начиная с работ Ленгмюра [2], привлекал внимание многих исследователей и имеет богатую историю. Современная теория столкновений в идеальной плазме основывается на работе Ландау [3], в которой интеграл столкновений Больцмана для кулоновской системы упрощается путем разложения по малому параметру энергии взаимодействия частиц. Расходимость получаемого интеграла столкновений устранялась путем ограничения верхнего предела интегрирования по прицельному параметру дебаевским радиусом. Обоснованность обрезания кулоновских сечений на дебаевском радиусе мотивируется эффектом экранирования в плазме и слабой логарифмической зависимостью частоты столкновений от радиуса обрезания. Подход к учету кулоновских столкновений, предложенный Балеску и Леннардом, основывается на понятии динамической поляризуемости плазмы. Но при этом значение диэлектрической проницаемости в области малых волновых чисел соответствует обрезанию на дебаевском радиусе. Таким образом, вводится эффективный радиус взаимодействия и процедура обрезания в неявном виде сохраняется.

На основе полученных в вычислительных экспериментах результатов можно говорить о большем, чем это обычно принято, значении для анализа процессов в плазме межчастичного расстояния между ионами. Этот факт уже отмечался в работе [4], где исследовалось затухание автокорреляционных функций микрополя. Недооценка этой характерной длины традиционно идет от кинетических теорий, построенных для газов систем с короткодействующим потенциалом. В кинетической теории для газов расстояние между частицами не играет никакой роли: имеет значение лишь размер атома и длина свободного пробега. А в плазме эта длина определяет характерный размер неоднородности электрического поля. Другая характерная длина в плазме, дебаевский радиус экранирования, часто подменяет среднее межчастичное расстояние. Поэтому изучение столкновений актуально не только с практической целью уточнения кинет!; ческих коэффициентов, но важно и для прояснения физической сущности кинетически теорий.

Подход Власова [5] позволяет естественным образом избежать проблемы расход! г мости, но его уравнение имеет дело с усредненной плотностью заряда, тем самым вводится масштаб пространственного усреднения, больший среднего межчастичного расстояния. При критике подхода на основе уравнения Больцмана для анализа плазмы Власов утверждает, что для плазмы учет последовательных парных взаимодействий совершенно не отражает реальности из-за невыполнимости условий парности столкновений. В то же время, отмечая, что в интеграле столкновений Больцмана надо учитывать столкновения с прицельным параметром до значений порядка половины межчастичного, он пишет, что существенную роль должны играть силы взаимодействия на расстояниях, больших расстояния между частицами, и действие этих сил не может быть учтено обычной схемой кинетического уравнения. Противоречия здесь нет. Это означает лишь то, что учет далеких взаимодействий надо проводить не по бинарной теории.

В работах Персико, Ленгмюра, Ландау, Давыдова, Дрювестейна, Спитцера и др. (обзор см. в [5]) исследовались столкновения в плазме, и в большинстве случаев при расчете частоты столкновений в качестве верхнего значения параметра обрезания предлагалось использовать радиус Дебая. Среднее межчастичное расстояние в качестве параметра обрезания выбирали Давыдов, Коулинг и Чандрасекар, который исследовал динамику

кает из-за неправильного применения решения задачи двух тел для описания далеких

расстояния. В настоящее время, следуя работам Ленгмюра, Ландау, Спитцера и др., в

т» иоттттг\^/,от/оп ттлттогп тт ттггл тч о /"• лг/^гт-т»"» »лргпт п/лтттт» x* и\/<11их ни ^/ислч/дишу/ч/ 1 ь

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

качестве верхнего предела обрезания используется дебаевский радиус. Борьба мнений здесь сосредоточена вокруг принципиально важных вопросов: 1) насколько правильно с помощью парных столкновений можно учесть множественность и немгновенность столкновений в плазме; 2) насколько важны далекие столкновения. Сторонники критического отношения к попыткам использования интеграла столкновений Больцмана для описания плазмы (Власов, Чандрасекар, Кога и др.) говорят о возможности учета столкновений в бинарном приближении только для прицельных параметров, меньших среднего межчастичного расстояния. С другой стороны, учитывая столкновения с прицельными параметрами до величины дебаевского радиуса, делается попытка учета далеких столкновений, которые безусловно важны. Но пока нет системы макроскопических уравнений, удовлетворительно описывающих динамику плазмы во всех режимах, то различные подходы - это методы, разрабатываемые для решения конкретной задачи и не имеющие общности [6].

Н. Н. Боголюбов в предисловии к работе [7] пишет, что в уравнении Больцмана имеется внутреннее противоречие. С одной стороны, делается предположение о стохастич-ности (гипотеза хаоса - 81оззгаЫап8а1г), при котором движение молекул трактуется как случайный процесс и вводится в рассмотрение статистический механизм бинарных соударений. С другой стороны, входящие в уравнение сечения получаются из решения уравнений динамики. В наибольшей мере это противоречие должно проявляться в системе из частиц с дальнодействующим потенциалом - плазме. Можно предположить, что причиной расходимости является нарушение критериев применимости уравнения Больцмана. В этом случае для устранения расходимости нет необходимости в привлечении эффекта экранирования.

Кога [6] развивает подход к устранению расходимости интеграла столкновений, основанный на выделении у каждой частицы ближайших соседей и учета взаимодействия с ними. Его трактовка причины расходимости интеграла столкновений Больцмана для плазмы представляется достаточно убедительной. Но представлены только качественные соображения, не подтвержденные какими-либо расчетами.

Вопрос о влиянии эффекта дебаевского экранирования возникает при исследовании практически важной задачи определения длины рассеяния в структурах без экрани-

плпотттгсг ттоттгмт»«,ог\ x» т/т *ттлттг'тгг\г»х1 1 ттттт iv ттг\тттгтттчптэг»лтттагт/ о v т^оттт тжтт/~»т^тто лгтпхюг\м/ттоотго ии/и^уНАУХУУ^ и ич/хушъиъиру/хллиитл ич/л^ и^у. хуъ.и.и инч/х м"1 ^ А V«»

что без экранирования длина свободного пробега равна нулю [8]. В то же время ясно, что это не так, ведь среднее поле, вычисленное Хольцмарком для системы случайно расположенных в пространстве тел, взаимодействующих по закону Кулона, является

конечной величиной.

Работа [9], опубликованная В.И. Коганом в 1960 году, содержала попытку устранить указанные противоречия путем прямого расчета влияния микрополя (см. также [10]). Полученный им результат интересен для понимания природы расходимости кулонов-ских столкновений и имеет некоторые важные приложения [11]. В упомянутой дискуссии ([12], см. также [13]) прозвучала мысль, что "физически прозрачная работа [8] ввиду достаточного объема заложенной в ней исходной информации имеет не меньше оснований считаться анализом из первопринципов, чем компьютерное моделирование [6]". Все это приводит к необходимости дополнительного анализа работы [9].

Исходная задача в [9] - это рассеяние пробной частицы с массой, зарядом и скоростью то, <?о, ^о в газе из N полевых частиц с массой и зарядом тп\, ^, беспорядочно распределенных по объему V и имеющих максвелловское распределение по скоростям с наивероятнейшей скоростью Полагалось также, что полевые частицы движутся по прямолинейным траекториям. Показано, что квадрат поперечного приращения импульса пробной частицы пропорционален корреляции электрического поля в моменты времени £ и I + т:

<

о

Здесь угловые скобки означают статистическое усреднение, Е±^) - проекция мгновенной напряженности микрополя на плоскость, нормальную к скорости пробной частицы. Корреляции микрополя в (1) в силу независимости частиц равны

и выражаются через коэффициент корреляции

Кг± =

где N - число полевых частиц, Е\±(1) - поле одной (любой) полевой частицы в месте нахождения пробной, усреднение в (2) проводится по максвелловскому распределению полевых частиц /(г>) и объему V.

Вычисление коэффициента корреляции дает следующий результат:

^-Ш^Ч-ЗМ1--'©]-1

Очевидно, что подстановка (4) в (2) и (1) приводит к бесконечному значению приращения поперечного импульса за любой конечный промежуток времени Таким образом, строгое рассмотрение в рамках исходных предположений приводит ко вполне ожидаемому результату - бесконечности. Если же в качестве нижнего предела интегрирования в (1) выбрать некоторую величину тТО|П, то получаем:

/(л л2\ ^ЯоЯ2^

\/ку0

Ч-ЗЬМ)-']^

В предельном случае неподвижных полевых частиц «1=0 получается

уоу тп

1 тпт

а в предельном случае неподвижной пробной частицы ь0 = 0 получается

С точностью до вида логарифма выражения

= дод^тр + тг)

Т~тп1П — — 1 ; V /

т07711^0!

где г>ог = (|г>о — В этом случае выражение (7) совпадает со строгим результатом, полученным другим способом для соответствующего случая неподвижных полевых частиц в [10].

В работе [9] не представлено аргументов в пользу выбора т^ в качестве нижнего предела интегрирования в выражении (1). Аналогия с корректным обрезанием прицельных параметров на соответствующей длине Ландау, которая приводит к точному ответу при вычислении вклада от близких столкновений, явно не проходит из-за различия физического смысла времени между корреляциями при вычислении интегралов (3) и временем сильного соударения. При любом ненулевом значении параметра обрезания убираются корреляции между близкими по времени точками для всех траекторий - как проходящих вблизи пробной частицы, так и для далеких пролетов. Для далеких пролетов при таком "обрезании" устраняются корреляции между полями в близкие моменты

времени (из-за спадающего характера зависимости от т они дают наибольший вклад в интеграл).

Логически правильнее было бы удаление из интегрирования близких траекторий (т.е. интегрирование по прицельным параметрам) или модификация потенциала взаимодействия на близких расстояниях. Но в этих случаях коэффициент корреляции зависит от частицы и равенство (2) неверно. На эту тему можно привести различные качественные соображения, в частности, такое. Если в качестве нижнего предела выбрать время, которое пробная частица за рассматриваемый промежуток времени проводит в области непрямолинейности траектории

г3

Tmin = rA,t = A i-§,

ri

где г0 = , г,- = (V/iV)1/3, то получаем значение кулоновского логарифма,

т0т iu01

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

Естественно, что флуктуации средней плотности полевых частиц приводят к дополнительному влиянию на частоту столкновений. Рассмотрим пример, когда влияние флуктуации велико [11]. В экспериментах по исследованию плазмы, образуемой фокусировкой сверхмощных и сверхкоротких импульсов лазерного излучения в газе, происходит сверхбыстрая фотоионизация газа. Периодичность движения электронов в сильном лазерном поле и малый размер пятна фокусировки могут менять характер электрон-ионных столкновений. Рассмотрим следующую модель. Пусть в начальный момент имеется полностью ионизованная плазма с плотностью ионов N, и электронов Ne = zN, ze - заряд ионов, — е - заряд электрона. Температуру ионов будем считать равной первоначальной температуре атомов, которая близка к комнатной температуре, а начальное распределение электронов будем полагать максвелловским с некоторой температурой Те, определяемой превышением поглощенной энергии квантов над потенциалом ионизации. Рассмотрим процесс воздействия на плазму электрического поля плоско поляризованной волны E(f) = (.£^,0,0), где

Ex

При отсутствии столкновений скорость и координаты электрона будут равны:

v(t) = v(0) + \е sina;£,

г(t) = г(0) + v(0)i + тЕ cos wt,

где векторы = —еЕо/тпш, ге = еЕ0/то;2 определяют скорость и амплитуду осцилля ций, Е0 = Е(0). Электрон-ионные столкновения приводят к хаотизации направленного движения во внешнем электрическом поле (1). Если скорость осцилляции электронов Ve значительно больше их тепловой скорости Ут = у/т] т, то частота столкновении электронов с ионами равна [16]:

16ze4Ne ( , VE\, к, и = ' ~

V Лт,п

Мощность столкновительного нагрева электронов в этом случае определена в работе [16] и равна:

цг = 8г2е4ЛГе7У,- Л . , УЕ\Л к.

тУЕ \ 2

Согласно [16], значения кт{п и ктах определяются обратным дебаевским радиусом и минимальным прицельным параметром, определяемым из условия применимости классической механики или теории возмущений для частицы, имеющей скорость УЕ. Куло-новский логарифм определяется следующим равенством:

Л = 1п^ = 1п-^-.

"-тт Ртхп

Логарифм отношения скорости осцилляции к тепловой скорости в (11) возник из-за расходимости частоты столкновений при обращении скорости направленного движения в ноль и обрезания столкновений со скоростями меньше тепловой. В случае круговой поляризации поля модуль скорости постоянен и такой логарифм не возникает [16]. В работах [17] для мощности столкновительного нагрева в сильном поле плоской волны получено похожее выражение:

тд/ 8г2е4ад- У. ктахУт /ЛА.

ж =

где значение ктах определяется из условия применимости классической механики ктах =

т

••"II '*•

В работе [18] для мощности столкновительного нагрева в сильном поле плоской волны получено приближение, в котором частота столкновений умножается на фактор, учитывающий отличие осциллятор ной скорости от тепловой:

е2Е2

v = vei

у2 \ -3/2

1 + Щ) ' - ^

4л/2^е4ЛГе

Где Л - обычное для плазмы значение кулоновского логарифма.

При воздействии на плазму электрического поля волны с круговой поляризацией Е(<) = (Ех, Еу, 0), где Ех{{) — ^ЕаСоьизЬ, Еу{1) = smшt, в отсутствие столкнове-

ний электроны будут двигаться по окружности радиуса ге/у/2 с постоянной скоростью Уе/\/2. В приближении мгновенных парных столкновений учитывается взаимодействие электронов с неподвижным ионом и не учитывается их взаимодействие между собой [19]. Тогда средняя сила "трения", действующая на электрон со стороны ионов, имеющих прицельный параметр меньше ртах, равна [19]:

_ 47гг2е47У,- А

где кулоновский логарифм Л = 1п ~ ln ¿^i, р± = ze2/mV2 - значение при-

цельного параметра, при котором частица отклоняется на прямой угол. Отметим, что рх здесь не является нижним пределом обрезания, т.к. учитываются все столкновения с р < Ртах Мощность столкновительного нагрева электронов получается умножением силы "трения" на скорость Ve!\J2 и плотность электронов:

w=4

т\е

Аппроксимация кулоновского логарифма. Кулоновский логарифм Л определяется через верхнюю и нижнюю границы прицельных параметров, учитываемых в интеграле столкновений Ландау [3], следующим образом: Л = \пртах/pmin. При учете всех возможных прицельных параметров 0 < р < +оо частота кулоновских столкновений логарифмически расходится как из-за близких, так и из-за далеких столкновений. Расходимость на нижнем пределе обусловлена нарушением предела применимости линейного разложения интеграла столкновений для близких соударений и устраняется следующим выбором минимального прицельного параметра:

ртхп = 1Мшах = Р± = ге2/тУ2,

что соответствует точному решению. Для больших температур, когда р± становится меньше длины волны де Бройля А = /г/ти, используется квантовый предел рт1П = А/4л" (наиболее подробное изложение см. в [10]). Устранение расходимости на верхнем пределе также производится путем выбора конечного значения ртах, отвечающего физически разумному механизму устранения далеких столкновений. Однако, сам выбор ртах не очевиден. Первоначально, большинство авторов полагало необходимым в качестве верхнего предела интегрирования использовать межчастичное расстояние:

о = П'-1/3

Углах — 5

но после вывода Ландау кинетического уравнения для плазмы общепринято использование радиуса Дебая: ртах = Го = {Те//^'ке2Ме)1^2.

Для кулоновского логарифма в плазме часто используется аппроксимация [20]:

А = 23.4 - (1/2)1пАГе + (3/2) 1п Те, при Те < 50 эВ,

А = 25.3 - (1/2)1пАГе + 1пТе, при Те > 50 эВ,

которая совпадает с предложенной Л. Спитцером [15] только в случае однозарядных ионов. Пренебрежение зависимостью от заряда ионов г оправдано для горячей плазмы ионов с небольшим зарядом ионов. В случае же холодной многозарядной плазмы с г >> 1, представляющей интерес для рентгеновских лазеров, отличие заряда ионов от единицы может приводить к существенному уменьшению кулоновского логарифма. Поэтому, для плазмы многозарядных ионов при оценке электрон-ионных столкновений следует полагать значение кулоновского логарифма равным:

Аг =23.4-

Аг = 25.3 - (1/2)1пАГе + 1пГе -\пг, при Те > 50г2 эВ.

Иттлпттп ТТ/.ТТЛТТТ шт/чт/.п ПТТ«ЛХТ,ЯТТТТЛ „„„ .т г, . п -----1---- А ( Л /0\ 1_ [ 1 I П ! А - С \

яьогдо, отралспис для и ии1 а^ифма л — у± / ) 11ц1тз/1« и /

через показатель неидеальности 6 = (z-\-l)z3e6Ne/T^ — 2е6Лге/Те3 при 2 = 1. Эта аппроксимация получена при учете экранирования как электронами, так и ионами, что справедливо только для плазмы частиц с близкими массами положительно и отрицательно

заряженных частиц (электрон-нозитронная, ионная плазма). Именно эта аппроксимация использовалась в работе [13] для анализа полученных в [1] результатов численного моделирования кулоновских столкновений электронов с бесконечно тяжелыми ионами. В результате было занижено полученное в расчетах [1] отклонение от теории. Следует отметить, что использованная в [1] методика численного моделирования, примененная для моделирования прямолинейно движущихся частиц, привела к результатам, точно совпадающим с аналитическими результатами [9, 10]. Следовательно, численные результаты [1] связаны не с методикой моделирования, как полагается в [13], а с немгновенным и множественным характером кулоновских столкновений.

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

Рассмотрим флуктуации плотности атомов в газе. Флуктуации плотности атомов ограничиваются из-за столкновений частиц. Величина флуктуации в реальном газе на чинает отличаться от флуктуации в идеальном газе, начиная с размеров порядка длины свободного пробега атомов Аа = 1 /сгаЛГа, где аа = -ксРа - газокинетическое сечение, Ма -плотность атомов, ¿а - диаметр атома. Характерное время релаксации та флуктуации в газе с температурой Та и массой атомов М определяется тепловой скоростью атомов (скоростью звука) с3 = (уТа/тпа)1^2, где 7 = 5/3 - показатель адиабаты идеального газа:

Та = К/Сз = 1/<7а^аС,.

Характерное время та, релаксации газовых флуктуаций плотности размера Аа в плазме с температурой электронов Те и зарядом ионов г определяется тепловой энергией приходящейся на один ион. Оно имеет величину порядка та, = Аа/ср; = 1/<та №аСр1, где среднемассовая скорость звука в плазме ср/ = [7(Та + гТе)/(М + гт)]1'2. За время действия сверхсильного и сверхкороткого лазерного импульса длительностью т>аз флуктуации плотности атомов в газе могут не успеть релаксировать ни за счет теплового движения атомов, ни за счет плазменных колебаний, т.к. обычно в экспериментах выполнено условие:

Т1а* « Таг « тс

а •

Таким образом, флуктуации плотности атомов в газе после сверхбыстрой фотоионизации становятся флуктуациями плотности ионов. При прямолинейном движении электронов под действием внешней силы (электрического поля лазерного излучения) влияние этих флуктуаций не будет экранировано. Эти флуктуации плотности ионов во время действия лазерного импульса могут увеличивать частоту электрон-ионных столкновений, если они значительно больше плазменных флуктуаций пространственного заряда (радиуса Дебая).

Сила динамического трения при прямолинейном движении пробной частицы с постоянной скоростью среди неподвижных полевых частиц логарифмически зависит от времени прошедшего с начала движения [9, 10]:

где величина гт,п определяется, как обычно, из условия применимости теории возмущений. Если прямолинейность движения обеспечивает внешняя сила, то частота столкновений будет определяться именно длиной участка прямолинейности движения. Радиус Дебая в этом случае не ограничивает действие кулоновских сил и не является верхним пределом возможных прицельных параметров. Поэтому, при учете движения в сильном поле (10), будем определять область возможных прицельных параметров, исходя из дли ны участка прямолинейности и среднего по полупериоду колебания значения квадрата

Кулоновский логарифм в этом случае определяется только частотой и напряженностью поля (10). Обозначив его Ае, получаем

^ =

4тг2г2е4Лг,-, Д*

тпу2 1п тт:,

скорости:

Ртах = 2гЕ = 2еЕ0/тш2, рт{п = 2ге2/тУ£.

АЕ = Ртах/Ртхп = ЩеЕо /2ГП2Ш4).

Мощность столкновительного нагрева электронов можно оценить, как произведение силы динамического трения на скорость и на плотность электронов:

\

\У =

При вычислении силы трения использовалось значение средней по полупериоду лазерного поля скорости (V) = Ve/у/2. Результаты компьютерного моделирования [21] показали хорошее совпадение именно с формулой (24), а не с аппроксимациями (12), (14) или (15).

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

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

Автор выражает благодарность коллегам в Институте общей физики им. А. М. Прохорова РАН за интерес к работе, и особенно В. И. Когану, А. А. Рухадзе и С. И. Яковленко за дискуссии. Автор благодарит также Российский фонд фундаментальных исследований (проект 02-02-16439) за финансовую поддержку работы.

ЛИТЕРАТУРА

[1] М а й о р о в С. А. Краткие сообщения по физике ФИ АН, N 9-10, 99 (1997).

[2] Langmuire L. Ргос. Nat. Acad., 14, 627 (1928).

[3] Л а н д а у Л. Д. ЖЭТФ, 7, 203 (1937); Phys. Z. der Sow. Union, 10, 154 (1936).

[4] M а й о p о в С. А., Ткачев А. Н., Яковленко С. И. ДАН СССР, 290, N 1, 106 (1988).

[5] В л а с о в А. А. Теория многих частиц. М.-Л., ГИТТЛ, 1950.

ífil К ora Т, Введение в кинетическою теорию стохастических процессов в газах, М., Наука, 1983.

[7] Боголюбов Н. Н. Проблемы динамической теории в статистической физике. М., Гостехиздат, 1946.

[8 [9 10 И

12

13

14

15

16

17

18

19

20

21 22

Форрестер А. Т. Интенсивные ионные пучки. М., Мир, 1992. Коган В. И. ДАН СССР, 135, N 6, 1374 (1960).

Сивухин Д. В. Вопросы теории плазмы, вып. 4, М., Госатомиздат, 1964. Майоров С. А. Краткие сообщения по физике ФИАН, N 7, 25 (1999); Физика плазмы, 27, N 4, 2001. Коган В. И. Частное сообщение, 1997.

Яковленко С. И. Краткие сообщения по физике ФИАН, N 7, 30 (1998). Cohen R. S., Spitzer L., R о u t 1 у P. Phys. Rev., 80, N 2, 230 (1950). Спитцер JI. Физика полностью ионизованного газа. М., ИЛ, 1957. Силин В. П. ЖЭТФ, 47, 2254 (1964).

Jones R. D. and Lee К. Phys. Fluids., 25(12), 2307 (1982). S h i e s s i n g e r L. and Wright J. Phys. Rev. A, 20, N 5, 1934 (1979). Трубников Б. А. В сб. Вопросы теории плазмы, вып. 1, М., Госатомиздат, 1963.

Брагинский С. И. В сб. Вопросы теории плазмы, т. 1, М., Госатомиздат, с.183, 1963.

Костюков И. Ю. Письма в ЖЭТФ, 73, вып. 7, 8, 438 (2001). Гордиенко С. Н. Физика плазмы, 26, N 6, 556 (2000).

Институт общей физики

им. А. М. Прохорова РАН Поступила в редакцию 6 февраля 2004 г.

\