An official website of the United States government
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 shock from the supersonic CGM–CGM collision. The resulting shock front extends from the LMC to beyond ∼R200,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, withvR ± 30–50 km s−1atR > 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 ∼108−9M⊙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.
more »
« less
The Clustering of Orbital Poles Induced by the LMC: Hints for the Origin of Planes of Satellites
Abstract A significant fraction of Milky Way (MW) satellites exhibit phase-space properties consistent with a coherent orbital plane. Using tailored N -body simulations of a spherical MW halo that recently captured a massive (1.8 × 10 11 M ⊙ ) LMC-like satellite, we identify the physical mechanisms that may enhance the clustering of orbital poles of objects orbiting the MW. The LMC deviates the orbital poles of MW dark matter particles from the present-day random distribution. Instead, the orbital poles of particles beyond R ≈ 50 kpc cluster near the present-day orbital pole of the LMC along a sinusoidal pattern across the sky. The density of orbital poles is enhanced near the LMC by a factor δ ρ max = 30% (50%) with respect to underdense regions and δ ρ iso = 15% (30%) relative to the isolated MW simulation (no LMC) between 50 and 150 kpc (150–300 kpc). The clustering appears after the LMC’s pericenter (≈50 Myr ago, 49 kpc) and lasts for at least 1 Gyr. Clustering occurs because of three effects: (1) the LMC shifts the velocity and position of the central density of the MW’s halo and disk; (2) the dark matter dynamical friction wake and collective response induced by the LMC change the kinematics of particles; (3) observations of particles selected within spatial planes suffer from a bias, such that measuring orbital poles in a great circle in the sky enhances the probability of their orbital poles being clustered. This scenario should be ubiquitous in hosts that recently captured a massive satellite (at least ≈1:10 mass ratio), causing the clustering of orbital poles of halo tracers.
more »
« less
- Award ID(s):
- 1941096
- PAR ID:
- 10334879
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 923
- Issue:
- 2
- ISSN:
- 0004-637X
- Page Range / eLocation ID:
- 140
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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.more » « less
-
ABSTRACT A variety of observational campaigns seek to test dark matter models by measuring dark matter subhaloes at low masses. Despite their predicted lack of stars, these subhaloes may be detectable through gravitational lensing or via their gravitational perturbations on stellar streams. To set measurable expectations for subhalo populations within Lambda cold dark matter, we examine 11 Milky Way (MW)-mass haloes from the FIRE-2 baryonic simulations, quantifying the counts and orbital fluxes for subhaloes with properties relevant to stellar stream interactions: masses down to $$10^{6}\, \text{M}_\odot$$, distances ≲50 kpc of the galactic centre, across z = 0 − 1 (tlookback = 0–8 Gyr). We provide fits to our results and their dependence on subhalo mass, distance, and lookback time, for use in (semi)analytical models. A typical MW-mass halo contains ≈16 subhaloes $$\gt 10^{7}\, \text{M}_\odot$$ (≈1 subhalo $$\gt 10^{8}\, \text{M}_\odot$$) within 50 kpc at z ≈ 0. We compare our results with dark matter-only versions of the same simulations: because they lack a central galaxy potential, they overpredict subhalo counts by 2–10×, more so at smaller distances. Subhalo counts around a given MW-mass galaxy declined over time, being ≈10× higher at z = 1 than at z ≈ 0. Subhaloes have nearly isotropic orbital velocity distributions at z ≈ 0. Across our simulations, we also identified 4 analogues of Large Magellanic Cloud satellite passages; these analogues enhance subhalo counts by 1.4–2.1 times, significantly increasing the expected subhalo population around the MW today. Our results imply an interaction rate of ∼5 per Gyr for a stream like GD-1, sufficient to make subhalo–stream interactions a promising method of measuring dark subhaloes.more » « less
-
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.more » « less
-
ABSTRACT Understanding the evolution of satellite galaxies of the Milky Way (MW) and M31 requires modelling their orbital histories across cosmic time. Many works that model satellite orbits incorrectly assume or approximate that the host halo gravitational potential is fixed in time and is spherically symmetric or axisymmetric. We rigorously benchmark the accuracy of such models against the FIRE-2 cosmological baryonic simulations of MW/M31-mass haloes. When a typical surviving satellite fell in ($$3.4\!-\!9.7\, \rm {Gyr}$$ ago), the host halo mass and radius were typically 26–86 per cent of their values today, respectively. Most of this mass growth of the host occurred at small distances, $$r\lesssim 50\, \rm {kpc}$$, opposite to dark matter only simulations, which experience almost no growth at small radii. We fit a near-exact axisymmetric gravitational potential to each host at z = 0 and backward integrate the orbits of satellites in this static potential, comparing against the true orbit histories in the simulations. Orbital energy and angular momentum are not well conserved throughout an orbital history, varying by 25 per cent from their current values already $$1.6\!-\!4.7\, \rm {Gyr}$$ ago. Most orbital properties are minimally biased, ≲10 per cent, when averaged across the satellite population as a whole. However, for a single satellite, the uncertainties are large: recent orbital properties, like the most recent pericentre distance, typically are ≈20 per cent uncertain, while earlier events, like the minimum pericentre or the infall time, are ≈40–80 per cent uncertain. Furthermore, these biases and uncertainties are lower limits, given that we use near-exact host mass profiles at z = 0.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3847/1538-4357/ac2c05
