Cosmological simulations of the same spiral galaxy: satellite properties, the role of baryonic physics and star formation history in shaping dark matter cores/cusps

We investigate the role of baryonic physics in shaping the population, structure, and internal dynamics of galactic subhalos using the Mochima suite of cosmological zoom-in simulations. A refined method is developed to identify bound subhalo material by isolating the local gravitational potential and applying multi-criteria phase-space selection. This approach enables a robust characterisation of subhalo properties across five baryonic runs with varying prescriptions for star formation, and supernova and protostellar feedback, as well as a dark matter-only baseline. At the population level, we find that the concentration of the central massive host halo, modulated by baryonic physics in the central disc galaxy, is a key predictor of subhalo survival. Subhalos with more massive stellar components exhibit deeper internal potentials and enhanced resilience to tidal disruption. At the structural level, we identify a broad diversity in inner dark matter profiles, consistent with observations of dwarf galaxies. We show that this diversity correlates with both star formation history and environmental interaction. In particular, galaxies that form most of their stars early tend to retain steep cusps, while those with extended or recent star formation exhibit oscillating inner slopes shaped by bursty feedback and tidal perturbations. These findings suggest that the so-called “diversity problem” may reflect the complex interplay between feedback history and gravitational environment, rather than a breakdown of cold dark matter predictions.