Metallicity offers a unique window into the baryonic history of the cosmos, being instrumental in probing evolutionary processes in galaxies between different cosmic environments. We aim to quantify the contribution of these environments to the scatter in the mass–metallicity relation (MZR) of galaxies. By analysing the galaxy distribution within the cosmic skeleton of the Horizon Run 5 cosmological hydrodynamical simulation at redshift z = 0.625, computed using a careful calibration of the T-ReX filament finder, we identify galaxies within three main environments: nodes, filaments, and voids. We also classify galaxies based on the dynamical state of the clusters and the length of the filaments in which they reside. We find that the cosmic environment significantly contributes to the scatter in the MZR; in particular, both the gas metallicity and its average relative standard deviation increase when considering denser large-scale environments. The difference in the average metallicity between galaxies within relaxed and unrelaxed clusters is ≈0.1dex, with both populations displaying positive residuals, δZg, from the averaged MZR. Moreover, the difference in metallicity between node and void galaxies accounts for $\approx 0.14 \, \text{dex}$ in the scatter of the MZR at stellar mass $M_{\star } \approx 10^{9.35}\, \text{M}_{\odot }$. Finally, both the average [O/Fe] in the gas and the galaxy gas fraction decrease when moving to higher large-scale densities in the simulation, suggesting that the cores of cosmic environments host – on average – older and more massive galaxies, whose enrichment is affected by a larger number of Type Ia Supernova events.
Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher.
Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?
Some links on this page may take you to non-federal websites. Their policies may differ from this site.
-
ABSTRACT -
ABSTRACT We derive empirical constraints on the nucleosynthetic yields of nitrogen by incorporating N enrichment into our previously developed and empirically tuned multizone galactic chemical evolution model. We adopt a metallicity-independent (‘primary’) N yield from massive stars and a metallicity-dependent (‘secondary’) N yield from AGB stars. In our model, galactic radial zones do not evolve along the observed [N/O]–[O/H] relation, but first increase in [O/H] at roughly constant [N/O], then move upward in [N/O] via secondary N production. By t ≈ 5 Gyr, the model approaches an equilibrium [N/O]–[O/H] relation, which traces the radial oxygen gradient. Reproducing the [N/O]–[O/H] trend observed in extragalactic systems constrains the ratio of IMF-averaged N yields to the IMF-averaged O yield of core-collapse supernovae. We find good agreement if we adopt $y_\text{N}^\text{CC}/y_\text{O}^\text{CC}=0.024$ and $y_\text{N}^\text{AGB}/y_\text{O}^\text{CC} = 0.062(Z/Z_\odot)$. For the theoretical AGB yields we consider, simple stellar populations release half their N after only ∼250 Myr. Our model reproduces the [N/O]–[O/H] relation found for Milky Way stars in the APOGEE survey, and it reproduces (though imperfectly) the trends of stellar [N/O] with age and [O/Fe]. The metallicity-dependent yield plays the dominant role in shaping the gas-phase [N/O]–[O/H] relation, but the AGB time-delay is required to match the stellar age and [O/Fe] trends. If we add ∼40 per cent oscillations to the star formation rate, the model reproduces the scatter in the gas phase [N/O]–[O/H] relation observed in external galaxies by MaNGA. We discuss implications of our results for theoretical models of N production by massive stars and AGB stars.
-
ABSTRACT Fluorine has many different potential sites and channels of production, making narrowing down a dominant site of fluorine production particularly challenging. In this work, we investigate which sources are the dominant contributors to the galactic fluorine by comparing chemical evolution models to observations of fluorine abundances in Milky Way stars covering a metallicity range of −2 < [Fe/H] < 0.4 and upper limits in the range of −3.4 < [Fe/H] < −2.3. In our models, we use a variety of stellar yield sets in order to explore the impact of varying both asymptotic giant branch (AGB) and massive star yields on the chemical evolution of fluorine. In particular, we investigate different prescriptions for initial rotational velocity in massive stars as well as a metallicity-dependent mix of rotational velocities. We find that the observed [F/O] and [F/Fe] abundance ratios at low metallicity and the increasing trend of [F/Ba] at [Fe/H] ≳ −1 can only be reproduced by chemical evolution models assuming, at all metallicities, a contribution from rapidly rotating massive stars with initial rotational velocities as high as 300 km s−1. A mix of rotational velocities may provide a more physical solution than the sole use of massive stars with vrot = 300 km s−1, which are predicted to overestimate the fluorine and average s-process elemental abundances at [Fe/H] ≳ −1. The contribution from AGB stars is predicted to start at [Fe/H] ≈ −1 and becomes increasingly important at high metallicity, being strictly coupled to the evolution of the nitrogen abundance. Finally, by using modern yield sets, we investigate the fluorine abundances of Wolf–Rayet winds, ruling them out as dominant contributors to the galactic fluorine.
-
Abstract We investigate the [X/Mg] abundances of 16 elements for 82,910 Galactic disk stars from GALAH+ DR3. We fit the median trends of low-Ia and high-Ia populations with a two-process model, which describes stellar abundances in terms of a prompt core-collapse and delayed Type-Ia supernova component. For each sample star, we fit the amplitudes of these two components and compute the residual Δ[X/H] abundances from this two-parameter fit. We find rms residuals ≲0.07 dex for well-measured elements and correlated residuals among some elements (such as Ba, Y, and Zn) that indicate common enrichment sources. From a detailed investigation of stars with large residuals, we infer that roughly 40% of the large deviations are physical and 60% are caused by problematic data such as unflagged binarity, poor wavelength solutions, and poor telluric subtraction. As one example of a population with distinctive abundance patterns, we identify 15 stars that have 0.3–0.6 dex enhancements of Na but normal abundances of other elements from O to Ni and positive average residuals of Cu, Zn, Y, and Ba. We measure the median elemental residuals of 14 open clusters, finding systematic ∼0.1–0.4 dex enhancements of O, Ca, K, Y, and Ba and ∼0.2 dex depletion of Cu in young clusters. Finally, we present a restricted three-process model where we add an asymptotic giant branch star (AGB) component to better fit Ba and Y. With the addition of the third process, we identify a population of stars, preferentially young, that have much higher AGB enrichment than expected from their SNIa enrichment.more » « less
-
ABSTRACT Using a sample of red giant stars from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) Data Release 16, we infer the conditional distribution $p([\alpha /{\rm Fe}]\, |\, [{\rm Fe}/{\rm H}])$ in the Milky Way disk for the α-elements Mg, O, Si, S, and Ca. In each bin of [Fe/H] and Galactocentric radius R, we model p([α/Fe]) as a sum of two Gaussians, representing ‘low-α’ and ‘high-α’ populations with scale heights $z_1=0.45\, {\rm kpc}$ and $z_2=0.95\, {\rm kpc}$, respectively. By accounting for age-dependent and z-dependent selection effects in APOGEE, we infer the [α/Fe] distributions that would be found for a fair sample of long-lived stars covering all z. Near the Solar circle, this distribution is bimodal at sub-solar [Fe/H], with the low-α and high-α peaks clearly separated by a minimum at intermediate [α/Fe]. In agreement with previous results, we find that the high-α population is more prominent at smaller R, lower [Fe/H], and larger |z|, and that the sequence separation is smaller for Si and Ca than for Mg, O, and S. We find significant intrinsic scatter in [α/Fe] at fixed [Fe/H] for both the low-α and high-α populations, typically ∼0.04-dex. The means, dispersions, and relative amplitudes of this two-Gaussian description, and the dependence of these parameters on R, [Fe/H], and α-element, provide a quantitative target for chemical evolution models and a test for hydrodynamic simulations of disk galaxy formation. We argue that explaining the observed bimodality will probably require one or more sharp transitions in the disk’s gas accretion, star formation, or outflow history in addition to radial mixing of stellar populations.more » « less
-
ABSTRACT We test the hypothesis that the observed first-peak (Sr, Y, Zr) and second-peak (Ba) s-process elemental abundances in low-metallicity Milky Way stars, and the abundances of the elements Mo and Ru, can be explained by a pervasive r-process contribution originating in neutrino-driven winds from highly magnetic and rapidly rotating proto-neutron stars (proto-NSs). We construct chemical evolution models that incorporate recent calculations of proto-NS yields in addition to contributions from asymptotic giant branch stars, Type Ia supernovae, and two alternative sets of yields for massive star winds and core-collapse supernovae. For non-rotating massive star yields from either set, models without proto-NS winds underpredict the observed s-process peak abundances by 0.3–$1\, \text{dex}$ at low metallicity, and they severely underpredict Mo and Ru at all metallicities. Models incorporating wind yields from proto-NSs with spin periods P ∼ 2–$5\, \text{ms}$ fit the observed trends for all these elements well. Alternatively, models omitting proto-NS winds but adopting yields of rapidly rotating massive stars, with vrot between 150 and $300\, \text{km}\, \text{s}^{-1}$, can explain the observed abundance levels reasonably well for [Fe/H] < −2. These models overpredict [Sr/Fe] and [Mo/Fe] at higher metallicities, but with a tuned dependence of vrot on stellar metallicity they might achieve an acceptable fit at all [Fe/H]. If many proto-NSs are born with strong magnetic fields and short spin periods, then their neutrino-driven winds provide a natural source for Sr, Y, Zr, Mo, Ru, and Ba in low-metallicity stellar populations. Conversely, spherical winds from unmagnetized proto-NSs overproduce the observed Sr, Y, and Zr abundances by a large factor.more » « 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] 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).more » « less
-
Abstract The SDSS-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey has obtained high-resolution spectra for thousands of red giant stars distributed among the massive satellite galaxies of the Milky Way (MW): the Large and Small Magellanic Clouds (LMC/SMC), the Sagittarius Dwarf Galaxy (Sgr), Fornax (Fnx), and the now fully disrupted Gaia Sausage/Enceladus (GSE) system. We present and analyze the APOGEE chemical abundance patterns of each galaxy to draw robust conclusions about their star formation histories, by quantifying the relative abundance trends of multiple elements (C, N, O, Mg, Al, Si, Ca, Fe, Ni, and Ce), as well as by fitting chemical evolution models to the [ α /Fe]–[Fe/H] abundance plane for each galaxy. Results show that the chemical signatures of the starburst in the Magellanic Clouds (MCs) observed by Nidever et al. in the α -element abundances extend to C+N, Al, and Ni, with the major burst in the SMC occurring some 3–4 Gyr before the burst in the LMC. We find that Sgr and Fnx also exhibit chemical abundance patterns suggestive of secondary star formation epochs, but these events were weaker and earlier (∼5–7 Gyr ago) than those observed in the MCs. There is no chemical evidence of a second starburst in GSE, but this galaxy shows the strongest initial star formation as compared to the other four galaxies. All dwarf galaxies had greater relative contributions of AGB stars to their enrichment than the MW. Comparing and contrasting these chemical patterns highlight the importance of galaxy environment on its chemical evolution.more » « less
-
Abstract We apply a novel statistical analysis to measurements of 16 elemental abundances in 34,410 Milky Way disk stars from the final data release (DR17) of APOGEE-2. Building on recent work, we fit median abundance ratio trends [X/Mg] versus [Mg/H] with a 2-process model, which decomposes abundance patterns into a “prompt” component tracing core-collapse supernovae and a “delayed” component tracing Type Ia supernovae. For each sample star, we fit the amplitudes of these two components, then compute the residuals Δ[X/H] from this two-parameter fit. The rms residuals range from ∼0.01–0.03 dex for the most precisely measured APOGEE abundances to ∼0.1 dex for Na, V, and Ce. The
correlations of residuals reveal a complex underlying structure, including a correlated element group comprised of Ca, Na, Al, K, Cr, and Ce and a separate group comprised of Ni, V, Mn, and Co. Selecting stars poorly fit by the 2-process model reveals a rich variety of physical outliers and sometimes subtle measurement errors. Residual abundances allow for the comparison of populations controlled for differences in metallicity and [α /Fe]. Relative to the main disk (R = 3–13 kpc), we find nearly identical abundance patterns in the outer disk (R = 15–17 kpc), 0.05–0.2 dex depressions of multiple elements in LMC and Gaia Sausage/Enceladus stars, and wild deviations (0.4–1 dex) of multiple elements inω Cen. The residual abundance analysis opens new opportunities for discovering chemically distinctive stars and stellar populations, for empirically constraining nucleosynthetic yields, and for testing chemical evolution models that include stochasticity in the production and redistribution of elements.