Search for: All records

Dwarf galaxies located in extremely underdense cosmic voids are excellent test beds for disentangling the effects of large-scale environments on galaxy formation and evolution. We present the first results of the Dwarfs in Void Environments Survey, which has obtained integral field spectroscopy for low-mass galaxies (M_⋆= 10⁷–10⁹M_⊙) located inside (N= 21) and outside (N= 9) cosmic voids using the Keck Cosmic Web Imager. Using measurements of stellar line-of-sight rotational velocityv_rotand velocity dispersionσ_⋆, we test the tidal stirring hypothesis, which posits that dwarf spheroidal galaxies are formed through tidal interactions with more massive host galaxies. We measure low values ofv_rot/σ_⋆≲ 2 for our sample of isolated dwarf galaxies, and we find no trend betweenv_rot/σ_⋆and the distance from a massive galaxy$d_{L^{\star }}$out to$d_{L^{\star }}\sim 10$Mpc. These suggest that dwarf galaxies can become dispersion-supported, “puffy” systems even in the absence of environmental effects like tidal interactions. We also find indications of an upward trend betweenv_rot/σ_⋆and galaxy stellar mass, perhaps implying that stellar disk formation depends on mass rather than environment. Although some of our conclusions may be slightly modified by systematic effects, our main result still holds: that isolated low-mass galaxies may form and remain as puffy systems rather than the dynamically cold disks predicted by classical galaxy formation theory.

 

Abstract Type Ia supernovae are critical for feedback and elemental enrichment in galaxies. Recent surveys like the All-Sky Automated Survey for Supernova (ASAS-SN) and the Dark Energy Survey (DES) find that the specific supernova Ia rate at z ∼ 0 may be ≲ 20 − 50 × higher in lower-mass galaxies than at Milky Way-mass. Independently, observations show that the close-binary fraction of solar-type Milky Way stars is higher at lower metallicity. Motivated by these observations, we use the FIRE-2 cosmological zoom-in simulations to explore the impact of metallicity-dependent rate models on galaxies of $M_* \sim 10^7\, \rm {M}_{\odot }-10^{11}\, \rm {M}_{\odot }$. First, we benchmark our simulated star-formation histories (SFHs) against observations, and show that assumed stellar mass functions play a major role in determining the degree of tension between observations and metallicity-independent rate models, potentially causing ASAS-SN and DES observations to agree more than might appear. Models in which the supernova Ia rate increases with decreasing metallicity ($\propto Z^{-0.5 \; \rm {to} \; -1}$) provide significantly better agreement with observations. Encouragingly, these rate increases (≳ 10 × in low-mass galaxies) do not significantly impact galaxy masses and morphologies, which remain largely unaffected except for our most extreme models. We explore implications for both [Fe/H] and [$\alpha /\rm {Fe}$] enrichment; metallicity-dependent rate models can improve agreement with the observed stellar mass-metallicity relations in low-mass galaxies. Our results demonstrate that a range of metallicity-dependent rate models are viable for galaxy formation and motivate future work. 

As they grow, galaxies can transition from irregular/spheroidal with ‘bursty’ star formation histories (SFHs), to discy with smooth SFHs. But even in simulations, the direct physical cause of such transitions remains unclear. We therefore explore this in a large suite of numerical experiments re-running portions of cosmological simulations with widely varied physics, further validated with existing FIRE simulations. We show that gas supply, cooling/thermodynamics, star formation model, Toomre scale, galaxy dynamical times, and feedback properties do not have a direct causal effect on these transitions. Rather, both the formation of discs and cessation of bursty star formation are driven by the gravitational potential, but in different ways. Disc formation is promoted when the mass profile becomes sufficiently centrally concentrated in shape (relative to circularization radii): we show that this provides a well-defined dynamical centre, ceases to support the global ‘breathing modes’ that can persist indefinitely in less-concentrated profiles and efficiently destroy discs, promotes orbit mixing to form a coherent angular momentum, and stabilizes the disc. Smooth SF is promoted by the potential or escape velocity Vesc (not circular velocity Vc) becoming sufficiently large at the radii of star formation that cool, mass-loaded (momentum-conserving) outflows are trapped/confined near the galaxy, as opposed to escaping after bursts. We discuss the detailed physics, how these conditions arise in cosmological contexts, their relation to other correlated phenomena (e.g. inner halo virialization, vertical disc ‘settling’), and observations.

 

