Как зарождается молния и почему она ветвится: физика диэлектрического пробоя и атмосферного электричества

Каждую секунду в атмосфере Земли вспыхивает от 40 до 50 молний — почти 4 миллиона высокоэнергетических разрядов в сутки. Спустя более 270 лет после знаменитого эксперимента Бенджамина Франклина с воздушным змеем (1752 г.) один из фундаментальных вопросов физики атмосферы долгое время оставался нерешённым: как именно зарождается молния внутри грозового облака? Метеозонды и самолёты-лаборатории внутри активных кучево-дождевых облаков регистрируют максимальные макроскопические электрические поля лишь в диапазоне \(0,1\text{–}0,4\,\text{МВ/м}\) — на порядок ниже классического порога диэлектрического пробоя сухого воздуха на уровне моря (\(E_{b0} \approx 3,0\,\text{МВ/м}\)), определяемого теорией лавин Таунсенда и законом Пашена. Почему молния никогда не бьёт по прямой, а образует сложное фрактальное дерево, подчиняющееся лапласовской неустойчивости роста? Почему отрицательный ступенчатый лидер (stepped leader) продвигается дискретными скачками по \(50\,\text{м}\) каждые \(50\,\text{мкс}\), а затем вызывает ослепительный обратный удар (return stroke), несущийся от земли вверх со скоростью в треть скорости света (\(v \approx c/3\)) и разогревающий канал до \(30\,000\,\text{К}\)? В этой статье мы последовательно выведем полную физическую картину молнии.

1. Электростатический генератор кучево-дождевого облака: безындукционная электризация и трипольная структура

Перед пробоем зрелое кучево-дождевое облако (Cumulonimbus) должно создать разность потенциалов \(10^8\text{–}10^9\,\text{В}\). В турбулентной смешанной фазовой зоне на высотах \(4\text{–}10\,\text{км}\) (температура от \(0^\circ\text{C}\) до \(-40^\circ\text{C}\)) мощные восходящие потоки (\(v_{\text{up}} \approx 10\text{–}30\,\text{м/с}\)) поддерживают сосуществование трёх видов гидрометеоров: переохлаждённых капель воды (\(d \approx 10\text{–}20\,\text{мкм}\)), растущих из пара ледяных кристаллов (\(20\text{–}200\,\text{мкм}\)) и миллиметровых частиц снежной крупы (graupel, \(1\text{–}5\,\text{мм}\)), быстро растущих за счёт обзернения (riming) — намерзания переохлаждённых капель при ударе.

Эксперименты Такахаси (Takahashi, 1978) и Сондерса (Saunders, 1991) доказали, что основной механизм — безындукционное разделение зарядов при упругих столкновениях частиц крупы с кристаллами льда в присутствии переохлаждённых капель (перенос \(10\text{–}100\,\text{фКл} = 10^{-14}\text{–}10^{-13}\,\text{Кл}\) за соударение). Микроскопическая природа эффекта связана с нанометровым квазижидким слоем (Quasi-Liquid Layer, QLL) на поверхности льда. Молекулы воды в решётке и QLL диссоциируют на протоны (\(\text{H}^+\)) и гидроксил-ионы (\(\text{OH}^-\)). Благодаря эстафетному механизму Гротгуса вдоль водородных связей подвижность протонов \(\text{H}^+\) на несколько порядков превышает подвижность \(\text{OH}^-\).

При контакте конкурируют два термодинамических градиента: (1) эффект температурного градиента: частица крупы нагревается скрытой теплотой замерзания переохлаждённых капель, поэтому из более тёплой крупы с высокой концентрацией ионов протоны \(\text{H}^+\) быстро диффундируют в холодный кристалл, заряжая крупу отрицательно (\(-\)), а кристалл положительно (\(+\)); (2) эффект толщины QLL (механизм Бейкера — Дэша): при более высокой температуре и высокой водности (\(\text{LWC}\)) более толстый слой QLL на крупе меняет знак переноса заряда. Граница между режимами определяет температуру инверсии (\(T_R \approx -15^\circ\text{C} \dots -20^\circ\text{C}\)):

Холодный режим (\(T < T_R\), высота \(7\text{–}10\,\text{км}\)): тяжёлая крупа заряжается отрицательно (\(-\)), лёгкие кристаллы льда — положительно (\(+\)).
Тёплый режим (\(T > T_R\), высота \(4\text{–}6\,\text{км}\)): инверсия знака — тяжёлая крупа заряжается положительно (\(+\)), лёгкие кристаллы — отрицательно (\(-\)).

