{"id":"dab9a303-7314-4e6d-a2c4-3044bcc42f61","arxiv_id":"2411.17094","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A halo-model formula for the 1D power spectrum of the 21-cm forest is presented and shown to match small-scale simulations built from the same gas and temperature assumptions.","lead":"This paper builds an analytical halo-model calculation for the one-dimensional power spectrum of the 21-cm forest, the absorption signal imprinted by neutral hydrogen on quasar spectra during the epoch of reionization. It could let future Square Kilometre Array observations extract dark-matter particle mass and early heating history without expensive simulations.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central mapping is unvalidated: Eq. (4) drops the velocity-gradient term and assumes T_S=T_K, and the Fig. 4 simulation uses the same Eq. (4), so agreement only checks internal consistency, not whether P21 predicts the observable forest.","rationale":"The reader's CONDITIONAL verdict is justified. I considered possible internal errors, including the use of a single concentration-mass relation across CDM and WDM, the projection formula in Eq. (6), and the treatment of sample variance in the forecast, but none is as load-bearing as the unvalidated brightness-temperature mapping. The halo-model integrals themselves are standard and likely correct; the problem is that the only validation target is built from the same approximations. The neglect of the H/(dv_parallel/dr_parallel) factor and of incomplete T_S coupling is asserted rather than tested, and if either is wrong the amplitude of P(k) changes directly, biasing both parameter forecasts. This is not a reason to reject the paper, since the model is explicit and testable, but it is a reason to require an external test before accepting the model as a description of the observable 21-cm forest. The recommendation remains conditional, pending a quantitative test of Eq. (4) or an independent comparison simulation.","tokens_in":18530,"tokens_out":12282,"duration_ms":125420,"concrete_test":"Recompute the 1D power spectrum from the same 2 Mpc simulation boxes used for Fig. 4, but with the full optical depth of Eq. (2) instead of Eq. (4): construct a velocity field from the same infall/NFW profiles (Hubble plus radial peculiar velocity, plus a thermal broadening term) and compute T_S from a Ly-alpha coupling model rather than setting T_S = T_K. Compare the resulting P(k) over 3 < k < 100 Mpc^-1 with the current solid curves. If the amplitude shifts by more than about 20% or the k-dependence changes, then the terms neglected in Eq. (4) are load-bearing and the m_W/f_X forecasts are not yet reliable.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Eqs. (1)-(2) contain the factor H/(dv∥/dr∥) and depend on T_S; Eq. (4) becomes the model's fundamental field only after dropping that factor, taking x_HI = 1, T_gamma >> T_S, and full Ly-alpha coupling T_S = T_K. These are not innocent assumptions: at k > 3 Mpc^-1, minihalo velocity dispersion and thermal broadening redistribute absorption in frequency, and incomplete Ly-alpha coupling places T_S between T_K and T_gamma, both changing the amplitude and slope of the 1D power spectrum. The simulation used for validation in Sec. IV is generated from the same prescriptions - the same NFW/infall density profiles, the same virial/adiabatic-plus-x-ray temperature prescription, and the same brightness-temperature mapping - so Fig. 4 confirms only that the analytic one-halo integral reproduces a Monte Carlo realization of the same sub-grid model. The conclusion's 'high consistency with simulation results' therefore does not validate the mapping from P21 to the actual observable spectrum, and the forecast errors for m_W and f_X assume that mapping is exact. The paper asserts that the redshift-space distortion effect is 'relatively weak' without a quantitative test.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops an analytical halo-model description of the 1D power spectrum of the 21-cm forest. Starting from the standard emission/absorption formalism, the authors simplify the brightness temperature to δTb ≈ T0(1+δ)/T_K, model the gas density and temperature profiles around halos, and compute the 1D power spectrum by projecting the 3D halo-model power spectrum. The model is compared against the small-scale simulation of Ref. [32] for CDM and WDM models with varying x-ray heating, showing visually good agreement on scales below 2 Mpc. The authors then use the model to forecast constraints on the warm dark matter mass m_W and the x-ray heating efficiency f_X from 10 radio-loud quasars observed with SKA-LOW for 100 hours each, obtaining σ_mW ≈ 1.3 keV and σ_fX ≈ 0.02.","tokens_in":18756,"tokens_out":6419,"duration_ms":57346,"significance":"The halo-model machinery is well suited to the scales of interest for the 21-cm forest, and the paper provides explicit closed-form expressions and a transparent forecasting pipeline, including a quasar population model and a noise estimate. If the underlying brightness-temperature mapping is accepted, this is a useful tool for quick parameter exploration and survey optimization. The main strength is the clear construction of the one-halo and two-halo terms and the demonstration that the one-halo term dominates at small scales, which physically maps the observable onto the halo mass function. However, the validation is internal: the simulation used for comparison is built from the same density and temperature prescriptions and the same brightness-temperature mapping, so the agreement checks the consistency of the analytic integration with a Monte Carlo realization of the same subgrid model rather than the accuracy of the physical mapping itself.","major_comments":[{"comment":"The simplification δTb ≈ T0(1+δ)/T_K drops the velocity-gradient factor H/(dv∥/dr∥) from Eq. (2) and assumes full Lyα coupling T_S = T_K. The validation simulation in Section IV is constructed from the same NFW/infall density profiles, the same virial/adiabatic-plus-x-ray temperature prescription, and the same brightness-temperature mapping. Therefore Figs. 4 and 5 demonstrate only that the analytic one-halo integral reproduces a Monte Carlo realization of the same subgrid model; they do not test whether Eq. (4) accurately represents the observable 21-cm forest. At k > 3 Mpc^-1, minihalo velocity dispersion and thermal broadening redistribute absorption in frequency, and incomplete Lyα coupling places T_S between T_K and T_γ, both of which would change the amplitude and slope of P(k∥). Please either validate Eq. (4) against a simulation that includes the velocity-gradient term and a more complete T_S treatment, or provide a quantitative estimate of the magnitude of these effects on P(k∥) and state the resulting systematic uncertainty in the forecast.","section":"§II.A, Eq. (4) and §IV"},{"comment":"The comparison between the halo model and the simulation is presented without error bars or any goodness-of-fit statistic. Since the claim that the model 'exhibits high consistency with the simulation results across a wide range of parameter spaces' is central to the paper, please add the simulation's sample variance (the text states that 500^2 lines of sight are averaged) and report a quantitative agreement metric, such as the fractional difference or χ² per degree of freedom, for each panel. This will allow readers to assess whether the visual agreement is genuine or simply the absence of error bars, and it will make the 'applicable across a wide range of parameters' claim falsifiable.","section":"§IV, Figs. 4 and 5"},{"comment":"The temperature profile is a free subgrid prescription: T_K equals the virial temperature inside halos and the maximum of the adiabatic and x-ray-heated temperatures outside, with f_X as a free parameter. The HERA constraint is used only to set the fiducial f_X = 0.1, not to validate the profile shape. Because the predicted power spectrum amplitude scales roughly as 1/T_K^2 through Eq. (4), a biased temperature profile would directly bias the m_W and f_X forecasts. Please add a sensitivity test that varies the temperature-profile shape (e.g., using the HERA 1σ bounds or a different x-ray heating model) and shows the impact on the forecast constraints, so that the systematic error is quantified rather than implicit.","section":"§II.B.3 and §III.B"}],"minor_comments":[{"comment":"The symbol 'M4' is used in the text but not defined; please define it explicitly as the mass at which the virial temperature equals 10^4 K.","section":"§II.B.1"},{"comment":"The normalization of the gas density profile in Eq. (16) is not specified; please state the boundary condition that fixes ρ_gc, since the profile is needed for the Fourier transform in Eq. (9).","section":"§II.B.2"},{"comment":"The survey area expression integrates sinθ dθ without specifying the integration limits; for a declination band from -86° to 34°, the area is obtained by integrating over the corresponding polar angle range, and the stated 10,000 deg^2 result should be verified with the actual limits.","section":"§III.A, Eq. (22)"},{"comment":"The phrase 'theoretically modeling' in the conclusion should be 'theoretically model' or 'theoretical modeling'; a typographical error.","section":"§V"},{"comment":"The legend repeats 'Halo model' and 'Simulation' for each color, which is redundant; please make the legend concise and associate the colors with f_X values in the caption.","section":"Fig. 4"}],"recommendation":"major_revision","confidential_remarks":"The validation simulation is not independent: Ref. [32] shares two authors with this manuscript and is built from the same density, temperature, and brightness-temperature prescriptions. This is not a conflict of interest, but it means the comparison in Figs. 4 and 5 is an internal consistency check rather than a physical validation. I would strongly encourage the authors to either use an independent simulation or explicitly reframe the comparison as a test of the halo-model implementation against a Monte Carlo realization of the same subgrid model. The paper's novelty claim of being the 'first' analytical halo-model description of the 21-cm forest should also be softened in light of prior halo-model 21-cm power spectrum work (Refs. [43-47]), unless the distinction is clearly articulated."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is the first analytic halo-model treatment of the 1D 21-cm forest power spectrum as far as I can tell. That is a real but modest step. Previous halo-model 21-cm work targeted the 3D power spectrum; here the projection in Eq. (6) and the one-halo-dominated small-scale limit are straightforward, and the SKA forecast pipeline is concrete. The math in Eqs. (5)-(8) is standard, and the curves in Figs. 4-5 track the simulation reasonably well by eye.\n\nWhat the paper does well: it writes down a complete, closed-form model with clearly stated ingredients (Press-Schechter mass function, NFW/infall density profiles, virial/adiabatic-plus-x-ray temperature prescription) and checks the one-halo term against small-scale simulations. The Bayesian forecast using 10 radio-loud quasars and the resulting sigma_mW ~ 1.3 keV, sigma_fX ~ 0.02 are internally consistent and give the community a useful first estimate of what SKA-LOW can do.\n\nThe soft spots are real, and the stress-test note lands. The biggest one: the validation in Fig. 4 is mostly internal consistency. The simulation uses the same NFW/infall density profiles, the same virial and x-ray temperature prescription, and the same brightness-temperature mapping, Eq. (4), as the analytic model; the simulation code comes from a paper sharing two authors. So agreement shows the halo-model integrals reproduce a Monte Carlo realization of the same sub-grid model. It does not validate Eq. (4) as a description of the actual 21-cm forest. Eq. (4) drops the velocity-gradient factor and assumes T_S = T_K, and the paper's claim that redshift-space distortions are 'relatively weak' is asserted, not tested. At k > 3 Mpc^-1, minihalo velocity dispersion and thermal broadening can move absorption in frequency, and incomplete Ly-alpha coupling would change amplitude and slope. The forecasts inherit those assumptions.\n\nA couple of smaller points: Figs. 4-5 have no error bars or residual/goodness-of-fit metric, and code/data are not public. For a methods paper, a public code release would let people check the internal-consistency issue quickly.\n\nNone of this is fatal. The paper is honest about its limitations in places, mentioning ionization, low-mass galaxy feedback, and Ly-alpha coupling as future work, and the analytical framework is sensible. The conclusion's 'high consistency with the simulation results' overreaches, but a revised version with independent validation, quantified residuals, and public code would be a solid contribution to 21-cm cosmology. My recommendation is to send it to peer review rather than desk reject, with major revision targeted at the validation strategy.","headline":"First analytic halo-model treatment of the 1D 21-cm forest power spectrum, with a useful SKA forecast, but the validation is mostly internal consistency and the forecast errors inherit unvalidated assumptions about spin temperature and velocity gradients.","tokens_in":19315,"tokens_out":3630,"would_cite":true,"duration_ms":33932,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper claims a halo model can analytically reproduce the 1D power spectrum of the 21-cm forest, with the one-halo term dominating below 2 Mpc, enabling SKA-LOW forecasts of warm dark matter and x-ray heating.","keywords":["21-cm forest","halo model","one-dimensional power spectrum","warm dark matter","x-ray heating","epoch of reionization","SKA-LOW","radio-loud quasars"],"falsifier":"Run a small-box radiative-transfer simulation that computes the spin temperature from the local Lyman-$\\alpha$ and x-ray fields instead of assuming $T_S = T_K$, and compare its 1D power spectrum at $k > 3$ Mpc$^{-1}$ with the halo-model prediction for the same halo population; a difference larger than the forecast measurement error $\\sigma_P$ would mean the quoted sensitivities are biased. A simpler observational check would be a SKA-LOW measurement toward a $z \\approx 8$ quasar that tests whether the predicted one-halo slope and amplitude are present.","tokens_in":18234,"feed_emoji":"📡","tokens_out":13031,"duration_ms":101216,"temperature":0.7,"pith_summary":"This paper tries to establish that the one-dimensional (1D) power spectrum of the 21-cm forest—the dense web of absorption lines that neutral hydrogen imprints on the spectra of high-redshift radio sources—can be computed analytically with a halo model, eliminating the need for costly small-box simulations. The model separates the signal into a one-halo term and a two-halo term; on scales below about 2 comoving Mpc the one-halo term dominates, so the observed spectrum is essentially the dark-matter halo mass function weighted by gas density and temperature profiles around halos. The authors compare this analytic spectrum against small-scale simulation results for both cold and warm dark matter and for several x-ray heating efficiencies, and report close agreement across that parameter space. If the model holds, it gives a fast, parameter-ready bridge from future Square Kilometre Array observations to two early-Universe unknowns: the mass of warm dark matter particles and the efficiency of x-ray heating.","feed_headline":"Halo model reproduces 21-cm forest spectrum; SKA can weigh dark matter","feed_subtitle":"New analytical model links absorption lines to halo counts, letting SKA-LOW probe dark matter and cosmic heating.","key_machinery":"The carrying object is the halo-model decomposition applied to the 21-cm brightness-temperature contrast, specifically to the profile $\\rho_{21}(r,z,M) = [1+\\delta(r)]/T_K(r)$ around a halo of mass $M$. A window function $W_{21}(k,z,M)$, the Fourier transform of $\\rho_{21}$, is integrated over the halo mass function $dn/dM$ to form the one-halo term, and the same window function weighted by halo bias is combined with the linear matter power spectrum for the two-halo term. The radial profiles come from an NFW dark-matter profile with gas in hydrostatic equilibrium inside the virial radius and an infall model outside; the temperature is the virial temperature inside halos and adiabatic plus x-ray heating outside, with $f_X$ as a free parameter. On scales below about 2 comoving Mpc the one-halo term dominates, making the predicted 1D power spectrum essentially a convolution of the halo mass function with these profiles.","core_discovery":"The paper's central claim is that the 1D power spectrum of the 21-cm forest along the line of sight is, up to a background-source prefactor $T_0^2(\\hat n,z)$, the projected 3D power spectrum of the field $[1+\\delta(z)]/T_K(z)$, where $\\delta$ is the gas overdensity and $T_K$ is the gas kinetic temperature. In the halo-model decomposition this becomes $P_{21}(k_\\parallel,z) = T_0^2 [P^{1h}_{21} + P^{2h}_{21}]$, with the one-halo term built from the Fourier-transformed halo profile $W_{21}(k,z,M)$ integrated over the halo mass function and the two-halo term built from the halo bias times the linear matter power spectrum. On the small scales the model targets, the one-halo term dominates, so the spectral shape is fixed by the halo mass function and the assumed density and temperature profiles. After validating the model against small-scale simulations for CDM and 6 keV WDM at $f_X = 0.01$, $0.1$, and $1$, the paper forecasts that SKA-LOW observations of ten radio-loud quasars, 100 hours each, would constrain a fiducial $m_W = 6$ keV to about $\\pm 1.3$ keV and a fiducial $f_X = 0.1$ to about $\\pm 0.02$.","pith_inferences":["Editorial extension: the one-halo dominance implies the 21-cm forest 1D power spectrum could in principle be inverted to measure the halo mass function over roughly $10^5$-$10^8\\,M_\\odot$, providing a small-scale structure census independent of galaxy surveys.","Editorial extension: the forecast's error bars assume the model's own temperature prescription in the likelihood; a self-consistent treatment of Lyman-$\\alpha$ coupling and inhomogeneous heating would provide a direct stress test of those sensitivities.","Editorial extension: the same halo-model pipeline should transfer to other small-scale dark-matter models with a cutoff, such as ultra-light axions, and to other one-dimensional statistics of the absorption spectra, like the variance or bispectrum."],"forward_implications":["The 1D 21-cm forest power spectrum on small scales can be evaluated directly from the halo mass function and profiles, bypassing expensive small-box simulations for survey design and parameter studies.","Because the one-halo term dominates below roughly 2 comoving Mpc, the observed small-scale spectrum traces the abundance of low-mass halos, so a deficit of power is a direct warm-dark-matter signature.","The forecast implies that SKA-LOW observations of ten bright radio-loud quasars could constrain the warm dark matter particle mass and the x-ray heating efficiency simultaneously.","Extending the model to larger scales requires adding ionization fluctuations and Lyman-$\\alpha$ coupling, a step the authors identify as the natural next development."],"supporting_citations":[{"why":"Supplies the halo-model formalism that splits the power spectrum into one-halo and two-halo terms.","marker":"[40]"},{"why":"Provides the small-scale 21-cm forest simulation results that the analytical model is compared against.","marker":"[32]"},{"why":"Gives the 21-cm optical depth and absorption signal formalism underlying the simplified brightness-temperature relation.","marker":"[14]"},{"why":"Supplies the original 21-cm optical depth formula used in Eq. (2).","marker":"[52]"},{"why":"Provides the warm dark matter halo mass function with free-streaming suppression used in the model.","marker":"[27]"},{"why":"Supplies the NFW dark-matter density profile assumed inside the virial radius.","marker":"[64]"},{"why":"Gives the hydrostatic-equilibrium gas density profile inside halos.","marker":"[65]"},{"why":"Supplies the IGM temperature evolution with Compton and x-ray heating.","marker":"[68]"},{"why":"Provides the updated radio-loud quasar distribution used to simulate SKA-LOW background sources.","marker":"[51]"},{"why":"Provides the 1D power spectrum projection and thermal-noise formulas used in the forecast.","marker":"[54]"}],"fun_headline_variants":["Halo model computes 21-cm forest spectrum; SKA to test warm dark matter","Analytical halo model turns 21-cm forest into dark matter probe","SKA-LOW could constrain warm dark matter via 21-cm forest","21-cm forest halo model forecasts dark matter mass and heating"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the spin temperature of neutral hydrogen is fully coupled to the gas kinetic temperature through an early Lyman-$\\alpha$ background and that the CMB dominates the background radiation, so the 21-cm brightness is just proportional to $[1+\\delta]/T_K$; if that coupling is incomplete or the assumed gas temperature is wrong, the predicted spectrum shifts, and because the validation simulation uses the same temperature and density prescriptions, it cannot detect that failure.","fun_headline_variants_meta":{"raw":{"variants":["Halo model computes 21-cm forest spectrum; SKA to test warm dark matter","Analytical halo model turns 21-cm forest into dark matter probe","SKA-LOW could constrain warm dark matter via 21-cm forest","21-cm forest halo model forecasts dark matter mass and heating"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000339,"raw_usage":{"total_tokens":1954,"prompt_tokens":1112,"completion_tokens":842,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":728,"completion_tokens_details":{"reasoning_tokens":760}},"tokens_in":728,"tokens_out":842,"duration_ms":8067,"temperature":1.0,"reasoning_tokens":760,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:35:53.438102+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a small-box radiative-transfer simulation that computes the spin temperature from the local Lyman-$\\alpha$ and x-ray fields instead of assuming $T_S = T_K$, and compare its 1D power spectrum at $k > 3$ Mpc$^{-1}$ with the halo-model prediction for the same halo population; a difference larger than the forecast measurement error $\\sigma_P$ would mean the quoted sensitivities are biased. A simpler observational check would be a SKA-LOW measurement toward a $z \\approx 8$ quasar that tests whether the predicted one-halo slope and amplitude are present.","supporting_citations":[],"review_version":1}