Pith. sign in

REVIEW 3 major objections 4 minor 43 references

DAMA/LIBRA's annual modulation signal is incompatible with the ANAIS-112 and COSINE-100 sodium-iodide experiments at 5.1σ when interpreted as dark-matter nuclear recoils, and reconciliation requires contrived fine-tuning.

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-04 11:20 UTC pith:5M6Y2HZO

load-bearing objection The full-spectrum comparison settles the qualitative question—DAMA can't be an ordinary recoil signal—but the '5.1σ' claim is a conditional point estimate that needs statistical validation and a less absolute abstract. the 3 major comments →

arxiv 2510.05216 v2 pith:5M6Y2HZO submitted 2025-10-06 hep-ph

DAMA/LIBRA and dark matter: decisive tension or contrived cancellation

classification hep-ph
keywords dark matterannual modulationDAMA/LIBRAANAIS-112COSINE-100sodium iodide detectorseffective field theorystatistical tension
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

Dark matter should imprint a yearly rise and fall in the rate of nuclear recoils in a detector as the Earth orbits the Sun. DAMA/LIBRA has reported such an annual modulation for more than two decades, but two newer experiments using the same sodium-iodide crystals, ANAIS-112 and COSINE-100, see no such signal. This paper asks whether the disagreement can be blamed on the detectors' different energy resolutions, quenching factors, and binning, rather than on the absence of dark matter. The authors find that, under any physically reasonable elastic-scattering dark-matter model, the tension between the datasets exceeds 5σ, and that models that reduce the tension require either many free parameters or unphysical opposite-sign modulation between sodium and iodine recoils. If true, this closes the loophole that DAMA's signal could be dark matter scattering that only shows up in certain NaI crystals.

Core claim

The paper's central claim is that the annual modulation signal reported by DAMA/LIBRA cannot be explained as dark-matter nuclear recoils once the full spectral information from ANAIS-112 and COSINE-100 is included. Using a dark-matter effective field theory with dimension-five and dimension-six operators and a standard Maxwellian halo, the combined fit to all three datasets is worse than the sum of the individual fits by δχ² = 51.5, corresponding to a 5.10σ tension. The authors also consider model-independent parameterisations of the recoil spectrum; in every case except those with many free parameters or with opposite-sign modulation on sodium and iodine, the tension remains above 5σ. Allow

What carries the argument

The argument is carried by a forward model that maps a single nuclear recoil spectrum for sodium and iodine into each experiment's measured modulation amplitudes, using the target mass fractions, measured quenching factors, detector energy resolution, and binning. On top of this sit two classes of generative models: a dark-matter effective field theory with dimension-five and -six operators (11 free parameters), and completely general parameterisations—binned recoil amplitudes or polynomial-exponential functions—with up to 40 free parameters. Tension between the datasets is quantified by the statistic δχ², the difference between the combined fit's χ² and the sum of the individual experiments

Load-bearing premise

The quoted 5.1σ tension assumes that the fit-quality difference δχ² follows a chi-squared distribution with as many degrees of freedom as the model has free parameters; if the null distribution is wider, the significance could be lower.

What would settle it

