Search for: All records

Abstract The large-scale morphology of Milky Way (MW)–mass dark matter halos is shaped by two key processes: filamentary accretion from the cosmic web and interactions with massive satellites. Disentangling their contributions is essential for understanding galaxy evolution and constructing accurate mass models of the MW. We analyze the time-dependent structure of MW-mass halos from zoomed cosmological-hydrodynamical simulations by decomposing their mass distribution into spherical harmonic expansions. We find that the dipole and quadrupole moments dominate the gravitational power spectrum, encoding key information about the halo’s shape and its interaction with the cosmic environment. While the dipole reflects transient perturbations from infalling satellites and damps on dynamical timescales, the quadrupole—linked to the halo’s triaxiality—is a persistent feature. We show that the quadrupole’s orientation aligns with the largest filaments, imprinting a long-lived memory on the halo’s morphology even in its inner regions (∼30 kpc). At the virial radius, the quadrupole distortion can reach 1–2 times the spherical density, highlighting the importance of environment in shaping MW-mass halos. Using multichannel singular spectrum analysis, we successfully disentangle the effects of satellite mergers and filamentary accretion on quadrupole. We find that, compared to isolated MW–LMC simulations that typically use a spherical halo, the LMC-mass satellite induces a quadrupolar response that is an order of magnitude larger in our cosmological halo. This highlights the need for models that incorporate the MW’s asymmetry and time evolution, with direct consequences for observable structures such as disk warps, the LMC-induced wake, and stellar tracers—particularly in the era of precision astrometry. 

Abstract The LMC’s stellar bar is offset from the outer disk center, tilted from the disk plane, and does not drive gas inflows. These properties are atypical of bars in gas-rich galaxies, yet the LMC bar’s strength and radius are similar to typical barred galaxies. UsingN-body hydrodynamic simulations, we show that the LMC’s unusual bar is explainable if there was a recent collision (impact parameter ≈2 kpc) between the LMC and SMC. Pre-collision, the simulated bar is centered and coplanar. Post-collision, the simulated bar is offset (≈1.5 kpc) and tilted (≈8 $\stackrel{\text{°}}{.}$6). The simulated bar offset reduces with time, and comparing with the observed offset (≈0.8 kpc) suggests the timing of the true collision to be 150–200 Myr ago. Then, 150 Myr post-collision, the LMC’s bar is centered with its dark matter (DM) halo, whereas the outer disk center is separated from the DM center by ≈1 kpc. The SMC collision produces a tilted-ring structure for the simulated LMC, consistent with observations. Post-collision, the simulated LMC bar’s pattern speed decreases by a factor of 2. We also provide a generalizable framework to quantitatively compare the LMC’s central gas distribution in different LMC–SMC interaction scenarios. We demonstrate that the SMC’s torques on the LMC’s bar during the collision are sufficient to explain the observed bar tilt, provided the SMC’s total mass within 2 kpc was (0.8–2.4) × 10⁹M_⊙. Therefore, the LMC bar’s tilt constrains the SMC’s pre-collision DM profile, and requires the SMC to be a DM-dominated galaxy. 

Abstract The infall of the Large Magellanic Cloud (LMC) into the Milky Way’s halo impacts the distribution of stars and dark matter (DM) in our Galaxy. Mapping the observational consequences of this encounter can inform us about the properties of both galaxies, details of their interaction, and possibly distinguish between different DM models.N-body simulations predict a localized overdensity trailing the LMC’s orbit both in baryonic and DM, known as the wake. We collected wide-field, deep near-infrared, and optical photometry using VIRCAM and DECam across four fields along the expected wake, covering the sky region expected to span most of its predicted density contrast. We identify over 400 stars comprising two different tracers, near main-sequence turnoff stars and red giants, which map the halo between 60 and 100 kpc, deriving stellar halo densities as a function of sky position and Galactocentric radius. We detect (1) a break in the halo radial density profile at 70 kpc not seen in northern halo studies and (2) a clear halo overdensity starting also at 70 kpc, with density contrast increasing steadily toward the expected current location of the wake. If this overdensity is the LMC wake, its peak density contrast is as pronounced as the most massive LMC model considered. Contamination from unidentified substructures may bias our wake detections, so wider-area surveys with similar depth are needed for confirmation. 

