Physical origin of the high-ionization temperature inhomogeneities causing the abundance discrepancy in H II regions, and the scope of the $t^2$ paradigm across ionized nebulae

Heavy-element abundances derived from collisionally excited lines (CELs) and recombination lines (RLs) in ionized nebulae disagree by factors of roughly 1.5–3 — the abundance discrepancy, quantified by the abundance discrepancy factor (ADF) — casting doubt on the absolute metallicity scale of H II regions, star-forming galaxies, and high-redshift systems. Méndez-Delgado et al. (Nature 618, 249, 2023) showed that in H II regions this discrepancy is caused by temperature inhomogeneities — described by Peimbert's $t^2$ parameter — concentrated in the high-ionization ($\mathrm{O}^{2+}$) gas, and gave empirical prescriptions correcting CEL-based abundances. This problem targets what that work left open: the physical mechanism generating the inhomogeneities (stellar feedback via winds and shocks and fluctuations in the hardness of the ionizing radiation are the candidates named in that work; turbulent/density structure is a further candidate raised by subsequent studies) has not been causally identified or quantitatively modeled; planetary nebulae do not follow the H II-region trends and their often extreme ADFs remain unexplained; the existence of a low-ionization abundance discrepancy in extragalactic H II regions is unverified; and the $t^2$-based corrections are still being validated on broad samples. A 2023–2024 exchange over the galaxy Mrk 71, where direct optical plus far-infrared [O III] diagnostics give $t^2\approx0.008$ versus the $\approx0.1$ required by the recombination lines, shows the paradigm itself is under active contest. As of mid-2026 the literature leaves all of these questions open.

Background

Ionized nebulae emit two families of lines from which chemical abundances are inferred. Collisionally excited lines (CELs), such as [O III] λ4363 and λ5007, are produced by electron-impact excitation of ions and depend exponentially on the electron temperature $T_e$; recombination lines (RLs), such as the O II V1 multiplet near λ4650, arise from electrons recombining with ions and are nearly temperature-insensitive. If the gas is thermally inhomogeneous, CEL-based abundances (which weight the hottest gas) are biased low relative to RL-based ones. Peimbert (1967) formalized this with the mean-square temperature fluctuation parameter $t^2=\langle(T_e-T_0)^2\rangle/T_0^2$, where $T_0$ is the volume-averaged electron temperature. The abundance discrepancy factor for an ion, ADF($X^{i+}$), is the ratio of the ionic abundance derived from RLs to that derived from CELs, commonly quoted as its logarithm in dex; it is ~0.2–0.5 dex in H II regions (zones of ionized interstellar gas around massive young stars) and can exceed 1 dex — reaching values above 10 in linear terms — in planetary nebulae (PNe), the ionized envelopes ejected by dying low- and intermediate-mass stars.

Because the widely used "direct method" for abundances — and the strong-line calibrations tied to it (empirical metallicity relations built on bright nebular line ratios, used when the faint temperature-sensitive lines are undetectable) — rests on CEL temperature diagnostics such as the ratio of the faint, temperature-sensitive auroral line [O III] λ4363 to the bright nebular line [O III] λ5007, the abundance discrepancy puts the absolute metallicity scale of the universe in doubt, from Galactic H II regions through the mass–metallicity relation of galaxies to high-redshift galaxies observed with the James Webb Space Telescope (JWST). Méndez-Delgado, Esteban, García-Rojas, Kreckel and Peimbert (Nature 618, 249, 2023) analyzed high-quality spectra of Galactic and extragalactic H II regions and found a tight relation between $t^2(\mathrm{O}^{2+})$ — the fluctuation parameter inferred by requiring RL and CEL abundances of $\mathrm{O}^{2+}$ to agree — and the temperature difference $\Delta T_e=T_e([\mathrm{O\,III}])-T_e([\mathrm{N\,II}])$ between the high- and low-ionization zones. They argued that temperature inhomogeneities confined to the high-ionization volume are the cause of the ADF in H II regions, discarding alternatives (chemically inhomogeneous gas, fluorescent excitation of the RLs, non-Maxwellian $\kappa$-distributed electrons, errors in recombination coefficients) for this object class, and provided empirical prescriptions (their Eqs. (3), (4), (5) and (7)) to correct CEL-based abundances even when only CELs are observed.