Simulate many noise-only realizations of the three experiments with no dark-matter modulation but with the same exposures and response functions, fit the same models to each realization, and record the distribution of δχ². If a value of 51.5 is not an extreme outlier in that simulated distribution, the claimed 5.1σ tension rests on an invalid statistical approximation.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • The same-technology loophole—that only pure NaI crystals can reveal DAMA's signal—is closed, since ANAIS-112 and COSINE-100 use the same target and still rule out the recoil interpretation.
  • No physically motivated elastic-scattering dark-matter model can explain the DAMA modulation, redirecting attention to inelastic scattering, electron interactions, or non-dark-matter explanations such as seasonal backgrounds.
  • The only generic fits that reduce the tension below 5σ require sodium and iodine recoil modulations to peak in opposite seasons, which forces the total unmodulated rate to be much larger than the observed modulation amplitude.
  • A future Southern-hemisphere NaI experiment could provide a decisive test by measuring both the modulation and the total rate in the same crystals.
  • Reanalyses that use only a few energy bins, as in previous combined analyses, underestimate the tension; using the full spectral information is necessary to see the incompatibility.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • A natural next step, which the paper does not take, is to apply the same forward model and δχ² test to future NaI datasets as they are released, giving an early model-independent compatibility check for any claimed modulation.
  • The opposite-sign sodium/iodine solutions imply a large unmodulated rate, which can already be checked against total-rate measurements from COSINE-100; the 'contrived cancellation' scenarios may therefore be falsifiable with data that already exist.
  • A Monte Carlo calibration of the tension statistic under the null hypothesis would pin down the exact significance; the qualitative conclusion, however, does not hinge on the precise p-value, because the tension persists across all physically motivated models and most generic parameterisations.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. This paper assesses whether the DAMA/LIBRA annual modulation signal can be reconciled with the null results from ANAIS-112 and COSINE-100 under the assumption of a common nuclear recoil signal from dark matter scattering. The authors perform fits to the modulation spectra using two approaches: a particle-physics EFT with dimension-5 and dimension-6 operators, and a set of model-independent parameterisations of the nuclear recoil spectrum (polynomial-exponential and binned amplitudes). They use the Maltoni–Schwetz tension statistic, δχ² of Eq. (5), to quantify the incompatibility between DAMA and the combined ANAIS+COSINE datasets. Their central finding is that, for the EFT analysis, the tension is 5.10σ with measured energy-dependent quenching factors, reducing to 4.63σ if DAMA's assumed constant quenching factors are used. For generic parameterisations, the tension remains above 5σ unless the number of fit parameters is large (more than 12 independent Na/I bins) or the sodium and iodine recoil amplitudes have opposite signs, which the authors argue is physically pathological. The paper concludes that the dark matter nuclear recoil interpretation of DAMA is ruled out at high significance by the combination of ANAIS-112 and COSINE-100.

Significance. If the central claim is correct, this paper closes a long-standing loophole in the DAMA interpretation: that differences in quenching, efficiency, resolution, and binning could allow a common nuclear recoil signal to be visible in DAMA but not in ANAIS-112 and COSINE-100. The analysis is careful and includes important robustness checks, such as repeating the EFT fit with DAMA's preferred constant quenching factors, which reduces the tension to 4.63σ but still disfavours the interpretation. A major strength is the use of the public GAMBIT/DDCalc ecosystem and the availability of the fitting code, which makes the analysis reproducible. The paper also goes beyond a single model by considering generic parameterisations, and it identifies the specific fine-tuned properties (opposite-sign Na/I modulation) required to reduce the tension. These are concrete and falsifiable predictions. However, the headline significance is a point estimate that does not include systematic uncertainties from quenching or halo parameters, and it relies on an asymptotic approximation for the tension statistic that is not tested for the high-dimensional generic fits.

major comments (3)
  1. [Quantifying the tension between the datasets, Eq. (5)] The paper assumes that δχ² follows a χ² distribution with n degrees of freedom, following Maltoni–Schwetz. This is only asymptotically valid for nested models with best-fit points in the interior of parameter space and enough data per bin. For the generic parameterisations in Table III, especially those with up to 20 free parameters fitted to a limited number of modulation bins, the effective number of degrees of freedom may be substantially smaller than n. Boundary effects (e.g., couplings or amplitudes driven to zero) and non-nested effective models can also distort the null distribution. The paper does not validate the statistic with Monte Carlo calibration or profile-likelihood methods. Since the quoted p-values (e.g., 3.37×10⁻⁷ for the EFT case) and the statement that 'more than 12 bins are required' rest entirely on this assumption, the tension values should be treated as condition
  2. [Table II and 'Detailed particle astrophysics test'] The fiducial 5.10σ tension assumes that DAMA's crystals have the same energy-dependent quenching factors measured for COSINE or ANAIS crystals (Refs. [22,32]). DAMA's quenching factors are unmeasured; the paper itself notes that using DAMA's assumed constant values (Q_Na=0.3, Q_I=0.09) reduces the tension to 4.63σ. No systematic uncertainty from quenching is propagated, nor are halo parameters or detector resolutions varied. Because the headline claim of 'decisive tension' sits near the 5σ threshold, this is a load-bearing point. The conclusion should either be softened to a robust lower bound (e.g., '>4.6σ under extreme quenching assumptions') or the analysis should profile over plausible quenching curves and report a range of tensions.
  3. [Table III and 'How generic is the tension?'] The conclusion that the tension only drops below 5σ with 'more than 12 independent Na/I bins' depends on the distributional assumption in Eq. (5). In the 20-bin case, the best-fit spectra require opposite-sign sodium and iodine modulation. The authors correctly note that this is pathological, but they do not impose the self-consistency constraint A0^T ≥ |A1^T| in the fits; they only state that this 'would place these scenarios under further strain'. This is acceptable as an argument, but it would strengthen the paper to show explicitly that imposing this constraint leaves the >5σ tension intact for all parameterisations considered.
