Search for: All records

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. Wemore »propose that the arms act as freeways channeling relatively metal poor gas radially inwards, and relatively enriched gas radially outwards.

« less

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’more »physics, including spectrally resolved cosmic rays and supermassive black hole growth and feedback.« less

We characterize the 3D spatial variations of [Fe/H], [Mg/H], and [Mg/Fe] in stars at the time of their formation, across 11 simulated Milky Way (MW)- and M31-mass galaxies in the FIRE-2 simulations, to inform initial conditions for chemical tagging. The overall scatter in [Fe/H] within a galaxy decreased with time until $\approx 7 \, \rm {Gyr}$ ago, after which it increased to today: this arises from a competition between a reduction of azimuthal scatter and a steepening of the radial gradient in abundance over time. The radial gradient is generally negative, and it steepened over time from an initially flat gradient $\gtrsim 12 \, \rm {Gyr}$ ago. The strength of the present-day abundance gradient does not correlate with when the disc ‘settled’; instead, it best correlates with the radial velocity dispersion within the galaxy. The strength of azimuthal variation is nearly independent of radius, and the 360 deg scatter decreased over time, from $\lesssim 0.17 \, \rm {dex}$ at $t_{\rm lb} = 11.6 \, \rm {Gyr}$ to $\sim 0.04 \, \rm {dex}$ at present-day. Consequently, stars at $t_{\rm lb} \gtrsim 8 \, \rm {Gyr}$ formed in a disc with primarily azimuthal scatter in abundances. All stars formed in amore »vertically homogeneous disc, Δ[Fe/H]$\le 0.02 \, \rm {dex}$ within $1 \, \rm {kpc}$ of the galactic mid-plane, with the exception of the young stars in the inner $\approx 4 \, \rm {kpc}$ at z ∼ 0. These results generally agree with our previous analysis of gas-phase elemental abundances, which reinforces the importance of cosmological disc evolution and azimuthal scatter in the context of stellar chemical tagging. We provide analytic fits to our results for use in chemical-tagging analyses.

« less

Abstract In the era of large-scale spectroscopic surveys in the Local Group, we can explore using chemical abundances of halo stars to study the star formation and chemical enrichment histories of the dwarf galaxy progenitors of the Milky Way (MW) and M31 stellar halos. In this paper, we investigate using the chemical abundance ratio distributions (CARDs) of seven stellar halos from the Latte suite of FIRE-2 simulations. We attempt to infer galaxies’ assembly histories by modeling the CARDs of the stellar halos of the Latte galaxies as a linear combination of template CARDs from disrupted dwarfs, with different stellar masses M ⋆ and quenching times t 100 . We present a method for constructing these templates using present-day dwarf galaxies. For four of the seven Latte halos studied in this work, we recover the mass spectrum of accreted dwarfs to a precision of <10%. For the fraction of mass accreted as a function of t 100 , we find the residuals of 20%–30% for five of the seven simulations. We discuss the failure modes of this method, which arise from the diversity of star formation and chemical enrichment histories that dwarf galaxies can take. These failure cases can be robustlymore »identified by the high model residuals. Although the CARDs modeling method does not successfully infer the assembly histories in these cases, the CARDs of these disrupted dwarfs contain signatures of their unusual formation histories. Our results are promising for using CARDs to learn more about the histories of the progenitors of the MW and M31 stellar halos.« less

The physics of magnetic fields (B) and cosmic rays (CRs) have recently been included in simulations of galaxy formation. However, significant uncertainties remain in how these components affect galaxy evolution. To understand their common observational tracers, we analyse the magnetic fields in a set of high-resolution, magnetohydrodynamic, cosmological simulations of Milky-Way-like galaxies from the FIRE-2 project. We compare mock observables of magnetic field tracers for simulations with and without CRs to observations of Zeeman splitting and rotation/dispersion measures. We find reasonable agreement between simulations and observations in both the neutral and the ionized interstellar medium (ISM). We find that the simulated galaxies with CRs show weaker ISM |B| fields on average compared to their magnetic-field-only counterparts. This is a manifestation of the effects of CRs in the diffuse, low density inner circumgalactic medium (CGM). We find that equipartition between magnetic and cosmic ray energy densities may be valid at large (> 1 kpc) scales for typical ISM densities of Milky-Way-like galaxies, but not in their haloes. Within the ISM, the magnetic fields in our simulated galaxies follow a power-law scaling with gas density. The scaling extends down to neutral hydrogen number densities < 300 cm−3, in contrast to observationally derived models,more »but consistent with the observational measurements. Finally, we generate synthetic rotation measure (RM) profiles for projections of the simulated galaxies and compare to observational constraints in the CGM. While consistent with upper limits, improved data are needed to detect the predicted CGM RMs at 10–200 kpc and better constrain theoretical predictions.

