Pith. sign in

REVIEW 3 major objections 7 minor 98 references

Ionized nebulae have no single density: each forbidden-line diagnostic reports only the part of a broad density distribution it is tuned to sense.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · grok-4.5

2026-07-10 14:28 UTC pith:N4LLDHM6

load-bearing objection The multi-diagnostic hierarchy is real and cleanly demonstrated; the ~10% dense-gas fraction is the only soft quantitative claim. the 3 major comments →

arxiv 2607.07973 v1 pith:N4LLDHM6 submitted 2026-07-08 astro-ph.GA

There is no single density: star-forming regions and galaxies hold more dense ionized gas than long assumed

classification astro-ph.GA
keywords H II regionselectron densityforbidden-line diagnosticsdensity distributionionized gasstar-forming galaxiesnebular spectroscopy
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

For decades, different density-sensitive emission lines drawn from the same ionized gas have returned electron densities that disagree by up to two orders of magnitude. The usual explanations—ionization stratification or imperfect atomic data—are incomplete. This paper shows that each diagnostic is intrinsically most sensitive at a particular density n_M set by atomic physics. When real nebulae contain gas spanning a wide density range, each line simply reports the part of that range it is tuned to. Across the Orion Nebula, local H II regions, nearby star-forming galaxies and high-redshift systems, the densities recovered by many independent diagnostics therefore follow a tight linear relation with their n_M values. The slope of that relation implies that a substantial fraction of the emission measure comes from dense gas (log n_e ≳ 4) that low-n_M diagnostics miss by construction. There is therefore no single electron density to measure; every mass, pressure, abundance and feedback estimate that rests on the single-density assumption inherits a predictable bias and must be revisited.

Core claim

The long-standing hierarchy among forbidden-line density diagnostics is neither ionization stratification nor atomic-data error. It is the natural result of folding a broad, unresolved electron-density probability distribution through the distinct atomic response kernels of each ratio. Median densities therefore obey the empirical relation log n_e,obs ≈ 0.29 log n_M + 1.52 from the Orion Nebula through local and high-redshift galaxies, revealing that ionized gas contains far more dense material than any single diagnostic implies.

What carries the argument

The maximum-sensitivity density n_M — the electron density at which the logarithmic derivative of a line ratio is largest — together with the linear relation that maps each diagnostic’s n_M onto the density it actually recovers when the gas spans a broad emission-measure distribution.

Load-bearing premise

That a single power-law emission-measure distribution with a slope near −1.3 adequately describes the true multi-phase density field of ionized nebulae.

What would settle it

A large sample of objects in which many independent density diagnostics spanning several decades in n_M return essentially the same density (slope near zero), or high-resolution spectroscopy of additional kinematically isolated dense knots that still show a hierarchy once the ambient nebula is cleanly removed.

Watch this falsifier — get emailed when new claim-graph text bears on it.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 7 minor

Summary. The paper argues that the long-standing disagreement among forbidden-line electron-density diagnostics is neither primarily ionization stratification nor atomic-data error, but the natural consequence of folding a broad unresolved density PDF through the distinct atomic response kernels of each ratio. Defining a maximum-sensitivity density n_M for each diagnostic, the authors show that median densities measured spaxel-by-spaxel across 1226 LVM spaxels in Orion obey log n_e,obs = (0.29 ± 0.05) log n_M + (1.52 ± 0.27). The same ordered hierarchy appears in DESIRED H II regions and local star-forming galaxies, and qualitatively in a small high-z sample. A minimal power-law emission-measure forward model (H_n ∝ n_e^β, β ≈ −1.3) reproduces the Orion slope without ionization stratification; a Bayesian fit implies ~10% of the emission measure above log n_e ≳ 4. The hierarchy vanishes in kinematically isolated HH 202-S/204 components, providing a narrow-PDF control. The authors conclude that nebulae have no single density and that masses, pressures, abundances and energetics based on the single-density assumption must be reconsidered.

Significance. If the central observational result holds, this is a high-impact clarification of a foundational tool in nebular spectroscopy. The multi-diagnostic Orion maps (16 ratios, same 1226 spaxels), the persistence of the hierarchy within single ionic species (S+, O+, Fe2+), the HH 202-S/204 narrow-PDF control, and the DESIRED extension are strong, falsifiable pieces of evidence that go well beyond re-stating known diagnostic scatter. The n_M framework and the simple linear relation provide a concrete, reusable diagnostic of unresolved density structure. The quantitative dense-gas fraction and the full suite of implied biases on masses/pressures/abundances are more model-dependent, but the qualitative claim that single-density analyses are systematically biased is well supported and would affect analyses from local H II regions to high-z galaxies.

