arxiv: 2605.01380 · v1 · submitted 2026-05-02 · 🌌 astro-ph.CO · astro-ph.GA

Recognition: unknown

A structural degeneracy explains reionization tensions and limits dark matter constraints

Zihan Wang , Huanyuan Shan

Authors on Pith no claims yet

Pith reviewed 2026-05-09 18:05 UTC · model grok-4.3

classification 🌌 astro-ph.CO astro-ph.GA

keywords reionizationionizing photonsescape fractionstar formation efficiencydegeneracydark matterThomson optical depthUV luminosity function

0 comments

The pith

The product of escape fraction and star formation efficiency creates an exact degeneracy in all reionization observations.

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

The paper demonstrates that the rate of ionizing photons is always proportional to the product of the escape fraction f_esc and the peak star formation efficiency f_*,0. Because of this, every major probe of reionization history from the optical depth to the neutral hydrogen fraction and the UV luminosity function responds only to that product. This makes it impossible to separately determine the two parameters from current data, explaining why different analyses report conflicting values. The authors confirm the effect persists even when galaxy structure changes substantially in dark matter simulations, as long as the effective emissivity stays the same.

Core claim

The ionizing photon rate is proportional to f_esc times f_*,0, rendering all reionization-history probes sensitive only to their product and creating an intrinsically non-invertible mapping between model parameters and observations. This structural degeneracy holds across Thomson optical depth, neutral hydrogen fraction, UV luminosity function, and quasar proximity zones.

What carries the argument

The algebraic degeneracy in the ionizing photon rate proportional to f_esc times f_*,0 that collapses separate constraints on escape fraction and star formation efficiency into a single effective emissivity.

If this is right

Published constraints on f_esc and f_*,0 from reionization are mutually inconsistent because they are degenerate.
Reionization observations cannot tightly constrain self-interacting dark matter or other models that alter halo structure.
Only observables that probe the spatial arrangement of ionized regions can separate the two parameters.
JWST observations of bright high-redshift galaxies are consistent with the degeneracy without requiring revisions to galaxy formation models.

Where Pith is reading between the lines

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

21 cm tomography or measurements of the ionized bubble size distribution could provide the missing topological information to break the degeneracy.
Similar product degeneracies may appear in other high-redshift observables that combine galaxy efficiency and photon escape.
Modeling efforts should prioritize the effective ionizing emissivity over independent values of f_esc and f_*,0.

Load-bearing premise

That the ionizing photon rate depends strictly on the product of escape fraction and star formation efficiency with no additional probe-specific compensating factors.

What would settle it

A simulation or observation showing that different pairs of f_esc and f_*,0 that give the same product produce measurably different reionization histories when structural effects are included.

read the original abstract

Over the past decade, reionization studies have yielded persistent factor-of-two-to-five disagreements in the inferred ionizing escape fraction $f_{\mathrm{esc}}$ and peak star formation efficiency $f_{*,0}$, compounded by JWST's discovery of unexpectedly bright $z>10$ galaxies. We show that this discrepancy arises from an algebraically exact structural degeneracy: the ionizing photon rate $\dot{n}_{\mathrm{ion}} \propto f_{\mathrm{esc}} \times f_{*,0}$ renders all reionization-history probes, including Thomson optical depth, neutral hydrogen fraction, UV luminosity function, and quasar proximity zones, sensitive only to their product, leading to an intrinsically non-invertible mapping between model parameters and observations. We demonstrate the robustness of this degeneracy using a large suite of N-body simulations of self-interacting dark matter haloes spanning $10^9$-$10^{11} M_\odot$. Despite substantial changes to galaxy-scale structure, observables remain indistinguishable once the effective ionizing emissivity is matched, severely limiting reionization-based dark matter probes. We identify that only observables sensitive to the spatial topology of ionized regions can break this degeneracy. Our results provide a unified explanation for the scatter among published constraints and establish a framework for interpreting reionization observations and their implications for early galaxy formation and dark matter.

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.

Desk Editor's Note private letter to a colleague

The paper correctly flags an algebraic degeneracy between f_esc and f_{*,0} for ionization history but incorrectly includes the UV luminosity function as one of the affected probes.

