The Milky Way Stellar Halo Is Twisted and Doubly Broken: Insights from DESI DR2 Milky Way Survey Observation

Using K giants from the second data release (DR2) of the Dark Energy Spectroscopic Instrument (DESI) Milky Way Survey, we measure the shape, orientation, radial profile, and density anisotropies of the Milky Way (MW) stellar halo over 8 kpc < rGC < 200 kpc. We identify a triaxial stellar halo (axis ratio 10:8:7), 43 degrees tilted from the disk, showing two break radii at similar to 16 and similar to 76 kpc, likely associated with Gaia-Sausage/Enceladus and the Large Magellanic Cloud (LMC), respectively. The inner stellar halo (<30 kpc) is oblate and aligned with the disk, whereas the outer stellar halo becomes prolate and perpendicular to the disk, consistent with the vast polar structure of MW satellites. The twisted halo may arise from the disk-halo angular momentum shift triggered by the infall of a massive satellite. The anisotropic density distribution of the stellar halo is also measured, with successful reidentification of the Hercules-Aquila Cloud North/South (HAC-N/S) overdensity and the Virgo overdensity (VOD). Break radii are found at 15 and 30 kpc for VOD and HAC-N/S, respectively. We identify the LMC transient density wake with a break radius at 60 kpc in the Pisces overdensity region. We also find new observational evidence of the LMC collective density wake, by showing a break radius at similar to 100 kpc in the northern Galactic cap with a clear density peak at 90 kpc. In the end, we found that more metal-poor halo stars are more radially extended. Our results provide important clues to the assembly and evolution of the MW stellar halo under the standard cosmic structure formation framework.