Pith. sign in

REVIEW 3 major objections 6 minor 33 references

The only difference between coherent and incoherent dark-matter scattering is crystal-momentum conservation, which weakens once multiple phonons are produced.

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 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.5

2026-07-31 05:31 UTC pith:3RNCQH4E

load-bearing objection Clean 1D proof that coherent vs incoherent is just crystal-momentum conservation, with useful numerics inside the model—but the abstract overclaims “validation for realistic 3D crystals.” the 3 major comments →

arxiv 2607.24933 v1 pith:3RNCQH4E submitted 2026-07-27 hep-ph hep-ex

Coherence from interference: a solvable model of sub-GeV dark matter-nucleus scattering

Lynn Lin , Tongyan Lin , Momei Fang This is my paper
classification hep-ph hep-ex
keywords sub-GeV dark matterdynamic structure factormultiphonon excitationsincoherent approximationcrystal momentum conservationdirect detectionphonon scattering
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.

Low-threshold dark-matter detectors must model how a crystal responds when a light particle deposits energy: at low momentum the response is a single collective phonon, while at high momentum it becomes an ordinary nuclear recoil. Computing the intermediate multiphonon regime exactly is prohibitively expensive in three dimensions. This paper solves a one-dimensional lattice of N atoms exactly and shows that the sole distinction between the full (coherent) response and the much cheaper incoherent approximation is a delta-function that enforces crystal-momentum conservation. Once two or more phonons are created that constraint becomes weak, so the incoherent formula already reproduces the exact rates. A simple hybrid that keeps the coherent single-phonon piece at long wavelength and switches to the incoherent formula elsewhere matches the exact cross-sections to within roughly thirty percent over most of the sub-GeV mass range. The result supplies a practical, validated shortcut for interpreting upcoming phonon-based searches.

Core claim

In an exactly solvable one-dimensional N-site lattice the full dynamic structure factor and the incoherent approximation differ only by the presence of a crystal-momentum-conserving delta function. That constraint becomes a weak restriction on the available phonon phase space once n greater than or equal to 2 phonons are produced, so the incoherent approximation (and a hybrid long-wavelength-plus-incoherent scheme) reproduces the exact scattering rates to within about 30 percent for both massive and massless mediators over most of the sub-GeV window.

What carries the argument

The n-phonon dynamic structure factor of the N-site chain (Eqs. 4.20 and 4.23), which isolates crystal-momentum conservation as the unique difference between the coherent and incoherent expressions and thereby explains why multiphonon rates can be computed without interference terms.

Load-bearing premise

That quantitative error estimates taken from a monoatomic harmonic one-dimensional chain with a hand-imposed infrared cutoff can be treated as a faithful stand-in for the isotropic three-dimensional structure factor of real crystals.

What would settle it

Compute the full multiphonon structure factor for a realistic three-dimensional monoatomic or diatomic crystal at intermediate q and compare the integrated rates against the hybrid Inc+LW prediction; a systematic discrepancy larger than 30 percent would falsify the claimed validation.

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

If this is right

  • Multiphonon contributions to sub-GeV dark-matter rates in crystals can be evaluated with the computationally cheap incoherent formula once n greater than or equal to 2.
  • A hybrid scheme that retains only the coherent long-wavelength single-phonon piece already yields percent-to-thirty-percent accuracy for experimental cross-section limits.
  • The same momentum-conservation argument supplies a first-principles justification for stitching single-phonon and nuclear-recoil calculations across the transition region.
  • Detector projections that previously relied on uncontrolled multiphonon approximations can now quote a controlled theoretical uncertainty.

