2023MNRAS.525.2472Y

Based on the Gaia DR3 RR Lyrae catalogue, we use two methods to fit the density profiles with an improved broken power law, and find that there are two break radii coinciding with the two apocentre pile-ups of high-eccentricity Gaia-Sausage-Enceladus (GSE) merger. Also, there is a break caused by the Sagittarius (Sgr) stream. Combining the positions of all breaks, we briefly analyse the metallicity and its dispersion as a function of r as well as its distribution in cylindrical coordinates. For the clean sample, the z-x ellipsoid axial ratio q in $36\, {\rm kpc}\, \lt \, r\, \lt \, 96\, {\rm kpc}$ becomes much smaller than that of the inner halo $(r\, \lt \, 36\, {\rm kpc})$, while the major axis has a large uncertainty in the region of $36-66\, {\rm kpc}$ and the one in the region of $66-96\, {\rm kpc}$ is obviously different from that dominated by the Hercules-Aquila Cloud (HAC) and the Virgo Overdensity (VOD) in the inner halo, which indicates that there is an over-density structure distributed at low-zenithal angles. Finally, we found that the over-density structure in the outer halo ($r\, \gt \, 50\, {\rm kpc}$) is shell-shaped and relatively metal-rich compared to the outer background halo. We conclude that the shells could be the apocentre pile-ups of the high-eccentricity GSE merger, which is supported by previous numerical simulations.