REVIEW 4 major objections 5 minor 48 references
A new technique for finding galaxies leaking Lyman-continuum radiation: [SII]-deficiency
T0 review · 4 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read Two of three [SII]-weak galaxies leak Lyman-continuum radiation
desk verdict A careful pilot with two real LyC detections, but the 'highly effective' claim outruns a sample of two and a DCO confound. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing quantity is the [SII]-deficiency $\Delta[\mathrm{SII}]$, defined as a galaxy's displacement in $\log_{10}([\mathrm{SII}]/\mathrm{H}\alpha)$ from a polynomial ridge-line fitted to the peak locus of SDSS DR12 star-forming galaxies in the plane of $\log_{10}([\mathrm{SII}]/\mathrm{H}\alpha)$ versus $\log_{10}([\mathrm{OIII}]/\mathrm{H}\beta)$. The physical mechanism is that [SII] (ionization potential 10.4 eV) forms mainly in the warm partially-ionized zone at the edge of a classical HII region, so an interstellar medium optically thin to ionizing radiation loses that zone and emits weaker [SII]. The observational machinery includes COS far-UV spectroscopy with a super-dark subtraction procedure inherited from earlier work, Starburst99 stellar-synthesis fits to the OVI and NV wind lines to establish that the far-UV light is stellar, and extinction corrections using two UV reddening laws.
What would settle it
Observe a sample of roughly twenty [SII]-weak star-forming galaxies at $z\sim0.3$ with HST/COS. If the fraction with direct Lyman-continuum detections is not significantly higher than in a control sample of starbursts with normal [SII] emission, then [SII]-deficiency is not an effective signpost. Alternatively, a photoionization-model grid that varies ionization parameter and density while keeping the gas optically thick to LyC would falsify the physical interpretation if it can reproduce the observed [SII] deficit.
Extended reading notes
Core claim
The central claim is that galaxies whose [SII] emission is weak relative to the locus of normal star-forming galaxies in the [SII]/H$\alpha$ versus [OIII]/H$\beta$ plane are likely to be leaking Lyman-continuum photons. The physical reason is that [SII] emission is produced mainly in the warm partially-ionized gas just outside the ionized region; if the interstellar medium is optically thin to ionizing radiation, that zone shrinks or disappears, so [SII] weakens. Two of three such galaxies observed with COS show direct flux below the Lyman edge, with relative escape fractions near unity but absolute escape fractions of only 3% to 4% because of internal dust. These [SII]-weak leakers have stellar masses $\sim10^{10.5}\,M_\odot$, near-solar metallicity, modest ionization parameters, and high star-formation rate per unit area, differing strongly from the Green Peas while sharing compactness and blueshifted Ly$\alpha$ emission with higher-redshift leakers. The paper concludes that [SII]-weakness is an effective and isotropic signpost, and that the five most [SII]-deficient galaxies observed so far are all confirmed leakers.
Load-bearing premise
The method assumes that a deficit of [SII] emission means the interstellar medium is transparent to ionizing photons, rather than being caused by a different nebular condition such as a lower ionization parameter, a harder radiation field, lower density, or lower metallicity.
Editorial extensions
If this is right
- [SII]-weakness can be measured from ground-based optical spectra, so large samples of candidate leakers can be assembled more cheaply than with UV observations.
- The ridge-line passes through the locus of $z\sim2$ to $3$ galaxies in existing surveys, indicating the technique can be exported to reionization-era galaxies and applied with JWST spectroscopy.
- The discovery of leakers that are massive, metal-rich, and dusty shows that LyC escape is not restricted to the extreme low-metallicity starbursts of the Green Pea type.
- All five of the most [SII]-deficient galaxies with LyC observations so far are confirmed leakers, suggesting the signpost may be a common property of LyC escape.
Reading between the lines
- A testable implication the paper leaves open: both leakers contain dominant central objects that produce nearly all the UV light, so the [SII] deficit may trace local optically-thin channels excavated by a compact starburst rather than a globally transparent interstellar medium; the signpost would still work for finding leakers, but its physical meaning would be more local.
- The $\Delta[\mathrm{SII}]$ metric is defined relative to a redshift-zero ridge-line; if the ridge-line evolves with metallicity or excitation, applying a fixed threshold at $z>2$ could bias the selection, and a redshift-calibrated definition would be a natural follow-up.
- A more direct test of the mechanism would combine [SII]-deficiency with independent optical-depth tracers (Ly$\alpha$ escape fraction, Ly$\beta$ residual intensity, or [OII]/[OIII] ratios) in one sample; a tight correlation would convert the signpost into a quantitative estimator of escape fraction.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This manuscript reports a pilot HST/COS program testing [SII]-deficiency as a signpost for Lyman-continuum leakage. The authors define Δ[SII] as the offset from a ridge line fitted to SDSS DR12 star-forming galaxies in the [SII]/Hα versus [OIII]/Hβ plane, select compact z~0.3 starbursts with Δ[SII] < -0.2 dex, and observe three targets with COS G140L after rejecting two objects whose UV spectra are quasar-dominated. Significant flux below the Lyman edge is detected in two galaxies (J0910 and J1432) and only an upper limit is obtained for J1242. Relative and absolute escape fractions are derived from Starburst99 models under several extinction prescriptions, and the leakers are compared with Green Peas, Lyman Break Analogs, and z~3 leakers. The paper concludes that [SII]-weakness is a highly effective way to identify LyC-leaking galaxies.
Significance. The two direct LyC detections are the main value of this paper. They are statistically strong, the super-dark subtraction and blank-sky comparisons are careful, and the escape-fraction analysis considers several SB99 models and two extinction laws. The validation is not circular: the LyC flux is measured directly from COS, while the selection ridge line is derived from SDSS and the intrinsic Lyman break from independent stellar-population models. However, the claim that [SII]-weakness is 'highly effective' is not yet established: the sample is tiny, the only non-detection does not satisfy the stated selection threshold, and both leakers contain dominant central objects while the non-leaker does not. The physical connection between [SII]-weakness and an optically thin ISM is plausible but is not independently tested here. If the conclusion is suitably qualified, the paper remains a useful proof-of-concept for a promising selection technique.
major comments (4)
- [Sec. 5.3] Both LyC detections (J0910 and J1432) contain dominant central objects that produce nearly all of the UV light, while the non-detection J1242 has a significant diffuse UV component. Since all three targets are [SII]-weak, the data cannot separate the hypothesis that LyC escape is caused by [SII]-weakness (density-bounded HII regions) from the hypothesis that escape is driven by DCO feedback and [SII]-weakness is a correlated byproduct. A control sample of compact DCO-bearing starbursts with normal [SII], or of [SII]-weak galaxies without DCOs, is needed before 'highly effective' can be claimed. Please either add such a test or explicitly restrict the conclusion to 'two [SII]-weak galaxies with DCOs show LyC escape.'
- [Sec. 3.1 and Table 2] J1242's measured Δ[SII] is -0.17 dex, below the stated selection threshold of ≥0.2 dex (the text notes this was caused by updating the SDSS sample from DR7 to DR12). Consequently the effective sample of targets satisfying the criterion is two galaxies, both of which are detected, and the paper's 'two out of three' statistic overstates the denominator. With two successes out of two (or three) targets, the binomial 95% confidence interval is extremely wide, so the quantitative conclusion of 'highly effective' is not supported by the pilot alone. Please state the effective sample explicitly, report a confidence interval, and temper the conclusion accordingly.
- [Sec. 4, Tables 4 and 5] The quoted relative escape fractions (93% and 80%) are close to the maximum allowed value, and Table 5 shows that the intrinsic Lyman-break amplitude changes by roughly a factor of two between the 10^7 yr and 10^8 yr constant-SFR models. The text itself notes that adopting the older burst age for J0910 and J1432 would increase fesc,rel by ~0.2 dex and push the values above unity. This means the reported escape fractions are strongly model-dependent. Please quote a systematic error budget that spans the full SB99 grid and/or present fesc as a range, rather than the point values in Table 4.
- [Sec. 1 and Sec. 2] The physical interpretation underlying the technique—that [SII]-weakness traces the absence of the partially ionized outer zone of a density-bounded HII region—is plausible but rests on the Pellegrini et al. (2012) models and is not independently tested here. Low ionization parameter, low electron density, or metallicity effects could also suppress [SII]/Hα, and any of these might accompany the DCOs seen in the two leakers. Please add an independent sanity check (e.g., photoionization modeling of the full SDSS line set) or state explicitly that the empirical correlation, not the physical mechanism, is what the pilot establishes.
minor comments (5)
- [Sec. 3.1] 'Starbust99' is a typo for 'Starburst99'.
- [Sec. 3.3] The choice of best-fit SB99 model is described as 'by eye' using the OVI and NV wind features; a quantitative goodness-of-fit criterion would improve reproducibility.
- [Table 4] The J1242 upper limits are described as 'inferred from a 3σ limit on dark fluxes'; it would be clearer to state that the limit is on the net LyC flux after dark subtraction.
- [Sec. 4, Eq. (5)] The notation F910−/F910+ is used before being defined in Eq. (5); please introduce it in Table 3 or in the text preceding it.
- [Figure 1 and Sec. 2] The grey uncertainty band of the ridge-line is mentioned but the bootstrap procedure is not described in enough detail to reproduce it; a brief description or reference would help.
Circularity Check
No circularity detected: the [SII]-deficiency selection is defined from SDSS data, while the LyC detections come from independent HST/COS measurements.
full rationale
The paper's central validation is self-contained rather than circular. The [SII]-deficiency ridge-line is fitted to the SDSS DR12 star-forming galaxy sample (Section 2, Eq. 1), not to the COS targets, so the selection criterion is not constructed from the measured LyC fluxes. The LyC signal is measured directly in the COS spectra below the Lyman edge (Section 4, Figure 2), with explicit dark-subtraction and sky-contamination checks, and the escape fractions are computed by comparing the observed Lyman break to Starburst99 model spectra selected by fits to stellar wind features (OVI and NV), not by fitting to the LyC flux itself. The non-detection in J1242 provides a genuine upper limit. Prior same-group work appears as methodological or contextual citations: the super-dark recipe (Leitherer et al. 2016), the Lyman Break Analog sample (Alexandroff et al. 2015), and the dominant-central-object interpretation (Heckman et al. 2011; Borthakur et al. 2014). These are tools or interpretive context, and the paper explicitly discusses the DCO alternative in Section 5.3 rather than using it to force the conclusion. The physical assumption that [SII]-weakness traces an optically-thin ISM, supported by Pellegrini et al. (2012) models, is a correctness/interpretation risk rather than a circular step, because the empirical LyC detections do not depend on that assumption. No equation or fitted parameter is re-used as its own prediction, and no load-bearing claim reduces by definition to its inputs.
Assumptions & free parameters
free parameters (3)
- Ridge-line polynomial coefficients =
9 coefficients: -0.487, 0.014, 0.028, -0.785, -3.870, 0.446, 8.696, 0.302, -6.623
- Internal extinction E(B-V)_int per galaxy =
J0910: 0.239; J1432: 0.243; J1242: 0.314 (Calzetti); 0.257, 0.252, 0.325 (Reddy)
- SB99 model burst age, metallicity, and rotation =
10^7 yr, solar Z, rotation for J0910/J1432; 10^8 yr for J1242
assumptions (3)
- domain assumption [SII] emission arises mostly in the warm, partially-ionized zone at the Stromgren edge; optically-thin HII regions lack this zone
- domain assumption Starburst99 stellar population synthesis models (Kroupa IMF, Geneva tracks) provide a reliable intrinsic Lyman break
- domain assumption Calzetti et al. (1999) and Reddy et al. (2015, 2016) reddening laws can be extrapolated to ~910 Å
Cite this review
Pith. "Pith review of A new technique for finding galaxies leaking Lyman-continuum radiation: [SII]-deficiency." pith.science (2026). https://pith.science/paper/ICXJDV3Y
@misc{pith2026190901368,
author = {Pith},
title = {Pith review of: A new technique for finding galaxies leaking Lyman-continuum radiation: [SII]-deficiency},
year = {2026},
howpublished = {\url{https://pith.science/paper/ICXJDV3Y}},
note = {Machine review of arXiv:1909.01368}
}
abstract
The source responsible for the reionization of the Universe is believed to be the population of star-forming galaxies at $z\sim6$ to 12. The biggest uncertainty concerns the fraction of Lyman-continuum photons that actually escape from the galaxies. In recent years, several relatively small samples of "leaky" galaxies have been uncovered, and clues have begun to emerge as to both the indirect signposts of leakiness and of the conditions/processes that enable the escape of ionizing radiation. In this paper we present the results of a pilot program aimed to test a new technique for finding leaky galaxies---using the weakness of the [SII] nebular emission-lines relative to typical star-forming galaxies as evidence that the interstellar medium is optically-thin to the Lyman continuum. We use the Cosmic Origins Spectrograph on the Hubble Space Telescope to detect significant emerging flux below the Lyman edge in two out of three [SII]-weak star-forming galaxies at $z\sim0.3$. We show that these galaxies differ markedly in their properties from the class of leaky "Green-Pea" galaxies at similar redshifts: our sample galaxies are more massive, more metal-rich, and less extreme in terms of their stellar population and the ionization state of the interstellar medium. Like the Green Peas, they have exceptionally high star-formation rates per unit area. They also share some properties with the known leaky galaxies at $z\sim3$, but are significantly dustier. Our results validate a new way to identify local laboratories for exploring the processes that made it possible for galaxies to reionize the Universe.
Figures
Figures from the paper (4 more)
Reference graph
Works this paper leans on
-
[1]
thebibliography [1] 20pt to REFERENCES 6pt =0pt 10pt plus 3pt =0pt =0pt =1pt plus 1pt =0pt =0pt -12pt =13pt plus 1pt =20pt =13pt plus 1pt \@M =10000 =-1.0em =0pt =0pt 0pt =0pt =1.0em @enumiv\@empty 10000 10000 `\.\@m \@noitemerr \@latex@warning Empty `thebibliography' environment \@ifnextchar \@reference \@latexerr Missing key on reference command Each re...
arXiv 2017
-
[2]
Alexandroff, R. M., Heckman, T. M., Borthakur, S., Overzier, R., & Leitherer, C. 2015, , 810, 104
work page 2015
-
[3]
P., Tollerud , E
Astropy Collaboration , Robitaille , T. P., Tollerud , E. J., et al. 2013, , 558, A33
2013
-
[4]
M., Sip o cz , B
Astropy Collaboration , Price-Whelan , A. M., Sip o cz , B. M., et al. 2018, , 156, 123
2018
-
[5]
Bennett , C. L., Larson , D., Weiland , J. L., et al. 2013, , 208, 20
work page 2013
-
[6]
M., Leitherer, C., & Overzier, R
Borthakur, S., Heckman, T. M., Leitherer, C., & Overzier, R. A. 2014, Science, 346, 216
work page 2014
-
[7]
J., Smit, R., Labb \' e , I., et al
Bouwens, R. J., Smit, R., Labb \' e , I., et al. 2016, , 831, 176
work page 2016
-
[8]
2011, in Secular Evolution of Galaxies, 419--458
Calzetti, D. 2011, in Secular Evolution of Galaxies, 419--458
work page 2011
Show all 48 references
-
[9]
C., et al
Calzetti, D., Armus, L., Bohlin, R. C., et al. 1999, , 533, 682
1999
-
[10]
2018, , 616, 30
Chisholm, J., Gazagnes, S., Schaerer, D., et al. 2018, , 616, 30
2018
-
[11]
2006, , 44, 415
Fan, X., Carilli, C., & Keating, B. 2006, , 44, 415
2006
-
[12]
A., Becker , R
Fan , X., Strauss , M. A., Becker , R. H., et al. 2006, , 132, 117
2006
-
[13]
C., Froning , C
Green , J. C., Froning , C. S., Osterman , S., et al. 2012, , 744, 60
2012
-
[14]
E., & Peterson , B
Gunn , J. E., & Peterson , B. A. 1965, , 142, 1633
1965
-
[15]
M., Alexandroff, R
Heckman, T. M., Alexandroff, R. M., Borthakur, S., Overzier, R., & Leitherer, C. 2015, , 809, 147
2015
-
[16]
M., Borthakur, S., Overzier, R., et al
Heckman, T. M., Borthakur, S., Overzier, R., et al. 2011, , 730, doi:10.1088/0004-637X/730/1/5
2011 doi
-
[17]
Hunter , J. D. 2007, Computing in Science Engineering, 9, 90
2007
-
[18]
I., Orlitov \' a , I., Schaerer, D., et al
Izotov, Y. I., Orlitov \' a , I., Schaerer, D., et al. 2016 a , , 529, 178
2016
-
[19]
I., Schaerer, D., Thuan, T
Izotov, Y. I., Schaerer, D., Thuan, T. X., et al. 2016 b , , 461, 3683
2016
-
[20]
I., Schaerer, D., Worseck, G., et al
Izotov, Y. I., Schaerer, D., Worseck, G., et al. 2018 a , , 474, 4514
2018
-
[21]
I., Worseck, G., Schaerer, D., et al
Izotov, Y. I., Worseck, G., Schaerer, D., et al. 2018 b , , 478, 4851
2018
-
[22]
E., Dowd , T., Oey , M
Jaskot , A. E., Dowd , T., Oey , M. S., Scarlata , C., & McKinney , J. 2019, arXiv e-prints, arXiv:1908.09763
2019 arXiv
-
[23]
2001, SciPy : Open source scientific tools for Python , , , [Online; accessed <today>]
Jones, E., Oliphant, T., Peterson, P., et al. 2001, SciPy : Open source scientific tools for Python , , , [Online; accessed <today>]
2001
-
[24]
C., & Evans, N
Kennicutt, R. C., & Evans, N. J. 2012, , 14, 3
2012
-
[25]
2001, , 322, 231
Kroupa, P. 2001, , 322, 231
2001
-
[26]
C., & Oey, M
Leitherer, C., Hernandez, S., Lee, J. C., & Oey, M. S. 2016, , 823, 64
2016
-
[27]
A., Bresolin, F., et al
Leitherer, C., Ortiz Ot \' a lvaro , P. A., Bresolin, F., et al. 2010, , 189, 309
2010
-
[28]
D., et al
Leitherer, C., Schaerer, D., Goldader, J. D., et al. 1999, , 123, 3
1999
-
[29]
2017, Astrophysics A & A, 614, 11
Marchi, F., Pentericci, L., Guaita, L., et al. 2017, Astrophysics A & A, 614, 11
2017
-
[30]
1990, , 28, 37
Mathis, J. 1990, , 28, 37
1990
-
[31]
2016, , doi:10.1146/annurev-astro-082214-122355
Mcquinn, M. 2016, , doi:10.1146/annurev-astro-082214-122355
2016 doi
-
[32]
A., Heckman, T
Overzier, R. A., Heckman, T. M., Tremonti, C., et al. 2009, , 706, 203
2009
-
[33]
W., Oey, M
Pellegrini, E. W., Oey, M. S., Winkler, P. F., et al. 2012, , 755, 40
2012
-
[34]
Pettini, M., & Pagel, B. E. J. 2004, , 348, L59
2004
-
[35]
2018 a , arXiv e-prints, arXiv:1807.06205
Planck Collaboration , Akrami , Y., Arroja , F., et al. 2018 a , arXiv e-prints, arXiv:1807.06205
2018 arXiv
-
[36]
2018 b , arXiv e-prints, arXiv:1807.06209
Planck Collaboration , Aghanim, N., Akrami, Y., et al. 2018 b , arXiv e-prints, arXiv:1807.06209
2018 arXiv
-
[37]
A., Pettini, M., Steidel, C
Reddy, N. A., Pettini, M., Steidel, C. C., et al. 2012, , 754, 25
2012
-
[38]
A., Steidel, C
Reddy, N. A., Steidel, C. C., Pettini, M., & Bogosavljevic, M. 2016, , arXiv:1606.00434
2016 arXiv
-
[39]
A., Kriek, M., Shapley, A
Reddy, N. A., Kriek, M., Shapley, A. E., et al. 2015, , 806, 259
2015
-
[40]
E., Ellis, R
Robertson, B. E., Ellis, R. S., Furlanetto, S. R., & Dunlop, J. S. 2015, , 802, L19
2015
-
[41]
2019, vrodgom/statmorph: v0.3.3, , , doi:10.5281/zenodo.2535876
Rodriguez-Gomez, V. 2019, vrodgom/statmorph: v0.3.3, , , doi:10.5281/zenodo.2535876
2019 doi
-
[42]
C., Bogosavljevi \' c , M., Shapley, A
Steidel, C. C., Bogosavljevi \' c , M., Shapley, A. E., et al. 2018, , 869, 123
2018
-
[43]
L., Steidel , C
Strom , A. L., Steidel , C. C., Rudie , G. C., Trainor , R. F., & Pettini , M. 2018, , 868, 117
2018
-
[44]
2013, , 431, 1383
Thomas , D., Steele , O., Maraston , C., et al. 2013, , 431, 1383
2013
-
[45]
C., & Varoquaux , G
van der Walt , S., Colbert , S. C., & Varoquaux , G. 2011, Computing in Science Engineering, 13, 22
2011
-
[46]
2018, , 476, L15
Vanzella , E., Nonino , M., Cupani , G., et al. 2018, , 476, L15
2018
-
[47]
B., Calura , F., et al
Vanzella , E., Caminha , G. B., Calura , F., et al. 2019, arXiv e-prints, arXiv:1904.07941
2019 arXiv
-
[48]
L., Eisenhardt , P
Wright , E. L., Eisenhardt , P. R. M., Mainzer , A. K., et al. 2010, , 140, 1868
2010
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.