Constrained cosmological simulations of the Local Group using Bayesian hierarchical field-level inference

TL;DR

This paper introduces a Bayesian hierarchical inference method to simulate the Local Group in detail, accurately constraining galaxy masses and structures within a cosmological framework.

Contribution

It extends the BORG algorithm with multi-resolution capabilities to resolve individual galaxies and incorporate diverse observational data for detailed Local Group modeling.

Findings

Estimated MW and M31 masses with high precision.

Revealed the nearly radial orbit of the M31-MW pair.

Showed the Local Group's structure aligns with the Supergalactic Plane.

Abstract

We present a novel approach based on Bayesian field-level inference capable of resolving individual galaxies within the Local Group (LG), enabling detailed studies of its structure and formation via posterior simulations. We extend the Bayesian Origin Reconstruction from Galaxies (BORG) algorithm with a multi-resolution approach, allowing us to reach smaller mass scales and apply observational constraints based on LG galaxies. Our updated data model simultaneously accounts for observations of mass tracers within the dark haloes of the Milky Way (MW) and M31, their observed separation and relative velocity, and the quiet surrounding Hubble flow represented through the positions and velocities of galaxies at distances from one to four Mpc. Our approach delivers representative posterior samples of $\Lambda$CDM realisations that are statistically and simultaneously consistent with all these observations, leading to significantly tighter mass constraints than found if the individual datasets are considered separately. In particular, we estimate the virial masses of the MW and M31 to be $\log_{10}(M_{200c}/M_\odot) = 12.07\pm0.08$ and $12.33\pm0.10$, respectively, their sum to be $\log_{10}(\Sigma M_{200c}/M_\odot)= 12.52\pm0.07$, and the enclosed mass within spheres of radius $R$ to be $\log_{10}(M(R)/M_\odot)= 12.71\pm0.06$ and $12.96\pm0.08$ for $R=1$ Mpc and 3 Mpc, respectively. The M31-MW orbit is nearly radial for most of our $\Lambda$CDM LG's, and most lie in a dark matter sheet that aligns approximately with the Supergalactic Plane, even though the surrounding density field was not used explicitly as a constraint. The approximate simulations employed in our inference are accurately reproduced by high-fidelity structure formation simulations, demonstrating the potential for future high-resolution, full-physics $\Lambda$CDM posterior simulations of LG look-alikes.