Where Pith is reading between the lines

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

  • The same weakening of momentum conservation should apply to optical branches and multi-atom unit cells, so the hybrid scheme is likely to remain accurate for polar targets once the long-wavelength optical matrix element is inserted by hand.
  • Because the one-dimensional density of states diverges at low frequency, the numerical errors quoted here are probably conservative upper bounds relative to three-dimensional crystals where soft modes are phase-space suppressed.
  • Anharmonic corrections, estimated small elsewhere, would first appear as a broadening of the multiphonon continuum rather than a revival of coherent interference, preserving the utility of the incoherent formula.

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 / 6 minor

Summary. The manuscript studies sub-GeV dark-matter scattering from a harmonic one-dimensional lattice for which the multiphonon dynamic structure factor can be calculated recursively and, in the relevant limits, analytically. A two-site model first illustrates how interference suppresses optical-mode excitation at small q and averages away as phonon multiplicity increases. For the N-site chain, the authors derive Eqs. (4.20) and (4.23), showing that the full and incoherent structure factors differ by enforcement of crystal-momentum conservation. Numerical structure factors and differential rates are then compared for full, incoherent, impulse, and hybrid incoherent-plus-long-wavelength prescriptions. The paper concludes that the hybrid scheme reproduces the exact model's cross sections to about 30% over most of the considered mass range and presents this as quantitative validation for realistic 3D crystals.

Significance. If the claims are restricted to the harmonic 1D model, this is a significant and useful benchmark: it gives an exact, non-fit derivation of the full-to-incoherent relation, exposes crystal-momentum conservation as the relevant constraint, supplies explicit multiphonon spectra and rates, and motivates a falsifiable hybrid approximation for expensive 3D calculations. Those analytic and numerical strengths merit publication. Their portability to realistic 3D targets, however, currently remains an assumption rather than a demonstrated quantitative result.

major comments (3)
  1. [Abstract; §5; §6] Abstract; §5 opening; §6. The claim of "quantitative validation" for realistic 3D crystals is not established by the calculations presented. All dR/dq, dR/dω, and σ(mχ) results insert the monoatomic, acoustic-only 1D S(q,ω) into the isotropic 3D rate in Eq. (5.1), while Eq. (4.23) is explicitly generalized only for a monoatomic isotropic lattice. Silicon has a two-atom basis, optical branches, anisotropic dispersions, and polarization-dependent matrix elements, and no full 3D benchmark is shown. Either restrict the conclusion to a 1D proof of principle or provide a 3D test using branch-resolved D(ω,k), basis structure, and polarization factors, demonstrating that the stated ~30% accuracy survives.
  2. [§4.1, Fig. 6; §5.1, Fig. 12; Eq. (5.1)] The 1 meV mode cutoff is said to represent the detector threshold, but Eq. (5.1) correctly applies ωth to the total deposited energy. Individual phonons below ωth can still contribute to detectable n≥2 events. Discarding all modes below 1 meV therefore changes the physical spectrum, and it does so differently in the full and incoherent calculations because they assign crystal momentum differently. The cutoff also regularizes a 1D IR divergence absent in 3D. The reported agreement and 0.2–1 MeV discrepancy window may consequently be regulator-dependent. Please separate the physical threshold from the IR regulator, include sub-threshold phonons in multiphonon states, and show convergence or sensitivity to several cutoffs.
  3. [Eqs. (4.20) and (4.23); §6] The identity is convincing for the 1D monoatomic Bravais lattice, but the unqualified statement that momentum conservation is the only difference between coherent and incoherent scattering does not follow for realistic crystals. In the general definition, dropping d≠d′ removes intra-unit-cell interference in addition to the inter-unit-cell sum that enforces q−Σiki=G; anisotropic q·e polarization weights also do not reduce to the scalar D(ω,k) used in Eq. (4.23). Please state the precise assumptions under which Eqs. (4.20) and (4.23) hold, and either derive the multi-atom/anisotropic generalization or qualify the claims about general 3D crystals.
