{"id":"04cb9b44-1a2e-4139-a275-b262e3d801b6","arxiv_id":"1908.10178","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"First IACT observations of dwarf irregular galaxy WLM find no gamma-ray excess and set 95% CL upper limits on dark matter annihilation cross section, reaching about 10^-21 cm^3/s for the tau+tau- channel at 1 TeV.","lead":"H.E.S.S. observed the dwarf irregular galaxy WLM for 19 hours, found no gamma-ray excess, and set the first imaging Cherenkov telescope upper limits on dark matter annihilation. The result adds a new target class to indirect dark matter searches.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (5.4) defines Junc as the maximizer of the J prior, making the claimed J-factor nuisance correction a tautology; the reported limits may not actually include J systematics as stated.","rationale":"The reader's weakest_assumption focused on the external reliability of the WLM J-factor from the coreNFW fit. My concern is closely related but more specific and internal: the paper claims to fold the J uncertainty into the upper limits via Eq. (5.4), but as written that equation is a tautology because the maximizing value of the log-normal in Eq. (5.3) is just the mean. This is not an attack on the null result: the absence of a significant excess is well supported by the event counts and significance. It is, however, a genuine gap in the derivation of the numerical upper limits, which are the central claim. If the actual H.E.S.S. analysis used a profile likelihood over J, the text does not describe it correctly, and the reported limits cannot be independently checked from the paper. Because the paper is a proceedings and the key equations are short, this should be fixable with a clarification or a corrected formula. That warrants a conditional acceptance rather than a rejection: the underlying experimental result is likely sound, but the published systematic treatment needs correction or explicit justification.","tokens_in":6131,"tokens_out":12507,"duration_ms":132740,"concrete_test":"Re-derive Eq. (5.4) from Eqs. (5.1)-(5.3) by profiling the full likelihood with respect to J. First, maximize Eq. (5.3) to verify that Junc = Jbar, making Eq. (5.4) an identity. Then recompute the 95% CL limits from Eq. (5.1) with J profiled over the log-normal term, using the ON/OFF counts in Table 1. If the profiled limits differ from the no-J limits by more than a few percent, the published J-systematic treatment is incorrect; if they do not differ, the 'conservative' claim in Sec. 6 should still be corrected or clarified.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central quantitative result is the set of 95% CL upper limits in Sec. 6, but their derivation in Sec. 5 is internally inconsistent. Equation (5.3) is a log-normal distribution in log10 J, so its maximum is at log10 J = log10 Jbar, i.e. Junc = Jbar. Inserting this into Eq. (5.4) gives <sigma_v>_95 = <sigma_v>_0, meaning the stated 'nuisance parameter on J' changes nothing. This contradicts the claim in Sec. 6 that the J uncertainties are included and make the limits more conservative. If the authors intended Junc to be a profile-likelihood value from Eq. (5.1), that is not what Eq. (5.3) describes; moreover, the standard profile-likelihood result does not reduce to a fixed rescaling by the mode of the prior. Because the upper limits are the paper's main deliverable, the treatment of the dominant astrophysical normalization is either absent or at least not transparently derived. The numerical effect may be modest because sigma_log10J = 0.037, but as written the text cannot justify the quoted limits' systematic coverage.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports the first imaging atmospheric Cherenkov telescope observations of the dwarf irregular galaxy WLM, carried out with H.E.S.S. for about 19 hours of live time in 2018. The analysis searches for a gamma-ray excess from dark matter annihilation, finds no significant excess (0.7 sigma), and derives 95% confidence level upper limits on the velocity-weighted annihilation cross section for the b-bbar, tau+tau-, W+W-, and Z+Z- channels. The limits reach about 1e-20 cm^3/s for the quark and boson channels and about 1e-21 cm^3/s for the tau+tau- channel at a dark matter mass of 1 TeV, improving on HAWC limits by an order of magnitude or more. The paper claims that the uncertainty on the J-factor is included as a log-normal nuisance parameter.","tokens_in":6364,"tokens_out":10236,"duration_ms":104021,"significance":"If the analysis is technically sound, this is a valuable first constraint on dark matter annihilation from a dwarf irregular galaxy, a new target class for ground-based gamma-ray indirect searches. The target selection is well motivated: WLM has a relatively high J-factor with a very small reported uncertainty (0.037 dex in log10 J), so the derived limits have a strong astrophysical normalization. The work uses standard H.E.S.S. multiple-OFF background estimation and a likelihood-ratio framework, and provides observed and expected limits with containment bands. The explicit statement of the non-detection and the comparison with HAWC limits are clear and useful. However, the treatment of the J-factor uncertainty as presented in Eqs. (5.3) and (5.4) is internally inconsistent and does not implement the claimed nuisance-parameter procedure, so the systematic coverage of the central result is not currently transparent.","major_comments":[{"comment":"The definition of J_unc as the value of J that maximizes Eq. (5.3) makes the J-factor correction nearly a tautology. Eq. (5.3) is a log-normal density in log10 J, so its maximum is at the mode, J_unc ≈ Jbar × 10^{-σ^2 ln 10}, which for σ = 0.037 differs from Jbar by less than 1%. Inserting this into Eq. (5.4) therefore rescales the limits by about 1%, which cannot account for the J-factor uncertainty in the way claimed in Section 6 ('the uncertainties on J as a nuisance parameter ... makes the derivation of the upper limits more conservative'). A correct profile-likelihood treatment would broaden the limits by a factor that depends on σ (roughly exp(k σ) for a one-sided 95% interval); the present equation does not do so. Please either present the full profile-likelihood calculation or revise the text to remove the claim that the J uncertainty is included in the quoted limits.","section":"Section 5, Eqs. (5.3)–(5.4) and Section 6"},{"comment":"The analysis inherits the J-factor log10 J = 16.6 ± 0.037 from the coreNFW fit of Ref. [5], which assumes WLM is dark-matter-dominated at all radii and uses a specific density profile. Because the limits scale directly with J via Eq. (5.4), a different dark matter profile assumption (e.g., a cuspier profile) or the effects of tidal stripping would shift the limits by the corresponding factor. The paper does not discuss or propagate this systematic uncertainty. At minimum, a sentence acknowledging this limitation and quoting the resulting range of limits would be needed to support the claim of a small systematic error.","section":"Section 3 and Eq. (5.4)"}],"minor_comments":[{"comment":"The definition of α is ambiguous: the text says α renormalizes the OFF region to the ON region, and Eq. (5.2) uses α N_Bi as the mean of the OFF Poisson, but the reported value α = 16.24 would imply the OFF region has 16.24 times the exposure of the ON region. Please clarify the definition (ON/OFF versus OFF/ON) and make the formulas and table consistent.","section":"Section 4, Table 1 and Eq. (5.2)"},{"comment":"The normalization of the log-normal density is missing a factor of ln 10 if σ is the standard deviation in log10 J; as written, the integral of Eq. (5.3) over J is not equal to unity.","section":"Section 5, Eq. (5.3)"},{"comment":"Reference [4] is malformed; it should read 'J. I. Read, G. Iorio, O. Agertz, et al.' rather than 'J. Read and I. Iorio, G., Agertz, et al.'.","section":"References"},{"comment":"The text states that the J-factor uncertainty is 'log10 J = 16.6 ± 0.037' but does not specify whether 0.037 is the standard deviation of the Gaussian fit in log10 J shown in Fig. 1. Please make this explicit.","section":"Section 3 and Figure 1"},{"comment":"The expected limits and 1σ/2σ bands are derived from only 100 Poisson realizations; with such a small sample, the containment bands will have noticeable statistical fluctuations. A larger set of realizations would give smoother bands.","section":"Section 6"},{"comment":"The phrase 'at 1 Mpc from the Milky Way' is imprecise; WLM is a galaxy at a distance of about 1 Mpc from the Milky Way (or from Earth), and the sentence should be reworded.","section":"Section 2"},{"comment":"In the right panel, the label 'log10 = 0.037' is missing a subscript J; it should read 'log10 J = 0.037'.","section":"Figure 1 caption"}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern about Eq. (5.4) is valid and lands on a load-bearing point of the paper. The authors should be asked to either provide the actual profile-likelihood calculation used to set the upper limits or explicitly state that the J-factor uncertainty has not been included. The numerical impact of the error is likely small because the quoted σ is only 0.037 dex, but the methodological inconsistency cannot be left as it stands in a journal publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Plain take: this is the first IACT constraint on dark matter annihilation from a dwarf irregular galaxy, and it's a clean null result from 19 hours of H.E.S.S. data on WLM. The limits are about a factor of 10–100 better than HAWC's for the same target class, but the way the J-factor uncertainty is claimed to be included in the limits is wrong. The numerical impact is small because log10 J has a quoted error of only 0.037 dex, so the final numbers are probably okay—but the text overstates what Eq. (5.4) actually does.\n\nWhat's new: WLM as a target, CT5 mono observations, no significant excess (0.7 sigma), and 95% CL upper limits on <sigma v> for bb, tau+tau, W+W, and Z+Z channels. The analysis uses standard H.E.S.S. multiple-OFF background estimation and a Poisson likelihood; expected limits and containment bands come from Poisson realizations. The paper is clearly written and the null result is well supported.\n\nThe soft spot is Sec. 5. Eq. (5.3) is a log-normal pdf on J. Its maximum is at Junc ~ Jbar (up to a tiny factor from the 1/J prefactor). Plugging that into Eq. (5.4) gives a rescaling of order unity—not a profile-likelihood treatment of a nuisance parameter. To actually include J systematics you would need to maximize the full likelihood in Eq. (5.1) over J for each test point, or marginalize over J. The paper says the uncertainties are included and make the limits more conservative; the equations don't support that. The effect is probably a few percent, so this is a presentation flaw rather than a decisively wrong result, but it should be fixed.\n\nA minor related point: the external J-factor from coreNFW fits has a very small quoted uncertainty. For a single galaxy that looks optimistic, and the paper doesn't discuss model dependence (core vs. cusp, tidal stripping). A caveat sentence would help.\n\nThis paper is for the indirect-dark-matter community, especially those planning CTA dwarf campaigns. It deserves a serious referee: the dataset is real, the analysis is mostly standard, and establishing a new target class is worth publishing. But the J-factor treatment needs rewriting before the limits can be taken as including systematic uncertainties. I'd recommend refereeing it, with a request to re-derive the limits with a proper profile likelihood or to clearly state that J systematics are not included.","headline":"First IACT search for DM in a dwarf irregular galaxy, with a clean null result and modest limits—but the J-factor nuisance treatment in Sec. 5 is not what it claims to be.","tokens_in":6888,"tokens_out":6200,"would_cite":false,"duration_ms":63613,"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":"After 19 hours of H.E.S.S. observations toward the dwarf irregular galaxy WLM, no significant gamma-ray excess is found, and the resulting 95% confidence upper limits on dark-matter annihilation reach about $10^{-20}$ cm$^3$ s$^{-1}$ for…","keywords":["dark matter annihilation","dwarf irregular galaxies","WLM","gamma-ray astronomy","imaging atmospheric Cherenkov telescope","J-factor","upper limits","indirect dark matter detection"],"falsifier":"Re-fit WLM's stellar and gas kinematics without the coreNFW assumption—for example with a cuspy NFW or tidally truncated profile—and recompute $\\log_{10}J$; if the result moves outside $16.6\\pm0.037$ by an amount larger than the quoted uncertainty, the reported limits scale proportionally and the claim that WLM has a precisely known J-factor fails.","tokens_in":5959,"feed_emoji":"🔭","tokens_out":12119,"duration_ms":106632,"temperature":0.7,"pith_summary":"This paper reports the first search for dark-matter annihilation gamma rays from a dwarf irregular galaxy with an imaging atmospheric Cherenkov telescope. H.E.S.S. observed WLM for 19 hours in 2018 and found no significant excess in the direction of the galaxy, with a signal significance of 0.7 $\\sigma$. From the null result, the paper derives 95% confidence upper limits on the velocity-weighted annihilation cross section for the $b\\bar{b}$, $W^+W^-$, and $Z^+Z^-$ channels of about $10^{-20}$ cm$^3$ s$^{-1}$, and about $10^{-21}$ cm$^3$ s$^{-1}$ for the $\\tau^+\\tau^-$ channel at a dark-matter mass of 1 TeV. These are the first IACT constraints on this new class of targets, and they improve on the previous limits by a factor of 10 to nearly 100. Dwarf irregulars matter because they are rotationally supported and dark-matter dominated at all radii, so their J-factors can be measured from kinematics with much smaller uncertainties than those of pressure-supported dwarfs.","feed_headline":"No dark-matter signal in 19 hours of WLM observations","feed_subtitle":"The first Cherenkov-telescope study of a dwarf irregular galaxy sets 95% limits near 10⁻²⁰ cm³/s at 1 TeV.","key_machinery":"The argument runs on three pieces. First, the predicted gamma-ray flux from dark-matter annihilation is factorized as $\\frac{1}{2}\\frac{\\langle\\sigma v\\rangle}{4\\pi m_\\chi^2}\\frac{d\\Phi^{\\mathrm{PP}}}{dE}J$, where $J$ is the line-of-sight integral of the dark-matter density squared over the 0.1-degree region of interest. Second, WLM's dark-matter distribution is described by a coreNFW profile, defined as the NFW profile modified by a factor $f^n(r)$ that flattens the inner density below a core radius $r_c$; an MCMC fit to stellar kinematics yields $\\log_{10}J = 16.6\\pm0.037$, giving a precisely known normalization. Third, the statistical machinery is a Poisson likelihood for ON/OFF counts with a log-normal likelihood for $J$, and the 95% limits are obtained from the profile likelihood through the asymptotic condition $\\mathrm{TS}=2.71$, using the scaling $\\langle\\sigma v\\rangle_{95\\%} = \\langle\\sigma v\\rangle_0 \\bar{J}/J_{\\mathrm{unc}}$ to absorb the J-factor uncertainty.","core_discovery":"The paper's central result is a null detection: in 18.6 hours of live time toward WLM, the ON region contains 1677 events against an OFF-normalized background expectation, corresponding to a gamma-ray excess of 31.2 events with a significance of 0.7 $\\sigma$. Treating the J-factor as a log-normal nuisance parameter with $\\log_{10} J(\\mathrm{GeV}^2\\,\\mathrm{cm}^{-5}) = 16.6 \\pm 0.037$ from the coreNFW fit, the analysis sets 95% CL upper limits on $\\langle\\sigma v\\rangle$ by a log-likelihood ratio test. At 1 TeV the observed limits reach $10^{-20}$ cm$^3$ s$^{-1}$ for $b\\bar{b}$, $W^+W^-$, and $Z^+Z^-$ and $10^{-21}$ cm$^3$ s$^{-1}$ for $\\tau^+\\tau^-$, improving on the previous limits by a factor of 10 to almost 100. The authors conclude that dwarf irregular galaxies, observed for the first time with an IACT, are viable targets for indirect dark-matter searches.","pith_inferences":["The same analysis chain could be applied to the other 35 known dwarf irregulars within 11 Mpc; a stacked analysis would improve sensitivity roughly as the square root of the number of galaxies, assuming comparable J-factor precision.","If future stellar-kinematic data confirm the small J-factor uncertainty, dwarf irregulars could serve as a cross-check on dwarf spheroidal results, whose dark-matter distributions typically carry much larger systematic uncertainties.","The non-detection at the $10^{-21}$ cm$^3$ s$^{-1}$ level in the leptonic channel remains above the canonical thermal-relic cross section, so the near-term payoff is methodological—demonstrating a low-background, well-calibrated target class—rather than a constraint on WIMP particle properties."],"forward_implications":["Dwarf irregular galaxies are now demonstrated to be viable dark-matter targets for imaging atmospheric Cherenkov telescopes, opening a new class of sources for the technique.","The 95% confidence upper limits on $\\langle\\sigma v\\rangle$ for WLM improve on the previous irregular-galaxy limits by a factor of 10 to almost 100 across the four annihilation channels.","The $\\tau^+\\tau^-$ channel is the most constraining, reaching about $10^{-21}$ cm$^3$ s$^{-1}$ at a dark-matter mass of 1 TeV.","Because the reported limits scale inversely with the J-factor, any future revision of WLM's dark-matter distribution directly rescales the constraints via the paper's Eq. (5.4)."],"supporting_citations":[{"why":"Establishes the baseline limits from prior irregular-galaxy observations that the WLM results improve on.","marker":"[1]"},{"why":"Supplies the per-channel gamma-ray spectra from Pythia used to compute the expected signal.","marker":"[2]"},{"why":"Describes the multiple-OFF background estimation method used to derive the ON/OFF counts.","marker":"[3]"},{"why":"Establishes WLM's rotation curve and dark-matter dominance used to justify the target and its mass model.","marker":"[4]"},{"why":"Provides the coreNFW dark-matter profile and the MCMC fit from which the J-factor is taken.","marker":"[5]"},{"why":"Provides the asymptotic profile-likelihood method used to set 95% CL upper limits.","marker":"[6]"},{"why":"Gives the J-factor nuisance scaling property used to include J uncertainty in the limits.","marker":"[7]"}],"fun_headline_variants":["No dark matter signal in first Cherenkov study of WLM dwarf galaxy","H.E.S.S. finds no dark matter in WLM, first for dwarf irregulars","WLM null result tightens dark matter limits from an IACT","No dark matter excess in WLM: first IACT limits for dwarf irregular","Dark matter bounds from WLM: no signal, limits improved 10-100x"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The limits rest on the assumption that WLM's dark-matter distribution is correctly given by the coreNFW profile fitted to its kinematics, with the quoted tiny uncertainty; if the inner profile is cuspier, tidally stripped, or otherwise different, the reported cross-section limits rescale by the same factor.","fun_headline_variants_meta":{"raw":{"variants":["No dark matter signal in first Cherenkov study of WLM dwarf galaxy","H.E.S.S. finds no dark matter in WLM, first for dwarf irregulars","WLM null result tightens dark matter limits from an IACT","No dark matter excess in WLM: first IACT limits for dwarf irregular","Dark matter bounds from WLM: no signal, limits improved 10-100x"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001593,"raw_usage":{"total_tokens":6366,"prompt_tokens":980,"completion_tokens":5386,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":596,"completion_tokens_details":{"reasoning_tokens":5282}},"tokens_in":596,"tokens_out":5386,"duration_ms":33879,"temperature":1.0,"reasoning_tokens":5282,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T10:49:27.579135+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-fit WLM's stellar and gas kinematics without the coreNFW assumption—for example with a cuspy NFW or tidally truncated profile—and recompute $\\log_{10}J$; if the result moves outside $16.6\\pm0.037$ by an amount larger than the quoted uncertainty, the reported limits scale proportionally and the claim that WLM has a precisely known J-factor fails.","supporting_citations":[{"cited_title":"Searching for TeV DM evidence from Dwarf Irregular Galaxies with the HAWC Observatory","cited_arxiv_id":"1708.04642","evidence_quote":"Establishes the baseline limits from prior irregular-galaxy observations that the WLM results improve on."},{"cited_title":"I., Iorio, G., Agertz, O., et al.\\ 2016, , 462, 3628","cited_arxiv_id":null,"evidence_quote":"Provides the coreNFW dark-matter profile and the MCMC fit from which the J-factor is taken."}],"review_version":1}