Модель calypso: новый инструмент для анализа аккреции вокруг двойных звезд
Астрофизики представили стохастическую эмуляторную модель calypso, предназначенную для изучения аккреционных потоков вокруг двойных звезд. Работа опубликована на сервере препринтов arXiv. Модель позволяет воспроизводить временные ряды аккреции с учетом внутренней случайности процесса, что особенно важно для анализа наблюдательных данных.
В основе calypso лежит метод главных компонент (PCA) и многомерное гауссово распределение. Коэффициенты разложения временных рядов моделируются как случайные величины, что позволяет эмулятору захватывать стохастичность аккреции и долгосрочную модуляцию из-за прецессии диска. Авторы также исследовали эпистемическую неопределенность, связанную с ограниченностью тренировочных данных, и пришли к выводу, что её включение не поддерживается имеющимися данными.
Для построения модели были использованы существующие наборы симуляций аккреции в двойных системах, а также проведены новые симуляции для заполнения пробелов в пространстве параметров (эксцентриситет орбиты e_b и отношение масс q_b). Для валидации эмулятора выполнено 13 тестовых симуляций. Разработчики опубликовали calypso в виде Python-пакета, устанавливаемого через pip, с открытым исходным кодом и подробной документацией.
Особенностью модели является вывод замкнутой гауссовой функции правдоподобия, которая позволяет напрямую оценивать параметры e_b и q_b по наблюдаемым временным рядам темпов аккреции. Это делает calypso удобным инструментом для интерпретации данных текущих и будущих обзоров транзиентов, таких как LSST и других.
Новый эмулятор снижает вычислительные затраты по сравнению с полномасштабными гидродинамическими симуляциями и может быть интегрирован в процессы быстрого анализа данных. Открытость кода способствует воспроизводимости и дальнейшему развитию инструмента.
Работа демонстрирует подход к моделированию стохастических процессов в астрофизике, который может быть применен и к другим типам аккреционных систем.
