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.
ABSTRACT
The ultra-faint dwarf (UFD) galaxy Reticulum 2 (Ret 2) was recently discovered in images obtained by the Dark Energy Survey. We have observed the four brightest red giants in Ret 2 at high ...spectral resolution using the Michigan/
Magellan
Fiber System. We present detailed abundances for as many as 20 elements per star, including 12 elements heavier than the Fe group. We confirm previous detection of high levels of
r
-process material in Ret 2 (mean Eu/Fe = +1.69 ± 0.05) found in three of these stars (mean Fe/H = −2.88 ± 0.10). The abundances closely match the
r
-process pattern found in the well-studied metal-poor halo star CS 22892–052. Such
r
-process-enhanced stars have not been found in any other UFD galaxy, though their existence has been predicted by at least one model. The fourth star in Ret 2 (Fe/H = −3.42 ± 0.20) contains only trace amounts of Sr (Sr/Fe = −1.73 ± 0.43) and no detectable heavier elements. One
r
-process enhanced star is also enhanced in C (natal C/Fe ≈ +1.1). This is only the third such star known, which suggests that the nucleosynthesis sites leading to C and
r
-process enhancements are decoupled. The
r
-process-deficient star is enhanced in Mg (Mg/Fe = +0.81 ± 0.14), and the other three stars show normal levels of
α
-enhancement (mean Mg/Fe = +0.34 ± 0.03). The abundances of other
α
and Fe-group elements closely resemble those in UFD galaxies and metal-poor halo stars, suggesting that the nucleosynthesis that led to the large
r
-process enhancements either produced no light elements or produced light-element abundance signatures indistinguishable from normal supernovae.
The velocity anisotropy parameter, β, is a measure of the kinematic state of orbits in the stellar halo, which holds promise for constraining the merger history of the Milky Way (MW). We determine ...global trends for β as a function of radius from three suites of simulations, including accretion-only and cosmological hydrodynamic simulations. We find that the two types of simulations are consistent and predict strong radial anisotropy ( ) for Galactocentric radii greater than 10 kpc. Previous observations of β for the MW's stellar halo claim a detection of an isotropic or tangential "dip" at r ∼ 20 kpc. Using the N-body+SPH simulations, we investigate the temporal persistence, population origin, and severity of "dips" in β. We find that dips in the in situ stellar halo are long-lived, while dips in the accreted stellar halo are short-lived and tied to the recent accretion of satellite material. We also find that a major merger as early as z ∼ 1 can result in a present-day low (isotropic to tangential) value of β over a broad range of radii and angles. While all of these mechanisms are plausible drivers for the β dip observed in the MW, each mechanism in the simulations has a unique metallicity signature associated with it, implying that future spectroscopic surveys could distinguish between them. Since an accurate knowledge of β(r) is required for measuring the mass of the MW halo, we note that significant transient dips in β could cause an overestimate of the halo's mass when using spherical Jeans equation modeling.
ABSTRACT
We study the kinematics of stars both at their formation and today within 14 Milky Way (MW)-mass galaxies from the FIRE-2 cosmological zoom-in simulations. We quantify the relative ...importance of cosmological disc settling and post-formation dynamical heating. We identify three eras: a Pre-Disc Era (typically ≳ 8 Gyr ago), when stars formed on dispersion-dominated orbits; an Early-Disc Era (≈8–4 Gyr ago), when stars started to form on rotation-dominated orbits but with high velocity dispersion, σform; and a Late-Disc Era (≲ 4 Gyr ago), when stars formed with low σform. σform increased with time during the Pre-Disc Era, peaking ≈8 Gyr ago, then decreased throughout the Early-Disc Era as the disc settled and remained low throughout the Late-Disc Era. By contrast, the dispersion measured today, σnow, increases monotonically with age because of stronger post-formation heating for Pre-Disc stars. Importantly, most of σnow was in place at formation, not added post-formation, for stars younger than ≈10 Gyr. We compare the evolution of the three velocity components: at all times, σR, form > σϕ, form > σZ, form. Post-formation heating primarily increased σR at ages ≲ 4 Gyr but acted nearly isotropically for older stars. The kinematics of young stars in FIRE-2 broadly agree with the range observed across the MW, M31, M33, and PHANGS-MUSE galaxies. The lookback time that the disc began to settle correlates with its dynamical state today: earlier-settling galaxies currently form colder discs. Including stellar cosmic-ray feedback does not significantly change disc rotational support at fixed stellar mass.
ABSTRACT
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.
Determining the velocity distribution of halo stars is essential for estimating the mass of the Milky Way and for inferring its formation history. Since the stellar halo is a dynamically hot system, ...the velocity distribution of halo stars is well described by the three-dimensional velocity dispersions or by the velocity anisotropy parameter . Direct measurements of consistently suggest β = 0.5-0.7 for nearby halo stars. In contrast, the value of β at large Galactocentric radius r is still controversial, since reliable proper motion data are available for only a handful of stars. In the last decade, several authors have tried to estimate β for distant halo stars by fitting the observed line-of-sight velocities at each radius with simple velocity distribution models (local fitting methods). Some results of local fitting methods imply at , which is inconsistent with recent predictions from cosmological simulations. Here we perform mock-catalog analyses to show that the estimates of β based on local fitting methods are reliable only at with the current sample size (∼103 stars at a given radius). As r increases, the line-of-sight velocity (corrected for the solar reflex motion) becomes increasingly closer to the Galactocentric radial velocity, so it becomes increasingly more difficult to estimate the tangential velocity dispersion from the line-of-sight velocity distribution. Our results suggest that the forthcoming Gaia data will be crucial for understanding the velocity distribution of halo stars at .
ABSTRACT
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, 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.
ABSTRACT
Pressure balance plays a central role in models of the interstellar medium (ISM), but whether and how pressure balance is realized in a realistic multiphase ISM is not yet well understood. ...We address this question by using a set of FIRE-2 cosmological zoom-in simulations of Milky Way-mass disc galaxies, in which a multiphase ISM is self-consistently shaped by gravity, cooling, and stellar feedback. We analyse how gravity determines the vertical pressure profile as well as how the total ISM pressure is partitioned between different phases and components (thermal, dispersion/turbulence, and bulk flows). We show that, on average and consistent with previous more idealized simulations, the total ISM pressure balances the weight of the overlying gas. Deviations from vertical pressure balance increase with increasing galactocentric radius and with decreasing averaging scale. The different phases are in rough total pressure equilibrium with one another, but with large deviations from thermal pressure equilibrium owing to kinetic support in the cold and warm phases, which dominate the total pressure near the mid-plane. Bulk flows (e.g. inflows and fountains) are important at a few disc scale heights, while thermal pressure from hot gas dominates at larger heights. Overall, the total mid-plane pressure is well-predicted by the weight of the disc gas and we show that it also scales linearly with the star formation rate surface density (ΣSFR). These results support the notion that the Kennicutt–Schmidt relation arises because ΣSFR and the gas surface density (Σg) are connected via the ISM mid-plane pressure.
ABSTRACT
Giant molecular clouds (GMCs) are well studied in the local Universe, however, exactly how their properties vary during galaxy evolution is poorly understood due to challenging resolution ...requirements, both observational and computational. We present the first time-dependent analysis of GMCs in a Milky Way-like galaxy and an Large Magellanic Cloud (LMC)-like dwarf galaxy of the FIRE-2 (Feedback In Realistic Environments) simulation suite, which have sufficient resolution to predict the bulk properties of GMCs in cosmological galaxy formation self-consistently. We show explicitly that the majority of star formation outside the galactic centre occurs within self-gravitating gas structures that have properties consistent with observed bound GMCs. We find that the typical cloud bulk properties such as mass and surface density do not vary more than a factor of 2 in any systematic way after the first Gyr of cosmic evolution within a given galaxy from its progenitor. While the median properties are constant, the tails of the distributions can briefly undergo drastic changes, which can produce very massive and dense self-gravitating gas clouds. Once the galaxy forms, we identify only two systematic trends in bulk properties over cosmic time: a steady increase in metallicity produced by previous stellar populations and a weak decrease in bulk cloud temperatures. With the exception of metallicity, we find no significant differences in cloud properties between the Milky Way-like and dwarf galaxies. These results have important implications for cosmological star and star cluster formation and put especially strong constraints on theories relating the stellar initial mass function to cloud properties.
ABSTRACT
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 (and 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.