read the letter

The paper's main contribution is showing that the ionizing photon rate scales with the product of escape fraction and peak star formation efficiency, which creates a non-invertible mapping for several reionization observables and accounts for the factor-of-two-to-five scatter in published constraints. The N-body simulations of self-interacting dark matter halos then demonstrate that structural changes do not alter the observables once effective emissivity is matched. That simulation check is a solid, independent test and adds real weight to the claim that reionization data have limited reach for dark matter constraints. The algebraic part follows directly from the definition of the rate and is exact on its own terms. The unified framing of the tensions is also helpful for people who work with these datasets. The clear soft spot is the inclusion of the UV luminosity function in the list of degenerate probes. Observed UV magnitudes at 1500 Å trace stellar continuum emission, which scales with f_{*,0} and the halo mass function under standard assumptions; f_esc only controls the fraction of Lyman-continuum photons that escape and does not change the UV continuum. Thomson optical depth, neutral fraction, and proximity zones integrate the escaped ionizing photons and are therefore sensitive only to the product, but the UV LF is not. This differential dependence means UV LF data can constrain f_{*,0} separately and breaks the claimed non-invertibility for that observable. The paper needs to drop or qualify the UV LF from the degeneracy statement. The simulations themselves look fine for the structural test they perform, but they do not address this modeling distinction. This work is for modelers who fit reionization histories and early galaxy populations. Readers focused on parameter degeneracies and the reach of reionization data for dark matter will find the algebraic and simulation pieces useful. It deserves a serious referee because the degeneracy for the ionization-sensitive probes is worth discussing and the simulation evidence is concrete, even though the UV LF claim requires revision.

Referee Report

2 major / 1 minor

Summary. The manuscript claims that an algebraically exact structural degeneracy, arising because the ionizing photon rate ṅ_ion ∝ f_esc × f_{*,0}, renders every reionization-history probe—including Thomson optical depth, neutral hydrogen fraction, UV luminosity function, and quasar proximity zones—sensitive exclusively to the product of these parameters. This produces a non-invertible mapping from model parameters to observations, explaining the persistent factor-of-two-to-five tensions in inferred f_esc and f_{*,0} values (including those highlighted by JWST bright z>10 galaxies). A large suite of N-body simulations of self-interacting dark matter halos (10^9–10^11 M_⊙) is used to show that substantial changes in galaxy-scale structure leave observables indistinguishable once the effective ionizing emissivity is fixed, thereby limiting the power of reionization data to constrain dark matter models. Only observables sensitive to the spatial topology of ionized regions can break the degeneracy.

Significance. If the degeneracy holds for the ionizing observables, the work supplies a unified algebraic explanation for the scatter among published reionization constraints and demonstrates that reionization-based limits on dark matter are weaker than commonly assumed. The N-body simulation suite is a concrete strength: it shows that once emissivity is matched, observables remain indistinguishable despite large structural variations, providing direct support for the robustness claim. The over-inclusion of the UV luminosity function in the list of degenerate probes, however, reduces the overall significance because that observable supplies an independent constraint on f_{*,0}.

major comments (2)

[Abstract] Abstract: the statement that 'all reionization-history probes, including ... UV luminosity function ... [are] sensitive only to their product' is incorrect. Observed UV magnitudes at ~1500 Å trace stellar continuum emission, which scales directly with f_{*,0} (and the halo mass function) under standard IMF and dust assumptions; f_esc affects only the Lyman-continuum leakage and leaves the UV continuum unchanged. In contrast, Thomson optical depth, x_HI, and proximity zones integrate the escaped rate ∝ f_esc × f_{*,0}. Consequently, UV LF data can constrain f_{*,0} separately and break the claimed universal non-invertibility.
[Abstract and simulation section] Abstract and simulation section: the N-body SIDM suite demonstrates that halo structural variations do not affect observables once emissivity is fixed, but this test assumes the emissivity (i.e., the product) has already been set and does not address the differential dependence of the UV LF on f_{*,0} alone. Because the central claim that the degeneracy 'limits dark matter constraints' rests on every listed probe being degenerate, the UV LF counter-example is load-bearing and requires correction of the scope of the argument.

