Килоновые, возникающие при слияниях нейтронных звёзд, производят тяжёлые элементы в ходе r-процесса и демонстрируют позднее инфракрасное излучение при температуре ниже 1000 К, которое невозможно объяснить атомными переходами. На примере событий AT2017gfo (GW170817) и AT2023vfi (GRB 230307A) показано, что в выбросах килоновых создаются условия для образования пыли из тугоплавких r-элементов — циркония, вольфрама и осмия. Кинетические расчёты на основе коэффициентов скоростей химических реакций для вольфрама как типичного представителя демонстрируют эффективную конденсацию пыли, особенно в медленных выбросах, что противоречит более ранним работам, основанным на классической теории нуклеации. Радиационно-гидродинамическое моделирование с учётом пыли успешно воспроизводит наблюдаемое инфракрасное излучение. Формирование и обилие пыли из r-процессных элементов сильно зависят от массы, состава и скорости выброса. Инфракрасное излучение такой пыли открывает новый метод диагностики синтеза тяжёлых элементов в слияниях нейтронных звёзд.
Слияния нейтронных звёзд, зарегистрированные по гравитационным волнам (событие GW170817), сопровождаются килоновыми — оптическими и инфракрасными транзиентами, питаемыми радиоактивным распадом свежесинтезированных элементов. Спектроскопия таких событий, особенно с помощью Джеймса Уэбба, позволила оценить массу выброшенного вещества и идентифицировать отдельные тяжёлые элементы. Однако поздние наблюдения обнаружили сильное инфракрасное излучение с температурой ниже 1000 K, которое не удаётся воспроизвести моделями, основанными только на атомных линиях. Ещё в середине XX века Фред Хойл и Маргарет Бербидж показали, что тяжёлые элементы синтезируются в звёздах, а сегодня детектирование гравитационных волн установками LIGO, созданными при ключевом участии Райнера Вайсса, открыло эру нейтронно-звёздной астрономии. Спектры килоновых характеризуются псевдоконтинуумом из-за доплеровского уширения множества линий.
Для проверки возможности пылеобразования исследователи рассчитали кинетическую эволюцию кластеров из тугоплавких элементов (используя вольфрам как типичный пример) в расширяющейся плазме выбросов. В отличие от классической теории нуклеации, учитывались реакции ассоциации и фрагментации для кластеров любых размеров, а также разрушение димеров быстрыми электронами от бета-распада. Затем с помощью численного моделирования переноса излучения в коде Sedona были синтезированы спектры запылённой килоновой, где микроскопические свойства пылинок (в рэлеевском пределе) определяли непрозрачность, а их температура рассчитывалась из баланса нагрев–охлаждение с учётом планковского закона излучения. Скорость роста кластеров задавалась молекулярно-динамическими коэффициентами для вольфрама, а образование устойчивых зёрен запускалось при падении температуры ниже ~1800 K.
Оказалось, что при температурах ниже ~1800 K (примерно через 10 дней после слияния) формирование пыли идёт очень эффективно в плотных, медленно расширяющихся слоях (со скоростями ≲0.1c). Масса пыли порядка 10^{-3} M⊙ достаточна для того, чтобы создать оптическую толщу около единицы в инфракрасном диапазоне. Внешние, быстрые слои (≳0.2c) остаются практически свободными от пыли из-за низкой плотности. Модельные спектры показывают почти чернотельный континуум с пиком на 5 мкм, что хорошо соответствует данным спектроскопии для AT2023vfi. Небольшой избыток излучения около 2.1 мкм объясняется вкладом линий от свободных от пыли внешних областей. Для сравнения, модель без пыли даёт чисто линейчатый, нетепловой спектр, совершенно непохожий на наблюдения. Таким образом, космическая пыль из тяжёлых элементов естественно решает загадку позднего ИК-излучения. Ключевой вывод: эффективность пылеобразования резко падает с ростом скорости выброса, что совпадает с наблюдаемым порогом ~0.1c для AT2017gfo и AT2023vfi.
Это открытие превращает инфракрасное свечение килоновых в чувствительный индикатор состава, массы и динамики выбросов. Эффективность пылеобразования сильно зависит от обилия тугоплавких элементов, которое, в свою очередь, определяется путями нуклеосинтеза. Так, для состава, похожего на звезду HD 222925 с усиленными третьим пиком r-процесса, пыли образуется много, а для состава, характерного для HD 122563 с ослабленными тяжёлыми элементами, — практически нет. Следовательно, наблюдаемое запыление AT2023vfi свидетельствует о богатом выходе элементов третьего пика, что согласуется с предсказаниями моделей слияний нейтронных звёзд. Наиболее вероятные источники медленных выбросов — аккреционные диски и их ветры, в которых и происходит основное пылеобразование.
Будущие наблюдения, особенно в среднем ИК-диапазоне за 5 мкм, позволят проследить переход от спектра с доминированием линий к пылевому континууму. Подобно тому, как в сверхновых II типа ловят момент конденсации пыли, временна́я спектроскопия килоновых в первые 10–20 дней может напрямую показать рождение тугоплавких пылинок. Кроме того, разнообразие условий в разных слияниях (масса, скорость, состав) должно порождать всё многообразие поздних спектров — от чисто линейчатых до полностью пылевых.
Результаты затрагивают не только астрофизику транзиентов, но и галактическую химическую эволюцию: образование пыли в килоновых может быть одним из источников межзвёздных тугоплавких пылинок в древних галактиках. Это также даёт новый метод оценки темпа производства тяжёлых элементов во Вселенной.
Важными шагами станут учёт более реалистичных смесей элементов (не только вольфрама) и включение процессов разрушения пыли β-частицами и ударными волнами в остатках килоновых. Также необходимо детальное неравновесное моделирование ионизации и нагрева газа в пылевой фазе.
Работа напрямую касается фундаментальной проблемы происхождения элементов тяжелее железа. Если килоновые действительно являются главными «фабриками» золота, платины и урана, то обнаружение и анализ их запылённых выбросов предоставляют уникальный инструмент для проверки теорий нуклеосинтеза в экстремальных условиях. Кроме того, образование пыли связывает ядерную астрофизику с физикой межзвёздной среды и процессами формирования планетных систем.
🎯 В выбросах килоновой может синтезироваться золота в сотни раз больше массы Земли, и часть его конденсируется в настоящие космические пылинки, которые некоторое время светятся, как крошечные раскалённые угольки.