{"id":"74fb64c8-98f8-47c0-8444-f03e875da1e4","arxiv_id":"1908.07911","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Shell-model calculations give two-neutrino double-beta-decay matrix elements of 0.0511 for calcium-48 and 0.0747 (or 0.0854) for zirconium-96, leading to larger predicted single-beta branching ratios than previous studies.","lead":"The paper computes the quantum-mechanical probabilities behind two kinds of radioactive decay in two rare nuclei, calcium-48 and zirconium-96, and combines them with measured half-lives to estimate how often the ordinary beta-decay path is taken. The estimated beta-decay branch is larger than earlier calculations suggested, so dedicated experiments might be able to detect it.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The branching predictions rest on one untested step: a gA extracted from the 2νββ Gamow-Teller channel is reused for the higher-rank axial operators of the 4th/6th-forbidden β branches; the paper gives only a plausibility argument, not a quantitative justification.","rationale":"The strongest part of the paper is the shell-model NME calculation itself: the 48Ca result with the full fp-space GXPF1A interaction and all 9470 intermediate 1+ states is close to the earlier value, and the extreme single-state dominance found for 96Zr agrees with the (3He,t) measurement of Thies et al. These are genuine, independent points of support. The claim that requires scrutiny is the branching prediction, which is exactly the part used for experimental motivation. The reader's weakest assumption identifies the transfer of gA from 2νββ to single β decay; I agree with that identification. I would sharpen it by noting that the dominant β branch is a unique forbidden transition whose axial operator is not the Gamow-Teller operator, so a single scalar gA is not automatically valid. The paper itself flags this only through the low-momentum-exchange sentence, without a quantitative test. The proposed multipole-dependent recalibration would settle the point. Because the branchings are nonlinear ratios of β and 2νββ rates, this is the step on which the headline numbers depend. The CONDITIONAL verdict remains appropriate: the calculation is not invalidated, but the transferability of gA must be demonstrated before the predicted branchings are used to plan dedicated experiments.","tokens_in":8159,"tokens_out":8217,"duration_ms":82237,"concrete_test":"Use the same shell-model wave functions to compute the 96Zr 5+ beta-decay shape factor with two different axial-operator normalizations: (i) the paper's single global gA=1.04, and (ii) an independent renormalization of the rank-5 unique-forbidden axial operator calibrated to measured forbidden beta-decay or charge-exchange strengths in the A≈100 mass region, while keeping the Gamow-Teller channel calibrated to the 2νββ half-life. If the resulting β branch differs from 18.4% by more than the quoted ±0.9%, the paper's consistency assumption is falsified; if it does not, the concern is resolved. A parallel test on 48Ca with gA=0.80 would show whether the effect is generic.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline branchings—7.5% for 48Ca and 18.4% for 96Zr—are controlled less by the shell-model NMEs themselves than by the effective axial coupling gA. The paper extracts gA^eff from Eq. (5) using the measured 2νββ half-life and the computed M2ν, then applies this same scalar gA to the β-decay shape factors. The only justification is the sentence after Table 1: 'This we consider to be a consistent approach since the 2νββ and β decays are low-momentum-exchange processes and thus the related axial couplings are expected to be quenched by a similar amount.' That is a plausible but untested assertion. The 2νββ operator is the allowed Gamow-Teller operator, whereas the dominant 5+ branch is a 4th-forbidden unique transition whose leading axial multipole has rank 5, and the 4+/6+ branches involve several interfering axial and vector multipoles. A single gA absorbed into the GT channel does not automatically renormalize those operators in the same way, especially inside a truncated model space. Because the branching ratio is λβ/(λβ+λ2νββ), even a 10–20% change in the β-decay partial half-life moves the 96Zr branching by several points, comparable to the quoted ±0.9%. This is an external-validity assumption rather than an internal inconsistency, but it is the load-bearing step for the experimental-motivation claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports large-space shell-model calculations of the two-neutrino double-beta-decay matrix elements M2ν for 48Ca and 96Zr, using the GXPF1A interaction in the full fp space for 48Ca and the glekpn interaction in an extended model space for 96Zr. All intermediate 1+ states are included. The computed M2ν values are 0.0511 for 48Ca and 0.0747 for 96Zr (0.0854 if the first 1+ state in 96Nb lies at 694.6 keV). The 96Zr calculation shows extreme single-state dominance, in agreement with the high-resolution charge-exchange experiment of Thies et al. Combining these matrix elements with measured 2νββ half-lives yields effective axial couplings geff_A = 0.80 ± 0.04 (48Ca) and 1.04^{+0.03}_{-0.02} (96Zr), which are then used to predict single-β branchings to the 4+, 5+, and 6+ states. The predicted total β branchings are 7.5 ± 2.8% for 48Ca and 18.4 ± 0.9% for 96Zr, substantially larger than earlier estimates, motivating dedicated experimental searches.","tokens_in":8512,"tokens_out":13257,"duration_ms":117583,"significance":"The 96Zr calculation is, as far as I am aware, the first large-space shell-model evaluation of M2ν for that nucleus, and the single-state-dominance result provides a theory-side confirmation of the experimental finding of Thies et al. The paper is transparent in its methodology: the NMEs are computed from the Hamiltonian rather than fitted to half-life data, the extraction of geff_A from Eq. (5) is a one-parameter step, and the gA dependence of the branchings is shown explicitly. The inclusion of all intermediate 1+ states in both nuclei is a clear improvement over earlier truncated calculations. However, the central experimental-motivation claim rests on an unquantified assumption about the renormalization of the axial coupling in highly forbidden β transitions, and the 96Zr NME carries an unquantified systematic uncertainty from the unmeasured 1+ spectrum. These issues do not affect the internal consistency of the shell-model matrix elements, but they do affect the reliability of the headline branching fractions.","major_comments":[{"comment":"The branchings are conditional on transferring a single geff_A extracted from the 2νββ Gamow-Teller channel to the forbidden β-decay shape factors. The 2νββ operator is an allowed Gamow-Teller operator, whereas the dominant 5+ branch is a 4th-forbidden unique transition whose leading axial multipole has rank five, and the 4+ and 6+ branches involve several interfering axial and vector multipoles. The sentence after Table 1 — 'This we consider to be a consistent approach since the 2νββ and β decays are low-momentum-exchange processes and thus the related axial couplings are expected to be quenched by a similar amount' — is a plausible physical statement but not a quantitative justification. Because the branching ratio is λβ/(λβ + λ2νββ), a 10–20% shift in the β partial half-life changes the 96Zr branching by several percentage points, comparable to the quoted ±0.9%. I request a sensitivity study in which the axial renormalization for the β operators is varied independently (for example, by rescaling the axial form factors in the range geff_A = 0.8–1.27 while keeping the 2νββ-derived value fixed) or, alternatively, a clear statement in the abstract and conclusion that the predicted branchings are model-dependent with respect to this assumption.","section":"Table 1 and the paragraph following Fig. 2"},{"comment":"The central value M2ν = 0.0747 for 96Zr is quoted without a systematic uncertainty reflecting the unmeasured 1+ spectrum in 96Nb. The paper itself shows that placing the lowest 1+ state at 694.6 keV changes M2ν to 0.0854, a 14% shift, and the 5+ branching from 18.4% to 16.7%. Since no 1+ states in 96Nb are known experimentally, Table 1 and the abstract should either adopt one scenario as the central value with the other included in the error budget, or present both values with equal prominence. As written, the precision implied by '0.0747' in Table 1 overstates what is known, and the same caveat propagates to the extracted geff_A.","section":"Table 1 and the paragraph following Fig. 2"},{"comment":"The quoted NMEs carry no estimate of model uncertainty. The 48Ca result differs by 5.5% from the earlier shell-model value of Horoi et al., and the 96Zr value changes by 14% under the alternative 1+ scenario, yet the errors quoted in Table 1 and in the branching fractions reflect only experimental half-life and Q-value uncertainties. The authors should add an explicit statement that M2ν and geff_A do not include uncertainties from the shell-model truncation, the choice of Hamiltonian, or the unknown 1+ energies, and, where possible, estimate the spread using the interaction/model-space variants already at hand.","section":"Abstract and Table 1"}],"minor_comments":[{"comment":"The uncertainty for the 96Zr total β branching is printed as '18.4 ± 0.09%' in the abstract and conclusion but as '18.4 ± 0.9%' in Section 3; the latter is consistent with the experimental half-life uncertainty and should be used throughout.","section":"Conclusion"},{"comment":"The denominator of Eq. (6) is ambiguous: the symbols E(1+_m), M_i, and the role of the excitation energy relative to the initial or final ground state should be defined explicitly, since the paper later discusses shifting the 1+ spectrum in 48Sc to the experimental 2200 keV.","section":"Eq. (6)"},{"comment":"The statement that 'the β decay might be up to 2.3 times faster than predicted by the older QRPA calculations in [7]' is not backed by an explicit comparison; please specify which half-life or branching ratio is compared and how the factor 2.3 is obtained.","section":"Section 3, comparison with Ref. [7]"},{"comment":"The gA-dependence figures are central to the argument, but the captions do not state the range of gA displayed or the meaning of the vertical lines; please add this information so the reader can judge the sensitivity at the extracted geff_A values.","section":"Figures 4–5 and 7–8"},{"comment":"Minor typographical issues: 'firs 1 + state' in the conclusion should read 'first 1+ state'.","section":"Conclusion"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is well within the scope of Physics Letters B and the shell-model calculations appear internally consistent. I found no evidence of circularity: the NMEs are not fitted to the half-life data, and the extraction of geff_A from Eq. (5) is transparent. The main issue is external validity rather than internal correctness: the headline branchings depend on transferring the 2νββ-derived axial quenching to forbidden β operators, and the 96Zr matrix element lacks an uncertainty from the unknown 1+ spectrum. These are fixable with additional sensitivity analyses and a more careful presentation of uncertainties, hence my recommendation of major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a serious shell-model paper, not a hype piece. The genuinely new thing is the first large-scale shell-model M2ν for 96Zr (0.0747, with 0.0854 under the 694.6 keV scenario) and the confirmation of extreme single-state dominance seen in the (3He,t) experiment. The 48Ca value is an update with a much larger intermediate-state sum (9470 states) and lands 5.5% below Horoi et al. The calculations are detailed enough to reproduce: NuShellX, standard interactions, cumulative NME plots, full intermediate sums. The gA extraction from Eq. (5) is transparent and the comparison with Barea et al. is a nice sanity check.\n\nThe soft spot is exactly the one the stress test flags: they extract a single gA_eff from the 2νββ Gamow-Teller matrix element and apply it unchanged to the 4th- and 6th-forbidden beta branches. The justification is plausible — low momentum exchange suggests similar quenching — but it is not quantitatively demonstrated. For the unique 5+ branch the half-life goes roughly as gA^-2, while for non-unique branches the shape factor has several interfering multipoles; there is no guarantee the same renormalization holds. Since the branchings depend on the ratio of beta to 2νββ rates, a 10–20% shift in the beta partial half-life moves the 96Zr branching by more than the quoted ±0.9%. This is an external-validity assumption, not an internal inconsistency, but it is the load-bearing step for the experimental-motivation claim.\n\nThere are also minor issues: M2ν is quoted without a model uncertainty; for 96Zr the intermediate 1+ energies are unknown and the result shifts by ~14% if the first state is at 694.6 keV; and the 48Ca branching uncertainty is dominated by the 5 keV Q-value error, which the paper itself acknowledges. None of these undercut the NME results themselves.\n\nWho should read this? People working on double beta decay matrix elements and beta-decay branchings in nuclei where single beta competes with 2νββ. It is a useful benchmark calculation and the 96Zr SSD result is a crisp confirmation of experiment. I would not use the branchings to plan an experiment without further testing the gA transfer, but the NMEs deserve a place in the literature.\n\nVerdict: worth sending to peer review. I would referee it myself. The branchings need to be flagged in the conclusions as conditional on the gA-transfer assumption, and the authors should be asked to state that clearly.","headline":"A serious shell-model paper with a genuinely new 96Zr NME and a clean SSD confirmation; the beta branchings are useful estimates but rest on an untested assumption about gA transfer.","tokens_in":9051,"tokens_out":2992,"would_cite":true,"duration_ms":29505,"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":"A consistent shell-model treatment now predicts measurable single-beta branches in 48Ca and 96Zr.","keywords":["double-beta decay","axial-vector coupling","48Ca","96Zr","shell model","nuclear matrix elements","single-beta decay branchings","single-state dominance"],"falsifier":"Measure the $5^+$ single-$\\beta$ branch of $^{96}$Zr: the paper predicts roughly $18.4\\%$ (or about $16.7\\%$ if the first $1^+$ state of $^{96}$Nb is at 694.6 keV), so a measured branch clearly outside this range would falsify the claim that the double-$\\beta$-derived coupling transfers to the single-$\\beta$ channel.","tokens_in":7938,"feed_emoji":"⚛️","tokens_out":13938,"duration_ms":108182,"temperature":0.7,"pith_summary":"The paper tries to show that one shell-model framework, with the same effective axial-vector coupling fixed by measured two-neutrino double-$\\beta$ half-lives, can describe both the double-$\\beta$ decay and the single-$\\beta$ decay branches of $^{48}$Ca and $^{96}$Zr. It reports $M_{2\\nu}=0.0511$ for $^{48}$Ca and the first large-scale shell-model $M_{2\\nu}=0.0747$ for $^{96}$Zr, with the latter showing extreme single-state dominance. These matrix elements, combined with measured half-lives, give effective axial couplings of about 0.80 and 1.04, which in turn predict total single-$\\beta$ branches of about 7.5% and 18.4%. The result matters because these are the only two known nuclei where single-$\\beta$ and double-$\\beta$ decay compete, and branches this large could be seen by dedicated experiments.","feed_headline":"One shell model now predicts measurable beta branches in 48Ca and 96Zr","feed_subtitle":"New matrix elements put single-beta branches at 7.5% for 48Ca and 18.4% for 96Zr, large enough for dedicated detectors.","key_machinery":"The load-bearing object is the interacting nuclear shell model expanded to large valence spaces, with all intermediate $1^+$ states included in the $2\\nu\\beta\\beta$ matrix-element sum (9470 states for $^{48}$Sc, 5894 for $^{96}$Nb). The $^{48}$Ca space is the full $fp$ shell; the $^{96}$Zr space adds the proton $f_{5/2}$, $p_{3/2}$, $p_{1/2}$, $g_{9/2}$ and neutron $g_{7/2}$, $d_{5/2}$, $d_{3/2}$, $s_{1/2}$ orbitals. The central identity is the matrix-element formula that sums Gamow-Teller transitions through intermediate $1^+$ states with energy denominators; inverting the half-life relation fixes the effective $g_A$, which is then reused in the single-$\\beta$ shape factors. The 'single-state dominance' pattern, where one first $1^+$ state supplies essentially the entire $M_{2\\nu}$, is what makes the $^{96}$Zr result particularly clean.","core_discovery":"The paper's central claim is that a consistent large-space shell-model treatment produces the $2\\nu\\beta\\beta$ matrix element and the competing single-$\\beta$ branches of $^{48}$Ca and $^{96}$Zr using one effective axial-vector coupling. For $^{48}$Ca the calculated matrix element is $M_{2\\nu}=0.0511$, 5.5% smaller than the earlier value of 0.0539; for $^{96}$Zr it gives $M_{2\\nu}=0.0747$, with essentially the whole value coming from a single intermediate $1^+$ state, confirming the extreme single-state dominance seen in charge-exchange data. If the first $1^+$ state of $^{96}$Nb lies at 694.6 keV, the $^{96}$Zr matrix element rises to 0.0854. Extracting the axial coupling from the measured half-lives gives $g_A=0.80\\pm0.04$ for $^{48}$Ca and $g_A=1.04^{+0.03}_{-0.02}$ for $^{96}$Zr, and using these same couplings to evaluate the $\\beta$ shape factors yields total single-$\\beta$ branchings of $(7.5\\pm2.8)\\%$ and $(18.4\\pm0.09)\\%$, both larger than earlier predictions.","pith_inferences":["An independent test: a dedicated experiment measuring the $5^+$ single-beta branch of $^{96}$Zr can invert the measured half-life to extract its own $g_A$; agreement with the double-beta-derived value would confirm the consistency assumption, and disagreement would break it.","Because the $^{48}$Ca branching uncertainty is dominated by the 5 keV Q-value error, a precise mass measurement of the $^{48}$Ca-$^{48}$Sc pair would tighten the $7.5\\%$ prediction without new nuclear-structure input.","The same consistent-coupling procedure could be applied to other double-beta emitters with competing single-beta branches, turning branching-ratio measurements into a systematic scan of axial quenching across the mass table."],"forward_implications":["The $^{96}$Zr $2\\nu\\beta\\beta$ matrix element becomes the first obtained in a large-scale shell-model space, and its extreme single-state dominance independently supports the earlier charge-exchange measurement.","The predicted single-beta branches of $7.5\\%$ and $18.4\\%$ are large enough that dedicated underground detectors could plausibly observe them, offering a new experimental handle on the axial-vector coupling.","The extracted effective couplings, about $0.80$ for $^{48}$Ca and about $1.04$ for $^{96}$Zr, give a consistent input for computing other low-momentum weak-transition rates in these model spaces.","If the first $1^+$ state of $^{96}$Nb is at 694.6 keV, the $^{96}$Zr matrix element becomes $0.0854$ and the $5^+$ single-beta branch drops to about $17\\%$, still within reach of detection."],"supporting_citations":[{"why":"Supplies the earlier 48Ca shell-model 2νββ matrix element of 0.0539 that this work revises.","marker":"[20]"},{"why":"Gives the previous 48Ca beta-decay branching-ratio study whose half-lives are updated here.","marker":"[6]"},{"why":"Provides the earlier smaller-space 96Zr shell-model beta-decay results used for comparison.","marker":"[8]"},{"why":"Provides the measured 48Ca 2νββ half-life used to extract the effective axial coupling.","marker":"[9]"},{"why":"Provides the measured 96Zr 2νββ half-life used to extract the effective axial coupling.","marker":"[10]"},{"why":"Supplies the phase-space integrals that turn the matrix elements into half-lives.","marker":"[25]"},{"why":"Reports the experimental extreme single-state dominance and the possible 694.6 keV first 1+ state in 96Nb.","marker":"[30]"},{"why":"Supplies the well-tested two-body interaction used for the 48Ca model space.","marker":"[21]"},{"why":"Supplies the well-tested two-body interaction used for the 96Zr model space.","marker":"[23]"}],"fun_headline_variants":["Shell model predicts 7.5% and 18.4% beta branches in 48Ca and 96Zr","Consistent large-space shell model boosts beta branch predictions for 48Ca and 96Zr","New shell-model matrix elements raise beta decay branches in 48Ca and 96Zr","One shell model explains both beta-beta and beta decay branches in 48Ca and 96Zr","Extreme single-state dominance confirmed in 96Zr shell-model matrix element"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The prediction stands on the assumption that the same quenched axial-vector coupling extracted from the measured double-beta half-lives also governs the highly forbidden single-beta transitions, because both are low-momentum-exchange processes.","fun_headline_variants_meta":{"raw":{"variants":["Shell model predicts 7.5% and 18.4% beta branches in 48Ca and 96Zr","Consistent large-space shell model boosts beta branch predictions for 48Ca and 96Zr","New shell-model matrix elements raise beta decay branches in 48Ca and 96Zr","One shell model explains both beta-beta and beta decay branches in 48Ca and 96Zr","Extreme single-state dominance confirmed in 96Zr shell-model matrix element"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000253,"raw_usage":{"total_tokens":1659,"prompt_tokens":1132,"completion_tokens":527,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":748,"completion_tokens_details":{"reasoning_tokens":409}},"tokens_in":748,"tokens_out":527,"duration_ms":50154,"temperature":1.0,"reasoning_tokens":409,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:54:33.464573+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the $5^+$ single-$\\beta$ branch of $^{96}$Zr: the paper predicts roughly $18.4\\%$ (or about $16.7\\%$ if the first $1^+$ state of $^{96}$Nb is at 694.6 keV), so a measured branch clearly outside this range would falsify the claim that the double-$\\beta$-derived coupling transfers to the single-$\\beta$ channel.","supporting_citations":[{"cited_title":"Horoi, S","cited_arxiv_id":null,"evidence_quote":"Supplies the earlier 48Ca shell-model 2νββ matrix element of 0.0539 that this work revises."},{"cited_title":"Haaranen, M","cited_arxiv_id":null,"evidence_quote":"Gives the previous 48Ca beta-decay branching-ratio study whose half-lives are updated here."},{"cited_title":"Alanssari, D","cited_arxiv_id":null,"evidence_quote":"Provides the earlier smaller-space 96Zr shell-model beta-decay results used for comparison."},{"cited_title":"Arnold, C","cited_arxiv_id":null,"evidence_quote":"Provides the measured 48Ca 2νββ half-life used to extract the effective axial coupling."},{"cited_title":"Argyriades, R","cited_arxiv_id":null,"evidence_quote":"Provides the measured 96Zr 2νββ half-life used to extract the effective axial coupling."},{"cited_title":"Neacsu, M","cited_arxiv_id":null,"evidence_quote":"Supplies the phase-space integrals that turn the matrix elements into half-lives."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports the experimental extreme single-state dominance and the possible 694.6 keV first 1+ state in 96Nb."},{"cited_title":"Honma, T","cited_arxiv_id":null,"evidence_quote":"Supplies the well-tested two-body interaction used for the 48Ca model space."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the well-tested two-body interaction used for the 96Zr model space."}],"review_version":1}