minor comments (1)

[Abstract] The phrase 'factor-of-two-to-five disagreements' would benefit from one or two explicit citations to the specific literature values being referenced.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful review and for identifying the overstatement regarding the UV luminosity function. We agree that this observable is not subject to the same degeneracy and will revise the abstract and manuscript to correct the scope of the claim. The core algebraic argument for reionization-history probes that depend on the escaped ionizing emissivity remains intact, as does the demonstration from the SIDM simulation suite that structural variations do not affect those observables once the product is fixed.

read point-by-point responses

Referee: [Abstract] Abstract: the statement that 'all reionization-history probes, including ... UV luminosity function ... [are] sensitive only to their product' is incorrect. Observed UV magnitudes at ~1500 Å trace stellar continuum emission, which scales directly with f_{*,0} (and the halo mass function) under standard IMF and dust assumptions; f_esc affects only the Lyman-continuum leakage and leaves the UV continuum unchanged. In contrast, Thomson optical depth, x_HI, and proximity zones integrate the escaped rate ∝ f_esc × f_{*,0}. Consequently, UV LF data can constrain f_{*,0} separately and break the claimed universal non-invertibility.

Authors: We agree with the referee. The UV continuum luminosity at 1500 Å is produced by stellar populations and is independent of the escape fraction f_esc under standard assumptions. Therefore the UV luminosity function provides a direct constraint on f_{*,0} (modulo the halo mass function, IMF, and dust) and is not degenerate with the product f_esc × f_{*,0}. The algebraic degeneracy applies strictly to observables that integrate the escaped ionizing photon rate. We will remove the UV luminosity function from the list of degenerate probes in the abstract and add a clarifying sentence in the introduction and discussion sections to delineate which observables are affected. revision: yes
Referee: [Abstract and simulation section] Abstract and simulation section: the N-body SIDM suite demonstrates that halo structural variations do not affect observables once emissivity is fixed, but this test assumes the emissivity (i.e., the product) has already been set and does not address the differential dependence of the UV LF on f_{*,0} alone. Because the central claim that the degeneracy 'limits dark matter constraints' rests on every listed probe being degenerate, the UV LF counter-example is load-bearing and requires correction of the scope of the argument.

Authors: The N-body simulations were designed to test whether galaxy-scale structural changes (in SIDM halos of 10^9–10^11 M_⊙) alter reionization observables when the effective ionizing emissivity is held fixed. They show that once the product f_esc × f_{*,0} is matched, quantities such as the Thomson optical depth, volume-averaged neutral fraction, and proximity-zone sizes remain indistinguishable. Because we are revising the manuscript to exclude the UV luminosity function from the degenerate set, the simulation results continue to support the claim that reionization-history data alone cannot tightly constrain dark-matter-induced structural variations. The UV LF will now be presented as an independent anchor on f_{*,0} that can be combined with the degenerate probes to tighten constraints. revision: partial

Circularity Check

0 steps flagged

Algebraic degeneracy follows from standard emissivity definition; derivation self-contained with independent simulation test

full rationale

The paper's core step states that the ionizing photon rate proportionality renders all listed probes (including UV LF) sensitive only to the product f_esc × f_*,0. This follows directly from the conventional definition of emissivity rather than any fitted parameter, self-referential assumption, or imported uniqueness theorem. The N-body SIDM suite supplies an independent structural check by holding effective emissivity fixed and verifying observable indistinguishability. No load-bearing self-citations, ansatz smuggling, or renaming of known results appear in the derivation chain. The debatable extension to UV LF (which traces stellar continuum independently) is a modeling scope issue, not a circular reduction of the claimed result to its inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the standard domain assumption that the ionizing photon rate is strictly proportional to the product of escape fraction and star formation efficiency. No free parameters are introduced for the degeneracy itself. No new physical entities are postulated. The simulations rely on existing N-body techniques for self-interacting dark matter.

axioms (1)

domain assumption The ionizing photon production rate is proportional to the product of the escape fraction f_esc and the peak star formation efficiency f_*,0
Invoked directly in the abstract as the basis for the structural degeneracy affecting all listed reionization probes.

