Гравитационные волны (ГВ), порождаемые ускоренно движущимися массивными объектами, в теории общей относительности описываются как возмущения метрики пространства-времени. В простейшей приближенной модели рассматривается линейная аппроксимация: волновое уравнение для малых возмущений метрики hμν на фоне плоского пространства:
□hμν = 0,
где □ — д’Аламбера оператор. В реальных астрофизических системах эта линейная модель оказывается недостаточной, поскольку источники гравитационных волн могут генерировать значительные нелинейные эффекты.
Нелинейность проявляется в нескольких аспектах:
Полные уравнения Эйнштейна:
$$ G_{\mu\nu} = \frac{8\pi G}{c^4} T_{\mu\nu}, $$
являются сильно нелинейными по метрике gμν. Для выделения волн вводят разложение:
gμν = ημν + hμν,
где ημν — фон плоского пространства, а hμν — малое возмущение. Вторая итерация разложения приводит к уравнению:
□hμν = Λμν(h, h) + ????(h3),
где Λμν — набор квадратичных по hμν членов. Эти квадратичные и кубические члены отвечают за самовоздействие волн, их нелинейное взаимодействие и генерацию вторичных гармоник.
Нелинейные эффекты оказывают прямое влияние на форму сигнала, который регистрируют детекторы типа LIGO и Virgo:
Для систем с малыми скоростями v ≪ c и слабым полем hμν ≪ 1 используют пост-Ньютоновские разложения. Амплитуда и фаза сигнала разворачиваются в ряд по параметру (v/c):
$$ h_{\mu\nu} = h_{\mu\nu}^{(0)} + \left(\frac{v}{c}\right)^2 h_{\mu\nu}^{(1)} + \left(\frac{v}{c}\right)^4 h_{\mu\nu}^{(2)} + \dots $$
Каждое последующее слагаемое учитывает все более сложные нелинейные эффекты, включая реакцию системы на излучение гравитационных волн.
Метод EOB заменяет сложное двухтельное взаимодействие на эффективное однокомпонентное движение в потенциальном поле, включающем нелинейные поправки. Это позволяет точно моделировать фазу и амплитуду до слияния и в момент него.
В тех случаях, когда амплитуды и скорости велики (v ∼ c), линейные и PN-методы теряют точность. Решение полных уравнений Эйнштейна на компьютере — единственный способ учесть все нелинейные эффекты. Используются методы разложения в 3+1 размерности, сеточные методы и адаптивные схемы времени и пространства.
Для анализа сигналов, зарегистрированных детекторами, часто объединяют PN-расчеты на ранних стадиях с нумерическими моделями на поздних стадиях слияния. Результатом является гибридный шаблон сигнала, учитывающий все важные нелинейные поправки.
Эффективный учет нелинейностей позволяет не только повысить точность анализа гравитационных волн, но и расширить возможности исследования физики слияния черных дыр, нейтронных звезд и других экстремальных астрофизических систем.