Direct-Method Inferred Metallicity

• Direct-method inferred metallicity is defined by directly measuring electron temperatures using auroral-to-nebular line ratios to calculate gas-phase and stellar metallicities.
• It underpins key scaling relations such as the mass–metallicity relation and metallicity gradients by providing robust, gold-standard abundance calibrations.
• Advances in spectral stacking and IFU mapping have enhanced its accuracy, though challenges like temperature structure variations and spectral contamination remain.

Direct‐method inferred metallicity refers to the determination of astrophysical object metallicity—typically the gas-phase oxygen abundance in galactic H II regions or the photospheric metallicity in stars—through the direct measurement of fundamental physical properties such as electron temperature ($T_e$) or via direct spectral modeling, rather than through empirical or model-dependent strong-line calibrations. In galaxies and star-forming regions, this approach relies on temperature-sensitive emission line ratios (“auroral” to “nebular” lines) to robustly determine ionic abundances. Direct methods yield metallicity measures that serve as gold-standard calibrators for the mass–metallicity relation, abundance gradients, and inform models of chemical evolution, stellar feedback, and baryon cycling across cosmic time.

1. Fundamental Principles and Measurement Procedures

The direct method, as established for nebular emission, determines the gas-phase metallicity (often quantified as $12 + \log(\mathrm{O/H})$) by first measuring electron temperatures in ionized regions using weak auroral lines and then calculating ionic and total abundances based on atomic physics. For the canonical two-zone model of an H II region:

• Electron temperature measurement: For the high-ionization zone, the auroral [O III] $\lambda4363$ line intensity divided by that of [O III] $\lambda\lambda4959,5007$ is used to derive $T_e(\mathrm{O\,III})$. For the low-ionization zone, often [O II] $\lambda\lambda7320,7330$ or [S II] $\lambda\lambda4068,4076$ relative to their nebular counterparts, or empirical $T_e(\mathrm{O\,II})$–$T_e(\mathrm{O\,III})$ relations are adopted:

$$T_e(\mathrm{O\,II}) = 0.7 \times T_e(\mathrm{O\,III}) + 3000\ \mathrm{K}$$

• Ionic abundance determination: With $12 + \log(\mathrm{O/H})$0 and electron density (from, e.g., [S II] $12 + \log(\mathrm{O/H})$1 ratio), ionic abundances are calculated using analytic expressions or codes like PyNeb:

$12 + \log(\mathrm{O/H})$2

summing low- and high-ionization zone contributions:

$12 + \log(\mathrm{O/H})$3

• Application: This process, requiring detection of faint auroral lines, is often applied to stacked spectra of galaxies (Andrews et al., 2012, Khoram et al., 2024), resolved mapping with high S/N (Cameron et al., 2021, Hamel-Bravo et al., 2024), and in the regime of extremely metal-poor galaxies to define robust abundance scaling relations (Nakajima et al., 2022, Hirschauer et al., 2021).

For late-type stars (e.g., M dwarfs), direct metallicity is inferred from high-resolution IR spectroscopy and detailed spectral modeling using synthetic spectra matched to observed atomic (and sometimes molecular) line profiles, fitting for $12 + \log(\mathrm{O/H})$4 through $12 + \log(\mathrm{O/H})$5 minimization (Lindgren et al., 2015, Lardo et al., 2015).

2. Applications to Scaling Relations and Gradients

Direct-method metallicities have been central to the empirical calibration of fundamental galaxy relations, most notably:

• Mass–metallicity relation (MZR): Using stacked SDSS spectra of $12 + \log(\mathrm{O/H})$6200,000 galaxies, a direct-method MZR is established from $12 + \log(\mathrm{O/H})$7 (Andrews et al., 2012). The relation exhibits a steep low-mass slope ($12 + \log(\mathrm{O/H})$8), a turnover at $12 + \log(\mathrm{O/H})$9, and saturates at $\lambda4363$0 for high mass.
• Fundamental metallicity relation (FMR): Incorporating SFR, direct-method metallicities reveal a strong monotonic secondary dependence—at fixed $\lambda4363$1, higher SFR galaxies are more metal-poor; parameterized as

$\lambda4363$2

• Metallicity gradients: Stacked MaNGA spectra enable direct-method gradients across six radial bins (Khoram et al., 2024), revealing trends: flatter or positive at low mass, steeply negative at intermediate mass, flattening at high mass. Legacy results from strong-line methods are confirmed for global trends, but direct methods challenge the universality and SFR dependence attributed to the FMR.
• Dwarfs and low-luminosity systems: Direct-method measurements at low $\lambda4363$3 document a flattening of the MZR slope, and the presence of a metallicity “floor”—dwarfs rapidly self-enrich to a minimum abundance (Hirschauer et al., 2021).

