Metallicity distributions of halo stars: do they trace the Galactic accretion history?

Abstract

The standard cosmological scenario predicts a hierarchical formation for galaxies. Many substructures have been found in the Galactic halo, usually identified as clumps in kinematic spaces, like the energy-angular momentum space ($E-L_z$), under the hypothesis that these quantities should be conserved during the interaction. If these clumps also feature different chemical abundances, such as the metallicity distribution function (MDF), these two arguments together (different kinematic and chemical properties) are often used to motivate their association with distinct and independent merger debris. The aim of this study is to explore to what extent we can couple kinematic characteristics and metallicities of stars in the Galactic halo to reconstruct the accretion history of the Milky Way (MW). In particular, we want to understand whether different clumps in the $E-L_z$ space with different MDFs should be associated with distinct merger debris. We analysed dissipationless, self-consistent, high-resolution N-body simulations of a MW-type galaxy accreting a satellite with a mass ratio of 1:10, with different orbital parameters and different metallicity gradients, which were assigned a posteriori. We confirm that accreted stars from a sim 1:10 mass ratio merger event redistribute in a wide range of $E$ and $L_z$, due to the dynamical friction process, and are thus not associated with a single region. Because satellite stars with different metallicities can be deposited in different regions of the $E-L_z$ space (on average the more metal-rich ones end up more gravitationally bound to the MW), this implies that a single accretion of sim 1:10 can manifest with different MDFs, in different regions of the $E-L_z$ space. Groups of stars with different $E$, $L_z$, and metallicities may be interpreted as originating from different satellite galaxies, but our analysis shows that these interpretations are not physically motivated. In fact, as we show, the coupling of kinematic information with MDFs to reconstruct the accretion history of the MW can bias the reconstructed merger tree towards increasing the number of past accretions and decreasing the masses of the progenitor galaxies.