Recognition: unknown
Double-Peaked Lyα Emission during Reionization Requires Nearby Voids and a Favorable Local Ionizing Background
Pith reviewed 2026-05-07 12:47 UTC · model grok-4.3
The pith
Double-peaked Lyα emission during reionization occurs only when sightlines pass through nearby underdense voids under a strong local ionizing background.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Double-peaked Lyα profiles are reproduced when sightlines from galaxies intersect underdense voids located a few comoving megaparsecs away. These voids become highly transmissive because neutral hydrogen density scales with the square of gas density in ionization equilibrium, dropping opacity sharply below 30 percent of mean density and allowing blueward photons to redshift out of resonance. No such cases appear at 60 percent global ionization fraction; at 80 percent the probability is approximately 3×10^{-3} and increases by a factor of about 10^4 for a tenfold rise in local ionizing background intensity.
What carries the argument
Lyα transmission calculation along sightlines from galaxies in the Cosmic Dawn III simulation, where voids with gas density ≲30 percent of cosmic mean permit resonant photons to pass due to quadratic suppression of neutral hydrogen.
If this is right
- The fraction of double-peaked Lyα emitters serves as a sensitive probe of the local ionizing background intensity in the late stages of reionization.
- Double-peaked profiles should be absent or extremely rare when the universe is only 60 percent ionized.
- The same underdense voids that enable double-peaked emission can also produce transmissive spikes in the Lyα forest.
- These configurations arise naturally from cosmic density fluctuations even in the vicinity of galaxies.
Where Pith is reading between the lines
- Surveys measuring the abundance of double-peaked emitters at different redshifts could map spatial variations in the ionizing background.
- If double-peaked emission proves more common than the simulated rate, it may indicate that current simulations underestimate the abundance or transmissivity of small-scale voids.
- Targeted observations around known double-peaked galaxies could directly search for the predicted nearby voids using other tracers such as the Lyα forest or 21 cm emission.
Load-bearing premise
The simulation accurately captures small-scale gas density fluctuations and ionization equilibrium around galaxies at redshifts above six, and the sightline transmission calculation identifies void effects without major numerical artifacts.
What would settle it
Finding double-peaked Lyα emission in a substantial fraction of galaxies at only 60 percent global ionization fraction, or observing no strong increase in occurrence with higher local ionizing background intensity, would contradict the mechanism.
Figures
read the original abstract
Several Lyman-alpha (Ly$\alpha$) emitters deep into the reionization era exhibit double-peaked Ly$\alpha$ emission profiles, raising the question of how the intergalactic medium can transmit photons blueward of the Ly$\alpha$ resonance at such high redshifts. To investigate this, we compute Ly$\alpha$ transmission along sightlines originating from galaxies in the Cosmic Dawn III simulation and identify cases that closely reproduce the observed double-peaked emission. In these cases, the sightlines intersect highly underdense voids located a few comoving megaparsecs from the source galaxy. These voids allow photons emitted blueward of Ly$\alpha$ to redshift through resonance without scattering while traversing them. The low opacity arises because the neutral hydrogen density scales with the square of the underlying gas density under ionization equilibrium, making sufficiently underdense regions with $\lesssim30~\%$ of cosmic mean density highly transmissive. Such voids naturally occur in the fluctuating cosmic density field, even in the vicinity of galaxies, and can also be associated with transmissive spikes in the Ly$\alpha$ forest. We find that the global probability of observing double-peaked emission is $\sim3\times10^{-3}$ during reionization at an 80\% global ionization fraction, while no cases are found at 60\% ionization. We also find that this probability depends sensitively on the local ionizing background intensity, increasing by $\sim10^{4}$ for a tenfold increase in intensity. These results suggest that the fraction of double-peaked Ly$\alpha$ emission in high-$z$ galaxies can serve as a sensitive probe of the ionizing background during the late stages of cosmic reionization.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript uses the Cosmic Dawn III simulation to compute Lyα transmission along sightlines from galaxies and identifies that double-peaked emission profiles during reionization arise when sightlines intersect nearby underdense voids (≲30% mean density a few cMpc away). These voids permit blueward photons to redshift through resonance with low opacity because n_HI scales as n_gas² under ionization equilibrium. The paper reports a global probability of ∼3×10^{-3} for observing such profiles at 80% global ionization fraction (none at 60%), with the probability increasing by ∼10^4 for a tenfold rise in local ionizing background intensity. It concludes that the fraction of double-peaked Lyα can probe the ionizing background in late reionization.
Significance. If the central results are robust, the work supplies a concrete physical mechanism linking observed double-peaked Lyα lines to the fluctuating density field and positions the incidence of such profiles as a sensitive diagnostic of local ionizing background intensity. The forward ray-tracing through an independent simulation (no fitted parameters) and the explicit link to transmissive Lyα-forest spikes are strengths that could be leveraged for future observational tests.
major comments (2)
- [Abstract / simulation and sightline selection] Abstract and simulation description: The headline probabilities (∼3×10^{-3} at 80% ionization, zero at 60%) and the 10^4 sensitivity to local background both rest on correctly counting sightlines that intersect voids with ≲30% mean density. No resolution or convergence study is presented for the tail of the gas density PDF at z>6 in Cosmic Dawn III, yet n_HI ∝ n_gas² makes transmission exponentially sensitive to these rare underdensities; shifts in the counted fraction by orders of magnitude are possible if the grid scale or sub-grid physics does not converge on this statistic.
- [Lyα transmission calculation] Ionization and transmission modeling: The assumption of ionization equilibrium with a uniform background plus a locally scaled intensity may break down near patchy ionization fronts; this directly affects whether the identified voids remain transmissive and therefore underpins both the reported probabilities and the claimed sensitivity to background intensity.
minor comments (2)
- [Abstract] Notation for the density threshold (≲30% of cosmic mean) should be tied explicitly to a figure or equation showing the transmission curve versus overdensity.
- [Results] Clarify the exact criteria used to classify a profile as 'double-peaked' (e.g., peak separation, flux ratio) and whether these match the observational samples cited.
Simulated Author's Rebuttal
We thank the referee for their constructive report and for highlighting the potential of double-peaked Lyα profiles as a probe of the ionizing background. We address each major comment below, clarifying the physical robustness of our results while acknowledging where additional discussion will strengthen the manuscript.
read point-by-point responses
-
Referee: [Abstract / simulation and sightline selection] Abstract and simulation description: The headline probabilities (∼3×10^{-3} at 80% ionization, zero at 60%) and the 10^4 sensitivity to local background both rest on correctly counting sightlines that intersect voids with ≲30% mean density. No resolution or convergence study is presented for the tail of the gas density PDF at z>6 in Cosmic Dawn III, yet n_HI ∝ n_gas² makes transmission exponentially sensitive to these rare underdensities; shifts in the counted fraction by orders of magnitude are possible if the grid scale or sub-grid physics does not converge on this statistic.
Authors: We agree that the absolute probabilities are sensitive to the accurate sampling of rare underdensities because of the n_HI ∝ n_gas² scaling. Cosmic Dawn III was run at high spatial resolution specifically to capture reionization-era structure, and the voids we identify are several cMpc in extent—well above the grid scale. The central physical result, however, is the existence of such voids in the fluctuating density field and the consequent sharp dependence of transmission on local background intensity; this scaling holds regardless of the precise normalization of the probability. We will revise the methods section to include a brief discussion of resolution, referencing convergence tests already published for Cosmic Dawn III, and to note the associated uncertainty in the absolute fraction while emphasizing that the relative sensitivity to background intensity is robust. revision: yes
-
Referee: [Lyα transmission calculation] Ionization and transmission modeling: The assumption of ionization equilibrium with a uniform background plus a locally scaled intensity may break down near patchy ionization fronts; this directly affects whether the identified voids remain transmissive and therefore underpins both the reported probabilities and the claimed sensitivity to background intensity.
Authors: The concern is well taken for sightlines that pass close to ionization fronts. In the cases that produce double-peaked profiles, however, the transmissive voids lie several comoving megaparsecs from the source galaxy, where the ionization field is smoother and the equilibrium assumption is more secure at the 60–80 % global ionization fractions we examine. The local background scaling is applied only to capture the enhanced ionization near the galaxy itself. We will expand the methods section to state the domain of validity of the equilibrium approximation explicitly, to discuss possible non-equilibrium effects near fronts, and to confirm that our selected sightlines avoid those immediate regions. revision: yes
Circularity Check
No significant circularity; results are direct simulation counts
full rationale
The paper computes Lyα transmission by ray-tracing through the independent Cosmic Dawn III simulation, identifies sightlines intersecting underdense voids, and reports the fraction of double-peaked profiles as a direct count (~3e-3 at 80% ionization). This is forward modeling with no fitted parameters renamed as predictions, no self-definitional loops, and no load-bearing self-citations or imported uniqueness theorems. The sensitivity to local ionizing background is likewise measured by varying the background in the same simulation setup. The derivation chain is self-contained against external simulation data.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Neutral hydrogen density scales with the square of the underlying gas density under ionization equilibrium
Reference graph
Works this paper leans on
-
[1]
2010, A&A, 523, A4
Baek, S., Semelin, B., Di Matteo, P., Revaz, Y., & Combes, F. 2010, A&A, 523, A4
2010
-
[2]
2001, PhR, 349, 125
Barkana, R., & Loeb, A. 2001, PhR, 349, 125
2001
-
[3]
D., & Bolton, J
Becker, G. D., & Bolton, J. S. 2013, MNRAS, 436, 1023
2013
-
[4]
D., Bolton, J
Becker, G. D., Bolton, J. S., & Lidz, A. 2015, PASA, 32, e045
2015
-
[5]
D., D’Aloisio, A., Christenson, H
Becker, G. D., D’Aloisio, A., Christenson, H. M., et al. 2021, MNRAS, 508, 1853
2021
-
[6]
2023, MNRAS, 523, 3749
Blaizot, J., Garel, T., Verhamme, A., et al. 2023, MNRAS, 523, 3749
2023
-
[7]
Bosman, S. E. I., Davies, F. B., Becker, G. D., et al. 2022, MNRAS, 514, 55
2022
-
[8]
2013, Reports on Progress in Physics, 76, 112901
Bromm, V. 2013, Reports on Progress in Physics, 76, 112901
2013
-
[9]
2020, A&A, 642, L16
Byrohl, C., & Gronke, M. 2020, A&A, 642, L16
2020
-
[10]
P., Becker, G
Calverley, A. P., Becker, G. D., Haehnelt, M. G., & Bolton, J. S. 2011, MNRAS, 412, 2543 D’Aloisio, A., McQuinn, M., Davies, F. B., & Furlanetto, S. R. 2018, MNRAS, 473, 560
2011
-
[11]
B., & Furlanetto, S
Davies, F. B., & Furlanetto, S. R. 2016, MNRAS, 460, 1328
2016
-
[12]
B., Bosman, S
Davies, F. B., Bosman, S. E. I., D’Odorico, V., et al. 2026, MNRAS, 545, staf1862
2026
-
[13]
E., et al
Dayal, P., Volonteri, M., Greene, J. E., et al. 2025, A&A, 697, A211
2025
-
[14]
2019, A&A, 622, A142
Deparis, N., Aubert, D., Ocvirk, P., Chardin, J., & Lewis, J. 2019, A&A, 622, A142
2019
-
[15]
L., Schmidt, L
DePoy, D. L., Schmidt, L. M., Ribeiro, R., et al. 2018, in Ground-based and Airborne Instrumentation for Astronomy VII, ed. C. J. Evans, L. Simard, & H. Takami, Vol. 10702, International Society for Optics and Photonics (SPIE), 107021X
2018
-
[16]
2014, PASA, 31, e040
Dijkstra, M. 2014, PASA, 31, e040
2014
-
[17]
Dijkstra, M., Mesinger, A., & Wyithe, J. S. B. 2011, MNRAS, 414, 2139
2011
-
[18]
B., Ciardi, B., Graziani, L., et al
Eide, M. B., Ciardi, B., Graziani, L., et al. 2020, MNRAS, 498, 6083
2020
-
[19]
B., & Hennawi, J
Eilers, A.-C., Davies, F. B., & Hennawi, J. F. 2018, ApJ, 864, 53
2018
-
[20]
J., & Stanway, E
Eldridge, J. J., & Stanway, E. R. 2022, ARA&A, 60, 455
2022
-
[21]
P., Bouwens, R
Endsley, R., Stark, D. P., Bouwens, R. J., et al. 2022, MNRAS, 517, 5642
2022
-
[22]
K., Steidel, C
Erb, D. K., Steidel, C. C., Trainor, R. F., et al. 2014, ApJ, 795, 33
2014
-
[23]
K., Strauss, M
Fan, X., Narayanan, V. K., Strauss, M. A., et al. 2002, AJ, 123, 1247
2002
-
[24]
A., Becker, R
Fan, X., Strauss, M. A., Becker, R. H., et al. 2006, AJ, 132, 117
2006
-
[25]
L., D’Aloisio, A., Paardekooper, J.-P., et al
Finkelstein, S. L., D’Aloisio, A., Paardekooper, J.-P., et al. 2019, ApJ, 879, 36
2019
-
[26]
2013, ApJ, 764, 41
Fragos, T., Lehmer, B., Tremmel, M., et al. 2013, ApJ, 764, 41
2013
-
[27]
R., Zaldarriaga, M., & Hernquist, L
Furlanetto, S. R., Zaldarriaga, M., & Hernquist, L. 2004, ApJ, 613, 1
2004
-
[28]
G., Davies, F
Gaikwad, P., Haehnelt, M. G., Davies, F. B., et al. 2023, MNRAS, 525, 4093 15
2023
-
[29]
2021, MNRAS, 504, 1902
Garel, T., Blaizot, J., Rosdahl, J., et al. 2021, MNRAS, 504, 1902
2021
-
[30]
Y., & Madau, P
Gnedin, N. Y., & Madau, P. 2022, Living Reviews in Computational Astrophysics, 8, 3
2022
-
[31]
2021, MNRAS, 508, 3697
Gronke, M., Ocvirk, P., Mason, C., et al. 2021, MNRAS, 508, 3697
2021
-
[32]
2018, MNRAS, 473, 227
Hassan, S., Dav´ e, R., Mitra, S., et al. 2018, MNRAS, 473, 227
2018
-
[33]
A., Scarlata, C., Verhamme, A., & Erb, D
Henry, A., Berg, D. A., Scarlata, C., Verhamme, A., & Erb, D. 2018, ApJ, 855, 96
2018
-
[34]
M., Cowie, L
Hu, E. M., Cowie, L. L., Songaila, A., et al. 2016, ApJL, 825, L7
2016
-
[35]
T., Mellema, G., Pen, U.-L., et al
Iliev, I. T., Mellema, G., Pen, U.-L., et al. 2006, MNRAS, 369, 1625
2006
-
[36]
T., Mellema, G., Shapiro, P
Iliev, I. T., Mellema, G., Shapiro, P. R., et al. 2012, MNRAS, 423, 2222
2012
-
[37]
2008, MNRAS, 391, 63
Pen, U.-L. 2008, MNRAS, 391, 63
2008
-
[38]
2013, MNRAS, 428, 1366
Jensen, H., Laursen, P., Mellema, G., et al. 2013, MNRAS, 428, 1366
2013
-
[39]
C., Bunker, A
Jones, G. C., Bunker, A. J., Saxena, A., et al. 2024, A&A, 683, A238
2024
-
[40]
L., Dickinson, M., et al
Jung, I., Finkelstein, S. L., Dickinson, M., et al. 2020, ApJ, 904, 144
2020
-
[41]
2025, ApJS, 278, 33
Kageura, Y., Ouchi, M., Nakane, M., et al. 2025, ApJS, 278, 33
2025
-
[42]
Laursen, P., Sommer-Larsen, J., & Razoumov, A. O. 2011, ApJ, 728, 52
2011
-
[43]
Lewis, J. S. W., Ocvirk, P., Dubois, Y., et al. 2023, MNRAS, 519, 5987
2023
-
[44]
Lewis, J. S. W., Ocvirk, P., Sorce, J. G., et al. 2022, MNRAS, 516, 3389
2022
-
[45]
Loeb, A., & Furlanetto, S. R. 2013, The First Galaxies in the Universe (Princeton University Press) Luki´ c, Z., Stark, C. W., Nugent, P., et al. 2015, MNRAS, 446, 3697
2013
-
[46]
2015, ApJL, 813, L8
Madau, P., & Haardt, F. 2015, ApJL, 813, L8
2015
-
[47]
2018, A&A, 619, A136
Matthee, J., Sobral, D., Gronke, M., et al. 2018, A&A, 619, A136
2018
-
[48]
D., Mesinger, A., & D’Odorico, V
McGreer, I. D., Mesinger, A., & D’Odorico, V. 2015, MNRAS, 447, 499
2015
-
[49]
2007, MNRAS, 377, 1043
McQuinn, M., Lidz, A., Zahn, O., et al. 2007, MNRAS, 377, 1043
2007
-
[50]
2015, MNRAS, 446, 566
Mesinger, A., Aykutalp, A., Vanzella, E., et al. 2015, MNRAS, 446, 566
2015
-
[51]
Mesinger, A., Ferrara, A., & Spiegel, D. S. 2013, MNRAS, 431, 621
2013
-
[52]
2021, MNRAS, 500, 558 Miralda-Escud´ e, J
Garel, T. 2021, MNRAS, 500, 558 Miralda-Escud´ e, J. 1998, ApJ, 501, 15
2021
-
[53]
2025, arXiv e-prints, arXiv:2510.18398, doi: 10.48550/arXiv.2510.18398
Mukherjee, T., Zafar, T., Nanayakkara, T., et al. 2025, arXiv e-prints, arXiv:2510.18398
-
[54]
P., Tacchella, S., Mason, C
Naidu, R. P., Tacchella, S., Mason, C. A., et al. 2020, ApJ, 892, 109
2020
-
[55]
The THESAN project: Lyman-alpha emitters as probes of ionized bubble sizes
Neyer, M., Smith, A., Vogelsberger, M., et al. 2025, arXiv e-prints, arXiv:2510.18946
work page internal anchor Pith review Pith/arXiv arXiv 2025
-
[56]
2024, MNRAS, 531, 2943
Neyer, M., Smith, A., Kannan, R., et al. 2024, MNRAS, 531, 2943
2024
-
[57]
2019, A&A, 626, A77
Ocvirk, P., Aubert, D., Chardin, J., Deparis, N., & Lewis, J. 2019, A&A, 626, A77
2019
-
[58]
Ocvirk, P., Lewis, J. S. W., Gillet, N., et al. 2021, MNRAS, 507, 6108
2021
-
[59]
R., et al
Ocvirk, P., Gillet, N., Shapiro, P. R., et al. 2016, MNRAS, 463, 1462
2016
-
[60]
G., et al
Ocvirk, P., Aubert, D., Sorce, J. G., et al. 2020, MNRAS, 496, 4087
2020
-
[61]
Ocvirk, P., Lewis, J. S. W., Conaboy, L., et al. 2025, A&A, 703, A98
2025
-
[62]
2007, ApJS, 170, 335
Page, L., Hinshaw, G., Komatsu, E., et al. 2007, ApJS, 170, 335
2007
-
[63]
2025, MNRAS, 544, 4246
Park, H., Song, H., Byrohl, C., et al. 2025, MNRAS, 544, 4246
2025
-
[64]
2021, ApJ, 922, 263
Park, H., Jung, I., Song, H., et al. 2021, ApJ, 922, 263
2021
-
[65]
S., & McLure, R
Parsa, S., Dunlop, J. S., & McLure, R. J. 2018, MNRAS, 474, 2904 Pell´ o, R., Puech, M., Prieto,´E., et al. 2024, in Ground-based and Airborne Instrumentation for Astronomy X, ed. J. J. Bryant, K. Motohara, & J. R. D
2018
-
[66]
Robertson, B. E. 2022, ARA&A, 60, 121
2022
-
[67]
E., Ellis, R
Robertson, B. E., Ellis, R. S., Furlanetto, S. R., & Dunlop, J. S. 2015, ApJL, 802, L19
2015
-
[68]
Santos, M. R. 2004, MNRAS, 349, 1137
2004
-
[69]
2022, MNRAS, 512, 3243
Smith, A., Kannan, R., Garaldi, E., et al. 2022, MNRAS, 512, 3243
2022
-
[70]
M., Barger, A
Songaila, A., Hu, E. M., Barger, A. J., et al. 2018, ApJ, 859, 91
2018
-
[71]
P., Ellis, R
Stark, D. P., Ellis, R. S., & Ouchi, M. 2011, ApJL, 728, L2
2011
-
[72]
P., et al
Torralba, A., Matthee, J., Naidu, R. P., et al. 2024, A&A, 689, A44
2024
-
[73]
2007, ApJ, 671, 1
Trac, H., & Cen, R. 2007, ApJ, 671, 1
2007
-
[74]
F., Steidel, C
Trainor, R. F., Steidel, C. C., Strom, A. L., & Rudie, G. C. 2015, ApJ, 809, 89
2015
-
[75]
2011, A&A, 536, A105 16Park et al
Vernet, J., Dekker, H., D’Odorico, S., et al. 2011, A&A, 536, A105 16Park et al
2011
-
[76]
A., Iliev, I
Watson, W. A., Iliev, I. T., D’Aloisio, A., et al. 2013, MNRAS, 433, 1230
2013
-
[77]
H., Haehnelt, M
Weinberger, L. H., Haehnelt, M. G., & Kulkarni, G. 2019, MNRAS, 485, 1350
2019
-
[78]
2023, MNRAS, 521, 4356
Xu, C., Smith, A., Borrow, J., et al. 2023, MNRAS, 521, 4356
2023
-
[79]
2016, ApJ, 820, 130
Yang, H., Malhotra, S., Gronke, M., et al. 2016, ApJ, 820, 130
2016
-
[80]
2020, ApJ, 904, 26
Yang, J., Wang, F., Fan, X., et al. 2020, ApJ, 904, 26
2020
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.