Физика плазмы, 2019, T. 45, № 4, стр. 346-352

О роли соударений с нейтралами в процессе модуляционного возбуждения пылевых звуковых возмущений в запыленной ионосфере

Н. Д. Борисов a, С. И. Копнин bc, Т. И. Морозова bd, С. И. Попель bc*

a Институт Земного магнетизма, ионосферы и распространения радиоволн
Москва, Троицк, Россия

b Институт космических исследований РАН
Москва, Россия

c Национальный исследовательский университет “Высшая школа экономики”
Москва, Россия

d Московский физико-технический институт (государственный университет)
Долгопрудный, Россия

* E-mail: popel@iki.rssi.ru

Поступила в редакцию 23.08.2018
После доработки 30.10.2018
Принята к публикации 30.10.2018

Полный текст (PDF)

Аннотация

Проведена оценка роли неупругих столковений электронов и ионов с нейтралами при развитии модуляционной неустойчивости электромагнитных волн с участием пылевых звуковых возмущений в плазме запыленной ионосферы, а также влияния столкновений электронов, ионов и пылевых частиц с нейтралами на проявления модуляционного взаимодействия в запыленной ионосфере. Показано, что влияние столкновений электронов и ионов с нейтралами при этом обычно менее существенно, чем влияние столкновений пылевых частиц с нейтралами. Продемонстрировано, что влияние модуляционной неустойчивости на распространение электромагнитных волн в плазме запыленной ионосферы наиболее существенно на высотах 100–120 км. Модуляционное взаимодействие в запыленной ионосфере важно для объяснения таких эффектов, как наземные наблюдения низкочастотных ионосферных радиошумов с частотами ниже 50 Гц, генерация инфразвуковых волн в ионосфере и возможность их детектирования у поверхности Земли, усиление интенсивности зеленого свечения ночного неба на длине волны 557.7 нм от слоя нижней ионосферы, расположенного на высотах 110–120 км, модуляционное возбуждение неоднородностей концентраций электронов и ионов в ионосфере на высотах 100–120 км. Объяснен факт отсутствия наблюдений низкочастотных ионосферных радиошумов во время таких явлений, как серебристые облака и полярные мезосферные радиоотражения, обусловленных присутствием пылевой плазмы на высотах 80–95 км. Показано, что данный факт связан с тем, что модуляционные процессы на этих высотах подавлены.

1. ВВЕДЕНИЕ

Интерес к описанию пылевой плазмы в ионосфере резко возрос в 2000-е годы, что было связано, во-первых, с разработкой к этому времени методов исследования пылевой (комплексной) плазмы, в том числе, и в природных системах [1, 3, 4]. Во-вторых, интерес к данной тематике обусловлен возможной связью плазменно-пылевых процессов в ионосфере с климатическими изменениями на Земле и, в частности, с процессами глобального потепления.

Пылевая ионосферная плазма состоит из электронов, нескольких сортов ионов, нейтралов и массивных отрицательно или положительно заряженных пылевых частиц и, как считают, присутствует на высотах от 80 до 120 км. Важным источником пылевых частиц на этих высотах являются метеорные потоки (например, Персеиды, Леониды, Геминиды, Ориониды и т. д.), порожденные выбросами из ядер короткопериодических комет. Частицы попадают в ионосферу в результате бомбардировки Земли метеоритами с их последующим сгоранием на высотах 80–120 км и конденсацией метеорного вещества [4]. Концентрация таких частиц варьируется в пределах 10–1000 см–3 в зависимости от сезона и метеоритной активности. На высотах 80–85 км в летней полярной ионосфере присутствуют тонкие (по сравнению со шкалой высот атмосферы) пылевые слои, известные как серебристые облака, а на высотах около 90–95 км аналогичные пылевые слои приводят к так называемым полярным мезосферным радиоотражениям – PMSE (Polar Mesosphere Summer Echoes).

Пылевые частицы оказывают существенное влияние на волновые процессы [5], протекающие в ионосфере и космическом пространстве. В настоящее время существенное развитие получили методы волновой диагностики пылевой ионосферной плазмы [6]. По существу, главные сведения о параметрах плазмы в ионосфере получены при помощи высокочастотных электромагнитных волн, применяемых для зондирования ионосферы с поверхности Земли или с борта космического аппарата.

Важным процессом, влияющим на распространение электромагнитных волн в плазме, является модуляционное взаимодействие [7]. В работе [8] проведено детальное исследование линейной стадии модуляционного взаимодействия – модуляционной неустойчивости электромагнитных волн при их распространении в запыленной ионосферной плазме во время высокоскоростных метеорных потоков таких, как Персеиды, Леониды, Геминиды и Ориониды. Оказалось [810], что модуляционное взаимодействие в запыленной ионосфере важно для объяснения таких эффектов, как:

1) наземные наблюдения низкочастотных ионосферных радиошумов с частотами ниже 50 Гц [11];

2) генерация инфразвуковых волн в ионосфере и возможность их детектирования у поверхности Земли;

3) усиление интенсивности зеленого свечения ночного неба на длине волны 557.7 нм от слоя нижней ионосферы, расположенного на высотах 110–120 км.

Кроме того [12], в условиях распространения в запыленной ионосфере достаточно интенсивных электромагнитных волн, например, при работе нагревных стендов типа HAARP, Tromsø и др. следует ожидать:

4) модуляционного возбуждения неоднородностей концентраций электронов и ионов в ионосфере.

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

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

Работа построена следующим образом. В разд. 2 вычисляются инкременты модуляционной неустойчивости. Влияние столкновений электронов, ионов и пылевых частиц с нейтралами на линейные свойства пылевых звуковых волн в ионосфере, а также возможные проявления модуляционного взаимодействия в запыленной ионосфере обсуждаются в разд. 3. В разд. 4 кратко сформулированы основные выводы исследования.

2. ИНКРЕМЕНТЫ МОДУЛЯЦИОННОЙ НЕУСТОЙЧИВОСТИ

В большинстве ситуаций детальное рассмотрение модуляционного взаимодействия оказывается возможным только на его линейной стадии, когда развивается модуляционная неустойчивость. Более того, аналитически сложно проанализировать случай неустойчивости волн с широким спектром. Этот случай подробно исследовался лишь для ленгмюровских [16] и нижнегибридных волн [17]. Данные исследования показали, что коротковолновая модуляционная неустойчивость волнового спектра (которая развивается при длинах волновых векторов модулирующих возмущений, намного больших, чем длина характерного волнового вектора в спектре) ведет себя, по существу, так же, как и коротковолновая модуляционная неустойчивость монохроматической волны. В этой связи, рассмотрим случай коротковолновой модуляционной неустойчивости электромагнитной волны в условиях запыленной ионосферы Земли.

Для описания модуляционной неустойчивости электромагнитных волн в запыленной ионосферной плазме используется стандартный подход [7]. В условиях диффузионного равновесия система основных уравнений имеет вид (ср. с [8, 9, 12])

(1)
$ - {{T}_{{e0}}}\nabla {{n}_{{e1}}} - {{n}_{{e0}}}\nabla {{T}_{{e1}}} + e{{n}_{{e0}}}{\mathbf{E}} - \frac{{{{n}_{{e0}}}{{e}^{2}}}}{{2{{m}_{e}}\omega _{0}^{2}}}\nabla {{\left| {\mathbf{E}} \right|}^{2}} = 0,$
(2)
$ - \kappa ({{T}_{{i0}}}\nabla {{n}_{{i1}}} + {{n}_{{i0}}}\nabla {{T}_{{i1}}}) - e{{n}_{{i0}}}{\mathbf{E}} = 0,$
(3)
$\frac{{{{\partial }^{2}}{{n}_{{d1}}}}}{{\partial {{t}^{2}}}} + {{\nu }_{{dn}}}\frac{{\partial {{n}_{{d1}}}}}{{\partial t}} = \frac{{{{n}_{{d0}}}{{q}_{{d0}}}\Delta \varphi }}{{{{m}_{d}}}},$
(4)
$\begin{gathered} \frac{3}{2}\frac{{\partial {{T}_{{e1}}}}}{{\partial t}} - \frac{{{{\chi }_{e}}}}{2}\Delta {{T}_{{e1}}} + {{{\bar {\nu }}}_{{ei}}}({{T}_{{i1}}} - {{T}_{{e1}}}) + \\ + \;{{{\bar {\nu }}}_{{en}}}{{T}_{{e1}}} - \frac{{{{T}_{{e0}}}}}{{{{n}_{{e0}}}}}\frac{{\partial {{n}_{{e1}}}}}{{\partial t}} = 0, \\ \end{gathered} $
(5)
$\begin{gathered} \frac{3}{2}\frac{{\partial {{T}_{{i1}}}}}{{\partial t}} - \frac{{{{\chi }_{i}}}}{2}\Delta {{T}_{{i1}}} + {{{\bar {\nu }}}_{{ei}}}({{T}_{{e1}}} - {{T}_{{i1}}}) + \\ + \;{{{\bar {\nu }}}_{{in}}}{{T}_{{i1}}} - \frac{{{{T}_{{i0}}}}}{{{{n}_{{i0}}}}}\frac{{\partial {{n}_{{i1}}}}}{{\partial t}} = \frac{{2{{\nu }_{e}}{{e}^{2}}}}{{{{m}_{e}}\omega _{0}^{2}}}{{\left| {\mathbf{E}} \right|}^{2}}, \\ \end{gathered} $
(6)
$\Delta \varphi = 4\pi \left( {{{n}_{{e1}}}e - {{n}_{{i1}}}e + {{q}_{{d0}}}{{n}_{{d1}}} + {{n}_{{d0}}}{{q}_{{d1}}}} \right),$
(7)
$\begin{gathered} \frac{{\partial {{q}_{{d1}}}}}{{\partial t}} = \mathop {\left( {\frac{{\partial ({{I}_{e}} + {{I}_{i}} + {{I}_{{ph}}})}}{{\partial {{q}_{d}}}}} \right)}\nolimits^{eq} {{q}_{{d1}}} + \mathop {\left( {\frac{{\partial {{I}_{e}}}}{{\partial {{n}_{e}}}}} \right)}\nolimits^{eq} {{n}_{{e1}}} + \\ + \;\mathop {\left( {\frac{{\partial {{I}_{e}}}}{{\partial {{T}_{e}}}}} \right)}\nolimits^{eq} {{T}_{{e1}}} + \mathop {\left( {\frac{{\partial {{I}_{i}}}}{{\partial {{n}_{i}}}}} \right)}\nolimits^{eq} {{n}_{{i1}}} + \mathop {\left( {\frac{{\partial {{I}_{i}}}}{{\partial {{T}_{i}}}}} \right)}\nolimits^{eq} {{T}_{{i1}}}. \\ \end{gathered} $

