{"id":"2072c359-f511-4202-bdf5-41244833b4e8","arxiv_id":"2505.03502","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"New pn-QRPA calculations give weak rates and the first electron capture cross sections for five N=Z waiting point nuclei, showing electron capture competes with positron decay under rp-process conditions.","lead":"This paper calculates nuclear shapes and weak interaction rates for five neutron-deficient waiting point nuclei used in rp-process nucleosynthesis models. It finds that electron capture is as important as positron decay under rp-process conditions, which strengthens the case for including it in X-ray burst network calculations.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Factor-two claim for 80Zr rests on a half-life-calibrated chi and conflicts with the paper's own conclusion naming 84Mo.","rationale":"Good-faith reading: the paper is a model-based calculation of weak rates; the central contribution is the claim of a factor-two excess for 80Zr and the significance of electron capture. The weakest point is not the models themselves but the calibration chain: chi is fitted to the experimental half-lives of these very nuclei and beta is borrowed from RMF; the same half-lives are used as validation. That is a self-recognized circular step in Section 3.1, where the paper states chi 'best reproduced the experimental half-lives'. A reader cannot tell whether the factor two is physics or fitting. Figure 8 shows half-life sensitivity to beta, but there is no equivalent for chi or for stellar rates, and no error bars are given, as the Reader also noted. The second issue is the abstract-versus-conclusion mismatch on which nucleus has the factor two: the Abstract and Section 4 say 80Zr, while Section 5 says 84Mo; this is an internal inconsistency that must be resolved. Both issues are addressable: a sensitivity sweep would answer the first, and a careful comparison of the data behind Figs. 15 and 16 would answer the second. The EC/positron claim in Table 3 is less problematic: even at T9 = 1.5, EC/beta+ ratios of roughly 0.2–0.4 are nontrivial, and the paper's core statement that EC should not be neglected is consistent with the table. Because the concerns are concrete and fixable, the CONDITIONAL verdict remains appropriate; no change in verdict is recommended.","tokens_in":23260,"tokens_out":5409,"duration_ms":50987,"concrete_test":"Recompute the pn-QRPA stellar rates for 80Zr and 84Mo at T9 = 1.5 GK and rho = 1e6 g/cm^3 for chi = 3.8/A, 4.2/A, and 4.6/A, keeping beta_RMF and all other parameters fixed, and repeat with beta_RMF scaled by ±20% at chi = 4.2/A. If the 80Zr-to-Sarriguren total rate ratio does not remain near 2 (and the 84Mo ratio near 1) across these variations, the factor-two claim is a calibration artifact. Separately, extract the total-rate curves from the data underlying Figs. 15 and 16; if the factor-two excess is attached to 84Mo rather than 80Zr, the Abstract's central claim is misattributed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim — that pn-QRPA total weak rates for 80Zr exceed Skyrme HF+BCS+QRPA rates by a factor of two under rp-process conditions — is not secured against the model's calibration. In Section 3.1 the particle-hole strength chi is 'concluded' as 4.2/A because it 'best reproduced the experimental half-lives' (Refs. 60, 61), and the deformation used in the Nilsson basis is taken from RMF (Ref. 59) rather than the IBM values the paper itself computed. The calculated half-lives are then presented in Section 4 (Fig. 7, percentage deviations) as validation. This is a circular verification: the same data used to fix chi and select beta are used to claim agreement with experiment. No sensitivity study is given for chi or for the RMF deformation, so it is unknown whether the factor-two excess for 80Zr is a robust model prediction or a by-product of the calibration. The concern is concrete: the factor-two statement in the Abstract and Section 4 would have to be confirmed when chi and beta are varied within uncertainties. Compounding this, the Abstract and Section 4 identify 80Zr as the factor-two nucleus, whereas Section 5 (Conclusions) states that 'total weak rates are a factor two bigger ... for 84Mo' and that the remaining nuclei, which would include 80Zr, agree. This internal contradiction directly affects which result is being claimed and must be resolved before the central claim can be evaluated.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript examines the structure and weak-interaction properties of the five N=Z rp-process waiting-point nuclei 80Zr, 84Mo, 88Ru, 92Pd, and 96Cd. The authors fit the IBM-1 Hamiltonian parameters to experimental level schemes (and to a shell-model calculation for 96Cd), compute potential energy surfaces and shape predictions, and then use a deformed Nilsson-basis pn-QRPA model to calculate GT+ strength distributions, terrestrial half-lives, stellar positron-decay and continuum electron-capture rates, electron-capture cross sections, and energy rates and emission probabilities of beta-delayed protons. The stellar rates are compared with Sarriguren's Skyrme HF+BCS+QRPA results; the abstract claims that the total weak rate for 80Zr is twice the Skyrme rate under rp-process conditions and that electron capture competes significantly with positron decay. The paper's headline result is clouded by an internal contradiction: the Abstract and Section 4 attribute the factor-two excess to 80Zr, while Section 5 attributes it to 84Mo; additionally, the particle-hole strength chi and the adopted deformations are calibrated on the same half-lives that are later presented as validation.","tokens_in":23602,"tokens_out":18102,"duration_ms":156600,"significance":"If the results hold, a factor-two change in the total weak rate for one of the heavier N=Z waiting-point nuclei and the finding that continuum electron capture reaches parity with positron decay at relevant densities (Table 3) are the kind of input that could matter for rp-process network calculations in Type I X-ray burst models; the electron-capture cross sections for these nuclei are presented as a new observable. The manuscript is transparent about its model choices, stating chi=4.2/A, kappa=0.1 MeV, the GT quenching factor f_q^2=0.6, and the RMF deformations explicitly, and the rate formalism follows the standard pn-QRPA equations (Eqs. 13-26). These are genuine strengths, as is the systematic comparison with a previous HF+BCS+QRPA calculation. The significance is capped, however, by the unresolved abstract-versus-conclusions contradiction and by the absence of a sensitivity analysis for the calibrated parameters, so the factor-two claim is not yet robustly established.","major_comments":[{"comment":"The central claim of the paper is internally inconsistent. The Abstract states that 'the calculated total weak rates are twice the Skyrme HF+BCS+QRPA rates for 80Zr,' and the Section 4 discussion of Fig. 15 repeats this for 80Zr, while the same section says that for 84Mo 'the total rates are in overall good agreement with Sarriguren's rates.' Section 5, however, concludes that 'Our total weak rates are a factor two bigger than the Skyrme HF+BCS+QRPA calculated rates for 84Mo' and that the remaining nuclei, which would include 80Zr, agree well. The manuscript must state consistently which nucleus or nuclei exhibit the factor-two excess and reconcile the Abstract, Section 4, and Section 5; as written, the abstract-level result cannot be evaluated.","section":"Abstract; §4 (Figs. 15-16); §5"},{"comment":"The validation of the rate calculation is partly circular, so the headline factor-two result is not yet secured. Section 3.1 fixes chi=4.2/A as the value that 'best reproduced the experimental half-lives' (Refs. 60, 61) and adopts the RMF deformations of Ref. 59; Section 4 then uses those same experimental half-lives in Fig. 7 as evidence of agreement, and Fig. 8 shows that the calculated half-lives depend strongly on the deformation parameter. Because the factor-two total-rate excess for 80Zr (Fig. 15) and the EC/β+ ratios of Table 3 come from the same calibrated calculation, a concrete sensitivity test is needed: please vary chi over a plausible range (for example 3.5-5)/A and beta within the RMF/IBM uncertainties, and report how the 80Zr excess and the Table 3 ratios change.","section":"§3.1; §4, Figs. 7-8"},{"comment":"The deformation input is load-bearing and conflicts with the paper's own structure calculation. Table 1 lists β_IBM=0.87 for 80Zr while the adopted RMF value is 0.437, and β_IBM=0.50 (prolate) for 84Mo while the adopted RMF value is -0.247 (oblate); these RMF values determine the Nilsson basis for all GT-strength and rate results. Since Section 4 (Fig. 8) demonstrates strong half-life sensitivity to beta, and since the GT strength distributions in Figs. 5-6 differ markedly from Sarriguren's, the claimed rate excess depends on this external input. Please either justify the RMF deformations quantitatively for these soft N=Z nuclei, or perform the rate calculation with the IBM-1 deformations (or a range of beta) and discuss the effect of the factor-of-3-5 difference between β_IBM and the geometric-model beta on the final rates.","section":"§3.1; Table 1"}],"minor_comments":[{"comment":"The text attributes the '(Ee − Q + Ei − Ef)^2 factor' to Eq. (14), but this factor appears in Eq. (13), the electron-capture cross-section formula; the cross-reference should be corrected.","section":"§4, paragraph on Fig. 9"},{"comment":"The caption refers to the 'pn-QRPA(N) model,' but the notation (N) is not defined anywhere in the text; the authors should either define it or remove it.","section":"Table 3"},{"comment":"The last term of Eq. (9) is missing a closing parenthesis: the printed 'δ(pf_2, p_i_2 < n^f_2c | ... >' should read 'δ(pf_2, p_i_2) < n^f_2c | ... >'; as printed, the equation cannot be parsed unambiguously.","section":"§3.1, Eq. (9)"},{"comment":"The shaded 'rp-process' temperature window in Figs. 15-19 is never defined quantitatively in the text or captions; stating the adopted T9 range would allow the reader to reconcile the abstract's claim that electron capture rates 'compete well' with positron-decay rates against Table 3, where the EC/β+ ratio reaches unity only at T9 ≈ 4.5-5.5 at ρ = 10^6 g/cm^3.","section":"§4, Figs. 10-19"},{"comment":"For 84Mo, the calculated half-life deviates by 71% from the measured value, yet the text concludes that the calculated half-lives are 'in good agreement' with experiment; given that Section 5 names 84Mo as the nucleus with the factor-two rate excess, the authors should soften this wording or explain the origin of the 84Mo deviation.","section":"§4, Fig. 7 paragraph"},{"comment":"The procedure of replacing theoretical levels with experimental levels when they lie within 0.5 MeV is described, but the manuscript does not state how many levels were replaced for each of the five nuclei; a brief statement would improve the reproducibility of the stellar-rate calculation.","section":"§3.1"},{"comment":"The authors acknowledge that 'Collective states cannot be treated in the current pn-QRPA model,' and they mitigate this by inserting experimental levels within 0.5 MeV. Since Section 4 classifies 84Mo as O(6)-like (γ-unstable) and 88Ru as close to the E(5) critical point, the static-deformation assumption is least reliable precisely for these soft nuclei; a short discussion of how shape fluctuations might affect the central rate comparison would strengthen the paper.","section":"§3.1; §4"}],"recommendation":"major_revision","confidential_remarks":"For the editor: the manuscript is within the scope of the journal and the work is honest in disclosing its parameter choices, but the headline claim is currently self-contradictory (Abstract and §4 name 80Zr, §5 names 84Mo), and the calibration-to-validation loop for chi and beta is a genuine robustness gap. I believe both issues are fixable within the manuscript's scope, so I recommend major revision rather than rejection. I do not see evidence of a citation or attribution problem; the comparison with Sarriguren's work is appropriately credited."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a serviceable continuation of the authors' pn-QRPA program, and it does add genuinely new numbers—first electron capture cross sections for 80Zr, 84Mo, 88Ru, 92Pd, 96Cd—but the headline result is internally inconsistent. The abstract and Section 4 say the total weak rates are twice the Skyrme HF+BCS+QRPA rates for 80Zr; the conclusion says the factor-two excess is for 84Mo and that the remaining nuclei, which would include 80Zr, agree. That is not a cosmetic slip; it decides what the paper actually claims. A referee should ask for a corrected statement before anything else.\n\nWhat is good: the IBM-1 fits and PES shapes look reasonable; the GT strength distributions are compared with Sarriguren's; the half-life comparison (deviations 21%, 71%, 17%, 2%, 17%) is a real check; and the electron capture cross sections and beta-delayed proton emission rates are new for this set. The general conclusion that continuum EC competes with positron decay under rp-process conditions is consistent with Sarriguren and supports including EC in network calculations.\n\nSoft spots, in order. First, the circularity concern is real: chi=4.2/A was chosen to best reproduce the experimental half-lives of these same nuclei, and those very half-lives are then presented as validation in Fig. 7. Second, the deformation is taken from RMF rather than the IBM values the paper itself computed; the decision is explained, but no sensitivity study is given for chi or beta. Third, there are no error bars on any of the rate tables. None of these, individually, kills the paper; they do mean the factor-two excess—for whichever nucleus it is—has not been shown to be robust.\n\nThe comparison with Sarriguren's rates is useful even though the two calculations differ in quenching and model space. The paper deserves a serious referee, but it should not be accepted as is. I would not cite the factor-two claim until the contradiction is resolved and a sensitivity test is added.","headline":"Useful new pn-QRPA rates for five N=Z waiting-point nuclei, undercut by an internal contradiction over which nucleus shows the factor-two excess and an unresolved calibration circularity.","tokens_in":24144,"tokens_out":3666,"would_cite":false,"duration_ms":34731,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["21.10.-k","21.60.Jz","21.60.Ev","21.60.Fw","23.40.-s","26.30.Jk","26.50.+x"],"model":"deepseek-v4-flash","headline":"The paper claims that the rp-process waiting-point nucleus $^{80}$Zr has stellar weak rates twice the previous QRPA values, with electron capture competing strongly with positron decay.","keywords":["nuclear structure","weak rates","IBM-1","pn-QRPA","rp-process","waiting point nuclei","N=Z nuclei","electron capture"],"falsifier":"Measure the Gamow-Teller strength distribution of $^{80}$Zr, for example with high-resolution charge-exchange or total-absorption $\\beta$-decay spectroscopy, and compare the total low-lying $B(GT^+)$ with the pn-QRPA prediction; if the measured distribution is closer to the Skyrme-QRPA one, the factor-two rate excess and the electron-capture dominance would not survive.","tokens_in":23051,"feed_emoji":"☢️","tokens_out":10536,"duration_ms":93194,"temperature":0.7,"pith_summary":"The paper argues that the weak-interaction rates governing the rp-process through the $N=Z$ waiting-point nuclei $^{80}$Zr, $^{84}$Mo, $^{88}$Ru, $^{92}$Pd, and $^{96}$Cd should be revised for $^{80}$Zr: under rp-process conditions the total stellar weak rate is about twice the rate from a previous Skyrme HF+BCS+QRPA calculation, while the other four nuclei agree within the usual tolerance. The same calculation finds that continuum electron capture is not a negligible correction: as density and temperature rise, capture rates grow until they exceed positron decay rates by up to a factor of five. If these results are right, rp-process network calculations for explosive hydrogen burning should include electron capture and use higher weak rates at $^{80}$Zr.","feed_headline":"For 80Zr, stellar weak rates are twice the previous QRPA values","feed_subtitle":"Electron capture rivals positron decay, so rp-process networks should include it.","key_machinery":"The engine is the proton-neutron quasiparticle random phase approximation (pn-QRPA), which generates Gamow-Teller strength distributions and the stellar rates from a deformed Nilsson basis; the basis uses quadrupole deformations from a relativistic mean-field calculation, with BCS pairing and a separable Gamow-Teller residual interaction. The particle-hole strength is fixed at $\\chi=4.2/A$ to reproduce terrestrial half-lives. The interacting boson model (IBM-1), with Hamiltonian $\\epsilon \\hat{n}_d + a_2 \\hat{Q}\\cdot\\hat{Q}$, is used only for energy levels and potential energy surfaces that classify each nucleus as spherical or deformed. The pn-QRPA's access to many (up to seven) oscillator shells is what lets the rate sum converge at high temperature.","core_discovery":"The central discovery is that the Gamow-Teller strength of the self-conjugate waiting-point nucleus $^{80}$Zr is much more fragmented and larger in the pn-QRPA model than in the Skyrme HF+BCS+QRPA model, and this extra low-lying strength doubles the computed stellar weak rate under rp-process conditions. For $^{84}$Mo, $^{88}$Ru, $^{92}$Pd, and $^{96}$Cd the total rates from the two models compare well, although the individual positron-decay and capture contributions differ. A separate result is that continuum electron capture is significant throughout the rp-process parameter range and becomes the dominant weak channel at $\\rho=10^7$ g cm$^{-3}$, up to five times the positron-decay rate. The paper takes this as evidence that electron capture belongs in rp-process weak-rate input.","pith_inferences":["If the factor-two excess for $^{80}$Zr is confirmed by direct strength data, rp-process network simulations would need rerunning; the most testable external consequences are X-ray burst light curves and the ash composition of the neutron-star crust.","Because the authors chose relativistic mean-field deformations instead of their own IBM-1 values, which are three to five times larger, one could recompute the same rates with the IBM-1 deformations to isolate how much of the $^{80}$Zr excess is shape-driven.","A shell-model diagonalization in the $A=80$--$100$ region would test whether the fragmented strength pattern is a robust physical feature or an artifact of the QRPA treatment; no such comparison is made in the paper.","The highest temperatures in the grid (up to 30 GK) go beyond rp-process conditions; those rates may be relevant to other astrophysical sites such as core-collapse supernovae or accretion disks, but the paper does not make that application."],"forward_implications":["If the $^{80}$Zr rate is indeed twice the Skyrme-QRPA value, rp-process network calculations that use the older rate will underestimate how quickly material clears the $^{80}$Zr waiting point, shifting the abundance flow toward heavier masses.","Continuum electron capture should be included in rp-process weak-rate input; at $\\rho=10^7$ g cm$^{-3}$ it exceeds positron decay by up to a factor of five, so omitting it changes the effective lifetime and heating of the flow.","The predicted electron-capture cross sections rise by roughly two orders of magnitude between $T=0.5$ and $1.0$ MeV because of thermal unblocking of the Gamow-Teller channel, so the temperature dependence of the thermal population matters, not just the ground-state strength.","The reported beta-delayed proton energy rates and emission probabilities for these five waiting-point nuclei provide new input that can alter the composition of X-ray burst ashes and, through them, models of neutron-star crusts."],"supporting_citations":[{"why":"Supplies the relativistic mean-field deformation parameters used in the Nilsson basis for the QRPA calculation.","marker":"(59)"},{"why":"Provides the pn-QRPA formalism for the Gamow-Teller matrix elements used in the rate calculation.","marker":"(38)"},{"why":"Establishes the particle-hole strength calibration and the beta-plus and electron-capture half-life systematics that set $\\chi=4.2/A$.","marker":"(40)"},{"why":"Supplies the Skyrme HF+BCS+QRPA comparison rates, half-lives, and Gamow-Teller strength distributions.","marker":"(66)"},{"why":"Supplies the AME2012 experimental masses and half-lives used for fitting and validation.","marker":"(60)"},{"why":"Supplies the experimental level data used to fit the IBM-1 Hamiltonian parameters for the known waiting-point nuclei.","marker":"(22)"},{"why":"Supplies the shell-model levels used to fit the IBM-1 parameters for the unknown nucleus $^{96}$Cd.","marker":"(23)"},{"why":"Previous work in the same mass region that fixes the present parameter choice and calculation procedure.","marker":"(21)"}],"fun_headline_variants":["80Zr weak rates double: electron capture matters in rp-process","Electron capture dominates weak rates for rp-process nuclei","80Zr weak rates double earlier QRPA; electron capture rivals decay","rp-process: 80Zr weak rates double, electron capture competes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The rates stand on adopting quadrupole deformations from a relativistic mean-field calculation rather than the authors' own IBM-1 values, and on tuning the particle-hole interaction strength to reproduce terrestrial half-lives of these same nuclei; if the true shapes of these $N=Z$ nuclei differ, or if half-life tuning does not carry over to stellar conditions, the factor-two excess for $^{80}$Zr would change.","fun_headline_variants_meta":{"raw":{"variants":["80Zr weak rates double: electron capture matters in rp-process","Electron capture dominates weak rates for rp-process nuclei","80Zr weak rates double earlier QRPA; electron capture rivals decay","rp-process: 80Zr weak rates double, electron capture competes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000844,"raw_usage":{"total_tokens":3715,"prompt_tokens":1022,"completion_tokens":2693,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":638,"completion_tokens_details":{"reasoning_tokens":2619}},"tokens_in":638,"tokens_out":2693,"duration_ms":18295,"temperature":1.0,"reasoning_tokens":2619,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:50:59.471329+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the Gamow-Teller strength distribution of $^{80}$Zr, for example with high-resolution charge-exchange or total-absorption $\\beta$-decay spectroscopy, and compare the total low-lying $B(GT^+)$ with the pn-QRPA prediction; if the measured distribution is closer to the Skyrme-QRPA one, the factor-two rate excess and the electron-capture dominance would not survive.","supporting_citations":[],"review_version":1}