Затем гравитационная сортировка по размерам разделяет заряды в масштабе всего облака: тяжёлая крупа (\(v_t \approx 5\text{–}10\,\text{м/с}\)) остаётся в средних и нижних слоях, а лёгкие кристаллы (\(v_t \ll 1\,\text{м/с}\)) уносятся восходящим потоком в верхнюю наковальню, формируя классическую электрическую трипольную структуру:

  flowchart TD
    subgraph Cloud["Электрическая трипольная структура кучево-дождевого облака"]
        UP["<strong>Верхняя основная положительная область</strong><br/>Высота 8–12 км (-30°C ... -50°C)<br/>Положительные кристаллы льда в наковальне (+40 ... +100 Кл)"]
        MN["<strong>Средняя основная отрицательная область</strong><br/>Высота 5–8 km (-10°C ... -25°C, T &lt; T_R)<br/>Отрицательная крупа и кристаллы тёплой зоны (-40 ... -100 Кл)"]
        LP["<strong>Нижний положительный центр заряда (LPCC)</strong><br/>Высота 3–5 км (0°C ... -10°C, T &gt; T_R)<br/>Положительная крупа ниже уровня инверсии (+5 ... +20 Кл)"]
    end
    Col1["Столкновения в холодной зоне (T &lt; -15°C)<br/>Крупа (-) и Кристалл (+)"] -->|"Восходящий поток уносит кристаллы"| UP
    Col1 -->|"Гравитация удерживает крупу"| MN
    Col2["Столкновения в тёплой зоне (T &gt; -15°C)<br/>Крупа (+) и Кристалл (-)"] -->|"Гравитация опускает крупу"| LP
    Col2 -->|"Восходящий поток уносит кристаллы"| MN
    MN -->|"Электростатическая индукция"| Ground["<strong>Проводящая поверхность Земли (индуцированный + заряд)</strong>"]

Плотность пространственного заряда \(\rho(\mathbf{r})\) определяет потенциал \(\phi(\mathbf{r})\) и электрическое поле \(\mathbf{E} = -\nabla\phi\) через уравнение Пуассона:

\[\nabla^2 \\phi(\mathbf{r}) = -\frac{\rho(\mathbf{r})}{\varepsilon_0}\]

Нижний положительный центр заряда (LPCC) усиливает направленное вниз поле у основания отрицательного слоя, работая как электростатическая катапульта для запуска отрицательных лидеров к земле. В зимних грозах (например, на побережье Хокурику в Японии) изотерма \(-15^\circ\text{C}\) опускается ниже \(2\,\text{км}\), и верхняя положительная область разряжается напрямую в землю в виде гигантских положительных суперболтов (+CG) с токами до \(200\text{–}300\,\text{кА}\).

2. Загадка инициации молнии: пробой Таунсенда против релятивистских лавин (RREA)

В классической физике газового разряда ударная ионизация электронами (первый коэффициент Таунсенда \(\alpha\)) конкурирует с прилипанием электронов к электроотрицательным молекулам кислорода (коэффициент \(\eta_a\)):

\[\frac{dN_e}{dx} = (\alpha - \eta_a) N_e \implies N_e(x) = N_0 \exp\left[ (\alpha - \eta_a) x \right]\]

Для экспоненциальной лавины Таунсенда необходимо \(\alpha > \eta_a\), что на уровне моря даёт порог \(E_{b0} \approx 3,0\,\text{МВ/м}\) (\(E/N \approx 120\,\text{Тд}\)). Согласно закону Пашена, пробивное поле убывает пропорционально относительной плотности воздуха \(n(z) = \exp(-z/H)\):

\[E_b(z) = E_{b0} \cdot n(z) \approx (3.0\,\\text{MV/m}) \cdot \exp(-z/H)\]

На высоте зарождения молний \(z \approx 6\,\text{км}\) (\(n \approx 0,5\)) классический пробой требует \(E_b(6\,\text{км}) \approx 1,5\,\text{МВ/м}\), тогда как измерения дают не более \(0,1\text{–}0,4\,\text{МВ/м}\). Этот парадокс разрешают два механизма:

Лавины релятивистских убегающих электронов (RREA), инициируемые космическими лучами: тепловые электроны останавливаются гигантским пиком неупругого торможения при \(100\text{–}200\,\text{эВ}\) (\(F_{\max} \approx 27\,\text{МэВ/м}\cdot n(z)\)). Однако по релятивистской формуле Бете — Блоха при кинетической энергии \(K_{\min} \approx 1,0\text{–}1,4\,\text{МэВ}\) сила трения падает до минимума \(F_{\min} \approx 0,218\,\text{МэВ/м}\cdot n(z)\). С учётом кулоновского рассеяния на ядрах пороговое поле RREA \(E_{\text{RREA}}\) (Гуревич, 1992; Дуайер, 2003) составляет:

\[E_{\\text{RREA}}(z) \approx 0.284\,\\text{MV/m} \cdot n(z), \qquad N_{\\text{RREA}}(L) = N_0 \exp\left(\frac{L}{\lambda_{\\text{RREA}}}\right)\]

На высоте \(6\,\text{км}\) порог равен \(E_{\text{RREA}} \approx 0,14\,\text{МВ/м}\), что в точности совпадает с данными зондов! Затравочные электроны (\(\sim 1\,\text{МэВ}\)) из широких атмосферных ливней (ШАЛ) от космических лучей переходят в режим убегания и размножаются за счёт мёллеровского рассеяния. 2. Усиление поля на гидрометеорах и быстрый положительный пробой (Fast Positive Breakdown): на остриях вытянутых ледяных кристаллов (\(l/r \gg 1\)) диэлектрическая поляризация усиливает локальное поле:

\[E_{\\text{tip}} \approx E_0 \frac{(l/r)^2}{\ln(2l/r) - 1}\]

Когда \(E_{\text{tip}}\) превышает \(E_b(z)\), стартуют холодные положительные коронные стримеры, способные неограниченно распространяться при фоновом поле выше порога стабильности \(E_{\text{st}}^{+} \approx 0,45\,\text{МВ/м}\cdot n(z)\). В 2016 г. Рисон и др. с помощью широкополосной УКВ-интерферометрии показали, что молнии начинаются со взрывного каскада положительных стримеров (\(v \approx 3\text{–}5 \times 10^7\,\text{м/с}\)) — быстрого положительного пробоя (FPB / узких биполярных событий NBE).

3. Почему молния ветвится: переход стример–лидер и \(\eta\)-лапласовский фрактальный рост

Одиночный стример — это холодная плазма (\(T_g \approx 300\,\text{К}\)), теряющая проводимость за десятки наносекунд из-за прилипания электронов к \(\text{O}_2\). Когда тысячи стримеров сливаются в общий стебель (stem), джоулев нагрев поднимает температуру газа выше критического порога \(T_{\text{crit}} \approx 1500\text{–}5000\,\text{К}\), запуская: (1) термический отлип электронов (\(\text{O}_2^- + \text{M} \to \text{O}_2 + \text{M} + \text{e}^-\)), (2) колебательно-поступательную (V-T) релаксацию \(\text{N}_2\) и (3) гидродинамическое расширение канала, снижающее плотность \(N\) и повышающее \(E/N\). Так рождается горячий проводящий лидер (leader), переносящий потенциал облака (\(-50\,\text{МВ}\)) к своей головке.

В 90% разрядов облако–земля (отрицательных, -CG) ступенчатый лидер (stepped leader) движется скачками, поскольку отрицательным стримерам нужно вдвое более сильное поле, чем положительным (\(E_{\text{st}}^{-} \approx 1,0\text{–}1,2\,\text{МВ/м}\cdot n\)). Впереди головки лидера зажигается изолированное плазменное ядро (объёмный стебель / space stem), из которого растёт двунаправленный лидер; когда его положительный конец смыкается с основным каналом, вспыхивает новая ступень длиной \(L_{\text{step}} \approx 10\text{–}50\,\text{м}\) каждые \(\tau_{\text{step}} \approx 20\text{–}50\,\text{мкс}\) (средняя скорость \(\sim 2 \times 10^5\,\text{м/с}\)).

Фрактальное ветвление объясняет модель диэлектрического пробоя (DBM) Нимейера, Пьетронеро и Висманна (1984). Проводящее дерево лидера (\(\phi = 0\)) находится в лапласовском поле, а вероятность пробоя \(P_i\) в граничной точке \(i\) пропорциональна \(\eta\)-й степени локального поля:

\[\nabla^2 \\phi(\mathbf{r}) = 0, \qquad P_i = \frac{|\nabla \\phi_i|^\eta}{\sum_{j \in \\text{boundary}} |\nabla \\phi_j|^\eta} \quad \left(|\nabla \\phi_i| > E_{\\text{crit}}\right)\]

Выступающие вершины сгущают эквипотенциали, а боковые ветви экранируются (неустойчивость Маллинза — Секерки), вызывая расщепление вершин. При \(\eta = 1\) модель изоморфна диффузионно-ограниченной агрегации (DLA, 2D \(D_f \approx 1,71\)), а нелинейный режим молнии (\(\eta \approx 2\)) даёт 2D-проекцию фрактальной размерности Хаусдорфа \(D_f \approx 1,35\text{–}1,55\).

