{"id":"45ef8c9e-5c8b-45f6-8d27-29eeadf9ebf2","arxiv_id":"2607.24933","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":4,"one_line_summary":"In an exactly solvable 1D lattice, coherent and incoherent DM-nucleus structure factors differ only by a crystal-momentum delta function that becomes unimportant for n≥2 phonons, validating hybrid Inc+LW rate calculations.","lead":"A solvable 1D crystal model shows that coherent and incoherent dark-matter scattering differ only by crystal-momentum conservation, which weakens as more phonons are made. That validates cheap multiphonon approximations used to interpret low-threshold direct-detection experiments.","discovery_kind":"extension","skeptic_critique":{"model":"moonshotai/kimi-k3","headline":"The exact 1D identity is sound, but the claimed quantitative validation for realistic crystals rests on a soft-phonon-dominated 1D DOS and an ad hoc 1 meV cutoff. That can distort the very multiphonon averaging being tested.","rationale":"This is the same load-bearing limitation identified by the reader, sharpened to the mechanism most likely to affect the quantitative conclusion: the 1D infrared density of states and cutoff alter the soft-phonon contribution precisely where momentum conservation is most restrictive. I do not see a problem with the analytic derivation of the coherent/incoherent difference, and the internal 1D numerics appear plausible. The concern is therefore not a correctness objection to the solvable model, but an external-validity objection to calling the results quantitative validation for realistic 3D crystals. A direct silicon n=2 comparison is feasible and tests the worst multiphonon case; higher orders should be more averaged if the paper’s mechanism is correct. The reader’s CONDITIONAL verdict already reflects this gap, so no verdict change is warranted. Released code would aid independent checking, but it is secondary to the dimensionality/realism question.","tokens_in":21755,"tokens_out":9293,"duration_ms":285384,"concrete_test":"Using tabulated silicon phonon eigenvalues and eigenvectors, compute the exact n=2 structure factor with crystal-momentum conservation and polarization factors, and the incoherent n=2 DOS convolution from the same spectrum. Direction-average q and insert both into Eq. (5.1) for a massive mediator at ωth=1 meV across the n=2 transition region mχ≈0.2–1 MeV, checking supercell convergence. If the n=2 contribution or integrated rate differs by more than ~30%, the 1D-based validation does not carry over; agreement within that tolerance would directly support the extrapolation at its least-averaged multiphonon order.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equations (4.20) and (4.23) convincingly isolate crystal-momentum conservation within the idealized harmonic model. The fragile step is the numerical portability of the resulting ~30% validation. In the 1D acoustic chain, D(ω) is finite as ω→0, so D(ω)/ω∼1/ω; in an isotropic 3D acoustic crystal, D(ω)∼ω² and D(ω)/ω∼ω. The authors remove all modes below ωcut=1 meV, although an experimental threshold applies to the total deposited energy, not to each constituent phonon. Soft modes therefore carry enhanced and regulator-dependent weight, especially for n=2 processes and for the forward-scattering-dominated massless-mediator case. Because the full and incoherent calculations assign different momenta to those soft modes, their agreement—and the location/size of the reported 0.2–1 MeV discrepancy window—could be sensitive to this 1D IR structure. Silicon’s optical branches and two-atom basis introduce additional untested interference effects. Thus the abstract’s claim of quantitative validation for realistic 3D crystals is not yet established, even though the 1D mechanism itself is well supported.","agreement_with_reader":"agree"},"referee_report":{"model":"moonshotai/kimi-k3","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.","tokens_in":22048,"tokens_out":5582,"duration_ms":37107,"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":[{"comment":"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.","section":"Abstract; §5; §6"},{"comment":"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.","section":"§4.1, Fig. 6; §5.1, Fig. 12; Eq. (5.1)"},{"comment":"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.","section":"Eqs. (4.20) and (4.23); §6"}],"minor_comments":[{"comment":"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.","section":"Eqs. (4.13) and (4.20)"},{"comment":"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ω̄).","section":"Table 1 and §5"},{"comment":"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.","section":"Eq. (2.1) and §5.3"},{"comment":"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.","section":"§4.3, §5, Appendix A"},{"comment":"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.","section":"Throughout"},{"comment":"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.","section":"Appendix B"}],"recommendation":"major_revision","confidential_remarks":"The exact 1D treatment and its physical interpretation are a good fit for the journal. My concern is mainly one of claim calibration: the manuscript currently presents a regulator-dependent 1D calculation as quantitative validation for realistic 3D crystals. A suitably qualified claim, together with an infrared-cutoff sensitivity study, may suffice; a genuine 3D comparison would make the practical case considerably stronger."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The useful core is simple and solid. They solve the multiphonon dynamic structure factor exactly on a 1D N-site harmonic chain and show, cleanly in Eqs. 4.20/4.23, that the only difference between the full (coherent) response and the incoherent approximation is the crystal-momentum delta. Once n≥2, that constraint is weak, so multiphonon rates average toward the incoherent result. The two-site warmup, the recursion in App. A, and the Inc+LW error budget inside the model are done carefully. That is real progress relative to Campbell-Deem et al. and the papers that already use the incoherent shortcut: we now have a transparent mechanism and a controlled 1D error estimate (~30% worst case, often much better).\n\nWhat the paper does less well is the leap in the abstract and Sec. 5. They treat the 1D S(q,ω) as a stand-in for an isotropic 3D crystal rate and call it “quantitative validation … in realistic 3D crystals.” The 1D acoustic DOS is finite as ω→0 (so D(ω)/ω ~ 1/ω), they impose a hard 1 meV cut on individual modes, and there are no optical branches or multi-atom cells. Soft-mode weight and the massless-mediator forward peak can be regulator-sensitive in ways that need not match 3D. The mechanism should survive; the precise percent-level error window and the 0.2–1 MeV bump may not. That is a soft spot on the claim, not on the 1D math.\n\nCitations look appropriate; circularity is low; no code is released, which is a minor practical minus. Significance is mid-subfield: it underwrites approximations already in use rather than changing what experiments search for.\n\nWho should read it: anyone computing or using multiphonon/incoherent rates for sub-GeV nuclear couplings. It deserves a serious referee. I would engage—cite the identity and the 1D error study, and discount the 3D-validation language until someone checks a real 3D DOS or a diatomic cell.","headline":"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.”","tokens_in":22518,"tokens_out":556,"would_cite":true,"duration_ms":16790,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"The only difference between coherent and incoherent dark-matter scattering is crystal-momentum conservation, which weakens once multiple phonons are produced.","keywords":["sub-GeV dark matter","dynamic structure factor","multiphonon excitations","incoherent approximation","crystal momentum conservation","direct detection","phonon scattering"],"falsifier":"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.","tokens_in":22386,"feed_emoji":"⚛️","tokens_out":936,"duration_ms":16530,"temperature":0.7,"pith_summary":"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.","feed_headline":"Coherent dark-matter scattering becomes incoherent after two phonons","feed_subtitle":"A solvable 1D lattice shows momentum conservation is the only difference, validating a cheap hybrid formula for detectors","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"fun_headline_variants":["1D lattice: DM scattering incoheres after two phonons","Momentum conservation alone splits coherent from incoherent DM scatter","Multiphonon phase space erases crystal momentum constraint","Exact 1D structure factor validates incoherent approx for sub-GeV DM","Hybrid formula matches exact rates once n≥2 phonons appear"],"cache_read_input_tokens":16512,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["1D lattice: DM scattering incoheres after two phonons","Momentum conservation alone splits coherent from incoherent DM scatter","Multiphonon phase space erases crystal momentum constraint","Exact 1D structure factor validates incoherent approx for sub-GeV DM","Hybrid formula matches exact rates once n≥2 phonons appear"]},"model":"grok-4.5","effort":"low","cost_usd":0.003212,"raw_usage":{"total_tokens":1061,"prompt_tokens":734,"num_sources_used":0,"completion_tokens":67,"cost_in_usd_ticks":32124000,"prompt_tokens_details":{"text_tokens":734,"audio_tokens":0,"image_tokens":0,"cached_tokens":128},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":260,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":734,"tokens_out":67,"duration_ms":5783,"temperature":1.0,"reasoning_tokens":260,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-31T05:31:10.452594+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"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.","supporting_citations":[],"review_version":1}