Здесь и далее все формулы и величины в СГС, ${{m}_{\xi }}$ – масса частиц сорта ξ; ${{n}_{\xi }}$ – концентрация частиц сорта ξ; $\xi = e,i,n,d$ для электронов, ионов, нейтралов и пылевых частиц соответственно; ${{T}_{{e(i)}}}$ – температура электронов (ионов), измеряемая в единицах энергии эрг, ${{T}_{n}}$ – температура нейтралов (поскольку температуру ионов и нейтралов на высотах 80–120 км можно считать одинаковой (см. например, [15], с. 15–16), в уравнениях (4), (5) не имеет смысла отдельно выделять температуру нейтралов, а следует использовать единое обозначение ${{T}_{i}}$); –e – заряд электрона, ионы предполагаются однозарядными, ${{q}_{d}}$, ${{m}_{d}}$, a – заряд, масса и радиус пылевой частицы соответственно; $\varphi $ – потенциал низкочастотного возмущения, ${\mathbf{E}}$ – электрическое поле электромагнитной волны накачки, ${{\omega }_{0}}$ – частота этой волны, $\kappa $ – показатель адиабаты (в случае изотермического процесса $\kappa = 1$, в случае адиабатического процесса $\kappa = 3$). Индекс “0” соответствует невозмущенным параметрам, индекс “1” – возмущенным величинам первого порядка малости, индекс “eq” обозначает, что соответствующая величина вычисляется при заряде пылевой частицы, равном его равновесному значению ${{q}_{{d0}}}$; ${{\bar {\nu }}_{{e(i)n}}} = {{\delta }_{{e(i)n}}}{{\nu }_{{e(i)n}}}$ – эффективная частота столкновений электронов (ионов) с нейтралами; ${{\nu }_{{\xi \eta }}}$ – частота столкновений частиц сорта $\xi = e,i,n,d$ с частицами сорта $\eta = e,i,n,d$, ${{\delta }_{{\xi \eta }}}$ – средняя доля энергии, передаваемая при однократном соударении частиц сорта $\xi = e,\;i,\;n,\;d$ с частицами сорта $\eta = e,\;i,\;n,\;d$; ${{\chi }_{e}} = 3.16{{T}_{e}}{\text{/}}({{m}_{e}}{{\nu }_{e}})$, ${{\chi }_{i}} = 3.9{{T}_{i}}{\text{/}}({{m}_{i}}{{\nu }_{i}})$ – электронный и ионный коэффициенты температуропроводности соответственно, ${{\nu }_{{e(i)}}} = {{\Sigma }_{{\xi = i(e),n,d}}}{{\nu }_{{e(i)\xi }}}$; $\Delta $ – оператор Лапласа; ${{I}_{{e(i)}}}$ – ток электронов (ионов) окружающей плазмы на пылевую частицу, ${{I}_{{ph}}}$ – ток фотоэлектронов (фототок). Фототок может возникать только в дневное время под действием солнечного излучения [18]. Однако фотоэффекту подвержены не все пылевые частицы, что обусловлено их составом. На высотах 80–120 км спектр солнечного излучения резко обрывается на длинах волн, порядка 170 нм, что соответствует энергии 7.3 эВ [19, 20]. Если, при этом, вещество, составляющее пылевую частицу, имеет работу выхода большую 7.3 эВ, то фотоэффект не происходит. Например, работа выхода чистого льда составляет около 8.9 эВ. Таким образом, пылевая частица, состоящая из чистого льда, не будет подвергаться фотоэффекту ни в дневное, ни, тем более, в ночное время.

Соударения электронов с ионами обычно носят упругий характер. В этом случае средняя доля энергии, теряемой электроном, составляет ${{\delta }_{{ie}}} = 2{{m}_{e}}{\text{/}}{{m}_{i}}$. Аналогичное значение имеет ${{\delta }_{{en}}}( = 2{{m}_{e}}{\text{/}}{{m}_{n}})$ и при упругих соударениях электронов с нейтралами. В ионосфере ${{\delta }_{{en}}}$, характеризующее упругие соударения, порядка (2–3) × 10–5. Упругие процессы учитывались при рассмотрении модуляционной неустойчивости [8]. Однако при взаимодействии электронов с нейтральными частицами важную роль играют и неупругие процессы, сопровождающиеся возбуждением ротационных, колебательных и оптических уровней. Поэтому для вычисления ${{\delta }_{{en}}}$ необходимо знать эффективные сечения всех неупругих процессов, существенных в ионосфере. Неупругие процессы могут оказывать существенное влияние на частоты столкновений. В настоящей работе используются частоты столкновений электронов и ионов с нейтралами для разных высот в ионосфере Земли из работы [15].

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

Основные выводы, касающиеся модуляционной неустойчивости электромагнитной волны в запыленной ионосфере, получаются при решении системы уравнений (1)–(7), которое осуществляется методами, описанными в работах [8, 12]. Предполагая, что низкочастотные возмущения в плазме во времени и пространстве изменяются как $exp( - i\Omega t + i{\mathbf{K}} \cdot {\mathbf{r}})$, где Ω и K – частота и волновой вектор, связанные с низкочастотными возмущениями, уравнения, описывающие эволюцию высокочастотного электромагнитного поля, можно записать в виде

