{"id":"781c8669-e2f9-4dc5-a18d-176d902874d9","arxiv_id":"1908.08884","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"No gamma-ray excess was found toward 31 dwarf irregular galaxies with HAWC, yielding 95% confidence upper limits on dark matter annihilation cross-sections and decay lifetimes in the 1 to 100 TeV mass range.","lead":"HAWC searched for gamma rays from 31 dwarf irregular galaxies and found no excess, so it set new limits on dark matter annihilation and decay for particle masses between 1 and 100 TeV. This shows that these small galaxies can serve as complementary targets for indirect dark matter searches, with limits comparable to classical dwarf spheroidals.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"J-factor systematics are not propagated in the HAWC dIrr limits, and the authors explicitly state no detailed kinematics were used; the headline numerical limits at 10 TeV are directly dependent on these unquantified factors.","rationale":"The paper's central claim is a set of exclusion limits, so the load-bearing input is the normalization of the expected dark-matter signal in each of the 31 targets. The J-/D-factors are computed by assuming Burkert halos with parameters from the Karukes & Salucci universal rotation curve, and the authors themselves state in the Conclusions that the analysis 'does not include detailed kinematics to compute J-factors and is therefore limited.' Because the limits scale almost directly with these factors, and because no systematic uncertainty on them is propagated, the numerical values in the abstract and Section 3 are not yet fully secure. I do not regard the Burkert/universal-rotation-curve choice as illegitimate; the issue is treating the resulting factors as exact constants while the paper simultaneously flags the absence of kinematics. The reader's weakest assumption is the same, so I agree with the conditional verdict. The remaining issues—the inconsistent 'sample of 30' in the Conclusions, the likelihood formula in Eq. (2.1) having a product of logarithms and an incorrectly signed exponent, and the mislabeled decay panels in the Fig. 3/4 captions—are less load-bearing than the J-factor systematics, but they reinforce the need for a revised version. The HAWC likelihood framework (Liff) is standard and the null result is plausible, so I see no reason to reject or to demand new data; the condition is a rigorous J-factor sensitivity study and correction of the presentation issues.","tokens_in":4576,"tokens_out":9700,"duration_ms":109314,"concrete_test":"Recompute the combined 95% CL limit at 10 TeV after replacing each fixed log10 J_i and log10 D_i in Table 1 with a nuisance parameter drawn from a log-normal prior with sigma = 0.3 dex, profiled over in the joint likelihood, and compare with the quoted 3.31e-22 cm^3/s and 6.16e25 s. Repeat the test with a coherent +0.3 and -0.3 dex shift of all J- and D-factors. If either shift moves the quoted limits by more than ~50%, the numerical constraints are dominated by unquantified J-factor assumptions and the conclusions must be revised to present limits as a band over the J-factor systematic.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The quantitative result—95% CL limits of 3.31e-22 cm^3/s (annihilation to tau+tau-) and 6.16e25 s (decay to tau+tau-) at 10 TeV—is controlled by the J- and D-factors of Table 1, because the signal expectation in every bin is proportional to these factors and the limits are approximately inversely proportional to them. The paper computes these factors with Clumpy using Burkert profiles whose parameters are taken from the Karukes & Salucci universal rotation curve, not from individual kinematic fits. The conclusions explicitly concede: 'This analysis does not include detailed kinematics to compute J-factors and is therefore limited.' No J-/D-factor uncertainty is propagated into the combined likelihood; dSph analyses in the same energy range typically profile over log-normal J-factor errors with sigma ~0.2-0.5 dex. The highest-weight galaxies (WLM, NGC 6822, IC 10, DDO 154) dominate the combined limit, so a factor of ~2 error in their normalization—plausible from distance, inclination, and Burkert-parameter scatter—shifts the quoted limits by a comparable factor. This is not an objection to the universal rotation curve as a model; it is an objection to treating derived J-factors as exact constants. The claim that dIrr limits are competitive with dSph limits should therefore be read as conditional on a J-factor normalization that is not yet demonstrated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript, an ICRC 2019 proceedings contribution from the HAWC Collaboration, reports a search for gamma-ray emission from dark matter annihilation and decay in 31 dwarf irregular (dIrr) galaxies within the HAWC field of view. The authors compute the annihilation J-factors and decay D-factors with CLUMPY, assuming Burkert dark matter profiles with parameters taken from the Karukes and Salucci universal rotation curve, and combine Poisson maximum-likelihood analyses over 1017 HAWC transits to set 95% CL limits on the annihilation cross-section and decay lifetime for WIMP masses between 1 and 100 TeV in the b bbar and tau+tau- channels. No significant excess is found. The headline limits at 10 TeV for the tau+tau- channel are 3.31e-22 cm^3/s for annihilation and 6.16e25 s for decay, and the combined limit is reported to improve on the best single-galaxy limit (DDO 154) above 10 TeV.","tokens_in":4900,"tokens_out":4765,"duration_ms":51307,"significance":"The paper's value is in extending TeV-scale indirect dark matter constraints to a new class of targets, dIrr galaxies, and in showing that a combined analysis of many low-astrophysical-factor objects can reach limits comparable to classical dSph galaxies for annihilation and to both classical and ultrafaint dSphs for decay. The analysis is not circular: the limits come from a null search with astrophysical factors fixed from external kinematic modeling, and the sample and numerical inputs are presented transparently in Table 1. The procedure follows prior HAWC dSph analyses, and the null result is plausible. However, the quantitative claims rest on two unquantified pillars: the absolute normalization of the J/D factors and the statistical significance of the null result. If those are supplied and demonstrated to be robust, the limits would be a useful addition to the HAWC dark matter program.","major_comments":[{"comment":"The conclusion explicitly states that \"This analysis does not include detailed kinematics to compute J-factors and is therefore limited,\" yet the paper quotes absolute exclusion limits that are directly proportional to the inverse of the J- and D-factors in Table 1. No uncertainty on these factors is computed or propagated, unlike dSph analyses that typically profile over log-normal J-factor errors with sigma ~0.2-0.5 dex. Because the combined limit is dominated by a few galaxies (WLM, NGC 6822, IC 10, DDO 154), a factor of about two normalization error shifts the headline limits by a comparable factor. The authors should propagate J/D-factor systematics or, at minimum, provide a sensitivity scan showing how the limits change under a range of halo-profile and distance assumptions.","section":"Section 4 and Table 1"},{"comment":"The text repeatedly states that \"no significant excess\" was found and that the significance was converted into an exclusion limit, but no test-statistic values, p-values, or expected-limit bands are reported, either for individual galaxies or for the combined analysis. Without these numbers the reader cannot determine whether the observed limits are consistent with the background-only expectation or are dominated by statistical fluctuations. The paper should report the observed and expected 95% CL limits with the associated significance (e.g., TS values or p-values) for each channel and for the combined analysis.","section":"Sections 2 and 3"},{"comment":"Equation (2.1) writes LSignal = product over i of ln[ (1/Ni!) (Bi+Si)^Ni exp(-(Bi+Si)) ], which is the product of logarithms and is not a likelihood; if LSignal is intended to be the log-likelihood, it should be a sum over bins. Because the limit-setting procedure and the null-hypothesis comparison both depend on this quantity, the notation must be corrected and the actual test statistic used by Liff should be specified, for example TS = -2 ln(LNull/LSignal).","section":"Equation (2.1)"}],"minor_comments":[{"comment":"The caption of Figure 4 labels the panels as \"Annihilation to b quarks\" and \"Annihilation to tau leptons,\" but the figure shows decay limits; the captions should say \"Decay\" instead.","section":"Figure 4"},{"comment":"The conclusions state that the local volume contains \"four times bigger than the sample of 30 we used here,\" but Table 1 lists 31 galaxies; the count should be corrected for consistency.","section":"Section 4"},{"comment":"The phrase \"ultra very high energies\" is redundant and should be replaced with a standard energy-band term such as \"very high energies.\"","section":"Abstract"},{"comment":"The caption of Table 1 does not define the assumed halo model or the angular integration aperture used for the J- and D-factors; a sentence stating the Burkert profile parameters and integration angle should be added.","section":"Table 1"},{"comment":"References [2] and [3] are cited for the claim that gamma rays from star-forming regions in dIrr galaxies are in the GeV range, but neither reference is about dIrr galaxies specifically; a more directly relevant citation would strengthen the background-free argument.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"This is a conference proceedings-style paper from ICRC 2019, and the methodological depth is intentionally limited. For a journal publication, the main risk is not the likelihood framework itself but the astrophysical normalization: the quoted limits are directly controlled by J/D factors that the authors themselves acknowledge are not based on detailed kinematics. I would urge the editor to require at least a quantitative treatment of J-factor systematics and the reporting of significances and expected limits before considering the paper for publication. The topic and results are in scope for the journal and are likely of interest to the astroparticle community, but the current manuscript does not yet meet the evidentiary standard for a refereed journal article."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is the first HAWC analysis of dwarf irregular galaxies for dark matter annihilation and decay, and for that reason alone it is worth a read. The analysis is a straightforward extension of the established HAWC dSph likelihood framework to 31 dIrr targets, with a null result and combined 95% limits. At 10 TeV in the tau+tau- channel, the quoted annihilation limit is 3.31e-22 cm^3/s and the decay limit is 6.16e25 s. The combined limit beats the best single galaxy above 10 TeV, and the comparison with previous dSph limits is honest: dIrr constraints are comparable to classical dSphs, not to ultra-faints. The paper is clearly written, the likelihood construction follows prior HAWC analyses, and the authors flag the main limitation themselves in the conclusions: no detailed kinematics were used to compute the J-factors.\n\nThe soft spot is real and it is the J-factors. The limits scale almost inversely with the astrophysical factors in Table 1, and those factors come from Burkert profiles fixed by the Karukes & Salucci universal rotation curve, not from individual kinematic fits. No systematic uncertainty on the J- or D-factors is propagated into the likelihood. The four heaviest-weight galaxies (WLM, NGC 6822, IC 10, DDO 154) drive the combined limit, so a factor-of-two error in their normalization shifts the quoted numbers by a comparable factor. This is not an objection to the universal rotation curve as a model; it is an objection to treating the derived factors as exact constants when the paper itself concedes the kinematics are missing. The dSph literature in the same energy range typically profiles over log-normal J-factor errors; this analysis does not.\n\nThere are also minor internal inconsistencies: the abstract and Table 1 say 31 galaxies, the conclusions say 30, and the caption of Figure 4 says \"Annihilation\" where it should say \"Decay\". These are cosmetic, and I would not hold up the paper over them, but a referee should ask for a clean pass.\n\nWho is this for? Model builders who want TeV-scale limits on a new target population, and HAWC collaboration members who want to know the status of the dIrr program. The circularity burden is low: the gamma-ray data are not fitted to the kinematics, the null result is plausible, and the self-citation (Karukes & Salucci 2017) is legitimate as the source of the rotation curve. Serious thinker: yes. The paper knows what it can and cannot claim.\n\nI would recommend engaging with it: send it to a serious referee, mostly for the J-factor systematic, and accept it as a conference proceedings after the number inconsistency is fixed and the limitations are stated where the limits are quoted rather than only in the conclusions.","headline":"First HAWC dIrr-galaxy dark matter limits, from a clean null search, but the J-factor normalizations are unquantified and the paper admits it.","tokens_in":5399,"tokens_out":699,"would_cite":false,"duration_ms":8095,"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":"31 dwarf irregular galaxies set new TeV dark matter limits","keywords":["dark matter","indirect detection","dwarf irregular galaxies","gamma-ray astronomy","WIMP annihilation","decay lifetime","TeV scale","HAWC"],"falsifier":"Measure the actual dark-matter distribution of DDO 154, the galaxy that drives the combined limit, from resolved rotation-curve or stellar-kinematic data and recompute its $J$-factor; if the new $J$-factor is, say, a factor of two below the value used here, the quoted 10 TeV cross-section limit would shift upward by the same factor. Alternatively, a TeV gamma-ray detection from any of the 31 galaxies that follows star-formation tracers rather than the $J$-factor would falsify the background-free assumption that the limits rely on.","tokens_in":4440,"feed_emoji":"🌌","tokens_out":8377,"duration_ms":78821,"temperature":0.7,"pith_summary":"This work argues that dwarf irregular (dIrr) galaxies, which are gas-rich and usually set aside because their star-forming regions produce gamma rays, can serve as clean, background-free targets for indirect dark matter searches at TeV energies. Using data from a wide-field ground-based gamma-ray observatory at 31 dIrr positions, the analysis finds no significant excess and converts the null result into 95% confidence exclusion limits on WIMP annihilation and decay between 1 and 100 TeV. The combined 31-galaxy limit outperforms the best single galaxy above 10 TeV, with the annihilation cross-section to tau pairs capped at $3.31\\times10^{-22}\\,\\mathrm{cm^3\\,s^{-1}}$ at 10 TeV and the decay lifetime at $6.16\\times10^{25}\\,\\mathrm{s}$. This matters because stacking many dark-matter-dominated dwarfs is a route to stronger TeV constraints without relying on a few exceptional targets.","feed_headline":"31 dwarf galaxies tighten TeV dark-matter limits","feed_subtitle":"A stacked analysis of 31 gas-rich dwarfs rules out annihilation and decay signals across 1–100 TeV.","key_machinery":"The load-bearing machinery is the joint likelihood over all 31 galaxies, in which the annihilation cross-section or decay lifetime is a single shared parameter while each galaxy contributes its own expected signal built from its astrophysical factor. The astrophysical factors, $J$ for annihilation and $D$ for decay, are computed from cored Burkert dark-matter profiles whose parameters come from a universal rotation-curve relation for dwarf disc galaxies, integrated with standard halo-model code. At TeV energies the suspected GeV-range emission from star-forming regions is neglected, so each galaxy is modelled as signal plus background counts, with the background taken from off-source data. The 95% limits are obtained by comparing the joint signal model with a null hypothesis.","core_discovery":"The central claim is that dIrr galaxies, despite their gas and star formation, can be treated as background-free at energies above about 1 TeV, and that a joint-likelihood analysis of 31 such galaxies in a wide-field survey yields meaningful TeV-scale constraints on dark matter. No statistically significant gamma-ray excess is found at any of the 31 positions. For a 10 TeV WIMP annihilating to $\\tau^+\\tau^-$, the 95% exclusion on the velocity-weighted cross-section is $\\langle\\sigma v\\rangle = 3.31\\times10^{-22}\\,\\mathrm{cm^3\\,s^{-1}}$; for decaying dark matter to the same channel, the exclusion on the lifetime is $\\tau_\\chi = 6.16\\times10^{25}\\,\\mathrm{s}$. The combined limit equals or exceeds the best individual galaxy (DDO 154) below 10 TeV and improves on it above 10 TeV, because the gain from stacking is offset at lower energies by the large spread in signal-to-background among the galaxies.","pith_inferences":["If the per-galaxy dark-matter content were recalibrated with detailed kinematic data, the $J$-factor systematics now dominating the limits could be reduced, potentially making the 31-galaxy stacking competitive with the best dwarf spheroidal targets for decay channels.","The background-free assumption above 1 TeV is a testable prediction: a future observation of gamma rays from a dIrr galaxy that tracks star-formation activity rather than dark-matter density would falsify it and force a reanalysis with diffuse emission models.","The crossover behaviour of the combined limit near 10 TeV suggests an optimization: weighting galaxies by their signal-to-background ratio rather than by $J$-factor alone might extend the stacking gain to lower masses.","The same joint-likelihood recipe could be transferred to other wide-field TeV instruments, where larger samples of gas-rich dwarfs would yield lifetime constraints complementary to those from dwarf spheroidals."],"forward_implications":["Dwarf irregulars can be added to dwarf spheroidal galaxies as a TeV dark-matter target population, roughly doubling the number of stacked targets available to wide-field observatories.","For decaying dark matter, the dIrr limits are comparable to those from classical and ultrafaint dwarf spheroidals, so a combined dIrr-plus-dSph lifetime constraint at TeV masses is feasible.","Because the local volume contains roughly four times more dIrr galaxies than the 31 used here, an enlarged sample should sharpen the combined annihilation and decay limits.","Above 10 TeV, stacking the full population beats the single best target, so future TeV analyses should prefer population stacking over selecting only the highest-$J$ galaxy."],"supporting_citations":[{"why":"supplies the universal rotation curve and dark-matter halo parameters from which the dIrr $J$- and $D$-factors are derived.","marker":"[1]"},{"why":"establishes that gamma-ray emission from star-forming galaxies is concentrated at GeV energies, supporting the background-free treatment above 1 TeV.","marker":"[2]"},{"why":"provides context for TeV emission from star-forming regions, used to argue that dIrr gamma-ray flux is negligible above 1 TeV.","marker":"[3]"},{"why":"is the code used to compute the annihilation and decay astrophysical factors for the cored halo profiles.","marker":"[4]"},{"why":"is used to compute the photon spectra for the $b\\bar{b}$ and $\\tau^+\\tau^-$ annihilation and decay channels.","marker":"[5]"},{"why":"is the earlier dwarf spheroidal analysis whose likelihood and analysis framework this work extends to dIrr galaxies.","marker":"[6]"},{"why":"is the high-level analysis framework used to compute significances and exclusion limits.","marker":"[7]"}],"fun_headline_variants":["HAWC stacks 31 dwarfs to bound dark-matter TeV signals","No gamma excess in 31 dwarfs: new dark-matter limits","Dark-matter constraints from 31 dwarf irregular galaxies","TeV dark-matter limits tightened by HAWC dwarf survey","31 gas-rich dwarfs yield tight TeV dark-matter bounds"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The analysis assumes that the dark-matter content of each of the 31 galaxies, computed from a standard cored halo profile with parameters taken from a universal rotation-curve relation rather than from each galaxy's own detailed kinematics, is accurate; if these inferred dark-matter amounts are too large, the quoted exclusion limits are too strong.","fun_headline_variants_meta":{"raw":{"variants":["HAWC stacks 31 dwarfs to bound dark-matter TeV signals","No gamma excess in 31 dwarfs: new dark-matter limits","Dark-matter constraints from 31 dwarf irregular galaxies","TeV dark-matter limits tightened by HAWC dwarf survey","31 gas-rich dwarfs yield tight TeV dark-matter bounds"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000146,"raw_usage":{"total_tokens":1153,"prompt_tokens":887,"completion_tokens":266,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":503,"completion_tokens_details":{"reasoning_tokens":176}},"tokens_in":503,"tokens_out":266,"duration_ms":3068,"temperature":1.0,"reasoning_tokens":176,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:26:07.834430+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the actual dark-matter distribution of DDO 154, the galaxy that drives the combined limit, from resolved rotation-curve or stellar-kinematic data and recompute its $J$-factor; if the new $J$-factor is, say, a factor of two below the value used here, the quoted 10 TeV cross-section limit would shift upward by the same factor. Alternatively, a TeV gamma-ray detection from any of the 31 galaxies that follows star-formation tracers rather than the $J$-factor would falsify the background-free assumption that the limits rely on.","supporting_citations":[{"cited_title":"Martin, Interstellar gamma-ray emission from cosmic rays in star forming galaxies, A&A 564, A61 (2014) [DOI: 10.1051/004-6361/201323329 ]","cited_arxiv_id":null,"evidence_quote":"establishes that gamma-ray emission from star-forming galaxies is concentrated at GeV energies, supporting the background-free treatment above 1 TeV."},{"cited_title":"H.E.S.S. observations of massive stellar clusters","cited_arxiv_id":"0906.2637","evidence_quote":"provides context for TeV emission from star-forming regions, used to argue that dIrr gamma-ray flux is negligible above 1 TeV."},{"cited_title":"Sjostrand, S","cited_arxiv_id":null,"evidence_quote":"is used to compute the photon spectra for the $b\\bar{b}$ and $\\tau^+\\tau^-$ annihilation and decay channels."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"is the high-level analysis framework used to compute significances and exclusion limits."}],"review_version":1}