On Cross-Correlating Line Intensity Maps from SPHEREx during Reionization
Pith reviewed 2026-06-27 15:54 UTC · model grok-4.3
The pith
Cross-correlating SPHEREx intensity maps for Hα and [OIII] at z=5 yields total signal-to-noise up to 99 after interloper masking.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
By simulating SPHEREx intensity maps for multiple emission lines during reionization and computing their cross-power spectra after accounting for interlopers and masking, the authors find that the highest signal-to-noise is 99 for Hα cross [OIII] at z=5 with no dust extinction, that intensity mapping can access galaxies with M < 4×10^{10} M⊙ below the 3σ direct detection threshold, but the majority of the observable signal is dominated by large directly detectable galaxies rather than smaller fainter ones, with marginal clustering detections possible at z=5 for that pair.
What carries the argument
Cross-power spectra of simulated multi-line intensity maps after modeling interloping lines from lower-redshift galaxies and masking bright interlopers.
If this is right
- Intensity mapping accesses galaxies with masses M<4×10^{10} M⊙ below the 3σ direct detection threshold of SPHEREx.
- Instrument noise is the dominant noise source on most scales, except for Lyα and [OII] on the largest scales where interlopers dominate.
- The observable cross-correlation signal is dominated by large directly detectable galaxies rather than the fainter population.
- Marginal detections of the clustering portion of the power spectrum are possible at z=5 for Hα × [OIII], but other lines or z>6 require more sensitive instruments.
Where Pith is reading between the lines
- Cross-correlations could provide a statistical route to map the contribution of faint galaxies to reionization without resolving them individually.
- If the dominance by bright galaxies holds, intensity mapping adds limited new information on the faint end beyond what direct detections already capture.
- Varying the dust model or adding more detailed radiative transfer effects could be tested against these baseline S/N predictions in future simulations.
Load-bearing premise
The simple dust extinction model calibrated on z<5 galaxy observations remains accurate at z≈5–7, and unmodeled systematics or interloper masking failures do not substantially alter the cross-power spectra.
What would settle it
An actual measurement of the Hα × [OIII] cross-power spectrum at z=5 showing total signal-to-noise well below 99 after realistic interloper masking and dust effects would falsify the highest detectability prediction.
read the original abstract
We have simulated Ly{\alpha}, H{\alpha}, H{\beta}, [OII], and [OIII] intensity maps which are observable by SPHEREx during cosmic reionization. We simulate these intensity maps including all significant sources of emission for each line, and include radiative transfer for the Ly{\alpha} intensity maps. We also include a simple model of dust extinction based on observations of galaxies at z<5. One of the main challenges of intensity mapping is interloping lines from galaxies at lower redshifts, which makes producing an auto-power spectrum challenging. We focus on cross-correlations between different lines, as this eliminates such foreground contamination of the signal. We have cross-correlated the simulated SPHEREx intensity maps to find the most observable cases. This includes modeling of interloping lines and masking bright interloping galaxies. Testing a range of cases motivated by observations, we find total signal-to-noise values up to 99 for the highest case of H{\alpha} cross-correlated with [OIII] at z=5 assuming no dust extinction. We also find cases which will not be detectable. We find that the dominant noise source in these intensity maps on most scales is from the instrument, except for Ly{\alpha} and [OII] and then only on the largest scales the interlopers are the dominant source. We find through intensity mapping we can probe galaxies with masses $M<4x10^{10}~M_{\odot}$ which are below the necessary luminosity for a 3{\sigma} signal-to-noise direct detection of galaxies by SPHEREx. However, the majority of our observable signal is dominated by large, directly detectable galaxies, rather than the smaller, fainter galaxies. We find marginal detections of the clustering portion of the power spectrum at z=5 for H{\alpha}x[OIII]. Detections of the clustering signal from other lines or at z>6 will require more sensitive instruments, such as the Cosmic Dawn Intensity Mapper.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper simulates Lyα, Hα, Hβ, [OII], and [OIII] intensity maps for SPHEREx during reionization, incorporating all major emission sources, radiative transfer for Lyα, and a simple dust extinction model calibrated on z<5 galaxies. It computes cross-power spectra between line pairs (with interloper modeling and bright-galaxy masking) to avoid foreground contamination, reports total S/N values reaching 99 for Hα × [OIII] at z=5 in the no-dust case, identifies instrument noise as dominant on most scales, and concludes that intensity mapping can access galaxies below the M=4×10^{10} M⊙ 3σ direct-detection threshold even though the observable signal is dominated by bright, directly detectable galaxies. Marginal clustering detections are found only for the z=5 Hα×[OIII] case.
Significance. If the modeling assumptions hold, the work provides concrete forecasts for the detectability of reionization-era line intensity mapping with SPHEREx via cross-correlations, quantifies the relative roles of instrument noise versus interlopers, and illustrates the complementarity of intensity mapping to direct galaxy detection. The forward-modeling framework with explicit interloper treatment is a strength for planning future observations such as those with CDIM.
major comments (3)
- [§3.3] §3.3 (Dust model): the extinction is taken from z<5 observations and applied at z=5–7 without additional validation or sensitivity tests; because the quoted peak S/N=99 is reported only for the no-dust case, any realistic increase in optical depth at reionization redshifts would proportionally lower the cross-power amplitude and total S/N, directly affecting the detectability conclusions.
- [§5] §5 (Mass threshold claim): the assertion that intensity mapping reaches M<4×10^{10} M⊙ galaxies below the SPHEREx 3σ direct-detection limit is not accompanied by a quantitative decomposition of the cross-power into luminosity bins; given the explicit statement that the observable signal is dominated by large, directly detectable galaxies, the incremental contribution from the faint population must be shown explicitly (e.g., via cumulative signal fractions) to support the claim.
- [§4.2] §4.2 (Interloper masking and noise dominance): the conclusion that instrument noise dominates on most scales (except Lyα and [OII] on the largest scales) rests on specific masking thresholds and noise modeling choices whose variation is not explored; altering the 3σ masking cut or the dust scaling would change the relative noise budgets and therefore the reported S/N values.
minor comments (2)
- [Figure 4] Figure 4 (or equivalent S/N summary plot): axis labels and line styles for the no-dust versus dust cases are not clearly distinguished in the caption, making it difficult to map the plotted curves to the tabulated values.
- Notation: the definition of the total S/N (summed over k-bins) is introduced without an explicit equation reference, forcing the reader to infer the precise summation from the text.
Simulated Author's Rebuttal
We thank the referee for their thorough review and constructive comments. We address each major comment point by point below, indicating revisions where the manuscript will be updated.
read point-by-point responses
-
Referee: [§3.3] §3.3 (Dust model): the extinction is taken from z<5 observations and applied at z=5–7 without additional validation or sensitivity tests; because the quoted peak S/N=99 is reported only for the no-dust case, any realistic increase in optical depth at reionization redshifts would proportionally lower the cross-power amplitude and total S/N, directly affecting the detectability conclusions.
Authors: We agree that the dust model is extrapolated from z<5 data without validation or sensitivity tests at z=5–7, and that the S/N=99 applies only to the no-dust case. In the revision we will add sensitivity tests scaling the extinction (e.g., by factors of 1.5–3) and recompute the cross-power spectra and total S/N for the key line pairs. This will quantify the impact on detectability conclusions. revision: yes
-
Referee: [§5] §5 (Mass threshold claim): the assertion that intensity mapping reaches M<4×10^{10} M⊙ galaxies below the SPHEREx 3σ direct-detection limit is not accompanied by a quantitative decomposition of the cross-power into luminosity bins; given the explicit statement that the observable signal is dominated by large, directly detectable galaxies, the incremental contribution from the faint population must be shown explicitly (e.g., via cumulative signal fractions) to support the claim.
Authors: We acknowledge that a quantitative decomposition by luminosity or mass bins is required to support the claim. Although the manuscript already states that the signal is dominated by bright galaxies, we will add cumulative signal-fraction plots and tables breaking down the cross-power contributions from galaxies below the 3σ direct-detection threshold in the revised version. revision: yes
-
Referee: [§4.2] §4.2 (Interloper masking and noise dominance): the conclusion that instrument noise dominates on most scales (except Lyα and [OII] on the largest scales) rests on specific masking thresholds and noise modeling choices whose variation is not explored; altering the 3σ masking cut or the dust scaling would change the relative noise budgets and therefore the reported S/N values.
Authors: We agree that the reported noise dominance and S/N values depend on the specific 3σ masking threshold and dust model, which were not varied. In the revision we will include additional tests or a dedicated discussion quantifying how changes in the masking cut and dust scaling alter the interloper versus instrument noise budgets and the resulting S/N. revision: yes
Circularity Check
No significant circularity in forward-modeled simulations
full rationale
The paper performs forward simulations of line intensity maps from assumed galaxy populations and instrument properties, applies an external dust model calibrated on z<5 observations, and computes cross-power spectra and S/N values directly from the resulting maps. No equations or steps reduce by construction to fitted inputs, self-citations, or renamed ansatzes; the central results (e.g., S/N up to 99 for Hα × [OIII]) are outputs of the simulation pipeline rather than tautological re-expressions of its inputs.
Axiom & Free-Parameter Ledger
free parameters (2)
- dust extinction scaling
- line luminosity functions and scalings
axioms (2)
- domain assumption Standard flat Lambda-CDM cosmology governs the distribution and evolution of galaxies at z~5-7.
- domain assumption All significant sources of line emission are included and radiative transfer for Lyα is modeled sufficiently accurately.
Reference graph
Works this paper leans on
-
[1]
DeBoer, A.R
D.R. DeBoer, A.R. Parsons, J.E. Aguirre, P. Alexander, Z.S. Ali, A.P. Beardsley et al., Hydrogen epoch of reionization array (hera),PASP129(2017) 045001
2017
- [2]
-
[3]
M. Santos, P. Bull, D. Alonso, S. Camera, P. Ferreira, G. Bernardi et al.,Cosmology from a SKA HI intensity mapping survey, inAdvancing Astrophysics with the Square Kilometre Array (AASKA14), p. 19, Apr., 2015, DOI [1501.03989]
Pith/arXiv arXiv 2015
-
[4]
Cleary, J
K.A. Cleary, J. Borowska, P.C. Breysse, M. Catha, D.T. Chung, S.E. Church et al.,Comap early science. i. overview,ApJ933(2022) 182
2022
- [5]
-
[6]
A. Cooray, T.-C. Chang, S. Unwin, M. Zemcov, A. Coffey, P. Morrissey et al.,Cosmic Dawn Intensity Mapper, inBulletin of the American Astronomical Society, vol. 51, p. 23, Sept., 2019, DOI [1903.03144]
Pith/arXiv arXiv 2019
-
[7]
Bernal and E.D
J.L. Bernal and E.D. Kovetz,Line-intensity mapping: theory review with a focus on star-formation lines,The Astronomy and Astrophysics Review30(2022) 5
2022
-
[8]
T.-C. Chang and A. Lidz,Line-Intensity Mapping,arXiv e-prints(2026) arXiv:2602.03011 [2602.03011]
arXiv 2026
-
[9]
Alibay, O.V
F. Alibay, O.V. Sindiy, P.A.T. Jansma, C.M. Reynerson, E.B. Rice, J. Rocca et al.,Spherex preliminary mission overview, in2023 IEEE Aerospace Conference, pp. 1–18, 2023, DOI
2023
-
[10]
E.D. Kovetz, M.P. Viero, A. Lidz, L. Newburgh, M. Rahman, E. Switzer et al.,Line-Intensity Mapping: 2017 Status Report,arXiv e-prints(2017) arXiv:1709.09066 [1709.09066]
Pith/arXiv arXiv 2017
-
[11]
R.J. Bouwens, P.A. Oesch, M. Stefanon, G. Illingworth, I. Labb´ e, N. Reddy et al.,New Determinations of the UV Luminosity Functions from z 9 to 2 Show a Remarkable Consistency with Halo Growth and a Constant Star Formation Efficiency,AJ162(2021) 47 [2102.07775]
arXiv 2021
-
[12]
E. Visbal and M. McQuinn,The Impact of Neutral Intergalactic Gas on LyαIntensity Mapping during Reionization,ApJL863(2018) L6 [1807.03370]
Pith/arXiv arXiv 2018
-
[13]
A.E. Ambrose, E. Visbal, M. Kulkarni and M. McQuinn,Radiative transfer simulations of Lyα intensity mapping during cosmic reionization including sources from galaxies and the intergalactic medium,Journal of Cosmology and Astroparticle Physics2025(2025) 080 [2502.18654]
arXiv 2025
-
[14]
A.R. Pullen, O. Dor´ e and J. Bock,Intensity Mapping across Cosmic Times with the LyαLine, ApJ786(2014) 111 [1309.2295]
Pith/arXiv arXiv 2014
-
[15]
et al.,Planck 2018 results - vi
Planck Collaboration, Aghanim, N., Akrami, Y., Ashdown, M., Aumont, J., Baccigalupi, C. et al.,Planck 2018 results - vi. cosmological parameters,A&A641(2020) A6. – 24 –
2018
-
[16]
A. Mesinger, S. Furlanetto and R. Cen,21CMFAST: a fast, seminumerical simulation of the high-redshift 21-cm signal,MNRAS411(2011) 955 [1003.3878]
Pith/arXiv arXiv 2011
-
[17]
J. Park, A. Mesinger, B. Greig and N. Gillet,Inferring the astrophysics of reionization and cosmic dawn from galaxy luminosity functions and the 21-cm signal,MNRAS484(2019) 933 [1809.08995]
Pith/arXiv arXiv 2019
-
[18]
S.R. Furlanetto, M. McQuinn and L. Hernquist,Characteristic scales during reionization, MNRAS365(2006) 115 [astro-ph/0507524]
Pith/arXiv arXiv 2006
-
[19]
P. Madau and M. Dickinson,Cosmic Star-Formation History,ARA&A52(2014) 415 [1403.0007]
Pith/arXiv arXiv 2014
-
[20]
S.R. Furlanetto, J. Mirocha, R.H. Mebane and G. Sun,A minimalist feedback-regulated model for galaxy formation during the epoch of reionization,MNRAS472(2017) 1576 [1611.01169]
Pith/arXiv arXiv 2017
-
[21]
J.Y.-C. Yeh, A. Smith, R. Kannan, E. Garaldi, M. Vogelsberger, J. Borrow et al.,The thesan project: ionizing escape fractions of reionization-era galaxies,MNRAS520(2023) 2757
2023
-
[22]
R. Barkana and A. Loeb,Detecting the Earliest Galaxies through Two New Sources of 21 Centimeter Fluctuations,ApJ626(2005) 1 [astro-ph/0410129]
Pith/arXiv arXiv 2005
-
[23]
H. Yajima, J.-H. Choi and K. Nagamine,Escape fraction of ionizing photons from high-redshift galaxies in cosmological SPH simulations,MNRAS412(2011) 411 [1002.3346]
Pith/arXiv arXiv 2011
-
[24]
A. Ferrara and A. Loeb,Escape fraction of the ionizing radiation from starburst galaxies at high redshifts,MNRAS431(2013) 2826 [1209.2123]
Pith/arXiv arXiv 2013
-
[25]
H. Xu, J.H. Wise, M.L. Norman, K. Ahn and B.W. O’Shea,Galaxy Properties and UV Escape Fractions during the Epoch of Reionization: Results from the Renaissance Simulations,ApJ 833(2016) 84 [1604.07842]
Pith/arXiv arXiv 2016
-
[26]
A. Hoag, M. Bradaˇ c, K. Huang, C. Mason, T. Treu, K.B. Schmidt et al.,Constraining the Neutral Fraction of Hydrogen in the IGM at Redshift 7.5,ApJ878(2019) 12 [1901.09001]
Pith/arXiv arXiv 2019
-
[27]
C.A. Mason, M. Trenti and T. Treu,The Galaxy UV Luminosity Function before the Epoch of Reionization,ApJ813(2015) 21 [1508.01204]
Pith/arXiv arXiv 2015
-
[28]
M. Hayes, D. Schaerer, G. ¨Ostlin, J.M. Mas-Hesse, H. Atek and D. Kunth,On the Redshift Evolution of the LyαEscape Fraction and the Dust Content of Galaxies,ApJ730(2011) 8 [1010.4796]
Pith/arXiv arXiv 2011
-
[29]
D. Schaerer,The transition from Population III to normal galaxies: Lyalpha and He II emission and the ionising properties of high redshift starburst galaxies,A&A397(2003) 527 [astro-ph/0210462]
Pith/arXiv arXiv 2003
-
[30]
J. Blaizot, T. Garel, A. Verhamme, H. Katz, T. Kimm, L. Michel-Dansac et al.,Simulating the diversity of shapes of the Lyman-αline,MNRAS523(2023) 3749 [2305.10047]
arXiv 2023
-
[31]
A. Smith, X. Ma, V. Bromm, S.L. Finkelstein, P.F. Hopkins, C.-A. Faucher-Gigu` ere et al.,The physics of Lymanαescape from high-redshift galaxies,MNRAS484(2019) 39 [1810.08185]
Pith/arXiv arXiv 2019
- [32]
-
[33]
A. D’Aloisio, M. McQuinn, O. Maupin, F.B. Davies, H. Trac, S. Fuller et al.,Heating of the Intergalactic Medium by Hydrogen Reionization,ApJ874(2019) 154 [1807.09282]
Pith/arXiv arXiv 2019
-
[34]
Santos, V
M.R. Santos, V. Bromm and M. Kamionkowski,The contribution of the first stars to the cosmic infrared background,MNRAS336(2002) 1082
2002
- [35]
-
[36]
I.T. Iliev, G. Mellema, P.R. Shapiro and U.-L. Pen,Self-regulated reionization,MNRAS376 (2007) 534 [astro-ph/0607517]
Pith/arXiv arXiv 2007
-
[37]
McQuinn, A
M. McQuinn, A. Lidz, O. Zahn, S. Dutta, L. Hernquist and M. Zaldarriaga,The morphology of H ii regions during reionization,MNRAS377(2007) 1043
2007
-
[38]
D’Aloisio, M
A. D’Aloisio, M. McQuinn, H. Trac, C. Cain and A. Mesinger,Hydrodynamic response of the intergalactic medium to reionization,ApJ898(2020) 149
2020
-
[39]
Osterbrock and G.J
D.E. Osterbrock and G.J. Ferland,Astrophysics of gaseous nebulae and active galactic nuclei (2006)
2006
-
[40]
K. Ahn, P.R. Shapiro, I.T. Iliev, G. Mellema and U.-L. Pen,The Inhomogeneous Background Of H 2-Dissociating Radiation During Cosmic Reionization,ApJ695(2009) 1430 [0807.2254]
Pith/arXiv arXiv 2009
-
[41]
Sheth and G
R.K. Sheth and G. Tormen,Large-scale bias and the peak background split,MNRAS308 (1999) 119 [https://academic.oup.com/mnras/article-pdf/308/1/119/18409158/308-1-119.pdf]
1999
-
[42]
R.C. Kennicutt, Jr.,Star Formation in Galaxies Along the Hubble Sequence,ARA&A36 (1998) 189 [astro-ph/9807187]
Pith/arXiv arXiv 1998
-
[43]
C. Ly, M.A. Malkan, N. Kashikawa, K. Shimasaku, M. Doi, T. Nagao et al.,The Luminosity Function and Star Formation Rate between Redshifts of 0.07 and 1.47 for Narrowband Emitters in the Subaru Deep Field,ApJ657(2007) 738 [astro-ph/0610846]
Pith/arXiv arXiv 2007
-
[44]
Y. Gong, A. Cooray, M.B. Silva, M. Zemcov, C. Feng, M.G. Santos et al.,Intensity Mapping of Hα, Hβ, [OII], and [OIII] Lines at z ¡ 5,ApJ835(2017) 273 [1610.09060]
Pith/arXiv arXiv 2017
-
[45]
Cheng and T.-C
Y.-T. Cheng and T.-C. Chang,Cosmic near-infrared background tomography with spherex using galaxy cross-correlations,ApJ925(2022) 136
2022
-
[46]
D. Calzetti, L. Armus, R.C. Bohlin, A.L. Kinney, J. Koornneef and T. Storchi-Bergmann,The Dust Content and Opacity of Actively Star-forming Galaxies,ApJ533(2000) 682 [astro-ph/9911459]
Pith/arXiv arXiv 2000
-
[47]
A.A. Khostovan, D. Sobral, B. Mobasher, P.N. Best, I. Smail, J.P. Stott et al.,Evolution of the Hβ+ [O III] and [O II] luminosity functions and the [O II] star formation history of the Universe up to z∼5 from HiZELS,MNRAS452(2015) 3948 [1503.00004]
Pith/arXiv arXiv 2015
- [48]
-
[49]
C.M. Hirata,Wouthuysen-Field coupling strength and application to high-redshift 21-cm radiation,MNRAS367(2006) 259 [astro-ph/0507102]
Pith/arXiv arXiv 2006
-
[50]
J.R. Pritchard and S.R. Furlanetto,Descending from on high: Lyman-series cascades and spin-kinetic temperature coupling in the 21-cm line,MNRAS367(2006) 1057 [astro-ph/0508381]
Pith/arXiv arXiv 2006
-
[51]
S.R. Furlanetto, J. Schaye, V. Springel and L. Hernquist,LyαEmission from Structure Formation,ApJ622(2005) 7 [astro-ph/0409736]
Pith/arXiv arXiv 2005
-
[52]
R.L. Sanders, A.E. Shapley, M.W. Topping, N.A. Reddy and G.B. Brammer,Excitation and Ionization Properties of Star-forming Galaxies at z = 2.0-9.3 with JWST/NIRSpec,ApJ955 (2023) 54 [2301.06696]
arXiv 2023
-
[53]
A.E. Shapley, R.L. Sanders, M.W. Topping, N.A. Reddy, D.A. Berg, R.J. Bouwens et al.,The AURORA Survey: A New Era of Emission-line Diagrams with JWST/NIRSpec,ApJ980 (2025) 242 [2407.00157]. – 26 –
arXiv 2025
-
[54]
M. Hayashi, D. Sobral, P.N. Best, I. Smail and T. Kodama,Calibrating [O II] star formation rates at z ¡ 1 from dual Hα-[O II] imaging from HiZELS,MNRAS430(2013) 1042 [1212.4905]
Pith/arXiv arXiv 2013
-
[55]
E. Visbal and A. Loeb,Measuring the 3D clustering of undetected galaxies through cross correlation of their cumulative flux fluctuations from multiple spectral lines,JCAP2010(2010) 016 [1008.3178]
Pith/arXiv arXiv 2010
-
[56]
CAMB: Code for Anisotropies in the Microwave Background
A. Lewis and A. Challinor, “CAMB: Code for Anisotropies in the Microwave Background.” Astrophysics Source Code Library, record ascl:1102.026, Feb., 2011
2011
- [57]
- [58]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.