4. Соединение, обратный удар при \(30\,000\,\text{К}\) и ударно-волновая механика грома

Когда ступенчатый лидер приближается к земле на дистанцию ориентировки (striking distance) \(r_s\), с вершин заземлённых объектов (молниеотводов, деревьев) стартуют встречные положительные лидеры. В электрогеометрической модели (EGM, метод катящейся сферы) связь \(r_s\) с пиковым током \(I_p\) (в \(\text{кА}\)) описывается формулой Лава (1973):

\[r_s = 10 \, I_p^{0.65} \quad [\\text{m}]\]

В момент контакта (\(20\text{–}100\,\text{м}\) над землёй) волна нейтрализации потенциала — обратный удар (return stroke) — устремляется от земли к облаку со скоростью \(v_{\text{RS}} \approx c/3 \approx 1,0 \times 10^8\,\text{м/с}\). За \(1\text{–}5\,\mu\text{s}\) пиковый ток \(\sim 30\,\text{кА}\) проходит через ядро диаметром \(1\text{–}3\,\text{см}\), разогревая плазму до \(T \approx 30\,000\,\text{К}\) (в пять раз горячее фотосферы Солнца). Последующие стреловидные лидеры (dart leaders) пробегают по остывающему главному каналу со скоростью \(\sim 10^7\,\text{м/с}\), вызывая мерцание молнии.

Поскольку микросекундный нагрев является строго изохорным, давление в канале скачком возрастает до \(10\text{–}100\,\text{атм}\), и плазменный цилиндр расширяется со сверхзвуковой скоростью \(10\text{–}20\,\text{Махов}\), порождая цилиндрическую взрывную волну Тейлора — Седова. Характерный радиус релаксации \(R_0\) и пиковая частота грома \(f_{\text{peak}}\) определяются погонной энергией \(E_l\) и атмосферным давлением \(P_0\):

\[R_0 = \sqrt{\frac{E_l}{\pi P_0}}, \qquad f_{\\text{peak}} \approx 0.63 \, c_s \sqrt{\frac{P_0}{E_l}} = \frac{0.63 \, c_s}{\sqrt{\pi} R_0} \approx 50\\text{–}150\,\\text{Hz}\]

На расстоянии нескольких километров N-волна превращается в долгие раскаты грома из-за разницы времён прихода звука от извилистого 3D-канала, квадратичного по частоте затухания в воздухе (\(\alpha_{\text{att}} \propto f^2\)) и тепловой рефракции звука вверх.

5. Транзиентные световые явления (TLE) и земные гамма-вспышки (TGF) с рождением антиматерии

Транзиентные световые явления (TLE):
- Красные спрайты (высота \(50\text{–}90\,\text{км}\)): после мощного удара \(+CG\) с большим изменением дипольного момента заряда (\(\text{CMC} \gtrsim 500\,\text{Кл}\cdot\text{км}\)) квазиэлектростатическое поле в мезосфере превышает экспоненциально падающий с плотностью воздуха порог пробоя \(E_b(z)\), зажигая гигантские комплексы холодных стримеров (красное свечение \(\text{N}_2(1\text{P})\) сверху, синие усики \(\text{N}_2^+(1\text{N})\) снизу).
- Синие джеты (\(15\text{–}50\,\text{км}\)) и Гигантские джеты (до \(90\,\text{км}\)): конусы положительных лидеров-стримеров, бьющие вверх из вершины облака со скоростью \(\sim 100\,\text{км/с}\).
- Эльфы (ELVES, \(85\text{–}100\,\text{км}\)): расширяющиеся со сверхсветовой фазовой скоростью красные кольца диаметром \(300\text{–}500\,\text{км}\) (\(<1\,\text{мс}\)) при пересечении сферического ЭМИ (EMP) обратного удара с нижней ионосферой.
Земные гамма-вспышки (TGF) и фотоядерное рождение антиматерии: лавины RREA ускоряют электроны до \(20\text{–}40\,\text{МэВ}\), излучая тормозные гамма-кванты. В 2017 г. Эното и др. (Enoto et al., Nature) доказали, что гамма-фотоны с энергией \(>10,55\,\text{МэВ}\) выбивают нейтроны из ядер \(^{14}\text{N}\); нестабильный изотоп \(^{13}\text{N}\) испытывает \(\beta^+\)-распад с испусканием позитронов антиматерии (\(\text{e}^+\)), аннигиляция которых с электронами воздуха даёт характерную гамма-линию \(0,511\,\text{МэВ}\):

