{"id":"e3727a03-1015-4ae5-8dbc-382f45e4fb0b","arxiv_id":"1909.00222","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A 10-hour MAGIC observation of the Galactic Halo, using same-night ON/OFF pairs at different Galactic latitudes, yields a 95% CL dark matter decay lifetime limit above 10^26 seconds for TeV masses in the b-bbar channel.","lead":"MAGIC telescopes used paired observations of two regions of the Milky Way halo to cancel ordinary gamma-ray backgrounds and search for dark matter decay. With 10 hours of data, they set a lower limit on the dark matter lifetime of 10^26 seconds for TeV-scale particles decaying to b quarks.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 10-h Galactic Halo limit assumes the residual Galactic diffuse γ-ray background cancels in the high-ΔJ ON/OFF pairs; this is asserted, not quantified, and the control-sample systematic does not cover it.","rationale":"The paper presents an ICRC proceedings describing a promising new ON/OFF-pair method for DM decay searches. The result is a limit, not a detection, so the main risk is an unmodeled background that biases the limit. The single most load-bearing point is the cancellation of the Galactic diffuse emission: the central equation ΔR = R_DM(φ1) − R_DM(φ2) hinges entirely on this, and the paper admits the anisotropy is only 'mild' without quantifying it. Because the high-ΔJ selection is geographically correlated with the Galactic Center, the diffuse residual is unlikely to be strictly zero, and its sign is not guaranteed. The σ_syst measurement, while carefully done on low-ΔJ and archival pairs, does not validate the foreground cancellation in the search configuration. This does not impugn the method's potential; rather, it means the quoted lifetime limit should be treated as conditional on foreground validation. The reader's weakest_assumption identifies exactly this issue, so I agree. A concrete test with a Fermi-LAT interstellar emission template for the actual pointing pairs would settle whether the concern lands. If the residual is small compared to the DM signal, the limit stands; if not, the analysis needs a foreground model or an expanded systematic. The verdict CONDITIONAL is appropriate: the paper should be accepted only with the foreground check demonstrated or explicitly flagged as an open systematic.","tokens_in":6810,"tokens_out":4599,"duration_ms":49973,"concrete_test":"For the 10 ON/OFF pairs used in Sec. 4, compute the expected γ-ray rate difference from the Galactic diffuse emission using the Fermi-LAT interstellar emission model (gll_iem_v06) or a comparable template, fold it with the MAGIC effective area, energy threshold (E > 60 GeV), and PSF, and compare |ΔR_Gal−γ| with the DM decay signal expected at the reported 95% CL limit for m_DM = 100 TeV, b-bbar channel. Repeat for the low-ΔJ control pairs. If |ΔR_Gal−γ| exceeds ~20% of the DM signal in any search pair, or if the control pairs bracket a different residual than the search pairs, the reported lifetime limit is biased and the analysis must include a foreground model or assign a dedicated systematic.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central result is a 95% CL lower limit on the DM lifetime from the difference between ON and OFF fields (Sec. 2). The derivation sets ΔR(φ1,φ2) = R_DM(φ1) − R_DM(φ2) by assuming R_Gal−γ(l1,b1) − R_Gal−γ(l2,b2) is negligible, justified only by the statement that the Galactic diffuse emission is 'mildly anisotropic for |b| > 10°' (Sec. 2). No foreground model, template map, or quantitative bound is given. The signal is proportional to ΔJ, and the selected high-ΔJ pairs have ON fields closer to the Galactic Center than OFF fields, so the Galactic diffuse gradient is expected to correlate with ΔJ. The residual could be positive or negative depending on local structures (Fermi Bubbles, Loop I, unresolved sources). If negative, the derived lower limit on τ_DM is artificially strong, and the claim that lifetimes below 10^26 s are excluded at 100 TeV could be false. The σ_syst = 4.8% is measured on low-ΔJ control pairs (ΔJ_syst ≈ 4% of the search average, Sec. 3) and on archival pairs with ΔJ ≈ 0; this does not probe the diffuse-background residual in the actual search configuration. Thus the load-bearing assumption is that the residual Galactic diffuse anisotropy is sub-dominant to the DM signal in the 10 high-ΔJ pairs used for the limit.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript reports a search for gamma-ray signals from dark-matter (DM) decay in the Milky Way halo using 20 hours of MAGIC observations taken in 2018. The method compares same-night ON/OFF pointings at different angular distances from the Galactic Center, with matched azimuth and zenith paths, and assumes that all non-DM diffuse components cancel in the ON-OFF difference. Ten hours with the lowest ΔJ are used to measure the systematic uncertainty σ_syst = 4.8%, and the remaining 10 hours with the highest ΔJ are used in a binned likelihood to set 95% CL lower limits on the DM lifetime for the b bbar channel; the strongest limit is τ_DM > 10^26 s at m_DM = 100 TeV, which the authors compare with a 200-hour Perseus observation.","tokens_in":7135,"tokens_out":4084,"duration_ms":40048,"significance":"If the result holds, the paper demonstrates that a dedicated matched-pair ON-OFF strategy can produce TeV-scale DM decay limits competitive with much longer exposures, and the methodology is of clear interest for future instruments such as CTA. The strengths of the paper include the careful matched azimuth/zenith data-taking, the cross-check of the systematic uncertainty against archival data, and the explicit demonstration that the J-factor is nearly model-independent in the selected sky region. The result is not a detection, and its main risk is statistical and systematic rather than the central limit itself; however, the foreground-cancellation assumption is load-bearing and is not quantitatively supported.","major_comments":[{"comment":"The central reduction ΔR(φ1,φ2) = R_DM(φ1) − R_DM(φ2) rests on the claim that R_Gal−γ(l1,b1) − R_Gal−γ(l2,b2) is negligible because the Galactic diffuse emission is 'mildly anisotropic' for |b| > 10°. This is asserted rather than demonstrated. Since the ON positions are closer to the Galactic Center than the OFF positions, the Galactic diffuse gradient is expected to correlate with ΔJ, so the residual could be of either sign and comparable to the DM signal. Please quantify this residual using a Galactic diffuse foreground model or Fermi-LAT template maps for the specific high-ΔJ pairs used in Sec. 4, and show that it is small compared to CΔJ. If a negative residual is possible, the derived lifetime limit could be artificially strong.","section":"Sec. 2, equation for R(l,b)"},{"comment":"The σ_syst = 4.8% used in the likelihood is measured on control samples with ΔJ_syst ≈ 4% of the search average and on archival pairs with ΔJ ≈ 0. These configurations do not reproduce the spatial foreground difference of the search pairs; they constrain the stability of the ON/OFF acceptance but not the residual Galactic diffuse anisotropy term in the equation for ΔR. The likelihood in Eq. (4.1) propagates σ_syst only into the κ uncertainty, so the final limit contains no contribution from the foreground-cancellation assumption. A quantitative bound on that residual is required before the reported limit can be considered robust.","section":"Sec. 3, systematic uncertainty"},{"comment":"The likelihood includes a term J(ΔJ|ΔJ_obs, σ_ΔJ) but sets σ_ΔJ = 0. The paper argues that the J-factor is robust because three halo profiles agree in the region of interest; nevertheless, the Milky Way halo parameters are not infinitely well known, and the choice σ_ΔJ = 0 should be justified with a quantitative estimate of the residual spread among models or a statement of which halo-parameter uncertainties were considered. As written, the limit omits any uncertainty in the astrophysical factor.","section":"Sec. 4, Eq. (4.1)"}],"minor_comments":[{"comment":"The lifetime limit is rendered as '10^26 s' in the abstract but as '1026 s' in the body text; please fix the typesetting of the superscript.","section":"Abstract and Sec. 4"},{"comment":"The right panel of Fig. 2 is discussed in terms of orange and magenta histograms, but the figure caption does not identify which histogram corresponds to which selection; the caption should be self-contained.","section":"Fig. 2"},{"comment":"The definition of R is given as an inline formula; please define it in a numbered equation and clarify that it is expressed as a percentage, and state how the quoted uncertainty on σ_tot was derived from the sample variance with the given number of nights.","section":"Sec. 3"},{"comment":"The notation τ_LL_DM is used without definition; please define it as the lower-limit lifetime or explain the subscript 'LL'.","section":"Sec. 4"},{"comment":"Reference [7] to the 3FHL catalog lacks volume and page information, and [10] should include the full collaboration name; please complete the bibliographic entries.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"This is a proceedings-style paper with preliminary results. The main technical risk is the unquantified cancellation of the Galactic diffuse foreground; if the authors can provide a quantitative estimate using Fermi diffuse templates or an independent cross-check on high-ΔJ fields, the central claim would be much more convincing. I do not see a novelty or attribution concern, and the paper is honest in labeling the results as preliminary. The appropriate next step is major revision with emphasis on the foreground residual and the σ_ΔJ assumption."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a preliminary ICRC proceedings, but the core idea is new and worth your attention. Instead of pointing at the GC or a cluster, the MAGIC authors pair same-night, same-azimuth/zenith ON/OFF observations at different angular distances from the GC, selecting pairs with large Delta-J. The isotropic components cancel in the subtraction, and the Galactic diffuse gamma-ray foreground is assumed to be negligible for |b|>10. Using 10 hours of the highest-Delta-J data, they derive a 95% CL lower limit on the DM lifetime that reaches ~10^26 s at 100 TeV for the bb channel, competitive with 200 hours on Perseus. The statistical treatment is standard (likelihood from Aleksic et al.), the cuts are optimized, the energy-independence of the systematic is checked, and no detection is claimed. Credit where due: the same-night paired strategy is a smart way to beat the usual IACT background problem for a diffuse target.\n\nThe soft spots are real. The load-bearing assumption is that the residual Galactic diffuse emission cancels in the high-Delta-J pairs. The paper states it is 'mildly anisotropic' but gives no foreground model, template map, or quantitative bound. Since the ON fields are closer to the GC than the OFF fields, the same spatial gradient that produces the DM signal could also produce a non-DM residual, with either sign. A negative residual would make the lifetime limit artificially strong. The systematic uncertainty is measured on low-Delta-J control pairs (Delta-J about 4% of the search average) and on archival pairs with Delta-J≈0; that does not probe the actual diffuse-background difference in the search configuration. The 4.8% is applied to the normalization kappa, but a foreground residual would affect the background term differently. For a conference proceedings these are acceptable caveats, but they matter if the limit is to be treated as final.\n\nThis is not a fatal flaw in the method, and the paper is honest about its limits. A follow-up with a template-based foreground estimate or with OFF positions that better isolate the DM gradient would answer the concern. I would send it to peer review — it deserves referee time — but with the expectation that the foreground assumption be quantified or the systematic study be redone on the actual search pairs. I'd cite it as a method note and might bring it to a reading group, but I would not yet quote the lifetime limit as a robust result.","headline":"A genuinely new same-night ON/OFF pair method for Galactic Halo dark matter decay searches gives an impressively competitive limit, but the diffuse-foreground cancellation is asserted rather than demonstrated.","tokens_in":7684,"tokens_out":1914,"would_cite":true,"duration_ms":20508,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["95.35.+d"],"model":"deepseek-v4-flash","headline":"Pairwise subtraction of matched MAGIC observations of the Galactic halo excludes dark-matter decay to bottom-quark pairs with lifetime below about $10^{26}$ s at a mass of 100 TeV.","keywords":["dark matter decay","Galactic halo","gamma-ray astronomy","Cherenkov telescopes","indirect dark matter search","lifetime limit","ON-OFF observation","J-factor"],"falsifier":"Compute the expected residual $\\Delta R_{\\rm Gal-\\gamma}$ between the selected ON and OFF fields using a standard model of interstellar gas and cosmic-ray distributions; if the predicted diffuse background difference is comparable to or larger than the 4.8% systematic uncertainty or the observed residual, then the derived lifetime limits would be biased.","tokens_in":6598,"feed_emoji":"🌌","tokens_out":7425,"duration_ms":62016,"temperature":0.7,"pith_summary":"This paper proposes that the gamma-ray decay signal of dark matter in the Milky Way's halo can be isolated by subtracting pairs of observations taken at different angular distances from the Galactic Centre under identical observational conditions. All isotropic diffuse components cancel in the subtraction, and the mildly anisotropic Galactic diffuse emission is argued to be negligible for the chosen fields, leaving a residual proportional to the difference in dark-matter J-factor. Applied to 10 hours of MAGIC 2018 data, the method sets a 95% confidence-level lower limit on the dark-matter lifetime of $\\tau_{\\rm DM} > 10^{26}$ s in the $\\mathrm{b}\\bar{\\mathrm{b}}$ decay channel at $m_{\\rm DM}=100$ TeV, competitive with limits from 200 hours on the Perseus cluster. A systematic uncertainty of 4.8%, measured with a separate 10-hour sample, is folded into the likelihood. The result matters because it shows that a modest investment of halo observing time can probe TeV-scale decaying dark matter as well as much longer cluster campaigns.","feed_headline":"10 hours rival 200 in dark-matter decay limits","feed_subtitle":"Matched halo sky pairs cancel diffuse backgrounds and set a TeV lifetime limit past 10^26 seconds.","key_machinery":"The load-bearing element is the matched ON/OFF pair: each ON observation at low Galactic longitude (high expected dark-matter flux) is paired with an OFF observation at higher longitude during the same night, following the same azimuth and zenith path, with both fields satisfying $|b|>10^\\circ$. The normalization factor $\\kappa$ between ON and OFF is derived from high-hadronness events outside the signal region, and the residual $R=2(N_{\\rm ON}-N_{\\rm OFF}/\\kappa)/(N_{\\rm ON}+N_{\\rm OFF}/\\kappa)$ is the test statistic whose width, measured on a low-$\\Delta J$ control sample, defines the systematic uncertainty. The signal term is proportional to $\\Delta J(\\phi_1,\\phi_2)$, the difference in J-factors, which is nearly model-independent for decay at $\\phi>10^\\circ$ across NFW, Einasto, and isothermal profiles.","core_discovery":"MAGIC's pairwise ON-OFF observation strategy for the Galactic halo can make the residual event rate between two sky positions depend only on the dark-matter decay contribution, $\\Delta R(\\phi_1,\\phi_2) = R_{\\rm DM}(\\phi_1)-R_{\\rm DM}(\\phi_2)$, once isotropic backgrounds are removed and the anisotropic Galactic gamma-ray background is assumed negligible. With 10 hours of observations selected for large $\\Delta J$, the analysis yields a 95% CL lower limit on the $\\mathrm{b}\\bar{\\mathrm{b}}$ decay lifetime of $10^{26}$ s at 100 TeV, and limits as constraining as the 200-hour Perseus cluster observation for masses up to 10 TeV. The paper also reports an energy-independent systematic uncertainty of $\\sigma_{\\rm syst}=4.8\\pm1.0\\%$ in the background-estimation procedure.","pith_inferences":["The same pairwise cancellation could be adapted to dark-matter annihilation, where the J-factor contrast between matched fields is stronger, though the brighter astrophysical foregrounds at small angular separation would need a quantitative foreground model.","The claim that the residual Galactic diffuse gamma-ray emission is negligible could be tested directly by measuring the ON-OFF residual as a function of interstellar gas column density in the selected fields; a correlation would indicate a foreground bias in the lifetime limit.","Comparable pairwise analyses by other instruments using different pointings could combine with MAGIC's limit to push the lifetime bound beyond $10^{26}$ s without requiring the full observation time of a deep cluster survey."],"forward_implications":["Ten hours of Galactic-halo observation with MAGIC exclude $\\mathrm{b}\\bar{\\mathrm{b}}$ dark-matter decay lifetimes below about $10^{26}$ s at 100 TeV at 95% confidence.","The pairwise subtraction method reaches the same sensitivity as a 200-hour Perseus cluster campaign for masses up to 10 TeV, using 5% of the observation time.","Extending the data set to 40 hours is projected to improve the lifetime constraints by at least a factor of two.","The method, including its measured 4.8% systematic uncertainty, transfers naturally to future wide-field Cherenkov observatories.","Because the decay J-factor is nearly identical for cusped and cored halo profiles at $\\phi>10^\\circ$, the limits are robust against the Milky Way halo modelling uncertainty."],"supporting_citations":[{"why":"Supplies the full likelihood formalism used to convert the ON-OFF counts into a dark-matter lifetime limit.","marker":"[8]"},{"why":"Provides the J-factor and decay gamma-ray flux calculations for the NFW, Einasto, and isothermal halo profiles.","marker":"[5]"},{"why":"The hadronness random-forest classifier that defines the signal and normalization regions.","marker":"[6]"},{"why":"The 3FHL catalog used to veto known gamma-ray sources inside the selected pointing fields.","marker":"[7]"},{"why":"The 200-hour Perseus cluster observation that serves as the benchmark sensitivity comparison.","marker":"[9]"},{"why":"The Fermi-LAT inner-halo lifetime limits shown for comparison at masses below 10 TeV.","marker":"[10]"}],"fun_headline_variants":["10 hours of halo data match 200 from Perseus","Matched halo sky pairs cancel diffuse backgrounds for DM decay","Novel ON-OFF method yields competitive dark matter decay limits","Dark matter decay lifetime limit above 10^26 seconds","Galactic halo observations set TeV DM decay limit in 10 hours"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The method assumes that the Galactic diffuse gamma-ray background is effectively the same in the ON and OFF fields, so that its difference is negligible compared with the dark-matter decay signal; the paper does not quantify this with a foreground emission model or template map.","fun_headline_variants_meta":{"raw":{"variants":["10 hours of halo data match 200 from Perseus","Matched halo sky pairs cancel diffuse backgrounds for DM decay","Novel ON-OFF method yields competitive dark matter decay limits","Dark matter decay lifetime limit above 10^26 seconds","Galactic halo observations set TeV DM decay limit in 10 hours"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000408,"raw_usage":{"total_tokens":2132,"prompt_tokens":973,"completion_tokens":1159,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":589,"completion_tokens_details":{"reasoning_tokens":1073}},"tokens_in":589,"tokens_out":1159,"duration_ms":39468,"temperature":1.0,"reasoning_tokens":1073,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:57:56.384110+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the expected residual $\\Delta R_{\\rm Gal-\\gamma}$ between the selected ON and OFF fields using a standard model of interstellar gas and cosmic-ray distributions; if the predicted diffuse background difference is comparable to or larger than the 4.8% systematic uncertainty or the observed residual, then the derived lifetime limits would be biased.","supporting_citations":[{"cited_title":"Optimized analysis method for indirect dark matter searches with Imaging Air Cherenkov Telescopes","cited_arxiv_id":"1209.5589","evidence_quote":"Supplies the full likelihood formalism used to convert the ON-OFF counts into a dark-matter lifetime limit."},{"cited_title":"Detecting Gamma-Ray Anisotropies from Decaying Dark Matter: Prospects for Fermi LAT","cited_arxiv_id":"0909.3514","evidence_quote":"Provides the J-factor and decay gamma-ray flux calculations for the NFW, Einasto, and isothermal halo profiles."},{"cited_title":"Acciari et al., Constraining Dark Matter lifetime with a deep gamma-ray survey of the Perseus Galaxy Cluster with MAGIC, Phys","cited_arxiv_id":null,"evidence_quote":"The 200-hour Perseus cluster observation that serves as the benchmark sensitivity comparison."},{"cited_title":"Constraints on the Galactic Halo Dark Matter from Fermi-LAT Diffuse Measurements","cited_arxiv_id":"1205.6474","evidence_quote":"The Fermi-LAT inner-halo lifetime limits shown for comparison at masses below 10 TeV."}],"review_version":1}