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

Автор: Александр Алексеенко

НАСЫЩЕНИЕ ВМЕСТО РЕЗОНАНСА: ОТ ЛАБОРАТОРНОЙ ГАЙКИ ДО ГРАВИТАЦИОННЫХ ВОЛН ПУЛЬСАРОВ

Из цикла «Рецепт хрупкого бессмертия»
Автор: Александр Алексеенко
Статус: рабочий протокол. Не проходила рецензирование. Требует проверки.


Пролог: не пугайтесь

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

  1. Есть гайка. Она крутится.
  2. Мы хотели узнать, не начнёт ли она крутиться сильнее, если внутри неё есть кое-что странное.
  3. Она не начала. А потом мы прыгнули от гайки к звёздам — и там кое-что нашли.

Всё остальное — детали. Поехали.


Кто такой лагранжиан и почему его все боятся

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

Это как рецепт борща. Можно рассказывать: «возьмите свёклу, потрите, добавьте капусту, варите двадцать минут, не забудьте сметану». А можно написать одну строчку: «Борщ = свёкла + капуста + сметана, всё вместе, варить до готовности». Второе — компактнее, но требует, чтобы вы знали, что такое «свёкла» и «до готовности».

Лагранжиан — это такая строчка. Если вы её понимаете — вы понимаете всю систему. Если нет — выглядит как заклинание. Ничего страшного. Дальше мы будем объяснять словами, а лагранжиан оставим в приложении, чтобы физики не обиделись.


Гайка, которая крутится не так, как должна

В прошлой лекции мы взяли гайку — обычную, шестигранную, из тех, что теряются в гараже. И задали ей странный вопрос: а что, если моменты инерции этой гайки зависят от локального поля α?

Стоп. Ещё проще.

Момент инерции — это «насколько трудно раскрутить предмет». У длинной палки он больше, у короткой — меньше. У гайки он зависит от того, как она расположена: вдоль одной оси крутить легче, вдоль другой — труднее.

Мы предположили: а что, если это «легче-труднее» меняется? Не потому, что гайка деформируется, а потому, что само пространство внутри неё немного другое — в нашей модели это описывается полем α.

И вот тут начинается интересное. Если моменты инерции зависят от α, а α зависит от того, как гайка крутится, — получается замкнутый круг:

Гайка крутится → поле α раскачивается → моменты инерции меняются → гайка крутится иначе → поле α раскачивается сильнее → …

Это как микрофон, поднесённый к колонке. Звук идёт в микрофон, усиливается, выходит из колонки, снова попадает в микрофон — и если вовремя не убавить громкость, получается тот самый визг.

В физике это называется параметрический резонанс. Мы искали именно его.


Что мы искали (и почему это было бы круто)

Если «микрофон-колонка» работают в резонансе — амплитуда растёт экспоненциально. Это значит: сначала чуть-чуть, потом сильнее, потом ещё сильнее, и очень быстро — очень сильно.

В нашей модели это означало бы: гайка начала крутиться чуть быстрее — и раскрутилась бы сама, без внешней силы. Только за счёт того, что поле α внутри неё «подстраивается».

Звучит как вечный двигатель? Не совсем. Энергия не берётся из ниоткуда — она переходит из поля α в движение. Но выглядит эффектно: система сама себя раскачивает.

Мы хотели это увидеть. Запустили численную симуляцию. Ждали.


Что мы нашли (и почему это не провал)

Резонанса не было.

Точнее, был, но не такой, как мы ждали. Вместо экспоненциального взрыва — линейный рост. Вместо «гайка раскручивается сама» — «гайка крутится чуть быстрее, но очень медленно, и поле α внутри неё накачивается до определённого уровня — и останавливается».

Это называется насыщение. И это, как ни странно, хорошая новость.

Почему? Потому что в реальном мире ничего не улетает в бесконечность. Если бы наша модель предсказывала бесконечный резонанс — это была бы модель, которая «взрывается» при малейшем возмущении. А реальные системы так себя не ведут. Они самоограничиваются: чем сильнее раскачиваются, тем сильнее что-то их тормозит.

В нашей модели этот «тормоз» встроен: чем больше α, тем сильнее меняются моменты инерции, тем сильнее сдвигается частота — и система сама уходит с резонанса. Как микрофон, который «слышит» колонку и автоматически убавляет громкость.

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

Пока честно: экспоненциального резонанса нет. Это устойчивый результат. А почему его нет — вопрос открытый.


Прыжок: от гайки к звёздам

А теперь — самое интересное. Мы подумали: а что, если ту же модель применить не к гайке, а к нейтронной звезде?

Нейтронная звезда — это то, что остаётся от массивной звезды после взрыва сверхновой. Размером с город, массой с Солнце. Одна чайная ложка её вещества весит как гора Эверест. И она крутится — иногда сотни оборотов в секунду.

Момент инерции у неё не как у гайки — на 44 порядка больше. Это как сравнивать вес пылинки и вес Луны. Но уравнения — те же. Если наша модель верна для гайки, она должна работать и для звезды.

И вот тут — первый сюрприз. При параметрах, которые для ядерной материи физичны (то есть не «мы подогнали, чтобы красиво», а «так должно быть при таких плотностях»), модель даёт наблюдаемый порядок для двух вещей, которые астрофизики видят прямо сейчас.


Глитчи: звёздные «икки»

Пульсар — это нейтронная звезда, которая мигает с очень точной периодичностью, как часы. Астрофизики используют её как эталон времени.

Но иногда — бац — и пульсар на мгновение ускоряется. Частота скачет на миллионную долю, потом восстанавливается. Это называется глитч.

