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

Почему физики пересматривают теорию слияния нейтронных звезд: новый инструмент для гравитационно-волновой астрономии

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

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

Почему старые модели ошибаются

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

Адиабатическая модель не учитывает ни резонансов, ни накопления энергии в колебаниях. Результат — систематическая ошибка. Для нынешних детекторов погрешность терпима. Но новые — Cosmic Explorer и Einstein Telescope — будут в десять раз чувствительнее. Если не исправить модель, они просто не смогут изучить внутреннюю структуру звезд.

Как работает динамический прилив: пошагово
1. Две нейтронные звезды сближаются, орбитальная частота растет.
2. Приливные силы меняются все быстрее.
3. Внутри звезды возбуждаются моды — типы пульсаций со своими частотами.
4. Когда частота прилива совпадает с частотой моды — резонанс, амплитуда резко растет.
5. Звезда деформируется несинхронно с приливом — возникает динамический отклик.
6. Этот отклик оставляет след в гравитационно-волновом сигнале — фазовый сдвиг.

Что сделали авторы: новая математика для ОТО

Группа физиков из университетов Иллинойса, Калифорнии, Монтаны и Института Тата (работа 2025 года) решила задачу, которая не поддавалась двадцать лет. В ньютоновской гравитации динамические приливы разлагали по модам — как сумму колебаний маятников. В общей теории относительности (ОТО) это не работало из-за трех препятствий: граничные условия для изолированной звезды не дают полного набора мод; в ОТО нельзя четко отделить собственное гравитационное поле от внешнего; само поле — динамическая переменная, его нельзя исключить.

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

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

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

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

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

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

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

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

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