{"id":"12edd236-bc24-4679-ae3d-10fca181ab05","arxiv_id":"2501.09789","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"CTAO's Galactic Plane Survey could detect the brightest Milky Way dark matter sub-halo at 5-sigma for TeV-scale WIMPs annihilating to b-quarks with cross section about 3e-25 cm^3/s, roughly ten times the canonical thermal value.","lead":"This paper simulates CTAO's planned Galactic Plane Survey to estimate how well it could detect gamma rays from annihilating dark matter clumps inside the Milky Way. It finds that a single bright clump could be detected at 5-sigma for annihilation cross sections near 3e-25 cm^3/s, about ten times above the standard thermal value, with the best chances slightly off the Galactic plane.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 3e-25 claim is set by the single most extreme subhalo among 10^10 SL17 Monte Carlo skies, not by a typical realization; the median brightest subhalo has ~100x lower J, moving the threshold to ~1e-23, and the outlier's listed NFW parameters are internally inconsistent.","rationale":"The reader's conditional verdict is well aligned with the paper's transparency, but the most load-bearing point is sharper than 'SL17 may be wrong.' The paper's own Table I shows a factor ~100 gap in Jtot between the global-max subhalo used for the abstract's headline (9.56e21 GeV^2/cm^5) and the median of per-realization brightest subhalos (8.94e19 resilient, 1.59e19 fragile). Since TS is proportional to flux, which is proportional to J*, the median-realization 5-sigma threshold is about 1e-23 cm^3/s, which is the paper's own population-level number. Thus the 3e-25 claim is an extreme-order-statistic ceiling over 10^10 Monte Carlo skies, and the paper does not convert it into a detection probability for the actual Milky Way. This is a framing and correctness issue in the abstract, not a fatal flaw, because the population-level statement is present and the paper explicitly says 'brightest from among various realizations.' A second, more technical anomaly supports caution: the top-row subhalo's stated rho_s, r_s, r_Delta, and M_Delta are inconsistent with the NFW mass formula by a factor ~4, and the implied J-factor by ~16, so the template used for the headline curve should be re-verified. Both issues are addressable and do not overturn the overall conditional acceptance, but they should be fixed before the abstract's 3e-25 number is taken as the GPS reach.","tokens_in":34502,"tokens_out":15594,"duration_ms":163440,"concrete_test":"Re-run the single-subhalo likelihood analysis of Sec. IV A 1 replacing the global-max subhalo with the median per-realization brightest subhalo template (Table I, row 7: M_Delta = 7.7e7 Msun, d = 12.3 kpc, Jtot = 8.94e19, theta_68 = 0.7 deg) and with the 95th-percentile realization; if the required for 5-sigma rises from 3e-25 to ~1e-23 cm^3/s or above, the abstract should state that 3e-25 is an extreme-order-statistic ceiling, not the expected MW reach. As a cross-check within the same run, recompute Jtot for the top-row template in CLUMPY from the stated rho_s, r_s, r_Delta and d; if the returned Jtot deviates from 9.56e21 GeV^2/cm^5 by more than ~20%, the template used for the headline curve is not the one documented in the text.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The abstract's headline 5-sigma sensitivity, ~ 3e-25 cm^3/s, is driven by the top row of Table I: the global maximum J-factor (Jtot = 9.56e21 GeV^2/cm^5) over the 10^10 SL17 realizations. The two 'median brightest subhalo' templates the authors themselves adopt for the population study (Table I, bottom rows) have Jtot = 8.94e19 (resilient) and 1.59e19 (fragile). Because the required cross section scales as 1/J at fixed mass and spectrum, the per-realization median threshold is ~1e-23 cm^3/s for the resilient case and ~5e-23 for the fragile case, matching the population-level numbers in the abstract. The 3e-25 number is therefore a 1-in-10^10 extreme-order-statistic ceiling, not the reach of the GPS for a typical Milky Way; the paper does not quantify the probability that the actual MW contains such an outlier. This matters because the entire single-object forecast, including the morphological-discrimination and profile-reconstruction results in Fig. 6, is built on this one non-representative object. In addition, the quoted properties of that object are not internally consistent for an NFW profile: with rho_s = 3.4e6 Msun/kpc^3, r_s = 2.36 kpc, r_Delta = 0.162 kpc, the enclosed mass M( of order 3 x 10^-25 cm^3/s at TeV masses, while the population-level reach is 10^-23 to 10^-22 cm^3/s. It also assesses morphological discrimination from point-like and Gaussian sources, angular profile reconstruction via annuli, cumulative diffuse emission from unresolved subhalos, and compares with main-halo and dwarf-spheroidal limits.","tokens_in":1848,"tokens_out":4592,"duration_ms":190626,"significance":"If the single-object reach were robust, the GPS would probe TeV WIMP annihilation cross sections an order of magnitude above the thermal value through Galactic substructure, complementing dwarf and Galactic-center searches. The paper has real strengths: it uses publicly available tools (CLUMPY, gammapy, CTAO prod5 IRFs), treats two tidal scenarios explicitly, tests three levels of instrumental systematics, and transparently labels the single-object case as the most optimistic. The population-level result (10^-23 to 10^-22 cm^3/s) is consistent with the median brightest subhalo J-factors and is the more defensible forecast. The quantitative single-object claim, however, is not yet reliable as printed because the benchmark template's tabulated NFW parameters are internally inconsistent and because the chosen object is an extreme order statistic over 10^10 realizations whose occurrence probability is not quantified.","major_comments":[{"comment":"The benchmark 'overall brightest' subhalo parameters in the first row of Table I are internally inconsistent for a truncated NFW profile. With rho_s = 3.4 x 10^6 M_sun/kpc^3, r_s = 2.36 kpc, and r_Delta = 0.162 kpc, the enclosed mass is M( ~ 3 x 10^-25 cm^3/s is obtained from the single realization with the largest J-factor among 10^10 SL17 Monte Carlo skies. Table I shows that the median brightest subhalos in the same model have J_tot = 8.94 x 10^19 (resilient) and 1.59 x 10^19 (fragile) GeV^2/cm^5, i.e., factors of 107 and 600 below the benchmark. Because the required cross section scales as 1/J (Eq. 1), the corresponding thresholds are about 3 x 10^-23 and 5 x 10^-23 cm^3/s, consistent with the population numbers in Sec. V A. The paper does not report the distribution of the brightest-subhalo J-factor over realizations, so the probability that the actual Milky Way contains such an extreme object is not quantified. The abstract and conclusions should either present this probability or explicitly frame the 3 x 10^-25 value as a 'best of 10^10 skies' upper envelope rather than as the GPS reach for a typical Milky Way; as written, the headline is easily misread as a robust single-object forecast.","section":"Abstract / Sec. II B 2 / Sec. V A"},{"comment":"The population-level detection thresholds in Sec. V A are computed by applying the re-scaling formula F(theta_68) = A sqrt(alpha theta_68^2 + sigma_PSF^2) to every subhalo in the realizations. The parameter alpha is fitted to three benchmark templates (point-like and the two median subhalos), with a reported minimum error of order 15%. Because this formula is used to classify the full population, including objects whose theta_68 lies outside the calibrated range and whose NFW concentrations may differ substantially from the benchmarks, the authors should either validate Eq. (8) against direct likelihood simulations for a wider grid of theta_68 and profile shapes, or explicitly state that the 15% figure is only demonstrated at the three calibration points. Without this, the population numbers inherit an unquantified systematic from this empirical scaling.","section":"Sec. III F / Eq. (8)"}],"minor_comments":[{"comment":"The abstract's 3 x 10^-25 cm^3/s value should be accompanied by the caveat that it is the no-systematics sensitivity; Sec. IV A 1 and Fig. 4 show that 3% systematics worsen the reach by roughly a factor of two and 10% systematics by an order of magnitude.","section":"Abstract / Sec. IV A 1"},{"comment":"The sentence \"the minimal error we can achieve with this method is of the order of 15%, corresponding to a position-independent value of alpha = 0.06\" is ambiguous: it should state whether alpha is fitted per position and then averaged, or fitted globally, and how the 15% error is defined and distributed over theta_68.","section":"Sec. III F / Eq. (8)"},{"comment":"Figure 1 shows the distribution of theta_68 for the brightest subhalo per realization, but the detection thresholds are driven by J-factor; showing the corresponding distribution of the maximum J-factor (and its percentiles) would directly address the probability content of the headline claim.","section":"Fig. 1 / Sec. II B 2"},{"comment":"The main-halo sensitivity curve is derived from a simplified setup with only CR and interstellar emission backgrounds and with the region |l| < 12 deg excluded; the text notes this, but the figure caption should carry the same caveat so the comparison with subhalo curves is not overinterpreted.","section":"Sec. VI / Fig. 10"},{"comment":"There are several typographical errors, including 'W design' in Sec. III F (should be 'We design'), 'SL-17' versus 'SL17' inconsistencies, and 'F ermi' in reference [77]; these should be corrected in a final proofread.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Hi — quick take before you spend an evening on it. This is the first CTAO Galactic Plane Survey forecast that treats DM subhalos as extended sources, and the population-level result is solid. But the 3e-25 headline is the absolute maximum over 10^10 simulated skies, not a typical Milky Way, and the top row of Table I has an internal NFW inconsistency I'd want fixed before quoting the single-object reach.\n\nWhat's new: previous CTAO subhalo forecasts used point-like sources and focused on the extragalactic survey or all-sky legacy data. This paper folds in the actual GPS exposure pattern, accounts for extended morphologies, and runs both the SL17 fragile and resilient population variants. The analysis is careful: public IRFs, gammapy, CLUMPY, explicit interstellar emission treatment, and a reasonable systematic-error survey. The morphological-discrimination and profile-reconstruction parts are genuinely useful, and the comparison with pulsar halos is a nice addition.\n\nWhere it gets shaky: the stress-test note checks out. The abstract's 3e-25 is the most extreme subhalo found across 10^10 realizations. The median brightest subhalo has a J-factor roughly 100 times lower, putting the per-realization threshold near 1e-23 — consistent with the population numbers the paper itself quotes. The abstract does not say \"most optimistic realization\" and does not say \"without systematics\"; both qualifications are buried in the text. With 3% systematics the best reach degrades by about a factor of two.\n\nThe Table I inconsistency is concrete. For the top-row template, the NFW mass enclosed within r_Delta = 0.162 kpc is ~1.2e6 Msun using the listed rho_s and r_s, not the stated M_Delta = 4.8e6 Msun. Either the parameters or the J-factor don't match the profile described in Sec. IIB2. That template drives the 3e-25 curve and the Fig. 6 results, so it's load-bearing, though it may be a typo or a mix-up between pre- and post-stripping values. It needs a direct check.\n\nThe matched-model injection is fine for a forecast, and the SL17 dependence is acknowledged. I'd send it out. Ask for the Table I fix, a clear statistical statement of what \"brightest\" means, and release of the subhalo catalogs or code. The population-level section alone is citable and will be useful to the CTAO community.","headline":"First GPS forecast for extended DM subhalos, but the headline sensitivity is an outlier and Table I has an internal NFW inconsistency.","tokens_in":35494,"tokens_out":6223,"would_cite":true,"duration_ms":66099,"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 CTAO Galactic Plane Survey could detect the brightest Milky Way dark-matter sub-halo at 5σ if TeV-scale WIMPs annihilate to b-bbar with a cross section of about 3×10^-25 cm^3/s—an order of magnitude above the thermal value.","keywords":["dark matter subhalos","indirect dark matter detection","Cherenkov Telescope Array Observatory","Galactic plane survey","gamma-ray astronomy","WIMP annihilation","J-factor"],"falsifier":"A concrete test: after the planned early-phase and long-term GPS exposures are accumulated, if no extended, multiwavelength-quiet source with the NFW-like surface-brightness profile modeled here is found and the derived upper limit on $\\langle\\sigma v\\rangle$ for TeV-mass annihilation into $b\\bar{b}$ remains above $\\sim 3\\times10^{-25}$ cm$^3$/s, the optimistic single-sub-halo claim fails. The population-level prediction would also be tested by searching for the predicted unresolved diffuse excess above $|b|>1^\\circ$ at $\\langle\\sigma v\\rangle\\approx 3\\times10^{-23}$ cm$^3$/s in the resilient scenario; its absence would disfavor the least-disrupted population model.","tokens_in":34245,"feed_emoji":"🔭","tokens_out":8937,"duration_ms":89999,"temperature":0.7,"pith_summary":"The paper asks whether the Cherenkov Telescope Array Observatory's Galactic Plane Survey—the deepest planned TeV survey of the inner five degrees of the plane—can find dark-matter sub-halos, the small bound clumps of dark matter that structure formation predicts but that have no luminous counterpart. It answers yes, under favorable but not exotic conditions: the brightest sub-halo in the Milky Way could be detected at $5\\sigma$ for an annihilation cross section of $\\langle\\sigma v\\rangle \\sim 3\\times10^{-25}$ cm$^3$/s for TeV-scale WIMPs annihilating to $b\\bar{b}$, roughly an order of magnitude above the thermal cross section and within reach of models with resonances or Sommerfeld enhancement. Detecting at least one member of the full population requires $\\langle\\sigma v\\rangle \\sim 10^{-23}$–$10^{-22}$ cm$^3$/s, depending on how resilient sub-halos are to tidal disruption. The study matters because sub-halos provide an indirect-detection channel that is independent of the uncertain dark-matter profile at the Galactic center, and because the crowded plane had previously been considered an unpromising place to look for them.","feed_headline":"Galactic plane survey could catch dark subhalos at 5σ","feed_subtitle":"A TeV WIMP clump needs only about ten times the thermal annihilation cross section to show up.","key_machinery":"Three pieces carry the argument. The first is the SL17 semi-analytic model of the Milky Way sub-halo population, in two variants: 'fragile,' where a sub-halo is destroyed once its tidal radius drops below its scale radius, and 'resilient,' where it survives until the tidal radius is below $0.01$ times the scale radius; from each variant $10^{10}$ population realizations are drawn, and the J-factor—the line-of-sight integral of the squared dark-matter density—of every sub-halo is computed numerically from an NFW profile. The second is a template-based binned Poisson likelihood analysis using simulated CTAO observations built from prod5-v0.16 instrument response functions, with the interstellar emission model 'Base-Max' as the astrophysical background; detection is defined by a test statistic $\\mathrm{TS}=25$ ($5\\sigma$). The third is a position-dependent rescaling law for extended-source sensitivity, $F(\\theta_{68})=A(x)\\sqrt{\\alpha\\theta_{68}^2+\\sigma_{\\mathrm{PSF}}^2}$ with $\\alpha=0.06$, which converts the survey's exposure map into an integrated-flux threshold for any sub-halo of angular extension $\\theta_{68}$, the angular radius containing 68% of the sub-halo's total J-factor.","core_discovery":"This paper establishes that the Cherenkov Telescope Array Observatory's Galactic Plane Survey can act as a discovery channel for dark-matter sub-halos, treating them as extended gamma-ray sources rather than point sources. Using templates for NFW sub-halos drawn from the SL17 semi-analytic population model in its fragile and resilient variants, the authors find that the single brightest sub-halo across many realizations would be detectable at $5\\sigma$ for an annihilation cross section $\\langle\\sigma v\\rangle \\approx 3\\times10^{-25}$ cm$^3$/s when a TeV-scale WIMP annihilates into $b\\bar{b}$, about an order of magnitude above the canonical thermal value. At the population level, detecting at least one sub-halo requires $\\langle\\sigma v\\rangle$ in the range $\\sim 10^{-23}$–$10^{-22}$ cm$^3$/s depending on the survival scenario. The paper also finds that the optimal sky region lies a few degrees above or below the Galactic plane, that the NFW surface-brightness profile allows morphological discrimination from point-like sources and Gaussian astrophysical sources, and that unresolved sub-halos may contribute a diffuse component at least comparable to the unresolved pulsar-wind-nebula population away from the plane.","pith_inferences":["The quoted single-sub-halo cross section is a statement about the luckiest realization; the model's realization scatter spans orders of magnitude, so a dedicated extreme-value treatment of the brightest-clump distribution would sharpen the discovery claim.","Because the SL17 model assumes cuspy NFW inner profiles, applying the same analysis to cored sub-halos or to alternative substructure models such as prompt cusps would directly test how sensitive the GPS reach is to the small-scale dark-matter density structure.","A null GPS search at the quoted population-level cross sections would be more than a constraint on particle physics: it would constrain the survival and concentration of Milky Way sub-halos, potentially favoring the fragile over the resilient scenario.","The morphological discrimination machinery could be turned around: a population of extended, TeV-bright sources with no counterparts, concentrated at $|b|\\approx 2$–$5^\\circ$, would be a distinctive signature of cold-dark-matter substructure that would be hard to mimic with astrophysical source classes."],"forward_implications":["A $5\\sigma$ detection of the brightest individual sub-halo is within reach at $\\langle\\sigma v\\rangle \\sim 3\\times10^{-25}$ cm$^3$/s for TeV-scale DM annihilating to $b\\bar{b}$ without added systematics; 3% systematics worsen this by roughly a factor of two.","For the full SL17 population, at least one sub-halo in the GPS band becomes detectable at $\\langle\\sigma v\\rangle \\sim 3.3\\times10^{-23}$ cm$^3$/s (resilient) and $\\sim 9.7\\times10^{-23}$ cm$^3$/s (fragile) for a 1 TeV WIMP.","The regions a few degrees above or below the plane are the most promising, because source confusion is lower there and the unresolved sub-halo population may dominate the diffuse emission at $|b|>1^\\circ$ in the resilient scenario.","Once detected, a sub-halo can be distinguished from a point-like or Gaussian source at only about 1.5 times the cross section needed for detection, but resolving its inner NFW-like profile into annuli requires about an order of magnitude higher cross section.","If no sub-halo is found, the GPS can still constrain $\\langle\\sigma v\\rangle$ at a level comparable to combined dwarf-spheroidal limits, providing a cross-check independent of Galactic-center profile uncertainties."],"supporting_citations":[{"why":"Supplies the SL17 semi-analytic sub-halo population model whose abundance, spatial distribution, and inner structure set every J-factor and detection threshold.","marker":"[105]"},{"why":"Defines the SL17 fragile and resilient disruption criteria used to generate the two population scenarios.","marker":"[110]"},{"why":"Provides the planned GPS observation schedule, exposure map, and survey simulation approach.","marker":"[85]"},{"why":"Supplies the template-based likelihood and angular-decomposition analysis framework adopted for extended sources.","marker":"[93]"},{"why":"Provides the interstellar emission model 'Base-Max' used as the astrophysical background along the plane.","marker":"[95]"},{"why":"Establishes the treatment of sub-halos as extended sources and the $\\theta_{68}$ angular-size convention.","marker":"[89]"},{"why":"Supplies the gamma-ray annihilation spectra for the $b\\bar{b}$ channel used in the flux calculations.","marker":"[41]"},{"why":"Baseline CTAO sub-halo sensitivity at high latitudes against which the GPS results are compared.","marker":"[84]"},{"why":"Supplies the statistical framework and the Galactic-center comparison sensitivity used for context.","marker":"[18]"}],"fun_headline_variants":["CTAO Galactic plane survey could catch dark subhalos at 5σ","5σ dark subhalo detection possible in CTAO plane survey","Dark subhalos may be visible a few degrees off the Galactic plane","CTAO survey: extended dark subhalos as gamma-ray sources"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Every detection threshold inherits the SL17 semi-analytic model's predictions for how many sub-halos exist, where they sit, and how concentrated their inner density profiles are; if the real Milky Way sub-halo population is more disrupted, less concentrated, or distributed differently, the quoted cross sections shift, potentially by orders of magnitude.","fun_headline_variants_meta":{"raw":{"variants":["CTAO Galactic plane survey could catch dark subhalos at 5σ","5σ dark subhalo detection possible in CTAO plane survey","Dark subhalos may be visible a few degrees off the Galactic plane","CTAO survey: extended dark subhalos as gamma-ray sources"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000248,"raw_usage":{"total_tokens":1668,"prompt_tokens":1186,"completion_tokens":482,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":802,"completion_tokens_details":{"reasoning_tokens":403}},"tokens_in":802,"tokens_out":482,"duration_ms":5945,"temperature":1.0,"reasoning_tokens":403,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T19:39:29.677223+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A concrete test: after the planned early-phase and long-term GPS exposures are accumulated, if no extended, multiwavelength-quiet source with the NFW-like surface-brightness profile modeled here is found and the derived upper limit on $\\langle\\sigma v\\rangle$ for TeV-mass annihilation into $b\\bar{b}$ remains above $\\sim 3\\times10^{-25}$ cm$^3$/s, the optimistic single-sub-halo claim fails. The population-level prediction would also be tested by searching for the predicted unresolved diffuse excess above $|b|>1^\\circ$ at $\\langle\\sigma v\\rangle\\approx 3\\times10^{-23}$ cm$^3$/s in the resilient scenario; its absence would disfavor the least-disrupted population model.","supporting_citations":[{"cited_title":"Modeling dark matter subhalos in a constrained galaxy: Global mass and boosted annihilation profiles","cited_arxiv_id":"1610.02233","evidence_quote":"Supplies the SL17 semi-analytic sub-halo population model whose abundance, spatial distribution, and inner structure set every J-factor and detection threshold."},{"cited_title":"A search for TeV gamma-ray emission from SNRs, pulsars and unidentified GeV sources in the Galactic plane in the longitude range between -2 deg and 85 deg","cited_arxiv_id":"astro-ph/0209360","evidence_quote":"Supplies the template-based likelihood and angular-decomposition analysis framework adopted for extended sources."},{"cited_title":"Dark matter substructure modelling and sensitivity of the Cherenkov Telescope Array to Galactic dark halos","cited_arxiv_id":"1606.04898","evidence_quote":"Baseline CTAO sub-halo sensitivity at high latitudes against which the GPS results are compared."}],"review_version":1}