Классическое объяснение: внутри звезды есть сверхтекучее ядро, которое накапливает угловой момент, а потом «отпускает» его в кору. Как если бы вы раскрутили юлу, а внутри неё сидел маленький волчок, который иногда «пинает» юлу изнутри.

Мы говорим: а может быть, это метрический сдвиг? Если поле α внутри коры звезды меняется на микроскопическую величину — момент инерции меняется, а поскольку момент импульса сохраняется, частота вращения скачет. И это даёт тот же порядок, что наблюдается: Δν/ν ~ 10⁻⁶.

Мы не говорим, что классика неправа. Мы говорим: есть второй способ получить тот же результат. И у него есть отличительный признак — см. следующий раздел.


Гравитационные волны: звёзды, которые «шепчут»

Если поле α внутри звезды колеблется, то форма звезды тоже колеблется — чуть-чуть, на микроны. А по общей теории относительности Эйнштейна, любое колеблющееся тело излучает гравитационные волны — рябь пространства-времени, которая распространяется со скоростью света.

Мы посчитали: на какой частоте звезда должна «шептать»? Ответ: на частоте m = 2Ω_prec — удвоенной частоте прецессии. Для PSR J0537-6910 это около 1.2 Гц (если кора мягкая) или 12.4 Гц (если жёсткая). Для Vela — 0.2 Гц или 2.2 Гц.

И вот тут — важный момент. Классическая сверхтекучесть объясняет глитч, но не предсказывает гравитационную волну на этой конкретной частоте. Наша модель предсказывает и то, и другое — и, главное, корреляцию: был глитч → появился пик на частоте m.

Это проверяемо. Не «когда-нибудь потом», а прямо сейчас: в данных LIGO, Einstein Telescope, Cosmic Explorer. Если пик есть — модель получает поддержку. Если нет — мы ограничиваем параметры ξ и δ сверху. Оба исхода — результат.


Где искать (и почему не там, где вы думаете)

Мы построили карту чувствительности: для каждого пульсара — частота m, для каждого детектора — диапазон. И вот что оказалось:

  • Vela — цель для Einstein Telescope (ET), не для LIGO. LIGO туда просто не смотрит.
  • PSR J0537-6910 — тоже цель для ET/CE. Только если кора очень жёсткая (δ ~ 0.1), частота попадает на нижний край LIGO (12.4 Гц).

И тут мы честно признаём ошибку: в первых оценках мы взяли расстояние до PSR J0537-6910 5 кпк и получили амплитуду ~3×10⁻²⁶ — «на пороге LIGO!». Потом выяснилось, что пульсар находится в Большом Магеллановом Облаке, 49.6 кпк — в десять раз дальше. Амплитуда упала в десять раз. Не LIGO, а ET. Это не «мы нашли сигнал» — это «мы говорим, где искать, и честно признаём, что ошиблись в расстоянии».


Честный финал: что мы НЕ доказали

Чтобы не быть как те учёные из мемов, которые «нашли пятно на фотографии Марса и объявили о жизни», перечислим что мы НЕ утверждаем:

  1. Мы не доказали, что поле α существует. Это гипотеза. Может, его нет — тогда всё это красивая математика.
  2. Мы не доказали, что глитчи — это метрический сдвиг. Мы показали, что это возможно и даёт правильный порядок. Но классика тоже даёт правильный порядок.
  3. Мы не открыли гравитационные волны от пульсаров. Мы сказали: вот частота, вот амплитуда, ищите. Если найдёте — отлично. Если нет — мы ограничим параметры.
  4. Мы не утверждаем, что насыщение в гайке — физика. Может, это артефакт метода. Следующий прогон покажет.
  5. Мы не знаем δ и ε. Они зависят от уравнения состояния ядерной материи, которое само по себе — открытая проблема. Поэтому мы даём диапазон, а не точку.

Что мы утверждаем:

  • Модель работает: уравнения интегрируются, энергия сохраняется, размерности сходятся.
  • В физичном лабораторном режиме экспоненциального резонанса нет — это устойчивый результат.
  • В астрофизическом режиме модель даёт наблюдаемый порядок для глитчей и гравитационных волн.
  • Есть разделительный тест: корреляция «глитч → пик на m». Это проверяемо.


Шутка на финал

Говорят, физики делятся на три категории:

  1. Те, кто понимает лагранжиан.
  2. Те, кто делает вид, что понимает.
  3. Те, кто пишет научно-популярные статьи.

Я отношусь к четвёртой: тем, кто пишет статьи с формулами, а потом объясняет их на пальцах — потому что сам когда-то был во второй категории.

Так что если вы дочитали до конца и всё ещё не знаете, кто такой лагранжиан, — добро пожаловать в клуб. Главное, что вы теперь знаете: гайка крутится, звезда мигает, а мы ищем, где они «шепчут». Остальное — детали.

От гайки до звёзд — 44 порядка. А от «не понимаю» до «примерно понимаю» — один спойлер.


Кристалл памяти обновлён: 01.11.2026. Проект «Рецепт хрупкого бессмертия». Версия для гуманитариев. Формулы — в основной статье, шутки — здесь.


Пролог: два полигона одной модели

В предыдущей лекции мы вывели минимальную метрическую модель динамики свободного твёрдого тела. Лагранжиан, уравнения Эйлера–Пуанкаре с динамическими моментами инерции, уравнение поля α — всё это дало проверяемое предсказание: параметрический резонанс на частоте 2Ω_prec, которого нет в классическом пределе ξ = 0.

Код был проверен. Энергия сохранялась с точностью 3.8×10⁻¹⁴. Отношение частот в диагностике дало f_α / f_ω₁ = 2.000 — точную подпись резонанса. Демонстрационный прогон при ξ = 0.1 показал рост max|ω₁| с насыщением. Всё сходилось.