(8)
${{\varepsilon }_{ \pm }}{{{\mathbf{E}}}_{ \pm }} - \frac{{{{c}^{2}}}}{{\omega _{ \pm }^{2}}}{{{\mathbf{k}}}_{ \pm }} \times ({{{\mathbf{k}}}_{ \pm }} \times {{{\mathbf{E}}}_{ \pm }}) = \frac{{{{n}_{{e1}}}}}{{{{n}_{{e0}}}}}\frac{{\omega _{{pe}}^{2}}}{{\omega _{ \pm }^{2}}}{{{\mathbf{E}}}_{{0 \pm }}},$
где ${{\omega }_{ \pm }} = \Omega \pm {{\omega }_{0}}$ и ${{{\mathbf{k}}}_{ \pm }} = {\mathbf{K}} \pm {{{\mathbf{k}}}_{0}}$; ${{{\mathbf{k}}}_{0}}$ – волновой вектор, соответствующий волне накачки; ${{\omega }_{{pe}}} = $ $ = \sqrt {4\pi {{n}_{{e0}}}{{e}^{2}}{\text{/}}{{m}_{e}}} $ – электронная плазменная частота; ${{\varepsilon }_{ \pm }} = 1 - \omega _{{pe}}^{2}{\text{/}}\omega _{ \pm }^{2}$ – высокочастотная диэлектрическая проницаемость ионосферной плазмы; ${{{\mathbf{E}}}_{ + }} = {\mathbf{E}}$, ${{{\mathbf{E}}}_{ - }} = {\mathbf{E}}{\text{*}}$, ${{{\mathbf{E}}}_{{0 + }}} = {{{\mathbf{E}}}_{0}}$, ${{{\mathbf{E}}}_{{0 - }}} = {\mathbf{E}}_{0}^{*}$. Здесь символ “*” означает комплексное сопряжение.

Основные выводы, касающиеся модуляционной неустойчивости электромагнитной волны в запыленной ионосфере, следующие. Возмущения на ионно-звуковых временных масштабах подавляются вследствие сильной амбиполярной диффузии [8, 10, 13]. Таким образом, модуляционная неустойчивость приводит к росту низкочастотных возмущений электрического поля, связанных с пылевой звуковой модой. Развитие модуляционной неустойчивости вызывается джоулевым нагревом, пондеромоторной силой, а также процессами, связанными с зарядкой пылевых частиц и их динамикой. При ${{q}_{{d0}}} < 0$ условие развития модуляционной неустойчивости следующее:

(9)
$\begin{gathered} \frac{{{{{\left| {{{{\mathbf{E}}}_{0}}} \right|}}^{2}}}}{{4\pi {{n}_{{e0}}}{{T}_{{e0}}}}} \gg \frac{3}{4}\frac{{{{\omega }_{\chi }} + {{{\bar {\nu }}}_{{en}}} + {{{\bar {\nu }}}_{{in}}}}}{{{{\nu }_{e}}}} \times \\ \times \;\frac{{\omega _{0}^{2}\nu _{e}^{2} + {{K}^{4}}{{c}^{4}}}}{{\omega _{{pe}}^{4}}}\frac{{\omega _{0}^{2}}}{{{{K}^{2}}{{c}^{2}}}}\frac{{\kappa + 4}}{{3\kappa }}. \\ \end{gathered} $

При ${{q}_{{d0}}} > 0$ модуляционная неустойчивость развивается, если

(10)
$\begin{gathered} \frac{{{{{\left| {{{{\mathbf{E}}}_{0}}} \right|}}^{2}}}}{{4\pi {{n}_{{e0}}}{{T}_{{e0}}}}} \gg {\text{max}}\left\{ {\frac{3}{8}\frac{{{{C}_{{Sd}}}K}}{{{{\nu }_{e}}}}\frac{{\omega _{0}^{2}\nu _{e}^{2} + {{K}^{4}}{{c}^{4}}}}{{{{K}^{2}}{{c}^{2}}}}\frac{{\omega _{0}^{2}}}{{\omega _{{pe}}^{4}}},} \right. \\ \left. {\frac{3}{8}\frac{{\mathop {\left( {{{\omega }_{{\chi e}}} + {{{\bar {\nu }}}_{{en}}}} \right)}\nolimits^3 }}{{{{\nu }_{e}}C_{{Sd}}^{2}{{K}^{2}}}}\frac{{\omega _{0}^{2}\nu _{e}^{2} + {{K}^{4}}{{c}^{4}}}}{{{{K}^{2}}{{c}^{2}}}}\frac{{\omega _{0}^{2}}}{{\omega _{{pe}}^{4}}}} \right\}. \\ \end{gathered} $

Здесь ${{{\mathbf{E}}}_{0}}$ – амплитуда электромагнитной волны накачки, ${{\omega }_{\chi }} = {{\omega }_{{\chi e}}} + {{\omega }_{{\chi i}}}$, ${{\omega }_{{\chi e}}} = {{\chi }_{e}}{{K}^{2}}{\text{/}}2$, ${{\omega }_{{\chi i}}} = {{\chi }_{i}}{{K}^{2}}{\text{/}}2$, $K = \left| {\mathbf{K}} \right|$ – длина волнового вектора модуляционных возмущений, c – скорость света, ${{C}_{{Sd}}} = $ $ = \left| {{{q}_{{d0}}}{\text{/}}e} \right|\sqrt {{{n}_{d}}{{T}_{e}}{\text{/}}{{n}_{e}}{{m}_{d}}} $ – скорость пылевого звука.

