Что на самом деле внутри нейтронной звезды? Как новая математическая модель раскрывает состав самых плотных объектов Вселенной
Гравитационно-волновая астрономия открыла человечеству совершенно новый способ изучения Вселенной. Фиксируя мельчайшие искажения пространства-времени, исходящие от слияний сверхплотных объектов, физики получают доступ к информации, которую невозможно добыть с помощью традиционных оптических или радиотелескопов. Одно из главных направлений в этой области — изучение двойных нейтронных звезд. Анализ их слияний дает уникальную возможность заглянуть внутрь этих объектов и понять, как ведет себя материя при плотностях, превышающих плотность атомного ядра.
Однако, чтобы точно расшифровать гравитационно-волновой сигнал и извлечь из него данные о внутреннем строении звезды, необходимо досконально понимать, как именно две звезды деформируют друг друга в процессе сближения. До недавнего времени математические модели, описывающие этот процесс, опирались на упрощенные классические теории. Новая работа группы физиков-теоретиков из Университета Иллинойса и Калифорнийского университета предлагает строгое аналитическое решение этой задачи в рамках полной Общей теории относительности, избавляя астрофизику от накопившихся неточностей.
Анатомия приливной деформации
Когда две нейтронные звезды вращаются вокруг общего центра масс, их гравитационные поля непрерывно воздействуют друг на друга. По мере сокращения орбиты это воздействие усиливается, что приводит к вытягиванию звезд вдоль оси, соединяющей их центры. Физики называют это явление приливной деформацией.
Степень этой деформации критически важна для науки. Она напрямую зависит от «уравнения состояния» нейтронной звезды — математической зависимости между давлением и плотностью в ее недрах. Если материя внутри звезды оказывает сильное сопротивление сжатию, звезда деформируется слабо. Если материя более податлива, деформация будет значительной. В астрофизике эта характеристика описывается параметрами, известными как числа Лява.
Приливная деформация оставляет четкий отпечаток на гравитационной волне, излучаемой системой. Энергия, которая могла бы пойти на поддержание скорости орбитального вращения, тратится на деформацию самой звезды и возбуждение внутренних движений ее материи. В результате орбита сокращается быстрее, и частота гравитационно-волнового сигнала меняется иначе, чем если бы звезды были абсолютно жесткими сферами.
Сложность заключается в том, что этот процесс не является статичным. В последние миллисекунды перед столкновением (в фазе так называемого позднего сближения) гравитационное поле партнера меняется настолько быстро, что звезда не успевает прийти в равновесие. Возникает динамический приливной отклик, который сопровождается возбуждением собственных резонансных частот внутри звезды.
Пределы ньютоновской физики и сложность ОТО
Исторически для описания приливных сил применялась ньютоновская физика. В рамках классической механики гравитация рассматривается как сила, действующая на массу, а деформацию звезды можно легко описать через возмущения плотности ее вещества. Чтобы учесть релятивистские эффекты — такие как влияние огромных скоростей и сильной гравитации — астрофизики используют пост-ньютоновский формализм. Этот метод добавляет к классическим уравнениям серию поправочных коэффициентов.
Пост-ньютоновский подход отлично описывает систему, когда звезды находятся на значительном расстоянии друг от друга. Но когда они сближаются вплотную, этот метод теряет свою предсказательную силу. Гравитационные поля становятся настолько экстремальными, что приближенные вычисления начинают расходиться с реальностью.
Казалось бы, логичным шагом должно стать использование Общей теории относительности в ее полном виде, без упрощений. Однако уравнения Эйнштейна представляют собой систему нелинейных дифференциальных уравнений в частных производных. В ОТО гравитация — это не сила, а искривление самого пространства-времени. В двойной системе пространство искривлено массами обеих звезд, и в математическом аппарате теории крайне сложно отделить гравитационное поле самой звезды от внешнего приливного поля, создаваемого ее партнером.
Прямое численное моделирование, при котором уравнения Эйнштейна решаются на суперкомпьютерах, позволяет получить конечный результат — формулу изменения сигнала. Но численные методы крайне ресурсоемки и, что более важно, они не дают физического понимания процесса. Численная симуляция не показывает, какие именно внутренние механизмы звезды привели к наблюдаемой деформации, и не позволяет разложить сложный процесс на базовые составляющие.
Деконструкция гравитационного отклика
Авторы нового исследования поставили перед собой цель создать физически прозрачную, аналитическую модель динамических приливов, которая бы полностью подчинялась законам Общей теории относительности, но при этом обладала математической структурой классической механики. Им удалось доказать, что сложнейшая динамика искривленного пространства-времени внутри деформируемой нейтронной звезды может быть сведена к набору независимых уравнений вынужденного гармонического осциллятора.
Для решения проблемы нелинейности уравнений Эйнштейна исследователи применили метод сращиваемых асимптотических разложений. Они разбили пространство вокруг двойной системы на несколько вложенных зон:
- Внутренняя зона: недра звезды и пространство в непосредственной близости от нее, где доминирует ее собственная сильная гравитация. Здесь уравнения Эйнштейна-Эйлера должны решаться в полной мере.
- Пост-ньютоновская зона: область на значительном удалении от звезды, где гравитационное поле партнера можно рассматривать как внешнее возмущение, а эффекты ОТО относительно слабы.
- Буферная зона: промежуточная область, в которой пересекаются границы применимости обеих моделей.
Ключевой математический прорыв работы заключается в том, что авторы смогли аналитически сшить решения уравнений на границе буферной зоны. Они доказали, что если сильное гравитационное поле внутри звезды корректно состыковать с внешним пост-ньютоновским полем в буферной зоне, то оператор, описывающий возмущения системы, становится самосопряженным.
В математической физике свойство самосопряженности гарантирует, что сколь бы сложной ни была система, ее реакцию на внешнее воздействие можно разложить на полную систему ортогональных мод. Говоря строгим языком, деформацию пространства-времени и материи нейтронной звезды можно представить как сумму независимых колебаний. Амплитуда каждого такого колебания вычисляется через уравнение, формально идентичное уравнению осциллятора, на который действует внешняя периодическая сила.
Релятивистский интеграл перекрытия
Следующим шагом стало определение природы этой внешней силы в контексте ОТО. В ньютоновской физике интенсивность возбуждения внутренних колебаний звезды определяется так называемым интегралом перекрытия. Он вычисляется на основе того, как внешнее гравитационное поле взаимодействует с распределением плотности массы внутри объекта.
В Общей теории относительности масса является не единственным источником гравитации. Согласно тензору энергии-импульса, давление жидкости, внутреннее напряжение материи и само искривление метрики пространства вносят вклад в гравитационное взаимодействие. Исследователи вывели принципиально новую формулу — релятивистский интеграл перекрытия.
Эта формула математически описывает, как сильное гравитационное поле нейтронной звезды влияет на ее же собственную приливную деформацию. Обобщенный интеграл состоит из нескольких компонентов:
- Связь внешнего приливного поля с градиентами энергии внутри звезды.
- Связь приливного поля с градиентами внутреннего давления.
- Взаимодействие приливного возмущения с векторным потенциалом гравитационного поля самой звезды (эффекты увлечения инерциальных систем отсчета).
Так, авторы создали полноценный математический аппарат. Он позволяет вычислять функцию приливного отклика нейтронной звезды через сумму мод колебаний, учитывая все нелинейные релятивистские поправки, которые классическая физика вынуждена игнорировать.
Значение для фундаментальной науки
Предложенная аналитическая модель появляется в критически важный для астрофизики момент. Существующие детекторы гравитационных волн (LIGO, Virgo, KAGRA) уже позволили зафиксировать слияния нейтронных звезд, но их чувствительность ограничена. В следующем десятилетии в строй вступят детекторы третьего поколения, такие как европейский Einstein Telescope и американский Cosmic Explorer. Их точность будет на порядки выше: они смогут улавливать гравитационные волны от слияний с огромным отношением сигнала к шуму.
При анализе данных такого качества главная проблема смещается от инструментальных ограничений к теоретическим. Если для расшифровки высокоточного сигнала от детекторов третьего поколения использовать устаревшие пост-ньютоновские модели, возникнет систематическая ошибка. Алгоритмы будут находить наилучшее совпадение между реальным сигналом и теоретическим шаблоном, но из-за несовершенства шаблона физики получат неверные значения массы, радиуса и жесткости нейтронной звезды. Это сделает невозможным точное определение уравнения состояния ядерной материи.
Новый формализм решает эту проблему превентивно. Сведение сложной динамики ОТО к набору независимых модальных уравнений позволяет создавать теоретические шаблоны гравитационных волн, которые являются одновременно вычислительно эффективными и релятивистски точными. Модель позволяет быстро и строго рассчитывать моменты резонансов, когда орбитальная частота двойной системы совпадает с собственными частотами колебаний звезды.
Работа исследователей доказывает, что даже в самых экстремальных режимах, где гравитация разрушает привычные физические концепции, система сохраняет внутреннюю математическую структуру, поддающуюся строгому анализу. Переведя язык тензорного исчисления Эйнштейна в плоскость модального анализа, физики получили надежный инструмент для расшифровки внутренней структуры объектов, состоящих из самой плотной материи во Вселенной.
Источник:arXiv