{"id":"e1c22967-dd08-40af-97ec-03d1cb4e69b4","arxiv_id":"2506.22659","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":5,"one_line_summary":"Using 1440 hours of neutron multiplicity data at 1166 m.w.e., the authors set 90% upper limits on dark matter-nucleus cross sections under the assumption that all dark matter mass energy converts to hadronic energy in a lead target.","lead":"A reanalysis of old underground neutron detector data finds no dark matter signal and sets upper limits on dark matter-matter cross sections between about 1e-45 and 2e-42 cm2. The result probes dark matter masses below 10 GeV, where conventional direct detection experiments lose sensitivity.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equation 14 inverts the target-area and A-scaling factors when converting the 90% CL likelihood bound into a cross section; the quoted limits do not follow from Eq. 13 as written.","rationale":"The reader's stated weakest assumption is the complete-hadronic-conversion hypothesis (Section VIII, f_average = 1), but the reader's rationale also flags Eq. 14. I focus on Eq. 14 because it is the unique algebraic link between the likelihood output and all quoted cross-section limits, and because the error is checkable analytically with the paper's own parameters. The null measurement itself is plausible: the detector calibration against 252Cf, the power-law description of multiplicities, and the Geant4 checks against Miyake fluxes and NESSI spallation data are all reasonable, and I am not disputing the absence of an observed excess. The problem is the final step: as printed, Eq. 14 cannot be derived from Eq. 13 and has the wrong dimensions. If the figure was produced with a different, correct formula, the manuscript needs a corrected equation and a reproducibility statement; if not, the numerical limits are not established. This warrants rejection of the current version rather than acceptance, because the central quantitative claim is not derivable from the equations provided. I partially agree with the reader because the reader also mentions the equation but identifies a different weakest assumption.","tokens_in":16052,"tokens_out":19056,"duration_ms":210832,"concrete_test":"Take a representative mass (e.g., mDM = 10 GeV), derive beta90 from Eq. 8 using the stated 1166 m.w.e. six-event spectrum and simulation p, and compute sigma90 two ways: (i) Eq. 14 as printed; (ii) the correct inversion sigma90 = beta90 * N1166 * Starget / (NDM * NPb * A^2 * epsilon), with NDM = 6.0e16 GeV/mDM, Starget = 900 cm^2, NPb = 8.9e26, A = 207, and the simulated acceptance epsilon at that mass. If curve (ii) differs from Fig. 17 by more than the likelihood uncertainty, the printed Eq. 14 is the cause and the limits must be re-derived; if Fig. 17 matches (ii), Eq. 14 is a typographical error and the central limits may survive. The dimensional check alone settles it: the correct inversion has units cm^2, while Eq. 14 as printed has units cm^-2.","verdict_should_be":"REJECT","load_bearing_attack":"Section X, Eq. 14 is the sole step that turns the beta upper limit into a cross section. From Eq. 13, Nevents = NDM * NPb * sigma_conversion / Starget, with NDM = rho*v*t*Starget/mDM from Eq. 12, so the target area cancels in the physical event rate. Solving for the per-nucleon SI cross section with sigma_conversion = A^2 * sigma_chiN and signal-count upper limit Nsig,UL = beta90 * N1166 gives sigma_chiN,90 = Nsig,UL * Starget / (NDM * NPb * A^2 * epsilon). Equation 14 instead writes sigma90 = N1166 / [beta90 * (Starget/A^2) * NDM * NPb * epsilon], placing Starget in the denominator and A^2 in the numerator, and uses N1166/beta90 rather than beta90 * N1166. Dimensionally Eq. 14 yields cm^-2, not cm^2. For Pb (A ~ 207, Starget = 900 cm^2) the misplacement changes the scale by A^4/(beta90^2 * Starget^2) ~ 10^3-10^4, so the curves in Fig. 17 cannot be reproduced from the printed equations. Since every quoted limit passes through this formula, the central numerical claim is unsupported as written; the SD version has the same 1/Starget inversion.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes neutron multiplicity data from the NMDS-II detector, a 305 kg lead target instrumented with 3He counters, operated at two underground depths (583 m.w.e. for 6504 hours and 1166 m.w.e. for 1440 hours). The authors model cosmic-ray-muon-induced neutron backgrounds with Geant4, find that both data and simulation multiplicity spectra are well described by a power law, and search for an excess at high multiplicities (n >= 5) at the deeper site, where only six events were observed. Two dark-matter interaction models, a spallation model and a pion fireball model, are used to predict the neutron multiplicity distribution from dark-matter energy deposition in the lead target. A likelihood fit with a signal fraction beta is performed, and upper limits on the dark-matter-nucleon cross section are derived. The paper claims spin-independent limits near 1e-45 cm^2 and spin-dependent limits near 2e-42 cm^2 for dark-matter masses from 300 MeV to 100 GeV, assuming full hadronic conversion of the dark-matter rest mass.","tokens_in":16392,"tokens_out":14730,"duration_ms":156060,"significance":"If the quoted limits are correct, the paper would provide an independent indirect search for dark-matter interactions in the sub-GeV to 100 GeV mass range using neutron multiplicity, a channel that is complementary to conventional direct detection. The identical detector operated at two depths is a useful feature for validating the muon-induced background model, and the simulation effort, including muon propagation from sea level and full shower modeling, is substantial. The paper also validates its spallation simulation against external experimental data. However, the central cross-section conversion equation is incorrect as printed, and the statistical definition of the upper limit is non-standard; until these issues are resolved, the numerical limits cannot be accepted. The paper does not compare its results with existing dark-matter limits, which would be important for context.","major_comments":[{"comment":"Equation (14) does not follow from Equations (12) and (13). From Eq. (13), the expected number of conversion events is Nevents = NDM * NPb * sigma_conversion / Starget, and substituting NDM = rho*v*t*Starget/mDM from Eq. (12) cancels Starget, giving Nevents = rho*v*t*NPb*sigma_conversion/mDM. Setting the 90% signal count to beta90*N1166 and solving for the per-nucleon spin-independent cross section yields sigma_chiN,90 = beta90*N1166*Starget / (NDM*NPb*A^2*epsilon). Equation (14) as displayed instead appears to write N1166 / [beta90 * (Starget/A^2) * NDM * NPb * epsilon], which inverts the beta factor (using N1166/beta90 rather than beta90*N1166) and places Starget in the denominator and A^2 in the numerator. Under this reading the expression has units of cm^-2, not cm^2, so it cannot be the correct conversion. Because every quoted limit in the abstract and Fig. 17 passes through this formula, the central numerical claim is unsupported as written. The authors must correct Eq. (14) and verify that the limits in Fig. 17 were computed with the corrected expression.","section":"Section X, Eq. (14)"},{"comment":"The significance measure in Eq. (9) is non-standard. With Ssig = beta*N1166 and Bback = (1-beta)*N1166, the expression Ssig/sqrt(Ssig+Bback) equals beta*sqrt(N1166), not beta/sqrt(N1166) as the printed equality suggests; this is not a likelihood-ratio test statistic and does not account for the number of multiplicity bins, the look-elsewhere effect, or systematic uncertainties in the background shape. The 90% upper limit beta90 is defined by integrating the normalized likelihood over beta from 0 to 90%, which is a Bayesian credible interval with an implicit uniform prior on beta, not a frequentist confidence interval. The paper should state this statistical construction explicitly, including the prior, and ideally provide frequentist coverage checks, since the entire limit rests on only six events.","section":"Section VII, Eq. (9), and Section X"},{"comment":"The paper acknowledges in Section VI C that the Geant4 simulation and the data disagree in the amplitude k of the multiplicity distribution, and it sets this disagreement aside. However, the dark-matter likelihood in Eq. (7) uses the simulated background shape with the index p fixed to the Geant4 value (2.31 at 1166 m.w.e.), while the data are consistent with p = 2.13 +/- 0.37. The uncertainty in p and the normalization is not propagated into beta90 or the cross-section limits. Given the small event sample, the authors should either fit p and beta simultaneously in the likelihood or explicitly propagate the uncertainty in the background shape into the final limits, rather than treating the simulated p as exact.","section":"Section VI C and Section VII"},{"comment":"The limits assume complete conversion of the dark-matter rest mass into hadronic energy in the lead target (f_average = 1). This assumption is stated, but the paper does not quantify how the limits scale with f_average or discuss the model dependence of the two bracketing hadronic models for dark-matter masses near threshold. The abstract states the assumption, but the conclusions should clearly flag that any model in which dark matter deposits only a fraction of its rest mass as hadronic energy, or produces different neutron multiplicities, would weaken the quoted limits. A simple scaling relation sigma ∝ 1/f_average, along with a short discussion of the validity range of the fireball model for masses below about 1 GeV, would make the caveat precise.","section":"Section VIII and Section X"}],"minor_comments":[{"comment":"The header contains a typo: 'Miiyake' should be 'Miyake'.","section":"Table I"},{"comment":"The displayed equality in Eq. (9) appears to have an algebraic slip: beta*N1166/sqrt(N1166) equals beta*sqrt(N1166), not beta/sqrt(N1166). Please correct the typo.","section":"Section VII, Eq. (9)"},{"comment":"There are numerous typographical errors, for example 'indepentently' in Section VI B, 'EXPERIMA T AL' and 'SIMUALTION' in the Section VI header, 'compered' in Section VIII, and 'curvture' in Appendix A. A careful proofread is needed.","section":"Throughout"},{"comment":"The paper does not compare its limits with existing dark-matter constraints, such as direct detection limits at higher masses or other neutron-based searches. A brief comparison, even in a figure inset, would help the reader judge the significance of the result.","section":"Section X / Conclusion"},{"comment":"The pion fireball model uses a Hagedorn temperature kT = 0.165 GeV taken from proton-nucleus interactions at 70 GeV; for dark-matter masses near the pion-production threshold (300 MeV to about 1 GeV), the applicability of this temperature and the Planck distribution should be discussed.","section":"Section VIII B"},{"comment":"The simulation uses Geant4-11.01 while the data were collected in 2001-2003; this is acceptable, but the paper should explicitly state that a modern Geant4 version was used to model the historical detector and that this version choice was validated against the 252Cf calibration.","section":"Section V A"}],"recommendation":"major_revision","confidential_remarks":"The core dataset and simulation work are potentially valuable, and the null result is plausible. However, the printed conversion formula in Eq. (14) is seriously incorrect, and the statistical method needs to be clarified. If the authors can confirm that their numerical calculations used the corrected formula and can present the limits with a properly defined and validated statistical procedure, the paper may be publishable. I would also suggest asking for a comparison with existing limits in the same mass range."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a look, but the central numbers are not reproducible from the printed equations. The paper reanalyzes archived NMDS-II data (2001-03) to set limits on dark-matter-induced neutron bursts in a 305 kg Pb target. That is genuinely new: no one else has used this detector for a DM search, and the neutron multiplicity signature combined with two extreme hadronization models (spallation vs pion fireball) is a sensible way to bracket the unknown final states. The null result is reported honestly, and the background work is diligent: muon propagation is validated against Miyake and mine measurements, and the Geant4 multiplicity shapes match the data at both depths in the slope p, even though the amplitude k differs (a mismatch they openly acknowledge and then set aside).\n\nThe soft spot is load-bearing. Eq. 14, which converts the 90% CL beta limit into a cross section, appears to invert the target area and A-scaling relative to Eq. 13. If the equation is read as written, the dimensions come out as cm^-2 rather than cm^2, and the resulting limits cannot match Fig. 17. The correct algebra from Eq. 13 gives sigma_chiN = (beta90 * N1166) * S_target / (NDM * N_Pb * A^2 * epsilon), with S_target in the numerator and A^2 in the denominator. Eq. 14 as printed puts S_target in the denominator and A^2 in the numerator, and uses N1166/beta90 instead of beta90 * N1166. This is not a cosmetic typo; it changes the scale by orders of magnitude. Every quoted limit passes through this step.\n\nThere are additional issues. The significance formula in Eq. 9 (theta = beta_peak / sqrt(N1166)) is not a valid significance for a Poisson search; it ignores the background-only variance. No systematic uncertainties are propagated into the limits. And the spin-dependent limit is derived with a simple A-scaling that is inappropriate for natural lead, whose dominant isotopes have zero net spin. The data set is also small: six events at 1166 m.w.e. carry the entire constraint. So while the analysis idea is sound and the paper is transparent about many of its choices, the quantitative limits as presented are not supported.\n\nFor a reader interested in low-mass DM search techniques, the paper is worth reading for the detector and background modeling sections. But I would not cite the limits until the normalization error is fixed and the SD treatment is redone. It deserves a serious referee: the underlying data and method could produce a valid constraint after major revision.","headline":"A genuinely new reanalysis of archived neutron-multiplicity data, but the printed limit equation inverts the target-area and A-scaling factors and the quoted cross sections do not follow.","tokens_in":16882,"tokens_out":3527,"would_cite":false,"duration_ms":36388,"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":"Neutron-multiplicity data from a lead cube at 1166 m.w.e. show no dark-matter signal, yielding 90% upper limits near $\\sim 10^{-45}\\,\\mathrm{cm}^2$ for spin-independent interactions from 300 MeV to 100 GeV.","keywords":["dark matter","indirect detection","neutron multiplicity","lead target","cross-section limits","cosmic-ray muons","low-mass dark matter","Monte Carlo simulation"],"falsifier":"Record the neutron multiplicity spectrum at 1166 m.w.e. for about ten times the current 1440-hour exposure: the no-signal claim is falsified if high-multiplicity events ($n\\ge 20$) appear at a rate incompatible with the simulated cosmic-ray-muon background. The complete-conversion assumption can be tested separately by measuring neutron yields from protons and pions of known energy in a 30 cm lead cube; if the true deposited-energy fraction is substantially below one, the reported limits must be scaled upward.","tokens_in":15818,"feed_emoji":"⚛️","tokens_out":12262,"duration_ms":123755,"temperature":0.7,"pith_summary":"The paper reanalyzes 1440 hours of NMDS-II neutron-multiplicity data taken with a 305 kg lead cube at 1166 meters-water-equivalent underground, searching for dark-matter interactions that produce bursts of neutrons. It finds no excess beyond cosmic-ray-muon-induced background, and derives 90% confidence upper limits on dark-matter-matter cross sections: spin-independent limits at roughly $\\sim 10^{-45}\\,\\mathrm{cm}^2$ and spin-dependent limits at roughly $2\\times10^{-42}\\,\\mathrm{cm}^2$, for dark-matter masses between 300 MeV and 100 GeV. The result matters because conventional direct-detection searches lose sensitivity below about 10 GeV, while a lead target amplifies low-mass dark-matter signals through high neutron multiplicity and coherent $A^2$ spin-independent scattering. The limits rest on the assumption that each dark-matter interaction converts its full rest-mass energy into hadronic energy inside the lead target.","feed_headline":"Lead-cube neutron search sets dark matter limits at 10^-45 cm^2","feed_subtitle":"1440 hours of neutron data at 1166 m.w.e. show no signal, constraining dark matter that direct searches miss.","key_machinery":"The carrying object is the NMDS-II detector: a 305 kg, 30 cm lead cube surrounded by 60 $^{3}$He proportional counters embedded in polyethylene, which records the number of neutrons produced per interaction. Lead acts as a high-gain neutron amplifier, and the analysis models the cosmic-ray background as a power-law multiplicity distribution $k\\,n^{-p}$, justified through a sample-space-reducing cascade mechanism in which each neutron emission shrinks the set of available nuclear states. A Monte Carlo simulation propagates muons from sea level through rock to the cavern and generates neutron multiplicity spectra from both muon showers and the two dark-matter models; the search then fits the 1166 m.w.e. data with the mixture $P(\\beta;n)=(1-\\beta)k\\,n^{-p}+\\beta S_{\\mathrm{DM}}(n)$, and the 90% integration of the resulting likelihood yields the cross-section limits.","core_discovery":"The central claim is that the NMDS-II data contain no resolvable dark-matter signal: the most likely dark-matter fraction in the likelihood fit peaks at a significance around $1.5\\sigma$, corresponding to about four excess events for deposited energies between 2 and 20 GeV, which the paper does not claim as a detection. Assuming complete hadronic energy conversion, the 90% upper limits are at the level of roughly $10^{-45}\\,\\mathrm{cm}^2$ for spin-independent interactions and $2\\times10^{-42}\\,\\mathrm{cm}^2$ for spin-dependent interactions across the dark-matter mass range 300 MeV to 100 GeV. The two deliberately extreme final-state models, a single-proton spallation and a pion fireball with a limiting temperature, produce nearly the same limits, so the bound is not sensitive to the assumed hadronization pattern.","pith_inferences":["Not tested in the paper: lowering the multiplicity threshold below five neutrons and using the 583 m.w.e. spectrum as a direct background template would extend the search to lower deposited energies, and the reach would then be limited by how well the power-law index is known rather than by the six-event sample.","Because the 90% upper limit scales roughly as the inverse square root of the number of background events, a year or more of running at 1166 m.w.e. with the same detector would push the spin-independent bound below $10^{-45}\\,\\mathrm{cm}^2$ without any hardware change.","A natural continuation would apply the same lead-cube multiplicity method to hidden-sector dark-matter models whose decay products include hadrons; the fireball simulation would need to be redone with the specific hidden-sector decay kinematics before quoting limits.","The power-law background extrapolation could be checked by comparing the 583 m.w.e. data and simulation at multiplicities above 30, where a break in the power law would signal that the background shape used at 1166 m.w.e. needs revision."],"forward_implications":["Spin-independent dark-matter-nucleon cross sections above about $10^{-45}\\,\\mathrm{cm}^2$ are excluded at 90% confidence for masses from 300 MeV to 100 GeV, under the complete-hadronic-conversion assumption.","The near-agreement of the spallation and fireball limits means the result is stable across very different hadronic final states, so it does not hinge on the details of how dark matter converts to pions or protons.","The $A^2$ coherent enhancement for spin-independent scattering is what lets a lead target set competitive limits at low dark-matter mass where direct-detection experiments lose sensitivity.","The 583 m.w.e. data set, with about 36 times more muon-induced events, serves as a high-statistics control that validates the cosmic-ray background shape used in the 1166 m.w.e. analysis.","If dark matter instead converts its energy mainly into leptons or photons, or deposits only a fraction of its rest mass in the target, these limits do not apply."],"supporting_citations":[{"why":"Supplies the wall-effect correction that brings the simulated $^{3}$He counter efficiency from 28% to the measured 23.2%, fixing the absolute normalization of every neutron count.","marker":"[15]"},{"why":"Provides the sea-level cosmic-ray muon energy and angular spectrum and the muon energy-loss parametrization used to propagate muons to the experimental depths.","marker":"[16]"},{"why":"Gives the empirical underground muon depth-intensity and angular formulas used to cross-check the simulated muon flux and angular distributions at both depths.","marker":"[17]"},{"why":"Supplies the measured underground muon flux data that set the flux normalization and validate the simulation's muon rate.","marker":"[18]"},{"why":"Establishes that sample-space-reducing cascades generate power-law distributions, the theoretical basis for fitting the cosmic-ray background as $k n^{-p}$.","marker":"[22]"},{"why":"Supports the claim that cascade power-law exponents are independent of nuclear structure and equal 2 under energy conservation, justifying the fixed background exponent in the likelihood.","marker":"[23]"},{"why":"Provides measured neutron multiplicities from protons on Pb, W, and Hg targets, the data against which the spallation dark-matter model is validated.","marker":"[28]"},{"why":"Supplies the pion energy distribution used to generate pion multiplicities in the fireball dark-matter model.","marker":"[30]"}],"fun_headline_variants":["Lead-cube neutron search sets dark matter limits","No dark matter neutrons: limits at 1e-45 cm^2","Neutron detector constrains dark matter cross sections","Dark matter neutron search yields null result, sharp limits"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The limits hold only if every dark-matter interaction deposits all of its rest-mass energy as hadrons inside the lead target; if the real interaction converts a fraction of that energy to neutrons, or emits mostly leptons or photons, the quoted cross-section bounds do not apply, and the six events at 1166 m.w.e. carry the entire statistical weight.","fun_headline_variants_meta":{"raw":{"variants":["Lead-cube neutron search sets dark matter limits","No dark matter neutrons: limits at 1e-45 cm^2","Neutron detector constrains dark matter cross sections","Dark matter neutron search yields null result, sharp limits"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000591,"raw_usage":{"total_tokens":2810,"prompt_tokens":1021,"completion_tokens":1789,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":637,"completion_tokens_details":{"reasoning_tokens":1722}},"tokens_in":637,"tokens_out":1789,"duration_ms":14264,"temperature":1.0,"reasoning_tokens":1722,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T22:02:40.729470+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Record the neutron multiplicity spectrum at 1166 m.w.e. for about ten times the current 1440-hour exposure: the no-signal claim is falsified if high-multiplicity events ($n\\ge 20$) appear at a rate incompatible with the simulated cosmic-ray-muon background. The complete-conversion assumption can be tested separately by measuring neutron yields from protons and pions of known energy in a 30 cm lead cube; if the true deposited-energy fraction is substantially below one, the reported limits must be scaled upward.","supporting_citations":[{"cited_title":"Shalev, Z","cited_arxiv_id":null,"evidence_quote":"Supplies the wall-effect correction that brings the simulated $^{3}$He counter efficiency from 28% to the measured 23.2%, fixing the absolute normalization of every neutron count."},{"cited_title":"displayed in Figure 11","cited_arxiv_id":null,"evidence_quote":"Provides the sea-level cosmic-ray muon energy and angular spectrum and the muon energy-loss parametrization used to propagate muons to the experimental depths."},{"cited_title":"For spin-independent cross sections, the sensitivity per nucleon requires dividing the factor of the atomic number A, twice","cited_arxiv_id":null,"evidence_quote":"Gives the empirical underground muon depth-intensity and angular formulas used to cross-check the simulated muon flux and angular distributions at both depths."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the measured underground muon flux data that set the flux normalization and validate the simulation's muon rate."},{"cited_title":"Kiryunin, H","cited_arxiv_id":null,"evidence_quote":"Establishes that sample-space-reducing cascades generate power-law distributions, the theoretical basis for fitting the cosmic-ray background as $k n^{-p}$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supports the claim that cascade power-law exponents are independent of nuclear structure and equal 2 under energy conservation, justifying the fixed background exponent in the likelihood."},{"cited_title":"Bauche-Arnoult and J","cited_arxiv_id":null,"evidence_quote":"Provides measured neutron multiplicities from protons on Pb, W, and Hg targets, the data against which the spallation dark-matter model is validated."},{"cited_title":"Letourneau, J","cited_arxiv_id":null,"evidence_quote":"Supplies the pion energy distribution used to generate pion multiplicities in the fireball dark-matter model."}],"review_version":1}