Приведем выражения для характерных зависимостей (от K) инкрементов модуляционной неустойчивости ${{\gamma }_{{mod}}}(K)$ и их максимальных значений Γ в наиболее важных с точки зрения описания процессов в ионосферной плазме случаях. При ${{q}_{{d0}}} < 0$ (например, в ночное время в отсутствие фотоэффекта), модуляционное возбуждение низкочастотных (пылевых звуковых) возмущений частоты Ω характеризуется следующими выражениями для инкрементов:

(11)
$\begin{gathered} {{\gamma }_{{mod}}}(K)\sim \\ \sim \mathop {\left( {\frac{{{{\nu }_{e}}{{c}^{2}}{{K}^{2}}}}{{{{\omega }_{\chi }}(K) + {{{\bar {\nu }}}_{{en}}} + {{{\bar {\nu }}}_{{in}}}}}} \right)}\nolimits^{1/2} \frac{{\omega _{{pe}}^{6}}}{{{{\nu }_{e}}\omega _{0}^{5}}}\mathop {\left( {\frac{{{{{\left| {{{{\mathbf{E}}}_{0}}} \right|}}^{2}}}}{{4\pi {{n}_{{e0}}}{{T}_{{e0}}}}}} \right)}\nolimits^{3/2} , \\ \end{gathered} $
(12)
$\Gamma \sim \mathop {\left( {\frac{{{{\nu }_{e}}{{c}^{2}}}}{{{{\chi }_{e}}}}} \right)}\nolimits^{1/2} \frac{{\omega _{{pe}}^{6}}}{{{{\nu }_{e}}\omega _{0}^{5}}}\mathop {\left( {\frac{{{{{\left| {{{{\mathbf{E}}}_{0}}} \right|}}^{2}}}}{{4\pi {{n}_{{e0}}}{{T}_{{e0}}}}}} \right)}\nolimits^{3/2} .$

При ${{q}_{{d0}}} > 0$ в случае ${{\omega }_{{\chi e}}} \gg \Omega \gg {{C}_{{Sd}}}K$ инкременты неустойчивости имеют вид

(13)
$\begin{gathered} {{\gamma }_{{mod}}}(K) \approx 2\sqrt 2 \mathop {\left( {{{\nu }_{e}}\frac{{C_{{Sd}}^{2}{{K}^{2}}}}{{{{\omega }_{{\chi e}}}(K) + {{{\bar {\nu }}}_{{en}}}}}} \right)}\nolimits^{1/2} \times \\ \times \;\frac{{\omega _{{pe}}^{2}}}{{{{\omega }_{0}}}}\frac{{Kc}}{{\sqrt {\omega _{0}^{2}\nu _{e}^{2} + {{K}^{4}}{{c}^{4}}} }}\mathop {\left( {\frac{{{{{\left| {{{{\mathbf{E}}}_{0}}} \right|}}^{2}}}}{{4\pi {{n}_{{e0}}}{{T}_{{e0}}}}}} \right)}\nolimits^{1/2} , \\ \end{gathered} $
(14)
$\Gamma \sim \mathop {\left( {{{\omega }_{0}}\frac{{C_{{Sd}}^{2}}}{{{{\chi }_{e}}}}} \right)}\nolimits^{1/2} \frac{{\omega _{{pe}}^{2}}}{{\omega _{0}^{2}}}\mathop {\left( {\frac{{{{{\left| {{{{\mathbf{E}}}_{0}}} \right|}}^{2}}}}{{4\pi {{n}_{{e0}}}{{T}_{{e0}}}}}} \right)}\nolimits^{1/2} .$

3. ЛИНЕЙНАЯ ДИСПЕРСИЯ И ОБСУЖДЕНИЕ РЕЗУЛЬТАТОВ

Из сравнения соотношений (12) и (14) с полученными ранее выражениями для максимальных инкрементов модуляционной неустойчивости [8] видно, что эти соотношения практически такие же как и те, которые были получены без учета неупругих столкновений электронов и ионов плазмы с нейтралами. Влияние таких столкновений на развитие модуляционной неустойчивости электромагнитных волн сказывается лишь для малых значений K (см. (11) и (13)). Так, например, на высоте 80 км это влияние является важным при $K < 7.5 \times {{10}^{{ - 3}}}$ см–1, на высоте 90 км – при K < $ < {{10}^{{ - 3}}}$ см–1, на высоте 100 км – при $K < 2.8 \times $ × 10‒4 см–1 днем и при $K < 2.6 \times {{10}^{{ - 4}}}$ см–1 ночью, на высоте 110 км – при $K < 6.1 \times {{10}^{{ - 5}}}$ см–1 днем и при $K < 5.6 \times {{10}^{{ - 5}}}$ см–1 ночью, на высоте 120 км – при $K < 2.8 \times {{10}^{{ - 5}}}$ см–1 днем и при $K < 2.6 \times $ × 10‒5 см–1 ночью. Характерным значением ${{K}_{{ch}}}$ длины волнового вектора возбуждаемых пылевых звуковых возмущений является ${{K}_{{ch}}} \sim {{\omega }_{d}}{\text{/}}{{C}_{{Sd}}}$, которое на высотах от 80 до 110 км изменяется в диапазоне от 0.03 до 3 см–1 (диапазон включает величины, полученные как для дневных, так и для ночных условий). Здесь ${{\omega }_{d}} = \sqrt {4\pi {{n}_{{d0}}}q_{{d0}}^{2}{\text{/}}{{m}_{d}}} $ – пылевая плазменная частота. Таким образом, для оценки влияния модуляционной неустойчивости на процессы, происходящие в запыленной ионосфере Земли, можно использовать полученные ранее [8] выражения для максимальных инкрементов модуляционной неустойчивости, выведенные без учета неупругих столкновений электронов и ионов плазмы с нейтралами. Однако при этом инкремент модуляционной неустойчивости должен превосходить декремент затухания пылевой звуковой волны. Иначе, развитие модуляционной неустойчивости происходить не может.