« less

We investigate stellar elemental abundance patterns at $z$ = 0 in eight low-mass ($M_{*}=10^{6}{-}10^{9}\ \text{M}_{\odot }$) galaxies in the Feedback in Realistic Environments cosmological simulations. Using magnesium (Mg) as a representative α-element, we explore stellar abundance patterns in magnesium-to-iron ([Mg/Fe]) versus iron-to-hydrogen ([Fe/H]), which follow an overall monotonic trend that evolved slowly over time. Additionally, we explore three notable secondary features in enrichment (in three different case-study galaxies) that arise from a galaxy merger or bursty star formation. First, we observe a secondary track with a lower [Mg/Fe] than the main trend. At $z$ = 0, stars from this track are predominantly found within 2–6 kpc of the centre; they were accreted in a 1:3 total-mass-ratio merger ∼0.4 Gyr ago. Second, we find a distinct elemental bimodality that forms following a strong burst in star formation in a galaxy at $t_{\text{lookback}}\, \sim 10$ Gyr. This burst quenched star formation for ∼0.66 Gyr, allowing Type Ia supernovae to enrich the system with iron (Fe) before star formation resumed. Third, we examine stripes in enrichment that run roughly orthogonal to the dominant [Mg/Fe] versus [Fe/H] trend; these stripes correspond to short bursts of star formation during which core-collapse supernovae enrich the surrounding medium with Mg (andmore »Fe) on short time-scales. If observed, these features would substantiate the utility of elemental abundances in revealing the assembly and star-formation histories of dwarf galaxies. We explore the observability of these features for upcoming spectroscopic studies. Our results show that precise measurements of elemental abundance patterns can reveal critical events in the formation histories of low-mass galaxies.

« less

ABSTRACT We develop a hybrid model of galactic chemical evolution that combines a multiring computation of chemical enrichment with a prescription for stellar migration and the vertical distribution of stellar populations informed by a cosmological hydrodynamic disc galaxy simulation. Our fiducial model adopts empirically motivated forms of the star formation law and star formation history, with a gradient in outflow mass loading tuned to reproduce the observed metallicity gradient. With this approach, the model reproduces many of the striking qualitative features of the Milky Way disc’s abundance structure: (i) the dependence of the [O/Fe]–[Fe/H] distribution on radius Rgal and mid-plane distance |z|; (ii) the changing shapes of the [O/H] and [Fe/H] distributions with Rgal and |z|; (iii) a broad distribution of [O/Fe] at sub-solar metallicity and changes in the [O/Fe] distribution with Rgal, |z|, and [Fe/H]; (iv) a tight correlation between [O/Fe] and stellar age for [O/Fe] > 0.1; (v) a population of young and intermediate-age α-enhanced stars caused by migration-induced variability in the Type Ia supernova rate; (vi) non-monotonic age–[O/H] and age–[Fe/H] relations, with large scatter and a median age of ∼4 Gyr near solar metallicity. Observationally motivated models with an enhanced star formation rate ∼2 Gyr ago improve agreement with the observed age–[Fe/H] and age–[O/H] relations, but worsen agreement with the observed age–[O/Fe]more »relation. None of our models predict an [O/Fe] distribution with the distinct bimodality seen in the observations, suggesting that more dramatic evolutionary pathways are required. All code and tables used for our models are publicly available through the Versatile Integrator for Chemical Evolution (VICE; https://pypi.org/project/vice).« less

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 abundancesmore »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.« less

ABSTRACT Kinematic studies of disc galaxies, using individual stars in the Milky Way or statistical studies of global disc kinematics over time, provide insight into how discs form and evolve. We use a high-resolution, cosmological zoom-simulation of a Milky Way-mass disc galaxy (h277) to tie together local disc kinematics and the evolution of the disc over time. The present-day stellar age–velocity relationship (AVR) of h277 is nearly identical to that of the analogous solar-neighbourhood measurement in the Milky Way. A crucial element of this success is the simulation’s dynamically cold multiphase ISM, which allows young stars to form with a low velocity dispersion (σbirth$\sim \!6 - 8 \ \mathrm{km\, s}^{-1}$) at late times. Older stars are born kinematically hotter (i.e. the disc settles over time in an ‘upside-down’ formation scenario), and are subsequently heated after birth. The disc also grows ‘inside-out’, and many of the older stars in the present-day solar neighbourhood are present because of radial mixing. We demonstrate that the evolution of σbirth in h277 can be explained by the same model used to describe the general decrease in velocity dispersion observed in disc galaxies from z ∼ 2–3 to the present-day, in which the disc evolves inmore »quasi-stable equilibrium and the ISM velocity dispersion decreases over time due to a decreasing gas fraction. Thus, our results tie together local observations of the Milky Way’s AVR with observed kinematics of high z disc galaxies.« less