How well can we unravel the accreted constituents of the Milky Way stellar halo? A test on cosmological hydrodynamical simulations

Abstract: Context. One of the primary goals of Galactic Archaeology is to reconstruct the Milky Way's accretion history. To achieve this, significant efforts have been dedicated to identifying signatures of past accretion events. In particular, the study of integrals-of-motion (IoM) space has proven to be highly insightful for uncovering these ancient mergers and understanding their impact on the Galaxy's evolution. Aims. This paper evaluates the effectiveness of a state-of-the-art method for detecting debris from accreted galaxies, by testing it on four Milky Way-like galaxies from the Auriga suite of cosmological magneto-hydrodynamical simulations. Methods. We employ the innovative method from L\"ovdal et al. (2022) to identify substructures in the integrals-of-motion space within the local stellar halos of the four simulated galaxies. This approach enables us to evaluate the method's performance by comparing the properties of the identified clusters with the known populations of accreted galaxies in the simulations. Additionally, we investigate whether incorporating chemical abundances and stellar age information can help to link distinct structures originating from the same accretion event. Results. This method is very effective in detecting debris from accretion events that occur less than 6-7 Gyr ago but struggles to detect most of the debris from older accretion. Furthermore, most of the detected structures suffer from significant contamination by in-situ stars. Our results also show that the method may also generate artificial detections. Conclusions. Our work show that the Milky Way's accretion history remains uncertain, and question the reality of some detected structures in the Solar vicinity.