pith-pipeline@v0.9.0 · 5533 in / 1407 out tokens · 85609 ms · 2026-05-09T18:05:24.644066+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

45 extracted references · 5 canonical work pages · 1 internal anchor

[1]

Aghanim, N.et al.Planck 2018 results. VI. Cosmological parameters.Astron. Astrophys.641, A6 (2020). [Erratum: Astron.Astrophys. 652, C4 (2021)]

2018
[2]

Umeda, H.et al.JWST Measurements of Neutral Hydrogen Fractions and Ionized Bubble Sizes at z = 7–12 Obtained with LyαDamping Wing Absorptions in 27 Bright Continuum Galaxies.apj971, 124 (2024)

2024
[3]

A.et al.Constraints on the z∼6−13 intergalactic medium from JWST spectroscopy of Lyman-alpha damping wings in galaxies.Astron

Mason, C. A.et al.Constraints on the z∼6−13 intergalactic medium from JWST spectroscopy of Lyman-alpha damping wings in galaxies.Astron. Astrophys.705, A114 (2026)

2026
[4]

& Mesinger, A

Greig, B. & Mesinger, A. Simultaneously constraining the astrophysics of reion- ization and the epoch of heating with 21CMMC.Mon. Not. Roy. Astron. Soc. 472, 2651–2669 (2017)

2017
[5]

T.et al.JWST PRIMER: a new multifield determination of the evolving galaxy UV luminosity function at redshifts z≃9 –15.mnras533, 7 3222–3237 (2024)

Donnan, C. T.et al.JWST PRIMER: a new multifield determination of the evolving galaxy UV luminosity function at redshifts z≃9 –15.mnras533, 7 3222–3237 (2024)

2024
[6]

Harikane, Y.et al.JWST, ALMA, and Keck Spectroscopic Constraints on the UV Luminosity Functions at z∼7–14: Clumpiness and Compactness of the Brightest Galaxies in the Early Universe.apj980, 138 (2025)

2025
[7]

J.et al.A Spectroscopic Analysis of the Ionizing Photon Production Efficiency in JADES and CEERS: Implications for the Ionizing Photon Budget

Pahl, A. J.et al.A Spectroscopic Analysis of the Ionizing Photon Production Efficiency in JADES and CEERS: Implications for the Ionizing Photon Budget. Astrophys. J.981, 134 (2025)

2025
[8]

mnras537, 3245–3264 (2025)

Begley, R.et al.The evolution of [O III] + Hβequivalent width from z≃3–8: implications for the production and escape of ionizing photons during reionization. mnras537, 3245–3264 (2025)

2025
[9]

Mascia, S.et al.Insights into the reionization epoch from cosmic-noon-C IV emitters in the VANDELS survey.aap674, A221 (2023)

2023
[10]

Unveiling the ISM properties of low-z Lyman-continuum emitters.aap663, A59 (2022)

Saldana-Lopez, A.et al.The Low-Redshift Lyman Continuum Survey. Unveiling the ISM properties of low-z Lyman-continuum emitters.aap663, A59 (2022)

2022
[11]

, keywords =

Naidu, R. P.et al.Two remarkably luminous galaxy candidates at z≈10–12 revealed by jwst.The Astrophysical Journal Letters940, L14 (2022). URL http://dx.doi.org/10.3847/2041-8213/ac9b22

work page doi:10.3847/2041-8213/ac9b22 2022
[12]

L.et al.Conditions for Reionizing the Universe with A Low Galaxy Ionizing Photon Escape Fraction.Astrophys

Finkelstein, S. L.et al.Conditions for Reionizing the Universe with A Low Galaxy Ionizing Photon Escape Fraction.Astrophys. J.879, 36 (2019)

2019
[13]

L.et al.A Long Time Ago in a Galaxy Far, Far Away: A Candidate z∼12 Galaxy in Early JWST CEERS Imaging.apjl940, L55 (2022)

Finkelstein, S. L.et al.A Long Time Ago in a Galaxy Far, Far Away: A Candidate z∼12 Galaxy in Early JWST CEERS Imaging.apjl940, L55 (2022)

2022
[14]

E., Ellis, R

