skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Spatially resolved gas-phase metallicity in FIRE-2 dwarfs: late-time evolution of metallicity relations in simulations with feedback and mergers
ABSTRACT We present an analysis of spatially resolved gas-phase metallicity relations in five dwarf galaxies ($$\rm \mathit{M}_{halo} \approx 10^{11}\, {\rm M}_\odot$$, $$\rm \mathit{M}_\star \approx 10^{8.8}{-}10^{9.6}\, {\rm M}_\odot$$) from the FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulation suite, which include an explicit model for sub-grid turbulent mixing of metals in gas, near z ≈ 0, over a period of 1.4 Gyr, and compare our findings with observations. While these dwarf galaxies represent a diverse sample, we find that all simulated galaxies match the observed mass–metallicity (MZR) and mass–metallicity gradient (MZGR) relations. We note that in all five galaxies, the metallicities are effectively identical between phases of the interstellar medium (ISM), with 95 $${{\ \rm per\ cent}}$$ of the gas being within ±0.1 dex between the cold and dense gas (T < 500 K and nH > 1 cm−3), ionized gas (near the H αT ≈ 104 K ridge-line), and nebular regions (ionized gas where the 10 Myr-averaged star formation rate is non-zero). We find that most of the scatter in relative metallicity between cold dense gas and ionized gas/nebular regions can be attributed to either local starburst events or metal-poor inflows. We also note the presence of a major merger in one of our galaxies, m11e, with a substantial impact on the metallicity distribution in the spatially resolved map, showing two strong metallicity peaks and triggering a starburst in the main galaxy.  more » « less
Award ID(s):
2107772 2108318 1715847 1713353
PAR ID:
10432524
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
515
Issue:
3
ISSN:
0035-8711
Page Range / eLocation ID:
p. 3555-3576
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT 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 (r1/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. 
    more » « less
  2. 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. 
    more » « less
  3. ABSTRACT We present spatially resolved kinematic measurements of stellar and ionized gas components of dwarf galaxies in the stellar mass range $$10^{8.5}\!-\!10^{10} \, \mathrm{M}_{\odot }$$, selected from Sloan Digital Sky Survey DR7 and DR8 and followed up with Keck/Low-Resolution Imaging Spectrometer spectroscopy. We study the potential effects of active galactic nuclei (AGNs) on Galaxy-wide gas kinematics by comparing rotation curves of 26 Galaxies containing AGNs, and 19 control Galaxies with no optical or infrared signs of AGNs. We find a strong association between AGN activity and disturbed gas kinematics in the host Galaxies. While star-forming Galaxies in this sample tend to have orderly gas discs that co-rotate with the stars, 73 per cent of the AGNs have disturbed gas. We find that 5 out of 45 Galaxies have gaseous components in counter-rotation with their stars, and all Galaxies exhibiting counter-rotation contain AGNs. Six out of seven isolated Galaxies with disturbed ionized gas host AGNs. At least three AGNs fall clearly below the stellar–halo mass relation, which could be interpreted as evidence for ongoing star formation suppression. Taken together, these results provide new evidence supporting the ability of AGN to influence gas kinematics and suppress star formation in dwarf galaxies. This further demonstrates the importance of including AGN as a feedback mechanism in galaxy formation models in the low-mass regime. 
    more » « less
  4. ABSTRACT We revisit the question of ‘hot mode’ versus ‘cold mode’ accretion on to galaxies using steady-state cooling flow solutions and idealized 3D hydrodynamic simulations. We demonstrate that for the hot accretion mode to exist, the cooling time is required to be longer than the free-fall time near the radius where the gas is rotationally supported, Rcirc, i.e. the existence of the hot mode depends on physical conditions at the galaxy scale rather than on physical conditions at the halo scale. When allowing for the depletion of the halo baryon fraction relative to the cosmic mean, the longer cooling times imply that a virialized gaseous halo may form in halo masses below the threshold of $$\sim 10^{12}\, {\rm M_{\odot }}$$ derived for baryon-complete haloes. We show that for any halo mass there is a maximum accretion rate for which the gas is virialized throughout the halo and can accrete via the hot mode of $${\dot{M}}_{\rm crit}\approx 0.7(v_{\rm c}/100\, \rm km\ s^{-1})^{5.4}(R_{\rm circ}/10\, {\rm kpc})(Z/\, {\rm Z_{\odot }})^{-0.9}\, {\rm M_{\odot }}\, {\rm yr}^{-1}$$, where Z and vc are the metallicity and circular velocity measured at Rcirc. For accretion rates $$\gtrsim {\dot{M}}_{\rm crit}$$ the volume-filling gas phase can in principle be ‘transonic’ – virialized in the outer halo but cool and free-falling near the galaxy. We compare $${\dot{M}}_{\rm crit}$$ to the average star formation rate (SFR) in haloes at 0 < z < 10 implied by the stellar-mass–halo-mass relation. For a plausible metallicity evolution with redshift, we find that $${\rm SFR}\lesssim {\dot{M}}_{\rm crit}$$ at most masses and redshifts, suggesting that the SFR of galaxies could be primarily sustained by the hot mode in halo masses well below the classic threshold of $$\sim 10^{12}\, {\rm M_{\odot }}$$. 
    more » « less
  5. ABSTRACT We analyse the rest-optical emission-line ratios of z ∼ 1.5 galaxies drawn from the Multi-Object Spectrometer for Infra-Red Exploration Deep Evolution Field (MOSDEF) survey. Using composite spectra, we investigate the mass–metallicity relation (MZR) at z ∼ 1.5 and measure its evolution to z = 0. When using gas-phase metallicities based on the N2 line ratio, we find that the MZR evolution from z ∼ 1.5 to z = 0 depends on stellar mass, evolving by $$\Delta \rm log(\rm O/H) \sim 0.25$$ dex at M*< $$10^{9.75}\, \mathrm{M}_{\odot }$$ down to $$\Delta \rm log(\rm O/H) \sim 0.05$$ at M* ≳ $$10^{10.5}\, \mathrm{M}_{\odot }$$. In contrast, the O3N2-based MZR shows a constant offset of $$\Delta \rm log(\rm O/H) \sim 0.30$$ across all masses, consistent with previous MOSDEF results based on independent metallicity indicators, and suggesting that O3N2 provides a more robust metallicity calibration for our z ∼ 1.5 sample. We investigated the secondary dependence of the MZR on star formation rate (SFR) by measuring correlated scatter about the mean M*-specific SFR and M*−$$\log (\rm O3N2)$$ relations. We find an anticorrelation between $$\log (\rm O/H)$$ and sSFR offsets, indicating the presence of a M*−SFR−Z relation, though with limited significance. Additionally, we find that our z ∼ 1.5 stacks lie along the z = 0 metallicity sequence at fixed μ = log (M*/M⊙) − 0.6 × $$\log (\rm SFR / M_{\odot } \, yr^{-1})$$ suggesting that the z ∼ 1.5 stacks can be described by the z = 0 fundamental metallicity relation (FMR). However, using different calibrations can shift the calculated metallicities off of the local FMR, indicating that appropriate calibrations are essential for understanding metallicity evolution with redshift. Finally, understanding how [N ii]/H α scales with galaxy properties is crucial to accurately describe the effects of blended [N ii] and H α on redshift and H α fiux measurements in future large surveys utilizing low-resolution spectra such as with Euclid and the Roman Space Telescope. 
    more » « less