minor comments (6)
  1. [Eqs. (4.13) and (4.20)] Clarify whether these deltas are Kronecker deltas at finite N or Dirac deltas in the continuum, including the associated factors of N and the finite-N broadening prescription used in the numerical plots.
  2. [Table 1 and §5] The stitching domains should be stated more explicitly. In particular, specify whether the n=1 contribution above qBZ is included through umklapp in the incoherent term or omitted, and confirm that there is no double counting at q=qBZ and q=2√(2mω̄).
  3. [Eq. (2.1) and §5.3] q0 is not defined when F̃(q)=q0²/q² is introduced, while §5.3 uses F̃(q)=(mχv0/q)². Please use one definition consistently.
  4. [§4.3, §5, Appendix A] Because the numerical error estimates are central, please collect N, maximum phonon order, q stitching, energy binning, cutoff, and averaging choices in one place. Supplying plotting scripts or machine-readable rate tables would substantially improve reproducibility.
  5. [Throughout] There are several typographical issues: "goal is use" should be "goal is to use" in §4.1; "rapidly calculations" should be "rapid calculations" in §4.4; "agree results" should be "agree with results" in §5.3; and "direct direction calculations" should be "direct detection calculations" in §6.
  6. [Appendix B] Figures 15 and 16 would benefit from complete captions identifying the mediator form factor, σp, N, threshold, and the meaning of any smoothing or averaging, rather than relying on the main-text definitions.

Circularity Check

0 steps flagged

No circularity: the coherent/incoherent identity and rate comparisons are derived from the lattice Hamiltonian and the definition of S(q,ω), not fitted or forced by self-citation.

full rationale

The paper’s load-bearing chain is self-contained. The dynamic structure factor is defined from the target operator OT(q) and the two-point correlation of displacements (Eqs. 2.2–2.8). The 2-site and N-site models are solved from their harmonic Hamiltonians; the full vs. incoherent difference is obtained by summing auto- and cross-correlations, which produces the crystal-momentum delta in Eqs. 4.12–4.13 and 4.20/4.23. That identity is an algebraic consequence of discrete translation invariance, not an input. Numerical S(q,ω), dR/dq, dR/dω, and σ(mχ) then compare the exact 1D object to the incoherent and Inc+LW schemes built from the same model. Prior work (including Ref. [13], which shares an author) supplies the incoherent baseline being tested and standard single-phonon formulas; it does not manufacture the agreement. External parameters (Si sound speed, 1 meV cutoff, Standard Halo Model) are not fitted to the target rates. Weaknesses of the 1D→3D portability are modeling assumptions, not circular reductions. Score 0; steps empty.

Axiom & Free-Parameter Ledger

4 free parameters · 5 axioms · 0 invented entities

The load-bearing content is a standard harmonic lattice plus the definition of the dynamic structure factor; almost everything else is either textbook solid-state physics or a modeling choice needed to turn the 1D solution into rate plots. No new particles or forces are introduced. The free parameters are conventional material and detector inputs, not fitted to manufacture the coherent/incoherent identity.

free parameters (4)
  • ω0 (acoustic scale) = ≈10 meV
    Set to cl/a ≈ 10 meV from silicon longitudinal sound speed so that energy scales match a real target; enters all dispersion and rate plots.
  • IR energy cutoff ωcut = 1 meV
    Hand-imposed 1 meV floor to regulate the 1D IR divergence of W(q) and D(ω)/ω; identified with detector threshold and changes low-q rates.
  • Two-site mode energies ω+, ω− = 20 meV, 60 meV
    Chosen as 20 and 60 meV to mimic Si acoustic/optical scales and to make frequencies commensurate for visualization; used only in the toy model of Sec. 3.
  • Finite N and multiphonon truncation / q-stitching boundaries = N~100 for rates; stitch at qBZ and 2√(2mω̄)
    N=50–2000 and the qBZ / 2√(2mω̄) boundaries in Table 1 are numerical choices that affect spike structure and where impulse approximation is glued on.
