pith. sign in

arxiv: 2606.28129 · v1 · pith:YLJI4VUPnew · submitted 2026-06-26 · 🌌 astro-ph.GA

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

classification 🌌 astro-ph.GA
keywords electron density[O II] doublet[S II] doubletstar-forming galaxiesDESInebular diagnosticsISM gas phasesdensity calibration
0
0 comments X

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.

The paper tests the common assumption that the [O II] and [S II] doublets can be treated as interchangeable electron-density diagnostics in galaxy spectra. Using the DESI DR1 Emission Line Catalog for star-forming galaxies with fiducial lines, it measures a median offset where [O II] gives higher densities. It supplies an empirical calibration relating the two: log n_e(OII) equals roughly 0.75 times log n_e(SII) plus 0.83. The size of the offset changes with stellar mass, star-formation rate, dust, a metallicity proxy, and an ionization proxy. The trends point to the two lines weighting different low-ionization gas phases within the same integrated spectrum.

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

These are editorial extensions of the paper, not claims the author makes directly.

  • 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

Figures reproduced from arXiv: 2606.28129 by Hu Zou, Shihong Liu, Yu Rong.

Figure 1
Figure 1. Figure 1: Same-galaxy comparison of DESI and SDSS spectra around the [O II]λλ3726, 3729 doublet. Black histograms show DESI coadd spectra, gray histograms show SDSS spectra, and orange curves show simple DESI local continuum-plus-doublet fits. For visual comparison only, each spectrum is locally continuum-subtracted and normalized by its [O II] peak. Vertical dotted lines indicate the observed-frame centers of the t… view at source ↗
Figure 2
Figure 2. Figure 2: The DESI [O II]–[S II] density plane for the fiducial σ(log ne) < 0.3 dex subset. Each cell is colored by the median value of the indicated galaxy property; green points show individual galaxies, black points with error bars show the binned median density relation, and the orange line is the median calibration fit. The black points are medians in equal-number bins of log ne(SII); the horizontal and vertica… view at source ↗
Figure 3
Figure 3. Figure 3: Galaxy-property trends with the [S II]-based electron density, shown as a control comparison for [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Direct correlations between galaxy properties and the [O II]–[S II] density offset. For each property X, we plot X against ∆ log ne = log ne(OII) − log ne(SII). Small green points show individual galaxies within the displayed plotting range; orange points and curves show equal-number binned medians. The title of each panel gives the Spearman rank correlation coefficient, rs, computed from the full fiducial… view at source ↗
Figure 5
Figure 5. Figure 5: Robustness of the [O II]–[S II] density offset. Panel (a) shows the median offset for the fiducial and stricter density-precision cuts; labels give the retained number of galaxies. Panel (b) shows the median offset after increasingly conservative [O II] ratio cuts. Panel (c) repeats the exercise for [S II] ratio cuts. Panel (d) shows the median offset after recomputing both diagnostics at different fixed e… view at source ↗
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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

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)
  1. [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.
  2. [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)
  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

2 responses · 0 unresolved

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
  1. 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

  2. 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

0 steps flagged

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

2 free parameters · 1 axioms · 0 invented entities

The result rests on an empirical fit whose coefficients are free parameters; it invokes the standard domain assumption that the [O II] and [S II] doublets function as density diagnostics via their line ratios.

free parameters (2)
  • slope = 0.752
    Coefficient in the binned-median linear relation between log n_e(OII) and log n_e(SII)
  • intercept = 0.832
    Coefficient in the binned-median linear relation between log n_e(OII) and log n_e(SII)
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.
    Invoked when treating the measured line ratios as direct density indicators.

pith-pipeline@v0.9.1-grok · 5861 in / 1649 out tokens · 66654 ms · 2026-06-29T03:26:35.480408+00:00 · methodology

discussion (0)

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

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Electron Densities of Typical Low-Mass Galaxies at z~2-7 from Stacked JWST/NIRSpec Spectra

    astro-ph.GA 2026-06 unverdicted novelty 7.0

    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

30 extracted references · cited by 1 Pith paper

  1. [1]

    K., Ag¨ueros, M

    Abazajian, K., Adelman-McCarthy, J. K., Ag¨ueros, M. A., et al. 2005, AJ, 129, 1755

  2. [2]

    A., Phillips, M

    Baldwin, J. A., Phillips, M. M., & Terlevich, R. 1981, PASP, 93, 5

  3. [3]

    J., Dopita, M

    Bian, F., Kewley, L. J., Dopita, M. A., & Juneau, S. 2016, ApJ, 822, 62

  4. [4]

    2008, MNRAS, 385, 769

    Brinchmann, J., Pettini, M., & Charlot, S. 2008, MNRAS, 385, 769

  5. [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

  6. [6]

    L., F¨orster Schreiber, N

    Davies, R. L., F¨orster Schreiber, N. M., Genzel, R., et al. 2021, ApJ, 909, 78

  7. [7]

    J., & Gupta, A

    Kaasinen, M., Bian, F., Groves, B., Kewley, L. J., & Gupta, A. 2017, MNRAS, 465, 3220

  8. [8]

    J., Bian, F., et al

    Kaasinen, M., Kewley, L. J., Bian, F., et al. 2018, MNRAS, 477, 5568

  9. [9]

    2003, MNRAS, 346, 1055

    Kauffmann, G., et al. 2003, MNRAS, 346, 1055

  10. [10]

    C., & Evans, N

    Kennicutt, R. C., & Evans, N. J. 2012, ARA&A, 50, 531

  11. [11]

    2001, ApJ, 556, 121

    Trevena, J. 2001, ApJ, 556, 121

  12. [12]

    J., Nicholls, D

    Kewley, L. J., Nicholls, D. C., Sutherland, R., et al. 2019, ApJ, 880, 16

  13. [13]

    Lacerda, E. A. D., Cid Fernandes, R., Couto, G. S., et al. 2018, MNRAS, 474, 3727

  14. [14]

    2025, ApJL, 979, L13

    Li, S., et al. 2025, ApJL, 979, L13

  15. [15]

    Luridiana, V ., Morisset, C., & Shaw, R. A. 2015, A&A, 573, A42

  16. [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

  17. [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)

  18. [18]

    2014, A&A, 561, A10

    Proxauf, B., ¨Ottl, S., & Kimeswenger, S. 2014, A&A, 561, A10

  19. [19]

    A., Sanders, R

    Reddy, N. A., Sanders, R. L., Shapley, A. E., et al. 2023, ApJ, 951, 56

  20. [20]

    L., Shapley, A

    Sanders, R. L., Shapley, A. E., Kriek, M., et al. 2016, ApJ, 816, 23

  21. [21]

    L., Shapley, A

    Sanders, R. L., Shapley, A. E., Zhang, K., Yan, R. 2017, ApJ, 850, 136

  22. [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 ...

  23. [23]

    C., Rudie, G

    Steidel, C. C., Rudie, G. C., Strom, A. L., et al. 2014, ApJ, 795, 165

  24. [24]

    L., Steidel, C

    Strom, A. L., Steidel, C. C., Rudie, G. C., et al. 2017, ApJ, 836, 164

  25. [25]

    A., Heckman, T

    Tremonti, C. A., Heckman, T. M., Kauffmann, G., et al. 2004, ApJ, 613, 898

  26. [26]

    W., Sanders, R

    Topping, M. W., Sanders, R. L., Shapley, A. E., et al. 2025, MNRAS, 541, 1707

  27. [27]

    Wang, W., Liu, X.-W., Zhang, Y ., & Barlow, M. J. 2004, A&A, 427, 873

  28. [28]

    G., Adelman, J., Anderson, J

    York, D. G., Adelman, J., Anderson, J. E., Jr., et al. 2000, AJ, 120, 1579

  29. [29]

    2017, MNRAS, 466, 3217

    Zhang, K., Yan, R., Bundy, K., et al. 2017, MNRAS, 466, 3217

  30. [30]

    2024, ApJ, 961, 173

    Zou, H., et al. 2024, ApJ, 961, 173