pith. machine review for the scientific record. sign in

arxiv: 2512.20825 · v2 · submitted 2025-12-23 · ✦ hep-ph

Recognition: 2 theorem links

· Lean Theorem

Irreducible Constraints on Hadronically Interacting Sub-GeV Dark Matter

Authors on Pith no claims yet

Pith reviewed 2026-05-16 19:52 UTC · model grok-4.3

classification ✦ hep-ph
keywords dark mattersub-GeV dark matterchiral effective theorybig bang nucleosynthesisfreeze-inmeson decaysnucleon scattering
0
0 comments X

The pith

Sub-GeV dark matter with hadronic couplings is ruled out above 10^{-36} cm² nucleon scattering by induced electromagnetic effects

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

The paper establishes conservative upper limits on dark-matter-nucleon scattering for masses below 1 GeV by working strictly inside low-energy chiral effective theory. Any dark matter that couples hadronically at leading order must develop electromagnetic couplings to photons or electrons at next-to-leading order. These induced interactions generate tight bounds from big bang nucleosynthesis and low-temperature freeze-in overproduction, while the leading hadronic operators are limited by meson decays. The combined constraints exclude both spin-independent and spin-dependent cross sections at or above 10^{-36} cm² across the entire keV to 100 MeV mass range, independent of any ultraviolet completion. These limits lie several orders of magnitude below existing astrophysical and cosmological bounds and directly shape the reach of planned low-mass direct detection experiments.

Core claim

Dark matter that interacts only hadronically at leading order in the chiral effective theory inevitably produces electromagnetic interactions at next-to-leading order. These electromagnetic couplings drive strong constraints from big bang nucleosynthesis and dark-matter overproduction via freeze-in at low temperatures, while the leading-order hadronic couplings face direct limits from meson decays. Taken together the limits rule out dark-matter-nucleon scattering cross sections of 10^{-36} cm² or larger for masses in the keV to 100 MeV window, without dependence on the details of the ultraviolet theory.

What carries the argument

Next-to-leading-order electromagnetic operators that are necessarily generated by leading-order hadronic dark-matter couplings inside chiral effective field theory

If this is right

  • Both spin-independent and spin-dependent dark-matter-nucleon cross sections are excluded above 10^{-36} cm² across the full keV-100 MeV range
  • Future low-mass direct detection experiments must reach sensitivities well below 10^{-36} cm² to retain discovery potential
  • The bounds hold independently of any specific high-energy completion of the dark-matter model
  • Cosmological and meson-decay limits dominate over prior astrophysical constraints by several orders of magnitude

Where Pith is reading between the lines

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

  • Models of light dark matter that rely on hadronic couplings to explain other anomalies would need additional suppression mechanisms to survive these limits
  • The effective-theory linking of hadronic and electromagnetic sectors at successive orders could be applied to other light-particle scenarios beyond dark matter
  • If a signal appears above the bound, it would require either breakdown of the chiral effective theory or a different leading interaction structure

Load-bearing premise

Dark matter interacts only hadronically at leading order in the chiral effective theory, with electromagnetic interactions appearing only at next-to-leading order, and the effective theory remains valid down to the momentum scales set by sub-GeV masses

What would settle it

A confirmed direct-detection signal of sub-GeV dark matter with nucleon scattering cross section above 10^{-36} cm², or a measurement of primordial light-element abundances that contradicts the predicted big-bang-nucleosynthesis effects from the induced electromagnetic couplings

Figures

Figures reproduced from arXiv: 2512.20825 by Avirup Ghosh, Matthew J. Dolan, Peter Cox.

