В поисках природы тёмной материи одна из самых интригующих гипотез уходит корнями в первые мгновения после Большого взрыва. Тогда из квантовых флуктуаций инфляционного поля рождались сгустки вещества, способные сколлапсировать в первичные чёрные дыры. Со времён Стивена Хокинга первичную плазму обычно считали идеальной жидкостью — текучей и лишённой внутреннего трения. Но реальность могла быть иной: вместо воды — густая патока, где каждое движение встречает сопротивление. Именно эту «вязкость» Вселенной и её роль в рождении чёрных дыр впервые удалось оценить в полноценном численном эксперименте.
Группа исследователей применила BSSN-формализм общей теории относительности вместе с уравнениями вязкой гидродинамики Эккарта, чтобы проследить судьбу сверхплотных областей в расширяющейся Вселенной. Результат нагляден: чем выше вязкость, тем труднее гравитации довести сгусток до схлопывания. Математически это выражается в эффективном уравнении состояния фоновой жидкости: \( w_{\text{eff}} = w - \epsilon \sqrt{24\pi} \). Параметр ε — мера «липкости» среды. Когда он растёт, давление падает, газ становится мягче, но одновременно возрастает внутреннее трение. Звуковые волны уже не могут легко стабилизировать возмущение, и для запуска коллапса требуется больший начальный перепад плотности. Сдвиг критической амплитуды µ_c подчиняется простой линейной зависимости: \(\mu_c - \mu_{c,0} = K \epsilon^\Gamma\), где Γ ≈ 1, а K порядка единицы. Это перекликается с классической неустойчивостью Джинса: диссипация замедляет сжатие, и гравитации приходится преодолевать дополнительное сопротивление — подобно тому, как Джеймс Джинс когда-то показал роль давления в противодействии сжатию.
Помимо роста порога, вязкость меняет и массу рождающихся чёрных дыр. Вблизи критического коллапса масса подчиняется скейлинговому закону \( M_{\text{PBH}} / M_H = k |\mu - \mu_c|^\gamma \). Вязкость увеличивает амплитуду k (например, с 8.30 до 10.22 для одного из профилей возмущений), слабо влияя на показатель γ. Иными словами, гравитационное сжатие в патоке рождает более массивные дыры при тех же относительных отклонениях от порога. Функция масс смещается в сторону больших значений, что напрямую затрагивает интерпретацию данных гравитационного микролинзирования и стохастического фона гравитационных волн.
Полученный результат — не просто уточнение. Он указывает, что для надёжных предсказаний обилия первичных чёрных дыр как кандидатов в тёмную материю необходимо учитывать неидеальность первичной плазмы. Более того, вязкость смягчает уравнение состояния фона, что, возможно, проливает свет на проблему хаббловского натяжения — загадочного расхождения в измерениях скорости расширения Вселенной. В перспективе моделирование с более реалистичной вязкостью второго порядка (теория Израэля–Стюарта) и включение анизотропии позволят вычислять спины и эллиптичности рождающихся чёрных дыр, прокладывая мост между микрофизикой коллапса и крупномасштабной космологией. Быть может, сама «липкость» ранней Вселенной окажется ключом к разгадке тёмной материи и судьбы космологического расширения.
Вязкость обычно ассоциируется с чем-то замедляющим, но здесь она играет роль стабилизатора: «вязкий мёд» космической среды требует большего избытка плотности, чтобы схлопнуться. Лишь в 2026 году численное моделирование позволило количественно оценить этот эффект — спустя полвека после первых работ Хокинга и Карра. Это означает, что даже неуловимые микроскопические свойства первичной плазмы могли определить, сколько тёмной материи мы видим сегодня.