axioms (5)
  • domain assumption Harmonic lattice dynamics (bilinear springs only); anharmonic corrections neglected
    Stated in Sec. 2; authors cite Ref. [34] that anharmonicity is small in most of phase space, but the entire exact solution uses the harmonic expansion.
  • domain assumption Dynamic structure factor factorization of the DM-target rate (Eq. 1.1 / 2.2)
    Standard in the light-DM direct-detection literature the paper builds on; assumed throughout.
  • ad hoc to paper 1D monoatomic periodic chain is a faithful enough proxy that its S(q,ω) can replace the isotropic 3D structure factor inside 3D rate integrals
    Explicit modeling step at the opening of Sec. 5; required for all claimed experimental-rate validations.
  • domain assumption Standard Halo Model velocity distribution with canonical vesc, v0, vE
    Sec. 5; conventional astrophysical input for rate plots, not derived here.
  • domain assumption Impulse approximation replaces the structure factor for q > 2√(2mω̄)
    Table 1 and Sec. 2; standard high-q limit, used uniformly for full, incoherent, and Inc+LW curves.

pith-pipeline@v1.2.0-grok45-kimik3 · 25528 in / 3409 out tokens · 78218 ms · 2026-07-31T05:31:10.452594+00:00 · methodology

0 comments
read the original abstract

How do dark matter-nucleus interactions transition from the regimes of coherent scattering, where single phonons are produced, to that of individual nuclear recoils? Answering this question relies on understanding multiphonon excitations. Multiphonons are important for interpreting low-threshold direct detection experiments, yet are computationally prohibitive to compute. In this paper, we employ a 1D $N$-site crystal lattice model where dark matter scattering can be computed exactly. We show that the only difference between coherent and incoherent scattering is that conservation of crystal momentum is enforced in coherent scattering. The momentum conservation constraint becomes less important as more phonons are produced, yielding the transition to incoherent scattering. Using numerical calculations of the 1D structure factor, we also obtain quantitative validation of using an incoherent approximation to compute sub-GeV dark matter scattering in realistic 3D crystals.

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

