Нейтронные звёзды перед столкновением деформируются не так, как считалось: физики обнаружили три релятивистских механизма, которых не было в моделях

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

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

Деформация влияет на движение системы и на гравитационные волны, которые она излучает. В 2017 году детекторы LIGO и Virgo зарегистрировали сигнал GW170817 — первое наблюдение слияния двух нейтронных звёзд. Из этого сигнала удалось извлечь ограничения на приливную деформируемость, а через неё — на уравнение состояния. Это стало одним из главных результатов гравитационно-волновой астрономии.

Но надёжность извлечённых параметров целиком определяется качеством теоретической модели, заложенной в анализ. И именно здесь обнаруживается серьёзная проблема.

Почему существующие модели недостаточны

Стандартный подход описывает приливной отклик нейтронной звезды в так называемом адиабатическом приближении. Оно предполагает, что внешнее гравитационное поле меняется медленно — настолько, что звезда успевает в каждый момент прийти в равновесие с текущей конфигурацией приливного поля. Форма звезды определяется мгновенным значением приливных сил, без учёта истории.

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

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

Для нынешних детекторов эта ошибка, вероятно, остаётся в пределах погрешности. Но детекторы следующего поколения — Cosmic Explorer и Einstein Telescope — будут на порядок чувствительнее. Систематическая ошибка адиабатической модели станет доминирующей. Без корректного описания динамических приливов новые инструменты окажутся бесполезны для изучения внутренней структуры нейтронных звёзд.

Что было известно из ньютоновской теории

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

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

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

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

Почему обобщение на ОТО не давалось

Перенос модовой картины в ОТО блокировался тремя трудностями.

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

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

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

Предпринимавшиеся ранее попытки решить задачу в ОТО давали только числа — конкретные значения приливного отклика, полученные прямым численным интегрированием. Но они не давали модовой картины: не позволяли разложить отклик на вклады отдельных мод, идентифицировать резонансы, построить аналитическую формулу.

Что сделали авторы

Группа физиков из университетов Иллинойса, Калифорнии, Монтаны и Института фундаментальных исследований Тата — Абхишек Хегаде, Джеймс Квон, Теджасви Венумадхав, Хан Юй и Николас Юнес — нашли способ обойти все три трудности. Их работа, опубликованная в июле 2025 года, содержит два основных результата.

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

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

Второй результат: приливная сила в ОТО. Авторы предложили способ продолжить приливное поле внутрь звезды — несмотря на то, что однозначного разделения не существует. Среди бесконечного множества возможных вариантов они выделили класс продолжений, при которых приливная составляющая играет роль внешней силы, возбуждающей колебания. При таком выборе амплитуды мод подчиняются уравнению гармонического осциллятора с внешней силой — в точности той же структуры, что и в ньютоновской теории.

Однако сама сила оказывается устроена существенно сложнее. В ньютоновской теории она определяется одним механизмом: связью между приливным полем и распределением плотности. В ОТО действуют четыре различных механизма.

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

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

Итоговая формула и её устройство

Авторы получили явное выражение для функции приливного отклика нейтронной звезды — зависимости мультипольного момента (характеризующего деформацию) от частоты и амплитуды приливного поля. Это выражение представляет собой сумму по модам, где каждый член содержит резонансный знаменатель (вклад резко возрастает при приближении частоты приливного поля к собственной частоте моды) и два обобщённых интеграла перекрытия. В ньютоновском пределе оба интеграла совпадают и сводятся к единственному классическому выражению. В полной ОТО они различны и содержат все четыре описанных выше механизма.

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

Практические приложения

Формализм открывает несколько направлений для непосредственного применения.

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

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

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

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

Что это меняет

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

До этой работы у физиков было два варианта. Либо использовать модовое разложение, но только в ньютоновском приближении — и мириться с неконтролируемой ошибкой в релятивистском режиме. Либо решать задачу в полной ОТО численно — и отказаться от модовой картины, от аналитических формул, от возможности строить быстрые шаблоны.

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

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

Источник:arXiv