An official website of the United States government
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
Spiral arms are metal freeways: azimuthal gas-phase metallicity variations in flocculent discs in the FIRE-2 cosmological zoom-in simulations
ABSTRACT We examine the azimuthal variations in gas-phase metallicity profiles in simulated Milky Way-mass disc galaxies from the Feedback in Realistic Environments (FIRE-2) cosmological zoom-in simulation suite, which includes a sub-grid turbulent metal mixing model. We produce spatially resolved maps of the discs at z ≈ 0 with pixel sizes ranging from 250 to 750 pc, analogous to modern integral field unit galaxy surveys, mapping the gas-phase metallicities in both the cold and dense gas and the ionized gas correlated with H ii regions. We report that the spiral arms alternate in a pattern of metal rich and metal poor relative to the median metallicity of the order of ≲0.1 dex, appearing generally in this sample of flocculent spirals. The pattern persists even in a simulation with different strengths of metal mixing, indicating that the pattern emerges from physics above the sub-grid scale. Local enrichment does not appear to be the dominant source of the azimuthal metallicity variations at z ≈ 0: there is no correlation with local star formation on these spatial scales. Rather, the arms are moving radially inwards and outwards relative to each other, carrying their local metallicity gradients with them radially before mixing into the larger-scale interstellar medium. We propose that the arms act as freeways channeling relatively metal poor gas radially inwards, and relatively enriched gas radially outwards.
more »
« less
- PAR ID:
- 10403164
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 521
- Issue:
- 3
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 3708-3726
- Size(s):
- p. 3708-3726
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The metallicity of galaxies, and its variation with galactocentric radius, provides key insights into the formation histories of galaxies and the physical processes driving their evolution. In this work, we analyze the radial metallicity gradients of star-forming galaxies in the EAGLE, Illustris, IllustrisTNG, and SIMBA cosmological simulations across broad mass (108.0M⊙≤M⋆ ≲ 1012.0M⊙) and redshift (0 ≤z≤ 8) ranges. We find that all simulations predict strong negative (i.e., radially decreasing) metallicity gradients at early cosmic times, likely due to their similar treatments of relatively smooth stellar feedback not providing sufficient mixing to quickly flatten gradients. The strongest redshift evolution occurs in galaxies with stellar masses of 1010.0–1011.0M⊙, while galaxies with stellar mass < 1010M⊙and >1011M⊙exhibit weaker redshift evolution. Our result of negative gradients at high redshift contrast with the many positive and flat gradients in the 1 < z < 4 observational literature. Atz > 6, the negative gradients observed with JWST and the Atacama Large Millimeter/submillimeter Array are flatter than those in simulations, albeit with closer agreement than at lower redshift. Overall, we suggest that these smooth stellar feedback galaxy simulations may not sufficiently mix their metal content radially, and that either stronger stellar feedback or additional subgrid turbulent metal diffusion models may be required to better reproduce observed metallicity gradients.more » « less
-
null (Ed.)ABSTRACT We use FIRE-2 simulations to examine 3D variations of gas-phase elemental abundances of [O/H], [Fe/H], and [N/H] in 11 MW and M31-mass galaxies across their formation histories at z ≤ 1.5 ($$t_{\rm lookback} \le 9.4 \, \rm {Gyr}$$), motivated by characterizing the initial conditions of stars for chemical tagging. Gas within $$1 \, \rm {kpc}$$ of the disc mid-plane is vertically homogeneous to $$\lesssim 0.008 \, \rm {dex}$$ at all z ≤ 1.5. We find negative radial gradients (metallicity decreases with galactocentric radius) at all times, which steepen over time from $$\approx \! -0.01 \, \rm {dex}\, \rm {kpc}^{-1}$$ at z = 1 ($$t_{\rm lookback} = 7.8 \, \rm {Gyr}$$) to $$\approx \! -0.03 \, \rm {dex}\, \rm {kpc}^{-1}$$ at z = 0, and which broadly agree with observations of the MW, M31, and nearby MW/M31-mass galaxies. Azimuthal variations at fixed radius are typically $$0.14 \, \rm {dex}$$ at z = 1, reducing to $$0.05 \, \rm {dex}$$ at z = 0. Thus, over time radial gradients become steeper while azimuthal variations become weaker (more homogeneous). As a result, azimuthal variations were larger than radial variations at z ≳ 0.8 ($$t_{\rm lookback} \gtrsim 6.9 \, \rm {Gyr}$$). Furthermore, elemental abundances are measurably homogeneous (to ≲0.05 dex) across a radial range of $$\Delta R \approx 3.5 \, \rm {kpc}$$ at z ≳ 1 and $$\Delta R \approx 1.7 \, \rm {kpc}$$ at z = 0. We also measure full distributions of elemental abundances, finding typically negatively skewed normal distributions at z ≳ 1 that evolve to typically Gaussian distributions by z = 0. Our results on gas abundances inform the initial conditions for stars, including the spatial and temporal scales for applying chemical tagging to understand stellar birth in the MW.more » « less
-
null (Ed.)ABSTRACT We present the measured gas-phase metal column densities in 155 sub-damped Ly α systems (subDLAs) with the aim to investigate the contribution of subDLAs to the chemical evolution of the Universe. The sample was identified within the absorber-blind XQ-100 quasar spectroscopic survey over the redshift range 2.4 ≤ zabs ≤ 4.3. Using all available column densities of the ionic species investigated (mainly C iv, Si ii, Mg ii, Si iv, Al ii, Fe ii, C ii, and O i; in order of decreasing detection frequency), we estimate the ionization-corrected gas-phase metallicity of each system using Markov chain Monte Carlo techniques to explore a large grid of cloudy ionization models. Without accounting for ionization and dust depletion effects, we find that the H i-weighted gas-phase metallicity evolution of subDLAs is consistent with damped Ly α systems (DLAs). When ionization corrections are included, subDLAs are systematically more metal poor than DLAs (between ≈0.5σ and ≈3σ significance) by up to ≈1.0 dex over the redshift range 3 ≤ zabs ≤ 4.3. The correlation of gas phase [Si/Fe] with metallicity in subDLAs appears to be consistent with that of DLAs, suggesting that the two classes of absorbers have a similar relative dust depletion pattern. As previously seen for Lyman limit systems, the gas phase [C/O] in subDLAs remains constantly solar for all metallicities indicating that both subDLAs and Lyman limit systems could trace carbon-rich ejecta, potentially in circumgalactic environments.more » « less
-
ABSTRACT We present predictions for the evolution of the galaxy dust-to-gas ratio (DGR) and dust-to-metal ratio (DTM) from z = 0 → 6, using a model for the production, growth, and destruction of dust grains implemented into the simba cosmological hydrodynamic galaxy formation simulation. In our model, dust forms in stellar ejecta, grows by the accretion of metals, and is destroyed by thermal sputtering and supernovae. Our simulation reproduces the observed dust mass function at z = 0, but modestly underpredicts the mass function by ∼×3 at z ∼ 1–2. The z = 0 DGR versus metallicity relationship shows a tight positive correlation for star-forming galaxies, while it is uncorrelated for quenched systems. There is little evolution in the DGR–metallicity relationship between z = 0 and 6. We use machine learning techniques to search for the galaxy physical properties that best correlate with the DGR and DTM. We find that the DGR is primarily correlated with the gas-phase metallicity, though correlations with the depletion time-scale, stellar mass, and gas fraction are non-negligible. We provide a crude fitting relationship for DGR and DTM versus the gas-phase metallicity, along with a public code package that estimates the DGR and DTM given a set of galaxy physical properties.more » « less