Figure 1
Figure 1. Figure 1: FIG. 1. Feynman diagrams that contribute to the DM ther [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Feynman diagrams contributing to DM thermalisa [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4 [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Constraints on the DM-nucleon spin-dependent cross-section [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Constraints on the DM-nucleon spin-dependent cross-section [PITH_FULL_IMAGE:figures/full_fig_p011_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Feynman diagrams contributing to [PITH_FULL_IMAGE:figures/full_fig_p017_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Feynman diagrams contributing to [PITH_FULL_IMAGE:figures/full_fig_p019_8.png] view at source ↗
read the original abstract

We derive conservative upper limits on the dark-matter--nucleon scattering cross-section for sub-GeV mass dark matter. Working exclusively within the low-energy chiral effective theory, we derive bounds that are independent of the details of the dark matter interactions in the UV. Dark matter that interacts only hadronically at leading order also inevitably interacts with photons or electrons at next-to-leading-order. We show that these electromagnetic interactions lead to strong constraints from big bang nucleosynthesis and over-production of dark matter via freeze-in at low temperatures, while the leading-order hadronic couplings face stringent constraints from meson decays. Combining these constraints, we rule out both spin-independent and spin-dependent dark-matter--nucleon scattering cross-sections $\gtrsim 10^{-36}\,{\rm cm}^2$ for dark matter masses in the keV - 100 MeV range. These bounds are several orders of magnitude stronger than the existing constraints from astrophysics and cosmology and have significant implications for future low-mass direct detection experiments.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper derives UV-independent upper limits on sub-GeV dark-matter–nucleon scattering cross sections by working in chiral effective theory. Dark matter is assumed to couple hadronically at leading order (LO), which generates electromagnetic couplings at next-to-leading order (NLO). LO hadronic operators are constrained by meson decays, while NLO electromagnetic operators are constrained by BBN and low-temperature freeze-in. Combining these yields the claim that both spin-independent and spin-dependent DM-nucleon cross sections ≳ 10^{-36} cm² are excluded for DM masses in the keV–100 MeV range.

Significance. If the central result holds, the work supplies conservative, model-independent bounds on light hadronic dark matter that are several orders of magnitude stronger than existing astrophysical and cosmological limits. The approach relies on standard chiral EFT power counting, BBN calculations, and freeze-in production, and the bounds would have direct implications for the sensitivity targets of future low-mass direct-detection experiments.

major comments (2)
  1. [§2] §2 (chiral Lagrangian): The central claim requires that NLO electromagnetic operators (DM-photon or DM-electron) cannot be parametrically suppressed relative to the LO hadronic couplings (DM-nucleon or DM-pion) while preserving the EFT validity. The manuscript states that electromagnetic interactions appear “only at next-to-leading order,” but does not demonstrate whether the NLO Wilson coefficients are fixed by the LO ones via matching or remain independent counterterms that can be tuned small. If independent tuning is allowed, the combined bound on the low-energy DM-nucleon cross section can be evaded. An explicit power-counting argument or matching calculation showing a lower bound on the NLO coefficients is needed.
  2. [§4–5] §4–5 (BBN and freeze-in): The quantitative translation from NLO electromagnetic couplings to the 10^{-36} cm² exclusion relies on specific freeze-in rates and BBN abundance limits. Without the explicit expressions for the NLO-induced DM production cross sections or the error propagation from the chiral scale, it is not possible to verify that the bound remains robust when the NLO coefficients are varied within the EFT uncertainty band.
minor comments (2)
  1. [Table 1] Table 1: the quoted meson-decay limits should include the precise branching-ratio inputs and the assumed form-factor parametrization used to convert them into bounds on the LO couplings.
  2. [Eq. (X)] Eq. (X) (cross-section formula): the relation between the low-energy DM-nucleon cross section and the chiral coefficients should be written explicitly, including the kinematic factors for both SI and SD cases.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We address each major point below and will incorporate clarifications to strengthen the presentation.

read point-by-point responses
  1. Referee: [§2] §2 (chiral Lagrangian): The central claim requires that NLO electromagnetic operators (DM-photon or DM-electron) cannot be parametrically suppressed relative to the LO hadronic couplings (DM-nucleon or DM-pion) while preserving the EFT validity. The manuscript states that electromagnetic interactions appear “only at next-to-leading order,” but does not demonstrate whether the NLO Wilson coefficients are fixed by the LO ones via matching or remain independent counterterms that can be tuned small. If independent tuning is allowed, the combined bound on the low-energy DM-nucleon cross section can be evaded. An explicit power-counting argument or matching calculation showing a lower bound on the NLO coefficients is needed.

    Authors: In the chiral EFT we employ, standard power counting ensures that LO hadronic operators (DM-pion and DM-nucleon) generate NLO electromagnetic operators through tree-level matching to the electromagnetic current and one-loop pion exchange, with the NLO Wilson coefficients fixed in terms of the LO couplings up to O(1) factors set by the chiral scale Λ_χ. Independent suppression of the NLO coefficients would require additional fine-tuning or new light degrees of freedom that lie outside the EFT validity range, which is inconsistent with our UV-independent approach. We will revise §2 to include an explicit schematic matching calculation and power-counting argument demonstrating the lower bound on the NLO coefficients. revision: yes

  2. Referee: [§4–5] §4–5 (BBN and freeze-in): The quantitative translation from NLO electromagnetic couplings to the 10^{-36} cm² exclusion relies on specific freeze-in rates and BBN abundance limits. Without the explicit expressions for the NLO-induced DM production cross sections or the error propagation from the chiral scale, it is not possible to verify that the bound remains robust when the NLO coefficients are varied within the EFT uncertainty band.

    Authors: We agree that explicit expressions improve verifiability. In the revised manuscript we will add the full analytic expressions for the NLO-induced DM production cross sections (both DM-photon and DM-electron channels) in §4, together with the freeze-in rate integrals. In §5 we will include a dedicated error-propagation analysis showing that O(1) variations in the NLO coefficients around the chiral-scale expectation leave the 10^{-36} cm² bound intact as a conservative order-of-magnitude limit. The BBN constraints are obtained from standard public codes with the induced electromagnetic interactions; we will reference the relevant input files and parameter ranges. revision: yes

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper's central bounds combine standard BBN and freeze-in limits on NLO electromagnetic operators with meson-decay limits on LO hadronic operators, all within chiral EFT. No equations reduce the final cross-section limit to a parameter fitted from the same data or to a self-citation chain. The claim that NLO EM interactions arise inevitably is presented as a derivation from the EFT structure rather than an assumption that loops back to the target result. The derivation remains self-contained against external benchmarks (BBN, meson lifetimes) without self-definitional or fitted-input reductions.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the validity of chiral effective theory for sub-GeV momentum transfers and the assumption that dark matter has no direct leading-order coupling to photons or electrons.

axioms (1)
  • domain assumption Chiral effective theory accurately captures low-energy hadronic interactions independent of UV completion
    Invoked to derive bounds that do not depend on high-energy details

pith-pipeline@v0.9.0 · 5472 in / 1304 out tokens · 34223 ms · 2026-05-16T19:52:17.869562+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

70 extracted references · 70 canonical work pages · 25 internal anchors

  1. [1]

    DM interacts only hadronically at leading order. 1

  2. [2]

    The mediator(s) between the DM and SM have mass greater than 2 O(100) MeV

  3. [3]

    (Motivated by aU(1) orZ 2 stabilising symmetry.)

    DM interacts with the SM via an operator that is bilinear in the DM fieldχ. (Motivated by aU(1) orZ 2 stabilising symmetry.)

  4. [4]

    The DM fieldχis a singlet underSU(3) c ×U(1) Q. Schematically, we consider DM interactions of the form L ⊃ cαβ Λn Oα SMOβ χ ,(1) whereO α SM are local, gauge-singlet operators consisting of quark and/or gluon fields and the operatorsO β χ are bi- linear inχand can be non-local. Thec αβ are dimension- less coefficients and Λ parametrises the scale of the U...

  5. [5]

    In fact, sub-MeV mass states that were in equilibrium at any time after the QCD phase transition but before the end of BBN are essentially excluded

    Big Bang Nucleosynthesis The success of SM BBN tightly constraints the abun- dances of any additional relativistic degrees of freedom during neutrino decoupling and nucleosynthesis. In fact, sub-MeV mass states that were in equilibrium at any time after the QCD phase transition but before the end of BBN are essentially excluded. More specifically, an elec...

  6. [6]

    Once again, this allows us to remain insen- sitive to the early cosmological history and the specifics of the DM production mechanism

    Irreducible DM abundance We do not require thatχsaturates the DM relic abun- dance (except when showing existing bounds from direct detection). Once again, this allows us to remain insen- sitive to the early cosmological history and the specifics of the DM production mechanism. We do, however, re- quire that any DM produced at low temperatures (below 10 M...

  7. [7]

    In particular, the decay rates ofπ,K, andBmesons to final states with miss- ing energy are all tightly constrained

    Meson decays Meson decays can provide strong constraints on the hadronic interactions of DM. In particular, the decay rates ofπ,K, andBmesons to final states with miss- ing energy are all tightly constrained. This is espe- cially true for the FCNC decaysK→π+ invisible, B→K+ invisible, andB→π+ invisible. In this work, we do not consider bounds from invisib...

  8. [8]

    First, there are those observables 5 that require the specification of a Lagrangian valid at higher energies, beyond a few hundred MeV

    Other constraints There are many other observables that can be used to constrain the DM parameter space, but which we do not consider here. First, there are those observables 5 that require the specification of a Lagrangian valid at higher energies, beyond a few hundred MeV. These in- clude searches for missing energy at the Large Hadron Collider, direct ...

  9. [9]

    The former process generally dominates, since the latter is suppressed by an additional factor ofαand the three-body phase space

    Vector operator For the vector operator, DM interacts with the SM plasma atT≲10 MeV via the processese +e− →χ¯χand/or γγ→γχ¯χ. The former process generally dominates, since the latter is suppressed by an additional factor ofαand the three-body phase space. However,γγ→γχ¯χbecomes the dominant process for the case of universal couplings, C V =1. In the foll...

  10. [10]

    Axial-vector operator For the axial-vector operator, the DM interacts with the SM plasma viaγγ→χ¯χ, which proceeds through the Feynman diagram shown in the left panel of fig. 4. The amplitude for this process is readily obtained from the DM-meson Lagrangian in eq. (14) and is given as follows: iMγγ→χ¯χ= iα πΛ2 ϵµναβ k1µk2νϵα(k1)ϵβ(k2) s−m 2 π0 +im π0Γπ0 P...

  11. [11]

    This process proceeds via the Feynman diagram shown in the left panel of fig

    Pseudoscalar operator For DM interactions with the quark pseudoscalar operator,γγ→χ¯χis the leading process for DM thermalisation and freeze-in production atT≲10 MeV. This process proceeds via the Feynman diagram shown in the left panel of fig. 4 and its amplitude can be readily obtained from the leading order Chiral Lagrangian in eq. (19): iMγγ→χ¯χ= iα π...

  12. [12]

    We provide the decay rates for the corresponding decays to DM,K + →π +χ¯χ andπ 0 →γχ¯χ, in the following subsections

    Vector operator For the vector operator, the relevant flavour physics constraints come from the upper-limits on the branching ratios ofK + →π +ν¯ν[28] andπ 0 →γν¯ν[32]. We provide the decay rates for the corresponding decays to DM,K + →π +χ¯χ andπ 0 →γχ¯χ, in the following subsections. a.K + →π +χ¯χ Under the assumption of flavour-diagonal Wilson coeffici...

  13. [13]

    C A u +C A d + 4 (C A u −C A d ) 1− 16B0L4 f2 X q mq + 16B0L5 f2 (C A u mu −C A d md) ! q2 12 −m 2 π q2 12 −m 2 π0 # +g S 8 X q C A q , (C12) CA π =g 8

    Axial-vector operator For the axial-vector operator, the relevant flavour physics constraints come from the measurements of, or upper- limits on, the branching ratios ofK + →π +ν¯ν[28],π 0 →ν¯ν[33] andη→ν¯ν[34]. In the following subsections we derive the decay rates for the corresponding decays to DM. 19 γ γ π0 π0, η ¯χ χ K + π+ χ ¯χ ¯χ χ ¯χ χ K + π+ π0 I...

  14. [14]

    There are also constraints from the decayK + →π +χ¯χ

    Pseudoscalar operator Similar to the axial-vector operator, for the pseudoscalar operator, the leading order Chiral Lagrangian (19) leads to the decaysπ 0 →χ¯χandη→χ¯χ. There are also constraints from the decayK + →π +χ¯χ. a.K + →π +χ¯χ The relevant terms in the ∆S= 1 Lagrangian (7) for the pseudoscalar operator are LLO ∆S=1 =− √ 2GF VudV ∗ usf2g8 (∂µK+)(...

  15. [15]

    Essiget al., inSnowmass 2021(2022) arXiv:2203.08297 [hep-ph]

    R. Essiget al., inSnowmass 2021(2022) arXiv:2203.08297 [hep-ph]

  16. [16]

    E. O. Nadler, V. Gluscevic, K. K. Boddy, and R. H. Wechsler, Astrophys. J. Lett.878, 32 (2019), [Erra- tum: Astrophys.J.Lett. 897, L46 (2020), Erratum: As- trophys.J. 897, L46 (2020)], arXiv:1904.10000 [astro- ph.CO]

  17. [17]

    E. O. Nadleret al.(DES), Phys. Rev. Lett.126, 091101 (2021), arXiv:2008.00022 [astro-ph.CO]

  18. [18]

    Maamari, V

    K. Maamari, V. Gluscevic, K. K. Boddy, E. O. Nadler, and R. H. Wechsler, Astrophys. J. Lett.907, L46 (2021), arXiv:2010.02936 [astro-ph.CO]

  19. [19]

    M. A. Buen-Abad, R. Essig, D. McKeen, and Y.-M. Zhong, Phys. Rept.961, 1 (2022), arXiv:2107.12377 [astro-ph.CO]

  20. [20]

    Constraining Dark Matter-Baryon Scattering with Linear Cosmology

    C. Dvorkin, K. Blum, and M. Kamionkowski, Phys. Rev. D89, 023519 (2014), arXiv:1311.2937 [astro-ph.CO]

  21. [21]

    W. L. Xu, C. Dvorkin, and A. Chael, Phys. Rev. D97, 103530 (2018), arXiv:1802.06788 [astro-ph.CO]

  22. [22]

    K. K. Rogers, C. Dvorkin, and H. V. Peiris, Phys. Rev. Lett.128, 171301 (2022), arXiv:2111.10386 [astro- ph.CO]

  23. [23]

    X.-l. Chen, S. Hannestad, and R. J. Scherrer, Phys. Rev. D65, 123515 (2002), arXiv:astro-ph/0202496

  24. [24]

    Constraints on scattering of keV--TeV dark matter with protons in the early Universe

    V. Gluscevic and K. K. Boddy, Phys. Rev. Lett.121, 081301 (2018), arXiv:1712.07133 [astro-ph.CO]

  25. [25]

    T. R. Slatyer and C.-L. Wu, Phys. Rev. D98, 023013 (2018), arXiv:1803.09734 [astro-ph.CO]

  26. [26]

    S. Das, R. Mondol, A. Singh, and R. Laha, (2025), arXiv:2511.02906 [hep-ph]

  27. [27]

    A. He, M. M. Ivanov, R. An, T. Driskell, and V. Glusce- vic, JCAP05(2025), 087, arXiv:2502.02636 [astro- ph.CO]

  28. [28]

    M. J. Dolan, F. Kahlhoefer, C. McCabe, and K. Schmidt- Hoberg, JHEP03(2015), 171, [Erratum: JHEP 07, 103 (2015)], arXiv:1412.5174 [hep-ph]

  29. [29]

    Probing Light Thermal Dark-Matter With a Higgs Portal Mediator

    G. Krnjaic, Phys. Rev. D94, 073009 (2016), arXiv:1512.04119 [hep-ph]

  30. [30]

    Light Dark Matter: Models and Constraints

    S. Knapen, T. Lin, and K. M. Zurek, Phys. Rev. D96, 115021 (2017), arXiv:1709.07882 [hep-ph]

  31. [31]

    S. Gori, S. Knapen, T. Lin, P. Munbodh, and B. Suter, Phys. Rev. D112, 075019 (2025), arXiv:2506.11191 [hep- ph]. 22

  32. [32]

    P. Cox, M. J. Dolan, and J. Wood, Phys. Rev. D112, 115021 (2025), arXiv:2408.12144 [hep-ph]

  33. [33]

    Gasser and H

    J. Gasser and H. Leutwyler, Nucl. Phys. B250, 465 (1985)

  34. [34]

    Chiral Effective Theory of Dark Matter Direct Detection

    F. Bishara, J. Brod, B. Grinstein, and J. Zupan, JCAP 02(2017), 009, arXiv:1611.00368 [hep-ph]

  35. [35]

    Matrix element analyses of dark matter scattering and annihilation

    J. Kumar and D. Marfatia, Phys. Rev. D88, 014035 (2013), arXiv:1305.1611 [hep-ph]

  36. [36]

    A Lower Bound on the Mass of Cold Thermal Dark Matter from Planck

    C. Boehm, M. J. Dolan, and C. McCabe, JCAP08 (2013), 041, arXiv:1303.6270 [hep-ph]

  37. [37]

    P. F. Depta, M. Hufnagel, K. Schmidt-Hoberg, and S. Wild, JCAP04, 029, arXiv:1901.06944 [hep-ph]

  38. [38]

    Sabti, J

    N. Sabti, J. Alvey, M. Escudero, M. Fairbairn, and D. Blas, JCAP01(2020), 004, arXiv:1910.01649 [hep- ph]

  39. [39]

    Sabti, J

    N. Sabti, J. Alvey, M. Escudero, M. Fairbairn, and D. Blas, JCAP08(2021), A01, arXiv:2107.11232 [hep- ph]

  40. [40]

    What is the lowest possible reheating temperature?

    S. Hannestad, Phys. Rev. D70, 043506 (2004), arXiv:astro-ph/0403291

  41. [41]

    Barbieri, T

    N. Barbieri, T. Brinckmann, S. Gariazzo, M. Lattanzi, S. Pastor, and O. Pisanti, Phys. Rev. Lett.135, 181003 (2025), arXiv:2501.01369 [astro-ph.CO]

  42. [42]

    Cortina Gilet al.(NA62), JHEP02(2025), 191, arXiv:2412.12015 [hep-ex]

    E. Cortina Gilet al.(NA62), JHEP02(2025), 191, arXiv:2412.12015 [hep-ex]

  43. [43]

    J. K. Ahnet al.(KOTO), Phys. Rev. Lett.134, 081802 (2025), arXiv:2411.11237 [hep-ex]

  44. [44]

    Pich and A

    A. Pich and A. Rodr´ ıguez-S´ anchez, JHEP06(2021), 005, arXiv:2102.09308 [hep-ph]

  45. [45]

    Kaon Decays in the Standard Model

    V. Cirigliano, G. Ecker, H. Neufeld, A. Pich, and J. Por- toles, Rev. Mod. Phys.84, 399 (2012), arXiv:1107.6001 [hep-ph]

  46. [46]

    Search for production of an invisible dark photon in $\pi^0$ decays

    E. Cortina Gilet al.(NA62), JHEP05(2019), 182, arXiv:1903.08767 [hep-ex]

  47. [47]

    Search forπ 0 decays to invisible particles,

    E. Cortina Gilet al.(NA62), JHEP02(2021), 201, arXiv:2010.07644 [hep-ex]

  48. [48]

    Y. M. Andreevet al.(NA64), Phys. Rev. Lett.133, 121803 (2024), arXiv:2406.01990 [hep-ex]

  49. [49]

    Ballesteros, M

    G. Ballesteros, M. A. G. Garcia, and M. Pierre, JCAP 03(2021), 101, arXiv:2011.13458 [hep-ph]

  50. [50]

    Stellar cooling bounds on new light particles: plasma mixing effects

    E. Hardy and R. Lasenby, JHEP02, 033, arXiv:1611.05852 [hep-ph]

  51. [51]

    Lella, P

    A. Lella, P. Carenza, G. Co’, G. Lucente, M. Giannotti, A. Mirizzi, and T. Rauscher, Phys. Rev. D109, 023001 (2024), arXiv:2306.01048 [hep-ph]

  52. [52]

    Bottaro, A

    S. Bottaro, A. Caputo, G. Raffelt, and E. Vitagliano, JCAP07, 071, arXiv:2303.00778 [hep-ph]

  53. [53]

    Hardy, A

    E. Hardy, A. Sokolov, and H. Stubbs, JHEP04, 013, arXiv:2410.17347 [hep-ph]

  54. [54]

    Pich, Rept

    A. Pich, Rept. Prog. Phys.58, 563 (1995), arXiv:hep- ph/9502366

  55. [55]

    Mesonic low-energy constants

    J. Bijnens and G. Ecker, Ann. Rev. Nucl. Part. Sci.64, 149 (2014), arXiv:1405.6488 [hep-ph]

  56. [56]

    Tools for model-independent bounds in direct dark matter searches

    M. Cirelli, E. Del Nobile, and P. Panci, JCAP10(2013), 019, arXiv:1307.5955 [hep-ph]

  57. [57]

    T. R. Slatyer, Phys. Rev. D93, 023527 (2016), arXiv:1506.03811 [hep-ph]

  58. [58]

    Aggarwalet al.(DAMIC-M), Phys

    K. Aggarwalet al.(DAMIC-M), Phys. Rev. Lett.135, 071002 (2025), arXiv:2503.14617 [hep-ex]

  59. [59]

    Adariet al.(SENSEI), Phys

    P. Adariet al.(SENSEI), Phys. Rev. Lett.134, 011804 (2025), arXiv:2312.13342 [astro-ph.CO]

  60. [60]

    Huanget al.(PandaX), Phys

    D. Huanget al.(PandaX), Phys. Rev. Lett.131, 191002 (2023), arXiv:2308.01540 [hep-ex]

  61. [61]

    Revisiting Scalar and Pseudoscalar Couplings with Nucleons

    H.-Y. Cheng and C.-W. Chiang, JHEP07(2012), 009, arXiv:1202.1292 [hep-ph]

  62. [62]

    Wang, K.-Y

    W. Wang, K.-Y. Wu, L. Wu, and B. Zhu, Nucl. Phys. B 983, 115907 (2022), arXiv:2112.06492 [hep-ph]

  63. [63]

    Di Luzio, G

    L. Di Luzio, G. Levati, P. Paradisi, and X. Ponce D´ ıaz, JHEP12(2023), 016, arXiv:2308.05215 [hep-ph]

  64. [64]

    Y. Kahn, G. Krnjaic, S. Mishra-Sharma, and T. M. P. Tait, JHEP05(2017), 002, arXiv:1609.09072 [hep-ph]

  65. [65]

    Axial Vector $Z'$ and Anomaly Cancellation

    A. Ismail, W.-Y. Keung, K.-H. Tsao, and J. Unwin, Nucl. Phys. B918, 220 (2017), arXiv:1609.02188 [hep-ph]

  66. [66]

    S. Ipek, D. McKeen, and A. E. Nelson, Phys. Rev. D90, 055021 (2014), arXiv:1404.3716 [hep-ph]

  67. [67]

    Pseudoscalar Portal Dark Matter

    A. Berlin, S. Gori, T. Lin, and L.-T. Wang, Phys. Rev. D92, 015005 (2015), arXiv:1502.06000 [hep-ph]

  68. [68]

    Sch´ erer and K

    H. Sch´ erer and K. Schutz, JHEP11(2024), 139, arXiv:2405.18466 [hep-ph]

  69. [69]

    Making dark matter out of light: freeze-in from plasma effects

    C. Dvorkin, T. Lin, and K. Schutz, Phys. Rev. D 99, 115009 (2019), [Erratum: Phys.Rev.D 105, 119901 (2022)], arXiv:1902.08623 [hep-ph]

  70. [70]

    Song, Phys

    C. Song, Phys. Lett. B329, 312 (1994)