minor comments (4)
  1. [Section 'How generic is the tension?'] Typo: 'spect ra' should be 'spectra'.
  2. [Table III] The notation is inconsistent: the text and table use 'polyxexp4', 'polyexp10', and 'polyxexp' interchangeably. Please unify to a single name, e.g., 'polyexp'.
  3. [Fig. 2 caption] '10 times 2 bins' is ambiguous; it would be clearer to write '10 bins per target' or '20 bins total'.
  4. [Eq. (3)] The integral over E_R extends to infinity, but the quenching factor Q_T(E_R) is only defined for finite E_R. For the resolution convolution this is harmless, but a note or a cutoff would avoid confusion.

Circularity Check

0 steps flagged

No significant circularity: the 5.1σ tension is a likelihood-ratio statistic computed from fits to independent experimental datasets, with no fitted quantity redefined as the target of the claim.

full rationale

The paper's derivation chain is self-contained and its central claim does not reduce to its inputs. The tension statistic in Eq. (5) is the Maltoni–Schwetz δχ² between separate best fits to DAMA/LIBRA and to ANAIS-112+COSINE-100, using free EFT or generic spectral parameters. No parameter is defined in terms of the tension outcome; the 'opposite-sign sodium/iodine' conclusion emerges from the best-fit recoil spectra rather than being imposed by construction. The quenching factors are taken from external measurements (Refs. [22] and [32]), and the paper explicitly recomputes the tension with DAMA's assumed constant quenching factors, reporting a reduced but still large 4.63σ tension. The GAMBIT/DDCalc/DirectDM software citations involve author overlap, but these are public, independently developed analysis tools whose likelihoods are tied to published experimental data; no load-bearing step invokes an unverified uniqueness theorem or an ansatz smuggled in via self-citation. The δχ²≈χ²_n approximation and the fixing of halo and resolution parameters are statistical and systematic assumptions that could affect the robustness of the quoted p-values, but they do not make the prediction equal to the fit inputs by construction. Therefore no circular step is identified.

Axiom & Free-Parameter Ledger

5 free parameters · 5 axioms · 0 invented entities

The central claim is conditional on the chosen model classes (EFT and generic recoil parameterisations), on the external quenching/resolution measurements used to map nuclear recoil energy to keVee, and on the asymptotic distribution of the tension statistic. The many fitted amplitudes in the generic parameterisations are themselves the 'fine-tuning' the paper quantifies; no new particles or forces are introduced.

free parameters (5)
  • EFT DM mass mχ = 17.1 GeV (all-data); 184 GeV (DAMA-only); 102 GeV (ANAIS+COSINE); Table I
    Fitted in GAMBIT/Diver scans over the modulation spectra.
  • EFT Wilson coefficients C_1^(5), C_2^(5), C_1^(6) ... C_4^(6) = Values in Table I, e.g., C_3^(6)=0.436 in DAMA-only fit
    Dimension-5 and dimension-6 operator coefficients; fitted to data.
  • Quark-coupling angles θ_1^(6)...θ_4^(6) = Values in Table I, e.g., θ_1=1.77π (all-data)
    Parameterise relative couplings to up- vs down-type quarks; fitted.
  • polyexp coefficients c_T^i and d_T (T=Na,I) = Not tabulated; models with 4, 6, 8, 10 free parameters
    Generic recoil-spectrum parameterisation in Eq. (15); fitted to each dataset combination.
  • Binned recoil amplitudes in 2N uniform ER bins = Not tabulated; scan from 2 to 20 independent Na/I amplitudes
    Model-independent binned spectra; the number of bins is the 'tuning' knob the paper scans.
