{"id":"330b4bfd-63a6-42e7-a371-b5e66ffd1f82","arxiv_id":"2607.15098","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Resonant excitation of nucleons into Δ(1232) during Earth passage is a non-negligible attenuation channel for boosted dark matter at E_χ ≈ 1–2 GeV, lowering the PandaX-4T upper bound on σ̄_n in the heavy-mediator regime.","lead":"This paper adds a missing nuclear reaction channel — excitation of the Δ(1232) nucleon resonance — to calculations of how fast, boosted dark matter is slowed or absorbed by Earth before reaching underground detectors. In the heavy-mediator regime the new channel is non-negligible near 1–2 GeV and shifts the PandaX-4T cross-section bound downward.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"In-medium Δ self-energy and unspecified QES/RES matching make the claimed O(10–50%) RES-induced boundary shift comparable to unquantified nuclear-model uncertainty.","rationale":"The paper is a legitimate extension of the Su/Wu/Zhu attenuation framework, and the derivation of the RES response is internally consistent, using standard nuclear machinery (NuWro spectral function, MAID2007 form factors, Rarita–Schwinger projector). I do not see an algebraic error in the central RES formula or an internal contradiction. The most load-bearing weakness is not the existence of the channel but the accuracy of its magnitude: the claimed exclusion-boundary shift is 10–50%, and nuclear-medium corrections to Δ production in nuclei are known to be of this size. If the free-width impulse-approximation treatment is off by tens of percent, the headline boundary shift is not quantitatively robust. This is exactly the condition the reader identified. The unspecified QES/RES and RES/DIS matching is a second aspect of the same uncertainty: without a stated W/Q² partition, the effective RES strength could be over- or under-counted. I also note the Sec. IV.A opening paragraph appears mislabeled (it describes the single-scattering model under the straight-line heading), but this is editorial and does not affect the physics. The concern is substantive yet addressable; it does not negate the qualitative claim that RES matters at GeV energies, so the reader's CONDITIONAL verdict remains appropriate.","tokens_in":12917,"tokens_out":15403,"duration_ms":174557,"concrete_test":"Recompute σ_w^RES for Fe at Eχ=1.5 GeV (m_V=10 GeV, m_χ=1 MeV) using Eq. (13) and Eq. (A15), replacing the free Breit–Wigner ΓΔ=117 MeV with a density-dependent Δ spectral function from an established nuclear model (e.g., Oset–Salcedo or GiBUU) that includes Pauli blocking and ΔN→NN absorption, evaluated along the same Earth trajectories. If the resulting σ_w^RES changes by less than ~20% relative to the free-width result, the Fig. 6 boundary shift is robust to in-medium effects; if it changes by more than ~50%, the claimed O(10–50%) shift is within nuclear-model uncertainty.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that including Δ(1232) RES lowers the PandaX-4T upper exclusion boundary by O(10–50%) at Eχ≈1–2 GeV—depends on the RES energy-loss and total cross sections computed in Eqs. (10)–(12)/(A15). These use a bound-nucleon spectral function in impulse approximation and a free Δ(1232) Breit–Wigner with ΓΔ=117 MeV. The paper explicitly states that the full W-dependence is kept only in this free Breit–Wigner (Sec. III, after Eq. 12), which is precisely where nuclear-medium effects enter. At densities relevant to Earth's interior, the Δ self-energy is known from neutrino/electron-nucleus studies to be modified at the tens-of-percent level: collisional broadening, Pauli blocking of Δ→πN decay, and ΔN→NN absorption all alter the inclusive Δ response. No such corrections are included or estimated. In addition, the paper does not specify the W/Q² matching used to separate RES from the QES and DIS channels taken from Ref. [60]; if the DIS prescription already includes low-W resonance contributions, adding Eq. (12) double-counts them. Because the reported RES-induced shift is of the same order as these unquantified nuclear-model uncertainties, the quantitative central claim—and the precise placement of the boundary in Fig. 6—is not yet robust, even though the existence of an additional RES attenuation channel is plausible.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper extends the Earth-attenuation treatment of boosted dark matter (BDM) by adding the resonant excitation of bound nucleons to the Δ(1232) resonance to the previously considered elastic, quasi-elastic, and deep-inelastic channels. The RES cross section is computed in the impulse approximation with a nuclear spectral function and a Breit-Wigner Δ width, and is then included in two attenuation models (continuous energy loss and single-scattering absorption) to compute the boosted-DM flux at the PandaX-4T detector. The authors report that, in the heavy-mediator regime, RES is sizable for E_χ ≈ 1–2 GeV in Fe and O, and that including it lowers the attenuation-induced upper boundary of the 90% C.L. exclusion region on the DM–nucleon cross section. The paper includes a detailed derivation of the nuclear response in Appendix A.","tokens_in":13256,"tokens_out":12415,"duration_ms":127694,"significance":"If the quantitative result is robust, the paper establishes that resonance production should be included in boosted-DM attenuation analyses for GeV-scale incoming energies, and it provides a concrete, detailed computation of the N→Δ(1232) response using a spectral-function impulse approximation and empirical form factors. The use of an external benchmark (PandaX-4T) and the explicit separation of the nuclear response into ES/QES/RES/DIS channels make the result directly usable in future BDM studies. The main value of the paper is phenomenological: it identifies a previously neglected channel and quantifies its effect on exclusion boundaries. The authors are transparent about the approximations used, and the derivation is largely self-contained.","major_comments":[{"comment":"The RES cross section is computed with the free Δ(1232) width Γ_Δ = 117 MeV and no in-medium modification of the resonance. At Earth-core densities, Δ self-energy effects — collisional broadening, Pauli blocking of Δ→πN decay, and ΔN→NN absorption — are known from neutrino/electron-nucleus studies to modify the inclusive Δ response at the tens-of-percent level. The reported RES-induced shift of the exclusion boundary in Fig. 6 is of the same order as these unquantified nuclear-model uncertainties. The qualitative existence of an additional RES attenuation channel is plausible, but the quantitative central claim is not yet robust. The authors should either model a density-dependent Δ width/self-energy or provide a sensitivity scan (e.g., varying Γ_Δ by ±30% or adding a collisional width) and show how the exclusion boundary shifts within that range.","section":"Sec. III, Eq. (12) and Appendix A (A15)"},{"comment":"The QES and DIS differential cross sections are taken from Ref. [60], but the paper does not specify the kinematic cuts that separate QES, RES, and DIS. If the DIS prescription in Ref. [60] already includes the low-W resonance region, adding the RES cross section of Eq. (12) double-counts part of the nuclear response. Since the attenuation calculation sums all four channels, this matching is load-bearing. The authors should state the W and Q² boundaries used for each channel and demonstrate that the four channels are disjoint and complete.","section":"Sec. III, Eq. (13) and the paragraph after Fig. 4"},{"comment":"The exclusion limits in Fig. 6 are derived from a 'likelihood-based analysis' with only a sentence describing the observed and expected event counts. To make the constraints reproducible and to assess how the attenuation-model differences propagate into the 90% C.L. boundary, the exact likelihood construction — e.g., Poisson likelihood with the 1356±43 background treated as a nuisance parameter, and the definition of the test statistic — should be given. This is especially important because the RES-induced ratio in the lower subpanels depends on the precise location of the upper boundary.","section":"Sec. V, Eq. (22) and likelihood analysis"}],"minor_comments":[{"comment":"The section is entitled 'Straight-line model' but the first sentence begins 'In the single-scattering model, any interaction...' This appears to be a copy-paste error and should be corrected to 'straight-line model'.","section":"Sec. IV A"},{"comment":"The symbol ⟨σν⟩ should presumably be ⟨σv⟩ (thermally averaged annihilation cross section). Please fix the notation.","section":"Eq. (2)"},{"comment":"The typesetting of the Gaussian prefactor is ambiguous due to the line break. As written, it appears to be 2/√(2π)σ0 mχ1, which integrates to 2, but the reader cannot tell whether a √2 is intended in the numerator. Please write the prefactor explicitly as 2/[√(2π)σ0 mχ1] (or state the intended normalization) so that the two-particle spectrum integrates to 2.","section":"Eq. (4)"},{"comment":"'is a Direct function' should read 'is a Dirac delta function' or simply 'delta function'.","section":"After Eq. (4)"},{"comment":"The caption contains a grammatically incomplete sentence: 'Together, The lower subpanels...' This should be rephrased.","section":"Fig. 6 caption"}],"recommendation":"major_revision","confidential_remarks":"The Gaussian normalization concern raised by the earlier reader appears to be a rendering artifact: the equation as typeset likely has the prefactor 2/√(2π)σ0 mχ1, which integrates to 2. The substantive risks are the in-medium Δ width and the unspecified QES/RES/DIS matching. Both are addressable with sensitivity tests and explicit definitions. I recommend major revision rather than rejection because the core idea — that the Δ(1232) channel contributes to BDM attenuation — is plausible and well-motivated; the paper just needs to make the quantitative claim robust against nuclear-model uncertainties."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, know this: the paper adds the Δ(1232) resonance channel to the boosted-DM Earth-attenuation framework of the same group's earlier work. The channel is genuinely missing, and the RES contribution is plausibly sizable at Eχ ~ 1–2 GeV for a heavy mediator. The derivation in App. A is internally consistent and uses standard nuclear machinery (impulse approximation, NuWro spectral function, MAID2007 form factors). The central result that RES lowers the PandaX-4T upper exclusion boundary by tens of percent is reasonable and new.\n\nThe main soft spots are three. First, Eq. (4) looks over-normalized: the Gaussian for the two-χ2 spectrum integrates to 2√2 rather than 2. That inflates the absolute flux by about 40%. The effect cancels in the RES-to-no-RES ratio, but it shifts the reported constraints. The authors should correct it. Second, the in-medium Δ dynamics are simply ignored: the calculation uses a free Breit–Wigner with ΓΔ = 117 MeV. From neutrino-nucleus work we expect collisional broadening, Pauli blocking, and absorption to modify the Δ self-energy at the tens-of-percent level. Since the claimed RES-induced boundary shift is also tens of percent, the exact size of the shift is not yet robust. This needs at least a discussion, ideally an estimate, before the number is trusted. Third, the paper doesn't say how the RES region is separated from the QES and DIS channels taken from ref [60]. If the DIS prescription already includes low-W resonance contributions, there's a double-counting risk. The matching should be spelled out.\n\nThere are minor issues too: the opening paragraph of Sec. IV.A mislabels the model as single-scattering in the straight-line section, and no code or data tables are provided.\n\nNone of this breaks the main point. The RES channel belongs in the attenuation framework, and the qualitative conclusion is almost certainly right. But the quantitative boundary shift carries unquantified nuclear-model uncertainty, and the flux normalization bug should be fixed.\n\nI'd send this to peer review. It's a serious, honest extension of an established program, and the referees can demand the fixes.","headline":"A genuine but incremental addition: the Δ(1232) resonance channel improves boosted-DM attenuation, yet the quoted boundary shift needs in-medium Δ treatment and a flux-normalization fix.","tokens_in":13834,"tokens_out":4817,"would_cite":true,"duration_ms":46981,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Resonant excitation of bound nucleons into the Δ(1232) baryon resonance gives a non-negligible contribution to Earth attenuation of GeV-scale boosted dark matter, and including it lowers the upper boundary of the PandaX-4T exclusion region.","keywords":["boosted dark matter","Earth attenuation","resonant scattering","Delta(1232)","dark photon mediator","PandaX-4T","nucleon spectral function","inelastic scattering"],"falsifier":"Compute the same resonant energy-loss cross section with an in-medium Δ self-energy and Pauli blocking of Δ→πN decay; if those corrections change the energy-transfer-weighted resonant cross section at Eχ ≈ 2 GeV by more than roughly the claimed shift of the exclusion boundary (tens of percent), the numerical conclusion would not survive. A cheaper cross-check is to compare the impulse-approximation response to electron- or neutrino-scattering data on iron and oxygen in the resonance region, which directly measures the N→Δ transition in the nuclear medium.","tokens_in":12721,"feed_emoji":"⚛️","tokens_out":4224,"duration_ms":48517,"temperature":0.7,"pith_summary":"This paper argues that boosted dark matter passing through the Earth loses energy and gets removed from the flux not only through elastic, quasi-elastic, and deep-inelastic scattering, but also through resonant excitation of nucleons into the Δ(1232) baryon state. In the heavy-mediator regime, this resonant channel is sizable for incoming dark-matter energies around 1–2 GeV in both iron and oxygen. Adding the channel to two Earth-propagation models reduces the detector-level flux near the spectral peak and moves the attenuation-induced upper boundary of the 90% confidence exclusion on the dark-matter–nucleon cross section to lower values. The claim matters because boosted-light-dark-matter searches rely on that upper boundary to close the allowed cross-section range, so omitting the resonance channel would place the boundary too high. The effect is larger in the straight-line energy-loss model than in the single-scattering absorption model because the resonant channel transfers a substantial fraction of the incoming energy.","feed_headline":"Baryon resonance lowers boosted-dark-matter limits","feed_subtitle":"Including Δ(1232) excitation in Earth attenuation shifts the PandaX-4T exclusion upper boundary at GeV energies.","key_machinery":"The central object is the nuclear resonance response tensor, built by folding the elementary N→Δ(1232) transition tensor with a nucleon spectral function under the impulse approximation, with the resonance's finite lifetime encoded in a Breit–Wigner distribution. This tensor converts the dark-photon-mediated dark-matter current into a nuclear excitation probability, and it is the mechanism that lets the paper add resonant scattering as a distinct channel lying between quasi-elastic and deep-inelastic scattering in energy transfer.","core_discovery":"The authors compute the dark-matter–nucleus resonant scattering cross section for excitation of the Δ(1232) resonance in the impulse approximation, folding the elementary electromagnetic N→Δ transition tensor with a nucleon spectral function and replacing the on-shell delta with a Breit–Wigner distribution of width ΓΔ ≈ 117 MeV. Using a heavy-mediator benchmark (mV′ = 10 GeV, mχ = 1 MeV), they find that the energy-transfer-weighted resonant cross section is comparable to the other inelastic channels at Eχ ≈ 1–2 GeV for Fe and O. Including this channel in both the straight-line and single-scattering Earth-attenuation models suppresses the boosted-DM flux reaching the detector and lowers the u","pith_inferences":["If the resonant channel is as sizable as claimed, the same Δ-excitation mechanism should also appear in other accelerated-dark-matter scenarios, such as cosmic-ray upscattered or atmospheric dark matter, and existing propagation calculations for those scenarios likely overestimate the flux at GeV energies.","The claimed shift of the exclusion boundary is comparable in size to plausible in-medium nuclear corrections, so the qualitative conclusion that a channel is missing is more robust than the exact numerical shift; a nuclear-model uncertainty band would clarify the strength of the constraint.","The stronger RES effect in the straight-line model suggests that a directional detector, whose trajectories traverse different Earth depths, could see a channel-dependent imprint and help discriminate attenuation models without relying on absolute flux normalization.","Because the resonance response is mediated by the electromagnetic current through kinetic mixing, precision pion-production data from neutrino or electron scattering on the same nuclear targets could calibrate the nuclear spectral function and Breit–Wigner treatment used here."],"forward_implications":["Boosted-dark-matter attenuation calculations at Eχ ≈ 1–2 GeV should include the Δ(1232) resonant channel; analyses that omit it will overestimate the flux reaching underground detectors and place the upper exclusion boundary too high.","The attenuation-induced upper boundary is lowered more in the straight-line continuous-energy-loss description than in the single-scattering absorption description, because the resonant channel transfers a large fraction of the incident energy.","The lower boundary of the exclusion region also changes: with inelastic channels included, it shifts upward by roughly a factor of two in the cases shown, except for the elastic-only straight-line case.","The relative importance of the resonant channel is specific to the heavy-mediator regime; for a light mediator, small-momentum-transfer elastic scattering dominates and the inelastic channels become less important.","The size of the effect is tied to the peak of the boosted-DM spectrum near 2 GeV, which sits at the kinematic threshold for Δ production, so the channel is most relevant for GeV-scale parent dark-matter annihilations."],"fun_headline_variants":["Resonant scattering tightens boosted-DM bounds","Δ(1232) resonance alters Earth attenuation for boosted DM","PandaX-4T limits shift with nuclear resonance inclusion","Heavy-mediator DM attenuation includes Δ(1232) channel","Boosted DM flux reduced by resonant nuclear scattering"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"That the nuclear environment can be ignored when computing the transition: the bound nucleon is described by a free-space spectral function and the produced Δ(1232) is assigned its free decay width, with no in-medium broadening, Pauli blocking, or final-state interactions included.","fun_headline_variants_meta":{"raw":{"variants":["Resonant scattering tightens boosted-DM bounds","Δ(1232) resonance alters Earth attenuation for boosted DM","PandaX-4T limits shift with nuclear resonance inclusion","Heavy-mediator DM attenuation includes Δ(1232) channel","Boosted DM flux reduced by resonant nuclear scattering"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000105,"raw_usage":{"total_tokens":840,"prompt_tokens":681,"completion_tokens":159,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":425,"completion_tokens_details":{"reasoning_tokens":86}},"tokens_in":425,"tokens_out":159,"duration_ms":2406,"temperature":1.0,"reasoning_tokens":86,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T00:10:31.112068+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the same resonant energy-loss cross section with an in-medium Δ self-energy and Pauli blocking of Δ→πN decay; if those corrections change the energy-transfer-weighted resonant cross section at Eχ ≈ 2 GeV by more than roughly the claimed shift of the exclusion boundary (tens of percent), the numerical conclusion would not survive. A cheaper cross-check is to compare the impulse-approximation response to electron- or neutrino-scattering data on iron and oxygen in the resonance region, which directly measures the N→Δ transition in the nuclear medium.","supporting_citations":[],"review_version":1}