Abstract The first infall of the LMC into the Milky Way (MW) represents a large and recent disruption to the MW circumgalactic medium (CGM). In this work, we use idealized, hydrodynamical simulations of an MW-like CGM embedded in a dark matter halo with an infalling LMC-like satellite initialized with its own CGM to understand how the encounter is shaping the global physical and kinematic properties of the MW CGM. First, we find that the LMC drives order-unity enhancements in MW CGM density, temperature, and pressure due to a $\mathcal{M}\approx 2$shock from the supersonic CGM–CGM collision. The resulting shock front extends from the LMC to beyond ∼R_200,MW, amplifying column densities, X-ray brightness, thermal Sunyaev–Zeldovich distortion, and potentially synchrotron emission from cosmic rays over large angular scales across the southern hemisphere. Second, the MW’s reflex motion relative to its outer halo induces a dipole in CGM radial velocities, withv_R ± 30–50 km s⁻¹atR > 50 kpc in the northern and southern hemispheres, respectively, consistent with measurements in the stellar halo. Finally, ram pressure strips most of the LMC’s CGM, leaving ∼10⁸⁻⁹M_⊙warm ionized gas along the past orbit of the LMC, moving at high radial and/or tangential velocities ∼50–100 kpc from the MW. Massive satellites like the LMC leave their mark on the CGM structure of their host galaxies, and signatures of such interactions may be observable in key all-sky tracers of the MW CGM and those of other massive galaxies. 

Abstract We present new measurements of the 2D geometry of the LMC’s stellar bar with precise astrometric observations of red clump stars in Gaia DR3. We develop a novel solution to tackle crowding-induced incompleteness in Gaia data sets with the Gaia BP-RP color excess. Utilizing the color excess information, we derive a 2D completeness map of the LMC’s disk. We find that incompleteness biases the bar measurements and induces large uncertainties. With the completeness-corrected 2D red clump map, we precisely measure the LMC bar’s properties using Fourier decomposition. The bar radius (semimajor axis) is $R_{\mathrm{bar}}={2.13}_{-0.04}^{+0.03}$kpc, and its position angle is 121.°26 ± 0.°21. The bar’s strength as quantified by the Fourier bi-symmetric amplitude isS_bar= 0.27, indicating that the LMC has a significant bar perturbation. We find the bar has an axis ratio of 0.54 ± 0.03, and is offset with respect to the center of the outer disk isophote atR≈ 5 kpc by 0.76 ± 0.01 kpc. These LMC bar properties agree with a hydrodynamic model where the SMC has undergone a recent direct collision with the LMC. We compare the LMC’s bar properties with other barred galaxies in the local Universe, and discover that the LMC is similar to other barred galaxies in terms of bar-galaxy scaling relations. We discuss how our completeness correction framework can be applied to other systems in the Local Group. 

