{"id":"c8df0cdc-4646-46e4-948c-783e8db0b0e8","arxiv_id":"2506.04074","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A simulation predicts that 500 hours of CTAO-North four-LST observations at large zenith angles can discriminate between dark-matter, millisecond-pulsar, and proton-induced models of the Galactic Center TeV source HESS J1745-290 by measuring the sharpness of its spectral cutoff.","lead":"This paper simulates future observations of the Milky Way's center with the four large telescopes of the CTAO-North array, estimating how quickly they could distinguish dark-matter, pulsar, and proton explanations of the TeV gamma-ray source HESS J1745-290. The forecast suggests about 500 hours of high-zenith-angle data could separate the dark-matter model from the others at 5-sigma, with a possible 3-sigma hint within the first year.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The forecast tests pure SEPL templates, not the physical emission models it claims to discriminate; if the true spectra yield intermediate beta values like the observed 1.53, the claimed rejection powers are not established.","rationale":"The reader's weakest assumption identified the same core issue: the forecast assumes the true spectrum is a single super-exponential power law whose sharpness parameter uniquely identifies the emission scenario, with no contamination from other gamma-ray sources. My stress-test sharpens this by pointing out that the mock data are generated from pure SEPL templates (Figure 1) rather than from the physical models described in Section 2, making the likelihood-ratio test a self-consistency check of SEPL fitting rather than a test of physical scenario discrimination. This is load-bearing because the paper's strongest claim explicitly maps beta=2 onto the DM-spike model, beta=0.6 onto hadronic emission, and beta=1 onto MSPs; if the true physical templates produce different effective beta values, the claimed 5-sigma rejections do not apply. The paper itself reports a best-fit beta of 1.53±0.28 from current data, which lies between the simulated values and is not tested, directly undermining the assumed one-to-one mapping. I agree with the reader that this warrants a conditional verdict, but I would not move to reject: the simulation methodology is standard, the use of public IRFs and injected systematics is transparent, and the paper is careful to state its assumptions even where they are idealized. The concrete test I propose would settle whether the concern lands by checking the forecast against the physical templates the paper actually discusses.","tokens_in":20,"tokens_out":6157,"duration_ms":122509,"concrete_test":"Regenerate the mock count cubes using the physical templates described in Section 2, namely the DM-spike plus CMZ multi-component model of Figure 2, the proton-injection plus pp-interaction spectrum from [46], and the MSP IC spectrum from [18], as true inputs while keeping the SEPL fitting model and the same IRFs and background treatment. Compare the resulting 68% containment bands of fitted beta at 120 h and 500 h, and the TS distributions for fixed beta=0.6, 1.0, and 2.0. If the DM+CMZ template yields a median fitted beta near 1.5, or if the TS values for distinguishing beta=2 drop below the claimed thresholds, the forecast's scenario-discrimination power is an artifact of the pure-SEPL template choice rather than a robust property of the physical models.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is that 500 h of four-LST observations can reject beta=2 (DM spike) at >5 sigma when the true spectrum is hadronic (beta=0.6) or pulsar-like (beta=1). This requires that the three scenarios map uniquely onto beta values of 2.0, 0.6, and 1.0, and that the mock data in Section 3.2 are generated from physically faithful templates. In fact, Section 3.2 generates count cubes from the pure SEPL 'exponential templates' shown in Figure 1 (beta=0.6, 1, 2), not from the physically motivated models of Section 2: the DM+CMZ multi-component model of Figure 2, the proton-injection plus pp-interaction spectrum from [46], or the MSP inverse-Compton spectrum from [18]. The likelihood-ratio test in Section 4 (Eq. 4.1) therefore measures how well a SEPL fit can reject an alternative SEPL, which is a self-consistency check of the fitting procedure rather than a test of discrimination among physical scenarios. Section 2.1's own combined fit gives beta=1.53±0.28, an intermediate value not included in the simulated grid; if the true input model were the DM+CMZ template, the reconstructed beta could be close to 1.5, collapsing the separation between scenarios. The field of view also contains known gamma-ray sources (Sgr A East, G359.95-0.04, CMZ diffuse emission), which the paper explicitly omits in Section 3.2 ('no additional gamma-ray source is presumed within the FoV'). Thus the headline rejection powers are conditional on a single-SEPL, empty-FoV assumption that is not justified by the current data or the source region.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper fits super-exponential power-law (SEPL, Eq. 1.1) models to archival H.E.S.S., MAGIC, and VERITAS spectra of HESS J1745-290, reporting a preferred cutoff sharpness beta = 1.53 +/- 0.28 (statistical only). It then simulates mock CTAO-North four-LST observations at large zenith angle (60 degrees) using public prod5 IRFs and gammapy, generating ~800 Poisson realizations for true SEPL templates with beta = 0.6, 1.0, and 2.0, with optional injected response systematics. Using Cash-statistic fits and a likelihood-ratio TS (Eq. 4.1), it finds that ~500 hours can reject beta = 2 when the true beta is 0.6 or 1.0 at more than 5 sigma (median), and that ~120 hours can give a roughly 3-sigma discrimination between beta = 2 and beta = 1. The paper concludes that CTAO-N LSTs can provide timely discrimination among DM-spike, hadronic, and MSP emission scenarios for the Galactic Center source.","tokens_in":14149,"tokens_out":6601,"duration_ms":63497,"significance":"The study is timely and methodologically solid in its technical core: it uses publicly available CTAO instrument response functions, a standard forward-folding framework (gammapy), many Monte Carlo realizations, and an explicit injection of response systematic uncertainties. The forecast that the sharpness of the multi-TeV cutoff can be measured well enough to reject beta = 2 at high significance would be a useful, falsifiable prediction for the upcoming four-LST array. However, the transfer of these sensitivities to the three physical scenarios is not yet established, because the simulated truths are pure SEPL templates rather than the multi-component models of Section 2, and the currently preferred beta = 1.53 lies outside the simulated grid. With additional simulations based on the physical templates and an intermediate-beta case, the paper could deliver on its main claim.","major_comments":[{"comment":"The mock data are generated from pure SEPL templates with beta = 0.6, 1.0, and 2.0, not from the physical models introduced in Section 2. The DM model of Section 2.2 is a multi-component (DM bbar/tau+tau plus CMZ-like diffuse) fit, the proton model of Section 2.3 is a pp-interaction spectrum, and the MSP model of Section 2.4 is an inverse-Compton spectrum; none is a pure SEPL. The likelihood-ratio test in Eq. (4.1) therefore measures how well a SEPL fit separates three artificial SEPL spectra, not how well it separates the physical emission scenarios claimed in the abstract. This is compounded by the fact that the joint fit in Section 2.1 gives beta = 1.53 +/- 0.28, an intermediate value not represented among the simulated truths; if the real spectrum is a blend (e.g., DM plus hadronic emission), the reconstructed beta could lie near 1.5 and the reported rejection powers would not transfer. Please regenerate the mock data from the physical templates (or at least add an intermediate-beta case such as beta = 1.5) and re-evaluate the TS distributions.","section":"§3.2 and §4"},{"comment":"The statement 'no additional gamma-ray source is presumed within the FoV' is not realistic for the Galactic Center field. Sgr A East, the pulsar-wind nebula candidate G359.95-0.04, and CMZ diffuse gamma-ray emission are all known sources in the region and are discussed in Sections 1 and 2. They are not removed by the cited background-estimation methods (RingBackground or Reflected-Region-Background), which only estimate the instrumental residual background. The forecast therefore applies to a source-free field of view; including at least a CMZ diffuse component and a check for Sgr A East contamination could change the recovered beta and the discrimination significances. I ask the authors to add such components to at least one mock setup and report the effect on the TS values.","section":"§3.2 (FoV assumptions)"},{"comment":"The significance claims are based on the median TS from roughly 800 realizations, as stated in the Figure 4 caption ('Lines within the bands show the TS median value'), rather than on the fraction of realizations exceeding the Wilks-theorem thresholds of TS = 9 and 25. A median TS of 9 means only about half of the realizations would yield a 3-sigma rejection, so the quoted 'rejected at 4.7 sigma' and 'preliminary 3-sigma hint' should be reported as median expectations, and the paper should give the distribution of TS values (e.g., the percentage of realizations above TS = 9 and 25). This is directly relevant to the abstract's claim of a first-year 3-sigma hint.","section":"§4, Figure 4"}],"minor_comments":[{"comment":"There are several typos: 'CT AO-N', 'LST s', 'T e V', 'senarios', and 'can well described' should be corrected.","section":"Abstract and title"},{"comment":"The sentence 'any parameter is unfrozen in neither the spatial components nor the background model' is grammatically ambiguous; please clarify which parameters are kept free in the fit.","section":"§3.2"},{"comment":"When reporting the SEPL fits at fixed beta, the cutoff energies are given as 'Ec = 7.8+1.5−1.1 TeV, 15.7+1.9−1.5 TeV, and 19.2+1.8−1.5 TeV' with inconsistent formatting; please standardize.","section":"§2.1 / Figure 1"},{"comment":"Refs. [35] and [38] appear to be the same GRAVITY mass-distribution paper with different journal formatting; please deduplicate or differentiate.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is suitable for JCAP. The technical implementation is sound, but the mapping from simulated SEPL templates to physical scenarios is the central gap. I would like the revision to include physical-template simulations; if the authors prefer to limit the scope to 'sensitivity to beta', they should rewrite the abstract and conclusions accordingly. No concerns about citation practices."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a read. This is the first quantitative forecast I know of for the four-LST CTAO-North sub-array observing HESS J1745-290 at large zenith angles, using public prod5 IRFs and a proper forward-folding simulation with injected systematics. The methodology is standard and transparent: mock count cubes, Poisson/Cash statistics, ~800 realizations per setup, and explicit treatment of energy-scale, effective-area, and resolution systematics. The main result — that ~500 h can reject beta=2 (DM-spike-like) against beta=0.6 or 1 at >5 sigma, and that a 3-sigma hint can come in 120 h — is a useful planning input for the near-term LST program.\n\nThe soft spots are real but not fatal. First, the mock true spectra are the pure SEPL templates with beta=0.6, 1, and 2, not the fuller physical models of Sec. 2 (the DM+CMZ multi-component fit, or the specific proton/MSP spectral shapes). The paper does argue that each scenario is reasonably approximated by those beta values, and that is defensible for a sensitivity study; the stress-test claim that this is merely a self-consistency check goes too far. It is a standard benchmark approach. But the authors' own combined fit to current data gives beta=1.53±0.28, an intermediate value they do not simulate. If the true spectrum sits between the benchmarks, the discrimination powers they quote will not transfer directly. They should add an intermediate-beta case, or at least show the beta containment as a function of true beta, to make the forecast robust.\n\nSecond, the FoV is assumed empty apart from the central point source and residual hadronic background. The GC region contains Sgr A East, G359.95-0.04, and CMZ diffuse emission. The paper states the assumption explicitly and uses background-subtraction methods, but does not quantify how such contaminants would bias the beta measurement. That is worth a caveat and maybe a simple check with a second source included.\n\nThird, minor: the text has a couple of typos (e.g., 'CT AO-N' in the title, 'senarios'), but nothing substantive.\n\nOverall: a solid, standard sensitivity study with a clear and important target. It deserves a serious referee. Recommend that the authors add an intermediate-beta simulation and a contamination caveat before publication. I would send it to review, not desk reject.","headline":"A solid, timely sensitivity forecast for CTAO-North's four-LST observations of the Galactic Center; headline claims hold for the assumed SEPL templates, but intermediate-beta and source-contamination cases are not tested.","tokens_in":14744,"tokens_out":4359,"would_cite":true,"duration_ms":44626,"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-North four-LST array can discriminate dark-matter, millisecond-pulsar, and hadronic explanations of the Galactic Center's multi-TeV source by measuring the cutoff sharpness parameter β.","keywords":["gamma ray experiments","dark matter experiments","CTAO-North","large zenith angle observations","Galactic Center","HESS J1745-290","super-exponential cutoff","millisecond pulsars"],"falsifier":"If 500 hours of actual four-LST data recover β with a 68% interval that overlaps the paper's current best fit of 1.53 ± 0.28 rather than separating the values 0.6, 1, and 2, the discrimination claim fails; equivalently, a detection of extended emission from Sgr A East within the point-source region would break the single-source assumption.","tokens_in":1700,"feed_emoji":"🔭","tokens_out":2646,"duration_ms":68926,"temperature":0.7,"pith_summary":"The paper argues that the shape of the energy cutoff in the gamma-ray spectrum of HESS J1745-290 carries enough information to identify the process producing the emission, and that the incoming four-LST array at CTAO-North can extract that shape within a few years. A super-exponential power law with sharpness parameter β is fitted to current H.E.S.S., MAGIC, and VERITAS data; the data prefer β = 1.53 ± 0.28, too broad to choose between β = 2 (dark-matter spike), β = 1 (millisecond pulsars), and β ≈ 0.6 (hadronic). Simulating large-zenith-angle observations with realistic instrument response, the paper finds that 120 hours can give a first 3σ hint and roughly 500 hours, possibly across several years, can reject the wrong scenarios at more than 5σ. A successful measurement would directly discriminate dark-matter annihilation from astrophysical emission at the center of the Milky Way.","feed_headline":"Four LSTs can reveal what powers the Galactic Center's TeV source","feed_subtitle":"A 500-hour campaign at 60° zenith can measure the cutoff sharpness that separates dark matter from pulsars and protons.","key_machinery":"The central object is the super-exponential cutoff power law $\\Phi(E) = \\Phi_0 (E/1\\,\\mathrm{TeV})^{-\\Gamma} \\exp[-(E/E_c)^\\beta]$, whose shape parameter $\\beta$ is the discriminator: $\\beta = 2$ for a dark-matter spike around Sgr A*, $\\beta = 1$ for inverse-Compton emission from millisecond pulsars, and $\\beta \\approx 0.6$ for proton-proton interactions. The observational mechanism is the large-zenith-angle technique: from CTAO-North the Galactic Center culminates near 58°, and at 60° zenith the four-LST array reaches an effective area of about $9.0\\times10^5$ m² at 10 TeV with energies up to about 100 TeV and an energy threshold of 400 GeV. The analysis pipeline forward-folds response-convolved models into count cubes and uses Cash statistics, with systematic uncertainties simulated by perturbing the effective area, energy scale, energy resolution, and background amplitude.","core_discovery":"The central claim is that the sharpness of the multi-TeV spectral cutoff can be measured well enough by four LSTs at CTAO-North to separate the three viable emission models of the Galactic Center source. With β fixed to 0.6, 1, or 2, and using the official 60° zenith instrument response functions, the paper simulates about 800 realizations per setup and computes likelihood-ratio test statistics between different assumed true models. It finds that after 500 hours a true β = 2 (DM spike) disfavors β = 0.6 and β = 1 at high significance, and a true β = 0.6 rejects β = 2 at about 9.7σ; after 120 hours, β = 2 can already be separated from β = 1 at about 3σ. Adding realistic response systematics widens the β containment bands by 20 to 40 percent but does not erase the discrimination. The paper also notes that the current combined fit gives β = 1.53 ± 0.28, an intermediate value that no pure scenario in the study reproduces.","pith_inferences":["The paper's own best fit of β = 1.53 ± 0.28 suggests that the real data may land between the three pure scenarios; an obvious extension is testing two-component models, such as a dark-matter contribution plus diffuse emission, with the same simulated likelihood machinery.","The same large-zenith-angle four-LST setup could be applied to other northern very-high-energy sources with cutoff features, where β is also theoretically constrained and could separate emission mechanisms.","A cross-check with CTAO-South, which views the Galactic Center at much lower zenith angles, would tighten systematic control and could test whether the large-zenith-angle energy-scale systematics bias the recovered β.","The result implies that CTAO-North can play a meaningful interim role in dark-matter searches at the Galactic Center before the more powerful Southern array is fully deployed, which may change near-term observation scheduling."],"forward_implications":["If the real spectrum after 500 hours has β near 2, the dark-matter spike interpretation would be strongly favored over pulsar and hadronic models.","If β is recovered near 1, the central stellar cluster's millisecond pulsars become the leading explanation for the multi-TeV emission.","If β is recovered near 0.6, the emission would support cosmic-ray proton interactions in the inner 10 parsecs of the Galactic Center.","A first-year 120-hour campaign can already provide about a 3σ hint, which motivates scheduling large-zenith-angle Galactic Center observations early in the four-LST array's operation.","Because systematic uncertainties broaden but do not erase the separation, the discrimination result appears robust at the level of the current instrument-response knowledge."],"supporting_citations":[{"why":"Establishes that the super-exponential cutoff shape parameter distinguishes annihilating dark matter (β > 1) from proton-induced emission (β < 1).","marker":"[19]"},{"why":"Provides the H.E.S.S. spectrum of HESS J1745-290 that anchors the current multi-TeV data and the cutoff measurement.","marker":"[6]"},{"why":"Supplies the MAGIC measurement of the Galactic Center region used in the combined spectral fit.","marker":"[23]"},{"why":"Supplies the VERITAS multi-TeV Galactic Center measurement used in the combined spectral fit.","marker":"[24]"},{"why":"Provides the earlier analysis of the VHE spectrum and future prospects, including the dark-matter interpretation and the SEPL framework.","marker":"[27]"},{"why":"Provides the public CTAO instrument response functions at 60° zenith used for all simulated four-LST observations.","marker":"[58]"},{"why":"Generates the millisecond-pulsar inverse-Compton spectrum that the paper identifies with β = 1.","marker":"[18]"},{"why":"Provides the parametrization of gamma-ray production in proton-proton interactions that yields β ≈ 0.6.","marker":"[47]"}],"fun_headline_variants":["LSTs to test dark matter vs pulsars at Galactic Center","500h with LSTs to resolve GC spectral cutoff","Four LSTs to separate GC emission models","CTAO-N LSTs may crack Galactic Center TeV puzzle"],"cache_read_input_tokens":16768,"weakest_assumption_plain":"The forecast assumes the true spectrum is one pure super-exponential power law whose sharpness β matches exactly one emission scenario, and that no other gamma-ray source contaminates the field of view.","fun_headline_variants_meta":{"raw":{"variants":["LSTs to test dark matter vs pulsars at Galactic Center","500h with LSTs to resolve GC spectral cutoff","Four LSTs to separate GC emission models","CTAO-N LSTs may crack Galactic Center TeV puzzle"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000355,"raw_usage":{"total_tokens":2032,"prompt_tokens":1150,"completion_tokens":882,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":766,"completion_tokens_details":{"reasoning_tokens":816}},"tokens_in":766,"tokens_out":882,"duration_ms":8673,"temperature":1.0,"reasoning_tokens":816,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T10:47:51.710584+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"If 500 hours of actual four-LST data recover β with a 68% interval that overlaps the paper's current best fit of 1.53 ± 0.28 rather than separating the values 0.6, 1, and 2, the discrimination claim fails; equivalently, a detection of extended emission from Sgr A East within the point-source region would break the single-source assumption.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes that the super-exponential cutoff shape parameter distinguishes annihilating dark matter (β > 1) from proton-induced emission (β < 1)."},{"cited_title":"Study of the very high energy gamma-ray spectrum from the Galactic Center and future prospects","cited_arxiv_id":"1610.10003","evidence_quote":"Provides the earlier analysis of the VHE spectrum and future prospects, including the dark-matter interpretation and the SEPL framework."},{"cited_title":"Gamma-rays from millisecond pulsar population within the central stellar cluster in the Galactic Center","cited_arxiv_id":"1306.4760","evidence_quote":"Generates the millisecond-pulsar inverse-Compton spectrum that the paper identifies with β = 1."}],"review_version":1}