Robertson, B. E., Ellis, R. S., Furlanetto, S. R. & Dunlop, J. S. Cosmic Reioniza- tion and Early Star-forming Galaxies: A Joint Analysis of New Constraints from Planck and the Hubble Space Telescope.apjl802, L19 (2015)

2015
[15]

Wyithe, J. S. B., Hopkins, A. M., Kistler, M. D., Y¨ uksel, H. & Beacom, J. F. Determining the escape fraction of ionizing photons during reionization with the GRB-derived star formation rate.mnras401, 2561–2571 (2010)

2010
[16]

& Shan, H

Wang, Z. & Shan, H. Breaking the UV Luminosity Function Degeneracy:Self- Interacting Dark Matter Constraints from Reionization Topology (2026)

2026
[17]

Reionization Topology as a Probe of Self-Interacting Dark Matter (2026)

Wang, Z. Reionization Topology as a Probe of Self-Interacting Dark Matter (2026)

2026
[18]

Koopmans, L. V. E.et al.The Cosmic Dawn and Epoch of Reionization with the Square Kilometre Array.PoSAASKA14, 001 (2015). 8

2015
[19]

URL http://dx.doi.org/10

Mellema, G.et al.Reionization and the cosmic dawn with the square kilometre array.Experimental Astronomy36, 235–318 (2013). URL http://dx.doi.org/10. 1007/s10686-013-9334-5

2013
[20]

Bullock, J. S. & Boylan-Kolchin, M. Small-Scale Challenges to the ΛCDM Paradigm.Ann. Rev. Astron. Astrophys.55, 343–387 (2017)

2017
[21]

Spergel, D. N. & Steinhardt, P. J. Observational evidence for selfinteracting cold dark matter.Phys. Rev. Lett.84, 3760–3763 (2000)

2000
[22]

& Yu, H.-B

Tulin, S. & Yu, H.-B. Dark Matter Self-interactions and Small Scale Structure. Phys. Rept.730, 1–57 (2018)

2018
[23]

& Yu, H.-B

Kaplinghat, M., Tulin, S. & Yu, H.-B. Dark Matter Halos as Particle Colliders: Unified Solution to Small-Scale Structure Puzzles from Dwarfs to Clusters.Phys. Rev. Lett.116, 041302 (2016)

2016
[24]

mnras515, 2386–2414 (2022)

Rosdahl, J.et al.LyC escape from SPHINX galaxies in the Epoch of Reionization. mnras515, 2386–2414 (2022)

2022
[25]

Yeh, J. Y. C.et al.The thesan project: ionizing escape fractions of reionization-era galaxies.Mon. Not. Roy. Astron. Soc.520, 2757–2780 (2023)

2023
[26]

Ma, X.et al.No missing photons for reionization: moderate ionizing photon escape fractions from the FIRE-2 simulations.Mon. Not. Roy. Astron. Soc.498, 2001–2017 (2020)

2001
[27]

Hopkins, P. F. A new class of accurate, mesh-free hydrodynamic simulation methods.Mon. Not. Roy. Astron. Soc.450, 53–110 (2015)

2015
[28]

Wang, Y.-Y., Yang, D., Lu, K., Tsai, Y.-L. S. & Fan, Y.-Z. Dark-matter-deficient Galaxies from Collisions: A New Probe of Bursty Feedback and Dark Matter Physics.Astrophys. J. Lett.997, L43 (2026)

2026
[29]

Bose, S.et al.The COpernicus COmplexio: Statistical Properties of Warm Dark Matter Haloes.Mon. Not. Roy. Astron. Soc.455, 318–333 (2016)

2016
[30]

& Huang, K.-W

Schive, H.-Y., Chiueh, T., Broadhurst, T. & Huang, K.-W. Contrasting Galaxy Formation from Quantum Wave Dark Matter,ψDM, with ΛCDM, using Planck and Hubble Data.Astrophys. J.818, 89 (2016)

2016
[31]

Carucci, I. P. & Corasaniti, P.-S. Cosmic Reionization History and Dark Matter Scenarios.Phys. Rev. D99, 023518 (2019)

2019
[32]

A.et al.Inferences on the timeline of reionization at z∼8 from the KMOS Lens-Amplified Spectroscopic Survey.Mon

