УДК 533.951
О ШИРИНЕ ЛИНИИ ВРМБ В ПЛАЗМЕ С ДВУМЯ СОРТАМИ СТОЛКНОВИТЕЛЬНЫХ ИОНОВ
К. Ю. Вагин, К. Н. Овчинников, В. П. Силин
В сравнительно широком интервале интенсивностей накачки получена простая формула для ширины смещения частоты поля рассеянного излучения относительно частоты поля накачки при ВРМБ, когда нелинейный закон имеет корневую зависимость от превышения интенсивности над порогом параметрической неустойчивости. Ширина линии ВРМБ может существенно превышать декремент затухания ионного звука.
I ' ' * *
В работах [1, 2] рассматривалась ширина линии спектра ВРМБ в полностью ионизованной плазме с ионами двух сортов. При этом рассматривалась бесстолкновительная плазма. Было показано, что уже на пороге ВРМБ неустойчивости возможно возникновение сравнительно широкой линии [1]. В работе [2] на примере углерод-водородной плазмы показано как изменяется ширина линии ВРМБ при изменении степени неизотер-мичности и соотношения концентраций ионных компонент плазмы. В отличие от работ [1, 2] излагаемое ниже посвящено теории ширины ВРМБ линии в плазме, столкновения ионов которой играют определяющую роль.
В 1994 году Эпперлейн, Шорт и Саймон [3] показали, что в теории звуковых волн в плазме с легкими и тяжелыми ионами, в пределе сильно столкновительной плазмы, когда спектр и затухание волн отвечают гидродинамическому пределу, важную роль играют эффекты, обусловленные различием отношений заряда к массе для ионов различных сортов. В работе [3] авторы использовали следствия кинетического уравнения с обычным интегралом столкновений Фоккера-Планка-Ландау. Результаты работы [3] относятся к такому пределу, когда длина свободного пробега ионов мала по сравнению с длиной звуковых волн, а частота звука мала по сравнению с частотами столкновений
ионов. Эти результаты интересны для теории вынужденного рассеяния Манделыптама-Бриллюэна (ВРМБ), в которой проявляется новая, выявленная в работе [3] диссипация.
В настоящем сообщении, посвященном теории ширины линии ВРМБ в плазме, имеющей два сорта ионов, легких и тяжелых, рассматривается сильностолкновительный предел для ионов, в то время как электроны будут считаться бесстолкновительными. При этом нас будет интересовать не только тот гидродинамический предел, которому была посвящена работа [3], но также и предел частот звука, больших ионных частот столкновений. При этом столкновительная диссипация ионов будет доминирующей. Для такого случая в работе [4] с помощью интеграла столкновений Батнагара-Гросса-Крука [5] было получено сравнительно простое выражение для продольной диэлектрической постоянной, которое мы используем далее.
Вслед за работой [4] при рассмотрении ионно-звуковых волн будем, с одной стороны, полностью пренебрегать вкладом электрон-ионных столкновений в диэлектрическую постоянную. То есть будем считать, что вклад электронов в продольную диэлектрическую постоянную имеет вид
где Г£>е - дебаевский радиус электронов, Уте = у к,дТе/те - тепловая скорость электронов. В (1.1) принято, что скорость ионного звука мала по сравнению с тепловой скоростью электронов. С другой стороны, вклад ионов в диэлектрическую постоянную к) будем считать столкновительным. При этом для интересующего нас случая плазмы с двумя сортами ионов, главным диссипативным вкладом от ионных столкновений, вслед за работой [4], будем считать вклад, определяющийся разностью отношений зарядов и масс ионов. Тогда имеем
Выражение (1.2) не учитывает столкновения между ионами одного сорта, и входящая в (1.2) частота и0 включает только вклады от столкновений ионов разного сорта.
(1.1)
бе^ш, к) =
ш2шь + и2ьу% . щ и>1 - ш2ь ш2(ш2 + 1/д) ш и2 + г/ц '
+ *
(1.2)
= ^12 + ^21
(1.3)
где
4 /—е^ЭДЛ / М2
1/12 = з^ {Квту )1м1(м1 + м3у 1/21 = ТЩ** (ы)
где кв - постоянная Больцмана. В выражении (1.2) использованы обычные обозначения для ионных ленгмюровских частот ионов в плазме с двумя сортами ионов
4 «4,+.4 = 4+ (1.5,
и обозначение
^ = 47ге ммТЖЖ"' (1-6)
Заметим, что благодаря соотношению
f^-^V (1.7)
VAfi Mj k '
(Z1N1 + Z2N2y
диссипативный вклад в (1.2) пропорционален квадрату разности отношений зарядов ионов к их массе (ср. [4]). В формулах (1.1), (1.2) к и ш - волновой вектор и частота звуковой волны, Na, Za, Ма - плотность, кратность ионизации и масса ионов сорта а соответственно. Формула (1.2) получена в предположениях Т, = Т\ = Т2 и в пренебрежении столкновениями одинаковых ионов между собой.
В низкочастотном гидродинамическом пределе, когда в соответствии с работой [4] выполнено условие
ш < min {i/i, v2) < v0, (1.8)
где
va = vaX + va2 (1.9)
частота столкновений ионов сорта а, учитывающая как столкновения ионов разных сортов (1.4), так и столкновения ионов одного сорта
4 e*Z¿N„A
■ з ^(^Fvk- ( '
выражение (1.2) принимает вид 34
M^-SLfi-i"^)- (lid
Ш2 \ U0
В таком гидродинамическом пределе формулы (1.1) и (1.11) позволяют с помощью дисперсионного уравнения
6ее(ш,к) + 6ф,к) = 0 (1.12)
записать необходимые для нашего рассмотрения следующие выражения для частоты UJs
Ws = kvs, (1-13)
и декремента затухания 7s звука
75 = 7е + 7». (1.14)
Здесь
VS—^LTDe (I-15) скорость звука в низкочастотном пределе (1.8),
7e = y/^/8(kv2s/vTe) (1.16)
бесстолкновительный декремент затухания Ландау низкочастотного звука на электронах,
7,- = - ш1)/й>1 (1.17)
декремент затухания низкочастотного звука, обусловленный столкновениями ионов разных сортов.
В высокочастотном пределе [4], когда справедливо
ш » max {1/1, i/2} > vo, (1-18)
ионный вклад (1.2) в диэлектрическую постоянную имеет вид
4( l-i*^). (1-19)
В этом случае уравнение (1.12) дает для u>s и декремента затухания 75 звука
us = kvs, (1.20)
7s = 7e + 7.-, (1.21)
где
vs = ULrDe (1.22)
скорость звука в высокочастотном пределе (1.18),
7е = V^/8(b|/<;Te) (1.23)
бесстолкновительный декремент затухания Ландау высокочастотного звука на электро-
7,- = Ы2)К - (1-24)
декремент затухания высокочастотного звука, обусловленный столкновениями ионов разных сортов.
Подчеркнем, что как в низкочастотном (1.8), так и в высокочастотном пределе (1.18) декремент затухания звука мал по сравнению с его частотой
7в «ые- (1-25)
ВРМБ неустойчивость. Используем постановку задачи, отвечающую развитию ВРМБ неустойчивости во времени. В основу такого описания ВРМБ положим следующее дисперсионное уравнение [6, 7]:
1 1
+
6ei(uj,k) 6ее(ш,к) ^[k-kcvs]!2
(2.1)
4(к - к0)2{с2(к - к0)2 - (ш - ы0)2е4г(и; - и*, к - к0)}' где ш0 и ко - частота и волновой вектор падающей электромагнитной волны; ыо — ш и ко — к - частота и волновой вектор рассеянной волны; и> и к - частота и волновой вектор участвующей в ВРМБ звуковой моды; — е~Е0/теш0 - скорость осцилляций электрона в поле падающей волны; - поперечная диэлектрическая проницаемость, определяемая следующим выражением:
е,г(ш - ш0, к - к0) = 1 - , (1 - ) ,
(ш — и>0у \ ш — ш0/
где
= А + = 4 /д~ е4пеЛ ЗД + ЗД
(2.2)
частота столкновений электронов с ионами. В уравнение (2.1) входит комплексная частота ш —* и; + ¿7, где 7 - временной инкремент звуковой моды, а значит и связанной с ней рассеянной волны. Везде далее ш - вещественная величина.
Используя представление о частотной расстройке резонанса, в соответствии с работами [1, 2] знаменатель в правой части уравнения (2.1) можно записать в обычном виде
(ш - ш0)2Ке[егт{ш - ш0, к - к0)] - с2(к - к0)2 = 2ш0А, (2.3)
где Д характеризует величину частотной расстройки параметрического резонанса.
В дальнейшем будем рассматривать развитие ВРМБ неустойчивости во времени с инкрементом
7<ы. (2.4)
Тогда, подставляя выражения (1.1) для 6ее и (1.11) или (1.19) для ¿6,- в уравнение (2.1) и используя (2.3), получим следующую систему уравнений для определения Д и 7, справедливую как в низкочастотном (1.8), так и в высокочастотном (1.18) пределах:
ш2 А
— - 1 = /7£Л, , ,—:-гг, (2.5)
Д2 + (7 + 7я)
7 + 75 т 7 + 7 е
2-= НЕ Л 2 , -тт. (2.6)
Ш5 Д2 + (7 + 7ег
Здесь 7е = ^ег^Ье/^о ~ декремент затухания поперечной электромагнитной волны, обусловленный ее поглощением при электрон-ионных столкновениях,
= ш0 |[к — ко, уд]12
4(к - к0)2< 1 "
безразмерная интенсивность накачки. Уравнения (2.5), (2.6) записаны в линейном приближении по малым согласно неравенствам (1.25) и (2.4) параметрам 75/0*5 и 7/и.
Уравнения (2.5) и (2.6) определяют следующим образом расстройку резонанса
А = ^ + (2'8)
Соответственно возникает следующее уравнение для инкремента 7:
7 + 7 е
1е
(а+ 2(7+ 75)
2а4(7 + 75) ^5
= I. (2.9)
Отсюда в частном случае 7 = 0 имеем уравнение для границы области неустойчивости:
1)2 + 275/а* = /. (2.10)
275/^5
Уравнение (2.10) описывает ширину ВРМБ линии. Оно дает зависимось частоты ш на границе области ВРМБ неустойчивости от интенсивности накачки 1. В линейном приближении по малому параметру 75/0^5 <С 1 из (2.10) имеем на пороге неустойчивости (ср. [6]):
/ = /,* = 2—, (2.11)
ш = ш3. (2.12)
Эти формулы позволяют для частоты на границе области ВРМБ неустойчивости записать следующее соотношение:
иь± = ± Цэу/!/!^ ~ 1- (2-13)
Формула (2.13) применима в весьма широкой области значений интенсивности полей накачки
1гк<1 75- (2.14)
Следует отметить, что, согласно уравнению (2.9), инкремент ВРМБ неустойчивости 7 как функция интенсивности накачки I положителен в области частот, удовлетворяющих неравенству и>ь- <ш< и>ь+.
Ширина линии ВРМБ, согласно формуле (2.12), определяется соотношением
До; = шь+ - шь- = 2^1/1^-1. (2.15)
При больших превышениях порога (/ Ith) ширина линии растет как квадратный корень из интенсивности поля накачки.
Для того, чтобы уточнить область применимости формулы (2.13), остановимся на выражении для инкремента. Чтобы не загромождать наше изложение, рассмотрим инкремент ВРМБ неустойчивости при ш = ujth = u>s, когда уравнение (2.9) дает
7 = ^ |\/(7s ~ 7Е)2 + 27Eusl ~ (7s + 7я)} • (2-16)
В пределе достаточно высокой надпороговости имеем
7 = VW/2. (2-17)
Согласно (2.17) видно, что использованное нами условие малости инкремента ВРМБ неустойчивости по сравнению с частотой ионного звука (2.4) сводится в рассматриваемом случае к неравенству
I < ushE, (2.18)
которое при 7е > 7s может быть более существенным по сравнению с правым неравенством формулы (2.14).
В заключение этого раздела обсудим как может влиять частотная расстройка (2.8) на величину волнового вектора к участвующей в ВРМБ звуковой моды. Обычно (см., например, [6, 7]) для определения к используют уравнение (2.3) с правой частью, равной нулю. Покажем, что учет конечной величины Д приводит лишь к малым поправкам в величине волнового вектора к. Так, на пороге ВРМБ неустойчивости относительное влияние частотной расстройки (2.8) на величину волнового вектора будет характеризоваться малой величиной (7e/u>0)(7s/<^s) «С 1, представляющей собой произведение малых отношений декрементов затухания к частотам соответственно электромагнитной и звуковой волн. При значительном превышении порога неустойчивости величина частотной расстройки (2.8) максимальна вблизи границы области неустойчивости и составляет Д = (7^/75)Ды/2. Тогда с учетом (2.15) в рамках наиболее жесткого из условий (2.14), (2.18) относительное влияние частотной расстройки на величину волнового вектора ограничено сверху величиной
(ws/wo) • min \JieIis, 7e/7s] < 1- (2-19)
Действительно, при < 75 условие (2.19) выполняется автоматически. Если же 1е > 75, ТО (ыя/^о)\JieIis < {шв/^о)\JlE~H~e ~ ^/(г>Те/с)(-)е/шо) < 1. Таким образом частотная расстройка приводит к малым поправкам в величине волнового вектора к.
Условия применимости нашего рассмотрения ограничиваются условиями примени мости исходных формул (1.1) и (1.2). В этой связи приближение бесстолкновительного описания электронов, приводящее к формуле (1.1), налагает на величину волнового вектора условие снизу
> 1, (2.20)
где £е, — ьте/ие{ - длина свободного пробега электронов относительно их столкновений с ионами. В свою очередь, формула (1.2) получена в работе [4] в старшем порядке разложения общего выражения для 5с, по малым параметрам
Рм2
к%)т° <1, а = 1,2 (2.21)
Ш1 + I/2
и условия ее применимости в низкочастотной области имеют вид
kvTa
ш < min va, max-- <С 1. (2.22)
В высокочастотной области формула (1.2) применима, если
ш maxi/a, max—— С 1. (2.23)
ш
Формула (1.2) также не учитывает не связанный с разложением по параметру (2.21)
бесстолкновительный диссипативный вклад ионов в диэлектрическую постоянную
<224>
существенный в высокочастотной области при kvxQ ^ иа- Такой подход оправдан в условиях, когда описываемая последним слагаемым (1.2) диссипация, обусловленная столкновениями ионов разных сортов, превосходит бесстолкновительный диссипативный вклад каждого сорта ионов (2.23)
Вклад от столкновений ионов одного сорта возникает при учете слагаемых следующего порядка малости по параметру (2.21) в разложении общего выражения для 6б, [4]. Таким образом неравенства (2.20), (2.22), (2.23) и (2.25) задают условия применимости нашего подхода, основанного на использовании выражений (1.1) для электронного и (1.2) для ионного вкладов в диэлектрическую постоянную.
Подводя итог всему вышеизложенному, можно утверждать, что для случая полностью ионизованной плазмы с двумя сортами ионов, столкновения которых существенны, и с бесстолкновительными электронами установлен простой аналитический закон нелинейного роста ширины линии рассеянного излучения при ВРМБ, который может приводить к тому, что ширина линии оказывается значительно большей декремента затухания ионного звука.
Работа выполнена при частичной поддержке CRDF (грант RP1-2268), РФФИ (проект 02-02-16047), а также в рамках программы государственной поддержки ведущих научных школ (проект 00-15-96720).
ЛИТЕРАТУРА
[1] К u z о г a I. V., К о z 1 о v М. V., М с К i п s t г i е С. J., et al., Phys. Lett., A, 296(1), 54 (2002).
[2] В а г и н К. Ю., Овчинников К. Н., Силин В. П., У р ю п и н С. А. Квантовая электроника, (2002), в печати.
[3] Е р р е г 1 е i п Е. М., S h о г t R. W., and S i m о n A. Phys. Rev., E, 49(3), 2480 (1990).
[4] В а г и н К. Ю., Овчинников К. Н., Силин В. П., Краткие сообщения по физике ФИАН, N 5, 18 (2002).
[5] Александров А. Ф., Богданкевич Л. С., Рухадзе А. А. Основы электродинамики плазмы. М., Высшая школа, 1988, с. 54.
[6] Г о р б у н о в Л. М. ЖЭТФ, 55, 2298 (1968).
[7] С и л и н В. П. Параметрическое воздействие излучения большой мощности на плазму. М., Наука, 1973, с. 147.
Поступила в редакцию 3 апреля 2002 г.