{"id":"0ca46c85-ed82-4131-95b0-de939f9b7102","arxiv_id":"1908.07255","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"IceCube's updated Earth-capture dark matter analysis projects up to a 3.8x sensitivity gain over its previous result at 1 TeV, with final limits on the real data set still to come.","lead":"The IceCube neutrino telescope reports projected sensitivities, not detections, for an eight-year search for dark matter that accumulates and annihilates inside Earth. The new event selection could improve the previous IceCube sensitivity by about 3.8 times at a 1 TeV dark matter mass.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Projected 3.8x sensitivity gain is computed from MC background with no systematics and an optimized BDT cut on the same MC; it is a best-case ceiling, not a robust expected sensitivity.","rationale":"The paper is a conference proceedings reporting projected sensitivities for an ongoing analysis, not measured limits. The methodology is standard: binned shape likelihood, Feldman-Cousins, and median pseudo-experiment sensitivity. The livetime increase alone accounts for a factor sqrt(8) ≈ 2.8; the BDT selection provides the remaining ~1.36. The main risk to the central claim is that the sensitivity is computed entirely from MC background simulations with no systematic uncertainties. The paper's own captions (Figs. 2 and 3) state 'only statistical uncertainties are included,' and Section 4 says only 10% of data validates data/MC agreement in the signal region. The BDT cut is optimized on the same MC used to evaluate the sensitivity, introducing a potential optimistic bias. These limitations do not invalidate the paper as a status report, but they mean the 3.8x factor should be interpreted as an idealized ceiling. The reader's conditional verdict is appropriate. A concrete systematic-uncertainty study could settle whether the projected improvement survives realistic background mismodeling.","tokens_in":4503,"tokens_out":12190,"duration_ms":124053,"concrete_test":"Recompute the median 90% upper limit with systematic nuisances: perturb the background rate by ±20% (Gaussian) and tilt the background zenith PDF by a 10% linear slope across the up-going region, regenerating pseudo-experiments from Eqs. 5.1-5.3. Additionally, re-evaluate the sensitivity at BDT cuts of 0.25 and 0.35 instead of the optimized 0.30, and, if possible, on an MC sample not used in the optimization scan. If the median limit worsens by more than ~30% or fluctuates strongly with the cut choice, the ~3.8x improvement claim is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim, a ~3.8x improvement in volumetric flux sensitivity over the previous IceCube analysis at mχ = 1 TeV (Section 5), is a median upper limit from 10^4 background-only pseudo-experiments. This estimate is only as valid as the simulated background defining both the total post-cut rate (0.19 mHz, Ntot = 43032) and the zenith-angle PDF B(bini) used in Eq. 5.2. The captions of Figs. 2 and 3 explicitly state that only statistical uncertainties are included, and Section 4 uses only a 10% data subset to validate data/MC agreement in the signal region. Because the BDT score cut was optimized for best sensitivity on this same MC (Section 5: 'A preliminary scan ... identified the region around 0.3 as the optimum cut'), the quoted median limit can be biased upward by unmodeled systematics in the background rate or zenith-angle shape, as well as by selection bias from optimizing the cut on the evaluation metric. A mismodeled background PDF near the Earth-center direction (cos θ ≈ -1) directly degrades the shape-only likelihood of Eq. 5.1. The 3.8x factor is therefore a best-case ceiling, not a robust expected performance.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports an updated sensitivity study for a search for dark matter annihilation in the center of the Earth using 8 years of IceCube data (2011-2018). The analysis employs a BDT-based event selection trained on two benchmark signals (mχ=50 GeV, τ+τ−; mχ=1 TeV, W+W−), a binned likelihood in the reconstructed zenith angle, and the Feldman-Cousins method to derive a 90% sensitivity defined as the median upper limit from 10^4 background-only pseudo-experiments. The central result is a projected volumetric neutrino-flux sensitivity that exceeds the previous IceCube result by a factor of ~3.8 at mχ=1 TeV in the W+W− channel, together with improved sensitivity on the annihilation rate compared to previous IceCube and ANTARES limits. The paper also outlines a future extension to an event-wise unbinned likelihood including reconstructed energy.","tokens_in":4646,"tokens_out":4755,"duration_ms":50442,"significance":"If the projected improvement is robust, this analysis would place IceCube in a competitive position for spin-independent dark matter-nucleon scattering constraints from Earth-capture signals, a complementary probe to direct-detection and solar-capture searches. The work uses standard, publicly documented simulation tools (CORSIKA, NuGen/GENIE, WimpSim), a well-established statistical framework (binned likelihood, Feldman-Cousins), and a large background-only pseudo-experiment ensemble, which provides a reproducible basis for the sensitivity estimate. These strengths are, however, partly offset by the issues identified in the major comments, so the quantitative significance of the claimed gain is currently not fully established.","major_comments":[{"comment":"The BDT score cut is set to 0.30 by a preliminary scan that optimizes the sensitivity on the same Monte Carlo used to compute the reported median limit. This in-sample optimization procedure introduces a selection bias that makes the quoted 3.8x improvement over the previous IceCube result a best-case ceiling rather than an unbiased expectation. The authors should either demonstrate that the sensitivity is stable under reasonable variations of the cut, use a validation sample that was not used for cut selection, or quantify the trial factor associated with the scan. Without this, the central claim is not fully robust.","section":"Section 5, central paragraph (after Fig. 4)"},{"comment":"The background estimate after the final BDT cut (0.19 mHz, Ntot=43032 events) and the zenith-angle background PDF B(bini) entering Eq. (5.1) are derived entirely from Monte Carlo (CORSIKA atmospheric muons, NuGen/GENIE neutrinos, and oscillation treatment), with the captions of Figs. 2 and 3 explicitly stating that only statistical uncertainties are included. Data/MC agreement in the signal region is validated with only a 10% data subset, which is then disregarded. Unmodeled systematics in the background rate or in the PDF shape, especially near the Earth-center direction (cos θ≈−1) where the signal peaks, directly propagate into the likelihood and would degrade the realized limits. The paper should propagate at least conservative systematic uncertainties on the normalization and shape of B(bini), or clearly label the reported sensitivity as a statistical-only upper bound.","section":"Section 4 and Figs. 2 and 3 captions"},{"comment":"The claimed improvement factor of ~3.8 compares a median expected sensitivity (from pseudo-experiments) with 'the previous IceCube sensitivity [2]'. If reference [2] actually reports an observed upper limit, as suggested by the phrase 'published upper limits on the spin-independent DM-nucleon scattering cross-section' in Sec. 1, then the comparison is not on equal footing: a median sensitivity should be compared with the previous analysis's median expected sensitivity, not with its observed limit. Please clarify the quantity taken from [2] and, if needed, recompute the improvement factor so that the comparison is apples-to-apples.","section":"Section 5, Fig. 5 and Sec. 1"}],"minor_comments":[{"comment":"The text states that DeepCore is installed at 'a depth of 1750 cm'; the intended unit is presumably meters (1750 m), consistent with the surrounding depths quoted in meters. Please correct this typo.","section":"Section 2, DeepCore description"},{"comment":"The sentence 'During the development of the event selection is was verified...' contains a typo ('is was' should be 'it was'), and the later sentence 'An initial series of cuts is applied on to reduce...' is grammatically awkward ('on to' should likely be 'to'). Please proofread these passages.","section":"Section 4, event selection paragraph"},{"comment":"The notation of the product index 'binmax ∏ bini=binmin' is nonstandard and slightly confusing; typically one writes ∏_{bini=binmin}^{binmax}. The authors may want to clarify that 'bini' is a bin label, not a running index.","section":"Section 5, Eq. (5.1)"},{"comment":"The conclusion says the new analysis 'could lead to world competitive limits' on the spin-independent cross-section, but the paper presents sensitivities, not observed limits. The wording should be adjusted to say that the analysis could produce world-competitive sensitivities or limits once applied to the full data.","section":"Section 6, Conclusion"},{"comment":"The paper would benefit from an explicit statement that all quoted results are median expected sensitivities and not observed limits, and from specifying how the 10% data validation subset is handled (the statement that it is 'disregarded for further analysis' is helpful but could be made more precise, e.g., that it is excluded from both the sensitivity estimate and any future limit calculation).","section":"General"}],"recommendation":"major_revision","confidential_remarks":"This is a proceedings paper from ICRC2019, so the review bar is appropriately lower than for a full journal article. The main scientific concern is that the headline sensitivity improvement is computed from Monte Carlo with only statistical uncertainties and a BDT cut tuned on the same Monte Carlo, which makes the quoted factor an optimistic ceiling. These issues are fixable within the scope of a proceedings paper by adding systematic uncertainties, demonstrating cut stability, and clarifying the comparison basis. I am not aware of any citation or novelty issues; the manuscript is a standard IceCube sensitivity study. I recommend major revision rather than rejection because the core methodology is standard and the data set is larger than in the previous analysis, so the qualitative conclusion of improved sensitivity is plausible; only the quantitative claim needs to be made more robust."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a conference proceedings paper from IceCube, not a final analysis. The headline number—a 3.8x better volumetric flux sensitivity at 1 TeV in the W+W- channel—is a median expected limit from 10^4 background-only pseudo-experiments. Treat it as a best-case ceiling, not a measured sensitivity. The paper is upfront about this: the analysis is labeled a sensitivity study, and the authors say they are moving to an unbinned energy-dependent likelihood.\n\nWhat it does well: the formalism is standard and cleanly presented (capture-annihilation equilibrium, flux conversion, Feldman-Cousins limits). The 8-year livetime and new BDT-based selection are a genuine extension over the 2017 IceCube result. The paper also gives credit where due: comparisons to previous IceCube and ANTARES limits, and explicit mention that only statistical uncertainties are included. For a status report, it is honest and useful.\n\nThe soft spots are real but not hidden. The BDT cut at 0.30 was chosen by scanning on the same MC that defines the sensitivity, which biases the quoted improvement upward. The background model is MC-only, with just a 10% data subset used to check data/MC agreement in the signal region; that is thin for a claim of factor-3.8. The stress-test note about the zenith-angle PDF being wrong near cos(theta) = -1 is a legitimate worry—the shape-only likelihood in Eq. 5.1 depends on it. These are exactly the caveats you'd expect from a proceedings paper, and the authors do not overstate the result.\n\nOne thing I'd push back on: the reader flags WimpSim being used both for signal PDFs and for flux-to-annihilation-rate conversion as circular. That's standard practice in this field—the conversion is a physics model, not an empirical fit—so I don't see it as a flaw here. The BDT optimization is the real issue, and it's a minor one because the paper calls it preliminary.\n\nWho is this for? Someone tracking IceCube's dark matter search program or wanting a reference for the current state of Earth-capture limits. It deserves a serious referee if submitted to a journal as a proceedings—not a desk reject—but the referee should insist on a clear statement that these are projected sensitivities, which the paper already makes. I'd engage with it as a useful status report, and I'd cite it for the sensitivity projection, but I wouldn't use the 3.8x number without the ceiling caveat.","headline":"Honest proceedings status report: the 3.8x sensitivity gain is a projected ceiling from MC with an optimized BDT cut, not a measured result, but the paper is transparent about that.","tokens_in":5262,"tokens_out":1574,"would_cite":true,"duration_ms":19773,"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 eight years of IceCube data, the projected sensitivity for dark-matter neutrinos from Earth's center improves by roughly 3.8.","keywords":["dark matter annihilation","Earth capture","IceCube","WIMP","neutrino telescope","boosted decision tree","sensitivity","spin-independent scattering"],"falsifier":"Unblind the full eight-year signal region and compare the observed post-cut event rate and zenith-angle distribution with the simulated 0.19 mHz background: any deviation beyond the quoted uncertainties would invalidate the projected median limit. Alternatively, compute the 90% upper limit directly from the observed data; if it is statistically consistent with the previous one-year IceCube limit rather than about 3.8 times more sensitive, the central claim is falsified.","tokens_in":4191,"feed_emoji":"🧊","tokens_out":12298,"duration_ms":114451,"temperature":0.7,"pith_summary":"Dark matter particles that scatter and lose energy inside Earth can accumulate at its center; annihilations there would produce a steady neutrino flux emerging from the planet's interior. This paper develops a search for that flux in eight years of IceCube data using a boosted decision tree, a machine-learning classifier, to separate signal-like neutrino events from atmospheric muons and neutrinos. The central result is a projected sensitivity, not a measured limit: for the $\\chi\\chi\\to W^+W^-$ channel at $m_\\chi=1$ TeV, the median expected 90% upper limit on the volumetric neutrino flux is about 3.8 times better than IceCube's previous one-year analysis, and the sensitivity on the annihilation rate improves on earlier IceCube and ANTARES limits. If the real data contain no signal, these sensitivities would translate into the strongest constraints yet from the Earth-capture channel on dark-matter-nucleon spin-independent scattering.","feed_headline":"IceCube data sharpen dark-matter search 3.8 times","feed_subtitle":"Machine-learning selection over 8 years would beat IceCube's old limit and rival ANTARES.","key_machinery":"The carrying object is the BDT classifier, which assigns every event a single signal-likeness score; one BDT is trained on a $m_\\chi=50$ GeV $\\tau^+\\tau^-$ benchmark and another on $m_\\chi=1$ TeV $W^+W^-$, with the high-mass model used for the reported sensitivities. After a cut on the score, a binned Poisson likelihood over reconstructed zenith angle compares data with a signal-plus-background model built from CORSIKA air-shower muons, NuGen/GENIE atmospheric neutrinos, and oscillation treatment below 100 GeV, while the signal is generated with WimpSim. The detector-level quantity that carries the result is the volumetric neutrino flux $\\Gamma_{\\nu\\to\\mu} = N_{\\rm sig}^{0.9}/(t_{\\rm live}V_{\\rm eff})$, which is converted through simulation to the annihilation rate $\\Gamma_A$, the physical variable compared with other experiments.","core_discovery":"On its own terms, the paper claims that a BDT-based event selection over eight years makes the Earth-center dark-matter channel a substantially more powerful neutrino search than the previous one-year IceCube analysis. For the benchmark channel $\\chi\\chi\\to W^+W^-$ with $m_\\chi = 1$ TeV, the median 90% sensitivity on the volumetric neutrino flux exceeds the earlier result by roughly a factor of 3.8, and the corresponding sensitivity on the annihilation rate is markedly better than both IceCube's previous upper limit and the ANTARES limit in the same channel. These numbers come from a binned likelihood in reconstructed zenith angle, a frequentist upper-limit prescription, and $10^4$ background-only pseudo-experiments, after a BDT score cut near 0.30 leaves a total background rate of 0.19 mHz. The paper presents this as the first sensitivity study of an ongoing analysis and expects further gains from moving to an event-wise unbinned likelihood that also uses reconstructed energy.","pith_inferences":["The quoted 3.8-fold gain is benchmarked at $m_\\chi=1$ TeV; a plausible extension, not stated in the paper, is that the gain varies with mass and may be smaller at low masses where the softer neutrino spectrum makes DeepCore veto performance the limiting factor.","Because the background estimate is validated on only 10% of the data, a natural robustness check is to compare the post-cut rate (0.19 mHz) and zenith distribution year by year; if the rate drifts, the median-limit projection degrades in a mass-dependent way.","The likelihood structure transfers directly to any fixed-direction source, so applying the same eight-year data set and BDT chain to the Sun would likely produce a comparable sensitivity gain for solar WIMP capture.","Adding reconstructed energy to the likelihood is a natural test of the method: the energy spectrum of secondary neutrinos differs sharply between $W^+W^-$ and $b\\bar b$ final states, so an energy-aware version of this analysis should separate channels more cleanly than the zenith-only version."],"forward_implications":["If the projected sensitivity becomes the observed result, IceCube's eight-year $W^+W^-$ search would surpass both its own one-year limit and ANTARES in that channel.","A null signal would convert the improved sensitivity into tighter upper limits on the spin-independent dark-matter-nucleon scattering cross-section, approaching the strongest constraints available from Earth-capture searches.","The same BDT and likelihood chain applies to other annihilation channels; the paper reports comparable selection performance for channels such as $b\\bar b$, so the method can be reused across final states.","Replacing the binned zenith-angle likelihood by the event-wise unbinned version with reconstructed energy is expected, as the paper states, to improve sensitivity further."],"supporting_citations":[{"why":"Supplies the previous one-year IceCube analysis whose sensitivity is the factor-3.8 comparison baseline, as well as the likelihood method being extended.","marker":"[2]"},{"why":"Gives the equilibrium formula relating the annihilation rate to the capture rate and annihilation cross-section, defining the target quantity of the search.","marker":"[4]"},{"why":"Provides the Earth capture rate as a function of WIMP mass that connects the flux sensitivity to the spin-independent scattering cross-section.","marker":"[7]"},{"why":"Simulates the atmospheric muon background from cosmic-ray air showers in the background estimate.","marker":"[8]"},{"why":"Supplies the treatment of neutrino oscillations inside Earth for energies below 100 GeV in the simulation.","marker":"[9]"},{"why":"WimpSim generates the neutrino energy spectra for each annihilation channel and mass used for signal simulation.","marker":"[10]"},{"why":"Provides the frequentist prescription used to convert the binned likelihood into a 90% upper limit on signal events.","marker":"[11]"},{"why":"ANTARES limits are the external comparison for the annihilation-rate sensitivity in the W+W- channel.","marker":"[12]"}],"fun_headline_variants":["IceCube 8-year data improves dark-matter sensitivity 3.8x","Dark matter in Earth: IceCube's 8-year hunt gains 3.8x","BDT selection lifts IceCube dark-matter sensitivity 3.8-fold","8 years of IceCube data sharpen Earth-core dark-matter search","IceCube's 8-year data improves Earth dark-matter probe 3.8x"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The projection assumes that the simulated atmospheric background reproduces the real detector data after the BDT cut well enough that the median expected limit from background-only pseudo-experiments is meaningful; the paper checks that agreement on only 10% of the data and shows only statistical uncertainties in the comparison plots.","fun_headline_variants_meta":{"raw":{"variants":["IceCube 8-year data improves dark-matter sensitivity 3.8x","Dark matter in Earth: IceCube's 8-year hunt gains 3.8x","BDT selection lifts IceCube dark-matter sensitivity 3.8-fold","8 years of IceCube data sharpen Earth-core dark-matter search","IceCube's 8-year data improves Earth dark-matter probe 3.8x"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000879,"raw_usage":{"total_tokens":3772,"prompt_tokens":889,"completion_tokens":2883,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":505,"completion_tokens_details":{"reasoning_tokens":2777}},"tokens_in":505,"tokens_out":2883,"duration_ms":19043,"temperature":1.0,"reasoning_tokens":2777,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:22:10.334049+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Unblind the full eight-year signal region and compare the observed post-cut event rate and zenith-angle distribution with the simulated 0.19 mHz background: any deviation beyond the quoted uncertainties would invalidate the projected median limit. Alternatively, compute the 90% upper limit directly from the observed data; if it is statistically consistent with the previous one-year IceCube limit rather than about 3.8 times more sensitive, the central claim is falsified.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the previous one-year IceCube analysis whose sensitivity is the factor-3.8 comparison baseline, as well as the likelihood method being extended."},{"cited_title":"Griest and D","cited_arxiv_id":null,"evidence_quote":"Gives the equilibrium formula relating the annihilation rate to the capture rate and annihilation cross-section, defining the target quantity of the search."},{"cited_title":"Sivertsson and J","cited_arxiv_id":null,"evidence_quote":"Provides the Earth capture rate as a function of WIMP mass that connects the flux sensitivity to the spin-independent scattering cross-section."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Simulates the atmospheric muon background from cosmic-ray air showers in the background estimate."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the treatment of neutrino oscillations inside Earth for energies below 100 GeV in the simulation."},{"cited_title":"Blennow, J","cited_arxiv_id":null,"evidence_quote":"WimpSim generates the neutrino energy spectra for each annihilation channel and mass used for signal simulation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the frequentist prescription used to convert the binned likelihood into a 90% upper limit on signal events."},{"cited_title":"Dark Univ","cited_arxiv_id":null,"evidence_quote":"ANTARES limits are the external comparison for the annihilation-rate sensitivity in the W+W- channel."}],"review_version":1}