{"id":"6fe59d2a-93fb-47fc-9213-e2eff179768c","arxiv_id":"2505.14258","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A phenomenological recipe maps 3D matter power to 1D Lyman-alpha flux power, adds warm dark matter and pressure cutoffs in quadrature, and matches simulations to within 5-20%.","lead":"This paper proposes a simple recipe for converting the 3D matter power spectrum into the 1D Lyman-alpha flux power spectrum, combining warm dark matter and gas pressure into one cutoff scale. If the recipe holds, researchers can use Lyman-alpha forest measurements to constrain dark matter properties without running a new simulation for every model.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Validation is self-referential: k_cut is extracted and tested with the same Eq. (2.7), and the observational analysis switches to a different fitting form (E.1), so the central mapping is not independently established.","rationale":"The reader's weakest assumption correctly identifies the absence of a first-principles justification for the projected-Gaussian ansatz. My concern is a sharper version of the same issue: the validation strategy is self-referential because the same functional form is used both to define and to extract k_cut, so the reported agreement does not independently confirm that Eq. (2.7) describes the flux power spectrum. In addition, the switch to Eq. (E.1) for observational data creates an unvalidated bridge between the theoretical k_cut and the observationally fitted k_cut. These points reinforce the reader's CONDITIONAL verdict rather than overturning it: the recipe is promising and honestly caveated, but the central claim would need an out-of-sample shape test and a demonstrated equivalence between the two fitting forms before the mapping could be considered robust. The proposed concrete test would settle whether the ansatz is genuinely predictive or merely a flexible fitting function.","tokens_in":16651,"tokens_out":4934,"duration_ms":52695,"concrete_test":"Take a WDM simulation with m = 3 keV and Puchwein thermal history, compute k_cut from the input mass and thermal history using Eqs. (A.3) and (B.1), then predict the full real-space flux power spectrum from Eq. (2.7) with k_cut fixed to that theoretical value, fitting only A, a, and xi, and compare residuals to the measured spectrum over the full k range. Repeat the same comparison using an alternative cutoff shape, such as the original Viel transfer function Eq. (A.1) or a Lorentzian, with the same fixed physical scale. Finally, refit the identical simulated spectra with Eq. (2.7) and with Eq. (E.1) and compare the extracted k_cut values; if the fixed-k_cut prediction is biased, or if the two fitting forms yield k_cut values differing by more than 5 percent, the claimed correspondence is form-dependent and the observational mapping in Sec. 4.3 is not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is that Eq. (2.7) correctly maps the 3D matter power spectrum to the 1D Lyman-alpha flux power spectrum, with k_cut carrying the physical WDM and pressure filtering scale. The load-bearing weakness is that the validation never tests Eq. (2.7) as a full predictive shape; it only tests whether the value of k_cut obtained by fitting simulated flux power spectra with Eq. (2.7) matches a theory value computed from Eq. (A.3) or Eq. (B.1). Because the fitting function is the same ansatz whose validity is at issue, the reported 5 percent, 20 percent, and 15 percent agreements measure internal consistency of the extraction procedure, not correctness of the ansatz. A wrong model with enough free parameters (A, a, k_cut, xi) can return a k_cut that tracks the true filtering scale while the full predicted flux power spectrum is biased. Moreover, the observational application in Sec. 4.3 abandons Eq. (2.7) for a different fitting form, Eq. (E.1), to account for thermal broadening and redshift-space distortions; no test establishes that k_cut from Eq. (E.1) equals the k_cut appearing in Eq. (2.7). The paper itself concedes in Sec. 5 that there is no rigorous theoretical justification, and no released code or data allow an independent re-analysis of the extraction and fitting steps.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a phenomenological mapping between the 3D matter power spectrum and the 1D Lyman-alpha flux power spectrum. Starting from the exact projection formula Eq. (2.1), it assumes the 3D spectrum is a power law times a combined Gaussian cutoff (Eq. 2.6), with WDM and pressure cutoffs combined in inverse quadrature (Eq. 2.5), and it absorbs the small-scale contribution into a constant term, leading to the fitting function Eq. (2.7). The free parameters A, a, kcut, and xi are fitted to simulated real-space flux power spectra, and the extracted kcut values are compared with theoretical predictions from Viel et al. (Eq. A.3) and GH98 (Eq. B.1). The paper reports agreement within 5% for WDM, 20% for pressure, and 15% for the combined case, and it demonstrates in Sec. 4.3 that a WDM simulation with a redesigned thermal history reproduces the flux power spectrum of a CDM simulation within about 5%. The paper explicitly acknowledges in Sec. 5 that there is no rigorous derivation of the ansatz.","tokens_in":17047,"tokens_out":8037,"duration_ms":73555,"significance":"If the recipe is correct, it provides a simple analytic route from observed Lyman-alpha flux power spectrum cutoffs to constraints on WDM mass and intergalactic-medium thermal history. The paper has several strengths: a broad simulation suite covering different thermal histories, WDM masses, and resolutions; the use of publicly available simulation and analysis codes; an honest statement of the phenomenological nature of the model; and a genuinely predictive test in Sec. 4.3, where a WDM model was designed analytically and only then run. The main weakness is that the validation is largely self-referential: the same functional form whose validity is at issue is used to extract the cutoff scale, and no fit uncertainties or full-shape residuals are reported. The significance of the claimed 5%, 20%, and 15% agreements therefore depends on the additional independent tests suggested below.","major_comments":[{"comment":"The headline agreements (5% for WDM, 20% for pressure, 15% for the sum) are reached by fitting the simulated real-space FPS with Eq. (2.7), the very ansatz whose validity is being tested, and then comparing the single fitted parameter kcut with theory values from Eqs. (A.3), (B.1), and (2.5). Because Eq. (2.7) has four free parameters (A, a, kcut, xi) and the fit range is chosen self-consistently as kmax = kcut (App. C), a functionally incorrect model could still return a kcut that tracks the true filtering scale while the full predicted FPS is biased. No error bars on kcut or goodness-of-fit statistics for the FPS fits are reported, so the precision of the comparison is not established. I recommend an independent full-shape test: fix kcut from the theoretical expressions, fix or fit A and a from the 3D power spectrum, and compare the complete predicted Delta1D(k) from Eq. (2.7) against the simulated FPS for k <= kcut, reporting residuals as a function of k and redshift.","section":"Sec. 4.2, Eq. (2.7), Figs. 5-7"},{"comment":"The observational analysis abandons the mapping defined by Eq. (2.7) and instead uses Eq. (E.1) to define kcut for both observations and simulations. No test is presented that the kcut from Eq. (E.1) equals the kcut that enters the theoretical combination rule Eq. (2.5) and the simulation validation. The Sec. 4.3 prediction chain (PuchweinLate thermal history, m = 3 keV) uses Eqs. (A.3), (B.1), and (2.5), i.e., the real-space cutoff, while the comparison in Figs. 8-10 uses kcut values from Eq. (E.1). If the two definitions differ by a redshift-dependent factor, the predicted degeneracy between the CDM and WDM simulations could be an artifact of the fitting forms. The authors should apply both fitting forms to simulations that include thermal broadening and redshift-space distortions and demonstrate the relation between the two kcut definitions.","section":"Sec. 4.3 and App. E, Eqs. (2.7) vs (E.1)"},{"comment":"The pressure-only comparison at z >~ 5.5 relies on a resolution-correction model, Eqs. (D.1)-(D.2), with a free parameter beta obtained from a linear fit across resolutions. The corrected filtering lengths are then compared with the GH98 theory without propagating the uncertainty of the extrapolation. Since at these redshifts the uncorrected filtering length is comparable to the resolution length (lF ~ 10-20 ckpc), the claimed 20% agreement is not yet robust. The authors should report uncertainties on the corrected lF and, where possible, validate the correction with a high-resolution run at the redshifts where the correction is largest.","section":"Sec. 4.2 and App. D, Eq. (D.2), Fig. 5 (left)"}],"minor_comments":[{"comment":"The dependence of the extracted kcut on the chosen mean flux normalization <F> = 0.5 is not quantified; a short test at, e.g., <F> = 0.3 and 0.7 would indicate the associated systematic.","section":"Sec. 3.1"},{"comment":"At z = 15 and z = 10 the 3D matter power spectrum no longer has a clean cutoff, yet kcut values from Eq. (2.6) are quoted; the text should state explicitly that these are biased by nonlinear power and are not used for the validation.","section":"Sec. 4.1, Fig. 4"},{"comment":"The algorithmic definition of i0 is ambiguous: 'the smallest i, when kmax,i = kN-i < kcut,i' does not specify whether kcut,i is from the same fit or the preceding one; a pseudocode listing or a worked example would remove the ambiguity.","section":"App. C"},{"comment":"The fitting pipeline is not released and the fitting function (2.7) has no documented uncertainty treatment; stating the availability of the code or providing a detailed numerical recipe would strengthen the paper.","section":"Reproducibility"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of JCAP and the authors are transparent about the phenomenological nature of the model. My main reservation is that the validation is self-referential and that the observational analysis changes the definition of the extracted cutoff without establishing equivalence to the theoretical quantity. These issues are fixable with additional analysis, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here is my read. The paper gives a practical recipe: treat the 1D flux power spectrum as the projection of a power-law 3D matter spectrum with a combined Gaussian cutoff, 1/k_cut^2 = 1/k_WDM^2 + 1/k_F^2, plus a constant small-scale term. What is new is not any single ingredient—those come from Lumsden, Viel, Gnedin & Hui—but the packaging and the validation. They run a broad suite of simulations covering different thermal histories, WDM masses, and resolutions, and show the fitted cutoff tracks theoretical values within 5% for WDM, 20% for pressure, and 15% when combined. The predictive test in Sec. 4.3 is the strongest part: they design a WDM plus late-reionization model that should mimic a CDM plus early-reionization model, run it, and find the flux power spectra agree within about 5%. That is a real result and the most convincing evidence the recipe works.\n\nThe soft spots are real but not fatal. The cutoff is extracted by fitting the same exponential form whose validity is at issue, so the reported agreements measure internal consistency of the extraction rather than an independent test of the model's shape. They do show the fits describe simulated spectra well over the fitted range, which is encouraging, but the range is limited to k below k_cut, so the shape above the cutoff is never tested. The observational extension switches to a different fitting form, Eq. (E.1), and no test establishes that k_cut from that form equals the k_cut appearing in Eq. (2.7). The paper openly concedes the lack of first-principles justification, which I appreciate, but that leaves the recipe on empirical ground only. There are also no error bars on the fitted cutoffs and no released code or data, limiting how far an independent reader can check the analysis.\n\nOverall, the paper is honest about its limits and the recipe is likely useful. It deserves a serious referee. I would ask for robustness checks: sensitivity of the fits to the chosen kmax, a direct comparison of the two fitting forms, error bars on the extracted cutoffs, and ideally a release of the fitting code. The central phenomenological claim seems supported as a practical mapping, even if the why is not yet understood.","headline":"A useful empirical recipe for converting Ly-alpha flux power spectrum cutoffs into WDM and pressure scales, but the validation is partly self-referential and the observational k_cut uses a different fitting form.","tokens_in":17529,"tokens_out":2525,"would_cite":true,"duration_ms":25499,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":null,"created_at":"2026-08-07T15:37:52.849020+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":null,"supporting_citations":[],"review_version":1}