axioms (5)
  • domain assumption Maxwellian halo velocity distribution with v0=240 km/s and vesc=528 km/s fixed to central values
    Used to compute the DM scattering rate in Eq. (1); halo parameters are not marginalised, only fixed.
  • domain assumption Dark matter is a Dirac fermion singlet with elastic scattering via EFT operators of dimension d≤6
    Section 'Detailed particle astrophysics test', Eqs. (6)–(14); excludes inelastic, exothermic, or non-nuclear-recoil models.
  • domain assumption Energy-dependent quenching factors measured on other NaI crystals apply to DAMA/LIBRA
    The main fits use quenching factors from Refs. [22,32] 'for all experiments'; DAMA's own crystals were never measured. A constant-quenching robustness check is performed.
  • domain assumption Detector resolution functions are as reported in DAMA, COSINE-100, and ANAIS-112 apparatus papers
    Resolution enters Eq. (3) and is taken from Refs. [29,30,31]; no uncertainty is propagated.
  • domain assumption The δχ² tension statistic follows a χ² distribution with n degrees of freedom
    Eq. (5) and surrounding text; the paper relies on the Maltoni–Schwetz asymptotic result without Monte-Carlo calibration for high-dimensional binned fits.

pith-pipeline@v1.3.0-alltime-deepseek · 11087 in / 15028 out tokens · 111728 ms · 2026-08-04T11:20:02.708708+00:00 · methodology

0 comments
read the original abstract

The ANAIS-112 and COSINE-100 experiments were constructed to test the long-observed dark matter-like annual modulation signal reported by DAMA. While they have reported null results in their annual modulation search, it remains possible that the combined effects of quenching, efficiency, resolution and binning could transform a common nuclear recoil rate into a signal that is visible in some detectors but not others. We assess the tension between DAMA/LIBRA and these latest experiments, under a range of hypotheses ranging from physical to general parameterisations of a common nuclear recoil input. We find that, in the most physically-motivated cases, the tension between DAMA and these other NaI experiments exceeds 5$\sigma$. Lowering the tension to reasonable values requires significant tuning, such as overfitting with large numbers of free parameters, and opposite-sign modulation between recoil signals on sodium versus iodine.

Figures

Figures reproduced from arXiv: 2510.05216 by Aaron C. Vincent, Anders Kvellestad, Felix Kahlhoefer, Giorgio Busoni, Jonathan M. Cornell, Lauren Street, Martin White, Masen Pitts, Will Handley.

