Лагранжиан плазмы: как метрическая теория выводит параметрическую модуляцию магнитных бурь

Лагранжиан плазмы: как метрическая теория выводит параметрическую модуляцию магнитных бурь

Проект «Рецепт хрупкого бессмертия» · блог-версия · 04.11.2026

Статус: вывод с гипотезами. Не проходила рецензирование. Требует проверки.

Пролог: откуда берётся формула бури

В статье «Магнитная буря 31 октября: пороговый срыв и метрические аномалии» мы использовали уравнение:

P̈_B + 2γṖ_B + ω₀²[1 + ξα(t)]P_B = F_CIR

Модуляция [1 + ξα(t)] была рабочей гипотезой. Она согласуется с данными, но не выведена из первых принципов. Рецензент сразу спросит: «Откуда формула? Почему ξα, а не ½ξα? Почему не ξ²α²?»

Эта статья — ответ. Мы выводим модуляцию из лагранжиана плазмы, с явными гипотезами и проверяемым количественным отличием. Главный результат: фактор ½ в модуляции. Феноменология даёт ξα, лагранжиан — ½ξα. Это дискриминатор: боковые полосы вдвое уже, чем в феноменологии. Проверяемо в данных магнитометров.

Но — с честными границами: ρ(α) — гипотеза, S_plasma — оценка, одна мода — упрощение. При физичном ξ = 10⁻⁴ эффект ничтожен — не видно ни в каком магнитометре. При демо ξ = 0.1 — виден в Шумане, не в Pc5.

Часть 1. Проблема: формула из воздуха

Уравнение бури выглядело так: колебания плазмы, затухание, собственная частота, и множитель [1 + ξα], который слегка меняет эту частоту. Поле α «дышит», и частота колебаний «дышит» вместе с ним.

Что не так. Множитель [1 + ξα] — это догадка. Мы её вставили, потому что она согласуется с данными. Но откуда она — не знали. А физика не любит формул из воздуха: если не можешь вывести — значит, не понимаешь.

Что нужно. Вывести этот множитель из более общего описания — из лагранжиана. Лагранжиан — это «энергетический портрет» системы: из него, как из семени, вырастают все уравнения движения. Если формула вырастает из лагранжиана — она не догадка, а следствие.

Ключевое предположение (гипотеза В.1). Мы уже знаем, что поле α меняет массу частиц по формуле «1 делить на корень из (1 + ξα)». Это подтверждено для фермионов. Мы предполагаем, что то же самое работает и для частиц плазмы — протонов, электронов, ионов. Это единственная точка отказа всей конструкции. Если она неверна — множитель ½ не следует из лагранжиана.

Запишем уравнение бури в явном виде.

P̈_B + 2γṖ_B + ω₀²[1 + ξα(t)]P_B = F_CIR

где P_B — давление плазмы (или смещение), γ — затухание, ω₀ — собственная частота, α(t) — поле α, F_CIR — внешняя сила (CIR — corotating interaction region). Модуляция [1 + ξα(t)] — вставка.

Ключевая гипотеза: конформная связь m(α) = m₀/√(1+ξα), подтверждённая для фермионов в fermion_mass.py, применима и к частицам плазмы. Это — предположение, не вывод из первых принципов. Но естественное: если α влияет на массу фермиона, почему не на массу частиц плазмы?

Часть 2. Лагранжиан Альвеновских колебаний

Представьте, что плазма — это газ из частиц. Плотность газа зависит от того, насколько тяжёлые частицы. Если поле α делает частицы легче, то при том же числе частиц в том же объёме газ становится менее плотным.

А частота Альвеновских волн (это колебания плазмы вдоль магнитных силовых линий) зависит от плотности: чем менее плотная плазма, тем быстрее волны. Суть в том, что частота связана с плотностью через квадратный корень.

Посмотрим, как это выглядит на уровне лагранжиана. Рассматриваем одну моду поверхностных Альвеновских волн (Pc5-пульсации) как осциллятор со смещением ξ(s, t) вдоль силовой линии:

L_plasma = ∫d³x [ ½ρ(α)ξ̇² − ½(B²/μ₀)(∂ξ/∂s)² ] + ∫d³x [ ½(∂α)² − V(α) ]

где ρ(α) — плотность плазмы, зависящая от α (ключевая гипотеза), B — магнитное поле, V(α) = ½m²α² — потенциал поля α (как в основной модели).

Границы гипотезы В.1:

• Квазинейтральность: n_e = n_i — сохраняется, если все массы масштабируются одинаково.

• Магнитное сжатие: B не зависит от α (предположение).

• Конечная температура: не учитывается (холодная плазма).

Часть 3. Откуда берётся ½

Смотрите цепочку:

1. Поле α делает массу частицы легче по формуле «1 делить на корень».

2. Плотность плазмы падает по той же формуле.

3. Скорость Альвеновских волн растёт — но уже по формуле «корень в четвёртой степени» (потому что скорость зависит от корня из плотности, а плотность — от корня).

4. Частота колебаний в квадрате растёт по формуле «корень» (потому что частота в квадрате зависит от скорости в квадрате).

И когда мы раскладываем этот «корень» в ряд (для малых изменений), там автоматически появляется ½. Это математическое свойство квадратного корня: √(1+x) ≈ 1 + x/2. Мы не выбирали ½ — оно вылезло само.

Если вернуться к формулам, цепочка выглядит так. При фиксированном числе частиц N и фиксированном объёме V:

ρ(α) = N·m(α)/V = ρ₀/√(1+ξα)

где ρ₀ = N·m₀/V. Тогда скорость Альвена:

v_A(α) = B/√(μ₀ρ(α)) = v_A0 · (1+ξα)^{1/4}

и частота:

ω_A²(α) = k²v_A²(α) = ω_A0² · (1+ξα)^{1/2}

При малых ξα (разложение Тейлора):

ω_A²(α) ≈ ω_A0² · (1 + ½ξα)

Фактор ½ — следствие ρ(α) ∝ (1+ξα)^{−1/2}, не вставка.

Это — не косметика. Это количественное уточнение: боковые полосы вдвое уже, чем в феноменологии.

Часть 4. Главное — не ½, а различие ½ против 1

Вот тут самое интересное. ½ против 1 — это не «вдвое уже». Это разные типы связи между полем α и материей.

Ширина боковой полосы в спектре — это отпечаток природы поля. Если мы когда-нибудь измерим эту ширину, мы узнаем не «какая модель верна», а «что такое α». Геометрическое поле или обычное?

Сравним два подхода в числах.

Фактор ½ — не косметика. Это количественное уточнение: боковые полосы вдвое уже, чем в феноменологии.

Часть 5. Числа: почему это пока не видно

Честно: при физичном значении силы связи (ξ = 10⁻⁴) эффект ничтожен. При демонстрационном значении (ξ = 0.1, в тысячу раз сильнее) — на грани.

Посмотрим на числа.

Честно:

• При физичном ξ = 1×10⁻⁴: эффект ничтожен (0.15 нГц для Pc5). Не видно ни в каком магнитометре.

• При демо ξ = 0.1: 0.15 мкГц для Pc5 (не видно), 0.4 мГц для Шумана (на грани — магнитометры чувствительны к ~мГц?).

В статье — диапазон, не одно число. «Наблюдаемое Δf ограничит ξ сверху».

Важное уточнение. Δf — это сдвиг несущей из-за DC-компоненты α, не расстояние между боковыми полосами. Боковые полосы — отдельный тест, требует m ≠ 2ω_A0 (не на резонансе). Разные интерпретации — разные числа, их нельзя путать.

Часть 6. Уравнения движения

До сих пор мы говорили о частоте. Но полная система — это два уравнения: одно для колебаний плазмы, другое для поля α. Они связаны: плазма раскачивает поле, поле меняет частоту плазмы.

Запишем оба уравнения движения.

Уравнение для ξ. Из вариации L_plasma по ξ:

ρ(α)ξ̈ + (B²/μ₀)(∂²ξ/∂s²) = F_CIR

Для одной моды (∂²ξ/∂s² → −k²ξ):

ξ̈ + ω_A0²(1 + ½ξα)ξ = F_CIR/ρ(α)

Это уже не феноменология — это вывод из лагранжиана. Модуляция ½ξα — следствие ρ(α).

Уравнение для α. Из вариации L_plasma по α:

□α + V'(α) = −∂ρ/∂α · ½ξ̇²

С ρ(α) = ρ₀(1+ξα)^{-1/2}:

∂ρ/∂α = −½ξρ₀(1+ξα)^{-3/2} = −½ξρ(α)/(1+ξα)

Тогда источник:

S_plasma = +¼ξρ(α)ξ̇²/(1+ξα) ≈ ¼ξρ₀ξ̇² (при малых ξα)

Уравнение: 

α̈ + m²α = S_plasma.

Проблема размерности. α̈ имеет размерность 1/с² (если α безразмерно). ρξ̇² — давление (Па = кг/(м·с²)). Это не 1/с². Решение: ввести явный коэффициент связи κ с размерностью 1/(кг·м):

S_plasma = κ · ¼ξρ₀ξ̇²

где κ — феноменологический параметр, аналог ξᵢ в твёрдом теле. Не выведен из первых принципов. Требует отдельного вывода (или экспериментального определения). Это честная граница, не провал.

Часть 7. Почему в плазме нет Aurora — мост к прогонам

В плазме условие резонанса (m ≈ 2Ω_prec) выполнено — это универсальное свойство уравнения Матьё. Но амплитуда не растёт. Почему?

Три причины — количественно:

1. Источник подавлен. κ = 1/M_α², при M_α ~ M_Pl это κ ~ 10⁻³⁸ ГэВ⁻². Плазменная плотность энергии, умноженная на 10⁻³⁸, не раскачает α до 0.1. Как если бы вы пытались раскачать качели, дуя на них.

2. Модуляция ничтожна. ε ≈ ½ξδα₀ ~ 10⁻³⁸. Уравнение Матьё с таким ε — гармонический осциллятор.

3. Кинетический член неканонический. ½M_α²(∂α)² с M_α ~ M_Pl подавляет амплитуду.

Мост к прогонам. Я запустил два прогона 48 часов и 132 часа. Они различаются по точности вычислений физического режима. В 48/132 источник механический — вращение тела, а не плазма. И он не подавлен. Поэтому Aurora возможна в принципе — в отличие от плазмы. Но там своя угроза — детюнинг: когда поле раскачивается, оно само сдвигает частоту, и система «уходит» с резонанса. Как радио, которое расстраивается, как только вы прибавляете громкость. Плазма объясняет, почему прогоны имеют шанс, и почему шанс не гарантирован.

Часть 8. Честные границы: где мы всё ещё работаем

Это не «окончательный» лагранжиан. Это более физичная модель, которая:

• Согласуется с сегодняшними результатами (резонанс в амплитуде, порог).

• Связывает бурю с твёрдым телом через общую структуру (параметрический резонанс).

• Даёт количественное предсказание (множитель ½).

Вот полная картина со статусами.

Проблема размерности (повторно, потому что важно): коэффициент κ с размерностью 1/(кг·м) введён, чтобы уравнение для α было размерностно согласовано. Но κ — феноменологический параметр, не выведенный из первых принципов. Это честная граница, не провал.

Часть 9. Что на данный момент нельзя сказать:

Часть 10. Карта достоверности — простыми словами

Гипотеза (не выведено):

• Поле α меняет массу частиц плазмы по той же формуле, что и для фермионов.

• Источник в уравнении для α имеет форму κ · ¼ξρ₀ξ̇² — с феноменологическим κ.

• Внешняя сила F_CIR не зависит от α.

• Одной моды достаточно для описания.

Следствие (если гипотеза верна):

• ω_A²(α) ≈ ω_A0²(1 + ½ξα) — множитель ½.