Mason, C. A.et al.Inferences on the timeline of reionization at z∼8 from the KMOS Lens-Amplified Spectroscopic Survey.Mon. Not. Roy. Astron. Soc.485, 3947–3969 (2019). 9

2019
[33]

& Hayes, M

Bruton, S., Lin, Y.-H., Scarlata, C. & Hayes, M. J. The Universe is at Most 88% Neutral at z = 10.6.apjl949, L40 (2023)

2023
[34]

BEACON: JWST NIRCam Pure-parallel Imaging Survey. III. Constraints on the UV LF and the Clustering of z~7-14 Galaxies

Kreilgaard, K. C.et al.Beacon: Jwst nircam pure-parallel imaging survey. iii. constraints on the uv lf and the clustering of z 7-14 galaxies (2026). URL https: //arxiv.org/abs/2604.17963. arXiv:2604.17963

work page internal anchor Pith review Pith/arXiv arXiv 2026
[35]

Rocha, M.et al.Cosmological Simulations with Self-Interacting Dark Matter I: Constant Density Cores and Substructure.Mon. Not. Roy. Astron. Soc.430, 81–104 (2013)

2013
[36]

H.et al.Sidm on fire: hydrodynamical self-interacting dark matter sim- ulations of low-mass dwarf galaxies.Monthly Notices of the Royal Astronomical Society472, 2945–2954 (2017)

Robles, V. H.et al.Sidm on fire: hydrodynamical self-interacting dark matter sim- ulations of low-mass dwarf galaxies.Monthly Notices of the Royal Astronomical Society472, 2945–2954 (2017). URL http://dx.doi.org/10.1093/mnras/stx2253

work page doi:10.1093/mnras/stx2253 2017
[37]

Meskhidze, H.et al.Comparing implementations of self-interacting dark matter in the GIZMO and AREPO codes.mnras513, 2600–2608 (2022)

2022
[38]

Dutton, A. A. & Macci` o, A. V. Cold dark matter haloes in the Planck era: evolution of structural parameters for Einasto and NFW profiles.Mon. Not. Roy. Astron. Soc.441, 3359–3374 (2014)

2014
[39]

F., Frenk, C

Navarro, J. F., Frenk, C. S. & White, S. D. M. A Universal density profile from hierarchical clustering.Astrophys. J.490, 493–508 (1997)

1997
[40]

galpy: A python library for galactic dynamics.The Astrophysical Journal Supplement Series216, 29 (2015)

Bovy, J. galpy: A python library for galactic dynamics.The Astrophysical Journal Supplement Series216, 29 (2015). URL http://dx.doi.org/10.1088/0067-0049/ 216/2/29

work page doi:10.1088/0067-0049/ 2015
[41]

Widrow, L. M. Distribution Functions for Cuspy Dark Matter Density Profiles. ApJS131, 39–46 (2000)

2000
[42]

A.et al.Self-interacting dark matter in the center of a Local Group dwarf galaxy and its satellites.arXiv e-printsarXiv:2510.05258 (2025)

Gutcke, T. A.et al.Self-interacting dark matter in the center of a Local Group dwarf galaxy and its satellites.arXiv e-printsarXiv:2510.05258 (2025)

work page arXiv 2025
[43]

M., Harness, A., Trenti, M

Shull, J. M., Harness, A., Trenti, M. & Smith, B. D. Critical Star Formation Rates for Reionization: Full Reionization Occurs at Redshift z≈7.apj747, 100 (2012)

2012
[44]

& Eke, V

Robertson, A., Massey, R. & Eke, V. Cosmic particle colliders: simulations of self-interacting dark matter with anisotropic scattering.Mon. Not. Roy. Astron. Soc.467, 4719–4730 (2017)

2017
[45]

Robertson, B. E. Galaxy Formation and Reionization: Key Unknowns and Expected Breakthroughs by the James Webb Space Telescope (2021). 10 Methods GIZMO simulation grid We run 230 dark-matter-onlyN-body simulations using GIZMO[27] with its SIDM module[35–37]. The SIDM implementation uses a Monte Carlo approach in which par- ticle pairs within a kernel smoot...

2021