Тогда мы задали два вопроса, и каждый потребовал своего полигона:

  1. Что будет в физичном лабораторном режиме?ξ = 1×10⁻⁴, I₂ = 1.1, m = 2Ω_prec = 2.0. Это Часть I — лабораторная гайка.
  2. Что будет, если масштабировать на 44 порядка? Нейтронная звезда, I ~ 10³⁸ кг·м², ξ ~ 0.5–1.0 (ядерная материя). Это Часть II — астрофизический полигон.

Мы провели длинный прогон в лабораторном режиме (60 ч настенного времени, 14.8% от t_end). Он дал устойчивый отрицательный результат — экспоненциального резонанса нет — и выявил ограничения метода: разреженную выборку и большие начальные условия. По результатам прогона код скорректирован для достижения более точного результата: плотная выборка, малые начальные возмущения, ограничение шага интегратора. Следующий прогон ответит на вопрос, была ли причина отсутствия резонанса физической или методической.

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


ЧАСТЬ I. НАСЫЩЕНИЕ ВМЕСТО РЕЗОНАНСА

1.1. В чём уникальность подхода

Стандартная теория параметрического резонанса работает с внешней периодической накачкой: уравнение Матьё, ẍ + ω₀²(1 + ε cos νt)x = 0, где ν — частота внешнего поля. Резонанс возникает, когда ν ≈ 2ω₀.

Наш случай другой. У нас нет внешней накачки. Поле α само является динамической переменной, которая:

  • возбуждается вращением через источник ½ΣξᵢIᵢ⁰ωᵢ²;
  • управляет моментами инерции через Iᵢ(α) = Iᵢ⁰(1 + ξᵢα);
  • влияет на вращение через уравнения Эйлера–Пуанкаре.

Это замкнутая петля: вращение → α → I(α) → вращение. Резонанс здесь не навязывается извне — он самовозбуждается, если частота источника (2Ω_prec) попадает в собственную частоту поля (m). Такая постановка — автономный параметрический резонанс с самосогласованной обратной связью — в литературе встречается редко.

Второе отличие: нелинейность встроена. В стандартной теории резонанс — это неустойчивость, амплитуда растёт до бесконечности. У нас связь Iᵢ(α) и источник ωᵢ² нелинейны, поэтому модель может ограничивать амплитуду — без трения, без обрезания. Может ли — это и есть вопрос, на который отвечает прогон.

Третье отличие: мы работаем в физичном режиме (ξ = 1×10⁻⁴), а не в усиленном (ξ = 0.1). Результат — не «демонстрация механизма», а граница применимости.


1.2. Что искали: линейное предсказание

Iᵢ(α) = Iᵢ⁰(1 + ξᵢα)

α̈ + m²α = ½ Σ ξᵢ Iᵢ⁰ ωᵢ²

При m = 2Ω_prec в линейной теории возникает параметрический резонанс: ω₁ растёт экспоненциально, exp(γt), где γ ≈ ξ·α·Ω_prec/2. Для физичных параметров γ ~ 4×10⁻⁷, t_e ~ 2.5×10⁶. За t ≈ 7×10⁸ это сотни времён t_e. Если бы резонанс был «чистым», ω₁ должна была вырасти на порядки. Этого не произошло.


1.3. Что показал пилотный прогон — и его ограничения

Физичный прогон за 60 ч настенного времени (14.8% от t_end):

ПараметрНачало (t ~ 10³)Конец (t ≈ 7.1×10⁸)Изменение
max|ω₁|0.36050.3624+0.52% (сопоставимо с разбросом ±1.2%)
max|α|1.56×10⁻³ (переходный)8.3×10⁻³DC-уровень
E_drift8.6×10⁻⁹6.5×10⁻⁴линейный

Что установлено устойчиво, а что — нет:

  1. Экспоненциального роста нет. Если бы резонанс был с γ ~ 4×10⁻⁷, за 7×10⁸ рост был бы e²⁸⁰ — его нельзя пропустить ни при какой выборке. Его нет. При этом разброс между чанками ±1.2% (0.3563–0.3624) означает, что тренд +0.5% сопоставим с шумом: линейность не установлена, установлена только граница γ < 4×10⁻¹¹.
  2. α накачивается до ~8×10⁻³ (DC-уровень). Это устойчиво: источник ½ΣξᵢIᵢ⁰ωᵢ² работает. Но назвать это «плато» или «насыщением» нельзя — три сценария (установившийся баланс, медленный рост, резонансная накачка) дают похожие числа в разреженной выборке.
  3. E_drift линейный и управляемый. 6.5×10⁻⁴ за 14.8% — дрейф интегратора, не физика. Меньше сигнала, но растёт с той же скоростью.


Рис. 1. физичный прогон за 60 ч (14.8% от t_end), усиление ×8. Четыре панели: max|ω₁| — колебания вокруг 0.3605–0.3624; рост max|ω₁| в % (серая тонкая — реальный, зелёная жирная — ×8); max|α| — накачивается до ~8.3×10⁻³; E_drift — линейный, 6.5×10⁻⁴. Подпись: α накачивается, ω₁ визуально растёт, но тренд сопоставим с разбросом между чанками (±1.2%). Рисунок показывает, что видели, а не что доказали.


1.4. Почему это не резонанс — и что это может быть

Два механизма могли бы останавливать рост раньше, чем экспонента успевает проявиться.

Первый: нелинейный сдвиг частоты. Ω_prec зависит от α, а α — динамическая переменная. Когда α колеблется, мгновенная частота прецессии Ω_prec(t) отличается от линейной. Резонанс «детюнится» — система уходит с резонансной частоты автоматически, как только амплитуда становится достаточной. Это классика нелинейных колебаний: резонанс существует, но его частота сдвинута на величину порядка q·Ω или q²·Ω, где q = ξα.

