A possibility of DM being multicomponent has a strong implication on resolving decades-long known cosmological problems on small scale. In addition to elastic scattering, the model allows for inelastic interactions, which can be characterized by a 'velocity kick' parameter. The simplest 2cDM model with cross section $0.01\lesssim\sigma/m<1\textrm{ cm}^{2}{ \rm g}^{-1}$ and the kick velocity $V_{k}\simeq 100\textrm{ km s}^{-1}$ has been shown to robustly resolve the missing satellites, core-cusp, and too-big-to-fail problems in $N$-body cosmological simulations tested on MW-like haloes of a virial mass $\sim5 \times 10^{11}$ M$_{\odot}$ (Paper I $\&$ II). With the aim of further constraining the parameter space available for the 2cDM model, we extend our analysis to dwarf and galaxy cluster haloes with their virial mass of $\sim 10^7 - 10^8$ and $\sim 10^{13} - 10^{14}$ M$_{\odot}$, respectively. We find $\sigma_{0} / m \gtrsim 0.1 \textrm{ cm}^{2}{\rm g}^{-1}$ is preferentially disfavored for both dwarfs and galaxy cluster haloes in comparison with observations, while $\sigma_{0} / m = 0.001 \textrm{ cm}^{2}{\rm g}^{-1}$ causes little perceptible difference from that of the CDM counterpart for most of the cross section's velocity dependence studied in this work. Our main result is that within the reasonable set of parameters the 2cDM model can successfully explain the observational trends seen in dwarf galaxy and galaxy cluster haloes and the model leaves us an open window for other possible alternative DM models.