Relation	Direct-method Slope/Trend	Turnover/Flattening	SFR/FMR Sensitivity
Mass–Metallicity (MZR)	$\lambda4363$4 at low mass	$\lambda4363$5	Strong, monotonic (α=0.66)
MZR (strong-line)	$\lambda4363$6 (shallower)	$\lambda4363$7	Weaker/complex/overlapping tracks

Direct-method scaling relations serve as benchmarks for both theoretical models (e.g., FIRE, IllustrisTNG (Gburek et al., 2022)) and extrapolations to high-redshift or reionization-era samples (Curti et al., 2022, Jones et al., 2020).

3. Comparison with Strong-line and Empirical Methods

There is systematic divergence between direct-method and strong-line (SEL) metallicity determinations:

• Steeper slope and lower turnover: The direct method yields a steeper low-mass slope and a lower mass turnover than seen in SEL MZRs (Andrews et al., 2012, Khoram et al., 2024).
• SFR dependence and FMR structure: Direct-method FMR shows a larger $\lambda4363$8 (0.66) and monotonic tracks with SFR, in contrast to overlapping or even reversed SFR trends at high mass in strong-line calibrations.
• Calibration of empirical indicators: Direct-method data underpin calibrations for strong-line indicators (e.g., $\lambda4363$9 O3N2). The $\lambda\lambda4959,5007$0 index emerges as the most robust, due to its sensitivity to both low- and high-ionization zone lines (Nakajima et al., 2022). Indicators that do not combine both stages (e.g., N2, O3N2) show increased scatter in the metal-poor regime ($\lambda\lambda4959,5007$1 dex for $\lambda\lambda4959,5007$2, larger for others).
• High-$\lambda\lambda4959,5007$3 and local analogues: At $\lambda\lambda4959,5007$4, strong-line calibrations must be carefully recalibrated; direct-method metallicities show that local calibrations often fail at matching high-$\lambda\lambda4959,5007$5 line ratios and metallicities (Curti et al., 2022, Sanders et al., 2019, Gburek et al., 2022).

4. Methodological Advances, Limitations, and Uncertainties

Advances

• Statistical power via stacking: Large stacking analyses enable detection of faint auroral lines across substantial bins of galaxy property space (mass, SFR, radius), greatly expanding the direct-method parameter regime (Andrews et al., 2012, Khoram et al., 2024).
• Non-parametric calibration frameworks: Genesis-metallicity leverages a 4D kernel density estimation across line ratios and EW(Hβ) for universal (redshift-invariant) metallicity estimation, and 5D calibration of low-ionization zone $\lambda\lambda4959,5007$6 estimators when auroral lines are not directly detected (Langeroodi et al., 2024).
• Resolved mapping: High S/N integral field spectrographs allow resolved $\lambda\lambda4959,5007$7 (and thus metallicity) mapping across galactic disks, outflows, and inflows on $\lambda\lambda4959,5007$8 pc scales (Cameron et al., 2021, Hamel-Bravo et al., 2024).

Limitations

• Temperature structure assumptions: Variations and underlying assumptions about $\lambda\lambda4959,5007$9 gradients within H II regions introduce substantial systematic uncertainties; classical $T_e(\mathrm{O\,III})$0–$T_e(\mathrm{O\,III})$1 conversions may bias gradients up to 0.5 dex, as shown via alternative auroral line diagnostics ([S II]) or empirical corrections (Cameron et al., 2020).
• Spectral contamination and ionization structure: Blending (e.g., [Fe II] $T_e(\mathrm{O\,III})$2 contaminating [O III] $T_e(\mathrm{O\,III})$3), non-thermal heating, and shocks may lead to anomalously high $T_e(\mathrm{O\,III})$4 at high metallicity, complicating abundance derivations (Khoram et al., 2024).
• Emission from diffuse ionized gas (DIG): Global spectra are subject to flux-weighting and DIG contamination, potentially biasing direct-method oxygen abundance by $T_e(\mathrm{O\,III})$5 dex (Sanders et al., 2017).
• Physical degeneracies in inversion models: For CGM and intergalactic absorbers, uncertainties in UV background shape and normalization propagate into factor-of-3 uncertainties (0.56 dex) in photoionization-corrected direct-method metallicity (Acharya et al., 2021).

5. Empirical and Physical Implications

Direct-method metallicity studies have yielded several key insights:

• Metal retention and galactic winds: Low-mass galaxies exhibit steeper MZR slopes and lower metallicities, reflecting the increased susceptibility of their shallow potential wells to metal-enriched outflows (Andrews et al., 2012, Hamel-Bravo et al., 2024). Observations of metal-rich outflows with $T_e(\mathrm{O\,III})$6 robustly support the scenario where galactic winds transport ISM metals to the CGM (Hamel-Bravo et al., 2024, Cameron et al., 2021).
• Baryon cycle and feedback: Simultaneous direct-method mapping of inflows (low-metallicity), outflows (enhanced metallicity), and the ISM reveals spatial segregation of baryon cycle processes (Cameron et al., 2021). The efficiency of metal ejection and mixing constrains chemo-dynamical models of galaxy evolution.
• Chemical evolution and enrichment sequence: The N/O versus O/H relation, with a plateau at low metallicity and a sharp upturn (transition to secondary nitrogen production), tightly correlates with stellar mass, signifying the link between integrated star formation history and enrichment (Andrews et al., 2012).
• Reionization and high-redshift assembly: Direct-method measurements at $T_e(\mathrm{O\,III})$7 using far-IR lines and JWST/NIRSpec-detected [O III] $T_e(\mathrm{O\,III})$8 record metallicities of $T_e(\mathrm{O\,III})$9 (Curti et al., 2022, Jones et al., 2020). The appearance of extreme deviations from the FMR in the most metal-poor, low-mass $\lambda\lambda7320,7330$0 galaxies implies departure from equilibrium between inflows, outflows, and enrichment during rapid assembly.
• Stellar metallicities from evolved stars: Direct spectroscopy of J-band atomic lines in RSGs and M dwarfs enables robust, model-dependent direct-method stellar metallicity measurements, cross-validated against binary companions (Lardo et al., 2015, Lindgren et al., 2015).

6. Calibration, Universal Estimators, and Prospects

The current consensus points toward direct-method calibrations as the authoritative basis for metallicity estimators:

• Universal calibrators: Genesis-metallicity’s non-parametric calibration of O2, O3, and EW(Hβ) in 4D achieves $\lambda\lambda7320,7330$1 dex accuracy, with strong redshift invariance from $\lambda\lambda7320,7330$2 to $\lambda\lambda7320,7330$3 (Langeroodi et al., 2024). The growing database of NIRSpec and ground-based direct-method measurements at $\lambda\lambda7320,7330$4 enables robust, universally applicable metallicity estimation.
• Physical parameter anchors: The integration of $\lambda\lambda7320,7330$5-based direct metallicities with measurements of UV continuum slope $\lambda\lambda7320,7330$6 and ionizing photon production efficiency $\lambda\lambda7320,7330$7 provides standardized benchmarks for interpreting high-redshift JWST spectra (Nakajima et al., 2022).
• Calibration of strong-line diagnostics: Direct-method samples (empirical and model-based) are necessary to recalibrate and interpret SEL indices, especially at high redshift where ISM conditions diverge from local templates (Nakajima et al., 2022, Curti et al., 2022).
• Spatial and temporal mapping: The advent of high spatial-resolution IFU spectroscopy enables mapping of metallicity gradients and small-scale mixing processes, crucial for tracing feedback-driven and accretion-driven enrichment. Extending such techniques to the early universe enables direct tracking of galaxy assembly and cosmic chemical evolution (Hamel-Bravo et al., 2024, Cameron et al., 2021, Khoram et al., 2024).
• Systematic uncertainties: Robust statistical cross-validation and awareness of limitations—temperature structure, line contamination, ionization corrections, UVB uncertainties—remain critical for accurate abundance determinations and their interpretation in the context of global galaxy evolution, the baryon cycle, and the cosmic metal budget.

7. Summary Table: Key Direct-Method Metallicty Paradigms

Domain	Diagnostic Lines/Indicators	Strengths	Key Uncertainties
H II regions (gas)	[O III] $\lambda\lambda7320,7330$8, [O II] $\lambda\lambda7320,7330$9	Physically motivated; robust to model bias	$\lambda\lambda4068,4076$0-structure, S/N, contamination
Stellar atmospheres	Atomic lines (J-band), spectral synthesis	Model-dependent yet precise for cool stars	Line lists, atmospheric models
Strong-line calibrators	$\lambda\lambda4068,4076$1, $\lambda\lambda4068,4076$2, etc.	Empirical, feasible at high-$\lambda\lambda4068,4076$3	Calibration dependency
Universal estimator	O2, O3, EW(Hβ); KDE in 4D	Universal $\lambda\lambda4068,4076$4 dex accuracy	Calibration sample, high-$\lambda\lambda4068,4076$5 ISM

Direct-method inferred metallicity thus constitutes the reference standard for galactic and extragalactic chemical abundance studies, underpinning the calibration of empirical diagnostics, enabling physically interpretable trends in galaxy evolution, and providing the pathway for robust metallicity mapping across cosmic time.