major comments (3)
  1. Methods, Bayesian constraints section and Extended Data Fig. 8: The quoted (10 +5/-4)% emission-measure fraction above log n_e ≳ 4 (and the precise β ≈ −1.3) is obtained only after adopting a pure power-law H_n and, in the preferred run, an informative n_rms prior. Other broad PDFs (log-normal + power-law tail, multi-component) can reproduce a slope ~0.3 with different dense fractions. The paper already notes that every viable model needs a dense component, but the abstract and main text present the ~10% figure as a firm implication. Please either (i) test at least one alternative functional form and report the range of dense fractions, or (ii) explicitly label the 10% result as power-law-model-dependent and move the quantitative claim out of the abstract-level summary.
  2. Fig. 2 and associated text (high-z SFGs): The high-redshift sample comprises only three systems (plus [O III] λ4363/λ5007-only measurements from Arellano-Córdova et al.). The statement that the relation “appears to extend” and that medians fall “above the Orion relation as expected for broader or denser distributions” is therefore under-constrained. Temper the high-z claim to “consistent with, but not yet a statistical demonstration of, the local hierarchy,” and avoid implying redshift evolution of the PDF width from the present sample.
  3. Implications paragraphs (thermal structure, masses, pressures, abundances, feedback): The qualitative direction of the bias is clear from the n_e,obs–n_M slope, but the manuscript asserts that masses, filling factors, abundances and outflow energetics “must be reconsidered” without even order-of-magnitude estimates of the bias for a typical [S II]-based analysis. A short worked example (e.g., how a power-law PDF with the Orion parameters shifts M_ion or P/k relative to a single n_e([S II]) assumption) would make the load-bearing claim quantitative rather than programmatic.
minor comments (7)
  1. Eq. (2) fit: State explicitly in the main text (not only Methods) that the fit uses the nine algebraically independent ratios and that the slope remains ~0.29–0.33 under alternative independent sets, so readers do not over-interpret the p-value from all 16 points.
  2. Methods, maximum-sensitivity density: For multi-peak Fe III and Fe II sensitivity curves, the peak-selection rule (highest peak above min critical density) is reasonable but should be illustrated for one diagnostic in Extended Data so the choice is reproducible.
  3. Analytic solution / Bayesian fit: The empirical contrast-factor scaling δ = 2 (n_M/10^3 cm^{-3})^{1/2} is convenient but ad hoc. Note the sensitivity of posterior medians to this choice, or replace it with diagnostic-specific δ from the actual critical densities where available.
  4. Fig. 1 caption: Clarify that grey bars are spatial 16th–84th percentiles across spaxels, not measurement errors; this is stated but easy to misread when comparing to Fig. 2 error bars (object-to-object percentiles).
  5. Extended Data Table 1: Mark which ratios enter the Orion linear fit versus which are displayed only, and note which are available only as blended sums in DESIRED, to avoid confusion when comparing Fig. 1 and Fig. 2.
  6. Planetary nebulae are mentioned as showing a shallower slope (Méndez-Delgado et al., in prep.). A one-sentence quantitative preview (slope/intercept range) would help readers judge whether the Orion calibration is H II-specific or more general.
  7. Typographical consistency: “Hiiregions” / “Hii regions” / “H II regions” appear in mixed forms; standardize to H II throughout.

Circularity Check

1 steps flagged

Empirical n_e,obs–n_M hierarchy is independent of models; power-law parameters are fitted to it (and the ~10% dense fraction is a posterior), but this is not load-bearing for the central claim or the HH narrow-PDF control.

specific steps
  1. fitted input called prediction [Forward model / Fig. 4 and Bayesian constraints section (Methods)]
    "Fig. 4 shows that this minimal model reproduces the observed Orion sequence with remarkable fidelity (m=0.29, b=1.50, r=0.84, compared to the observed m=0.29, b=1.52), with no tuning of the slope to the data. ... the observed Orion relation is reproduced by distributions with slope β≈−1.3 spanning nearly five orders of magnitude in density ... Integrating the posterior emission-measure distribution, one finds that (10+5−4)% of the recombination-line emission originates from gas with log10 ne ≳4"

    β (and n_min) are free parameters of the assumed power-law H_n; they are chosen/fitted so that the synthetic diagnostics recover the already-measured slope m≈0.29 of the Orion sequence. The model therefore matches its target by construction for that functional form, and the quoted ~10% dense-gas fraction is a direct posterior of the same fit rather than an independent prediction. The paper’s phrasing “with no tuning” overstates independence, but the step is not load-bearing for the empirical hierarchy or the HH control.