Второй: обратная связь через ωᵢ². Источник в уравнении α растёт как квадрат амплитуды. Когда ω₁ увеличивается, источник усиливается, α растёт, моменты инерции меняются сильнее — и это ещё больше сдвигает частоту. Система могла бы найти новый баланс: предельный цикл, где накачка и диссипация уравновешиваются.

В фазовом пространстве это могло бы быть не спиралью наружу, а замкнутой траекторией. Это был бы предельный цикл. Но чтобы это утверждать, нужна плотная выборка и квазилинейный режим — их в этом прогоне нет.


Рис. 2.Сравнение (12ч Vs 60 ч). За 5× времени рост max|ω₁| вырос в ~4.9× (+0.107% → +0.524%). Это не экспоненциальный рост (тот дал бы e⁵ ≈ 148×). Но линейный ли он — не установлено: абсолютные значения сопоставимы с разбросом. Устойчиво только: рост не экспоненциальный.


Вот какие выводы мы можем сделать:

Пилотный прогон при ξ = 1×10⁻⁴, I₂ = 1.1, m = 2Ω_prec (14.8% от t_end) не показал экспоненциального роста ω₁. Разброс между чанками (±1.2%) ограничивает инкремент γ < 4×10⁻¹¹ — на три порядка ниже линейной оценки. Энергия сохраняется линейно (E_drift ~ 6.5×10⁻⁴). α накачивается до ~8×10⁻³ — источник работает. Однако разреженная выборка (n_pts_chunk = 500 на dt = 10⁶) и большие начальные условия (ω₁(0)=0.3) не позволяют отличить медленный рост от шума и «плато» от резонансной накачки. Количественные утверждения о линейном росте и предельном цикле не делаются — они проверяются в следующем прогоне.

Это не «модель не работает». Это модель работает, но результат требует уточнения.


1.5. Почему это важно

В линейной теории резонанс — катастрофа: амплитуда растёт, пока система не разрушится. В реальности этого не происходит, потому что нелинейности всегда ограничивают рост. Наша модель, по-видимому, воспроизводит это из первых принципов: не потому, что мы добавили трение, а потому, что связь Iᵢ(α) и уравнение поля сами создают обратную связь.

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


Рис. 3. 16 ч, точный график. Виден провал −1.312% — это биения, а не «резонанс кончился». Огибающая колеблется, но тренд вверх. Если смотреть только на max|ω₁| в моменте, провал можно принять за остановку роста.


1.6. Что дальше в Части I

Отрицательный результат при m = 2Ω_prec не означает, что резонанса нет нигде. Он означает, что при этом m и этих ξ система не показывает экспоненциального роста. Два способа проверить, где резонанс есть:

Сканер по m. Короткие симуляции при m в окрестности 2Ω_prec (±10%, 21 точка), смотреть ENV-рост. Если пик есть при m ≠ 2Ω_prec — нелинейный сдвиг сдвинул резонанс. Если нет — резонанс подавлен во всей окрестности.

Усиленный режим. ξ = 1×10⁻² (в 100 раз сильнее). Резонанс должен проявиться за минуты. Это не подгонка — стандартная практика: сначала показать механизм в усиленном режиме, потом искать в физичном.

Скорректированный код уже готов. По результатам длинного прогона мы изменили три вещи: (1) плотная выборка (dt_eval = 1.0 вместо 500 точек на чанк 10⁶) — теперь max|ω₁| это реальный максимум, а не выборочная точка; (2) малые начальные (ω₁ = 0.01, ω₂ = 0.005) — квазилинейный режим, где резонанс должен проявиться как экспонента, если он есть; (3) max_step = 0.3 — интегратор не залипает на жёстких участках, данные непрерывны. Следующий прогон ответит на вопрос: была ли причина отсутствия резонанса физической или методической.


ЧАСТЬ II. АСТРОФИЗИЧЕСКИЙ ПОЛИГОН

2.1. Зачем прыгать на 44 порядка

В Части I мы остановились на границе: физичный лабораторный прогон не показал экспоненциального резонанса, но причина требует проверки. Независимо от этого, если уравнения верны, они должны работать не только для гайки с I ~ 10⁻⁶ кг·м², но и для нейтронной звезды с I ~ 10³⁸ кг·м². Это масштабирование на 44 порядка.

И вот что оказывается: при параметрах, которые для ядерной материи физичны (ξ ~ 0.5–1.0, сильная нелинейность EOS), модель даёт наблюдаемый порядок для двух вещей:

  1. Глитчи пульсаров — скачки частоты на Δν/ν ~ 10⁻⁶. Классика объясняет их сверхтекучестью. Мы говорим: это может быть метрический сдвиг инерции, и он даёт тот же порядок.
  2. Непрерывные гравитационные волны на частоте m = 2Ω_prec. Классика их не предсказывает. Мы предсказываем — и говорим, где искать: не «пятно», а пик на m, коррелирующий с глитчами.

Это не «объяснение глитчей». Это наблюдательная программа: вот частота, вот амплитуда, вот разделительный тест. И честная граница: δ и ε мы не знаем точно, поэтому даём диапазон, а не точку.


2.2. Масштабирование: от лаборатории к звезде

Уравнения — те же. Меняются числа:

ПараметрГайка (Часть I)Нейтронная звездаПорядок
I10⁻⁶ кг·м²10³⁸ кг·м²10⁴⁴
ω₃10² рад/с10³ рад/с10
Ω_prec2π · Гцдоли Гц—
m = 2Ω_prec~2 рад/с0.1–100 рад/с—
ξ1×10⁻⁴ (физ.) / 0.1 (демо)0.5–1.0 (ядерная материя)—