33 extracted references · 30 linked inside Pith

  1. [3]

    Kahn and T

    Y. Kahn and T. Lin,Searches for light dark matter using condensed matter systems,Rept. Prog. Phys.85(2022) 066901 [2108.03239]

    Pith/arXiv arXiv 2022

  2. [4]

    Zurek,Dark Matter Candidates of a Very Low Mass,Ann

    K.M. Zurek,Dark Matter Candidates of a Very Low Mass,Ann. Rev. Nucl. Part. Sci.74 (2024) 287 [2401.03025]. [5]TESSERACTcollaboration,First Limits on Light Dark Matter Interactions in a Low Threshold Two Channel Athermal Phonon Detector from the TESSERACT Collaboration, 2503.03683. – 27 – Figure 16.dR/dωresult forN= 100. [6]CRESSTcollaboration,First result...

    Pith/arXiv arXiv 2024

  3. [7]

    Knapen, T

    S. Knapen, T. Lin, M. Pyle and K.M. Zurek,Detection of Light Dark Matter With Optical Phonons in Polar Materials,Phys. Lett.B785(2018) 386 [1712.06598]

    Pith/arXiv arXiv 2018

  4. [8]

    Griffin, S

    S. Griffin, S. Knapen, T. Lin and K.M. Zurek,Directional Detection of Light Dark Matter with Polar Materials,Phys. Rev.D98(2018) 115034 [1807.10291]

    Pith/arXiv arXiv 2018

  5. [9]

    Trickle, Z

    T. Trickle, Z. Zhang, K.M. Zurek, K. Inzani and S. Griffin,Multi-Channel Direct Detection of Light Dark Matter: Theoretical Framework,JHEP03(2020) 036 [1910.08092]

    Pith/arXiv arXiv 2020

  6. [10]

    Knapen, T

    S. Knapen, T. Lin and K.M. Zurek,Light Dark Matter in Superfluid Helium: Detection with Multi-excitation Production,Phys. Rev.D95(2017) 056019 [1611.06228]

    Pith/arXiv arXiv 2017

  7. [11]

    Campbell-Deem, P

    B. Campbell-Deem, P. Cox, S. Knapen, T. Lin and T. Melia,Multiphonon excitations from dark matter scattering in crystals,Phys. Rev. D101(2020) 036006 [1911.03482]

    Pith/arXiv arXiv 2020

  8. [12]

    Y. Kahn, G. Krnjaic and B. Mandava,Dark Matter Detection with Bound Nuclear Targets: The Poisson Phonon Tail,Phys. Rev. Lett.127(2021) 081804 [2011.09477]

    Pith/arXiv arXiv 2021

  9. [13]

    Campbell-Deem, S

    B. Campbell-Deem, S. Knapen, T. Lin and E. Villarama,Dark matter direct detection from the single phonon to the nuclear recoil regime,Phys. Rev. D106(2022) 036019 [2205.02250]

    Pith/arXiv arXiv 2022

  10. [14]

    S. Gori, S. Knapen, T. Lin, P. Munbodh and B. Suter,Spin-dependent scattering of sub-GeV dark matter: Models and constraints,Phys. Rev. D112(2025) 075019 [2506.11191]

    Pith/arXiv arXiv 2025

  11. [15]

    Schutz and K.M

    K. Schutz and K.M. Zurek,Detectability of Light Dark Matter with Superfluid Helium,Phys. Rev. Lett.117(2016) 121302 [1604.08206]

    Pith/arXiv arXiv 2016

  12. [16]

    Caputo, A

    A. Caputo, A. Esposito and A.D. Polosa,Light Dark Matter and Superfluid He-4 from EFT, J. Phys. Conf. Ser.1468(2020) 012060 [1911.07867]. – 28 –

    Pith/arXiv arXiv 2020

  13. [17]

    Acanfora, A

    F. Acanfora, A. Esposito and A.D. Polosa,Sub-GeV Dark Matter in Superfluid He-4: an Effective Theory Approach,Eur. Phys. J.C79(2019) 549 [1902.02361]

    Pith/arXiv arXiv 2019

  14. [18]

    Caputo, A

    A. Caputo, A. Esposito, F. Piccinini, A.D. Polosa and G. Rossi,Directional detection of light dark matter from three-phonon events in superfluid4He,Phys. Rev. D103(2021) 055017 [2012.01432]

    Pith/arXiv arXiv 2021

  15. [19]

    Knapen, J

    S. Knapen, J. Kozaczuk and T. Lin,Migdal Effect in Semiconductors,Phys. Rev. Lett.127 (2021) 081805 [2011.09496]

    Pith/arXiv arXiv 2021

  16. [20]

    Berghaus, A

    K.V. Berghaus, A. Esposito, R. Essig and M. Sholapurkar,The Migdal effect in semiconductors for dark matter with masses below∼100 MeV,JHEP01(2023) 023 [2210.06490]

    Pith/arXiv arXiv 2023

  17. [21]

    Berghaus, R

    K.V. Berghaus, R. Essig and M.H. McDuffie,The Migdal effect in Semiconductors for the Effective Field Theory of Dark Matter Direct Detection,2603.13052

    arXiv

  18. [22]

    M. Ibe, W. Nakano, Y. Shoji and K. Suzuki,Migdal Effect in Dark Matter Direct Detection Experiments,JHEP03(2018) 194 [1707.07258]

    Pith/arXiv arXiv 2018

  19. [23]

    Berghaus, R

    K.V. Berghaus, R. Essig, Y. Hochberg, Y. Shoji and M. Sholapurkar,Phonon background from gamma rays in sub-GeV dark matter detectors,Phys. Rev. D106(2022) 023026 [2112.09702]

    Pith/arXiv arXiv 2022

  20. [24]

    P. Cox, T. Melia and S. Rajendran,Dark matter phonon coupling,Phys. Rev. D100(2019) 055011 [1905.05575]

    Pith/arXiv arXiv 2019

  21. [25]

    Griffin, K

    S.M. Griffin, K. Inzani, T. Trickle, Z. Zhang and K.M. Zurek,Multichannel direct detection of light dark matter: Target comparison,Phys. Rev.D101(2020) 055004 [1910.10716]

    Pith/arXiv arXiv 2020

  22. [26]

    Caputo, A

    A. Caputo, A. Esposito, E. Geoffray, A.D. Polosa and S. Sun,Dark Matter, Dark Photon and Superfluid He-4 from Effective Field Theory,Phys. Lett. B802(2020) 135258 [1911.04511]

    Pith/arXiv arXiv 2020

  23. [27]

    Griffin, Y

    S.M. Griffin, Y. Hochberg, K. Inzani, N. Kurinsky, T. Lin and T.C. Yu,Silicon carbide detectors for sub-GeV dark matter,Phys. Rev. D103(2021) 075002 [2008.08560]

    Pith/arXiv arXiv 2021

  24. [28]

    Trickle, Z

    T. Trickle, Z. Zhang and K.M. Zurek,Effective field theory of dark matter direct detection with collective excitations,Phys. Rev. D105(2022) 015001 [2009.13534]

    Pith/arXiv arXiv 2022

  25. [29]

    Taufertshöfer, M

    N. Taufertshöfer, M. Garcia-Sciveres and S.M. Griffin,Proposal for broad-range directional detection of light dark matter in cryogenic ice,Phys. Rev. D110(2024) 103552 [2301.04778]

    Pith/arXiv arXiv 2024

  26. [30]

    Ashour and S.M

    O.A. Ashour and S.M. Griffin,Pressure-Tunable Targets for Light Dark Matter Direct Detection: The Case of Solid Helium,2409.02439

    Pith/arXiv arXiv

  27. [31]

    X.-X. Li, N. Valsan and Z. Zhang,Astrophysical Uncertainties in Sub-GeV Dark Matter Detection via Single Phonon Excitations,2606.04091

    Pith/arXiv arXiv

  28. [32]

    Coskuner, T

    A. Coskuner, T. Trickle, Z. Zhang and K.M. Zurek,Directional detectability of dark matter with single phonon excitations: Target comparison,Phys. Rev. D105(2022) 015010 [2102.09567]

    Pith/arXiv arXiv 2022

  29. [33]

    Knapen, J

    S. Knapen, J. Kozaczuk and T. Lin,DarkELF: A python package for dark matter scattering in dielectric targets,Phys. Rev. D105(2022) 015014 [2104.12786]

    Pith/arXiv arXiv 2022

  30. [34]

    Lin, C.-H

    T. Lin, C.-H. Shen, M. Sholapurkar and E. Villarama,Anharmonic effects in nuclear recoils from sub-GeV dark matter,Phys. Rev. D109(2024) 095020 [2309.10839]. – 29 –

    Pith/arXiv arXiv 2024

  31. [35]

    Stratman and T

    C. Stratman and T. Lin,Daily modulation of low-energy nuclear recoils from sub-GeV dark matter,Phys. Rev. D111(2025) 043030 [2411.03433]

    Pith/arXiv arXiv 2025

  32. [36]

    Squires,Introduction to the Theory of Thermal Neutron Scattering, Dover Publications, Inc., Mineola, New York (1996)

    G.L. Squires,Introduction to the Theory of Thermal Neutron Scattering, Dover Publications, Inc., Mineola, New York (1996)

    1996

  33. [37]

    Schober,An introduction to the theory of nuclear neutron scattering in condensed matter, Journal of Neutron Research17(2014) 109

    H. Schober,An introduction to the theory of nuclear neutron scattering in condensed matter, Journal of Neutron Research17(2014) 109. – 30 –

    2014