Uncovering the first-infall history of the LMC through its dynamical impact in the Milky Way halo

Yanjun Sheng, Yuan-Sen Ting, Xiang-Xiang Xue, Jiang Chang, Hao Tian

Published 2024-04-13Version 1

The gravitational interactions between the LMC and the Milky Way can give rise to dynamical perturbations in the MW halo, leading to a biased distribution of stellar density and other kinematic signals. These disequilibrium phenomena exhibit variations under different parameter combinations of the MW-LMC model. In this work, we run 50 high-resolution N-body simulations spanning different masses and halo shapes of the Milky Way and LMC and investigate how the LMC-induced perturbations evolve with these model parameters. We measure the magnitude of kinematic perturbations from the mean velocities of simulated halo stars and identify a discontinuity between the first-infall and second-passage scenarios of the LMC's orbital history. We demonstrate that, due to the short dynamical times of the Galactic inner halo, the reduced perturbation magnitude in the second-passage scenario is mainly a result of the LMC's second infall into the MW, which starts at a much lower velocity relative to the inner halo compared to the first-infall scenario. Using a subset of $\sim 1200$ RR Lyrae stars located in the outer halo ($50 \leq R_{\mathrm{GC}} < 100$ kpc), which are selected from a larger sample of 135,873 RR Lyrae stars with precise distance estimates from Gaia, we find the mean latitudinal velocity ($v_{b}$) in the heliocentric frame to be $\langle v_{b} \rangle = 30.8 \pm 4.0$ km/s. The observation contradicts the second-passage scenario and supports the first-infall scenario with a massive LMC ($\sim 2.1 \times 10^{11} \mathrm{M}_{\odot}$) at infall, an oblate MW halo with a virial mass $M_{200} < 1.4 \times 10^{12} \mathrm{M}_{\odot}$ and a flattening parameter $q > 0.7$.