Search

We present a suite of baryonic cosmological zoom-in simulations of self-interacting dark matter (SIDM) haloes within the 'Feedback In Realistic Environment' (FIRE) project. The three simulated haloes have virial masses of ∼ 1012 M⊙ at z = 0, and we study velocity-independent self-interaction cross sections of 1 and 10 cm2, g-1. We study star formation rates and the shape of dark matter density profiles of the parent haloes in both cold dark matter (CDM) and SIDM models. Galaxies formed in the SIDM haloes have higher star formation rates at z ≤ 1, resulting in more massive galaxies compared to the CDM simulations. While both CDM and SIDM simulations show diverse shape of the dark matter density profiles, the SIDM haloes can reach higher and more steep central densities within few kpcs compared to the CDM haloes. We identify a correlation between the build-up of the stars within the half-mass radii of the galaxies and the growth in the central dark matter densities. The thermalization process in the SIDM haloes is enhanced in the presence of a dense stellar component. Hence, SIDM haloes with highly concentrated baryonic profiles are predicted to have higher central dark matter densities than the CDM haloes. Overall, the SIDM haloes are more responsive to the presence of a massive baryonic distribution than their CDM counterparts.

We perform cosmological hydrodynamical simulations to study the formation of proto-globular cluster candidates in progenitors of present-day dwarf galaxies |$(M_{\rm vir} \approx 10^{10}\, {\rm M}_\odot$| at z = 0) as part of the 'Feedback in Realistic Environment' (FIRE) project. Compact (r 1/2 < 30 pc), relatively massive (0.5 × 105 ≲ M ⋆/M⊙ ≲ 5 × 105), self-bound stellar clusters form at 11 ≳ z ≳ 5 in progenitors with |$M_{\rm vir} \approx 10^9\, \mathrm{M}_{\odot }$|⁠. Cluster formation is triggered when at least |$10^7\, \mathrm{M}_{\odot }$| of dense, turbulent gas reaches |$\Sigma _{\rm gas} \approx 10^4\, {\rm M}_\odot \, {\rm pc}^{-2}$| as a result of the compressive effects of supernova feedback or from cloud–cloud collisions. The clusters can survive for |$2-3\, {\rm Gyr}$|⁠ ; absent numerical effects, they could possibly survive substantially longer, perhaps to z = 0. The longest lived clusters are those that form at significant distance – several hundreds of pc – from their host galaxy. We therefore predict that globular clusters forming in progenitors of present-day dwarf galaxies will be offset from any pre-existing stars within their host dark matter haloes as opposed to deeply embedded within a well-defined galaxy. Properties of the nascent clusters are consistent with observations of some of the faintest and most compact high-redshift sources in Hubble Space Telescope lensing fields and are at the edge of what will be detectable as point sources in deep imaging of non-lensed fields with JWST. By contrast, the star clusters' host galaxies will remain undetectable. [ABSTRACT FROM AUTHOR]

Increasingly, uncertainties in predictions from galaxy formation simulations (at sub-Milky Way masses) are dominated by uncertainties in stellar evolution inputs. In this paper, we present the full set of updates from the Feedback In Realistic Environment (FIRE)-2 version of the FIRE project code, to the next version, FIRE-3. While the transition from FIRE-1 to FIRE-2 focused on improving numerical methods, here we update the stellar evolution tracks used to determine stellar feedback inputs, e.g. stellar mass-loss (O/B and AGB), spectra (luminosities and ionization rates), and supernova rates (core-collapse and Ia), as well as detailed mass-dependent yields. We also update the low-temperature cooling and chemistry, to enable improved accuracy at |$T \lesssim 10^{4}\,$| K and densities |$n\gg 1\, {\rm cm^{-3}}$|⁠ , and the meta-galactic ionizing background. All of these synthesize newer empirical constraints on these quantities and updated stellar evolution and yield models from a number of groups, addressing different aspects of stellar evolution. To make the updated models as accessible as possible, we provide fitting functions for all of the relevant updated tracks, yields, etc, in a form specifically designed so they can be directly 'plugged in' to existing galaxy formation simulations. We also summarize the default FIRE-3 implementations of 'optional' physics, including spectrally resolved cosmic rays and supermassive black hole growth and feedback. [ABSTRACT FROM AUTHOR]

Self-interacting dark matter (SIDM) models offer one way to reconcile inconsistencies between observations and predictions from collisionless cold dark matter (CDM) models on dwarf-galaxy scales. In order to incorporate the effects of both baryonic and SIDM interactions, we study a suite of cosmological-baryonic simulations of Milky-Way (MW)-mass galaxies from the Feedback in Realistic Environments (FIRE-2) project where we vary the SIDM self-interaction cross-section σ/ m. We compare the shape of the main dark matter (DM) halo at redshift z  = 0 predicted by SIDM simulations (at σ/ m  = 0.1, 1, and 10 cm2 g−1) with CDM simulations using the same initial conditions. In the presence of baryonic feedback effects, we find that SIDM models do not produce the large differences in the inner structure of MW-mass galaxies predicted by SIDM-only models. However, we do find that the radius where the shape of the total mass distribution begins to differ from that of the stellar mass distribution is dependent on σ/ m. This transition could potentially be used to set limits on the SIDM cross-section in the MW. [ABSTRACT FROM AUTHOR]

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-1, whose diversity has been found to be greater than predicted by the Λ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 a 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 haloes with compact baryonic distributions versus underdense haloes with extended baryonic discs. 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, emphasizing the need for a better understanding of the frequency of extreme outliers present in current observational samples. [ABSTRACT FROM AUTHOR]