More Like this

1. The MOSDEF-LRIS Survey: The Interplay Between Massive Stars and Ionized Gas in High-Redshift Star-Forming Galaxies1

   https://doi.org/10.1093/mnras/staa1410  

   Topping, Michael W ; Shapley, Alice E ; Reddy, Naveen A ; Sanders, Ryan L ; Coil, Alison L ; Kriek, Mariska ; Mobasher, Bahram ; Siana, Brian ( July 2020 , Monthly Notices of the Royal Astronomical Society)

   Abstract We present a joint analysis of rest-UV and rest-optical spectra obtained using Keck/LRIS and Keck/MOSFIRE for a sample of 62 star-forming galaxies at z ∼ 2.3. We divide our sample into two bins based on their location in the [OIII]5007/Hβ vs. [NII]6584/Hα BPT diagram, and perform the first differential study of the rest-UV properties of massive ionizing stars as a function of rest-optical emission-line ratios. Fitting BPASS stellar population synthesis models, including nebular continuum emission, to our rest-UV composite spectra, we find that high-redshift galaxies offset towards higher [OIII]λ5007/Hβ and [NII]λ6584/Hα have younger ages ($\log (\textrm {~Age/yr})=7.20^{+0.57}_{-0.20}$) and lower stellar metallicities ($Z_*=0.0010^{+0.0011}_{-0.0003}$) resulting in a harder ionizing spectrum, compared to the galaxies in our sample that lie on the local BPT star-forming sequence ($\log (\textrm {Age/yr})=8.57^{+0.88}_{-0.84}$, $Z_*=0.0019^{+0.0006}_{-0.0006}$). Additionally, we find that the offset galaxies have an ionization parameter of $\log (U)=-3.04^{+0.06}_{-0.11}$ and nebular metallicity of ($12+\log (\textrm {~O/H})=8.40^{+0.06}_{-0.07}$), and the non-offset galaxies have an ionization parameter of $\log (U)=-3.11^{+0.08}_{-0.08}$ and nebular metallicity of $12+\log (\textrm {~O/H})=8.30^{+0.05}_{-0.06}$. The stellar and nebular metallicities derived for our sample imply that the galaxies offset from the local BPT relation are more α-enhanced ($7.28^{+2.52}_{-2.82}\textrm {~O/Fe}_{\odot }$) compared to those consistent with the local sequencemore »($3.04^{+0.95}_{-0.54}\textrm {~O/Fe}_{\odot }$). However, even galaxies that are entirely consistent with the local nebular excitation sequence appear to be α-enhanced – in contrast with typical local systems. Such differences must be considered when estimating gas-phase oxygen abundances at high redshift based on strong emission-line ratios. Specifically, a similarity in the location of high-redshift and local galaxies in the BPT diagram may not be indicative of a similarity in their physical properties.« less