Abstract We describe a public data release of the FIRE-2 cosmological zoom-in simulations of galaxy formation (available at http://flathub.flatironinstitute.org/fire ) from the Feedback In Realistic Environments (FIRE) project. FIRE-2 simulations achieve parsec-scale resolution to explicitly model the multiphase interstellar medium while implementing direct models for stellar evolution and feedback, including stellar winds, core-collapse and Type Ia supernovae, radiation pressure, photoionization, and photoelectric heating. We release complete snapshots from three suites of simulations. The first comprises 20 simulations that zoom in on 14 Milky Way (MW)–mass galaxies, five SMC/LMC-mass galaxies, and four lower-mass galaxies including one ultrafaint; we release 39 snapshots across z = 0–10. The second comprises four massive galaxies, with 19 snapshots across z = 1–10. Finally, a high-redshift suite comprises 22 simulations, with 11 snapshots across z = 5–10. Each simulation also includes dozens of resolved lower-mass (satellite) galaxies in its zoom-in region. Snapshots include all stored properties for all dark matter, gas, and star particles, including 11 elemental abundances for stars and gas, and formation times (ages) of star particles. We also release accompanying (sub)halo catalogs, which include galaxy properties and member star particles. For the simulations to z = 0, including all MW-mass galaxies, we release the formation coordinates and an “ex situ” flag for all star particles, pointers to track particles across snapshots, catalogs of stellar streams, and multipole basis expansions for the halo mass distributions. We describe publicly available python packages for reading and analyzing these simulations. 

ABSTRACT 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 We explore the effect of including progenitor mass- and metallicity-dependent yields, supernova rates and energetics on variations in elemental abundance ratios (particularly [α/Fe]) in dwarf galaxies. To understand how the scatter and overall trends in [α/Fe] are affected by including variable metal yields from a discretely sampled initial mass function, we run FIRE simulations of a dwarf galaxy (M⋆(z = 0$) \sim 10^6\rm \, M_{\odot })$ using nucleosynthetic yields from the NuGrid data base that depend on the stellar progenitor mass and metallicity. While NuGrid exhibits lower aggregate α-element production than default FIRE yields, we find that its explicit mass dependence, even when including turbulent metal diffusion, substantially widens the intrinsic scatter in the simulated [Fe/H]-[α/Fe] – a phenomenon visible in some observations of dwarf galaxies. 

Extended, old, and round stellar halos appear to be ubiquitous around high-mass dwarf galaxies (10^8.5<M_⋆/M_⊙< 10^9.6) in the observed universe. However, it is unlikely that these dwarfs have undergone a sufficient number of minor mergers to form stellar halos that are composed of predominantly accreted stars. Here, we demonstrate that FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulations are capable of producing dwarf galaxies with realistic structures, including both a thick disk and round stellar halo. Crucially, these stellar halos are formed in situ, largely via the outward migration of disk stars. However, there also exists a large population of “nondisky” dwarfs in FIRE-2 that lack a well-defined disk/halo and do not resemble the observed dwarf population. These nondisky dwarfs tend to be either more gas-poor or to have burstier recent star formation histories than the disky dwarfs, suggesting that star formation feedback may be preventing disk formation. Both classes of dwarfs underscore the power of a galaxy’s intrinsic shape—which is a direct quantification of the distribution of the galaxy’s stellar content—to interrogate the feedback implementation in simulated galaxies.

 

ABSTRACT We analyse the cold dark matter density profiles of 54 galaxy haloes simulated with Feedback In Realistic Environments (FIRE)-2 galaxy formation physics, each resolved within $0.5{{\ \rm per\ cent}}$ of the halo virial radius. These haloes contain galaxies with masses that range from ultrafaint dwarfs ($M_\star \simeq 10^{4.5}\, \mathrm{M}_{\odot }$) to the largest spirals ($M_\star \simeq 10^{11}\, \mathrm{M}_{\odot }$) and have density profiles that are both cored and cuspy. We characterize our results using a new, analytic density profile that extends the standard two-parameter Einasto form to allow for a pronounced constant density core in the resolved innermost radius. With one additional core-radius parameter, rc, this three-parameter core-Einasto profile is able to characterize our feedback-impacted dark matter haloes more accurately than other three-parameter profiles proposed in the literature. To enable comparisons with observations, we provide fitting functions for rc and other profile parameters as a function of both M⋆ and M⋆/Mhalo. In agreement with past studies, we find that dark matter core formation is most efficient at the characteristic stellar-to-halo mass ratio M⋆/Mhalo ≃ 5 × 10−3, or $M_{\star } \sim 10^9 \, \mathrm{M}_{\odot }$, with cores that are roughly the size of the galaxy half-light radius, rc ≃ 1−5 kpc. Furthermore, we find no evidence for core formation at radii $\gtrsim 100\ \rm pc$ in galaxies with M⋆/Mhalo < 5 × 10−4 or $M_\star \lesssim 10^6 \, \mathrm{M}_{\odot }$. For Milky Way-size galaxies, baryonic contraction often makes haloes significantly more concentrated and dense at the stellar half-light radius than DMO runs. However, even at the Milky Way scale, FIRE-2 galaxy formation still produces small dark matter cores of ≃ 0.5−2 kpc in size. Recent evidence for a ∼2 kpc core in the Milky Way’s dark matter halo is consistent with this expectation.