Анализ звёздных движений помог определить структуру гало тёмной материи вокруг Млечного Пути

Звёздные потоки, формирующиеся при разрушении приливными силами карликовых галактик и шаровых скоплений, уже давно используются как инструменты для изучения гравитационного поля галактик. Их конфигурация непосредственно связана с распределением массы, главным образом определяя структуру гало тёмной материи. Однако при наблюдении за внешними галактиками данные о потоках обычно сводятся к двумерным фотометрическим изображениям, что затрудняет точное определение параметров гало по каждому объекту.

В рамках новой работы был предложен и реализован иерархический байесовский метод, который позволяет обойти это принципиальное ограничение. Вместо оценки формы гало для каждого потока по отдельности, исследователи применили байесовский подход, позволяющий объединить все имеющиеся данные и выявить общую тенденцию. Суть метода заключается в том, что каждый поток предоставляет лишь нечёткое и неоднозначное представление о гало, однако, объединив данные десятков потоков, можно получить стабильную картину. Байесовский анализ формально учитывает все возможные варианты, определяя, какие конфигурации наиболее вероятны для всей совокупности потоков. Это позволяет преодолеть фундаментальную неопределённость, связанную с отдельными наблюдениями, и сделать надёжные выводы о типичной форме гало тёмной материи.

Авторы провели обширный байесовский анализ, который потребовал сотен тысяч симуляций для каждого из нескольких десятков потоков данных.

Моделирование основано на методе «разбрызгивания частиц». Исходный поток перемещается в гравитационном поле родительской галактики, а звёзды отрываются от него вблизи точек Лагранжа L1 и L2. Затем частицы эволюционируют как объекты, не обладающие массой. Благодаря независимой динамике частиц, вычисления можно эффективно распараллеливать, что и реализовано в программном пакете StreaMAX. Этот подход обеспечил ускорение прямого моделирования звёздных потоков на несколько порядков.

Не пропустите:  Презентация полноприводного внедорожника Jetour T1 в России не состоится

Для сопоставления трёхмерных симуляций с реальными наблюдениями, поток проецируется на небесную сферу, откуда извлекается только центральный трек — «линия», определяющая геометрию потока. Данные о ширине и плотности намеренно исключаются, так как текущая модель не позволяет точно воспроизвести эти параметры. В ходе анализа рассматриваются исключительно геометрические ограничения, обусловленные формой трека.

Для определения формы и ориентации эллипса используется нестандартный метод параметризации: вместо углов вводится трёхмерный вектор, указывающий направление оси сплюснутости, а его длина служит параметром. Это позволяет задать априорное предположение о сплюснутости и избежать мультимодальности, свойственной параметрам, определяемым углами. С учетом симметрий, каждая модель потока описывается 13 параметрами.

По результатам изучения отдельных потоков установлено, что даже при точном определении исходных параметров апостериорные распределения сохраняют широкую и многомодальную структуру. Причина кроется в базовых физических неоднозначностях, в частности, в невозможности независимо установить масштаб массы гало и типичные скорости звёзд без данных о кинематике, а также в вырождениях, обусловленных проецированием трёхмерной структуры на двумерное изображение. Вследствие этого, гало, имеющие вытянутую и сплюснутую форму, и наблюдаемые под разными углами, могут демонстрировать практически идентичные траектории.

Не пропустите:  Представлен новый Chevrolet Bolt

Неопределенность, присущая этим потокам, делает их оптимальными для иерархического анализа. Авторы, сформировав выборку из 35 потоков, воспроизводящую каталог STRRINGS, использовали байесовскую модель, включающую в себя распределение сплюснутости гало в популяции. Применение метода перевзвешивания позволило исключить необходимость повторного всестороннего анализа каждого потока и значительно уменьшить объем вычислений.

Использование иерархической модели позволило с высокой степенью точности определить параметры исходных распределений для трех модельных популяций: сплюснутой, сферической и вытянутой. Мультимодальность, свойственная анализу отдельных потоков, была устранена при рассмотрении популяций в целом. Незначительное систематическое смещение оценок в сторону сферической формы обусловлено геометрическими особенностями и погрешностями ориентации, однако не мешает надёжному дифференцировании различных морфологий гало.

Авторы представляют эту работу как демонстрацию принципиальной возможности подхода и отмечают, что даже фотометрические данные, полученные при наблюдении за звёздными потоками, позволяют установить жёсткие ограничения на форму гало тёмной материи при использовании популяционного анализа. В перспективе планируется применение данного метода к фактическим каталогам наблюдений и усовершенствование модели.

Похожие статьи