REVIEW 3 major objections 4 minor 23 references
DarkSide-50's multi-scatter search excludes composite dark nuclei from 10^7 to 10^13 GeV/c^2
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · deepseek-v4-flash
2026-08-03 02:50 UTC pith:MD3DTOHZ
load-bearing objection First dual-phase argon UHDM search with a solid pipeline, but Eq. (7) uses the wrong mass in the flux normalization — if the code follows it, the exclusions are off by orders of magnitude. the 3 major comments →
First Nuclear Ultra-Heavy Dark Matter Search in Argon Time Projection Chambers with the DarkSide-50 Experiment
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
For the first time, a dual-phase liquid argon TPC has been used to search for nuclear ultra-heavy dark matter, composite objects made of many dark nucleons that form when the dark sector has no long-range force. Using the 532-day low-radioactivity argon run of DarkSide-50, the analysis applies quality cuts that select multi-scatter event topologies, accounts for energy loss in the Earth's overburden with a continuous slowing-down model, and derives exclusion limits on the UHDM-nucleon cross section for four dark-nucleon masses. The central result is a 90% C.L. excluded band in total-mass versus cross-section space, bounded below by the saturated-overburden cross section and above by the geom
What carries the argument
The signal is a series of S1 scintillation pulses from successive elastic scatters of a UHDM nucleus as it crosses the TPC. The recoil spectrum is shaped by the dark-nucleus form factor Fχ = 3j₁(qRχ)/(qRχ), a spherical top-hat density, and by the geometric radius Rχ = (9πMχ/4mχ⁴)^{1/3} obtained by assuming dark-nucleon forces equivalent to Standard Model ones. Overburden energy loss is computed with a continuous-slowing-down approximation, valid because the UHDM mass vastly exceeds the target mass, and the search is bounded from below by the saturated-overburden cross section σs. The analysis uses only S1, since S2 drift times are comparable to the UHDM transit time.
Load-bearing premise
The search assumes a dark nucleus has the same force structure as ordinary strong interactions, giving radius Rχ = (9πMχ/4mχ⁴)^{1/3} and a top-hat density; if the dark binding force differs, the entire exclusion band shifts in cross section and mass.
What would settle it
A single observed event in the multi-scatter channel with an S1 train of two or more pulses separated by microseconds and total energy in the 100–8000 PE window would directly contradict the null result. Alternatively, a dedicated calculation showing that the dark-nucleus radius is not set by Standard-Model-like forces would undermine the cross-section axis of the limits.
If this is right
- A liquid argon TPC can detect composite UHDM through the multi-scatter topology, where a single crossing produces many S1 pulses.
- The excluded band spans UHDM masses from roughly 10⁷ to 10¹³ GeV/c² for all four dark-nucleon masses considered.
- The limits apply only for the assumed dark-nucleus radius; changing mχ shifts the band, so future searches can map the internal structure of UHDM.
- The S1-only multi-scatter analysis with overburden energy loss is directly applicable to larger argon TPCs.
- Cross sections above the geometric limit are already ruled out by self-interaction constraints, so this search closes the remaining window.
Where Pith is reading between the lines
- If the dark binding force differs from the Standard Model strong force, the radius formula and hence the excluded band would shift by orders of magnitude; the paper's limits quantify a specific dark-sector model, not a universal bound.
- The same multi-scatter analysis could be applied to other noble-liquid detectors, since the signature depends only on a track of energy deposits, not on argon's specific response.
- A testable extension is to re-analyze the same dataset with alternative dark-nucleus form factors (e.g., Gaussian) to gauge model dependence.
- The dark-nucleon mass mχ acts as an internal-structure parameter; future experiments could treat it as a continuous variable rather than sampling a few fixed values.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports the first search for nuclear ultra-heavy dark matter (UHDM) in a dual-phase liquid argon TPC, using the 532-day low-radioactivity argon campaign of DarkSide-50. UHDM candidates are modeled as composite dark nuclei of total mass Mχ and constituent dark-nucleon mass mχ; the analysis includes overburden energy loss with the Verne package, a multi-scatter signal simulation with LAr response from ARIS measurements, and a profile-likelihood fit to radiogenic backgrounds. The central result is a set of 90% CL exclusion curves in the (Mχ, σχ,n) plane for mχ = 10, 50, 100, and 500 GeV/c² (Fig. 9).
Significance. If the numerical results are correct, this is the first direct UHDM constraint from a dual-phase LAr TPC and a useful complement to existing searches by DEAP-3600 and LZ. Strengths of the manuscript are the coherent use of an external theoretical signal model, the published DarkSide-50 background model with small nuisance pulls (all |pull| < 0.03), and the explicit treatment of multi-scatter event topologies. However, the central quantitative claim depends on the event-rate normalization in Eq. (7), which appears to be internally inconsistent with the composite-mass formalism used everywhere else in the paper. Because this issue directly affects the reported exclusion curves and p-values, the manuscript cannot be accepted in its present form.
major comments (3)
- [§IV A, Eq. (7)] The event-rate formula N = T(ρχ/mχ)∫|v|f d³v dA uses the constituent dark-nucleon mass mχ rather than the composite UHDM mass Mχ. The incident flux of composite dark-matter nuclei is nχ⟨v⟩ = ρχ⟨v⟩/Mχ, so the printed expression overestimates the event rate by Mχ/mχ, which is ≈10⁵–10¹² over the parameter range considered. Since the right-hand, exposure-limited edge of Fig. 9 and the p-value calculation for μ=1 depend directly on this normalization, the quoted exclusions are not justified from the formula as written. The authors must either state explicitly that the computation uses Mχ (correcting Eq. (7) as a typographical error) or recompute the limits with the correct density.
- [§II, Eqs. (1)–(3); §IV A, Eq. (7)] The manuscript is internally inconsistent about which mass governs UHDM kinematics. Equation (1) and the surrounding text treat σχ,n as a whole-nucleus–nucleon cross section and use reduced masses μχ,N and μχ,n that should be evaluated with the composite mass Mχ, while Eq. (7) uses mχ for the flux normalization. If the reduced masses in Eq. (1) are instead evaluated with mχ, the A² coherence enhancement and the form-factor suppression shown in Fig. 5 would also be affected. Please define explicitly which mass enters each reduced mass, the number density, and the kinetic-energy-loss formalism, and confirm that the signal Monte Carlo uses the same convention.
- [§IV A, detector response] The treatment of S1 for multi-scatter events is not fully specified. The text states that S1 is evaluated using the sum of all energy deposits in a TPC pass, but the LAr response parameterized by L_eff(E_R) is nonlinear; summing individual recoil energies before applying the quenching model differs from summing the S1 produced by each individual scatter. Because the ROI boundaries (100–8000 PE) and the acceptance map in Fig. 6 are derived from S1, the magnitude of this approximation should be quantified, or the simulation should be described more precisely.
minor comments (4)
- [§I] Typographical error: “time projection champer” should be “time projection chamber.”
- [§III] The sentence “masses below ∼Mχ = 10¹⁰ GeV/cm²” mixes a mass with a dimension of inverse length; this should be GeV/c², and the phrasing should be clarified.
- [§IV A] The sentence “No events had more than one S1 pulse” appears to contradict the preceding statement that successive energy deposits are not expected to merge into a single S1 pulse. Please clarify whether the toy simulation tested full multi-scatter events or individual energy deposits.
- [References] Reference [11] contains a malformed author entry (“Tailby.”); this should be corrected.
Circularity Check
No significant circularity: the UHDM exclusion is a data-driven analysis using an external signal model and measured detector response.
full rationale
The claimed derivation chain is not circular. The UHDM signal model is taken from external theory papers (Hardy et al. [14], Coskuner et al. [15], Butcher et al. [16]), with the form factor (Eq. 2) and radius relation (Eq. 3) explicitly adopted from those external works rather than derived from DarkSide-50 data or from the authors' prior results. Overburden transport uses the external Verne toolkit [18,19]; the liquid-argon S1 response uses measured ARIS/DarkSide scintillation calibrations [21,22]; the background model uses measured radioactive activities [8] and is validated by a background-only fit with nuisance pulls consistent with zero; and the limit-setting uses the standard Cowan et al. likelihood [23]. No parameter is fitted to UHDM candidates and then renamed as a prediction, and no equation reduces to another by construction. The self-citations in refs. [6–9,21,22] provide detector calibration, background, and dataset descriptions; they are external measurements and do not smuggle in the UHDM signal claim. The apparent use of mχ rather than Mχ in Eq. (7) is a potential unit/consistency issue that could affect the exposure-limited boundary, but it is not a circular reduction of the prediction to its inputs and cannot be scored as circularity without confirmation of the implemented analysis code.
Axiom & Free-Parameter Ledger
free parameters (4)
- dark nucleon mass mχ =
10, 50, 100, 500 GeV/c² (scan points, not fitted)
- dark-matter mass density ρχ =
0.3 GeV/cm³ (standard halo value, unstated in text)
- halo velocity parameters =
v0 = 220 km/s, Maxwellian distribution (implicit in Fig. 4)
- nuclear recoil scintillation response (Leff, Fq) =
ARIS fits (refs. [21,22])
axioms (7)
- domain assumption Composite dark-matter nuclei with up to 10²⁰ constituents form via Big-Bang Dark Nucleosynthesis in a dark sector without a long-range force (refs. [14,15]).
- domain assumption Dark nucleus radius Rχ = (9πMχ/4mχ⁴)^{1/3} assuming SM-equivalent dark-nucleon forces (Eq. 3).
- domain assumption Uniform-density spherical top-hat form factor Fχ = 3j₁(qRχ)/(qRχ) (Eq. 2).
- domain assumption Standard halo: ρχ = 0.3 GeV/cm³, v0 = 220 km/s, Maxwellian speed distribution with no galactic-parameter variation.
- domain assumption Continuous energy-loss formalism (Verne) accurately describes UHDM passage through Earth's overburden for Mχ ≫ mA.
- domain assumption LAr scintillation response models (Eqs. 5–6) with ARIS quenching describe the merged S1 of a multi-deposit track.
- ad hoc to paper Trigger and data-quality selection have ~100% efficiency for signal events.
read the original abstract
We report the first search for nuclear ultra-heavy dark matter (UHDM) in a dual-phase liquid argon time projection chamber using the DarkSide-50 experiment. Unlike conventional weakly interacting massive particles (WIMPs), nuclear UHDM candidates may be composed of many dark nucleons and scatter numerous times while passing through the detector. Accounting for energy loss through the Earth's overburden, we apply selection criteria optimized for multi-scatter event topologies using the 532-day low-radiation campaign of the DarkSide-50 detector. Excluded limits on the UHDM-nucleon scattering cross section for dark nucleon masses of $m_\chi = 10, 50, 100, 500 \, \mathrm{GeV/c^2}$ are presented.
Figures
Reference graph
Works this paper leans on
-
[1]
S. M. Faber and J. S. Gallagher, Masses and Mass-To- Light Ratios of Galaxies, Annual Review of Astronomy and Astrophysics17, 135 (1979)
1979
-
[2]
Okamoto, S
T. Okamoto, S. Yachi, and A. Habe, Galaxy Distribu- tion in Clusters of Galaxies, Symposium - International Astronomical Union183, 262 (1999)
1999
-
[3]
Planck Collaboration, Planck 2018 results. VI. Cosmo- logical parameters, Astronomy & Astrophysics641, A6 (2020), arXiv:1807.06209 [astro-ph.CO]
Pith/arXiv arXiv 2018
-
[4]
D. N. Spergel, L. Verde, H. V. Peiris,et al., First-Year Wilkinson Microwave Anisotropy Probe ( WMAP ) Ob- servations: Determination of Cosmological Parameters, The Astrophysical Journal Supplement Series148, 175 (2003)
2003
-
[5]
M. W. Goodman and E. Witten, Detectability of cer- tain dark-matter candidates, Physical Review D31, 3059 (1985)
1985
-
[6]
Agnes, T
P. Agnes, T. Alexander, A. Alton,et al., First results from the darkside-50 dark matter experiment at labora- tori nazionali del gran sasso, Physics Letters B743, 456 (2015)
2015
-
[7]
Agnes, L
P. Agnes, L. Agostino, A. I.F.M,et al.(DarkSide Collab- oration), Results from the first use of low radioactivity argon in a dark matter search, Phys. Rev. D93, 081101 (2016)
2016
-
[8]
P. Agnes, I. F. M. Albuquerque, T. Alexander,et al., DarkSide-50 532-day dark matter search with low- radioactivity argon, Physical Review D98, 102006 (2018), 1802.07198
Pith/arXiv arXiv 2018
-
[9]
P. Agneset al.(DarkSide), The veto system of the DarkSide-50 experiment, JINST11(03), P03016, arXiv:1512.07896 [physics.ins-det]
-
[10]
P. Adhikari, R. Ajaj, M. Alp ´ ızar-Venegas,et al.(DEAP Collaboration), First Direct Detection Constraints on Planck-Scale Mass Dark Matter with Multiple-Scatter Signatures Using the DEAP-3600 Detector, Physical Re- view Letters128, 011801 (2022), 2108.09405
Pith/arXiv arXiv 2022
-
[11]
E. Aprile, J. Aalbers, F. Agostini,et al.(The XENON collaboration), Projected WIMP sensitivity of the XENONnT dark matter experiment, Journal of Cosmology and Astroparticle Physics2020(11), 031, 2007.08796
Pith/arXiv arXiv 2007
-
[12]
Ebadi, A
R. Ebadi, A. Mathur, E. H. Tanin, and Tailby., Ultra- heavy dark matter search with electron microscopy of geological quartz, Phys. Rev. D104, 015041 (2021)
2021
-
[13]
Aalbers, D
J. Aalbers, D. S. Akerib, A. K. Al Musalhi,et al.(LZ Col- laboration), New constraints on ultraheavy dark matter from the lz experiment, Phys. Rev. D109, 112010 (2024)
2024
-
[14]
Hardy, R
E. Hardy, R. Lasenby, J. March-Russell, and S. M. West, Big bang synthesis of nuclear dark matter, Journal of High Energy Physics2015, 11 (2015)
2015
-
[15]
Coskuner, D
A. Coskuner, D. M. Grabowska, S. Knapen, and K. M. Zurek, Direct detection of bound states of asymmetric dark matter, Physical Review D100, 035025 (2019)
2019
-
[16]
Butcher, R
A. Butcher, R. Kirk, J. Monroe, and S. M. West, Can tonne-scale direct detection experiments discover nuclear dark matter?, Journal of Cosmology and Astroparticle Physics2017(10), 035
-
[17]
R. H. Helm, Inelastic and Elastic Scattering of 187-Mev Electrons from Selected Even-Even Nuclei, Physical Re- view104, 1466 (1956)
1956
-
[18]
B. J. Kavanagh, Earth scattering of superheavy dark matter: Updated constraints from detectors old and new, Physical Review D97, 123013 (2018)
2018
-
[19]
B. J. Kavanagh, Verne, Astrophysics Source Code Li- brary (2016), record ascl:1802.005
2016
-
[20]
J. Bramante, B. Broerman, R. F. Lang, and N. Raj, Satu- rated overburden scattering and the multiscatter frontier: Discovering dark matter at the Planck mass and beyond, Physical Review D98, 083516 (2018), 1803.08044
Pith/arXiv arXiv 2018
-
[21]
Agnes, J
P. Agnes, J. Dawson, S. De Cecco,et al.(The ARIS Collaboration), Measurement of the liquid argon energy response to nuclear and electronic recoils, Phys. Rev. D 97, 112005 (2018)
2018
-
[22]
Agnes, I
P. Agnes, I. Albuquerque, T. Alexander,et al., Sim- ulation of argon response and light detection in the DarkSide-50 dual phase TPC, Journal of Instrumenta- tion12(10), P10015
-
[23]
Cowan, K
G. Cowan, K. Cranmer, E. Gross, and O. Vitells, Asymp- totic formulae for likelihood-based tests of new physics, The European Physical Journal C71, 1554 (2011)
2011
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.