Анализ распределения эффективного спина (χeff) двойных черных дыр в каталоге GWTC-5.0 с помощью иерархического байесовского вывода выявил переход при массе главной черной дыры. Популяция моделируется смесью двух спиновых компонент с граничной массой, извлеченной из данных: 15,2(+4,3/−3,6) M⊙. Анализ искусственных каталогов показал, что такой переход вряд ли вызван флуктуациями конечной выборки. Ниже переходной массы распределение χeff узкое, с пиком при малом положительном значении, но со значительной поддержкой отрицательных спинов; выше — распределение шире, пик смещен к χeff≈0, а поддержка положительных и отрицательных значений примерно одинакова. Эти результаты указывают на статистически различные популяции: доминирующая группа с массами ~10 M⊙, вероятно, формируется иначе. Низкомассовая популяция согласуется с образованием из массивных кратных звезд в поле — либо через эволюцию изолированных двойных с большими начальными толчками черных дыр (≳100 км/с) ниже граничной массы, либо через динамическую эволюцию иерархических тройных систем; стандартные модели изолированной двойной эволюции не воспроизводят поддержку отрицательных χeff.
Слияния бинарных чёрных дыр, предсказанных ещё Карлом Шварцшильдом, — главный источник гравитационных волн, регистрируемых детекторами, такими как LIGO. Происхождение этих систем остаётся одной из центральных проблем астрофизики. Эффективный спин χeff — ключевая наблюдаемая, несущая информацию о массах, вращении и ориентации осей компонентов. Разные каналы формирования предсказывают различный характер спинов — от упорядоченных при изолированной эволюции до хаотических в плотных скоплениях. Понимание этих различий необходимо для реконструкции эволюции массивных звёзд и динамики звёздных популяций в галактиках.
В работе использованы данные 259 достоверных событий из каталога GWTC-5.0, зарегистрированных с помощью гравитационно-волновых обсерваторий. Популяционный анализ проводился методом иерархического байесовского вывода, моделируя распределение χeff как смесь двух компонент, разделённых массой перехода ˜m, которая выводилась непосредственно из данных. Компоненты описывались непараметрически с помощью гауссовских процессов, что позволяло избежать априорных предположений о форме. Для проверки устойчивости применялась также параметрическая модель. Чтобы исключить артефакты, было проведено моделирование 100 искусственных каталогов с плавным распределением спина, не зависящим от массы.
Обнаружен выраженный переход при массе ˜m = 15,2+4,3–3,6 M⊙. Ниже этого значения распределение χeff узкое (ширина около 0,05) с пиком около +0,03, но с заметной долей отрицательных значений. Выше перехода распределение становится значительно шире (∼0,12), пик смещается к нулю, а положительные и отрицательные спины встречаются почти одинаково часто. Различие между двумя популяциями статистически значимо на уровне >99% для медиан. Моделирование показало, что такой переход не может возникнуть из-за случайных флуктуаций при независимом от массы распределении спина. Важно, что результат не зависит от учёта отдельных необычных событий, таких как GW241011 и GW241110.
Результаты показывают, что доминирующая популяция слияний с массами около 10 M⊙ статистически отлична от остальных. Это ставит под сомнение универсальность сценариев изолированной двойной эволюции, которые в стандартном виде не способны воспроизвести долю отрицательных χeff. Для объяснения требуются либо большие натальные «пинки» (толчки) при коллапсе в чёрные дыры (≳100 км/с), что связано с физикой сверхновых, либо динамические эффекты в тройных системах. Таким образом, происхождение низкомассивных чёрных дыр оказывается тесно связанным с финальными стадиями эволюции массивных звёзд, где важную роль играют водород и гелий в оболочках, а также синтез элементов вплоть до углерода.
В будущем, с накоплением данных от текущих и планируемых детекторов (LIGO, Virgo, KAGRA, а также третьего поколения, таких как Einstein Telescope и Cosmic Explorer), появится возможность уточнить тонкую структуру перехода и изучить зависимость от красного смещения. Улучшенное моделирование коллапса звёзд и гидродинамики сверхновых поможет прояснить механизмы рождения чёрных дыр с большими импульсами. Также ожидается прогресс в прямых численных симуляциях тройных и кратных систем, что позволит сопоставить предсказания с данными.
Работа затрагивает широкий круг областей: от физики гравитационных волн и астрофизики компактных объектов до теории эволюции массивных звёзд и динамической астрономии. Она также стимулирует исследования в области ядерной астрофизики, где критически важны процессы при коллапсе ядра.
Дальнейшие шаги включают детальное сравнение предсказаний популяционного синтеза тройных систем с расширенной выборкой событий, а также поиск дополнительных маркеров, таких как распределение массовых отношений и эксцентриситета. Планируется более тщательное исследование зависимости от металличности и среды.
Обнаруженная спиновая граница напрямую связана с нерешёнными вопросами: каковы типичные скорости натальных толчков чёрных дыр? Как часто массивные звёзды образуются в тройных и более сложных системах? Насколько эффективно приливное взаимодействие синхронизирует вращение в тесных парах? Ответы на эти вопросы заполнят пробелы в нашем понимании финальных стадий звёздной эволюции и природы популяции гравитационно-волновых источников.
🎯 Масса Солнца составляет около 2×10³⁰ кг, а типичная чёрная дыра с массой 15 M⊙ имеет радиус горизонта событий всего около 45 километров — меньше, чем расстояние от Москвы до Серпухова.
🎬 В романе «Контакт» Карла Сагана и одноимённом фильме сигнал из системы Веги содержал зашифрованные гравитационно-волновые паттерны. Сегодня мы сами расшифровываем «сигналы» из глубин Вселенной, чтобы понять, как устроены фабрики чёрных дыр.