• Ширина боковых полос вдвое уже, чем в феноменологии.

• При физичном ξ эффект не виден; при демо ξ виден в Шумане, не в Pc5.

Проверяемо (если чувствительность вырастет):

• Сдвиг несущей частоты Pc5/Шумана во время бури (INTERMAGNET, GOES).

• Боковые полосы на ω_A0 ± m (при m ≠ 2ω_A0).

• Пороговое поведение при α_crit.

• Резонанс в амплитуде (⟨P_α⟩ ≈ 0) — раскачка амплитуды без накопления.

Ограничение: при физичном ξ = 1×10⁻⁴ ни один из этих тестов не виден в существующих данных. Нужны либо усиленные ξ (демо), либо более чувствительные магнитометры, либо другой канал (не Pc5).

Эпилог: шаг вперёд, не финал

В предыдущей статье: ω²[1+ξα(t)] — феноменология.

Теперь: ω_A²(α) = ω_A0²(1+ξα)^{1/2} ≈ ω_A0²(1 + ½ξα) — вывод из ρ(α).

Дискриминатор: боковые полосы вдвое уже (½ξα vs ξα). Проверяемо в данных магнитометров.

Границы: ρ(α) — гипотеза, S_plasma — оценка с κ, одна мода — упрощение. При физичном ξ=10⁻⁴ эффект не виден; при демо ξ=0.1 — виден в Шумане, не в Pc5.

Это шаг вперёд — от феноменологии к лагранжиану, с честными границами и проверяемым уточнением. Не «окончательный» лагранжиан — более физичная модель. И если фактор ½ подтвердится — это ограничит ξ сверху или подтвердит конформную связь для плазмы. Оба продвигаемый вперёд — результат.

Главный вывод — одной фразой

Мы вывели множитель ½ в формуле для магнитной бури из конформной связи — это не вставка, а следствие квадратного корня. Главное — не ½, а различие ½ против 1: это дискриминатор типа связи (геометрическое поле vs обычное). Но ρ(α) — гипотеза, источник — оценка с κ, при физичном ξ эффект не виден. Это шаг вперёд от феноменологии к лагранжиану, с честными границами и проверяемым уточнением.

Что осталось

1. Численная проверка — решить систему для ξ и α с κ. Пока не сделано.

2. Наблюдательная проверка — боковые полосы в Pc5/Шумане. Нужны более чувствительные магнитометры или больше событий.

3. Теоретическая проверка — вывести ρ(α) для плазмы из первых принципов. Пока не сделано.

Каждый шаг — реальный, не «симуляционный бег».

Ссылки на материалы цикла

1. «Магнитная буря 31 октября: пороговый срыв и метрические аномалии» — феноменологическая модель, которую это приложение уточняет.
2. «Скомканный лист: почему мы видим только мюоны и электроны» — конформная связь m(α), рекурсия опровергнута.
3. «Скомканный лист II: лёгкие электроны» — Aurora-режим, четыре состояния, порог ω₁/ω₃ > 0.29.
4. «Почему NICER не увидел сдвиг массы в глитчах» — три механизма маскировки, честное признание.
5. «Метрическая модель динамики твёрдого тела: параметрический резонанс как проверяемое следствие» — лагранжиан, уравнения Эйлера–Пуанкаре.
6. «Вселенная, которая дышит: уравнения Фридмана как мост между ОТО и метрической теорией» — математическое ядро.
7. «Квантование метрики граничными условиями: эффект Казимира, лестница состояний и рубцевание» — почему спектр дискретен.

«Рецепт хрупкого бессмертия» — цикл статей. Настоящая статья предлагает концепцию к рассмотрению. Не проходила рецензирование. Это гипотеза, которая требует проверки. Фактор ½ — количественное предсказание. ρ(α) — гипотеза. S_plasma — оценка. Проверка — INTERMAGNET, GOES. При физичном ξ эффект не виден.

«Физический» ≠ «доказанный». Меньше произвола, больше проверяемых следствий. Флаг стоит: точка на карте, а не доказательство.