k, H 0.
Подставим H из (IX.2.1) в (IX.2.2) и раскроем двойное векторное произведение
2D k2E − k k, E , c
или в компонентах
2 |
− k |
2 |
k k |
E |
0. |
c |
|
Эта система уравнений имеет нетривиальные решения, если ее детерминант обращается в ноль:
det |
2 |
− k |
2 |
k k |
0. |
c |
|
(IX.2.4)
(IX.2.5)
(IX.2.6)
Аналогичная ситуация возникает в кристаллооптике. Там уравнение (IX.2.6) называют уравнением Френеля.
Вычислим в холодной магнитоактивной плазме, т.е. будем считать, что kvTe. В этом случае можно пренебречь такими кинетическими эффектами, как затухание Ландау, и производить вычисления в гидродинамическом приближении.
Разобьем поляризацию на электронную и ионную части, P Pe Pi. Из уравнения (IX.1.3) для каждой из них имеем:
− i Pj 4 jj 4 ZjnjVj, |
(IX.2.7) |
откуда находим вектор поляризации
196
P ∑Pj; |
Pj 4 i |
Zjenj |
Vj. |
(IX.2.8) |
|
j
Здесь индекс j означает сорт частиц, j e, i. Zj означает знак заряда частицы. Мы будем
рассматривать |
простую |
|
|
электрон-протонную |
||||||||||
плазму: Zi 1, Ze |
−1. Vj |
|
_ коллективная скорость |
|||||||||||
частиц. |
|
|
|
|
|
|
|
|
|
|
|
|
||
Найдем скорости из уравнений движения |
|
|
|
|||||||||||
|
∂Vj |
|
Zje |
1 |
|
|
|
|
||||||
|
|
|
|
|
|
mj |
E |
c |
Vj, H , |
|
|
(IX.2.9) |
||
|
∂t |
|
|
|
|
|||||||||
или, в фурье-представлении, |
|
|
|
|||||||||||
− i Vj |
|
|
Zje |
|
E Hj Vj, h . |
|
|
(IX.2.10) |
||||||
|
mj |
|
|
ZjeH |
||||||||||
Здесь введены ларморовские частоты Hj |
|
|
||||||||||||
mjc |
||||||||||||||
и единичный вектор h H/H, направленный вдоль магнитного поля.
Пусть система координат выбрана так, что ось z направлена вдоль магнитного поля. Тогда систему (IX.2.10) можно переписать в компонентах так:
i Vjx |
Zje |
Ex HjVjy, |
(IX.2.11) |
mj |
|||
i Vjy |
Zje |
Ey − HjVjx, |
(IX.2.12) |
mj |
197
i Vjz |
Zje |
Ez. |
(IX.2.13) |
mj |
Уравнение (IX.2.13) содержит только z-компоненту скорости, а (IX.2.11) и (IX.2.12),
напротив, ее не содержат. Таким образом, она отщепляется. Выразим скорости через электрическое поле с помощью (IX.2.11) _ (IX.2.13) и подставим (IX.2.8). Используя условие квазинейтральности ni ne ≡ n, получаем:
Px i ∑ |
pj |
|
i Ey HjEx |
, |
(IX.2.14) |
||||
|
|
|
|
2 |
2 |
||||
j |
|
|
|
|
|
− Hj |
|
|
|
|
|
|
j |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Py i ∑ |
pj |
|
HjEy i Ex |
, |
(IX.2.15) |
||||
|
|
|
|
2 |
2 |
||||
j |
|
|
|
|
|
− Hj |
|
|
|
|
|
|
j |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Pz − ∑ |
pj |
|
|
2Ez. |
|
(IX.2.16) |
|||
|
|
|
|
||||||
|
|
|
|
j |
|
|
|
|
|
Пользуясь этими соотношениями и связью между P и E (IX.1.9), (IX.1.10), находим компоненты тензоров и :
хх хх − ∑ |
|
|
pj2 |
, |
xx yy |
(IX.2.17) |
|
2 |
2 |
||||
j |
|
− Hj |
|
|
||
|
|
|
|
|
|
|
198
1 −∑ |
|
|
|
pj2 |
|
|
, |
|
|
|
|
|
||||||||
|
|
2 |
2 |
|
|
|
|
|
||||||||||||
|
|
|
|
|
j |
|
|
|
− Hj |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
zz − ∑ |
|
pj2 |
; zz |
1 |
− |
|
∑ |
|
pj2 |
. |
(IX.2.18) |
|||||||||
|
2 |
|
|
2 |
||||||||||||||||
j |
|
|
|
|
|
|
|
|
|
|
|
|
j |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
В этом выражении обычно можно пренебречь |
|
|||||||||||||||||||
членом, пропорциональным |
pi2 , |
|
и пользоваться |
|
||||||||||||||||
выражением zz 1 |
− |
|
p2е |
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
xy − yx |
i |
∑ |
|
|
|
pj2 Hj |
|
, |
|
(IX.2.19) |
||||||||||
|
|
|
2 |
|
2 |
|
|
|
||||||||||||
|
|
|
|
|
|
|
j |
|
|
|
− Hj |
|
|
|||||||
xy xy − yx, |
|
|
|
|
|
|
|
|||||||||||||
zx |
zx |
xz |
zx |
|
|
|
(IX.2.20) |
|||||||||||||
zy |
zy |
yz |
yz 0. |
|
|
|
|
|||||||||||||
Введя |
|
xy ≡ ig, |
|
|
|
|
|
|
обозначения |
|
||||||||||
xx ≡ , zz ≡ , |
|
|
|
|
запишем |
тензор |
|
|||||||||||||
диэлектрической проницаемости в виде
199
|
|
|
ig |
0 |
|
|
|
|
−ig |
|
0 |
. |
(IX.2.21) |
00
Влитературе иногда также встречаются обозначения: ≡ ; ≡ .
Заметим, что все компоненты этого тензора не
зависят от волнового вектора.
Итак, теперь мы можем выразить все коэффициенты в дисперсионном уравнении (IX.2.6) через и k. После простых, но громоздких вычислений можно получить:
A |
|
kc |
4 B |
|
kc |
2 C 0. |
(IX.2.22) |
|
|
|
|
||||
|
|
|
|
|
|||
Здесь |
sin2 cos2 , B − 1 cos |
||||||
A k−2 k k |
|||||||
|
|
|
|
|
угол |
между |
|
C 2 − g2 , kH0 |
− |
||||||
направлениями |
|
волнового |
вектора и магнитного |
||||
поля.
Так как функции A, B и C, так же, как и тензор, не зависят от волнового вектора, решение биквадратного уравнения (IX.2.22) позволяет найти зависимость фазовой скорости /k от частоты. Уравнение (IX.2.22) биквадратное, поэтому имеются два решения для квадрата фазовой скорости, т.е. два типа волн.
В общем случае выражение для фазовой скорости очень громоздко, поэтому рассмотрим некоторые частные случаи.
200