Частота прецессии: Ω_prec² = (I₃−I₁)(I₃−I₂)ω₃²/(I₁I₂).

Для нейтронной звезды δ = (I₃−I₁)/I₁ — относительная разность моментов инерции (кора vs сверхтекучее ядро). Тогда:

Ω_prec ≈ δ · ω₃, m = 2Ω_prec ≈ 2δ · ω₃

Вся неопределённость — в δ. Для ядерной материи δ может быть от 10⁻³ до 10⁻¹. Значит, m для одного и того же пульсара — от 0.1 до 12 Гц, и попадает в разные детекторы.

Рис. 4.«Карта чувствительности»: по горизонтали δ (10⁻³–10⁻¹), по вертикали m = 2Ω_prec для Vela, Crab, PSR J0537-6910. Поверх — зоны LIGO (10–1000 Гц), Einstein Telescope (1–10⁴ Гц), Cosmic Explorer (0.1–10⁴ Гц). Видно: Vela и PSR J0537-6910 — цели для ET/CE, не для LIGO. Только при δ ~ 0.1 PSR J0537-6910 попадает на нижний край LIGO (m ~ 12.4 Гц).

Вывод из карты: мы не можем «подогнать» δ под LIGO. Мы должны честно сказать: «Если кора жёсткая (δ ~ 0.1), смотрим LIGO. Если мягкая (δ ~ 10⁻²), смотрим ET».


2.3. Глитчи как метрический сдвиг инерции

Что наблюдается. Vela: Δν/ν ~ 10⁻⁶ (до 3×10⁻⁶). PSR J0537-6910: Δν/ν ~ 0.36×10⁻⁶. Crab: 10⁻⁸–10⁻⁶.

Что говорит классика. Сверхтекучее ядро накапливает угловой момент, затем «отпускает» его в кору — кора резко ускоряется. Это работает, но не предсказывает гравитационного сигнала на конкретной частоте.

Что говорим мы. Если в коре происходит локальный метрический сдвиг Δα, момент инерции меняется на ΔI/I ≈ ξ·Δα. Момент импульса L = Iω сохраняется, значит:

Δω/ω ≈ −ΔI/I ≈ −ξ·Δα

Чтобы получить Δν/ν ~ 10⁻⁶: при ξ ~ 0.1 нужно Δα ~ 10⁻⁵; при ξ ~ 1 нужно Δα ~ 10⁻⁶.

Это физично. Δα ~ 10⁻⁶ — микроскопический сдвиг метрики, не «взрыв». Перераспределение сверхтекучей жидкости или напряжение в коре могут его вызвать.

Вот какие выводы мы можем сделать:

При ξ ~ 0.5–1.0 (ядерная материя, сильная нелинейность EOS) наблюдаемые глитчи Δν/ν ~ 10⁻⁶ соответствуют метрическому сдвигу Δα ~ 10⁻⁶–10⁻⁵. Это не «объяснение глитчей», а порядковая оценка: модель способна дать наблюдаемый порядок при физичных параметрах. Отличие от классики — не в Δν/ν, а в сопровождающем гравитационном сигнале.


2.4. Гравитационные волны: m = 2Ω_prec

Если поле α(t) внутри пульсара колеблется, квадрупольный момент Qᵢⱼ ~ Iᵢⱼ становится периодической функцией времени. По ОТО, переменный квадруполь излучает гравитационные волны. Частота — m ≈ 2Ω_prec.

Амплитуда на Земле: h₀ = (4G/c⁴) · (I₃ · ε · ω² / r)

где ε = (I₃−I₁)/I₃ — эллиптичность, r — расстояние. В нашей модели ε связано с ξ и Δα: при ξ ~ 1 и Δα ~ 10⁻⁶ достижимо ε ~ 10⁻⁵; при ξ ~ 0.1 — ε ~ 10⁻⁶.

Vela (PSR B0833-45)

  • ω ≈ 70.4 рад/с, I ~ 10³⁸ кг·м², r ≈ 287 пк = 8.86×10¹⁸ м
  • δ = 10⁻²: m ≈ 0.22 Гц (ET/CE)
  • δ = 10⁻¹: m ≈ 2.2 Гц (ET/CE, не LIGO)
  • ε = 10⁻⁶: h₀ ≈ 4.6×10⁻²⁸ (ET)
  • ε = 10⁻⁵: h₀ ≈ 4.6×10⁻²⁷ (ET, на грани)

Vela — цель для Einstein Telescope, не для LIGO.

PSR J0537-6910

  • ω ≈ 392 рад/с, I ~ 10³⁸ кг·м², r ≈ 49.6 кпк = 1.53×10²¹ м (в 173 раза дальше Vela)
  • δ = 10⁻²: m ≈ 1.24 Гц (ET/CE)
  • δ = 10⁻¹: m ≈ 12.4 Гц (LIGO, нижний край)
  • ε = 10⁻⁶: h₀ ≈ 8.2×10⁻²⁹ (ET)
  • ε = 10⁻⁵: h₀ ≈ 8.2×10⁻²⁸ (ET, на грани)

Важная поправка: в первых оценках мы взяли r = 5 кпк и получили h₀ ~ 3×10⁻²⁶. Это была ошибка: PSR J0537-6910 находится в Большом Магеллановом Облаке, r ≈ 49.6 кпк. Амплитуда падает в 10 раз, и цель — ET/CE, не LIGO.

Вот какие выводы мы можем сделать:

Для Vela и PSR J0537-6910 при физичных ξ ~ 0.5–1.0 и δ ~ 10⁻²–10⁻¹ модель предсказывает непрерывные гравитационные волны на частоте m = 2Ω_prec с амплитудой h₀ ~ 10⁻²⁸–10⁻²⁷. Это на уровне чувствительности Einstein Telescope и Cosmic Explorer, но ниже порога LIGO (кроме PSR J0537-6910 при δ ~ 0.1, где m ~ 12.4 Гц — нижний край LIGO). Это не «мы нашли сигнал», а «мы говорим, где искать»: пик на m, а не пятно.

Рис. 5.Огибающая h(t) после глитча (безразмерная симуляция PSR J0537-6910, δ=0.01, ξ=0.1). Сигнал не чистый синус, а биения — α и ω обмениваются энергией. Длительность переходного процесса — несколько периодов релаксации α. В безразмерной модели это ~10 периодов α; в физических единицах для PSR J0537-6910 при ξ = 0.1 это секунды (t_end ≈ 12 с), а не дни. Для дней–недель нужен ξ ~ 10⁻⁸ — другой режим, который в этой симуляции не показан. Форма огибающей — предсказание модели: после глитча сигнал нарастает, затем затухает по мере релаксации α.


2.5. Разделительная диагностика

Классические модели сверхтекучести тоже дают Δν/ν ~ 10⁻⁶. Как отличить?

ПризнакКлассика (сверхтекучесть)Метрика (наша модель)
Глитч Δν/νДа, ~10⁻⁶Да, ~10⁻⁶ (при ξΔα ~ 10⁻⁶)
Зависимость от температурыДа (сверхтекучесть чувствительна к T)Нет (зависит от геометрии/напряжений)
Непрерывная волна на 2Ω_precНетДа
Корреляция «глитч → пик на m»Не предсказываетПредсказывает

Вот какие выводы мы можем сделать:

Классические модели сверхтекучести объясняют глитчи, но не предсказывают непрерывный гравитационный спектр на частоте 2Ω_prec. Наша модель предсказывает: сразу после глитча на частоте m должен появиться когерентный сигнал, который затухает по мере релаксации α(t). Это проверяемое предсказание: искать не «шум», а корреляцию — «был глитч, появился пик на m». Если LIGO/ET не видят пика на m после глитча, это ограничивает ξ и δ сверху.


Рис. 6. Схема разделительной диагностики (не данные, не симуляция — иллюстрация предсказания). Панель 1: ν(t) с глитчами — обе модели объясняют. Панель 2: классика (сверхтекучесть) — пика на m нет. Панель 3: метрика — после каждого глитча пик на m, затухающий по мере релаксации α(t). Разделительный тест: корреляция «глитч → пик на m» в данных LIGO/ET/CE. Схема показывает, что искать, не утверждая, что сигнал уже найден.


2.6. Честные границы

Первое: δ не известен точно. δ зависит от EOS. Диапазон 10⁻³–10⁻¹, а не точка.

Второе: ε не известен точно. ε зависит от ξ и Δα. Диапазон 10⁻⁶–10⁻⁵, а не точка.

Третье: ξ ~ 1 — гипотеза. При α ~ 1 инерция может менять знак — риск нефизичности. ξ = 0.5 — компромисс, помечен как гипотеза.

Четвёртое: это порядковые оценки, не доказательство. Мы говорим, где искать, а не «мы нашли».

Пятое: численная симуляция — безразмерная. Для пульсара мы работаем в безразмерных единицах, чтобы не было переполнения float64. Амплитуда h(t) — грубая оценка, не точное значение.


ОБЩИЙ ИТОГ

Часть I показала: в физичном лабораторном режиме (ξ = 1×10⁻⁴) экспоненциального резонанса нет. Инкремент ограничен γ < 4×10⁻¹¹. Причина — нелинейный сдвиг или большие начальные условия — проверяется следующим прогоном. Код уже скорректирован: плотная выборка, малые начальные, max_step = 0.3.

Часть II показывает: при астрофизических параметрах (ξ ~ 0.5–1.0, δ ~ 10⁻²–10⁻¹) та же модель даёт наблюдаемый порядок для:

  1. Глитчей Δν/ν ~ 10⁻⁶ — при ξΔα ~ 10⁻⁶, что физично.
  2. Непрерывных гравитационных волн на m = 2Ω_prec с h₀ ~ 10⁻²⁸–10⁻²⁷ — на уровне ET/CE, ниже порога LIGO (кроме PSR J0537-6910 при δ ~ 0.1, m ~ 12.4 Гц — нижний край LIGO).

Разделительный тест: классика (сверхтекучесть) даёт глитч, но не даёт пика на m. Метрика даёт оба и предсказывает корреляцию. Это проверяемо: искать в данных LIGO/ET не «пятно», а пик на m после глитча.

Что это меняет для цикла. Модель — кандидат на объяснение астрофизических аномалий. Не «объяснение» — кандидат. В лаборатории — устойчивый отрицательный результат, который ограничивает инкремент и ставит вопрос о причине. В астрофизике — наблюдательная программа. Одно без другого неполно.

Отрицательный результат — это не провал. Это граница, за которой начинается следующий слой.


КАРТА ДОСТОВЕРНОСТИ