Abstract Stellar streams in the Milky Way are promising detectors of low-mass dark matter (DM) subhalos predicted by ΛCDM. Passing subhalos induce perturbations in streams that indicate the presence of the subhalos. Understanding how known DM-dominated satellites impact streams is a crucial step toward using stream perturbations to constrain the properties of dark perturbers. Here, we cross-match a Gaia Early Data Release 3 and SEGUE member catalog of the Cetus-Palca stream (CPS) with H3 for additional radial velocity measurements and fit the orbit of the CPS using this six-dimensional (6D) data. We demonstrate for the first time that the ultra-faint dwarf Segue 2 had a recent (77 ± 5 Myr ago) close flyby (within the stream's 2σwidth) with the CPS. This interaction enables constraints on Segue 2’s mass and density profile at larger radii ( $O\left(1\right)$kpc) than are probed by its stars ( $O\left(10\right)$pc). While Segue 2 is not expected to strongly affect the portion of the stream covered by our 6D data, we predict that if Segue 2’s mass within ∼ 6 kpc is 5 × 10⁹M_⊙, the CPS's velocity dispersion will be ∼ 40 km s⁻¹larger atϕ₁ > 20° than atϕ₁ < 0°. If no such heating is detected, Segue 2’s mass cannot exceed 10⁹M_⊙within ∼ 6 kpc. The proper motion distribution of the CPS near the impact site is mildly sensitive to the shape of Segue 2’s density profile. This study presents a critical test for frameworks designed to constrain properties of dark subhalos from stream perturbations. 

Abstract The merger timescales of isolated low-mass pairs (10⁸<M_*< 5 × 10⁹M_⊙) on cosmologically motivated orbits have not yet been studied in detail, though isolated high-mass pairs (5 × 10⁹<M_*< 10¹¹M_⊙) have been studied extensively. It is common to apply the same separation criteria and expected merger timescales of high-mass pairs to low-mass systems, however, it is unclear if their merger timescales are similar, or if they evolve similarly with redshift. We use the Illustris TNG100 simulation to quantify the merger timescales of isolated low-mass and high-mass major pairs as a function of cosmic time, and explore how different selection criteria impact the mass and redshift dependence of merger timescales. In particular, we present a physically motivated framework for selecting pairs via a scaled separation criterion, wherein pair separations are scaled by the virial radius of the primary’s Friends-of-Friends (FoF) group halo (r_sep< 1R_vir). Applying these scaled separation criteria yields equivalent merger timescales for both mass scales at all redshifts. Alternatively, static physical separation selections applied equivalently to all galaxy pairs at all redshifts lead to a difference in merger rate of up to ∼1 Gyr between low- and high-mass pairs, particularly forr_sep< 150 kpc. As a result, applying the same merger timescales to physical-separation-selected pairs will lead to a bias that systematically overpredicts low-mass galaxy merger rates. 

Abstract We study the effect of the Large Magellanic Cloud (LMC) on the dark matter (DM) distribution in the Solar neighborhood, utilizing the Auriga magneto-hydrodynamical simulations of Milky Way (MW) analogues that have an LMC-like system. We extract the local DM velocity distribution at different times during the orbit of the LMC around the MW in the simulations. As found in previous idealized simulations of the MW-LMC system, we find that the DM particles in the Solar neighborhood originating from the LMC analogue dominate the high speed tail of the local DM speed distribution. Furthermore, the native DM particles of the MW in the Solar region are boosted to higher speeds as a result of a response to the LMC's motion.We simulate the signals expected in near future xenon, germanium, and silicon direct detection experiments, considering DM interactions with target nuclei or electrons. We find that the presence of the LMC causes a considerable shift in the expected direct detection exclusion limits towards smaller cross sections and DM masses, with the effect being more prominent for low mass DM. Hence, our study shows, for the first time, that the LMC's influence on the local DM distribution is significant even in fully cosmological MW analogues. 

Abstract The interaction between the supersonic motion of the Large Magellanic Cloud (LMC) and the circumgalactic medium (CGM) is expected to result in a bow shock that leads the LMC’s gaseous disk. In this letter, we use hydrodynamic simulations of the LMC’s recent infall to predict the extent of this shock and its effect on the Milky Way’s (MW) CGM. The simulations clearly predict the existence of an asymmetric shock with a present-day standoff radius of ∼6.7 kpc and a transverse diameter of ∼30 kpc. Over the past 500 Myr, ∼8% of the MW’s CGM in the southern hemisphere should have interacted with the shock front. This interaction may have had the effect of smoothing over inhomogeneities and increasing mixing in the MW CGM. We find observational evidence of the existence of the bow shock in recent Hαmaps of the LMC, providing a potential explanation for the envelope of ionized gas surrounding the LMC. Furthermore, the interaction of the bow shock with the MW CGM may also explain the observations of ionized gas surrounding the Magellanic Stream. Using recent orbital histories of MW satellites, we find that many satellites have likely interacted with the LMC shock. Additionally, the dwarf galaxy Ret2 is currently sitting inside the shock, which may impact the interpretation of the reported gamma-ray excess in Ret2. This work highlights how bow shocks associated with infalling satellites are an underexplored yet potentially very important dynamical mixing process in the circumgalactic and intracluster media. 

Abstract Low-mass galaxy pair fractions are understudied, and it is unclear whether low-mass pair fractions evolve in the same way as more massive systems over cosmic time. In the era of JWST, Roman, and Rubin, selecting galaxy pairs in a self-consistent way will be critical to connect observed pair fractions to cosmological merger rates across all mass scales and redshifts. Utilizing the Illustris TNG100 simulation, we create a sample of physically associated low-mass (10⁸<M_*< 5 × 10⁹M_⊙) and high-mass (5 × 10⁹<M_*< 10¹¹M_⊙) pairs betweenz= 0 and 4.2. The low-mass pair fraction increases fromz= 0 to 2.5, while the high-mass pair fraction peaks atz= 0 and is constant or slightly decreasing atz> 1. Atz= 0 the low-mass major (1:4 mass ratio) pair fraction is 4× lower than high-mass pairs, consistent with findings for cosmological merger rates. We show that separation limits that vary with the mass and redshift of the system, such as scaling by the virial radius of the host halo (r_sep< 1R_vir), are critical for recovering pair fraction differences between low-mass and high-mass systems. Alternatively, static physical separation limits applied equivalently to all galaxy pairs do not recover the differences between low- and high-mass pair fractions, even up to separations of 300 kpc. Finally, we place isolated mass analogs of Local Group galaxy pairs, i.e., Milky Way (MW)–M31, MW–LMC, LMC–SMC, in a cosmological context, showing that isolated analogs of LMC–SMC-mass pairs and low-separation (<50 kpc) MW–LMC-mass pairs are 2–3× more common atz≳ 2–3.