{"id":"9fe0d807-42b2-415c-8b1f-cd6cb123a047","arxiv_id":"2512.20560","paper_version":3,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A non-factorized amplitude treatment with a fitted sign-changing nuclear transition density reduces the predicted νe-71Ga capture rate by ~20%, absorbing the gallium anomaly without new physics.","lead":"This paper argues that the long-standing 'gallium anomaly' — a 20% shortfall in neutrino captures measured by gallium experiments — may come from an error in the standard calculation of the capture rate, not from new physics. The authors show that allowing the nuclear and electron wave functions to mix inside the nucleus, with a plausible sign-changing nuclear response, lowers the predicted rate enough to explain the discrepancy.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Anomaly resolution rests on a fitted nodal transition density whose shape (node position/depth) is neither independently derived nor robust; the paper proves existence, not that the physical density has this node.","rationale":"The reader's weakest_assumption identifies precisely the load-bearing concern: the physical transition density must have a node with fitted parameters, but this is not independently established. I agree. The paper's methodological contribution—relaxing leptonic factorization and using the full transition density—is physically sound, and the half-life constraint plus the 2pF benchmark are legitimate internal consistency checks. However, the central quantitative result is underdetermined: the fit includes the anomaly cross sections as input data, and the nodal shape is not predicted by nuclear theory. The paper explicitly frames the result as a 'proof of principle' and concedes non-uniqueness, so the appropriate verdict is CONDITIONAL rather than rejection. The proposed ab initio calculation or sensitivity scan would settle whether the node is physical. No internal mathematical inconsistency or author misconduct is alleged; the concern is about the evidentiary weight of a fit against a small dataset. Therefore the reader's CONDITIONAL verdict remains unchanged.","tokens_in":11250,"tokens_out":13311,"duration_ms":161260,"concrete_test":"Independently compute the 71Ga→71Ge Gamow-Teller transition density ρ_TD(r) using a large-scale shell model (e.g., GXPF1A or PFSDG-U) or an ab initio method with chiral NN+3N forces as in Refs. [67,72], then evaluate Eqs. (8)-(11) to predict σ_gs for 51Cr and 37Ar. If the predicted cross sections are not within the reduced ~4.4/5.1×10^-45 cm² range (or the density is positive-definite/no node), the proposed resolution is falsified. A practical fallback: re-run the combined fit with the DG node position fixed at r_b ± 0.5 fm; if σ_gs shifts by more than ~10%, the mechanism is fine-tuned rather than robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim depends on the physical 71Ga→71Ge Gamow-Teller transition density having a sign-changing radial profile with the fitted parameters (Table I, SOG/DG/DG2). This node is load-bearing: without it the SG fit fails and the 2pF benchmark fails. But the nodal parameters are free parameters optimized against the two anomaly cross sections themselves (Eq. 13); with 5-6 parameters and only three constraints (two σ_gs and t1/2), good χ² is essentially guaranteed and the ~20% reduction can be a fitting artifact. The paper concedes the densities are 'not unique' and that a first-principles calculation remains challenging [67]. Moreover, because the electron and neutrino wavefunctions are nearly constant over the nuclear volume, a 20% suppression requires delicate cancellation between positive and negative lobes, making the result highly sensitive to the exact node position, width, and depth; no fit uncertainties or sensitivity analysis are reported. Reproducing the 71Ge half-life is a genuine but weak constraint—it fixes a weighted integral of ρ_TD, not its radial shape.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper revisits the theoretical treatment of the ν_e + 71Ga → 71Ge + e^- cross section and argues that the standard detailed-balance approach, combined with the factorization of leptonic and nuclear wave functions, is biased. It introduces a non-factorized amplitude formalism in which the transition is expressed through a weak transition density ρ_TD(r) (Eqs. 7–11), evaluated with exact Dirac-Hartree-Fock-Slater lepton wave functions. Since ρ_TD is not known from first principles, the authors fit several phenomenological parameterizations (SG, SOG, DG, DG2) simultaneously to the measured 71Ge half-life and to the GALLEX/SAGE/BEST experimental ground-state cross sections (Eqs. 13–14). With sign-changing SOG/DG/DG2 densities they reproduce both quantities and obtain σ_gs ≈ 4.42–4.45 ×10^-45 cm^2 (51Cr) and 5.13–5.17 ×10^-45 cm^2 (37Ar), a ~20% reduction relative to previous estimates, and conclude that the gallium anomaly can be resolved without invoking sterile neutrinos. The manuscript explicitly acknowledges that the adopted densities are not unique and that this is a proof of principle.","tokens_in":11439,"tokens_out":6320,"duration_ms":69430,"significance":"If the required transition density were independently established, this would be a significant resolution of a long-standing 5σ anomaly and would weaken the sterile-neutrino interpretation of the gallium data. The non-factorized amplitude treatment, the use of exact lepton wave functions, and the imposition of the precisely measured 71Ge half-life as a constraint are genuine improvements, and the paper is transparent about its model dependence. However, in its present form the quantitative result is an existence proof: the reduction is obtained by fitting the same experimental cross sections that define the anomaly, with free nodal-density parameters, and no nuclear-structure calculation or independent observable fixes the node. The paper's value therefore lies mainly in identifying an accurate weak transition density as a concrete path forward, rather than in demonstrating that the path is realized in 71Ga.","major_comments":[{"comment":"The claimed solution is fitted rather than predicted. The χ²_IBD defined in Eq. (13) is minimized against the very same experimental σ_gs values that define the gallium anomaly. The SOG/DG/DG2 entries in Table I achieve σ_gs close to σ_exp because the parameters are optimized to do exactly that; with five or six free parameters and only two cross-section data points, a near-perfect χ² is essentially guaranteed. The 71Ge half-life is an independent constraint, but it is a single weighted integral of ρ_TD and is reproduced by all parameterizations, including the failing SG and 2pF benchmarks. Thus Table I demonstrates the existence of densities that remove the anomaly, not that the physical density has this property. The abstract's wording 'we find that the revised cross-section can be significantly reduced' overstates the predictive content; at present the reduction is an imposed outcome","section":"Fit results, Eq. (13), Table I"},{"comment":"The resolution is entirely tied to the sign-changing (nodal) form of ρ_TD. The single-lobe SG and the 2pF charge-density benchmark fail to solve the anomaly, while SOG/DG/DG2 succeed precisely because they allow a node. No independent evidence is presented that the 71Ga→71Ge Gamow-Teller transition density actually has a node with the fitted parameters. The cited ab initio studies [67,72,73] concern weak densities in general or in other transitions and do not provide this specific density; the text itself concedes that 'a first-principles calculation with controlled uncertainties remains challenging [67]'. The half-life constraint cannot determine the node position, width, or depth. Because the ~20% suppression relies on cancellation between positive and negative lobes, small changes in the node parameters could substantially alter or eliminate the effect. This is a load-bearing assumpti","section":"Transition density parametrization, Eq. (12), Fig. 1"},{"comment":"No uncertainties are reported for the best-fit transition-density parameters or for the resulting cross-section reduction. Given that the suppression arises from a delicate cancellation between positive and negative contributions to the integrals in Eqs. (8), (9), and (11), the stability of the result under parameter variations should be demonstrated. A covariance matrix, a profile-likelihood scan over σ_gs, or a scan over node position and amplitude is required to assess whether the ~20% reduction is robust or an artifact of the chosen functional forms. This is particularly important because the fit is underdetermined relative to its parameters, as noted above.","section":"Table I, Conclusions"},{"comment":"The title and abstract attribute the solution to moving beyond the leptonic wave-function factorization. However, the numerical comparison in Table I shows that the non-factorized formalism alone, when combined with a conventional positive-definite density (SG or 2pF), does not solve the anomaly; the suppression appears only when an ad hoc sign-changing ρ_TD is inserted. Thus the essential ingredient is the assumed nodal density, not the removal of the factorization. The framing should be corrected to avoid implying that the improved lepton treatment itself produces the ~20% reduction. This distinction matters for how readers interpret the mechanism.","section":"Abstract and Beyond the factorization scheme"}],"minor_comments":[{"comment":"In the SOG row, the entries '0.0050.0003' appear to be concatenated without a separator. The table caption should define χ²_IBD and χ²_EC explicitly, since the subscripts may be confused with the source labels.","section":"Table I"},{"comment":"The curves are hard to identify because SG and 2pF are scaled by a factor of 10 and the line styles may be difficult to distinguish in print. Adding a legend or explicit labels, including the unscaled curves, would improve readability.","section":"Fig. 1"},{"comment":"The notation p_j^e E_j^e in the sum over electron energies is confusing. Define p_e and E_e before use and clarify that the sum runs over the electron-energy branches of the source.","section":"Eq. (1)"},{"comment":"There is a typo: 'ab-initiostudies' should read 'ab initio studies'. Also, in the same paragraph, 'Aρ_0(c,a)' should be defined explicitly so the normalization of the 2pF density is unambiguous.","section":"Transition density parametrization, text after Eq. (15)"}],"recommendation":"major_revision","confidential_remarks":"The formal machinery and the critique of detailed-balance factorization are sound, and the paper is honestly labeled as a proof of principle. The main issue is that the quantitative resolution is an artifact of fitting the anomaly data with ad hoc nodal densities; no independent input fixes the node, and no sensitivity analysis is provided. I would be willing to reconsider after the authors either (i) obtain or cite a credible nuclear-structure calculation of ρ_TD for this transition, (ii) demonstrate robustness under parameter variations and report uncertainties, or (iii) significantly reframe the claim as a proposal for the data needed to settle the anomaly rather than a solution. The paper may also benefit from a clearer separation of the effect of the improved lepton wave functions from the effect of the assumed density shape."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nYou should read this if you care about the gallium anomaly. The paper does two things worth noticing. First, it drops the usual factorization of leptonic and nuclear parts in the 71Ga(νe,e−)71Ge cross-section and treats the full amplitude with a radial Gamow-Teller transition density. That critique of detailed balance is correct: the standard treatment assumes the electron and neutrino wave functions are constant over the nucleus and maps the IBD matrix element onto EC, which is only an approximation. The formalism (Eqs. 7–11) is adapted cleanly from beta-decay literature and is applied correctly. Second, the 71Ge half-life constraint is real and has teeth; a no-node single Gaussian fails to solve the anomaly even while matching t1/2 via nuisance pulls, and the 2pF charge-density benchmark fails too. Good control checks.\n\nWhere it gets soft: the resolution is not a prediction. The SOG, DG and DG2 parameters are fit by minimizing χ² against the experimental cross-sections (Eq. 13). With five or six shape parameters and essentially two cross-section data points plus the half-life, a good fit is almost guaranteed. The suppression is \"intimately connected\" to a node in the transition density, but that node is a fit ingredient, not something independently derived. The paper is honest about this—it calls the result a proof of principle and says the densities are not unique—but it means the ~20% reduction is a fitted outcome. No fit uncertainties are reported, no sensitivity scan over node position or depth, and no check of what a 20% reduction does to the gallium solar-neutrino rate or BEST's inner/outer counting ratio. Those are addressable, and until they are done the anomaly cannot be said to be resolved.\n\nWho gets value: neutrino phenomenologists working on the gallium anomaly, and nuclear physicists interested in weak transition densities. The paper deserves a serious referee. The methodological point is solid and the proof-of-principle claim is honestly scoped. A referee should demand an independent (even schematic) transition density from nuclear theory, error bars on the fit, and the solar/BEST consistency checks before accepting the conclusion as more than suggestive.\n\nMy recommendation: send it to peer review. It may not close the anomaly, but the framework is legitimate and the question it raises—whether detailed-balance factorization biases the cross-section—is worth the community's attention.","headline":"A legitimate challenge to the detailed-balance factorization that gets concrete numbers only by fitting a nodal transition density to the anomaly itself—proof of principle, not proof of the node.","tokens_in":12138,"tokens_out":2450,"would_cite":true,"duration_ms":25425,"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":"The gallium anomaly disappears once the 71Ga→71Ge neutrino-capture cross-section is computed with a sign-changing Gamow-Teller transition density instead of the factorized detailed-balance approximation.","keywords":["gallium anomaly","71Ga neutrino capture","inverse beta decay","detailed balance","Gamow-Teller transition density","sterile neutrino","lepton wave functions","electron capture"],"falsifier":"Compute the 71Ga→71Ge Gamow-Teller transition density from first principles with controlled uncertainties. If the radial density is positive-definite (no node), the ~20% cross-section suppression disappears and the gallium anomaly returns at full strength; the authors themselves note that such a calculation is not yet available.","tokens_in":10954,"feed_emoji":"⚛️","tokens_out":7357,"duration_ms":75073,"temperature":0.7,"pith_summary":"For thirty years, gallium-based neutrino-source experiments have seen about 20% fewer neutrino captures than predicted, a deficit now beyond 5σ. This paper argues that the deficit is a calculation artifact: the standard cross-section splits the neutrino-capture amplitude into a nuclear matrix element times a leptonic factor, using detailed balance with electron capture. That factorization is only valid if the lepton wave functions are constant over the nucleus; dropping it couples the lepton wave functions to the nuclear Gamow-Teller transition density. Using realistic electron wave functions and data-driven transition densities that reproduce the precisely measured 71Ge half-life, the authors find the ground-state cross-section drops by about 20%, into agreement with GALLEX, SAGE, and BEST. If right, the gallium anomaly is explained by nuclear and leptonic structure, and no sterile neutrinos are needed.","feed_headline":"Gallium neutrino deficit traced to a sign-changing nuclear density","feed_subtitle":"Computing the 71Ga capture with an un-factorized Gamow-Teller density matches source experiments, so no sterile neutrinos are needed.","key_machinery":"The key object is the Gamow-Teller weak transition density ρ_TD(r) = Ψ*(71Ge) Ĥ_GT Ψ(71Ga), the radial nuclear overlap weighted by the charge-exchange operator. The paper replaces the factorized amplitude ψ_e*(r0)ψ_ν(r0) M_nuc with integrals of ρ_TD(r) against exact Dirac-Hartree-Fock-Slater electron radial components g_κ, f_κ and the spherical Bessel functions j_0(qr), j_1(qr). The same density must reproduce the 71Ge electron-capture rate, which provides the tight experimental anchor: any density that solves the gallium anomaly must also give the measured half-life. The load-bearing feature is the node, which produces partial cancellation between positive and negative radial lobes in the i","core_discovery":"The central claim is that the ~20% gallium deficit is not new physics but the failure of the detailed-balance/factorization approximation. When the leptonic radial functions are kept inside the nuclear integral, the inverse-beta-decay amplitude depends on the weak transition density ρ_TD(r) rather than on a single matrix element times ψ(r0). The authors introduce phenomenological densities—single Gaussian, sum of Gaussians, and double Gaussians—fit simultaneously to the precisely measured 71Ge electron-capture half-life (11.465 d) and to the experimental 51Cr and 37Ar cross-sections. The best-fit sign-changing (nodal) densities reduce the predicted ground-state cross-section from roughly 5.3","pith_inferences":["A testable extension follows from the mechanism: because the node suppresses the cross-section through a momentum-dependent cancellation, the correction is source-energy dependent, whereas a sterile-neutrino deficit would be flat after phase-space normalization; a radioactive-source campaign at two well-separated neutrino energies could discriminate the two.","The result implicitly calls for recomputing other gallium-based rates, notably the solar-neutrino absorption cross-section on 71Ga, in the same un-factorized scheme; a ~20% reduction there would feed directly into solar-model comparisons, though the paper does not address this.","The mechanism's sensitivity to the node makes a specific prediction: high-resolution charge-exchange measurements of the 71Ga→71Ge Gamow-Teller response should show a sign-changing radial form or fragmented strength; future data fixing the node's position and depth would test the fitted densities against the same two χ² constraints.","If the anomaly is indeed nuclear in origin, the sterile-neutrino interpretation of the other short-baseline anomalies is weakened only insofar as those calculations share the same detailed-balance input; the paper does not extend its conclusion to reactor or accelerator anomalies."],"forward_implications":["If the reduced cross-sections are correct, the ~5σ gallium anomaly becomes statistically consistent, removing the strongest short-baseline hint for sterile neutrinos and aligning gallium data with recent accelerator and tritium-based bounds.","The same de-factorized treatment applies to other low-energy charged-current neutrino capture processes, so cross-sections calculated under the detailed-balance approximation may carry similar biases.","The fitted nodal densities are phenomenologically viable because they reproduce the precisely measured 71Ge half-life to within a fraction of a day, so the solution is not bought at the expense of a measured decay rate.","The paper's results are grounded in the ground-state contribution; with the adopted excited-state subtraction (5.3% for 51Cr and 5.8% for 37Ar), the source-averaged experimental cross-sections become consistent with the reduced theory.","The paper explicitly urges a dedicated nuclear-structure effort to compute or measure ρ_TD with controlled uncertainties, since a first-principles calculation with reliable error bars does not yet exist."],"fun_headline_variants":["Gallium anomaly resolved by nodal density, no new physics","Sign-changing density fixes 5σ gallium deficit","Beyond factorization: gallium deficit is just nuclear structure","Un-factorized wave functions erase gallium anomaly"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The physical 71Ga→71Ge Gamow-Teller transition density really has a radial sign-changing node with roughly the fitted position, width, and depth; no first-principles calculation or direct measurement establishes this, and the no-node parametrizations fail to resolve the anomaly.","fun_headline_variants_meta":{"raw":{"variants":["Gallium anomaly resolved by nodal density, no new physics","Sign-changing density fixes 5σ gallium deficit","Beyond factorization: gallium deficit is just nuclear structure","Un-factorized wave functions erase gallium anomaly"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000232,"raw_usage":{"total_tokens":1306,"prompt_tokens":707,"completion_tokens":599,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":451,"completion_tokens_details":{"reasoning_tokens":536}},"tokens_in":451,"tokens_out":599,"duration_ms":6970,"temperature":1.0,"reasoning_tokens":536,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T14:24:12.721266+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the 71Ga→71Ge Gamow-Teller transition density from first principles with controlled uncertainties. If the radial density is positive-definite (no node), the ~20% cross-section suppression disappears and the gallium anomaly returns at full strength; the authors themselves note that such a calculation is not yet available.","supporting_citations":[],"review_version":1}