Problem Statement

Establish the physical origin of the temperature inhomogeneities that cause the abundance discrepancy in H II regions, and fix the scope of that explanation across photoionized nebulae. Concretely: (1) identify the mechanism(s) — candidates include stellar feedback through winds and shocks, fluctuations in the hardness of the ionizing stellar radiation, and turbulent or density structure of the gas — that generate temperature inhomogeneities concentrated in the high-ionization ($\mathrm{O}^{2+}$) volume, and demonstrate quantitatively, through causal physical modeling confronted with observation, that they reproduce the observed relation between $t^2(\mathrm{O}^{2+})$ and $\Delta T_e=T_e([\mathrm{O\,III}])-T_e([\mathrm{N\,II}])$, the observed ADF($\mathrm{O}^{2+}$) magnitudes and their increase toward low metallicity, and the spatially resolved ADF behavior now being mapped in nearby H II regions, together with the direct multi-wavelength constraints on $t^2$ (e.g., optical auroral lines combined with far-infrared [O III] 52/88 μm fine-structure lines, as applied to Mrk 71, where the same optical plus far-infrared data are contested between $t^2\approx0.008$ and $t^2\approx0.1$; a proposed mechanism must predict both of these $t^2$ diagnostics and state which of the two readings it supports, and identify the observation that would discriminate between them, rather than adopting one of the two values in advance); (2) determine whether the temperature-inhomogeneity mechanism identified in (1) or a distinct cause explains the abundance discrepancy in planetary nebulae — current evidence shows that planetary nebulae do not follow the H II-region trends and that the extreme-ADF objects remain unexplained — the candidates for a distinct cause being chemically distinct inclusions, binary-central-star gas components, and deviations from Case B recombination (the standard nebular assumption that photons from recombinations directly to the ground state are re-absorbed locally rather than escaping as observable lines), and identify the observation that would discriminate among these possibilities; and (3) establish observationally whether a significant low-ionization abundance discrepancy, such as ADF($\mathrm{O}^{+}$), exists in integrated spectra of extragalactic H II regions. For component (1), the mechanism must generate these relations rather than merely be compatible with them, and "reproduce" requires agreement, within the measurement uncertainties quoted in the cited works, with the slope and scatter of the $t^2(\mathrm{O}^{2+})$–$\Delta T_e$ relation, with the ADF($\mathrm{O}^{2+}$) magnitudes and their increase toward low metallicity, and with the spatially resolved ADF behavior, together with a determinate outcome for the contested Mrk 71 $t^2$ diagnostics as specified above. (1), (2) and (3) are jointly posed as a single objective — a physically grounded account of which nebular abundances are correct.

References