УтверждениеСтатус
Код реализует уравнения корректноПроверено:E_drift ~ 10⁻¹⁰, отношение частот 2.000
При ξ = 0.1, I₁ ≠ I₂ резонанс развивается и насыщается (демо)Численно подтверждено
При ξ = 1×10⁻⁴, I₂ = 1.1, m = 2Ω_prec экспоненциального резонанса нет за 14.8% времениЧисленно подтверждено (устойчивый отрицательный результат)
Инкремент ограничен γ < 4×10⁻¹¹Численно подтверждено (из разброса между чанками ±1.2%)
Энергия сохраняется линейно, E_drift ~ 6.5×10⁻⁴ за 14.8%Численно подтверждено
α накачивается до ~8×10⁻³ (DC-уровень)Численно подтверждено
ω₁ растёт линейно+0.52%НЕ установлено — тренд меньше разброса, требует плотной выборки
Насыщение — физическое свойство моделиНЕ установлено, гипотеза, проверяется следующим прогоном
Причина отсутствия резонанса — большие начальные условияГипотеза, проверяется следующим прогоном
Причина отсутствия резонанса — нелинейный сдвиг частотыГипотеза, проверяется сканером
Уравнения масштабируются на 44 порядкаПроверено: безразмерная симуляция, E_drift ~ 10⁻¹⁰
Глитч Δν/ν ~ 10⁻⁶ соответствует Δα ~ 10⁻⁶–10⁻⁵ при ξ ~ 0.5–1.0Оценка, физична, требует проверки EOS
m = 2Ω_prec для Vela: 0.22–2.2 ГцОценка, зависит от δ
m = 2Ω_prec для PSR J0537-6910: 1.24–12.4 ГцОценка, зависит от δ
h₀ для Vela: 4.6×10⁻²⁸–4.6×10⁻²⁷Оценка, зависит от ε, ξ
h₀ для PSR J0537-6910: 8.2×10⁻²⁹–8.2×10⁻²⁸Оценка, r = 49.6 кпк
Vela и PSR J0537-6910 — цели для ET/CE, не LIGOОценка, следует из карты чувствительности
Классика не предсказывает CW на 2Ω_precУтверждение, требует проверки литературы
Корреляция «глитч → пик на m» — разделительный тестПредсказание, проверяемо в данных LIGO/ET
ξ ~ 1 для ядерной материиГипотеза, риск смены знака инерции при α ~ 1
Реальность поля αГипотеза, требует эксперимента

ПРИЛОЖЕНИЕ. ФОРМУЛЫ ЦИКЛА

A.1. Лагранжиан и уравнения движения

Лагранжиан: L = ½ Σ Iᵢ(α) ωᵢ² + ½ α̇² − ½ m² α²

Моменты инерции: Iᵢ(α) = Iᵢ⁰(1 + ξᵢα)

Связь ξ₃:ξ₃ = (ξ₁I₁⁰ + ξ₂I₂⁰) / I₃⁰

Уравнения Эйлера–Пуанкаре:Iᵢ(α) ω̇ᵢ + ξᵢIᵢ⁰ α̇ ωᵢ + (ω × I(α)ω)ᵢ = 0

Уравнение поля:α̈ + m²α = ½ Σ ξᵢ Iᵢ⁰ ωᵢ²

Энергия:E = ½ Σ Iᵢ(α) ωᵢ² + ½ α̇² + ½ m² α² = const


A.2. Частота прецессии и резонанс

Ω_prec² = (I₃−I₁)(I₃−I₂)ω₃² / (I₁I₂)

m ≈ 2Ω_prec (главный резонанс, n = 1)

Подпись резонанса:f_α / f_ω₁ = 2.000

Критическое поле инверсии:α_crit = (I₂⁰−I₁⁰) / (ξ(I₁⁰+I₂⁰))


A.3. Демо-режим и физичный режим (Часть I)

Демо (ξ = 0.1, ω₃ = 10, I₂ = 1.2):Ω_prec = 10, m = 20. Резонанс развивается, насыщается. E_drift ~ 10⁻⁸ за t = 10⁵.

Физичный (ξ = 1×10⁻⁴, ω₃ = 1, I₂ = 1.1):Ω_prec = 1, m = 2. Экспоненциального резонанса нет за 14.8% от t_end; инкремент ограничен γ < 4×10⁻¹¹. α накачивается до ~8×10⁻³ (DC-уровень). E_drift ~ 6.5×10⁻⁴ за 14.8%, линейный. Количественные утверждения о линейном росте ω₁ и плато не установлены (разреженная выборка).

Инкремент: γ ≈ ξ·α·Ω_prec / 2. Время развития:t_e ≈ 1/γ.

Скорректированный код: dt_eval = 1.0, ω₁(0)=0.01, ω₂(0)=0.005, max_step=0.3, n_pts_chunk = 10⁶.


A.4. Астрофизический полигон (Часть II)

δ = (I₃−I₁)/I₁

Ω_prec ≈ δ·ω₃, m = 2Ω_prec ≈ 2δ·ω₃

Связь глитча: Δω/ω ≈ −ΔI/I ≈ −ξ·Δα

Амплитуда: h₀ = (4G/c⁴) · (I₃·ε·ω² / r)

Эллиптичность: ε = (I₃−I₁)/I₃


A.5. Численные значения для пульсаров

ПульсарP (с)ω₃ (рад/с)I (кг·м²)rΔν/ν
Vela (B0833-45)0.08970.410³⁸287 пк~10⁻⁶ (до 3×10⁻⁶)
Crab (B0531+21)0.03319010³⁸2 кпк10⁻⁸–10⁻⁶
PSR J0537-69100.01639210³⁸49.6 кпк0.36×10⁻⁶


A.6. Карта чувствительности

Пульсарδ = 10⁻² → mδ = 10⁻¹ → mДетектор
Vela0.22 Гц2.2 ГцET/CE
PSR J0537-69101.24 Гц12.4 ГцET/CE (LIGO при δ ~ 0.1)
Пульсарh₀ (ε = 10⁻⁶)h₀ (ε = 10⁻⁵)Детектор
Vela4.6×10⁻²⁸4.6×10⁻²⁷ET
PSR J0537-69108.2×10⁻²⁹8.2×10⁻²⁸ET