2. The MOSDEF survey: a comprehensive analysis of the rest-optical emission-line properties of z ∼ 2.3 star-forming galaxies

   https://doi.org/10.1093/mnras/stab119  

   Runco, Jordan N ; Shapley, Alice E ; Sanders, Ryan L ; Topping, Michael W ; Kriek, Mariska ; Reddy, Naveen A ; Coil, Alison L ; Mobasher, Bahram ; Siana, Brian ; Freeman, William R ; et al ( February 2021 , Monthly Notices of the Royal Astronomical Society)

   ABSTRACT We analyse the rest-optical emission-line spectra of z ∼ 2.3 star-forming galaxies in the complete MOSFIRE Deep Evolution Field (MOSDEF) survey. In investigating the origin of the well-known offset between the sequences of high-redshift and local galaxies in the [O iii]λ5008/Hβ versus [N ii]λ6585/Hα (‘[N ii] BPT’) diagram, we define two populations of z ∼ 2.3 MOSDEF galaxies. These include the high population that is offset towards higher [O iii]λ5008/Hβ and/or [N ii]λ6585/Hα with respect to the local SDSS sequence and the low population that overlaps the SDSS sequence. These two groups are also segregated within the [O  iii]λ5008/Hβ versus [S ii]λλ6718,6733/Hα and the [O iii]λλ4960,5008/[O ii ]λλ3727,3730 (O32) versus ([O  iii]λλ4960,5008+[O ii]λλ3727,3730)/Hβ (R23) diagrams, which suggests qualitatively that star-forming regions in the more offset galaxies are characterized by harder ionizing spectra at fixed nebular oxygen abundance. We also investigate many galaxy properties of the split sample and find that the high sample is on average smaller in size and less massive, but has higher specific star formation rate (SFR) and SFR surface density values and is slightly younger compared to the low population. From Cloudy+BPASS photoionization models, we estimate that the high population has a lower stellar metallicity (i.e. harder ionizing spectrum) but slightly higher nebular metallicity and higher ionizationmore »parameter compared to the low population. While the high population is more α-enhanced (i.e. higher α/Fe) than the low population, both samples are significantly more α-enhanced compared to local star-forming galaxies with similar rest-optical line ratios. These differences must be accounted for in all high-redshift star-forming galaxies – not only those ‘offset’ from local excitation sequences.« less
3. The MOSDEF survey: direct-method metallicities and ISM conditions at z ∼ 1.5–3.5

   https://doi.org/10.1093/mnras/stz3032  

   Sanders, Ryan L. ; Shapley, Alice E. ; Reddy, Naveen A. ; Kriek, Mariska ; Siana, Brian ; Coil, Alison L. ; Mobasher, Bahram ; Shivaei, Irene ; Freeman, William R. ; Azadi, Mojegan ; et al ( November 2019 , Monthly Notices of the Royal Astronomical Society)

   We present detections of [O iii] λ4363 and direct-method metallicities for star-forming galaxies at z = 1.7–3.6. We combine new measurements from the MOSFIRE Deep Evolution Field (MOSDEF) survey with literature sources to construct a sample of 18 galaxies with direct-method metallicities at z > 1, spanning 7.5 < 12+log(O/H) < 8.2 and log(M*/M⊙) = 7–10. We find that strong-line calibrations based on local analogues of high-redshift galaxies reliably reproduce the metallicity of the z > 1 sample on average. We construct the first mass–metallicity relation at z > 1 based purely on direct-method O/H, finding a slope that is consistent with strong-line results. Direct-method O/H evolves by ≲0.1 dex at fixed M* and star formation rate from z ∼ 0 to 2.2. We employ photoionization models to constrain the ionization parameter and ionizing spectrum in the high-redshift sample. Stellar models with supersolar O/Fe and binary evolution of massive stars are required to reproduce the observed strong-line ratios. We find that the z > 1 sample falls on the z ∼ 0 relation between ionization parameter and O/H, suggesting no evolution of this relation from z ∼ 0 to z ∼ 2. These results suggest that the offset of the strong-line ratios of this sample from local excitation sequences is driven primarilymore »by a harder ionizing spectrum at fixed nebular metallicity compared to what is typical at z ∼ 0, naturally explained by supersolar O/Fe at high redshift caused by rapid formation time-scales. Given the extreme nature of our z > 1 sample, the implications for representative z ∼ 2 galaxy samples at ∼1010 M⊙ are unclear, but similarities to z > 6 galaxies suggest that these conclusions can be extended to galaxies in the epoch of reionization.

   « less
