Калькулятор столкновения чёрных дыр
Рассчитайте время слияния, энергию гравитационных волн и финальные свойства системы из двух чёрных дыр.
Введите массы двух чёрных дыр, их начальное орбитальное расстояние, эксцентриситет и угол наклона, чтобы вычислить параметры спирального сближения и слияния по формулам общей теории относительности Питерса.
Калькулятор столкновения чёрных дыр
Рассчитайте время слияния, энергию гравитационных волн и финальные свойства системы из двух чёрных дыр.
О калькуляторе столкновения чёрных дыр
Когда две чёрные дыры образуют гравитационно связанную бинарную систему, они постепенно сближаются, теряя энергию и угловой момент за счёт излучения гравитационных волн. Этот процесс, называемый inspiral, следует математической схеме, впервые выведенной Питером Питерсом в 1964 году на основе общей теории относительности Эйнштейна. Калькулятор столкновения чёрных дыр реализует формулу Питерса, чтобы оценить, сколько времени займёт слияние системы, сколько энергии она излучит и как будет выглядеть финальный объект.
Самая важная производная величина в астрономии гравитационных волн — chirp-масса, определяемая как M_chirp = (m₁ m₂)^(3/5) / (m₁+m₂)^(1/5). Chirp-масса задаёт скорость роста частоты гравитационной волны — характерный «свип» от низких частот к высоким, благодаря которому сигнал и получил своё название. Вместе с симметричным отношением масс η = m₁m₂/(m₁+m₂)² chirp-масса содержит всю информацию, необходимую для расчёта динамики inspiral в ведущем порядке.
Формула Питерса (1964) для круговой орбиты имеет вид T = (5/256) × c⁵ a⁴ / (G³ m₁ m₂ M_total), где a — начальная большая полуось, G — гравитационная постоянная, а c — скорость света. Для эксцентричной орбиты время слияния уменьшается примерно на множитель (1−e²)^(7/2), что означает: чем выше эксцентриситет, тем эффективнее система излучает энергию и тем быстрее сливается по сравнению с круговой орбитой при том же начальном расстоянии. Это приближение становится всё точнее при e ≲ 0.6.
Энергия гравитационных волн, излучённая во время inspiral, оценивается примерно в 5% от энергии приведённой массы (μc²), что согласуется с численнорелятивистскими моделями бинарных систем сравнимых масс. Оставшаяся масса образует финальную чёрную дыру, радиус Шварцшильда которой равен r_s = 2 G M_final / c². Пиковая частота гравитационных волн на самой внутренней устойчивой круговой орбите (ISCO) равна f_peak = c³ / (π × 6√6 × G × M_total) и обозначает переход от inspiral к plunge и ringdown — самый громкий и наиболее заметный момент слияния.
Первое прямое обнаружение гравитационных волн, GW150914, было сделано LIGO 14 сентября 2015 года. Оно возникло при слиянии двух чёрных дыр массой примерно 36 и 29 солнечных масс на расстоянии 1,3 миллиарда световых лет. Событие излучило около 3 солнечных масс энергии в виде гравитационных волн за долю секунды, на мгновение превзойдя по гравитационной светимости всю наблюдаемую Вселенную. С тех пор коллаборация LIGO–Virgo–KAGRA обнаружила более 90 событий бинарных слияний, превратив астрономию гравитационных волн в точную науку.
Примеры столкновения чёрных дыр
В таблице ниже показаны типичные системы из двух чёрных дыр и их ключевые параметры слияния.
| Параметры системы | Ключевые результаты | Событие / контекст |
|---|---|---|
| m₁=36 M☉, m₂=29 M☉, a=10,000,000 km, e=0 | T_merge ≈ 94.4 Myr, M_chirp ≈ 28.1 M☉, f_peak ≈ 67.6 Hz | Похоже на GW150914 (LIGO, 2015) |
| m₁=1000 M☉, m₂=800 M☉, a=100,000,000 km, e=0.3 | T_merge ≈ 32.0 Myr, M_chirp ≈ 778 M☉, f_peak ≈ 2.44 Hz | Бинарная система чёрных дыр промежуточной массы |
| m₁=20 M☉, m₂=20 M☉, a=5,000,000 km, e=0 | T_merge ≈ 25.0 Myr, M_chirp ≈ 17.4 M☉, f_peak ≈ 110 Hz | Звёздная бинарная система равных масс |
Как использовать калькулятор столкновения чёрных дыр
- Введите массу каждой чёрной дыры в солнечных массах (M☉). Звёздные чёрные дыры обычно имеют массу от ~3 до ~100 M☉; сверхмассивные чёрные дыры могут превышать 10⁹ M☉.
- Введите начальное орбитальное расстояние в километрах. Это начальная большая полуось бинарной орбиты.
- Задайте орбитальный эксцентриситет от 0 (круговая орбита) до 0.99 (почти радиальная). Большинство событий, зарегистрированных LIGO, входят в полосу детектора уже почти на круговых орбитах.
- Введите угол наклона в градусах (0° = face-on, 90° = edge-on). Он влияет на амплитуду гравитационных волн, наблюдаемую с Земли, но не на время слияния.
- Нажмите Рассчитать, чтобы увидеть chirp-массу, оценочное время слияния, энергию гравитационных волн, финальную массу чёрной дыры, радиус Шварцшильда и пиковую частоту ГВ.
Часто задаваемые вопросы
Что такое chirp-масса и почему она важна?
Chirp-масса M_chirp = (m₁m₂)^(3/5)/(m₁+m₂)^(1/5) — самый важный отдельный параметр для обнаружения гравитационных волн. Она задаёт скорость роста частоты гравитационной волны (chirp rate), поэтому астрономы могут очень точно измерять chirp-массу по форме сигнала ещё до того, как известны отдельные массы.
Насколько точна оценка времени слияния?
Формула Питерса, используемая здесь, точна для ранней фазы inspiral, когда расстояние намного больше радиуса Шварцшильда. Поправка на эксцентриситет (1−e²)^(7/2) — это приближение, которое хорошо работает при e ≲ 0.6. Для сильно эксцентричных орбит или очень компактных систем для точных оценок нужна численная относительность.
Почему более высокий эксцентриситет ускоряет слияние?
В точке наибольшего сближения (перицентре) эксцентричной орбиты тела движутся быстрее всего и находятся ближе всего друг к другу, что резко увеличивает мощность излучаемых гравитационных волн. За один оборот излучается больше энергии, орбита быстрее сжимается, и время слияния становится меньше по сравнению с круговой орбитой при том же среднем расстоянии.
Что такое ISCO и почему она задаёт пиковую частоту ГВ?
Самая внутренняя устойчивая круговая орбита (ISCO) — это граница, внутри которой вокруг чёрной дыры Шварцшильда (без вращения) не существует устойчивых круговых орбит. Когда inspiral достигает этой точки, меньший объект быстро проваливается внутрь. Орбитальная частота на ISCO, удвоенная для получения частоты гравитационной волны, соответствует самой высокой частоте inspiral-сигнала и началу ringdown.
Сколько энергии излучается в виде гравитационных волн?
Для слияний чёрных дыр сравнимых масс численные модели показывают, что примерно 4–8% общей массы излучается в виде гравитационных волн. Калькулятор использует около 5% энергии приведённой массы как грубую оценку. Для GW150914 около 3 солнечных масс (≈5% от общей массы) были преобразованы в энергию гравитационных волн за долю секунды.
Можно ли использовать этот калькулятор для слияний нейтронных звёзд?
Формула inspiral одинаково применима к двойным системам нейтронная звезда–нейтронная звезда (BNS) и нейтронная звезда–чёрная дыра (NSBH). Однако для BNS-ивентов приливное разрушение и уравнение состояния нейтронной звезды вносят поправки, которые здесь не учитываются. Для оценки порядка величины его можно использовать, но для точных результатов BNS нужны постньютоновские или численнорелятивистские модели формы сигнала.