{"id":"085a362b-edfa-4305-932f-114eabf9189c","arxiv_id":"1908.08302","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Shell-model shape factors for 36 first-forbidden beta decays shift predicted reactor antineutrino spectra by 4-5% in the 4-7 MeV range, making forbidden transitions a necessary ingredient in reactor anomaly analyses.","lead":"The paper calculates the shapes of 36 'forbidden' nuclear beta decays that dominate reactor antineutrino spectra above 4 MeV. It finds these shapes differ strongly from the simple 'allowed' approximations used in standard reactor flux models, adding several-percent corrections and large uncertainty to predictions of the reactor neutrino anomaly and the 5 MeV bump.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Combined-solution claim rests on an uncomputed 2.4%/MeV allowed-shape slope; the forbidden-only bump (~4.5%) does not by itself close the reactor anomaly or shoulder.","rationale":"The reader's weakest assumption identifies exactly the load-bearing soft spot: the quantitative size of the forbidden-induced spectral bump is measured against an allowed-shape baseline that the paper itself leaves incomplete, and the combined-solution claim depends on an uncomputed 2.4%/MeV slope. I agree with this assessment. The paper has real independent support: a sizeable shell-model calculation, explicit shape factors for 36 transitions, consistency with half-lives for most cases, and partial agreement with Fang and Brown for 140Xe. The qualitative conclusion that first-forbidden transitions are important for reactor antineutrino spectra is credible and likely correct. The conditional verdict is therefore appropriate rather than a rejection. My stress-test does not move the verdict: the central quantitative and combined-solution claims should be accepted only with the caveat that the allowed-shape baseline must be computed before the 8.5% and 'solved at the same time' statements are taken at face value.","tokens_in":44436,"tokens_out":5361,"duration_ms":60949,"concrete_test":"Replace the ad hoc 2.4%/MeV slope used in Fig. 15 with the yield-weighted average of the full allowed shape factors computed from Eqs. (22)-(24) for the dominant allowed transitions in the same ENDF+ENSDF data set, using the same shell-model Hamiltonians. If that average slope is below about 1%/MeV in the 4-7 MeV window, the forbidden-induced bump remains near 4.5% and the combined-solution claim fails; if it approaches 2.4%/MeV with a defensible uncertainty band, the claim is supported. Report the slope as a function of energy and its bin-to-bin correlations.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing quantitative assertion is that a proper treatment of forbidden transitions can simultaneously address the normalization anomaly and the spectral shoulder. That assertion is made in Secs. VII B 1 and VIII A, but it is calibrated against an allowed-shape baseline the paper itself never computes. In Sec. II D the authors list the missing induced-tensor term (Eq. 23), the full weak-magnetism expression (Eq. 22), and the Lambda-prime finite-size correction (Eq. 24), then state that these 'will be investigated in a future work.' The forbidden-only bump in Fig. 14 reaches only about 4.5% when compared to the simplified weak-magnetism slope of 0.67%/MeV (Eq. 20). The advertised ~8.5% bump in Fig. 15 is obtained by replacing that baseline with an ad hoc slope of 2.4%/MeV. If the true yield-weighted average allowed slope in the 4-7 MeV region is close to 0.67%/MeV, forbidden transitions contribute only roughly half of the observed shoulder and cannot resolve the rate anomaly; if the slope is closer to 2.4%/MeV, the combined-solution conclusion becomes plausible. The central summary sentence in Sec. IX that an increased allowed slope 'can yield a solution' is therefore a proposal resting on an uncomputed quantity, not a demonstrated result. A secondary concern, noted by the paper itself, is that individual shape factors are validated mostly against half-lives, with 86Br and 89Br not reproduced and the 136Te slope sensitive to phase convention and gA/gV; this mainly affects the absolute magnitude of the forbidden-only bump but does not remove the primary dependence on the allowed baseline.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents microscopic shell-model calculations of the dominant first-forbidden beta transitions contributing to reactor antineutrino spectra above 4 MeV, using the Behrens-Bühring formalism with full Coulomb corrections. The authors calculate shape factors for 36 transitions, find strong deviations from the allowed approximation, and use these to construct both a direct numerical correction and a Monte Carlo parametrization covering the remaining forbidden branches in summation calculations. They then compare the resulting electron and antineutrino spectra with the ILL data and with the Huber-Mueller model, concluding that a correct treatment of forbidden transitions is indispensable for both the normalization anomaly and the spectral shoulder, and proposing that a larger allowed-shape slope (2.4%/MeV) could resolve both simultaneously when combined with their forbidden corrections.","tokens_in":44816,"tokens_out":2779,"duration_ms":30175,"significance":"If the main results hold, this is a significant contribution to the reactor antineutrino problem. The work applies a standard and sophisticated formalism (Behrens-Bühring with shell-model matrix elements) to a much larger set of forbidden transitions than previous microscopic studies, and it provides a transparent uncertainty treatment by varying gA and epsilon_MEC. The comparison with the independent ILL electron data through a composite approach is a strength, as is the favorable reproduction of the 140Xe shape factor against Fang and Brown. The paper also ships a detailed parametrization procedure and Monte Carlo code logic, which is useful for future summation calculations. However, the central quantitative claim about simultaneously solving the anomaly and shoulder depends on an uncomputed baseline slope for allowed transitions, and the paper itself identifies several transitions (86Br, 89Br, 136Te) where the calculations are not fully validated. The significance of the paper is therefore real but conditional on resolving these load-bearing uncertainties.","major_comments":[{"comment":"The claim that a proper treatment of forbidden transitions can simultaneously solve the normalization anomaly and spectral shoulder rests on the 2.4%/MeV slope of the allowed shape factor, which is never computed. In Sec. II D the authors list the missing induced-tensor term (Eq. 23), the full weak-magnetism expression (Eq. 22), and the Lambda-prime finite-size correction (Eq. 24), and state that these will be investigated in future work. The forbidden-only bump in Fig. 14 reaches only about 4.5% when compared with the simplified 0.67%/MeV weak-magnetism slope, whereas Fig. 15 shows an 8.5% bump only after replacing that baseline with an ad hoc 2.4%/MeV slope. Since the true yield-weighted average allowed slope in the 4-7 MeV region is not computed, the central conclusion in Sec. IX that an increased allowed slope 'can yield a solution' is a proposal rather than a demonstrated result. This is the most load-bearing quantitative assertion of the paper and needs to be either supported by an explicit calculation of the allowed shape factors for the dominant allowed branches or softened accordingly.","section":"Sec. VII B 1, Fig. 15, and Sec. II D"},{"comment":"The reliability of the numerical shape factors is not uniform across the 36 transitions, and the paper acknowledges unresolved discrepancies. The half-lives of 86Br and 89Br are not reproduced even with the larger jj45pna model space (Sec. IV C), and the 136Te shape-factor slope is sensitive to the phase convention and to the gA/gV ratio (Sec. IV E). Since 86Br (Delta J = 1) and 89Br contribute to the highest-energy part of the spectrum and show some of the largest deviations in Fig. 3, these unresolved cases could disproportionately affect the magnitude of the bump in the 4-7 MeV region. The paper should quantify how the final cumulative spectral changes in Figs. 14 and 16 would shift if the problematic transitions were assigned conservative alternative shape factors (e.g., the allowed shape or the unique-first-forbidden shape), or should explicitly state which transitions dominate the bump and show that they are among the validated cases.","section":"Sec. IV C, IV D, and IV E; Table I"},{"comment":"The parametrized shape-factor distributions are validated against the same numerical shape factors to which the polynomial of Eq. (26) was fitted, so Fig. 10 is not an independent test of the parametrization's predictive power. The paper then applies this parametrization to roughly 1600 forbidden branches in the database, assuming the 36 selected transitions are representative, an assumption stated explicitly in Sec. VI B 1. The authors themselves note in Sec. VII B 2 that the 2-4 MeV enhancement in the parametrized results 'could be a true verifiable feature or a limitation of our parametrization' and in Sec. VIII A that the increased spread in shape factors 'is possibly a limitation of our current approach.' Given that the parametrization uncertainty is comparable to or larger than the main systematic uncertainties in the Huber-Mueller model (Appendix), the paper should provide a more direct test of representativeness, for example by using a leave-one-out cross-validation on the 36 transitions or by comparing the parametrized shape-factor distributions to the few available experimental shape-factor measurements beyond half-lives.","section":"Sec. VI B, Fig. 10, and Sec. VII B 2"}],"minor_comments":[{"comment":"The abstract uses the phrase 'ab initio electron cumulative spectra,' but the calculations are shell-model calculations with effective interactions, not ab initio; please rephrase to avoid confusion.","section":"Abstract and Sec. IV C"},{"comment":"The sentence 'This has for many isotopes resulted in a correction of branching ratios to high-lying states which had previously gone eluded due to the pandemonium effect' contains a typo; 'gone eluded' should be 'escaped detection' or 'been missed.'","section":"Sec. IV C"},{"comment":"The symbol for antineutrino momentum in Eq. (18) is written as pν, but elsewhere the paper uses p with W0 - W; please define pν and p_e explicitly in the text preceding Eq. (18) for clarity.","section":"Sec. II C 2, Eq. (18)"},{"comment":"In Table III, the row labels φ and R_IBD are not defined in the caption or in the text before the table; please define them as the integrated antineutrino flux and the inverse-beta-decay-weighted flux, respectively.","section":"Table III and Sec. VII C"},{"comment":"The caption of Fig. 11 mentions 'ENDF Q' but the legend in the figure shows 'ENDF Q' and the text refers to the Qβ approximation; please make the label consistent and explain the acronym in the caption.","section":"Sec. VII A and Fig. 11"},{"comment":"The figure caption labels the curves as '0.67%' and '2.4%' but does not state that these are slopes in units of %/MeV; please specify the units explicitly.","section":"Sec. VIII A, Fig. 18"}],"recommendation":"major_revision","confidential_remarks":"The paper is an extended version of the authors' earlier rapid communication (Ref. [18]) and includes additional transitions and an improved parametrization. The main novelty relative to that work should be described more explicitly in the introduction to avoid any perception of duplicate publication. The central claim about resolving the reactor anomaly and shoulder simultaneously is not fully supported by the present calculations, since it relies on an uncomputed allowed-shape slope; that is a fixable issue if the authors either compute the allowed shape factors or reframe the conclusion as conditional. The parametrization's lack of external validation is a more fundamental concern, but it does not invalidate the numerical shape-factor calculations themselves, which are a useful contribution to the field."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know up front. This is the most complete microscopic treatment of first-forbidden transitions in reactor antineutrino spectra so far, with 36 shell-model shape factors and full Behrens-Bühring Coulomb corrections. And the paper's biggest claim—that these forbidden transitions, combined with a steeper allowed shape slope, can resolve both the reactor normalization anomaly and the 5 MeV shoulder—is a proposal, not a demonstrated result.\n\nWhat's genuinely new: the extension from 29 to 36 transitions, the explicit shape factors for a large fraction of the flux above 4 MeV, and the Monte Carlo parametrization that turns the scatter in those shape factors into a correlated spectral uncertainty. The 140Xe comparison to Fang and Brown is favorable, which is a good sanity check. The point that forbidden transitions dominate the flux in the 4–8 MeV window is well made and is the real takeaway.\n\nThe soft spots are mostly ones the authors flag themselves. The forbidden-only bump reaches about 4.5% relative to the weak-magnetism slope of 0.67%/MeV; the advertised ~8.5% bump comes from replacing that with an ad hoc 2.4%/MeV slope. The paper lists missing induced-tensor and Lambda-prime corrections in Sec. IID and says they will be investigated later, so the combined-solution conclusion rests on a quantity that is not computed. Half-lives for 86Br and 89Br are not reproduced, the parametrization is validated only in-sample, the 2–4 MeV increase is acknowledged as possibly an artifact, and the uncertainty band is a bound treated as 1 sigma. Those are real limitations, but they do not undercut the basic finding that forbidden shapes deviate strongly from allowed approximations and affect antineutrino spectra by several percent.\n\nThe citation pattern is appropriate; this builds on earlier work by the same group and by Fang and Brown, and the formalism references are standard. The paper is written honestly, with limitations stated rather than hidden.\n\nBottom line: this is a useful paper for anyone working on reactor flux predictions or sterile neutrino fits. It deserves peer review. A referee should push on whether the allowed shape slope can be pinned down, because that controls the magnitude of the effect and the headline conclusion. I would send it out.","headline":"A solid shell-model calculation of 36 first-forbidden shape factors that shows a few-percent effect on reactor antineutrino spectra, but the headline claim that it solves both the rate anomaly and the spectral shoulder depends on an uncomputed allowed-shape slope.","tokens_in":45354,"tokens_out":1552,"would_cite":true,"duration_ms":17298,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Microscopic treatment of first-forbidden beta decays yields antineutrino spectral changes of several percent, with a bump in the 4–7 MeV range.","keywords":["first-forbidden beta decay","reactor antineutrino anomaly","5 MeV spectral bump","nuclear shell model","beta spectrum shape factor","Coulomb corrections","summation method","spectral shoulder"],"falsifier":"A high-statistics measurement of the beta spectrum of a dominant fission fragment such as 92Rb or 96Y, converted into an extracted shape factor, would settle the claim: if the measured shape factor is flat instead of the calculated downward-sloping curve, the predicted several-percent antineutrino bump shrinks or disappears.","tokens_in":44171,"feed_emoji":"⚛️","tokens_out":7961,"duration_ms":78265,"temperature":0.7,"pith_summary":"Reactor antineutrino spectrum predictions have generally treated beta decays that change nuclear spin or parity as though they had the same spectral shape as fully allowed decays. This paper computes the dominant first-forbidden transitions above 4 MeV with the nuclear shell model and a fully Coulomb-corrected beta-decay formalism, and finds that their energy-dependent shape factors deviate strongly from that approximation. The cumulative antineutrino spectrum develops a bump of up to about 4.5% between roughly 4 and 7 MeV, and up to about 8.5% when the comparison baseline has a steeper allowed slope. Since forbidden transitions carry roughly half the electron flux in the 4–8 MeV window, the paper concludes that neither the reactor flux normalization anomaly nor the so-called 5 MeV spectral bump can be understood without them.","feed_headline":"Forbidden decays reshape reactor neutrino spectra by up to 4.5%","feed_subtitle":"Shell-model shapes for 36 transitions add a bump in the 4–7 MeV range and enlarge the flux-anomaly uncertainty.","key_machinery":"The load-bearing object is the $\\beta$-decay shape factor $C(Z,W)$, the energy-dependent multiplier in the spectrum formula that contains all nuclear-structure information; in the Behrens–Bühring formalism it is built from Coulomb functions such as $\\lambda_k$, lepton phase-space factors, and nuclear form factors obtained from shell-model wave functions. The argument works by computing these shape factors fully for 36 dominant first-forbidden transitions instead of setting them to $C=1$ or to the standard weak-magnetism term, and then comparing cumulative electron and antineutrino spectra built with the calculated shapes against spectra built with those allowed approximations. The comparison is what converts a small change in electron spectra into a several-percent change in antineutrino spectra, because the steeply falling cumulative flux amplifies any reweighting of the high-energy tail.","core_discovery":"On its own terms, the paper establishes that first-forbidden transitions are not a small correction in reactor antineutrino spectra but a dominant component above 4 MeV, and that their true shapes change the predicted flux in exactly the region where experiments see an excess. For 36 transitions selected by their fission yield and branching ratio, the paper evaluates the full shape factor in the Behrens–Bühring formalism without the usual $\\xi$-approximation, including Coulomb distortion of the electron wave function and finite-size effects. The calculated shape factors for many pseudovector transitions slope downward with energy, while some pseudoscalar transitions show strong quadratic behaviour from higher-order operators, and unique transitions require the Coulomb function $\\lambda_2$ for percent-level accuracy. When these shapes are folded into summation and composite spectra matched to measured electron data, the electron spectrum changes by only about 1–2%, but the antineutrino spectrum rises by several percent in the 4–7 MeV window, with the exact magnitude depending on the allowed shape to which the comparison is made.","pith_inferences":["A direct shell-model survey of the allowed-shape slopes for the dominant fission fragments would decide whether the 2.4%/MeV baseline is physical; the paper motivates it but does not compute it.","High-statistics beta-spectrum measurements of individual dominant fragments (for example 92Rb or 96Y) could extract their shape factors and anchor or falsify the parametrization before it is applied to reactor data.","The heavy-tailed, multimodal parameter distributions suggest that older Gaussian-style uncertainty propagation for forbidden transitions underestimates the high-energy spectral error; the Monte Carlo approach used here may be the safer template for future summation calculations.","If sterile-neutrino searches continue to rely on reactor spectral predictions, the size of this forbidden-transition uncertainty will matter as much as detector systematics in the 4–7 MeV region."],"forward_implications":["Spectral-shape analyses that assume forbidden decays are allowed will mis-estimate the high-energy antineutrino flux by several percent, so the 5 MeV excess cannot be interpreted without a forbidden-transition correction.","The forbidden-transition contribution adds a correlated spectral uncertainty comparable to or larger than the leading systematic uncertainties of the standard summation model, so previous significance estimates for the rate anomaly were missing a major error term.","With the 36 numerical shape factors alone, the predicted inverse-beta-decay rate rises by about 0.8(5)%; including the parametrized shapes for all other forbidden branches raises the shift to about 2.3(13)%, which shifts but does not remove the anomaly.","If the true average slope of allowed shapes is near 2.4% per MeV rather than the standard 0.67% per MeV, the forbidden-transition bump reaches the size of the observed shoulder, offering a single mechanism for both the rate and shape discrepancies.","Because the effect is nearly identical for the main fission actinides, small differences in reactor fuel composition do not change the predicted distortion, making the correction robust for current and near-future experiments."],"supporting_citations":[{"why":"provides the baseline allowed-shape model and weak-magnetism term that the calculated forbidden shapes are compared against","marker":"[14]"},{"why":"provides the original summation-method reactor spectrum model whose treatment of forbidden transitions is corrected here","marker":"[15]"},{"why":"supplies the only previous microscopic shape-factor calculations for 136Te and 140Xe, which are compared and extended","marker":"[17]"},{"why":"gives the full beta-spectrum correction formulas used to build the spectra around the shape factors","marker":"[21]"},{"why":"provides the Behrens–Bühring formalism used to evaluate the shape factors without approximations","marker":"[23]"},{"why":"compiles the dominant beta branches above 4 MeV that determine which 36 transitions are calculated","marker":"[45]"},{"why":"supplies the measured reactor electron spectra used for normalization in the composite approach","marker":"[47]"},{"why":"established earlier that forbidden transitions can have a significant influence on the reactor anomaly, motivating the microscopic calculation","marker":"[16]"}],"fun_headline_variants":["Forbidden transitions dominate reactor antineutrino spectra","First-forbidden decays skew reactor antineutrino flux","Forbidden transitions reshape reactor anomaly spectra","Forbidden decays alter reactor antineutrino shapes","Forbidden transitions cause reactor antineutrino excess"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central quantitative claim assumes that the 'allowed' spectral shape used as the comparison is the right baseline, even though the paper itself lists missing corrections to that baseline; a different allowed slope would change both the size and the shape of the predicted bump.","fun_headline_variants_meta":{"raw":{"variants":["Forbidden transitions dominate reactor antineutrino spectra","First-forbidden decays skew reactor antineutrino flux","Forbidden transitions reshape reactor anomaly spectra","Forbidden decays alter reactor antineutrino shapes","Forbidden transitions cause reactor antineutrino excess"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000404,"raw_usage":{"total_tokens":2081,"prompt_tokens":902,"completion_tokens":1179,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":518,"completion_tokens_details":{"reasoning_tokens":1105}},"tokens_in":518,"tokens_out":1179,"duration_ms":10301,"temperature":1.0,"reasoning_tokens":1105,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:43:57.471359+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A high-statistics measurement of the beta spectrum of a dominant fission fragment such as 92Rb or 96Y, converted into an extracted shape factor, would settle the claim: if the measured shape factor is flat instead of the calculated downward-sloping curve, the predicted several-percent antineutrino bump shrinks or disappears.","supporting_citations":[{"cited_title":"The 5 MeV bump - a nuclear whodunit mystery","cited_arxiv_id":"1609.03910","evidence_quote":"provides the baseline allowed-shape model and weak-magnetism term that the calculated forbidden shapes are compared against"},{"cited_title":"Hayen, J","cited_arxiv_id":null,"evidence_quote":"provides the Behrens–Bühring formalism used to evaluate the shape factors without approximations"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"compiles the dominant beta branches above 4 MeV that determine which 36 transitions are calculated"},{"cited_title":"Herman and A","cited_arxiv_id":null,"evidence_quote":"supplies the measured reactor electron spectra used for normalization in the composite approach"}],"review_version":1}