The LBT $Y_{\rm p}$ Project II: MODS Spectra, Physical Conditions, and Oxygen Abundances in Local Metal-Poor Nebulae

TL;DR

This work addresses the determination of the primordial helium mass fraction $Y_{ m p}$ by assembling a homogeneous, high‑fidelity dataset of low‑metallicity nebulae with direct temperature and density diagnostics. Using deep rest‑frame optical spectroscopy from MODS on the Large Binocular Telescope, the authors measure $T_e$ and $n_e$ across multiple ionization zones for 62 galaxies, derive direct O/H abundances with median uncertainty ~4%, and quantify the relationships between ionic abundances and physical conditions. A core result is the empirical Te[S III]–Te[O III] scaling, with a tight $T_e$–$T_e$ relation that improves abundance inferences in metal-poor environments, and the discovery of systematic density stratification with higher densities in high‑ionization zones. The dataset also reveals widespread narrow He II emission and very high ionization, reinforcing the need for multizone diagnostics; overall, the study provides a robust, uniform basis for estimating $Y_{ m p}$ and for broader ISM abundance analyses relevant to JWST‑era observations, with implications for chemical evolution and primordial nucleosynthesis constraints.

Abstract

Empirically measuring the primordial He mass fraction, $Y_{\rm p}$, requires a significant number of low-metallicity nebulae with direct constraints on He/H and O/H abundances. This technique requires high-fidelity measurements of the gas-phase physical conditions, namely the electron temperature ($T_e$) and density ($n_e$). To this end, we present deep rest-optical spectroscopy for a sample of 62 low-metallicity ($\lesssim$ 20% solar O/H) galaxies acquired using the Multi-Object Double Spectrographs (MODS) on the Large Binocular Telescope (LBT) as part of the LBT $Y_{\rm p}$ Project. We discuss new fitting methods that recover the intensity of up to 61 H and He recombination lines, of which, up to 26 will be used to determine gas-phase He abundances, and we examine the emission line properties of the LBT $Y_{\rm p}$ Project sample. We assess different scaling relations in the low-metallicity interstellar medium (ISM), finding that $n_e$[Ar IV] measured in 31 targets is systematically larger than $n_e$[S II] or $n_e$[O II]. The larger densities are insufficient to significantly bias $T_e$[O III] or the O/H abundance. $T_e$[S III] and $T_e$[O III] are strongly correlated over a range of $\sim$10$^4$ K with very low scatter, and we calibrate new $T_e$[S III]-$T_e$[O III] scaling relations for use in other low-metallicity environments. We examine different $T_e$ measured in the low-ionization gas, finding significant scatter compared to $T_e$[O III]. The precision direct O/H derived in this analysis (median uncertainty $\sim$4%) are consistent with prior literature measurements, albeit with relatively large scatter. These data provide a key component necessary to empirically measure $Y_{\rm p}$ and the abundance patterns of other elements in the ISM.

Figures (9)

• Figure 1: Combined MODS spectrum of the low-metallicity galaxy HS0811+4913. Faint $T_e$- and $n_e$-sensitive emission lines are highlighted in the subplots. The $T_e$-sensitive auroral line of [O3]$\lambda$4363 is clearly visible without magnification. The spectrum shows a distinct Balmer and Paschen jump, along with numerous intense and faint He1 lines crucial for the determination of He$^+$/H$^+$.
• Figure 2: Examples of emission line fits in the galaxy HS0811+4913. The left two panels uses a combination of a linear continuum plus Gaussian emission line profiles to fit the H14, H13, and [O2]$\lambda\lambda$3726,29 (first panel) and [O3]$\lambda$5007 and He1$\lambda$5016 (second panel). The 1D spectrum is plotted in black, the full fit is provided in orange, individual narrow emission lines are in alternating blue and red. The standard Gaussian fits have normalized residuals (lower panels) that can exceed 20%. The right two panels plot new fits using super-Gaussian profiles for the emission lines. This provides significantly better fits to the blend of H1 and [O2] lines in the third panel. The inclusion of a two-component broad wing on [O3]$\lambda$5007 (fourth panel, shown in pink and purple) provides a reliable fit to both the strong [O3] line and the faint He1 line.
• Figure 3: Emission line diagrams proposed by bald1981 and veil1987: log(I([O3]$\lambda$5007)/I(H$\beta$)) vs. log(I([N2]$\lambda$6584)/I(H$\alpha$)) (Left), log(I([S2]$\lambda$6717,31)/I(H$\alpha$)) (Center), and log(I([O1]$\lambda$6300)/I(H$\alpha$)) (Right). Line intensity ratios measured in the LBT $Y_{\rm p}$ Project sample are plotted as blue diamonds. Line ratios measured in relatively metal-rich H2 regions from the CHAOS project are plotted as purple circles for reference. The kewl2001 and kauf2003 demarcations proposed to distinguish between stellar and AGN ionizing sources are plotted as black dashed and green dotted lines, respectively. The LBT $Y_{\rm p}$ Project targets are characterized by highly-ionized, low-metallicity ISM and have line ratios consistent with extreme stellar ionization.
• Figure 4: Additional line ratio diagrams to highlight the highly-ionized ISM of the LBT $Y_{\rm p}$ Project targets. The shape and color of each point is consistent with the labeling from Figure \ref{['fig:bptfigs']}. Left: log(O$_{32}$) vs. log(R$_{23}$), the ionization-excitation diagram. Center: log(I([Ne3]$\lambda$3869)/I([O2]$\lambda$3727)) vs. log(I([Ne3]) + I([O2])/I(H$\delta$)), an alternative ionization-excitation diagram. Right: log(I(He2$\lambda$4686)/I(H$\beta$)) vs. log(I([N2]$\lambda$6584)/I(H$\alpha$)), a diagram to highlight nebular He2 emission. The red dot-dashed line indicates line ratios from 10% AGN ionization shir2012. All diagrams highlight the high ionization and excitation conditions observed in the targeted sample, while the line ratios remain consistent with predominantly stellar ionization.
• Figure 5: Direct density trends measured in the LBT $Y_{\rm p}$ Project sample. $n_e$[S2] is plotted against $n_e$[O2] in the low-ionization zone (Left), $n_e$[Cl3] in the intermediate-ionization zone (Center), and $n_e$[Ar4] in the high-ionization zone (Right). The dotted black line represents equivalent densities. The number of nebulae with simultaneous $n_e$ is provided in the bottom right of each panel. The [S2] and [O2] densities scatter around the one-to-one line, but $n_e$[Cl3] and $n_e$[Ar4] are systematically higher than the low-ionization zone densities.