A.7. Константы

  • G = 6.674×10⁻¹¹ м³/(кг·с²)
  • c = 2.998×10⁸ м/с
  • 4G/c⁴ = 3.3×10⁻⁷⁷ с²/(кг·м²)
  • 1 пк = 3.086×10¹⁶ м
  • 1 кпк = 3.086×10¹⁹ м


A.8. Следующие шаги (пункты 1–7)

1. Закрыть пик. Скан m вокруг 2Ω_prec, ±10%, 21 точка, ENV-рост. Пик/фон > 2 — сигнал.

2. Карта режимов. 2D скан (ξ, m), область, где ENV-рост > 5%.

3. Связь с формулой X. Проверить, зависит ли X ≈ (m_p/(m_e·α²))^(1/3) от ξ и m.

4. Нелинейности и диссипация. Заменить линейную связь на нелинейную, добавить γ·α̇. Проверить, выживает ли резонанс.

5. Тензорная метрика. Iᵢⱼ = Iᵢⱼ⁰ + hᵢⱼ, вывести Iᵢ(α) из геометрии. Следующий цикл.

6. Термодинамика. Добавить γ·α̇ и производство энтропии. Связь с LENR. После эксперимента.

7. Квантование. Уравнение Шрёдингера для α с потенциалом от кинетической энергии. «Метрическая спектроскопия». Долгосрочно.

Кристалл памяти обновлён: 01.11.2026. Проект «Рецепт хрупкого бессмертия». 


Рабочая модель "Чебурашка и Джанибеков". Обновлено 02.11.2026

Представьте себе в космосе летит огромная нейтронная суперзвезда — Чебурашка. И представим, что Джанибеков получил богатырскую силушку и решил, что гайка это слишком мелко, пора попробовать раскрутить звезду. И вот он берёт нейтронного Чебурашку и раскручивает его. Если бы Джанибеков ограничился бы гайкой, то в какой-то момент по гайке бы начала пробегать волна, и если бы мы покрасили её в синий, то увидели бы что атомы оказавшиеся внутри волны идущей по материи уходят в пространственный карман. Если бы сверхлёгкие чёрные дыры существовали, то эффект был бы аналогичным.

Но из-за этого же эффекта на нейтронной звезде возникает глитч. Глитч это «чих» Звезды.

Чебурашка такого, конечно, не ожидал, и потому, оказавшись в инометрическом иномирье (столкнувшись с новым для себя миром) звезда как детектор начинает неистово чихать, всё чаще переворачиваясь. Как мы знаем, из предыдущих статей: «Чих» Чебурашки, как детектора, означает, что он попал в область инометрического поля. А мы помним, что по нему бегут синие волны, которые погружают его атомы в инометрику.

С каждым новым чихом звезда союзмультфильма переводит энергию из вращения в переворачивание. И потому эффект переворачивания, а это и есть эффект Джанибекова, возникает все чаще. Но поскольку звезда обладает своим карманом — пузырём, остановить вращение Чебурашки с помощью трения о реальность, так как такой объект сам является источником реальности, не получается. Во всяком случае быстро. Поэтому он переворачивается все чаще и чаще.

Так вот возвращаясь к микрофону. Если перевести это учащение переворачиваний в звук, то частота была бы всё выше и выше.

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

Таким образом Чебурашка сам себя переворачивает всё сильнее с помощью чихания, и снижает скорость своего вращения со временем. Волны по нему бегут всё медленнее, а переворачивание идёт всё быстрее. И энергию, что получил от Джанибекова, звезда тратит на то, чтобы переворачиваться всё чаще.

Если в этот момент внутри звезды есть инометрические атомы, то они станут светиться всё ярче. Что означает, что нейтронный Чебурашка станет светиться со временем только ярче, пока не потратит весь их запас. Так бы всё и было, если бы мы не учли, что он раскручен.

Но опять же если вернуться к нашей системе дней вещества, то мы увидим, что «чем быстрее вращается объект, тем быстрее его атомы переходят в инометрику». Дело в том, что из-за изменения размера волн их стало больше на поверхности звезды. И они начинают наслаиваться друг на друга. Система начинает насыщаться и настаиваться. Процесс настаивания — это новая тонкая структура. Поскольку мы выразили её в виде волн. То со стороны выглядит, так словно вода поднялась, возникает неспокойный океан волн. Это и есть процесс насыщения. Сама же звезда наоборот сжимается в более малый размер под действием поля. Поэтому инометрические атомы, попав под влияние синих волн, получают комфортные для себя условия. Что позволяет им испускать свет без перехода в другое состояние. Без возврата в Т1, они так и остались в Т2, просто стали светиться чуть ярче. Если же атомов Т2 инометрики окажется в какой-то момент недостаточно, то к ним присоединятся атомы Т1. Так как попадут под влияние океана из многочисленных волн и перейдут в состояние инометрики Т2.

Поэтому Джанибеков просто раскрутил, а дальше процесс идет сам. И со временем Чебурашка светится всё ярче.

Подробнее о волнах и океане вы можете прочитать: Эффект Джанибекова: Восемь состояний, Оумуамуа и лестница вещества.

А найти чихающего чебурашку вы можете здесь:  Реальность как ракушка: витки, неваляшки и детектор инометрики

+13
57

0 комментариев, по

531 0 298
Мероприятия

Список действующих конкурсов, марафонов и игр, организованных пользователями Author.Today.

Хотите добавить сюда ещё одну ссылку? Напишите об этом администрации.

Полезное

Подборка полезных статей и записей в блогах.

Хотите добавить сюда ещё одну ссылку? Напишите об этом администрации.

Наверх Вниз