We perform a series of controlled N-body simulations to study realizations of the recently discovered Antlia 2 galaxy in cold dark matter (CDM) and self-interacting dark matter (SIDM) scenarios. Our simulations contain six benchmark models, where we vary the initial halo concentration and the self-scattering cross section. We adopt well-motivated initial stellar and halo masses, and our fiducial orbit has a small pericenter. After evolving in the Milky Way's tidal field, the simulated galaxies experience significant mass loss and their stellar distributions expand accordingly. These tidal effects are more prominent if the initial halo concentration is lower and if the self-scattering cross section is larger. Our results show that Antlia 2-like galaxies could be realized in CDM if the halo concentration is low and the stellar distribution is diffuse at the infall time, while these conditions could be relaxed in SIDM. We also find all the simulated galaxies predict approximately the same stellar velocity dispersion after imposing selection criteria for stellar particles. This has important implications for testing dark matter models using tidally disturbed systems.

We study evolution of self-interacting dark matter (SIDM) subhalos in the Milky Way (MW) tidal field. The interaction between the subhalos and the MW's tides lead to more diverse dark matter distribution in the inner region, compared to their cold dark matter counterparts. We test this scenario with two MW satellite galaxies, Draco and Fornax, opposite extremes in the inner dark matter content, and find that they can be accommodated within the SIDM model proposed to explain the diverse rotation curves of spiral galaxies in the field.

We study halo mass functions with the Press-Schechter formalism for interacting dark matter models, where matter power spectra are damped due to dark acoustic oscillations in the early universe. After adopting a smooth window function, we calibrate the analytical model with numerical simulations from the "effective theory of structure formation" (ETHOS) project and fix the model parameters in the high mass regime, $M_{\rm h}\gtrsim3\times10^{10}\;{\rm M}_{\odot}$. We also perform high-resolution cosmological simulations with halo masses down to $M_{\rm h}\sim10^8\;{\rm M}_{\odot}$ to cover a wide mass range for comparison. Although the model is calibrated with ETHOS1 and CDM simulations for high halo masses at redshift $z=0$, it successfully reproduces simulations for two other ETHOS models in the low mass regime at low and high redshifts. As an application, we compare the cumulative number density of haloes to that of observed galaxies at $z=6$, and find the interacting dark matter models with a kinetic decoupling temperature below $0.5\ \rm{keV}$ is disfavored. We also perform the abundance-matching analysis and derive the stellar-halo mass relation for these models at $z=4$. Suppression in halo abundance leads to less massive haloes that host observed galaxies in the stellar mass range $M_*\simeq 10^5-10^7\ {\rm M}_{\odot}$.

We explore a scenario where metal poor globular clusters (GCs) form at the centres of their own dark matter halos in the early universe before reionization. This hypothesis leads to predictions about the abundance, distribution and kinematics of GCs today that we explore using cosmological N-body simulations and analytical modelling. We find that selecting the massive tail of collapsed objects at $z\gtrsim 9$ as GC formation sites leads to four main predictions: i) a highly clustered population of GCs around galaxies today, ii) a natural scaling between number of GCs and halo virial mass that follows roughly the observed trend, iii) a very low number of free floating GCs outside massive halos and iv) GCs should be embedded within massive and extended dark matter (sub)halos. We find that the strongest constraint to the model is given by the combination of (i) and (ii): a mass cut to tagged GCs halos which accounts for the number density of metal poor GCs today predicts a radial distribution that is too extended compared to recent observations. On the other hand, a mass cut sufficient to match the observed half number radius could only explain 60% of the metal poor population. In all cases, observations favour early redshifts for GC formation ($z\geq 15$).

We employ isolated N-body simulations to study the response of self-interacting dark matter (SIDM) halos in the presence of the baryonic potentials. Dark matter self-interactions lead to kinematic thermalization in the inner halo, resulting in a tight correlation between the dark matter and baryon distributions. A deep baryonic potential shortens the phase of SIDM core expansion and triggers core contraction. This effect can be further enhanced by a large self-scattering cross section. We find the final SIDM density profile is sensitive to the baryonic concentration and the strength of dark matter self-interactions. Assuming a spherical initial halo, we also study evolution of the SIDM halo shape together with the density profile. The halo shape at later epochs deviates from spherical symmetry due to the influence of the non-spherical disc potential, and its significance depends on the baryonic contribution to the total gravitational potential, relative to the dark matter one. In addition, we construct a multi-component model for the Milky Way, including an SIDM halo, a stellar disc and a bulge, and show it is consistent with observations from stellar kinematics and streams.

Galactic rotation curves are a fundamental constraint for any cosmological model. We use controlled N-body simulations of galaxies to study the gravitational effect of baryons in a scenario with collisionless cold dark matter (CDM) versus one with a self-interacting dark matter (SIDM) component. In particular, we examine the inner profiles of the rotation curves in the velocity range Vmax = 30-250 km/s, whose diversity has been found to be greater than predicted by the Lambda-CDM scenario. We find that the scatter in the observed rotation curves exceeds that predicted by dark matter only mass-concentration relations in either the CDM nor SIDM models. Allowing for realistic baryonic content and spatial distributions, however, helps create a large variety of rotation curve shapes; which is in better agreement with observations in the case of self-interactions due to the characteristic cored profiles being more accommodating to the slowly rising rotation curves than CDM. We find individual fits to model two of the most remarkable outliers of similar Vmax, UGC 5721 and IC 2574, the former a cusp-like rotation curve and the latter a seemingly 8 kpc cored profile. This diversity in SIDM arises as permutations of overly concentrated halos with compact baryonic distributions versus underdense halos with extended baryonic disks. The SIDM solution is promising and its feasibility ultimately depends on the sampling of the halo mass-concentration relation and its interplay with the baryonic profiles, emphasising the need for a better understanding of the frequency of extreme outliers present in current observational samples.