Quantifying the Milky Way, LMC and their interaction using all-sky kinematics of outer halo stars

TL;DR

This study uses all-sky stellar kinematics data and simulations to measure the Milky Way and LMC masses, reflex motion, and halo anisotropy, highlighting the importance of including the LMC in models.

Contribution

It introduces a simulation-based inference method to jointly constrain MW and LMC masses, reflex motion, and halo anisotropy from outer halo star data.

Findings

Reflex velocity magnitude is approximately 39 km/s.

Enclosed MW mass within 50 kpc is about 3.36 x 10^{11} solar masses.

LMC mass within 50 kpc is roughly 8.76 x 10^{10} solar masses.

Abstract

The recent pericentric passage of the Large Magellanic Cloud (LMC) has dislodged the Milky Way's (MW) centre of mass, inducing dynamical disequilibrium, the reflex motion, in the kinematics of outer stellar halo stars. Using data out to $160 \, \rm kpc$ from the combined H3+SEGUE+MagE outer halo survey, we constrain the mass of the MW and LMC, as well as the resulting reflex motion and the stellar halo velocity anisotropy. Using a suite of 32,000 rigid MW--LMC simulations, each with a MW stellar halo evolved to the present day in the combined MW--LMC potential, we perform Simulation Based Inference by training a neural posterior estimator on the means and dispersions of the radial and tangential velocities of stars from the combined H3+SEGUE+MagE outer halo sample. Relative to halo stars at $100 \, \rm kpc$, we find the magnitude of the reflex velocity to be $v_{\rm travel} = 38.6^{+8.3}_{-7.8}\,\rm km \, s^{-1}$. Simultaneously, we determine the enclosed MW mass, $M_{\rm MW}(< 50 \, \rm kpc) = 3.36 \pm 0.15 \times 10^{11}\, \rm M_{\odot}$ and the enclosed LMC mass, $M_{\rm LMC}(< 50 \, \rm kpc) = 8.76^{+1.94}_{-1.77} \times 10^{10}\, \rm M_{\odot}$. Our results suggest that the total LMC mass must be at least $\sim20\%$ that of the MW. The velocity anisotropy prior to the LMC's infall is constrained to be $\beta_0 = 0.68 \pm 0.02$. Finally, we demonstrate that neglecting the LMC in models biases the estimated MW mass to prefer more massive values.