\[\gamma + {{}}^{14}_{7}\\text{N} \longrightarrow {{}}^{13}_{7}\\text{N} + \\text{n} \quad (E_\gamma > 10.55\,\\text{MeV})\]\[{{}}^{13}_{7}\\text{N} \longrightarrow {{}}^{13}_{6}\\text{C} + \\text{e}^+ + \nu_e, \qquad \\text{e}^+ + \\text{e}^- \longrightarrow 2\gamma \quad (E_\gamma = 0.511\,\\text{MeV})\]

6. Сводное физическое сравнение и современные методы наблюдений

Явление / Режим	Высота / Пространственный масштаб	Длительность	Скорость распространения \(v\)	Характерный ток, энергия или температура	Основной физический механизм
Безындукционная зарядка	\(4\text{–}10\,\text{км}\) (смешанная зона)	\(10\text{–}30\,\text{мин}\) цикл конвекции	Восходящий поток \(10\text{–}30\,\text{м/с}\)	\(\sim 10^{-14}\,\text{Кл}\)/удар; \(\pm 40\text{–}100\,\text{Кл}\) в облаке	Соударения крупы и кристаллов, диффузия \(\text{H}^+\) в QLL, \(T_R \approx -15^\circ\text{C}\)
Быстрый положительный пробой (FPB)	\(5\text{–}12\,\text{км}\); длина \(0,3\text{–}1\,\text{км}\)	\(10\text{–}20\,\mu\text{s}\)	\(3 \times 10^7 \text{–} 5 \times 10^7\,\text{м/с}\)	Мультимегаваттный УКВ-радиовсплеск	Лавина RREA + усиление поля на остриях ледяных кристаллов
Отрицательный ступенчатый лидер	Облако–земля (шаги \(10\text{–}50\,\text{м}\))	\(10\text{–}30\,\text{мс}\) (паузы \(20\text{–}50\,\mu\text{s}\))	Средняя \(2 \times 10^5\,\text{м/с}\)	\(100\text{–}1000\,\text{А}\); ядро \(T \approx 5000\,\text{К}\)	Термический отлип (\(>1500\,\text{К}\)), space stems, фрактал DBM
Обратный удар (Return Stroke)	Земля–облако (диаметр ядра \(1\text{–}3\,\text{см}\))	Фронт \(1\text{–}5\,\mu\text{s}\); спад \(\sim 100\,\mu\text{s}\)	\(v \approx c/3 \approx 1,0 \times 10^8\,\text{м/с}\) (вверх)	Медиана \(30\,\text{кА}\) (+CG до \(300\,\text{кА}\)); \(T \approx 30\,000\,\text{К}\)	Волна нейтрализации, изохорный нагрев (\(10\text{–}100\,\text{атм}\)), цилиндрическая ударная волна
Спрайты / Джеты / Эльфы	\(15\text{–}100\,\text{км}\) (стратосфера–ионосфера)	\(<1\,\text{мс}\) (ELVES) до \(300\,\text{мс}\) (джеты)	\(10^5\,\text{м/с}\) до фазовой скорости \(>c\)	Холодная плазма (\(T_g \approx 200\text{–}300\,\text{К}\))	Квазиэлектростатический пробой, дисбаланс заряда, сферический ЭМИ
Земная гамма-вспышка (TGF)	\(10\text{–}15\,\text{км}\) (\(2\text{–}4\,\text{км}\) зимой)	\(<1\,\text{мс}\) (свечение \(\text{e}^+\) \(\sim 10\,\text{мин}\))	Релятивистские электроны \(v \approx 0,99c\)	Энергия фотонов до \(20\text{–}40\,\text{МэВ}\)	Тормозное излучение RREA, реакция \(^{14}\text{N}(\gamma,\text{n}){}^{13}\text{N}\) и аннигиляция \(\text{e}^+\text{e}^-\) (\(0,511\,\text{МэВ}\))

Сегодня системы 3D VHF Lightning Mapping Arrays (LMA) и широкополосные цифровые интерферометры реконструируют трёхмерную динамику молний внутри облаков с метровой точностью. В 2023 г. на вершине горы Зентис (\(2502\,\text{м}\), Швейцария) европейский консорциум с помощью тераваттного фемтосекундного лазера создал ионизированный канал за счёт лазерной филаментации и впервые в истории направил природную молнию вдоль лазерного луча на расстояние \(50\,\text{м}\) (Houard et al., 2023). В планетарном масштабе около 2000 одновременно действующих гроз питают глобальную электрическую цепь Вильсона, поддерживая потенциал ионосферы \(+250\,\text{кВ}\), поле хорошей погоды (\(100\text{–}150\,\text{В/м}\)) и резонансы Шумана (\(7,83\,\text{Гц}\)).