A DESI Calibration of the [O II]--[S II] Electron-density Offset in Integrated Star-forming Galaxies
Pith reviewed 2026-06-29 03:26 UTC · model grok-4.3
The pith
Star-forming galaxies show [O II] electron densities 0.228 dex higher than [S II] on average.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
For star-forming galaxies with fiducial emission lines, [O II] yields systematically higher electron densities than [S II], with a median offset of 0.228 dex. The binned median calibration is log n_e(OII)=(0.752)log n_e(SII) +(0.832). The offset is larger in galaxies with higher stellar mass, Hα star-formation rate, dust attenuation, and N2, and smaller in galaxies with higher log O32. These trends are consistent with [O II] and [S II] sampling different low-ionization gas phases in integrated spectra, with [S II] more strongly weighted toward lower-density diffuse or outer gas.
What carries the argument
The empirical binned-median calibration between [O II]- and [S II]-derived electron densities, together with its trends against galaxy properties.
If this is right
- ISM pressure and nebular density studies must apply the calibration when comparing or combining [O II]- and [S II]-based values.
- Measurements of gas density evolution across redshift cannot mix the two diagnostics without the offset correction.
- Galaxy samples selected by different emission-line availability will carry systematic density biases unless calibrated.
- The phase weighting implies that [S II] densities better trace diffuse outer gas while [O II] includes denser inner regions.
Where Pith is reading between the lines
- High-redshift surveys that rely on [O II] because [S II] shifts out of the observed frame may report systematically higher densities than local [S II] studies.
- The mass and metallicity trends suggest the calibration slope could change in metal-poor or low-mass systems not well sampled by DESI.
- Resolved spectroscopy of individual galaxies could test whether the integrated offset disappears when both lines are measured in the same spatial aperture.
Load-bearing premise
The observed offset and trends arise from genuine differences in the low-ionization gas phases sampled by the two doublets rather than from survey selection, line-fitting systematics, or aperture effects.
What would settle it
Repeating the measurement in an independent catalog with identical aperture coverage and independent line fitting that finds a median offset consistent with zero would falsify the claimed density difference.
Figures
read the original abstract
The [O II]$\lambda\lambda3726,3729$ and [S II]$\lambda\lambda6716,6731$ doublets are widely used as low-ionization electron-density diagnostics in galaxy spectra and are often treated as interchangeable when only one of them is accessible. We test this assumption using the DESI DR1 Emission Line Catalog. For star-forming galaxies with fiducial emission lines, [O II] yields systematically higher electron densities than [S II], with a median offset of 0.228 dex. The binned median calibration is $\log n_e({\rm OII})=(0.752^{+0.182}_{-0.097})\log n_e({\rm SII}) +(0.832^{+0.231}_{-0.422})$. The offset is larger in galaxies with higher stellar mass, H$\alpha$ star-formation rate, dust attenuation, and $N2\equiv\log([{\rm NII}]\lambda6583/{\rm H}\alpha)$, an empirical gas-phase metallicity proxy, and smaller in galaxies with higher \(\log O_{32}\equiv\log\{[{\rm OIII}]\lambda5007/ ([{\rm OII}]\lambda3726+\lambda3729)\}\), an ionization proxy; no significant trend is found with specific star-formation rate. These trends are consistent with [O II] and [S II] sampling different low-ionization gas phases in integrated spectra, with [S II] more strongly weighted toward lower-density diffuse or outer gas. Our results show that [O II]- and [S II]-based densities should not be mixed without empirical calibration in studies of ISM pressure, nebular density, and their evolution across galaxy samples and redshift.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript uses the DESI DR1 Emission Line Catalog to derive an empirical calibration between electron densities measured from the [O II] λλ3726,3729 and [S II] λλ6716,6731 doublets in star-forming galaxies. It reports a median offset of 0.228 dex with [O II] yielding higher values, provides the binned-median relation log n_e(OII)=(0.752^{+0.182}_{-0.097}) log n_e(SII) + (0.832^{+0.231}_{-0.422}), and identifies trends of the offset with stellar mass, Hα SFR, dust attenuation, N2, and O32 (but not sSFR), interpreting these as evidence that the doublets sample different low-ionization gas phases.
Significance. If the offset and trends are shown to be physical rather than catalog artifacts, the result supplies a practical, data-driven calibration for mixing [O II]- and [S II]-based densities in studies of ISM pressure and nebular conditions across redshift. The large DESI sample enables statistically robust binned trends that go beyond a simple median offset.
major comments (2)
- [Abstract] Abstract (and § on sample selection, not numbered here): the central claim that the 0.228 dex offset and its correlations with M*, SFR, A_V, N2, and O32 reflect genuine differences in low-ionization phases requires that the DESI DR1 Emission Line Catalog measurements are free of differential biases in doublet deblending, sky subtraction, fiber aperture weighting, and star-forming selection cuts. No explicit null tests against fitting variants or selection variants are described; if such systematics dominate, the physical interpretation and the non-interchangeability conclusion do not hold.
- [Abstract] The binned calibration relation (abstract): the reported asymmetric uncertainties on slope and intercept are derived from the binned medians, but without the explicit binning scheme, number of galaxies per bin, or propagation of individual line-ratio uncertainties into the fit, it is unclear whether the quoted errors fully capture the scatter or selection effects that could correlate with the binned galaxy properties.
minor comments (1)
- [Abstract] The abstract uses 'fiducial emission lines' without definition; a brief parenthetical or reference to the catalog paper would improve clarity for readers.
Simulated Author's Rebuttal
We thank the referee for their careful review and constructive comments, which highlight important aspects of robustness and clarity. We address each major comment below and will revise the manuscript accordingly to strengthen the presentation and support for our conclusions.
read point-by-point responses
-
Referee: [Abstract] Abstract (and § on sample selection, not numbered here): the central claim that the 0.228 dex offset and its correlations with M*, SFR, A_V, N2, and O32 reflect genuine differences in low-ionization phases requires that the DESI DR1 Emission Line Catalog measurements are free of differential biases in doublet deblending, sky subtraction, fiber aperture weighting, and star-forming selection cuts. No explicit null tests against fitting variants or selection variants are described; if such systematics dominate, the physical interpretation and the non-interchangeability conclusion do not hold.
Authors: We agree that explicit tests for differential systematics are important to support the physical interpretation. The manuscript describes the DESI DR1 emission-line catalog, fitting pipeline, and star-forming selection (based on BPT diagram and line S/N cuts) but does not present dedicated null tests varying deblending, sky subtraction, or selection. We will add a new subsection in the methods or results detailing robustness checks, including: repeating the analysis with stricter S/N thresholds, alternative star-forming cuts, and examination of sky-subtraction residuals via negative-flux statistics. These additions will quantify any impact on the offset and trends. revision: yes
-
Referee: [Abstract] The binned calibration relation (abstract): the reported asymmetric uncertainties on slope and intercept are derived from the binned medians, but without the explicit binning scheme, number of galaxies per bin, or propagation of individual line-ratio uncertainties into the fit, it is unclear whether the quoted errors fully capture the scatter or selection effects that could correlate with the binned galaxy properties.
Authors: We thank the referee for noting the need for greater transparency on the fit. The calibration uses logarithmic bins in log n_e(SII) with a minimum occupancy of ~100 galaxies per bin; asymmetric uncertainties come from bootstrap resampling of the binned medians (16th/84th percentiles). The manuscript does not explicitly state the bin edges, per-bin counts, or discuss propagation of individual line-ratio errors. We will revise the methods section and relevant figure caption to include the binning details, report galaxy counts per bin, and clarify that the quoted errors reflect median scatter rather than formal propagation of flux uncertainties, while noting any limitations this introduces for correlated selection effects. revision: yes
Circularity Check
Empirical calibration from DESI data shows no circularity
full rationale
The paper measures electron densities from [O II] and [S II] doublets directly in the DESI DR1 Emission Line Catalog for star-forming galaxies, reports a median offset of 0.228 dex, and fits the binned-median relation log n_e(OII) = (0.752^{+0.182}_{-0.097}) log n_e(SII) + (0.832^{+0.231}_{-0.422}) as an empirical calibration. No derivation chain reduces any claimed result to a prior fit, self-citation, or input by construction; the trends with mass, SFR, dust, N2, and O32 are also observational measurements. The work is self-contained against external benchmarks with no load-bearing self-citations or ansatzes invoked.
Axiom & Free-Parameter Ledger
free parameters (2)
- slope =
0.752
- intercept =
0.832
axioms (1)
- domain assumption The [O II] λλ3726,3729 and [S II] λλ6716,6731 doublets serve as reliable low-ionization electron-density diagnostics through their flux ratios.
Forward citations
Cited by 1 Pith paper
-
Electron Densities of Typical Low-Mass Galaxies at z~2-7 from Stacked JWST/NIRSpec Spectra
Stacked spectra give n_e ~100-150 cm^{-3} at z=2-5 rising to 381 cm^{-3} at z=5-7 for low-mass galaxies, with fitted evolution n_e = n_e0 [(1+z)/(1+2.3)]^alpha.
Reference graph
Works this paper leans on
-
[1]
K., Ag¨ueros, M
Abazajian, K., Adelman-McCarthy, J. K., Ag¨ueros, M. A., et al. 2005, AJ, 129, 1755
2005
-
[2]
A., Phillips, M
Baldwin, J. A., Phillips, M. M., & Terlevich, R. 1981, PASP, 93, 5
1981
-
[3]
J., Dopita, M
Bian, F., Kewley, L. J., Dopita, M. A., & Juneau, S. 2016, ApJ, 822, 62
2016
-
[4]
2008, MNRAS, 385, 769
Brinchmann, J., Pettini, M., & Charlot, S. 2008, MNRAS, 385, 769
2008
-
[5]
A., Clayton, G
Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245 DESI Collaboration et al. 2026, AJ, 171, 285 DESI Collaboration et al. 2022, AJ, 164, 207
1989
-
[6]
L., F¨orster Schreiber, N
Davies, R. L., F¨orster Schreiber, N. M., Genzel, R., et al. 2021, ApJ, 909, 78
2021
-
[7]
J., & Gupta, A
Kaasinen, M., Bian, F., Groves, B., Kewley, L. J., & Gupta, A. 2017, MNRAS, 465, 3220
2017
-
[8]
J., Bian, F., et al
Kaasinen, M., Kewley, L. J., Bian, F., et al. 2018, MNRAS, 477, 5568
2018
-
[9]
2003, MNRAS, 346, 1055
Kauffmann, G., et al. 2003, MNRAS, 346, 1055
2003
-
[10]
C., & Evans, N
Kennicutt, R. C., & Evans, N. J. 2012, ARA&A, 50, 531
2012
-
[11]
2001, ApJ, 556, 121
Trevena, J. 2001, ApJ, 556, 121
2001
-
[12]
J., Nicholls, D
Kewley, L. J., Nicholls, D. C., Sutherland, R., et al. 2019, ApJ, 880, 16
2019
-
[13]
Lacerda, E. A. D., Cid Fernandes, R., Couto, G. S., et al. 2018, MNRAS, 474, 3727
2018
-
[14]
2025, ApJL, 979, L13
Li, S., et al. 2025, ApJL, 979, L13
2025
-
[15]
Luridiana, V ., Morisset, C., & Shaw, R. A. 2015, A&A, 573, A42
2015
-
[16]
2021, MNRAS, 508, 1582 M´endez-Delgado, J
Mannucci, F., Belfiore, F., Curti, M., et al. 2021, MNRAS, 508, 1582 M´endez-Delgado, J. E., et al. 2023, MNRAS, 523, 2952
2021
-
[17]
E., & Ferland, G
Osterbrock, D. E., & Ferland, G. J. 2006, Astrophysics of Gaseous Nebulae and Active Galactic Nuclei, 2nd ed. (Sausalito, CA: University Science Books)
2006
-
[18]
2014, A&A, 561, A10
Proxauf, B., ¨Ottl, S., & Kimeswenger, S. 2014, A&A, 561, A10
2014
-
[19]
A., Sanders, R
Reddy, N. A., Sanders, R. L., Shapley, A. E., et al. 2023, ApJ, 951, 56
2023
-
[20]
L., Shapley, A
Sanders, R. L., Shapley, A. E., Kriek, M., et al. 2016, ApJ, 816, 23
2016
-
[21]
L., Shapley, A
Sanders, R. L., Shapley, A. E., Zhang, K., Yan, R. 2017, ApJ, 850, 136
2017
-
[22]
C., et al
Shimakawa, R., Kodama, T., Steidel, C. C., et al. 2015, MNRAS, 451, 1284 8 RONG ET AL. < 0.30 < 0.25 < 0.20 σ(log ne) cut [dex] 0.100 0.125 0.150 0.175 0.200 0.225 0.250 0.275median Δlog ne N=1938 N=700 N=162 a Stricter density precision 0.50-1.50 0.70-1.45 0.80-1.40 0.90-1.35 0.95-1.30 allowed [O II] λ3729/λ3726 0.075 0.100 0.125 0.150 0.175 0.200 0.225 ...
2015
-
[23]
C., Rudie, G
Steidel, C. C., Rudie, G. C., Strom, A. L., et al. 2014, ApJ, 795, 165
2014
-
[24]
L., Steidel, C
Strom, A. L., Steidel, C. C., Rudie, G. C., et al. 2017, ApJ, 836, 164
2017
-
[25]
A., Heckman, T
Tremonti, C. A., Heckman, T. M., Kauffmann, G., et al. 2004, ApJ, 613, 898
2004
-
[26]
W., Sanders, R
Topping, M. W., Sanders, R. L., Shapley, A. E., et al. 2025, MNRAS, 541, 1707
2025
-
[27]
Wang, W., Liu, X.-W., Zhang, Y ., & Barlow, M. J. 2004, A&A, 427, 873
2004
-
[28]
G., Adelman, J., Anderson, J
York, D. G., Adelman, J., Anderson, J. E., Jr., et al. 2000, AJ, 120, 1579
2000
-
[29]
2017, MNRAS, 466, 3217
Zhang, K., Yan, R., Bundy, K., et al. 2017, MNRAS, 466, 3217
2017
-
[30]
2024, ApJ, 961, 173
Zou, H., et al. 2024, ApJ, 961, 173
2024
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.