4. The BPT Diagram in Cosmological Galaxy Formation Simulations: Understanding the Physics Driving Offsets at High Redshift

   https://doi.org/10.3847/1538-4357/ac43b8  

   Garg, Prerak ; Narayanan, Desika ; Byler, Nell ; Sanders, Ryan L. ; Shapley, Alice E. ; Strom, Allison L. ; Davé, Romeel ; Hirschmann, Michaela ; Lovell, Christopher C. ; Otter, Justin ; et al ( February 2022 , The Astrophysical Journal)

   The Baldwin, Philips, & Terlevich diagram of [Oiii]/Hβversus [Nii]/Hα(hereafter N2-BPT) has long been used as a tool for classifying galaxies based on the dominant source of ionizing radiation. Recent observations have demonstrated that galaxies atz∼ 2 reside offset from local galaxies in the N2-BPT space. In this paper, we conduct a series of controlled numerical experiments to understand the potential physical processes driving this offset. We model nebular line emission in a large sample of galaxies, taken from thesimbacosmological hydrodynamic galaxy formation simulation, using thecloudyphotoionization code to compute the nebular line luminosities from Hiiregions. We find that the observed shift toward higher [Oiii]/Hβand [Nii]/Hαvalues at high redshift arises from sample selection: when we consider only the most massive galaxiesM_*∼ 10^10–11M_⊙, the offset naturally appears, due to their high metallicities. We predict that deeper observations that probe lower-mass galaxies will reveal galaxies that lie on a locus comparable toz∼ 0 observations. Even when accounting for samples-selection effects, we find that there is a subtle mismatch between simulations and observations. To resolve this discrepancy, we investigate the impact of varying ionization parameters, Hiiregion densities, gas-phase abundance patterns, and increasing radiation field hardness on N2-BPT diagrams. We find that either decreasing themore »ionization parameter or increasing the N/O ratio of galaxies at fixed O/H can move galaxies along a self-similar arc in N2-BPT space that is occupied by high-redshift galaxies.

   « less
5. Infrared Spectral Energy Distributions and Dust Masses of Sub-solar Metallicity Galaxies at z ∼ 2.3

   https://doi.org/10.3847/1538-4357/ac54a9  

   Shivaei, Irene ; Popping, Gergö ; Rieke, George ; Reddy, Naveen ; Pope, Alexandra ; Kennicutt, Robert ; Mobasher, Bahram ; Coil, Alison ; Fudamoto, Yoshinobu ; Kriek, Mariska ; et al ( March 2022 , The Astrophysical Journal)

   Abstract We present results from Atacama Large Millimeter/submillimeter Array (ALMA) 1.2 mm continuum observations of a sample of 27 star-forming galaxies at z = 2.1–2.5 from the MOSFIRE Deep Evolution Field survey with metallicity and star formation rate measurements from optical emission lines. Using stacks of Spitzer, Herschel, and ALMA photometry (rest frame ∼8–400 μ m), we examine the infrared (IR) spectral energy distributions (SED) of z ∼ 2.3 subsolar-metallicity (∼0.5 Z ⊙ ) luminous infrared galaxies (LIRGs). We find that the data agree well with an average template of higher-luminosity local low-metallicity dwarf galaxies (reduced χ 2 = 1.8). When compared with the commonly used templates for solar-metallicity local galaxies or high-redshift LIRGs and ultraluminous IR galaxies, even in the most favorable case (with reduced χ 2 = 2.8), the templates are rejected at >98% confidence. The broader and hotter IR SED of both the local dwarfs and high-redshift subsolar-metallicity galaxies may result from different grain properties or a harder/more intense ionizing radiation field that increases the dust temperature. The obscured star formation rate (SFR) indicated by the far-IR emission of the subsolar-metallicity galaxies is only ∼60% of the total SFR, considerably lower than that of the local LIRGsmore »with ∼96%–97% obscured fractions. Due to the evolving IR SED shape, the local LIRG templates fit to mid-IR data overestimate the Rayleigh–Jeans tail measurements by a factor of 2–20. These templates underestimate IR luminosities if fit to the observed ALMA fluxes by >0.4 dex. At a given stellar mass or metallicity, dust masses at z ∼ 2.3 are an order of magnitude higher than z ∼ 0. Given the predicted molecular gas fractions, the observed z ∼ 2.3 dust-to-stellar mass ratios suggest lower dust-to-molecular gas masses than in local galaxies with similar metallicities.« less