full rationale

The load-bearing result is the ordered sequence of median densities versus atomic n_M (Eq. 2, Fig. 1) measured spaxel-by-spaxel on new LVM data, its persistence across DESIRED objects and high-z systems without refitting (Fig. 2), and its disappearance in the kinematically isolated HH 202-S/204 components (Extended Data Fig. 9). n_M itself is defined solely from atomic response functions (Eq. 1) and is independent of the nebular data. The power-law emission-measure model (Eq. 3) and its Bayesian fit (Methods; Extended Data Fig. 8) are used only to interpret the already-observed slope and to quote a dense-gas emission-measure fraction; they do not generate the hierarchy. The narrow-PDF validation (Extended Data Fig. 7) and the HH control supply independent falsification tests that do not rely on the fitted β or n_min. No equation reduces the central claim to a definitional identity or to a self-citation uniqueness theorem. The single mild circularity is the presentation of a tuned model as reproducing the data “with no tuning.” This is proportionate to score 2, not higher.

Axiom & Free-Parameter Ledger

4 free parameters · 5 axioms · 0 invented entities

The observational hierarchy stands on standard atomic physics and public spectroscopy. The quantitative claim that a dense tail contributes ~10 % of the emission measure rests on three free parameters of an assumed power-law PDF, an isothermal temperature field, and the absence of ionization–density correlations in the forward model. No new physical entities are postulated; n_M is a derived atomic quantity.

free parameters (4)
  • power-law slope β = ≈ −1.3 (forward model); −1.41^{+0.13}_{-0.15} (informative prior)
    Fitted to the Orion sequence (posterior β ≈ −1.3 to −1.4); controls the slope of the diagnostic hierarchy and the mass fraction in the dense tail.
  • lower density bound log10 n_min = 1.47^{+0.44}_{-1.22} (uninformative); 1.74^{+0.27}_{-0.32} (informative)
    Fitted; sets the intercept of the observed relation and the characteristic bulk density.
  • upper density bound log10 n_max (or width) = >5.2 (68 % credibility, informative prior)
    Only a lower limit is constrained by the data; the high-density tail is effectively unbounded and contributes the dense-gas fraction.
  • empirical contrast-factor scaling δ = 2 (n_M / 10^3 cm^{-3})^{1/2} = scaling coefficient 2
    Adopted for the analytic two-level solution to keep the Bayesian forward model fast; calibrated to the range of real diagnostics rather than computed ion-by-ion.
axioms (5)
  • ad hoc to paper A power-law emission-measure distribution H_n ∝ n_e^β is an adequate minimal description of the unresolved density field.
    Chosen as the simplest continuous form with few free parameters; other functional forms (log-normal + power-law tail, multi-phase) are not explored for the quantitative fractions.
  • domain assumption Electron temperature is essentially isothermal (t^{2} = 5×10^{-5}) so that the diagnostic hierarchy arises solely from density structure.
    Explicitly imposed in the forward model (Methods) to isolate the density effect; real nebulae have temperature variations that could couple to density.
  • domain assumption No ionization stratification or spatial density–ionization correlations are required to produce the hierarchy.
    Built into the forward model; supported by the persistence of the hierarchy within single ionic species (S+, O+, Fe2+), but not proven to be zero in real nebulae.
  • domain assumption Standard collisional-radiative atomic data (PyNeb data set of Méndez-Delgado et al. 2023) correctly describe the emissivities of all 16 diagnostics.
    Assumed throughout; the continuous monotonic sequence across independent ions and calculations is used to argue that atomic-data errors cannot produce the observed hierarchy.
  • standard math The maximum-sensitivity density n_M is well-defined as arg max |d log R / d log n_e| and, for multi-peak Fe ions, is the highest-sensitivity peak above the minimum critical density.
    Definition given in equation (1) and Methods; the selection rule for multi-peak cases is a reasonable but non-unique choice.

pith-pipeline@v1.1.0-grok45 · 31660 in / 3520 out tokens · 37929 ms · 2026-07-10T14:28:11.601015+00:00 · methodology

0 comments
read the original abstract