При $\omega \ll {{\omega }_{\chi }}$ линейное дисперсионное соотношение пылевых звуковых волн определяется выражением

(15)
$\omega \left( k \right) = {\text{Re}}\omega \left( k \right) + i{\text{Im}}\omega \left( k \right),$
где

(16)
${\text{Re}}\omega \left( k \right) = \sqrt {\frac{{C_{{Sd}}^{2}{{k}^{2}}}}{{1 + {{k}^{2}}\lambda _{D}^{2}}} - \frac{{\nu _{{dn}}^{2}}}{4}} ,$
(17)
${\text{Im}}\omega \left( k \right) = - \frac{{{{\nu }_{{dn}}}}}{2} - \frac{{\mathop {\left( {{\text{Re}}\omega \left( k \right)} \right)}\nolimits^2 }}{2}\left( {\frac{{{{\nu }_{{en}}}}}{{{{k}^{2}}v_{{Te}}^{2}}} + \frac{{{{\nu }_{{in}}}}}{{{{k}^{2}}v_{{Ti}}^{2}}}} \right).$

Здесь ${{\nu }_{{dn}}} = \left( {4{\text{/}}3} \right)\pi {{a}^{2}}\sqrt {8{{T}_{{n0}}}{\text{/}}\pi {{m}_{n}}} {{n}_{n}}\left( {{{m}_{n}}{\text{/}}{{m}_{d}}} \right)$ – эффективная частота столкновений пылевых частиц и нейтралов, $\lambda _{D}^{{ - 2}} = \lambda _{{De}}^{{ - 2}} + \lambda _{{Di}}^{{ - 2}}$, ${{\lambda }_{{De(i)}}} = \sqrt {{{T}_{{e(i)}}}{\text{/}}4\pi {{n}_{{e(i)}}}{{e}^{2}}} $ – дебаевский радиус электронов (ионов), ${{v}_{{Te(i,d)}}}$ – тепловая скорость электронов (ионов, пылевых частиц). Выражение (17), полученное в предположениях

(18)
$\begin{gathered} \omega \gg {{\nu }_{{dn}}},k{{v}_{{Td}}};\quad {{\nu }_{{en}}} \gg \omega ,\;k{{v}_{{Te}}};\quad \omega {{\nu }_{{en}}} \ll {{k}^{2}}v_{{Te}}^{2}; \\ {{\nu }_{{in}}} \gg \omega ,k{{v}_{{Ti}}};\quad \omega {{\nu }_{{in}}} \ll {{k}^{2}}v_{{Ti}}^{2} \\ \end{gathered} $
(ср. с (4.3.20)–(4.3.21) из [21]), учитывает столкновения с нейтралами пылевых частиц, электронов и ионов. Условия (18) обычно выполнены для пылевых звуковых волн в ионосфере Земли. Далее, слагаемые (17), характеризующие вклад столкновений электронов и ионов с нейтралами являются функциями волнового вектора k. Вклад этих слагаемых в условиях запыленной ионосферной плазмы обычно не превосходит вклада от ${{\nu }_{{dn}}}{\text{/}}2$. Таким образом, условие развития модуляционной неустойчивости имеет вид: $\Gamma > {{\nu }_{{dn}}}{\text{/}}2$. Пороговое значение амплитуды поля электромагнитной волны $\left| {{{{\mathbf{E}}}_{0}}} \right|$, характеризующее развитие модуляционной неустойчивости, определяется из условия
$\Gamma = {{\nu }_{{dn}}}{\text{/}}2,$
где $\Gamma $ определяется соотношениями (12) и (14) в случае соответственно отрицательных и положительных зарядов пылевых частиц. С другой стороны, в выражении для действительной части частоты (16) в подкоренном выражении присутствует слагаемое $\nu _{{dn}}^{2}{\text{/}}4$. Таким образом, для существования распространяющейся пылевой звуковой волны необходимо, чтобы ${{\omega }_{d}} \approx {{C}_{{Sd}}}{\text{/}}{{\lambda }_{D}} > {{\nu }_{{dn}}}{\text{/}}2$.

