{"id":"594e1f13-3fd5-4857-b631-a6879771bbe2","arxiv_id":"2504.17478","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A modeled five-year CTAO campaign is projected to rule out local primordial black hole burst rates above about 36 per cubic parsec per year, roughly ten times better than LHAASO.","lead":"This paper calculates how well the future CTAO gamma-ray observatory could detect the final gamma-ray flash of evaporating primordial black holes. If the projection holds, a five-year campaign could set the strongest limit yet on how often these bursts occur near Earth, roughly ten times better than current detectors.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The projected 5-year CTAO limit treats S in Eqs. (2.8) and (2.14) as 5 calendar years; using realistic IACT livetime (~0.5–1 yr) weakens the headline 36 pc⁻³ yr⁻¹ to ~180–360 and puts the claimed order-of-magnitude improvement over LHAASO at risk.","rationale":"I reviewed the paper in good faith. The central claim is a prospective 99% CL upper limit on the local PBH burst rate density from a 5-year CTAO campaign, with the South array at τ = 100 s giving 36.26 pc⁻³ yr⁻¹ and an order-of-magnitude improvement over LHAASO. The formalism is standard and internally arithmetically consistent: the Table 1 values follow from Eqs. (2.4)–(2.14) given the adopted effective area, background rates, FOV, and S. My stress-test focused on the least secure condition supporting the headline number: the meaning of S. The manuscript neither defines S nor applies the duty-cycle correction that IACT observations require. Since the limit is inversely proportional to S, treating a 5-year campaign as S = 5 yr rather than as ~0.5–1 yr of livetime overestimates sensitivity by roughly an order of magnitude. This directly threatens the core comparative claim against LHAASO and HAWC. The paper's own treatment of τ = 10⁵ s bursts shows awareness that CTAO cannot observe continuously, making the omission for the main S more surprising. I also noted a secondary inconsistency between the written FOV formula and the numerical volume used in Table 1; while the table appears to use a realistic small FOV, the text's Eq. (2.12) describes a full-sky annulus, so the manuscript must be corrected to avoid a factor-of-10³ discrepancy in the search volume. These issues are fixable, but they are load-bearing because they change the headline number by an order of magnitude and may invert the comparison with LHAASO. The reader's weakest assumption identifies exactly the S/livetime ambiguity; I agree. I recommend no change to the reader's CONDITIONAL verdict: the paper is a useful prospective estimate with a plausible method, but it requires clarification of S and correction/reconciliation of the FOV equation before the quoted limit and improvement claim can be accepted.","tokens_in":136,"tokens_out":16184,"duration_ms":259086,"concrete_test":"Recompute Table 1 with S equal to the actual CTAO livetime: adopt a representative dark-time duty cycle (e.g., ~1500 h yr⁻¹ for CTAO South, about 0.17 yr per calendar year) over 5 calendar years, and compare the resulting τ = 100 s South limit with the LHAASO limit from Ref. [31] at the same τ. If the quotient drops below ~10, the abstract's 'improvement of one order of magnitude' fails, confirming that S was treated as calendar time. As a first step, ask the authors to state explicitly whether S is livetime or calendar time and to provide the livetime assumed in Eqs. (2.8) and (2.14).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper never defines whether the search duration S in Eqs. (2.8) and (2.14) is calendar time or telescope livetime. The abstract and Table 1 use '3-yr' and '5-yr' campaign labels, and the Table values scale linearly with S, so the headline 36.26 pc⁻³ yr⁻¹ for the South array at τ = 100 s is directly proportional to S = 5 yr. CTAO, as a pointed IACT, can only accumulate exposure during dark, clear nights; the typical usable dark-time duty cycle is ~10–20%, giving roughly 0.5–1.0 yr of livetime in a 5-year campaign. If the correct S is 0.5 yr, the limit becomes ~363 pc⁻³ yr⁻¹; with S = 1 yr, ~181 pc⁻³ yr⁻¹. This is an order of magnitude worse and removes the claimed 'one to two orders of magnitude' advantage over LHAASO, whose wide-field duty cycle is nearly continuous. The paper's own §4 acknowledges that 'telescopes cannot observe continuously over such long timescales' when correcting τ = 10⁵ s bursts, but no analogous duty-cycle factor is applied to S. A secondary internal inconsistency: Eq. (2.12) defines FOV(θᵢ) = 2π(cosθᵢ,min − cosθᵢ,max), a full-sky zenith annulus, not the ~4–8° instantaneous CTAO field of view; the numerical Table 1 is consistent with a ~7.5° single-pointing FOV, so either the written formula or the implementation is wrong, and this must be reconciled before the volume calculation can be trusted.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper projects the sensitivity of the Cherenkov Telescope Array Observatory (CTAO) to the local burst rate density of evaporating primordial black holes (PBHs). Using the standard evaporation model, a Poisson detection criterion with trial-factor correction, and public CTAO instrument response functions, the authors compute the maximum detectable distance and the 99% CL upper limit on the local PBH burst rate density for the CTAO South and North arrays for remaining lifetimes from 10^-3 s to 100 s (Table 1), with additional approximate corrections for longer lifetimes up to 10^5 s. They find a headline 5-year limit for the South array of rho_dot < 36.26 pc^-3 yr^-1 at tau = 100 s, claim an improvement of one to two orders of magnitude over LHAASO and HAWC, and propose an observational strategy favoring long (>10^3 s) exposures of the same sky region.","tokens_in":12156,"tokens_out":12345,"duration_ms":110863,"significance":"If the calculation is correct, this is a useful and falsifiable projection for CTAO: a five-year campaign with no detected bursts would exclude rho_dot > 36.26 pc^-3 yr^-1 at 99% confidence, assuming the quoted effective areas, background rates, and exposure duration are accurate. The paper's main strengths are its transparency (the derivation from the source spectrum to the upper limit is spelled out in Section 2), its use of public CTAO performance data, and the absence of any parameter fitted to the target sensitivity; the trial-factor and Poisson formulæ follow the established literature. However, the two ambiguities discussed below (the definition of S and the field-of-view formula) must be resolved before the headline number can be taken at face value. The projected improvement over existing wide-field instruments is interesting but currently rests on an unverified duty-cycle assumption.","major_comments":[{"comment":"The written definition of the field of view as FOV(theta_i) = 2*pi*(cos theta_i,min - cos theta_i,max) describes a full-sky zenith annulus, appropriate for an all-sky monitor, not the ~4-8 degree instantaneous field of view of a pointed IACT. The numerical results in Table 1, however, are consistent with a single pointing cone of ~7.5 degree diameter; for the South tau=100 s row, the volume implied by r_max=1.782 pc, S=5 yr, and UL=36.26 pc^-3 yr^-1 gives Omega_eff ~ 0.013 sr. The paper must reconcile this discrepancy by replacing Eq. (2.12) with the actual pointing-cone solid angle (or by defining the zenith-band binning and survey strategy explicitly); as written, the formula would overestimate the search volume by orders of magnitude and the derivation is not reproducible as stated.","section":"Sec. 2.2, Eq. (2.12) and Table 1"},{"comment":"The paper never states whether the search duration S is calendar time or effective livetime. Because CTAO is a pointed IACT, a five-year calendar campaign yields only roughly 0.5-1.0 yr of on-sky livetime (a duty cycle of 10-20%). Since UL_99 = 4.6/(V S), the headline 36.26 pc^-3 yr^-1 scales linearly with S; with S=0.5 yr it becomes ~360 pc^-3 yr^-1, which would remove the claimed one-to-two-order-of-magnitude improvement over LHAASO. The paper's own Section 4 applies a nighttime-visibility correction for tau=10^5 s bursts, but no analogous duty-cycle factor is applied to S; this internal inconsistency must be addressed, and the reported limits should be labeled as calendar-time or livetime projections.","section":"Sec. 2.2, Eqs. (2.8) and (2.14)"},{"comment":"The claim that tau=100 s yields the best constraint is not supported by the numbers shown. Section 4 states that without the approximate correction, the strongest constraint would come from tau=10^3 s, yet the corrected limits for tau=10^3-10^5 s are not tabulated, so the reader cannot verify that 36.26 pc^-3 yr^-1 is the minimum after all corrections. The corrected values should be listed (or plotted with numerical labels) and the conclusion should identify the exact lifetime, integration-window, and correction assumptions that produce the quoted best limit.","section":"Sec. 4 and Table 1 / Conclusions"}],"minor_comments":[{"comment":"The sentence 'Section 4 outlines the expected sensitivities of CTAO' is incorrect: the sensitivities are presented in Section 3, while Section 4 is the Discussion.","section":"Sec. 1 (last paragraph)"},{"comment":"The abstract states an improvement of 'one order of magnitude' over LHAASO, while the Conclusions claim 'one to two orders of magnitude'; these statements should be harmonized after the duty-cycle question is settled.","section":"Abstract and Conclusions"},{"comment":"The phrase '50-hour photon effective area of CTAO after gamma/hadron separation' is ambiguous: the 50-hour point-source sensitivity is not an effective area per se, so the authors should clarify whether they used the point-source effective-area curve or the differential sensitivity, and how the 50-hour normalization relates to the assumed S.","section":"Sec. 3"},{"comment":"The '50% correction' for tau=10^3 s and 10^4 s is not clearly derived: if the integration window is Delta t = 10^2 s, the fraction of photons collected from a tau=10^3 s burst should be about 10%, not 50%; the correction scheme should be defined in terms of the search-window function.","section":"Sec. 4"},{"comment":"The trial-factor expression uses (theta_fov/theta_res)^2 but does not specify whether theta_res is the point-spread-function resolution or the angular bin size; this should be stated for reproducibility.","section":"Eq. (2.8)"},{"comment":"The units of S are described as '3-yr' and '5-yr' but the equations use seconds; the text should state the unit conversion explicitly (e.g., 5 yr = 1.58 x 10^8 s).","section":"Table 1 and text"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a timely topic and the analytic framework is a reasonable starting point, but the two ambiguities in the core calculation (definition of S and the field-of-view formula) must be fixed before the headline number can be endorsed. If the authors clarify that S is livetime and correct Eq. (2.12), the revised manuscript may well be suitable for publication; if the duty-cycle effect substantially weakens the claimed improvement over LHAASO, the paper's central claim will need to be reframed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper computes CTAO's prospective sensitivity to local PBH burst rate density using the standard Ukwatta-type formalism and the public prod5 response functions. That is genuinely new: previous CTAO statements were either qualitative or based on MAGIC extrapolation. The specific South/North limits and the 10^3 s staring recommendation are useful additions. The arithmetic in the simple parts checks out, and the headline 36.26 pc^-3 yr^-1 follows from the stated equations if you take S=5 yr as the exposure.\n\nThe problems are mostly definitional. First, S in Eqs. (2.8) and (2.14) is never defined as calendar time or livetime. CTAO is a pointed IACT; a five-year calendar campaign yields maybe 0.5–1 yr of on-sky time. The headline limit scales linearly with S, so using livetime moves 36.26 to roughly 180–360 pc^-3 yr^-1, which erodes the claimed advantage over LHAASO. The paper's own §4 acknowledges that 'telescopes cannot observe continuously over such long timescales' when applying corrections for tau=10^5 s, but no analogous duty cycle is applied to S.\n\nSecond, the written FOV formula (2.12) gives a zenith annulus (2π(cosθmin−cosθmax)), but the numbers in Table 1 are consistent with a single circular FOV of about 7.5° diameter. Either the equation or the implementation is wrong; this needs reconciliation before the volume calculation is trustworthy.\n\nThird, the corrections for long-duration bursts are ad hoc (50% fudge factors for tau≥10^3 s). That affects the proposed observing strategy, not the headline 100 s limit, but it should be stated as such.\n\nThe paper provides no code or data, but the method is clear enough to reimplement. The citation pattern is fine; the only self-citation is the LHAASO comparison, not the CTAO derivation.\n\nThe core idea is sensible and the topic is of interest to the gamma-ray and PBH communities. The paper deserves a serious referee, but it needs a major revision to define S explicitly, fix the FOV formula, and propagate a realistic duty cycle. As written, the headline claim is probably optimistic by an order of magnitude.\n\nMy take: send to peer review, but make sure the referee checks the exposure definition. I would not cite this in its current form.\n\nBest","headline":"Useful CTAO projection undermined by an unstated definition of the search duration and a FOV formula that contradicts the numbers.","tokens_in":12809,"tokens_out":12232,"would_cite":false,"duration_ms":113236,"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":"Five years of no bursts would cap CTAO's primordial black hole burst-rate limit at 36.26 per cubic parsec per year.","keywords":["primordial black holes","Hawking radiation","gamma-ray bursts","CTAO","local burst rate density","very-high-energy gamma rays","IACT","dark matter"],"falsifier":"Recompute the 99% upper limit from Eq. (2.14) using $S$ equal to the realistic on-sky livetime of a five-year CTAO campaign rather than five calendar years, keeping everything else fixed; if the resulting limit exceeds the LHAASO value the paper quotes, the claimed one-order-of-magnitude improvement is not realized.","tokens_in":11548,"feed_emoji":"🔭","tokens_out":7780,"duration_ms":69053,"temperature":0.7,"pith_summary":"This paper calculates how well the Cherenkov Telescope Array Observatory could detect the final gamma-ray burst of an evaporating primordial black hole, and what a null result would mean. The authors argue that a five-year CTAO campaign with no detected bursts would rule out local PBH burst rate densities above about 36.26 bursts per cubic parsec per year at 99% confidence, using the South array for bursts with a remaining lifetime of 100 seconds. That limit is one to two orders of magnitude tighter than the LHAASO limit and the current HAWC constraints, and several times tighter than the projected SWGO sensitivity. The paper also identifies bursts with remaining lifetimes around 1000 seconds as CTAO's sweet spot and recommends dwelling on the same sky region for at least 1000 seconds. The result matters because a detection would be direct evidence that primordial black holes exist, while a tighter limit constrains the small-scale primordial density fluctuations that could produce them.","feed_headline":"CTAO would beat LHAASO 10x on primordial black hole burst limits","feed_subtitle":"No bursts in five years would cap the local primordial-black-hole burst rate at 36.26 per cubic parsec per year.","key_machinery":"The load-bearing identity is the Poisson no-detection formula $\\mathrm{UL}_{99} = 4.6/(V S)$: with no bursts seen, the 99% upper limit on the local burst rate density is 4.6 divided by the product of the detectable volume $V$ and the search duration $S$. The volume itself is built from Eq. (2.10), where the maximum detectable distance $r_{\\max}$ follows from comparing the expected burst photon count $\\mu(r,\\theta_i,\\tau)$ against the background through a detection threshold corrected for $N_t = (S/\\tau)(\\theta_{\\rm fov}/\\theta_{\\rm res})^2$ trials. These two formulas turn the detector's effective area, background rate, field of view, angular resolution, and campaign length directly into a rate-density limit, which is why the result is sensitive to the assumed exposure time.","core_discovery":"The paper's central claim is that CTAO's sensitivity to primordial black hole evaporation bursts is set by a volume-time trade-off: its field of view is small, so the maximum detectable distance is tiny, about 1.8 pc for the 100-second bursts, but the upper limit on the local burst rate density scales inversely with that volume times the five-year search duration. Combining the standard Hawking spectrum for the burst, CTAO's effective area and background rates, a Poisson trial correction for blind searches, and the requirement that a burst be bright enough to pass a 5-$\\sigma$ threshold, the authors obtain a 99% confidence upper limit of $\\dot{\\rho} < 36.26\\,\\mathrm{pc}^{-3}\\,\\mathrm{yr}^{-1}$ for the South array at $\\tau = 100$ s. They further show that for bursts shorter than about 0.01 s, LHAASO remains more sensitive, while for longer bursts CTAO wins by one to two orders of magnitude; the North array is weaker by roughly an order of magnitude. Based on a correction for bursts longer than the integration window, the paper concludes that the optimal CTAO strategy is to integrate on one sky region for more than 1000 seconds.","pith_inferences":["A direct extension the paper leaves implicit: if the five-year campaign is expressed in usable on-sky time rather than calendar time, the 36.26 figure scales linearly upward; for a pointed telescope observing only on dark, clear nights, the limit would weaken substantially, potentially eroding the claimed advantage over LHAASO.","The same volume-time identity implies the biggest sensitivity lever is the maximum detectable distance, because it enters cubed; improving the low-energy threshold or effective area would therefore pay off faster than extending the campaign.","The trial-count formula treats time bins of length equal to the burst lifetime; using shorter or adaptive binning may reduce the trials penalty for short bursts and is a testable optimization the paper mentions only as future work."],"forward_implications":["If no bursts are seen, the 99% upper limits in Table 1 become direct experimental targets: 36.26 bursts per cubic parsec per year for the South array at 100 seconds, rising to 6755 for 1-millisecond bursts.","For bursts with remaining lifetimes longer than about 0.01 seconds, CTAO's limits beat LHAASO's by one to two orders of magnitude; for shorter bursts, LHAASO stays ahead.","The optimal search strategy is to point CTAO at a single sky region for at least 1000 seconds, since that is the burst duration at which CTAO's constraints are strongest.","Any detection or improved limit on the local burst rate density translates into a constraint on the amplitude of primordial density fluctuations on scales far smaller than those probed by the cosmic microwave background."],"supporting_citations":[{"why":"Supplies the Hawking radiation temperature-mass relation that makes black hole evaporation a physical process.","marker":"[11]"},{"why":"Provides the time-integrated photon spectrum from an evaporating black hole used to compute expected burst counts.","marker":"[36]"},{"why":"Establishes the blind-search trial correction and detection threshold logic that sets the minimum detectable burst flux.","marker":"[14]"},{"why":"Provides CTAO background particle rates and effective area inputs for the sensitivity calculation.","marker":"[37]"},{"why":"Provides CTAO's effective area and field-of-view parameters used to compute maximum detectable distance.","marker":"[38]"},{"why":"Introduces the upper-limit formula and the burst search framework adopted here.","marker":"[26]"},{"why":"Gives the HAWC local burst rate density constraint that the paper claims to improve by one to two orders of magnitude.","marker":"[29]"},{"why":"Gives the LHAASO prospective sensitivity that is the paper's main comparison baseline.","marker":"[31]"},{"why":"Provides the SWGO projected sensitivity that CTAO is claimed to beat by several times.","marker":"[30]"}],"fun_headline_variants":["For long PBH bursts, CTAO beats LHAASO by 10x","CTAO beats LHAASO 10x on PBH bursts, but only with long stares","CTAO could exclude PBH bursts at a rate 10x tighter than LHAASO","Staring at one sky spot for 1000s: CTAO's PBH burst edge"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes the five-year campaign duration is the exposure time $S$ that enters the volume-time product; for a pointed telescope that only records data on dark, clear nights, that assumption is the load-bearing one, and a shorter real livetime would weaken the quoted limit almost in proportion.","fun_headline_variants_meta":{"raw":{"variants":["For long PBH bursts, CTAO beats LHAASO by 10x","CTAO beats LHAASO 10x on PBH bursts, but only with long stares","CTAO could exclude PBH bursts at a rate 10x tighter than LHAASO","Staring at one sky spot for 1000s: CTAO's PBH burst edge"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001033,"raw_usage":{"total_tokens":4379,"prompt_tokens":1006,"completion_tokens":3373,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":622,"completion_tokens_details":{"reasoning_tokens":3273}},"tokens_in":622,"tokens_out":3373,"duration_ms":22190,"temperature":1.0,"reasoning_tokens":3273,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T10:40:28.220976+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the 99% upper limit from Eq. (2.14) using $S$ equal to the realistic on-sky livetime of a five-year CTAO campaign rather than five calendar years, keeping everything else fixed; if the resulting limit exceeds the LHAASO value the paper quotes, the claimed one-order-of-magnitude improvement is not realized.","supporting_citations":[{"cited_title":"Cta observatory performance","cited_arxiv_id":null,"evidence_quote":"Provides CTAO background particle rates and effective area inputs for the sensitivity calculation."},{"cited_title":"The technology behind the next generation very high-energy gamma-ray detector","cited_arxiv_id":null,"evidence_quote":"Provides CTAO's effective area and field-of-view parameters used to compute maximum detectable distance."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the SWGO projected sensitivity that CTAO is claimed to beat by several times."}],"review_version":1}