1. Méndez-Delgado, J. Eduardo; Esteban, César; García-Rojas, Jorge; Kreckel, Kathryn; Peimbert, Manuel (2023). Temperature inhomogeneities cause the abundance discrepancy in H ii regions. Nature 618, 249–251. DOI 10.1038/s41586-023-05956-2, https://doi.org/10.1038/s41586-023-05956-2 (preprint arXiv:2305.11578, https://arxiv.org/abs/2305.11578)
2. Peimbert, Manuel (1967). Temperature Determinations of H II Regions. The Astrophysical Journal 150, 825. DOI 10.1086/149385, https://doi.org/10.1086/149385
3. Chen, Yuguang; Jones, Tucker; Sanders, Ryan; Fadda, Dario; Sutter, Jessica; Minchin, Robert; Huntzinger, Erin; Senchyna, Peter; Stark, Daniel; Spilker, Justin; Weiner, Benjamin; Roberts-Borsani, Guido (2023). Accurate oxygen abundance of interstellar gas in Mrk 71 from optical and infrared spectra. Nature Astronomy 7(7), 771–778. DOI 10.1038/s41550-023-01953-7, https://doi.org/10.1038/s41550-023-01953-7 (preprint arXiv:2304.09898, https://arxiv.org/abs/2304.09898)
4. Méndez-Delgado, J. Eduardo; Esteban, César; García-Rojas, Jorge; Kreckel, Kathryn; Peimbert, Manuel (2024). Effects of density and temperature variations on the metallicity of Mrk 71. Nature Astronomy 8(3), 275–277. DOI 10.1038/s41550-024-02198-8, https://doi.org/10.1038/s41550-024-02198-8 (preprint arXiv:2310.01197, https://arxiv.org/abs/2310.01197)
5. Chen, Yuguang; Jones, Tucker; Sanders, Ryan; Fadda, Dario; Sutter, Jessica; Minchin, Robert; Huntzinger, Erin; Senchyna, Peter; Spilker, Justin; Weiner, Benjamin; Roberts-Borsani, Guido (2024). Reply to: Effects of density and temperature variations on the metallicity of Mrk71. Nature Astronomy 8(3), 278–279. DOI 10.1038/s41550-024-02199-7, https://doi.org/10.1038/s41550-024-02199-7 (preprint arXiv:2402.03105, https://arxiv.org/abs/2402.03105)
6. Rickards Vaught, Ryan J.; Sandstrom, Karin M.; Belfiore, Francesco; Kreckel, Kathryn; Méndez-Delgado, J. Eduardo; Emsellem, Eric; Groves, Brent; Blanc, Guillermo A.; Dale, Daniel A.; Egorov, Oleg V.; Glover, Simon C. O.; Grasha, Kathryn; Klessen, Ralf S.; Neumann, Justus; Williams, Thomas G. (2024). Investigating the Drivers of Electron Temperature Variations in H ii Regions with Keck-KCWI and VLT-MUSE. The Astrophysical Journal 966(1), 130. DOI 10.3847/1538-4357/ad303c, https://doi.org/10.3847/1538-4357/ad303c (preprint arXiv:2309.17440, https://arxiv.org/abs/2309.17440)
7. Méndez-Delgado, J. E.; Skillman, E. D.; Aver, E.; Morisset, C.; Esteban, C.; García-Rojas, J.; Kreckel, K.; Rogers, N. S. J.; Rosales-Ortega, F. F.; Arellano-Córdova, K. Z.; Flury, S. R.; Reyes-Rodríguez, E.; Orte-García, M.; Tan, S. (2025). Generalized Te([O III])–Te(He I) Discrepancies in Ionized Nebulae: Possible Evidence of Case B Deviations and Temperature Inhomogeneities. The Astrophysical Journal 986(1), 74. DOI 10.3847/1538-4357/adc67a, https://doi.org/10.3847/1538-4357/adc67a (preprint arXiv:2410.17381, https://arxiv.org/abs/2410.17381)
8. Singh, Amrita; Blanc, Guillermo A.; Kumari, Nimisha; Méndez-Delgado, J. E.; Sánchez, Sebastián F.; Morisset, Christophe; Congiu, Enrico; Kreckel, Kathryn; Roman-Lopes, Alexandre; Egorov, Oleg; Drory, Niv; Sankrit, Ravi; Mejía-Narváez, Alfredo; Egorova, Evgeniya; Jones, Amy M.; Bizyaev, Dmitry; Sattler, Natascha; Johnston, Evelyn J.; Minniti, Dante; Zermeño, Rodolfo de J.; Fernández-Trincado, José G.; Kollmeier, Juna A. (2026). SDSS-V LVM: A Spatially Resolved Study of the Physical Conditions and the Chemical Abundance Discrepancy in the Lagoon Nebula (M 8). The Astrophysical Journal 1001, 238. DOI 10.3847/1538-4357/ae563a, https://doi.org/10.3847/1538-4357/ae563a (preprint arXiv:2603.22382, https://arxiv.org/abs/2603.22382)
9. Méndez-Delgado, J. E.; García-Rojas, J. (2023). The abundance discrepancy in ionized nebulae: which are the correct abundances? Proceedings of the International Astronomical Union 19, 370–378. DOI 10.1017/s174392132400005x, https://doi.org/10.1017/s174392132400005x (preprint arXiv:2311.10280, https://arxiv.org/abs/2311.10280)
10. Gómez-Llanos, V.; García-Rojas, J.; Morisset, C.; Monteiro, H.; Jones, D.; Wesson, R.; Boffin, H. M. J.; Corradi, R. L. M. (2024). MUSE spectroscopy of the high abundance discrepancy planetary nebula NGC 6153. Astronomy & Astrophysics 689, A228. DOI 10.1051/0004-6361/202450822, https://doi.org/10.1051/0004-6361/202450822 (preprint arXiv:2407.06385, https://arxiv.org/abs/2407.06385)
11. Kwok, Sun; Balick, Bruce; Chu, You-Hua; Hrivnak, Bruce J.; López, Alberto; Parker, Quentin; Sahai, Raghvendra; Zijlstra, Albert (2026). Current Unsolved Problems in Planetary Nebulae Research. Galaxies 14(2), 30. DOI 10.3390/galaxies14020030, https://doi.org/10.3390/galaxies14020030 (preprint arXiv:2604.22999, https://arxiv.org/abs/2604.22999)
12. Rosales-Ortega, F. F.; Méndez-Delgado, J. E.; Guerrero-González, J. U.; Esteban, C.; García-Rojas, J.; Arellano-Córdova, K. Z.; Lugo-Aranda, A. Z.; Espíndola-Camacho, O.; López-Gutiérrez, J. C.; Martínez-Rivero, L. E.; Morisset, C.; Orte-García, M.; Reyes-Rodríguez, E.; Toribio San Cipriano, L.; Kreckel, K.; Egorov, O.; Zinchenko, I. A.; Sánchez, S. F.; Vílchez, J. M. (2026). The DESIRED strong-line calibrations: I. New empirical metallicity relations for the local and high-redshift universe. arXiv:2604.16273, https://arxiv.org/abs/2604.16273
13. Xing, Larrance; Choustnikov, Nicholas; Katz, Harley; Cameron, Alex J. (2026). Inferring Interstellar Medium Density, Temperature, and Metallicity from Turbulent H II Regions. Open Journal of Astrophysics 9. DOI 10.33232/001c.156223, https://doi.org/10.33232/001c.156223 (preprint arXiv:2602.02949, https://arxiv.org/abs/2602.02949)
14. Nemer, Ahmad; Mendez-Delgado, J. E.; Sattler, Natascha; Blanc, Guillermo A.; Singh, Amrita; Kreckel, Kathryn; Gelfand, Joseph D.; Drory, Niv (2024). Analyzing the Abundance Discrepancy Problem in HII Regions with Photoionization Modeling. arXiv:2410.20819, https://arxiv.org/abs/2410.20819
15. Habjan, Eric; Faesi, Christopher; Kreckel, Kathryn; Méndez-Delgado, J. Eduardo; Belfiore, Francesco; Vaught, Ryan J.; Groves, Brent; Scheuermann, Fabian; Williams, Thomas G.; Klessen, Ralf S.; Amiri, Amirnezam; Grasha, Kathryn; Glover, Simon (2026). Toward Unbiased Abundance Measurements in Inhomogeneous H II Regions. arXiv:2607.05295, https://arxiv.org/abs/2607.05295
16. Méndez-Delgado, J. Eduardo; Morisset, Christophe; Henney, William J.; Drory, Niv; Egorov, Oleg V.; Sánchez, Sebastián F.; Kollmeier, Juna A.; Kreckel, Kathryn; Blanc, Guillermo; Stasińska, Grażyna; Johnston, Evelyn J.; Ibarra-Medel, Héctor J.; Mejía-Narváez, Alfredo J.; Esteban, César; García-Rojas, Jorge; Singh, Amrita; Katkov, Ivan Yu.; Skillman, Evan D.; Orozco-Duarte, Rogelio; Zinchenko, Igor A.; Lugo-Aranda, Alejandra Z.; Wofford, Aida; Glover, Simon C. O.; Egorova, Evgeniya; Zermeño, Rodolfo de J.; Castañeda-Carlos, Lesly C.; Liang, Fu-Heng; Sattler, Natascha; Cruz-González, Irene; Brownstein, Joel R.; Hilder, Thomas; Schneider, Donald P. (2026). There is no single density: star-forming regions and galaxies hold more dense ionized gas than long assumed. arXiv:2607.07973, https://arxiv.org/abs/2607.07973
17. Amayo, Alexia; Carigi, Leticia; Méndez-Delgado, José Eduardo; García-Rojas, Jorge; Reyes-Rodríguez, Elena (2026). A comparative study of O, Ne, Cl, and Ar in Hii regions and PNe of the Galactic disk: Temporal evolution of radial gradients? arXiv:2606.19061, https://arxiv.org/abs/2606.19061