{"id":"3f7a9e1b-7fbb-432b-b19b-1433090db177","arxiv_id":"2509.14534","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A first-principles QFT treatment of dark matter scattering off atomic electrons shows that free-electron factorization can fail and that relativistic Dirac wave functions reduce the xenon atomic factor by 30-50% relative to non-relativistic Schrödinger wave functions.","lead":"This paper argues that the standard shortcut for dark matter scattering off bound electrons, which factors the calculation into a free-electron piece times an atomic form factor, can produce negative cross sections. The authors build a quantum field theory treatment of bound and ionized electron states and show that relativistic corrections shrink the predicted xenon ionization rate by roughly 30 to 50 percent.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 30–50% suppression claim mixes relativity with a per-shell effective-charge recalibration: Dirac and Schrödinger wavefunctions are compared with different Z_eff values (e.g., 6.87 vs 6.5 for 5s), so the reduction may be partly a model-calibration artifact.","rationale":"The formal parts of the paper are valuable: the QFT derivation in Sec. 3, the non-relativistic factorization in Sec. 4, and the demonstration that the free-electron factorization can become negative are all coherent and independently checkable. The central quantitative claim, however, depends on a comparison between relativistic and non-relativistic hydrogen-like wavefunctions in which the effective charge is recalibrated separately for each model. Comparing Dirac and Schrödinger solutions with different Z_eff does not isolate relativity; it conflates relativistic dynamics with a different screening/calibration assumption. The reader's weakest assumption identified the hydrogen-like single-particle model as load-bearing; my concern is more specific: even within that model, the two sides of the comparison are not evaluated at the same potential. A same-Z_eff recomputation would settle the issue cleanly. The existing figure-text inconsistency in Fig. 4 and the absence of numerical uncertainty estimates reinforce the need for this check. The verdict is unchanged: conditional acceptance pending a same-potential comparison and a clarified figure.","tokens_in":33847,"tokens_out":10401,"duration_ms":103506,"concrete_test":"Recompute R(Tr) for Xe 1s and 5s using one common fixed effective charge for both the Dirac and Schrödinger radial integrals, e.g., Z_eff = 50.4 for 1s and Z_eff = 6.5 for 5s, with the same continuum final-state Z_eff as in Ref. [61]. If R_rel/R_nr remains 30–50% over the displayed Tr range, the suppression is genuinely relativistic; if it shrinks toward 10–20% (or near 1 for 5s), the headline claim is dominated by the Z_eff recalibration. A useful cross-check is to repeat the 1s comparison with the Dirac-Hartree-Fock/RPA potentials of Refs. [55,68], which removes the single-particle screening ambiguity entirely.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The quantitative headline is the 30–50% reduction in the K-factor phase-space ratio R(Tr) when Dirac wavefunctions replace Schrödinger wavefunctions for xenon. In the hydrogen-like implementation, the effective charge is not held fixed across the two calculations. Sec. 5.2.1 fixes Z_nκ separately for the Dirac bound state through Eq. (5.13), while the non-relativistic radial functions use a Schrödinger-calibrated Z_nl chosen to reproduce the same empirical binding energy; Sec. 5.3 explicitly quotes Z_eff = 6.87 (Dirac) vs 6.5 (Schrödinger) for 5s, with correspondingly different energies. The K-factor comparison in Sec. 5.4 inherits these different potentials. The authors themselves note (Fig. 3, right panel, and Sec. 5.4) that the suppression scales with Z_eff, so a ~6% charge mismatch for 5s—and a comparable mismatch for inner shells—can produce an apparent 30–50% reduction that is not purely relativistic. Because the abstract claim is precisely this numerical reduction, the paper does not currently establish it. The figure-text mismatch (Fig. 4 caption says n=3 while the text describes 1s and 2s) further prevents reconstruction of which Z_eff grid was used.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper argues that the standard factorization of DM-bound-electron scattering into a free-electron matrix element times an atomic form factor is inconsistent, because imposing free-electron kinematics on bound electrons can produce negative squared matrix elements (Sec. 2). It then develops a QFT formalism in the Furry picture, quantizing electrons in the atomic Coulomb potential (Sec. 3), derives the cross section, recovers the non-relativistic factorization in the appropriate limit (Sec. 4), and computes scalar-interaction atomic K-factors for xenon using Dirac and Schrödinger hydrogenic wave functions (Sec. 5). The authors report a 30%–50% reduction in the phase-space ratio R(T_r) when relativistic wave functions are used, and conclude that a relativistic treatment is necessary.","tokens_in":34091,"tokens_out":9631,"duration_ms":89783,"significance":"If the quantitative claim is established, the paper provides a cleaner theoretical basis for DM-electron scattering calculations and a practical method for evaluating relativistic atomic form factors; the formal QFT derivation, the Wigner-Eckart reduction of the K-factor, and the analytic transformation of the radial integrals are useful contributions. The paper is also explicit about the negative-cross-section pathology of the factorized treatment, which is an important caution for the field. However, the headline 30%–50% reduction currently rests on a model-dependent comparison that is not fully controlled, so the significance is somewhat tempered until the quantitative claim is isolated and validated.","major_comments":[{"comment":"The comparison supporting the 30%–50% reduction is not controlled with respect to the effective charge. The Dirac and Schrödinger bound states are computed with different Z_eff values (e.g., Z_nκ=6.87 vs Z_nl=6.5 for 5s, quoted in Sec. 5.3), and Fig. 3 (right) shows that the relativistic phase shift grows with Z_eff. If both Z values are calibrated to the same empirical binding energy, the comparison uses different Coulomb potentials and therefore conflates the relativistic treatment of the electron with the different effective charge required by the hydrogenic model. Please provide a control calculation at identical Z_eff for the two wave-function sets, or explicitly quantify how much of the reported 30%–50% reduction is attributable to the Z_eff difference rather than to Dirac vs Schrödinger dynamics.","section":"Secs. 5.3–5.4, Eq. (5.30), Fig. 4"},{"comment":"The headline reduction is extracted from the K-factor-only ratio R(T_r)=∫|q|K(ΔE,q)d|q|/(2m_e T_r), not from the full differential cross section in Eq. (5.2), which contains |M_χ^{SS}|^2 D_SS^2 inside the q integral. Since the DM spinor factor and the mediator propagator depend on q, the reduction in dσ/dT_r is not necessarily 30%–50% unless those factors are effectively q-independent over the integration range. The abstract and conclusion generalize beyond what is actually computed. Please compute the full cross-section ratio for a representative mediator mass, or state explicitly the conditions under which it reduces to the K-factor ratio.","section":"Sec. 5.4, Eq. (5.2), definition of R(T_r)"},{"comment":"The central quantitative figures are internally inconsistent. The Fig. 4 right-panel caption says n=3, l=0–3, while the text describes the panel as showing 1s and 2s sub-shells with varying Z_eff; the Fig. 2 caption quotes T_r=24 keV while the text says T_r=50 keV. Because these figures carry the paper's main quantitative claim, the captions and the text must be reconciled and the exact states and Z_eff values used in each panel must be specified unambiguously.","section":"Fig. 4 and Fig. 2"},{"comment":"The calculation assumes that each xenon electron is a single electron in a hydrogen-like Coulomb potential with a shell-dependent Z_eff, and it does not include electron-electron screening beyond the effective charge, final-state relaxation, or many-body correlations. If the 30%–50% reduction is intended as a xenon-specific prediction, the paper should justify this approximation for inner shells and compare with existing many-body or RPA calculations (e.g., Refs. [38,68]); if the result is meant as an illustration of the formalism, that limitation should be stated prominently in the abstract and conclusion.","section":"Secs. 4.2 and 5.2, hydrogenic approximation"}],"minor_comments":[{"comment":"The spinor trace is written as Tr[\\bar u(me)u(me)\\bar u(me)u(me)], while the subsequent text switches to the projection u(me)\\bar u(me) and γ^0. Please align the notation and define the spinor normalization convention explicitly.","section":"Eq. (5.29a)"},{"comment":"There is a typo: 'Comloub' should be 'Coulomb' in the heading of Sec. 5.2.","section":"Sec. 5.2 heading"},{"comment":"The phrase 'our precious work [58]' should read 'our previous work [58]'. There are also typos such as 'qunatum' near Eq. (5.30b).","section":"Sec. 5.4"},{"comment":"The color assignment for the ionized wave functions is inconsistent between the caption (T_r=24 keV) and the main text (T_r=50 keV); please unify.","section":"Fig. 2 caption"},{"comment":"No numerical tables or code are provided for the K-factors or the phase-space ratios. Given the nontrivial radial integrations, a small table of R(T_r) values or a code repository would substantially improve reproducibility.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The formal QFT derivation appears internally consistent and the paper addresses a real issue in the factorized treatment. However, the quantitative 30%–50% claim is currently supported by a comparison that mixes effective-charge calibration with the relativistic wave-function effect, and the central figures contain caption/text mismatches that should have been caught before submission. I would be comfortable with a revised version that isolates the relativistic effect, computes or bounds the full cross-section ratio, and corrects the figures."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the formal part is the real contribution. The claim that the usual free-electron factorization can produce negative cross section is correct and worth having; the Furry-picture quantization and the Wigner-Eckart reduction of the K-factor are worked out carefully. I would send this to a referee. But I would not let the abstract keep the 30–50% number until the authors isolate relativity from model recalibration.\n\nThe negative-cross-section pathology in Sec. 2 is real. They take the free-electron matrix element, impose bound-electron kinematics, and show a large region where |M|^2 goes negative. Anyone using the standard factorization should know this. The new formalism in Secs. 3–4 is also mostly sound: quantize the electron in the atomic potential, keep only the energy-conserving delta, use the ionized-state phase space, and recover the usual atomic factor in the non-relativistic limit. The analytic trick for the oscillatory radial integrals using the integral representation of 1F1 is useful and seems correct as far as I can check.\n\nThe relativistic suppression itself is not entirely new—Roberts et al. and Bloch et al. already reported reduced rates with Dirac wave functions—but the explicit demonstration of the pathology and the systematic reduction of the integrals are new elements. Credit where it is due.\n\nNow the soft spot. The headline number is computed as a ratio of K-factors, not from the full differential cross section, and the relativistic and non-relativistic wave functions are not evaluated in the same potential. For 5s they quote Z_eff = 6.87 for Dirac and 6.5 for Schrödinger, with correspondingly different binding energies. They choose effective charges to match atomic data, so the comparison mixes a real relativistic phase shift with a 5–6% change in the potential. The paper itself notes that the suppression grows with Z_eff. That means the “relativistic effect” is not cleanly separated from recalibration. A clean calculation would use identical Z_eff for both wave functions, or at least show the R(Tr) comparison at fixed Z_eff. Until that is done, 30–50% is a model-dependent estimate, not a demonstrated relativistic correction.\n\nThere is also a mismatch in Fig. 4: the caption says n=3 for the right panel while the text describes 1s and 2s. That is minor, but it makes reconstruction harder. There is no code or data, and no uncertainty estimate, which is consistent with the paper being a formalism paper, but it does mean the reader cannot independently check the numerical integrals.\n\nWho is this for? Direct detection theorists, especially those who use QEdark/DarkARC-type atomic factors for sub-GeV DM. They should read Sec. 2 and the formalism; they should be careful with the suppression estimate. I would accept for peer review, but insist on the fixed-Z_eff comparison and the full differential cross section before publication.","headline":"The QFT reframing is worth taking seriously, but the headline 30–50% relativistic suppression is not established because the two calculations use different effective charges.","tokens_in":34653,"tokens_out":2935,"would_cite":true,"duration_ms":29244,"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":"Treating xenon electrons as relativistic bound states lowers dark-matter scattering rates by 30–50%.","keywords":["dark matter electron scattering","atomic form factor","bound electrons","relativistic effects","Dirac equation","xenon direct detection","sub-GeV dark matter","K-factor"],"falsifier":"Compute the inner-shell K-factor for xenon using self-consistent many-body Dirac wave functions that include screening and ionic relaxation; if the 30–50% reduction relative to the Schrödinger result shrinks or vanishes, the claimed relativistic effect is an artifact of the single-particle hydrogen-like model rather than a robust prediction.","tokens_in":33615,"feed_emoji":"⚛️","tokens_out":13699,"duration_ms":112562,"temperature":0.7,"pith_summary":"This paper argues that the standard factorization of dark-matter–electron scattering—a free-electron matrix element times an atomic form factor—is not self-consistent, because it imposes free-electron kinematics and phase space on electrons that are bound and off-shell, and in part of parameter space it even yields negative differential cross sections. Starting from quantum field theory, the authors quantize the electron directly in the atomic Coulomb potential, so the initial and final electron states are exact bound and ionized solutions rather than plane waves, and derive a cross section whose atomic content lives in a relativistic K-factor. They evaluate that K-factor for scalar dark-matter interactions in xenon with both Schrödinger and Dirac wave functions. The relativistic treatment reduces the phase-space ratio and differential cross section by 30–50%, an effect that grows with the effective nuclear charge and comes mainly from the amplitude and Coulomb phase of the ionized final-state wave function. A reliable sub-GeV dark-matter search therefore needs both a consistent bound-state formalism and relativistic atomic wave functions.","feed_headline":"Relativistic bound electrons cut dark matter rates by 30–50%","feed_subtitle":"Dirac wave functions for atomic electrons fix a cross-section formula that can otherwise go negative.","key_machinery":"The object that carries the calculation is the atomic K-factor $K^S_{n\\kappa}(q,\\Delta E)$, the sum over final ionized states of the squared inner product of the initial bound-state and final ionized-state electron wave functions with the momentum-transfer plane wave. In the non-relativistic limit this reduces to a sum over overlaps of Schrödinger wave functions, while the relativistic version uses four-component Dirac spinors and is simplified with the Wigner–Eckart theorem into a double sum over $\\kappa'$ and $L$, with radial integrals $R_{PP}$ and $R_{QQ}$ over the large and small components. The companion quantity is the phase-space ratio $R(T_r) = \\int |q|\\, K(q,\\Delta E)\\,d|q| / (2m_e T_r)$, which measures how far the atomic response is from a free electron at rest and is the quantity that drops by 30–50%. The numerical bottleneck, highly oscillatory radial integrals coming from the confluent hypergeometric function in the ionized wave function, is removed by switching to the integral representation of that function, which turns each radial integral into a compact sum of Gauss hypergeometric functions.","core_discovery":"The central claim is that $|M|^2$ for dark-matter scattering off a bound electron does not factorize as $|M_{\\rm free}|^2$ times an atomic form factor. Because the initial electron has negative binding energy and a momentum distribution, the free-electron on-shell dispersion and phase space do not apply; imposing them can make $|M|^2$ negative, for example for an N-shell xenon electron with $m_\\chi = 10$ keV, $T_\\chi = 10$ keV, and recoil energies between 1.5 and 5 keV. A consistent quantum field theory treatment, with electron fields quantized in the bound-state picture and final ionized states normalized by recoil energy, gives a differential cross section integrated over momentum transfer, with the atomic effects encoded in a relativistic K-factor built from Dirac wave functions. For scalar interactions in xenon, replacing the non-relativistic Schrödinger wave functions by Dirac wave functions suppresses the phase-space ratio $R(T_r)$ by 30–50% for the $s$-shell states; the suppression grows with the effective nuclear charge and is driven mainly by the reduced amplitude and altered Coulomb phase of the final ionized wave function, not by a change in the bound-state wave function.","pith_inferences":["Beyond the paper, the growth of the suppression with effective nuclear charge suggests that even heavier targets such as gold or lead would show a larger relativistic reduction in inner-shell atomic factors; a cross-element scan with the same formalism would test this scaling.","Beyond the paper, only scalar interactions are demonstrated; for vector or axial-vector couplings the electron bilinear does not reduce to a single scalar overlap, so the relativistic correction could differ in size or sign, and applying the same Wigner–Eckart machinery to those couplings is a natural next step.","Beyond the paper, a self-consistent many-body treatment of xenon could shift the quoted 30–50% number, since screening and ionic relaxation affect exactly the inner-shell wave functions where the relativistic effect is largest.","Beyond the paper, the same phase-shift mechanism should be visible in neutrino–electron scattering in liquid xenon, where the predicted recoil spectrum near threshold would fall below the Schrödinger-based prediction by a measurable factor."],"forward_implications":["Existing XENON-style exclusion limits on sub-GeV dark matter that use non-relativistic Schrödinger form factors would loosen by roughly the same 30–50%, because the true bound-electron cross section is smaller than assumed.","The unphysical negative differential cross sections that appear for fast incoming dark matter, such as cosmic-ray-boosted dark matter, disappear in the consistent formalism, so the full parameter space becomes usable.","Because the suppression grows with effective nuclear charge, inner-shell ionization in high-$Z$ targets like xenon is affected most, while outer shells and lighter atoms remain closer to the non-relativistic prediction.","The phase-space ratio returns to unity at high recoil energy, so relativistic atomic corrections matter mainly in the low-energy window relevant to sub-GeV dark matter rather than for very energetic electron recoils.","The same bound-state formalism applies to other neutral projectiles that ionize atomic electrons, notably neutrinos, so predicted solar-neutrino and supernova-neutrino electron-recoil rates in xenon inherit the same relativistic correction."],"supporting_citations":[{"why":"Introduces the factorization of dark-matter–bound-electron scattering into a free-electron matrix element and an atomic form factor; this is the scheme the paper argues is inconsistent.","marker":"[22]"},{"why":"Applies that factorization to derive XENON10 limits, establishing the non-relativistic phenomenology whose rates the relativistic K-factor would lower by 30–50%.","marker":"[48]"},{"why":"Extends the atomic form-factor calculation to semiconductor and xenon targets and supplies the Schrödinger wave-function ingredients the paper compares against.","marker":"[50]"},{"why":"Previous work by two of the authors defining the atomic K-factor and phase-space treatment that the relativistic formalism generalizes.","marker":"[58]"},{"why":"Earlier accurate relativistic calculations of dark-matter–electron scattering that set the context for the suppression found here.","marker":"[53]"},{"why":"Uses the full Dirac equation for electron recoils in direct detection, one of the results the paper builds on for the relativistic reduction.","marker":"[55]"},{"why":"Provides the general atomic-response formalism and the convention that the ionized final-state electron keeps the effective charge of its initial shell, used in the radial wave functions.","marker":"[61]"},{"why":"A quantum field theory treatment of low-energy neutrino scattering by bound electrons, whose Furry-picture quantization of bound electron states this paper builds on.","marker":"[89]"},{"why":"The quantum field theory textbook whose scattering-theory conventions underlie the wave-packet, phase-space, and cross-section derivations.","marker":"[91]"},{"why":"The special-functions source for the integral representation of the confluent hypergeometric function that makes the radial integrals numerically tractable.","marker":"[107]"}],"fun_headline_variants":["Atomic bound electrons break dark matter scattering assumptions","Old dark matter formula goes negative; Dirac equation fixes it","Dirac wave functions cut dark matter rates by up to 50%"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes each xenon electron is a single particle in a hydrogen-like Coulomb potential with a shell-dependent effective charge, ignoring electron–electron screening, correlations, and the rearrangement of the ion left behind; if those many-body effects change the inner-shell wave functions, the size of the claimed relativistic reduction could be different.","fun_headline_variants_meta":{"raw":{"variants":["Atomic bound electrons break dark matter scattering assumptions","Old dark matter formula goes negative; Dirac equation fixes it","Dirac wave functions cut dark matter rates by up to 50%"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001623,"raw_usage":{"total_tokens":6453,"prompt_tokens":937,"completion_tokens":5516,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":553,"completion_tokens_details":{"reasoning_tokens":5463}},"tokens_in":553,"tokens_out":5516,"duration_ms":35367,"temperature":1.0,"reasoning_tokens":5463,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T15:50:37.090289+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the inner-shell K-factor for xenon using self-consistent many-body Dirac wave functions that include screening and ionic relaxation; if the 30–50% reduction relative to the Schrödinger result shrinks or vanishes, the claimed relativistic effect is an artifact of the single-particle hydrogen-like model rather than a robust prediction.","supporting_citations":[{"cited_title":"Solar Active-Sterile Neutrino Conversion with Atomic Effects at Dark Matter Direct Detection Experiments","cited_arxiv_id":"2112.05560","evidence_quote":"Previous work by two of the authors defining the atomic K-factor and phase-space treatment that the relativistic formalism generalizes."},{"cited_title":"Scattering of low energy neutrinos and antineutrinos by neon and argon","cited_arxiv_id":"2209.01732","evidence_quote":"A quantum field theory treatment of low-energy neutrino scattering by bound electrons, whose Furry-picture quantization of bound electron states this paper builds on."},{"cited_title":"Andrews, Richard Askey, and Ranjan Roy, Special Functions","cited_arxiv_id":null,"evidence_quote":"The special-functions source for the integral representation of the confluent hypergeometric function that makes the radial integrals numerically tractable."}],"review_version":2}