Ionized gas fills star-forming regions and galaxies, and nearly everything we know about its temperature, pressure, mass, and composition is inferred from its emission lines [1-3]. The electron density is needed for all of these, yet a longstanding puzzle has resisted explanation: different density-sensitive lines, applied to the same gas, return values that disagree by up to two orders of magnitude. This is usually attributed either to each line tracing a physically distinct ionization zone or to imperfect atomic data [4-7]. Here we show that the disagreement is neither a flaw in the atomic data nor an ionization-stratification effect, but something more fundamental. Each diagnostic is tuned to a particular density, and when a nebula contains gas across a wide range of densities, as real nebulae do, each line reports the part of that range it is most sensitive to. The diagnostics do not measure a representative average density; they respond to different parts of a broad density distribution. This resolves the discrepancy with a simple relation between the density each line returns and the density it is most sensitive to, a relation that holds from individual H II regions to whole galaxies, near and far, and reveals that ionized nebulae contain far more dense gas than any one diagnostic implies. A nebula has no single electron density to measure, but a broad density distribution, and the masses, pressures, abundances and energetics built on the single-density assumption must be reconsidered, from nearby star-forming regions to galaxies across cosmic time.

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

98 extracted references · 98 canonical work pages · 2 internal anchors

  1. [1]

    Seaton, M. J. & Osterbrock, D. E. Relative [O II] Intensities in Gaseous Nebulae. ApJ125, 66 (1957)

  2. [2]

    Planetary Nebulae II

    Peimbert, M. Planetary Nebulae II. Electron Temperatures and Electron Densities.Boletin de los Observatorios Tonantzintla y Tacubaya6, 29–37 (1971)

  3. [3]

    Osterbrock, D. E. & Ferland, G. J.Astrophysics of gaseous nebulae and active galactic nuclei(2006)

  4. [4]

    & Rodr´ ıguez, M

    Juan de Dios, L. & Rodr´ ıguez, M. Atomic data and the density structures of planetary nebulae. MNRAS507, 5331–5339 (2021)

  5. [5]

    A.et al.Characterizing Extreme Emission-line Galaxies

    Berg, D. A.et al.Characterizing Extreme Emission-line Galaxies. I. A Four-zone Ionization Model for Very High-ionization Emission. ApJ922, 170 (2021)

  6. [6]

    Exploring UV Diagnostics of the Interstellar Medium in Local High-z Analogs at the Dawn of the JWST Era

    Mingozzi, M.et al.CLASSY IV. Exploring UV Diagnostics of the Interstellar Medium in Local High-z Analogs at the Dawn of the JWST Era. ApJ939, 110 (2022)

  7. [7]

    ApJ995, 204 (2025)

    Martinez, Z.et al.Under Pressure: Decoding the Effect of High Densities on Derived Nebular Properties. ApJ995, 204 (2025)

  8. [8]

    Rubin, R. H. The Effect of Density Variations on Elemental Abundance Ratios in Gaseous Nebulae. ApJS69, 897 (1989)

  9. [9]

    Copetti, M. V. F. & Writzl, B. C. Study of electron density in planetary nebulae. A comparison of different density indicators. A&A382, 282–290 (2002). 21 5h36m00s 35m30s 00s 34m30s 00s 5°10' 20' 30' 40' RA DEC [S II] 6717/ 6731 1.0 1.5 2.0 2.5 3.0 3.5 4.0 log ne [cm 3] 5h36m00s 35m30s 00s 34m30s 00s 5°10' 20' 30' 40' RA DEC [O II] 3727/ 3729 1.0 1.5 2.0 ...

  10. [10]

    & Barlow, M

    Wang, W., Liu, X.-W., Zhang, Y. & Barlow, M. J. A reexamination of electron density diagnostics for ionized gaseous nebulae. A&A427, 873–886 (2004)

  11. [11]

    J.et al.Investigating the Drivers of Electron Temperature Variations in H II Regions with Keck-KCWI and VLT-MUSE

    Rickards Vaught, R. J.et al.Investigating the Drivers of Electron Temperature Variations in H II Regions with Keck-KCWI and VLT-MUSE. ApJ966, 130 (2024)

  12. [12]

    E., Esteban, C., Garc´ ıa-Rojas, J., Kreckel, K

    M´ endez-Delgado, J. E., Esteban, C., Garc´ ıa-Rojas, J., Kreckel, K. & Peimbert, M. Effects of density and temperature variations on the metallicity of Mrk 71. Nature Astronomy8, 275–277 (2024)

  13. [13]

    J., Nicholls, D

    Kewley, L. J., Nicholls, D. C. & Sutherland, R. S. Understanding Galaxy Evolution Through Emission Lines. ARA&A57, 511–570 (2019)

  14. [14]

    Steidel, C. C.et al.Strong Nebular Line Ratios in the Spectra of z ˜2-3 Star 23 2 3 4 5 6 log ne [cm 3] 0 1 2 3 4Probability density Input density PDF [S II] 6717/6731 [O II] 3727/3729 [Cl III] 5538/5518 [Ar IV] 4740/4711 [Fe III] 4658/4702 C III] 1909/1907 Si III] 1882/1892 [N IV] 1483/1487 Fig. 6 Distributions of inferred electron densities from the for...

  15. [15]

    L.et al.The KMOS 3D Survey: Investigating the Origin of the Elevated Electron Densities in Star-forming Galaxies at 1≲z≲3

    Davies, R. L.et al.The KMOS 3D Survey: Investigating the Origin of the Elevated Electron Densities in Star-forming Galaxies at 1≲z≲3. ApJ909, 78 (2021)

  16. [16]

    The impact of spectra quality on nebular abundances

    Rodr´ ıguez, M. The impact of spectra quality on nebular abundances. MNRAS 495, 1016–1034 (2020)

  17. [17]

    Filippenko, A. V. New evidence for photoionization as the dominant excitation mechanism in liners. ApJ289, 475–489 (1985)

  18. [18]

    & Verner, D

    Baldwin, J., Ferland, G., Korista, K. & Verner, D. Locally Optimally Emitting Clouds and the Origin of Quasar Emission Lines. ApJ455, L119 (1995)

  19. [19]

    W., Korista, K

    Ferguson, J. W., Korista, K. T., Baldwin, J. A. & Ferland, G. J. Locally Optimally Emitting Clouds and the Narrow Emission Lines in Seyfert Galaxies. ApJ487, 122–141 (1997)

  20. [20]

    O’Dell, C. R. The Orion Nebula and its Associated Population. ARA&A39, 99–136 (2001)

  21. [21]

    7 Narrow-PDF model validation.Probability density distributions of log 10 ne recovered by each diagnostic from a forward model identical to that of Extended Data Fig

    Drory, N.et al.The SDSS-V Local Volume Mapper (LVM): Scientific Motivation 24 1.5 2.0 2.5 3.0 3.5 4.0 log ne [cm 3] 0 2 4 6 8 10Probability density Input density PDF [S II] 6717/6731 [O II] 3727/3729 [Cl III] 5538/5518 [Ar IV] 4740/4711 [Fe III] 4658/4702 C III] 1909/1907 Si III] 1882/1892 [N IV] 1483/1487 Fig. 7 Narrow-PDF model validation.Probability de...

  22. [22]

    A.et al.Sloan Digital Sky Survey

    Kollmeier, J. A.et al.Sloan Digital Sky Survey. V. Pioneering Panoptic Spectroscopy. AJ171, 52 (2026)

  23. [23]

    A&A689, A352 (2024)

    Kreckel, K.et al.SDSS-V Local Volume Mapper (LVM): A glimpse into Orion. A&A689, A352 (2024)

  24. [24]

    E.et al.Density biases and temperature relations for DESIRED H II regions

    M´ endez-Delgado, J. E.et al.Density biases and temperature relations for DESIRED H II regions. MNRAS523, 2952–2973 (2023)

  25. [25]

    Z.et al.A self-consistent direct method for chemical abundances in high-z galaxies with JWST

    Arellano-C´ ordova, K. Z.et al.A self-consistent direct method for chemical abundances in high-z galaxies with JWST. MNRAS547, stag380 (2026)

  26. [26]

    & Shaw, R

    Luridiana, V., Morisset, C. & Shaw, R. A. PyNeb: a new tool for analyzing emission lines. I. Code description and validation of results. A&A573, A42 (2015)

  27. [27]

    Garc´ ıa-V´ azquez, J., Henney, W. J. & Casta˜ neda, H. O. Turbulence in compact to giant H II regions. MNRAS523, 4202–4226 (2023)

  28. [28]

    Schneider, N.et al.Detection of two power-law tails in the probability distribution 25 = 1.32+0.30 0.18 1.5 0.0 1.5 log10 (nmin) log10 (nmin) = 1.47+0.44 1.22 2 1 0 4 6 8 10 log10 (nmax) 1.5 0.0 1.5log10 (nmin) 4 6 8 10 log10 (nmax) log10 (nmax) = 6.25+2.25 1.68 Uninformative prior: [ 2, 0.5] log10 (nmin) [ 3, 3.5] log10 (nmax/nmin) [0.1, 8] a = 1.41+0.13...

  29. [29]

    H.-H., Burkhart, B., Goodman, A

    Chen, H. H.-H., Burkhart, B., Goodman, A. & Collins, D. C. The Anatomy of the Column Density Probability Distribution Function (N-PDF). ApJ859, 162 (2018)

  30. [30]

    & Klessen, R

    Federrath, C. & Klessen, R. S. On the Star Formation Efficiency of Turbulent Magnetized Clouds. ApJ763, 51 (2013)

  31. [31]

    Klessen, R. S. & Glover, S. C. O. Physical Processes in the Interstellar Medium. Saas-Fee Advanced Course43, 85 (2016)

  32. [32]

    & Preite-Martinez, A

    Panagia, N. & Preite-Martinez, A. Forbidden lines from inhomogeneous nebulae: the case of a density distribution N∝r −2. A&A43, 93–100 (1975)

  33. [33]

    Aller, L. H. & Epps, H. W. Electron densities in gaseous nebulae. ApJ204, 445 (1976). 26 3.0 3.5 4.0 4.5 5.0 5.5 6.0 6.5 7.0 log n [cm 3] 3.0 3.5 4.0 4.5 5.0 5.5inferred log ne [cm 3] HH 202-S HH 204 1 2 34 56 7 8 9 10 11 Diagnostic

  34. [34]

    [S II] 4070/(6717+6731)

  35. [35]

    [O II] (7319+7320)/(3727+3729)

  36. [36]

    [Fe III] 4702/4881 10.[Fe III] 4658/4881 11.[Fe III] 4658/4702 Fig. 9 Observational counterpart of the narrow-PDF test: the hierarchy vanishes in a kinematically isolated dense component.Inferred electron density as a function of the maximum-sensitivity densityn M for the photoionized Herbig-Haro objects HH 202-S (black) and HH 204 (blue) in the Orion Neb...

  37. [37]

    Viegas, S. M. & Clegg, R. E. S. Density condensations in planetary nebulae and the electron temperature. MNRAS271, 993–998 (1994)

  38. [38]

    & Szczerba, R

    Stasi´ nska, G. & Szczerba, R. The temperature structure of dusty planetary nebulae. A&A379, 1024–1038 (2001)

  39. [39]

    Temperature Determinations of H II Regions

    Peimbert, M. Temperature Determinations of H II Regions. ApJ150, 825 (1967)

  40. [40]

    E., Esteban, C., Garc´ ıa-Rojas, J., Kreckel, K

    M´ endez-Delgado, J. E., Esteban, C., Garc´ ıa-Rojas, J., Kreckel, K. & Peimbert, M. Temperature inhomogeneities cause the abundance discrepancy in H II regions. Nature618, 249–251 (2023)

  41. [41]

    M.et al.The SDSS-V Local Volume Mapper Telescope System

    Herbst, T. M.et al.The SDSS-V Local Volume Mapper Telescope System. AJ 168, 267 (2024)

  42. [42]

    F.et al.The SDSS-V Local Volume Mapper: Data Analysis Pipeline

    S´ anchez, S. F.et al.The SDSS-V Local Volume Mapper: Data Analysis Pipeline. AJ169, 52 (2025)

  43. [43]

    E.et al.Gradients of chemical abundances in the Milky Way from H II regions: distances derived from Gaia EDR3 parallaxes and temperature inhomogeneities

    M´ endez-Delgado, J. E.et al.Gradients of chemical abundances in the Milky Way from H II regions: distances derived from Gaia EDR3 parallaxes and temperature inhomogeneities. MNRAS510, 4436–4455 (2022)

  44. [44]

    Orozco-Duarte, R.et al.SDSS-V LVM: Revealing the Physical and Chemical 27 Structure of the Helix Nebula.arXiv e-printsarXiv:2606.11394 (2026)

  45. [45]

    A&A706, A81 (2026)

    Sattler, N.et al.SDSS-V LVM: Resolving physical conditions in the Trifid Nebula. A&A706, A81 (2026)

  46. [46]

    A., Clayton, G

    Cardelli, J. A., Clayton, G. C. & Mathis, J. S. The Relationship between Infrared, Optical, and Ultraviolet Extinction. ApJ345, 245 (1989)

  47. [47]

    Blagrave, K. P. M.et al.Deviations from He I Case B Recombination Theory and Extinction Corrections in the Orion Nebula. ApJ655, 299–315 (2007)

  48. [48]

    Storey, P. J. & Hummer, D. G. Recombination line intensities for hydrogenic ions-IV. Total recombination coefficients and machine-readable tables for Z=1 to

  49. [49]

    MNRAS272, 41–48 (1995)

  50. [50]

    J.et al.Pumping up the [N I] Nebular Lines

    Ferland, G. J.et al.Pumping up the [N I] Nebular Lines. ApJ757, 79 (2012)

  51. [51]

    Fluorescence of [Fe II] in H II regions

    Rodr´ ıguez, M. Fluorescence of [Fe II] in H II regions. A&A348, 222–226 (1999)

  52. [52]

    M., Verner, D

    Verner, E. M., Verner, D. A., Baldwin, J. A., Ferland, G. J. & Martin, P. G. Continuum Pumping of [Fe II] in the Orion Nebula. ApJ543, 831–839 (2000)

  53. [53]

    A.et al.Physical Conditions in Low Ionization Regions of the Orion Nebula

    Baldwin, J. A.et al.Physical Conditions in Low Ionization Regions of the Orion Nebula. ApJ468, L115 (1996)

  54. [54]

    E.et al.Photoionized Herbig-Haro Objects in the Orion Nebula through Deep High Spectral Resolution Spectroscopy

    M´ endez-Delgado, J. E.et al.Photoionized Herbig-Haro Objects in the Orion Nebula through Deep High Spectral Resolution Spectroscopy. II. HH 204. ApJ 918, 27 (2021)

  55. [55]

    E., Bautista, M., Garc´ ıa-Rojas, J

    Mendoza, C., M´ endez-Delgado, J. E., Bautista, M., Garc´ ıa-Rojas, J. & Moris- set, C. Atomic Data Assessment with PyNeb: Radiative and Electron Impact Excitation Rates for [Fe II] and [Fe III].Atoms11, 63 (2023)

  56. [56]

    A., Skillman, E

    Berg, D. A., Skillman, E. D., Henry, R. B. C., Erb, D. K. & Carigi, L. Carbon and Oxygen Abundances in Low Metallicity Dwarf Galaxies. ApJ827, 126 (2016)

  57. [57]

    M.et al.Hubble Space Telescope Observations of the Dusty Small Magellanic Cloud H II Region N88A

    Kurt, C. M.et al.Hubble Space Telescope Observations of the Dusty Small Magellanic Cloud H II Region N88A. ApJ518, 246–261 (1999)

  58. [58]

    The Sunburst Arc with JWST: III. An Abundance of Direct Chemical Abundances

    Welch, B.et al.The Sunburst Arc with JWST: III. An Abundance of Direct Chemical Abundances.arXiv e-printsarXiv:2405.06631 (2024)

  59. [59]

    W.et al.Metal-poor star formation at z ¿ 6 with JWST: new insight into hard radiation fields and nitrogen enrichment on 20 pc scales

    Topping, M. W.et al.Metal-poor star formation at z ¿ 6 with JWST: new insight into hard radiation fields and nitrogen enrichment on 20 pc scales. MNRAS529, 3301–3322 (2024)

  60. [60]

    E.et al.Gas-phase Fe/O and Fe/N abundances in star- forming regions: Relations between nucleosynthesis, metallicity, and dust

    M´ endez-Delgado, J. E.et al.Gas-phase Fe/O and Fe/N abundances in star- forming regions: Relations between nucleosynthesis, metallicity, and dust. A&A 28 690, A248 (2024)

  61. [61]

    ApJS153, 501–522 (2004)

    Garc´ ıa-Rojas, J.et al.Chemical Abundances of the Galactic H II Region NGC 3576 Derived from Very Large Telescope Echelle Spectrophotometry. ApJS153, 501–522 (2004)

  62. [62]

    MNRAS362, 301–312 (2005)

    Garc´ ıa-Rojas, J.et al.Deep echelle spectrophotometry of S 311, a Galactic HII region located outside the solar circle. MNRAS362, 301–312 (2005)

  63. [63]

    MNRAS368, 253–279 (2006)

    Garc´ ıa-Rojas, J.et al.Faint emission lines in the Galactic HII regions M16, M20 and NGC 3603 ∗. MNRAS368, 253–279 (2006)

  64. [64]

    A revision based on deep VLT echelle spectrophotometry

    Garc´ ıa-Rojas, J.et al.The chemical composition of the Galactic H II regions M8 and M17. A revision based on deep VLT echelle spectrophotometry. Rev. Mexicana Astron. Astrofis.43, 3–31 (2007)

  65. [65]

    & Rodr´ ıguez, M

    Esteban, C., Peimbert, M., Torres-Peimbert, S. & Rodr´ ıguez, M. Optical Recom- bination Lines of Heavy Elements in Giant Extragalactic H II Regions. ApJ581, 241–257 (2002)

  66. [66]

    ApJ700, 654–678 (2009)

    Esteban, C.et al.Keck HIRES Spectroscopy of Extragalactic H II Regions: C and O Abundances from Recombination Lines. ApJ700, 654–678 (2009)

  67. [67]

    MNRAS433, 382–393 (2013)

    Esteban, C.et al.NGC 2579 and the carbon and oxygen abundance gradients beyond the solar circle. MNRAS433, 382–393 (2013)

  68. [68]

    & Toribio San Cipriano, L

    Esteban, C., Fang, X., Garc´ ıa-Rojas, J. & Toribio San Cipriano, L. The radial abundance gradient of oxygen towards the Galactic anti-centre. MNRAS471, 987–1004 (2017)

  69. [69]

    & Garc´ ıa-Rojas, J

    Esteban, C. & Garc´ ıa-Rojas, J. Revisiting the radial abundance gradients of nitrogen and oxygen of the Milky Way. MNRAS478, 2315–2336 (2018)

  70. [70]

    & Toribio San Cipriano, L

    Esteban, C., Bresolin, F., Garc´ ıa-Rojas, J. & Toribio San Cipriano, L. Carbon, nitrogen, and oxygen abundance gradients in M101 and M31. MNRAS491, 2137–2155 (2020)

  71. [71]

    & Fierro, J

    Torres-Peimbert, S., Peimbert, M. & Fierro, J. Physical Conditions of H II Regions in M101 and the Pregalactic Helium Abundance. ApJ345, 186 (1989)

  72. [72]

    & Bresolin, F

    Zurita, A. & Bresolin, F. The chemical abundance in M31 from H II regions. MNRAS427, 1463–1481 (2012)

  73. [73]

    & Toribio San Cipriano, L

    Dom´ ınguez-Guzm´ an, G., Rodr´ ıguez, M., Garc´ ıa-Rojas, J., Esteban, C. & Toribio San Cipriano, L. The homogeneity of chemical abundances in H II regions of the Magellanic Clouds. MNRAS517, 4497–4514 (2022). 29

  74. [74]

    & Sixtos, A

    Valerdi, M., Peimbert, A., Peimbert, M. & Sixtos, A. Determination of the Pri- mordial Helium Abundance Based on NGC 346, an H II Region of the Small Magellanic Cloud. ApJ876, 98 (2019)

  75. [75]

    A., Peimbert, A., Peimbert, M

    Pe˜ na-Guerrero, M. A., Peimbert, A., Peimbert, M. & Ruiz, M. T. Analysis of Two Small Magellanic Cloud H II Regions Considering Thermal Inhomogeneities: Implications for the Determinations of Extragalactic Chemical Abundances. ApJ 746, 115 (2012)

  76. [76]

    The Chemical Composition of the 30 Doradus Nebula Derived from Very Large Telescope Echelle Spectrophotometry

    Peimbert, A. The Chemical Composition of the 30 Doradus Nebula Derived from Very Large Telescope Echelle Spectrophotometry. ApJ584, 735–750 (2003)

  77. [77]

    R., Esteban, C., Garc´ ıa-Rojas, J., Peimbert, M

    L´ opez-S´ anchez,´A. R., Esteban, C., Garc´ ıa-Rojas, J., Peimbert, M. & Rodr´ ıguez, M. The Localized Chemical Pollution in NGC 5253 Revisited: Results from Deep Echelle Spectrophotometry. ApJ656, 168–185 (2007)

  78. [78]

    Implications for chemical evolution

    Esteban, C.et al.Carbon and oxygen abundances from recombination lines in low-metallicity star-forming galaxies. Implications for chemical evolution. MNRAS443, 624–647 (2014)

  79. [79]

    F., P´ erez-Montero, E., D´ ıaz,´A

    H¨ agele, G. F., P´ erez-Montero, E., D´ ıaz,´A. I., Terlevich, E. & Terlevich, R. The temperature and ionization structure of the emitting gas in HII galaxies: impli- cations for the accuracy of abundance determinations. MNRAS372, 293–312 (2006)

  80. [80]

    F.et al.Precision abundance analysis of bright HII galaxies

    H¨ agele, G. F.et al.Precision abundance analysis of bright HII galaxies. MNRAS 383, 209–229 (2008)

Showing first 80 references.