{"id":"ad104abf-18fb-4ae7-a36e-f1bbe60332bc","arxiv_id":"2504.13705","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"Using 116 CHIME baseband non-repeaters, the energy function is Schechter-like at low redshift, but the high-redshift slope and rate evolution are too uncertain to choose between star-formation and old-stellar progenitors.","lead":"The authors measure how many non-repeating fast radio bursts exist at different energies and redshifts using 116 CHIME telescope bursts with precise baseband data. The rates are too uncertain to say whether FRB sources follow star formation or old stellar populations.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The Vmax calculation depends on an unspecified Flimit and on Hashimoto et al. weights rather than the newly fitted selection functions, so the corrected densities and inferred gamma slopes are not reproducible.","rationale":"The paper's goal is straightforward: to use improved baseband fluence measurements and the Vmax method to derive non-repeater energy functions and their redshift evolution. The low-redshift result is a plausible extension of earlier work, and the authors are appropriately cautious about small samples at high redshift. The load-bearing problem is not the existence of selection effects but the reproducibility of the corrected number densities. The reader already identified the unspecified Flimit in Eq. 15 and the dependence on Hashimoto et al. assumed distributions; I agree those are serious. I add a sharper observation: the selection functions that the paper fits in Eqs. 5-6 appear to be disconnected from the weight functions actually used in Eq. 17, because Eq. 17 explicitly cites Hashimoto et al. (2022) for wi(DM) and wi(F). If the authors intended their own selection functions, the omission is a straightforward reproducibility gap; if they used the Hashimoto weights literally, the new selection-function analysis does not enter the final result at all. Either way, the high-z gamma about -2 in the abstract rests on an undocumented step. A two-value rerun of the pipeline is a cheap, decisive test: if the fitted gamma remains stable under Flimit choices and under the two weight conventions, the central claim survives; if it shifts, the manuscript needs a stated Flimit, a clear definition of wi, a table of free-fit high-z slopes, and a discussion of the post-hoc point exclusion. I therefore keep the reader's conditional verdict rather than moving to accept or reject; the paper is not yet reproducible, but the core idea is testable and the gaps are fixable.","tokens_in":10318,"tokens_out":7415,"duration_ms":80222,"concrete_test":"Recompute all Vmax values and the three redshift-bin energy functions with Flimit = 0.5 Jy ms and with Flimit = 5 Jy ms, using Eqs. 14-18 and taking the selection functions from Eqs. 5-6 as wi = 1/si; if the best-fit gamma for redshift bin 3 shifts by more than its 1-sigma MC uncertainty relative to the curve plotted in Fig. 5, the reported high-redshift slope is not robust to the unspecified detection limit.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the low-redshift energy function has gamma = -1.66 and high-redshift regions have gamma about -2 rests entirely on the corrected number densities in Eq. 18, and those densities inherit two unstated or inconsistent inputs. First, Vmax in Eq. 15 contains Flimit, the CHIME fluence detection limit, but no numerical value is supplied anywhere in the paper. Because dL,max scales as (E/Flimit)^1/2, every source's Vmax, every rho value, and every fitted Schechter parameter depends on this unreported constant. The effect is not a global normalization: low-energy points, which dominate the high-redshift slope, are the most sensitive to Flimit, so a different choice can change the shape of the energy function. Second, although Section 2 derives updated selection functions s(DM) and s(Fnu) in Eqs. 5-6, the scaling used in Eq. 17 is explicitly written with weights wi(DM) and wi(F) derived from T. Hashimoto et al. (2022). If the published equations are literal, the new selection functions are never applied to the density correction, and the results inherit Hashimoto et al.'s assumed intrinsic DM lognormal and fluence power-law. If instead the authors intended w = 1/s with their own Eqs. 5-6, the text does not say so. Either way the corrected densities are not reproducible from the manuscript. The high-redshift gamma about -2 is additionally presented without tabulated free-fit parameters and after post-hoc exclusion of the lowest-energy point in redshift bins 2 and 3 (Section 3.1). The caveats in Section 4 do not mention these issues.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses 116 CHIME/FRB baseband non-repeaters with updated fluence measurements, estimates pseudo-redshifts from dispersion measures via Bayesian analysis, and applies the V_max method to derive the energy function in three redshift bins (0.04-0.44, 0.44-0.94, 0.94-2.50). The authors fit Schechter functions, claiming a low-redshift slope gamma = -1.66 (+0.33, -0.23), a steeper high-redshift slope gamma ~ -2, and an ambiguous volumetric-rate evolution that does not clearly track star-formation or stellar-mass density. Selection functions for DM and fluence are constructed by assuming intrinsic distributions from Hashimoto et al. (2022), and volumetric rates are compared with SFRD and SMD after arbitrary scaling.","tokens_in":10778,"tokens_out":3785,"duration_ms":34951,"significance":"If the central claims held, the paper would provide updated energy-function parameters based on the better-calibrated CHIME baseband fluences, which would be a useful contribution to FRB population studies. The manuscript explicitly attempts to account for selection effects, uses Monte Carlo and nested-sampling error estimates, and compares with SFRD/SMD evolution. However, the load-bearing results currently rest on an unreported detection limit, a weighting scheme that appears to bypass the newly fitted selection functions, and post-hoc exclusions that are not tabulated or justified. These issues mean the claimed energy-function slopes and rate evolution are not yet established from this analysis.","major_comments":[{"comment":"The value of Flimit is never specified anywhere in the manuscript. Since d_L,max is obtained from (1+z)E / (4 pi d_L,max^2 Delta nu) >= Flimit, every V_max in Eq. (14), every number density rho_i in Eq. (16), and hence every fitted Schechter parameter in Table 1 depends on this unreported constant. This is not a harmless overall normalization: for a fixed energy, the V_max of low-energy sources is most sensitive to Flimit, so changing Flimit can alter the shape of the energy function, not just its amplitude. Please state the adopted Flimit in physical units and show that the reported gamma values are robust to plausible variations of this threshold.","section":"Section 2.2, Eq. (15)"},{"comment":"The corrected number density is written explicitly with wi(DM) and wi(F) 'derived from T. Hashimoto et al. (2022)', while Section 2 states that the newly fitted s(DM) and s(Fnu) in Eqs. (5)-(6) are 'utilized for correcting the FRB number densities'. These statements are contradictory. If Eq. (17) is literal, the new selection functions are never applied, and the results inherit Hashimoto et al.'s assumed intrinsic DM log-normal and fluence power-law distributions; if the intended relation is w_i = 1/s with the authors' own Eqs. (5)-(6), the manuscript does not say so. Either way, the corrected densities, and therefore the energy-function slopes, are not reproducible from the text as written. Please make the weighting scheme explicit and rerun the analysis with the selection functions derived in this paper.","section":"Section 2.3, Eq. (17)"},{"comment":"The abstract's claim of a high-redshift slope gamma ~ -2 is not supported by the reported fits: Table 1 fixes gamma = -1.66 for bins 2 and 3, and the 'unfixed' fits shown in Figure 5 are described only qualitatively, with the lowest-energy points in both bins discarded post hoc. No best-fit values, uncertainties, or goodness-of-fit statistics are given for the unfixed fits. Because those low-energy points dominate the high-redshift slope, the exclusion rule must be stated in advance (e.g., a completeness or V_max criterion) and the free-fit parameters must be tabulated, including a test of whether gamma ~ -2 persists when the excluded points are retained or when Flimit is varied.","section":"Section 3.1, Table 1 and Figure 5"},{"comment":"The abstract states that the V_max method 'allows us to measure redshift evolution without prior assumptions', but the analysis assumes that the intrinsic DM and fluence distributions of the baseband sample are the same as those adopted for CHIME/FRB Catalog 1 by Hashimoto et al. (2022), as encoded in Eqs. (2)-(3). These assumed distributions enter the selection functions and hence the corrected number densities, so the energy-function normalization and slopes are conditional on those priors. Please qualify this claim and demonstrate robustness of the results to alternative choices of the intrinsic DM parameters (mu0, sigma) and fluence index alpha, for example by rerunning the analysis over a grid of these parameters.","section":"Section 2, Eqs. (1)-(3), and Abstract"}],"minor_comments":[{"comment":"There are typographical and grammatical errors, including 'Energy F unction' in the title and 'Analyse of the energy function' in the abstract; these should be corrected before resubmission.","section":"Title and Abstract"},{"comment":"Equation (3) as typeset, 'P(Fnu) ∝ −α (Fnu/Fnu,0)^{α−1}', appears to have a typographical issue: a power-law probability distribution is usually written P(F) ∝ F^{α−1} (or with a positive coefficient for a decreasing distribution), and the minus sign and placement of alpha are confusing. Please correct the expression.","section":"Section 2, Eq. (3)"},{"comment":"The exposure time is given as t_obs = 0.59 yr, attributed to the CHIME/FRB Catalog 1 observation, but the baseband data span 2018 December 9 to 2019 July 1. Please clarify whether this same exposure is appropriate for the baseband subset and whether the subset's live-time fraction is accounted for.","section":"Section 2.3, Eq. (16)"},{"comment":"The text says the energy functions are poorly constrained in bins 2 and 3 when the slope is fixed, but the main conclusion about the high-redshift slope is then drawn from the unfixed fits. Please make the caveat explicit in the abstract and conclusions, since the high-redshift result is the most uncertain part of the analysis.","section":"Section 3.1"},{"comment":"The manuscript does not state whether the derived energy-function data points, the fitted selection functions, or the analysis code will be made available. Given the reproducibility issues in Eqs. (15) and (17), a data/code availability statement would be valuable.","section":"Data availability"}],"recommendation":"major_revision","confidential_remarks":"The paper's methodology is closely modeled on Hashimoto et al. (2022) but uses baseband fluences; this is an incremental but potentially useful contribution if the analysis can be made reproducible. The two most serious technical problems are the unspecified Flimit in Eq. (15) and the apparent mismatch between the newly derived selection functions and the weights in Eq. (17); both directly affect the paper's central quantitative claims and must be fixed before the manuscript can be considered for publication. I would also encourage the editor to ask for the unfixed high-redshift fit parameters to be tabulated, because the abstract's gamma ~ -2 claim currently rests on fits that are not reported."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Punchline: This is a legitimate update of the Hashimoto et al. Vmax analysis to the CHIME baseband catalog, but the headline high-redshift slope (γ≈-2) is not reproducible from the text, and the low-redshift fit inherits an unspecified detection limit. The paper deserves a serious referee, but needs major revision.\n\nWhat's actually new: they use the 116 non-repeaters with baseband fluences and better localization from CHIME/FRB Collaboration 2024, which is a real improvement over Catalog 1. Applying the Vmax method to this sample is a reasonable step, and the low-redshift Schechter fit with γ=-1.66 is consistent with earlier work. They also re-fit selection functions for DM and fluence (Eqs 5-6) rather than just importing them, and they are honest about the ambiguity in the rate evolution.\n\nSoft spots, in order of importance. First, the CHIME detection limit Flimit in Eq. 15 is never given a number. Because Vmax is ∝ (E/Flimit)^{3/2} at fixed distance, every corrected density depends on this constant. It is not just a global normalization: low-energy sources, which drive the high-z slope, are the most sensitive, so different choices can change the shape of the energy function. Second, Eq. 17 says the weights wi(DM) and wi(F) are taken from Hashimoto et al. 2022, not from the new selection functions derived in Section 2. If that is literal, the new s(DM) and s(Fν) are never used. If the authors intended w=1/s, they do not say so. Either way, the corrected densities are not reproducible. Third, the high-z slope is presented after first fixing γ=-1.66, then unfixing after discarding the lowest-energy points in bins 2 and 3, but the unfixed fits are not tabulated. The abstract's \"γ≈-2\" is therefore not a well-defined measurement. Fourth, the selection functions themselves are built by assuming the intrinsic DM and fluence distributions are the same as in Hashimoto et al., which puts a prior on the very quantity being measured. That is a known limitation but should be stated plainly. Finally, only statistical errors are propagated; systematics from the DM host distribution, the pseudo-redshift prior, and Flimit are not included.\n\nThese are not fatal flaws. The low-redshift result is plausible and the dataset is worth examining for FRB population studies. But as it stands, the paper's central high-redshift claim is not supported by the written methods. I would send it to peer review, but insist that the authors specify Flimit, clarify the weighting scheme, tabulate the unfixed fits, and add systematic error estimates.","headline":"A worthwhile Vmax update to CHIME baseband data, but the high-redshift slope is not reproducible until Flimit and the actual weighting scheme are specified.","tokens_in":11243,"tokens_out":3995,"would_cite":false,"duration_ms":32272,"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":"Baseband CHIME data put the low-redshift energy-function slope of non-repeaters at -1.66.","keywords":["fast radio bursts","non-repeaters","energy function","volumetric rate","Vmax method","dispersion measure","redshift evolution","baseband data"],"falsifier":"Recompute the $V_\\mathrm{max}$ densities with intrinsic DM and fluence distributions measured directly from a complete, injection-calibrated baseband sample, or from localized events with known host DM; if the low-redshift slope then lies outside $\\gamma=-1.66^{+0.33}_{-0.23}$, or the high-redshift slope no longer settles near $-2$, the paper's central claim is falsified.","tokens_in":10144,"feed_emoji":"📡","tokens_out":8464,"duration_ms":68570,"temperature":0.7,"pith_summary":"This paper aims to establish the energy function and redshift-dependent volumetric rate of non-repeating fast radio bursts from the 116 CHIME/FRB baseband detections that survive its cuts, using the $V_\\mathrm{max}$ method without assuming how the rate evolves with redshift. It finds a Schechter-like low-redshift energy function in the bin $0.04