В табл. 1 представлены величины, характеризующие развитие модуляционной неустойчивости электромагнитных волн, а также распространение и затухание пылевых звуковых возмущений на разных высотах над Землей, при ${{n}_{d}} = 100$ см–3, $a = 0.1$ мкм, ${{\omega }_{0}} = 2\pi \times {{10}^{6}}$ с–1, а также параметров ионосферы на высотах 80–120 км ночью и днем [15].

Таблица 1.

Величины, характеризующие развитие модуляционной неустойчивости электромагнитных волн, а также распространение и затухание пылевых звуковых возмущений на разных высотах

  Ночь День
высота, км $\Gamma \mathop {\left( {\frac{{{\text{|}}{{{\mathbf{E}}}_{0}}{{{\text{|}}}^{2}}}}{{4\pi {{n}_{{e0}}}{{T}_{{e0}}}}}} \right)}\nolimits^{ - 3/2} $, с–1 $\frac{{{{C}_{{Sd}}}}}{{{{\lambda }_{D}}}}$, с–1 $\frac{{{{\nu }_{{dn}}}}}{2}$, с–1 $\Gamma \mathop {\left( {\frac{{{\text{|}}{{{\mathbf{E}}}_{0}}{{{\text{|}}}^{2}}}}{{4\pi {{n}_{{e0}}}{{T}_{{e0}}}}}} \right)}\nolimits^{ - 1/2} $, с–1 $\frac{{{{C}_{{Sd}}}}}{{{{\lambda }_{D}}}}$, с–1 $\frac{{{{\nu }_{{dn}}}}}{2}$, с–1
80 10.6 0.38 12.9 2.8 60.2 12.9
90 2.2 × 103 0.74 1.8 3.1 61.0 1.8
100 1.7 × 107 1.45 0.5 5.4 62.3 0.5
110 5.1 × 107 1.8 0.1 3.4 63.0 0.1
120 7.0 × 107 2.5 0.05 2.5 61.9 0.05

Из табл. 1 видно, что для используемых параметров развитие модуляционной неустойчивости затруднено на высотах 80 и 90 км днем и 80 км ночью, поскольку нарушается соотношение $\Gamma \geqslant {{\nu }_{{dn}}}{\text{/}}2$ и, кроме того, выражения для инкрементов модуляционной неустойчивости справедливы в условиях слабой нелинейности ${{\left| {{{{\mathbf{E}}}_{0}}} \right|}^{2}}{\text{/}}4\pi {{n}_{{e0}}}{{T}_{{e0}}} \ll 1$. Далее, ночью на высотах 80 и 90 км пылевые звуковые волны не могут распространяться, так как нарушается неравенство ${{C}_{{Sd}}}{\text{/}}{{\lambda }_{D}} > {{\nu }_{{dn}}}{\text{/}}2$. Следовательно, на высотах 80 и 90 км невозможно ожидать развития пылевой звуковой турбулентности (днем – из-за невозможности развития модуляционной неустойчивости, за счет которой пылевая звуковая турбулентность могла бы возбуждаться, ночью – из-за невозможности распространения пылевых звуковых волн).

Таким образом, влияние модуляционной неустойчивости на распространение электромагнитных волн в плазме запыленной ионосферы наиболее существенно на высотах 100–120 км. На меньших высотах оно может быть подавлено. Соответственно, модуляционное взаимодействие в запыленной ионосфере важно для объяснения таких эффектов, как наземные наблюдения низкочастотных ионосферных радиошумов с частотами ниже 50 Гц [11], поскольку генерация низкочастотных радиошумов за счет модуляционного взаимодействия вполне может происходить на высотах 100–120 км. Кроме того, при этом пороговые значения амплитуды поля электромагнитной волны $\left| {{{{\mathbf{E}}}_{0}}} \right|$, характеризующие генерацию наблюдаемых низкочастотных радиошумов, определяются из условия (19), что соответствует данным наблюдений [8]. Возможно, отсутствие наблюдений низкочастотных ионосферных радиошумов с частотами ниже 50 Гц во время таких явлений, как серебристые облака и полярные мезосферные радиоотражения связано с тем, что генерация радиошумов за счет модуляционных процессов возможна на высотах 100–120 км, а серебристые облака и полярные мезосферные радиоотражения обусловлены присутствием пылевой плазмы на высотах 80–95 км. Соответственно, корреляция наблюдений низкочастотных ионосферных радиошумов [11] с метеорными потоками Персеиды, Леониды, Геминиды, Ориониды обусловлена возможностью генерации пылевой плазмы на высотах 100–120 км вследствие взаимодействия этих потоков с ионосферой.

Далее, тот факт, что влияние модуляционной неустойчивости в плазме запыленной ионосферы существенно на высотах 100–120 км, никоим образом не сказывается на выводах о том, что этот эффект приводит к генерации инфразвуковых волн в ионосфере (только с учетом эффекта на высотах 100–120 км, а не 80–120 км), детектирование которых оказывается возможным у поверхности Земли [9], а также на объяснении эффекта усиления интенсивности зеленого свечения ночного неба на длине волны 557.7 нм, поскольку этот эффект происходит в слое нижней ионосферы, расположенном на высотах 110–120 км. Что же касается модуляционного возбуждения неоднородностей концентраций электронов и ионов в ионосфере [12], то данный эффект, например, при работе нагревных стендов типа HAARP, Tromsø и др. также может происходить на высотах 100–120 км.

4. ЗАКЛЮЧЕНИЕ