We study various implementations of supernova feedback model and present the results of our `Osaka feedback model' using isolated galaxy simulations performed by the smoothed particle hydrodynamics (SPH) code {\small GADGET-3}. Our model is a modified version of Stinson et al.'s work, and we newly add the momentum kick for SN feedback rather than only thermal feedback. We incorporate the physical properties of SN remnants from the results of Chevalier and McKee \& Ostriker, such as the effective radius of SN bubble and the remnant life-time, in the form of Sedov-Taylor (ST)-like solutions with the effect of radiative cooling. Our model utilizes the local, physical parameters such as density and temperature of the ISM rather than galactic or halo properties to determine the galactic wind velocity or mass-loading factor. The Osaka model succeeds in self-regulating star formation, and naturally produces galactic outflow with variable velocities depending on the local environment and available SN energy as a function of time.An important addition to our previous work by Aoyama et al. is the implementation of the {\small CELib} chemistry library which allows us to deal with the time-dependent input of energy and metal yields for type Ia \& II supernovae (SNe) and asymptotic giant branch (AGB) stars. As initial tests of our model, we apply it to isolated galaxy simulations, and examine various galactic properties and compare with observational data including metal abundances.

Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly-formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low density region, and between more diffusive and less diffusive schemes in the high density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.

The multicomponent dark matter model with self-scattering and inter-conversions of species into one another is an alternative dark matter paradigm that is capable of resolving the long-standing problems of $\Lambda$CDM cosmology at small scales. In this paper, we have studied in detail the properties of dark matter halos with $M \sim 4-5 \times10^{11} M_{\odot}$ obtained in $N$-body cosmological simulations with the simplest two-component (2cDM) model. A large set of velocity-dependent cross-section prescriptions for elastic scattering and mass conversions, $\sigma_s(v)\propto v^{a_s}$ and $\sigma_c(v)\propto v^{a_c}$, has been explored and the results were compared with observational data. The results demonstrate that self-interactions with the cross-section per particle mass evaluated at $v=100$ km s$^{-1}$ being in the range of $0.01\lesssim \sigma_0/m\lesssim 1$ cm$^2$g$^{-1}$ robustly suppress central cusps, thus resolving the core-cusp problem. The core radii are controlled by the values of $\sigma_0/m$ and the DM cross-section's velocity-dependent power-law indices $(a_s,a_c)$, but are largely insensitive to the species' mass degeneracy. These values are in full agreement with those resolving the substructure and too-big-to-fail problems. We have also studied the evolution of halos in the 2cDM model with cosmic time.

Multicomponent dark matter with self-interactions, which allows for inter-conversions of species into one another, is a promising paradigm that is known to successfully and simultaneously resolve major problems of the conventional $\Lambda$CDM cosmology at galactic and sub-galactic scales. In this paper, we present $N$-body simulations of the simplest two-component (2cDM) model aimed at studying the distribution of dark matter halos with masses $M\lesssim10^{12}M_\odot$. In particular, we investigate how the maximum circular velocity function of the halos is affected by the velocity dependence of the self-interaction cross-sections, $\sigma(v)\propto v^a$, and compare them with available observational data. The results demonstrate that the 2cDM paradigm with the range of self-interaction cross-section per particle mass (evaluated at $v=100$ km s$^{-1}$) of $0.01\lesssim \sigma_0/m\lesssim 1 $ cm$^2$g$^{-1}$ and the mass degeneracy $\Delta m/m\sim 10^{-7}-10^{-8}$ is robustly resolving the substructure and too-big-to-fail problems by suppressing the substructure having small maximum circular velocities, $V_{\rm max}\lesssim100$ km s$^{-1}$. We also discuss the disagreement between the radial distribution of dwarfs in a host halo observed in the Local Group and simulated with CDM. This can be considered as one more small-scale problem of CDM. We demonstrate that such a disagreement is alleviated in 2cDM. Finally, the computed matter power-spectra of the 2cDM structure indicate the model's consistency with the existing Ly-$\alpha$ forest constraints.

We perform smoothed particle hydrodynamics (SPH) simulations of an isolated galaxy with a new treatment for dust formation and destruction. To this aim, we treat dust and metal production self-consistently with star formation and supernova feedback. For dust, we consider a simplified model of grain size distribution by representing the entire range of grain sizes with large and small grains. We include dust production in stellar ejecta, dust destruction by supernova (SN) shocks, grain growth by accretion and coagulation, and grain disruption by shattering. We find that the assumption of fixed dust-to-metal mass ratio becomes no longer valid when the galaxy is older than 0.2 Gyr, at which point the grain growth by accretion starts to contribute to the nonlinear rise of dust-to-gas ratio. As expected in our previous one-zone model, shattering triggers grain growth by accretion since it increases the total surface area of grains. Coagulation becomes significant when the galaxy age is greater than $\sim$ 1 Gyr: at this epoch the abundance of small grains becomes high enough to raise the coagulation rate of small grains. We further compare the radial profiles of dust-to-gas ratio $(\mathcal{D})$ and dust-to-metal ratio $(\mathcal{D}/Z)$ (i.e., depletion) at various ages with observational data. We find that our simulations broadly reproduce the radial gradients of dust-to-gas ratio and depletion. In the early epoch ($\lesssim 0.3$ Gyr), the radial gradient of $\mathcal{D}$ follows the metallicity gradient with $\mathcal{D}/Z$ determined by the dust condensation efficiency in stellar ejecta, while the $\mathcal{D}$ gradient is steeper than the $Z$ gradient at the later epochs because of grain growth by accretion. The framework developed in this paper is applicable to any SPH-based galaxy evolution simulations including cosmological ones.

We introduce the AGORA project, a comprehensive numerical study of well-resolved galaxies within the LCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of 8 galaxies with halo masses M_vir ~= 1e10, 1e11, 1e12, and 1e13 Msun at z=0 and two different ("violent" and "quiescent") assembly histories. The numerical techniques and implementations used in this project include the smoothed particle hydrodynamics codes GADGET and GASOLINE, and the adaptive mesh refinement codes ART, ENZO, and RAMSES. The codes will share common initial conditions and common astrophysics packages including UV background, metal-dependent radiative cooling, metal and energy yields of supernovae, and stellar initial mass function. These are described in detail in the present paper. Subgrid star formation and feedback prescriptions will be tuned to provide a realistic interstellar and circumgalactic medium using a non-cosmological disk galaxy simulation. Cosmological runs will be systematically compared with each other using a common analysis toolkit, and validated against observations to verify that the solutions are robust - i.e., that the astrophysical assumptions are responsible for any success, rather than artifacts of particular implementations. The goals of the AGORA project are, broadly speaking, to raise the realism and predictive power of galaxy simulations and the understanding of the feedback processes that regulate galaxy "metabolism." The proof-of-concept dark matter-only test of the formation of a galactic halo with a z=0 mass of M_vir ~= 1.7e11 Msun by 9 different versions of the participating codes is also presented to validate the infrastructure of the project.