Figure 1
Figure 1. Figure 1: FIG. 1. Comparison of experimental data to the predicted rates at the best-fit points of selected fits. [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Best-fit nuclear recoil spectra necessary to pro [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

43 extracted references · 29 linked inside Pith

  1. [1]

    Bernabeiet

    R. Bernabeiet. al.,Dark matter search,Riv. Nuovo Cim.26N1(2003) 1–73, [astro-ph/0307403]

  2. [2]

    Bernabeiet

    R. Bernabeiet. al.,Dark Matter: DAMA/LIBRA and its perspectives,SciPost Phys. Proc.12(2023) 025, [arXiv:2209.00882]

  3. [3]

    Kelso, P

    C. Kelso, P. Sandick, and C. Savage,Lowering the Threshold in the DAMA Dark Matter Search,JCAP09 (2013) 022, [arXiv:1306.1858]

  4. [4]

    S. Baum, K. Freese, and C. Kelso,Dark Matter implications of DAMA/LIBRA-phase2 results,Phys. Lett. B789(2019) 262–269, [arXiv:1804.01231]

  5. [5]

    Agnese, T

    SuperCDMS Collaboration: R. Agnese, T. Aramaki, et. al.,Results from the super cryogenic dark matter search experiment at soudan,Phys. Rev. Lett.120 (2018) 061802

  6. [6]

    Agnese, T

    R. Agnese, T. Aralis,et. al.,First dark matter constraints from a supercdms single-charge sensitive detector,Phys. Rev. Lett.121(2018) 051301

  7. [7]

    Agnes, I

    DarkSide Collaboration: P. Agnes, I. F. M. Albuquerque,et. al.,Low-mass dark matter search with the darkside-50 experiment,Phys. Rev. Lett.121(2018) 081307

  8. [8]

    Amole, M

    PICO Collaboration: C. Amole, M. Ardid,et. al.,Dark matter search results from the complete exposure of the pico-60c 3f8 bubble chamber,Phys. Rev. D100(2019) 022001

  9. [9]

    Aalbers, D

    LUX-ZEPLIN Collaboration: J. Aalbers, D. S. Akerib, et. al.,First dark matter search results from the lux-zeplin (lz) experiment,Phys. Rev. Lett.131(2023) 041002

  10. [10]

    Aprile, K

    XENON Collaboration: E. Aprile, K. Abe,et. al.,First dark matter search with nuclear recoils from the xenonnt experiment,Phys. Rev. Lett.131(2023) 041003

  11. [11]

    B. Ahmedet. al.,The NAIAD experiment for WIMP searches at Boulby mine and recent results,Astropart. Phys.19(2003) 691–702, [hep-ex/0301039]

  12. [12]

    S. C. Kimet. al.,New Limits on Interactions between Weakly Interacting Massive Particles and Nucleons Obtained with CsI(Tl) Crystal Detectors,Phys. Rev. Lett.108(2012) 181301, [arXiv:1204.2646]

  13. [13]

    Barbosa de Souzaet

    DM-Ice: E. Barbosa de Souzaet. al.,First search for a dark matter annual modulation signal with NaI(Tl) in the Southern Hemisphere by DM-Ice17,Phys. Rev. D 95(2017) 032006, [arXiv:1602.05939]

  14. [14]

    Amar´ eet

    J. Amar´ eet. al.,First Results on Dark Matter Annual Modulation from the ANAIS-112 Experiment,Phys. Rev. Lett.123(2019) 031301, [arXiv:1903.03973]

  15. [15]

    J. Amareet. al.,Annual modulation results from three-year exposure of ANAIS-112,Phys. Rev. D103 (2021) 102005, [arXiv:2103.01175]

  16. [16]

    Coarasaet

    I. Coarasaet. al.,ANAIS–112 three years data: a sensitive model independent negative test of the DAMA/LIBRA dark matter signal,Commun. Phys.7 (2024) 345, [arXiv:2404.17348]

  17. [17]

    Amar´ eet

    J. Amar´ eet. al.,Towards a Robust Model-Independent Test of the DAMA/LIBRA Dark Matter Signal: ANAIS-112 Results with Six Years of Data,Phys. Rev. Lett.135(2025) 051001, [arXiv:2502.01542]

  18. [18]

    Adhikariet

    G. Adhikariet. al.,An experiment to search for dark-matter interactions using sodium iodide detectors, Nature564(2018) 83–86, [arXiv:1906.01791]. [Erratum: Nature 566, E2 (2019)]

  19. [19]

    Adhikariet

    COSINE-100: G. Adhikariet. al.,Strong constraints from COSINE-100 on the DAMA dark matter results using the same sodium iodide target,Sci. Adv.7(2021) abk2699, [arXiv:2104.03537]

  20. [20]

    Adhikariet

    COSINE-100: G. Adhikariet. al.,Search for a Dark Matter-Induced Annual Modulation Signal in NaI(Tl) with the COSINE-100 Experiment,Phys. Rev. Lett.123 (2019) 031302, [arXiv:1903.10098]

  21. [21]

    Adhikariet

    COSINE-100: G. Adhikariet. al.,Three-year annual modulation search with COSINE-100,Phys. Rev. D106 (2022) 052005, [arXiv:2111.08863]

  22. [22]

    Carlinet

    N. Carlinet. al.,COSINE-100 Full Dataset Challenges the Annual Modulation Signal of DAMA/LIBRA, arXiv:2409.13226

  23. [23]

    Carlinet

    COSINE-100, ANAIS-112: N. Carlinet. al.,Combined Annual Modulation Dark Matter Search with COSINE-100 and ANAIS-112,arXiv:2503.19559

  24. [24]

    Bernabeiet

    R. Bernabeiet. al.,Further results from 7 DAMA/Libra-phase2 and perspectives,Nucl. Phys. Atom. Energy22(2021) 329–342

  25. [25]

    A. K. Drukier, K. Freese, and D. N. Spergel,Detecting Cold Dark Matter Candidates,Phys. Rev. D33(1986) 3495–3508

  26. [26]

    D. N. Spergel,The Motion of the Earth and the Detection of Wimps,Phys. Rev. D37(1988) 1353

  27. [27]

    Freese, J

    K. Freese, J. A. Frieman, and A. Gould,Signal Modulation in Cold Dark Matter Detection,Phys. Rev. D37(1988) 3388–3405

  28. [28]

    Bernabeiet

    R. Bernabeiet. al.,Recent Results from DAMA/LIBRA and Comparisons,Moscow Univ. Phys. Bull.77(2022) 291–300

  29. [29]

    Bernabeiet

    DAMA: R. Bernabeiet. al.,The DAMA/LIBRA apparatus,Nucl. Instrum. Meth. A592(2008) 297–315, [arXiv:0804.2738]

  30. [30]

    S. Kang, S. Scopel, and G. Tomar,Probing DAMA/LIBRA data in the full parameter space of WIMP effective models of inelastic scattering,Phys. Rev. D99(2019) 103019, [arXiv:1902.09121]

  31. [31]

    Amar´ eet

    J. Amar´ eet. al.,Performance of ANAIS-112 experiment after the first year of data taking,Eur. Phys. J. C79(2019) 228, [arXiv:1812.01472]

  32. [32]

    S. H. Leeet. al.,Measurements of low-energy nuclear recoil quenching factors for Na and I recoils in the NaI(Tl) scintillator,Phys. Rev. C110(2024) 014614, [arXiv:2402.15122]

  33. [33]

    Maltoni and T

    M. Maltoni and T. Schwetz,Testing the statistical compatibility of independent data sets,Phys. Rev. D68 (2003) 033020, [hep-ph/0304176]

  34. [34]

    M. J. Reidet. al.,Trigonometric Parallaxes of High Mass Star Forming Regions: the Structure and Kinematics of the Milky Way,Astrophys. J.783(2014) 130, [arXiv:1401.5377]

  35. [35]

    A. J. Deason, A. Fattahi,et. al.,The local high-velocity tail and the galactic escape speed,MNRAS485(2019) 3514–3526, [arXiv:1901.02016]

  36. [36]

    Bishara, J

    F. Bishara, J. Brod, B. Grinstein, and J. Zupan, DirectDM: a tool for dark matter direct detection, arXiv:1708.02678

  37. [37]

    J. Brod, A. Gootjes-Dreesbach, M. Tammaro, and J. Zupan,Effective Field Theory for Dark Matter Direct Detection up to Dimension Seven,JHEP10(2018) 065, [arXiv:1710.10218]. [38]GAMBITCollaboration: P. Athron, C. Bal´ azs,et. al., GAMBIT: The Global and Modular Beyond-the-Standard-Model Inference Tool, Eur. Phys. J. C77(2017) 784, [arXiv:1705.07908]. Addend...

  38. [40]

    Bringmann et

    GAMBIT Dark Matter Workgroup: T. Bringmann et. al.,DarkBit: A GAMBIT module for computing dark matter observables and likelihoods,Eur. Phys. J. C 77(2017) 831, [arXiv:1705.07920]. [41]GAMBITScanner Workgroup: G. D. Martinez, J. McKay,et. al.,Comparison of statistical sampling methods with ScannerBit, the GAMBIT scanning module,Eur. Phys. J. C77(2017) 761,...

  39. [42]

    Bishara, J

    F. Bishara, J. Brod, B. Grinstein, and J. Zupan,Chiral Effective Theory of Dark Matter Direct Detection, JCAP1702(2017) 009, [arXiv:1611.00368]

  40. [43]

    Bishara, J

    F. Bishara, J. Brod, B. Grinstein, and J. Zupan,From quarks to nucleons in dark matter direct detection, JHEP11(2017) 059, [arXiv:1707.06998]

  41. [44]

    Athronet

    GAMBIT: P. Athronet. al.,Global analyses of Higgs portal singlet dark matter models using GAMBIT,Eur. Phys. J. C79(2019) 38, [arXiv:1808.10465]

  42. [45]

    Bernabeiet

    R. Bernabeiet. al.,New limits on WIMP search with large-mass low-radioactivity NaI(Tl) set-up at Gran Sasso,Phys. Lett. B389(1996) 757–766

  43. [46]

    Barberioet

    SABRE: E. Barberioet. al.,Simulation and background characterisation of the SABRE South experiment: SABRE South Collaboration,Eur. Phys. J. C83(2023) 878, [arXiv:2205.13849]. APPENDIX We present additional details on the results of our fits in table III. 8 TABLE III.χ 2 values, p-values and tension statistics obtained by maximizing individual experimental...