Итак, проведена оценка роли неупругих столковений электронов и ионов с нейтралами при развитии модуляционной неустойчивости электромагнитных волн с участием пылевых звуковых возмущений в плазме запыленной ионосферы, а также влияния столкновений электронов, ионов и пылевых частиц с нейтралами на проявления модуляционного взаимодействия в запыленной ионосфере. Оказывается, что влияние столковений электронов и ионов с нейтралами при этом обычно менее существенно, чем влияние столкновений пылевых частиц с нейтралами. В частности, именно столкновения пылевых частиц с нейтралами определяют условия развития модуляционной неустойчивости электромагнитных волн и распространения пылевых звуковых волн в запыленной ионосфере. Влияние модуляционной неустойчивости на распространение электромагнитных волн в плазме запыленной ионосферы наиболее существенно на высотах 100–120 км. Корреляция наблюдений низкочастотных ионосферных радиошумов с частотами ниже 50 Гц [11] с метеорными потоками Персеиды, Леониды, Геминиды, Ориониды обусловлена возможностью генерации пылевой плазмы на высотах 100–120 км вследствие взаимодействия этих потоков с ионосферой и развития в ней модуляционного взаимодействия. Отсутствие наблюдений низкочастотных ионосферных радиошумов во время таких явлений, как серебристые облака и полярные мезосферные радиоотражения, обусловленных присутствием пылевой плазмы на высотах 80–95 км, связано с тем, что модуляционные процессы на этих высотах подавлены. Модуляционное взаимодействие электромагнитных волн в плазме запыленной ионосферы приводит к генерации инфразвуковых волн на высотах 100–120 км, детектирование которых оказывается возможным у поверхности Земли. Кроме того, оно позволяет объяснить эффект усиления интенсивности зеленого свечения ночного неба на длине волны 557.7 нм, происходящий в слое нижней ионосферы, расположенном на высотах 110–120 км. Модуляционное возбуждение неоднородностей концентраций электронов и ионов в ионосфере [12] при работе нагревных стендов типа HAARP, Tromsø и др. также может происходить на высотах 100–120 км.

Список литературы

  1. Shukla P.K., Mamun A.A. Introduction to Dusty Plasmas Physics. Bristol/Philadelphia: Instit. Phys. Publ., 2002.

  2. Tsytovich V.N., Morfill G.E., Vladimirov S.V., Thomas H. Elementary Physics of Complex Plasmas. Berlin/Heidelberg: Springer, 2008.

  3. Fortov V.E., Ivlev A.V., Khrapak S.A., Khrapak A.G., Morfill G.E. // Phys. Reports. 2005. V. 421. P. 1.

  4. Turco R.P., Toon O.B., Whitten R.C., Keesee R.G., Hollenbach D. // Planet. Space Sci. 1982. V. 30. P. 1147.

  5. Verheest F. // Waves in Dusty Space Plasmas. Dordrecht: Kluwer Acad. Publ., 2000.

  6. Buonsanto M.J., Foster J.C., Sipler D.P. // J. Geophys. Res. 1992. V. 97. P. 1225.

  7. Vladimirov S.V., Tsytovich V.N., Popel S.I., Khakimov F.Kh. // Modulational Interactions in Plasmas. Dordrecht: Kluwer Acad. Publ., 1995.

  8. Копнин С.И., Попель С.И., Ю Минг // Физика плазмы. 2007. Т. 33. С. 323.

  9. Копнин С.И., Попель С.И. // Физика плазмы. 2008. Т. 34. С. 517.

  10. Kopnin S.I., Popel S.I., Yu M.Y. // Phys. Plasmas. 2009. V. 16. P. 063705.

  11. Мусатенко C.И., Мусатенко Ю.С., Курочка Е.В., Ласточкин А.В., Чолий B.Я., Максименко О.И., Слипченко А.С. // Геомагнетизм и аэрономия. 2006. Т. 46. С. 182.

  12. Копнин С.И., Попель С.И., Mорозова Т.И. // Физика плазмы. 2015. Т. 41. С. 188.

  13. Lie-Svendsen Ш., Blix T.A., Hoppe U.-P., Thrane E.V. // J. Geophys. Res. 2003. V. 108. P. 8442.

  14. Popel S.I., Kopnin S.I., Kosarev I.N., Yu M.Y. // Adv. Space Res. 2006. V. 37. P. 414.

  15. Гуревич А.В., Шварцбург А.Б. Нелинейная теория распространения радиоволн в ионосфере. М.: -Наука, 1973.

  16. Popel S.I., Tsytovich V.N., Vladimirov S.V. // Phys. Plasmas. 1994. V. 1. P. 2176.

  17. Попель С.И. // Физика плазмы. 1998. Т. 24. С. 1093.

  18. Копнин С.И., Моржакова А.А., Попель С.И., Шук-ла П.К. // Физика плазмы. 2011. Т. 37. С. 745.

  19. Клумов Б.А., Морфилл Г.Е., Попель С.И. // ЖЭТФ. 2005. Т. 127. С. 171.

  20. Дубинский А.Ю., Попель С.И. // Письма в ЖЭТФ. 2012. Т. 96. С. 22.

  21. Александров А.Ф., Рухадзе А.А. Лекции по электродинамике плазмоподобных сред. М.: Изд-во Моск. ун-та. Физ. факультет, 1999.

Дополнительные материалы отсутствуют.