Untangling stellar components of galaxies: Evaluation of dynamical decomposition methods in simulated galaxies with GalaxyChop

Abstract

Galaxy formation is intrinsically connected to the distinct evolutionary processes of disk and spheroidal systems, which are the fundamental stellar components of galaxies. Understanding the mutual, dynamical interplay and co-evolution of these components requires a detailed dynamical analysis to allow for disentanglement of these systems. We introduce JEHistogram, a new method for the dynamical decomposition of simulated galaxies into disk and spheroidal stellar components that utilizes the angular momentum and energy of star particles. We evaluate its performance against five previously established methods using a sample of equilibrium galaxies with stellar masses in the range 1010 ≤ Mgal/M⊙ ≤ 1012. Our assessment involves several metrics, including the completeness and purity of stellar particle classification, scale lengths, mass density profiles, velocity dispersion, and rotational velocity profiles. While all methods approximate the properties of the original components, such as mass fractions and density or velocity profiles, JEHistogram demonstrates a better accuracy, particularly in the inner regions of galaxies where component overlap complicates separation. Additionally, we apply JEHistogram to a Milky Way-like galaxy from the IllustrisTNG cosmological simulations, showcasing its capability to derive properties such as the size, mass, velocity, color, and age of dynamically defined disk and spheroidal components. All the dynamical decomposition methods we analyzed are publicly accessible through the Python package GalaxyChop.