Search for: All records

ABSTRACT Self-interacting dark matter (SIDM) models have received great attention over the past decade as solutions to the small-scale puzzles of astrophysics. Though there are different implementations of dark matter (DM) self-interactions in N-body codes of structure formation, there has not been a systematic study to compare the predictions of these different implementations. We investigate the implementation of dark matter self-interactions in two simulation codes:gizmo and arepo. We begin with identical initial conditions for an isolated 1010 M⊙ dark matter halo and investigate the evolution of the density and velocity dispersion profiles in gizmo and arepo for SIDM cross-section over mass of 1, 5, and 50 $\rm cm^2\, g^{-1}$. Our tests are restricted to the core expansion phase, where the core density decreases and core radius increases with time. We find better than 30 per cent agreement between the codes for the density profile in this phase of evolution, with the agreement improving at higher resolution. We find that varying code-specific SIDM parameters changes the central halo density by less than 10 per cent outside of the convergence radius. We argue that SIDM core formation is robust across the two different schemes and conclude that these codes can reliably differentiate between cross-sections of 1, 5, and 50 $\rm cm^2\,more »g^{-1}$, but finer distinctions would require further investigation.« less

ABSTRACT We explore the origin of stellar metallicity gradients in simulated and observed dwarf galaxies. We use FIRE-2 cosmological baryonic zoom-in simulations of 26 isolated galaxies as well as existing observational data for 10 Local Group dwarf galaxies. Our simulated galaxies have stellar masses between 105.5 and 108.6 M⊙. Whilst gas-phase metallicty gradients are generally weak in our simulated galaxies, we find that stellar metallicity gradients are common, with central regions tending to be more metal-rich than the outer parts. The strength of the gradient is correlated with galaxy-wide median stellar age, such that galaxies with younger stellar populations have flatter gradients. Stellar metallicty gradients are set by two competing processes: (1) the steady ‘puffing’ of old, metal-poor stars by feedback-driven potential fluctuations and (2) the accretion of extended, metal-rich gas at late times, which fuels late-time metal-rich star formation. If recent star formation dominates, then extended, metal-rich star formation washes out pre-existing gradients from the ‘puffing’ process. We use published results from ten Local Group dwarf galaxies to show that a similar relationship between age and stellar metallicity-gradient strength exists among real dwarfs. This suggests that observed stellar metallicity gradients may be driven largely by the baryon/feedback cycle rather thanmore »by external environmental effects.« less

ABSTRACT Understanding the rate at which stars form is central to studies of galaxy formation. Observationally, the star formation rates (SFRs) of galaxies are measured using the luminosity in different frequency bands, often under the assumption of a time-steady SFR in the recent past. We use star formation histories (SFHs) extracted from cosmological simulations of star-forming galaxies from the FIRE project to analyse the time-scales to which the H α and far-ultraviolet (FUV) continuum SFR indicators are sensitive. In these simulations, the SFRs are highly time variable for all galaxies at high redshift, and continue to be bursty to z = 0 in dwarf galaxies. When FIRE SFHs are partitioned into their bursty and time-steady phases, the best-fitting FUV time-scale fluctuates from its ∼10 Myr value when the SFR is time-steady to ≳100 Myr immediately following particularly extreme bursts of star formation during the bursty phase. On the other hand, the best-fitting averaging time-scale for H α is generally insensitive to the SFR variability in the FIRE simulations and remains ∼5 Myr at all times. These time-scales are shorter than the 100 and 10 Myr time-scales sometimes assumed in the literature for FUV and H α, respectively, because while the FUV emission persists for stellar populations oldermore »than 100 Myr, the time-dependent luminosities are strongly dominated by younger stars. Our results confirm that the ratio of SFRs inferred using H α versus FUV can be used to probe the burstiness of star formation in galaxies.« less

We use FIRE-2 zoom cosmological simulations of Milky Way size Galaxy haloes to calculate astrophysical J-factors for dark matter annihilation and indirect detection studies. In addition to velocity-independent (s-wave) annihilation cross-sections 〈σv〉, we also calculate effective J-factors for velocity-dependent models, where the annihilation cross-section is either p-wave (∝ v2/c2) or d-wave (∝ v4/c4). We use 12 pairs of simulations, each run with dark matter-only (DMO) physics and FIRE-2 physics. We observe FIRE runs produce central dark matter velocity dispersions that are systematically larger than in DMO runs by factors of ∼2.5–4. They also have a larger range of central (∼400 pc) dark matter densities than the DMO runs (ρFIRE/ρDMO ≃ 0.5–3) owing to the competing effects of baryonic contraction and feedback. At 3 deg from the Galactic Centre, FIRE J-factors are 3–60 (p-wave) and 10–500 (d-wave) times higher than in the DMO runs. The change in s-wave signal at 3 deg is more modest and can be higher or lower (∼0.3–7), though the shape of the emission profile is flatter (less peaked towards the Galactic Centre) and more circular on the sky in FIRE runs. Our results for s-wave are broadly consistent with the range of assumptions in most indirect detection studies. We observemore »p-wave J-factors that are significantly enhanced compared to most past estimates. We find that thermal models with p-wave annihilation may be within range of detection in the near future.

« less