{"id":"b79562be-7853-4b92-abc8-83fed8180354","arxiv_id":"2511.21567","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Adding microscopic beta-decay rates to supernova simulations makes the pre-bounce antineutrino signal orders of magnitude brighter than positron-capture-only models.","lead":"This paper ran core-collapse supernova simulations with a new microscopic beta-decay rate table, finding much stronger pre-bounce antineutrino emission than models with only positron capture. Future neutrino detectors could use this enhanced signal to probe the last seconds of a massive star's life.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"FT-QRPA β-decay rates are unbenchmarked and exceed shell-model rates by up to 3 dex; the claimed 4-dex emissivity enhancement could be an artifact of overcounted de-excitation strength.","rationale":"The reader's weakest_assumption identified the FT-QRPA rate library as the load-bearing premise, and my independent reading agrees. The paper's own Supplemental Material is explicit that FT-QRPA rates exceed shell-model rates by up to three orders of magnitude, and this discrepancy is not tied to any measured quantity. Because the enhanced antineutrino emissivity and luminosity are direct linear functions of these rates, an overcount in the de-excitation channel undermines the central numerical claim. The lack of a public rate table compounds the issue, as no independent reproduction is possible. The inconsistency between the abstract and introduction further indicates the enhancement factors are not well constrained. However, the qualitative conclusion—that beta decay on neutron-rich nuclei, especially around N=50, can produce high-energy antineutrinos and appreciably contribute to the pre-bounce signal—does not hinge on the exact factor and is supported by the physical argument that these nuclei have larger Q-values. Thus the appropriate verdict remains CONDITIONAL, consistent with the reader.","tokens_in":12659,"tokens_out":6920,"duration_ms":57870,"concrete_test":"Isolate the de-excitation contribution by recomputing the FT-QRPA β-decay rate for a representative set of pf-shell and N≈50 nuclei (e.g., 76Fe, 78Ni, 82Ge) with the parent ground state only, and then with the full finite-temperature ensemble; compare both to shell-model rates and to any available experimental β-strength distributions (e.g., from total absorption spectroscopy). If the full-ensemble rate exceeds the ground-state rate by orders of magnitude in the temperature range 5–12 GK, verify the transition strengths against shell-model or experimental data. Additionally, re-run the GR1D simulation with the de-excitation contributions set to zero to see whether the factor-of-10^4 emissivity enhancement and the 11-MeV tail persist; if they disappear, the central claim is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim—'antineutrino emissivity increases by more than 4 orders of magnitude'—rests entirely on the FT-QRPA β-decay rate table generated by the authors. The paper provides no experimental validation of these rates. Supplemental Fig. S2 shows that where FT-QRPA can be compared with the shell-model rates of Refs. [26,47], FT-QRPA rates are systematically larger by up to three orders of magnitude; the authors attribute this to de-excitations from highly excited parent states (Supplemental Material, 'Comparison of FT-QRPA and shell-model β decay rates'). No measured β-decay half-lives or β-strength functions are presented to support this de-excitation component. Since the emissivity in Eq. (4) sums antineutrino contributions weighted by NSE abundances, an overcount of the de-excitation strength would propagate linearly into the emissivity and luminosity, potentially collapsing the claimed enhancement to the much smaller shell-model enhancement. The abstract (two orders of magnitude/factor of 50) and introduction (4 orders/3 orders) also report inconsistent enhancement factors, suggesting the headline numbers have not been stabilized. The rate table is only 'available from the authors upon request,' not publicly archived, which impedes independent verification. Until the FT-QRPA β-decay rates are benchmarked or the de-excitation contribution is isolated and validated, the specific enhancement factors should be treated as conditional.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This Letter presents, to the authors' knowledge, the first core-collapse supernova (CCSNe) simulation that incorporates global β-decay rates from a microscopic theory. The rates are computed with a relativistic energy density functional (D3C*) and finite-temperature quasiparticle random-phase approximation (FT-QRPA) for 902 nuclei from Z=20 to 40, including allowed Gamow-Teller and first-forbidden transitions, with full temperature and density coverage. The antineutrino spectral function is parameterized with an effective Q-value ansatz (Eq. 3), and the rates are implemented in the open-source code GR1D using NSE composition from the SFHo equation of state. The authors compare models with only positron capture on free neutrons (e+ + n baseline), with shell-model rates (limited to the pf-shell), and with the new FT-QRPA rates. They report that β decay enhances the pre-bounce antineutrino emissivity by more than four orders of magnitude and the luminosity by about three orders of magnitude relative to the e+ + n baseline, extending the antineutrino spectrum to about 11 MeV, while leaving the collapse dynamics essentially unchanged.","tokens_in":13052,"tokens_out":6858,"duration_ms":62464,"significance":"If the FT-QRPA rate table is reliable, the paper is a substantial step forward: it provides a global microscopic β-decay rate library for CCSNe simulations extending to neutron-rich nuclei near N=50, and it quantifies a potentially observable pre-bounce antineutrino signal. The qualitative direction—that β decay enhances antineutrino emission—is supported by both the FT-QRPA and the shell-model rates, and the use of a self-consistent framework for electron capture and β decay is a strength. However, the quantitative enhancement factors are not stable (the abstract, introduction, and conclusions disagree by up to two orders of magnitude) and rest on FT-QRPA de-excitation strength that exceeds shell-model rates by up to three orders of magnitude without experimental validation. The paper is clearly written and uses open-source tools (GR1D, NuLib), but the central numerical claims require additional benchmarking and reconciliation before they can be accepted.","major_comments":[{"comment":"The enhancement factors are internally inconsistent. The abstract reports antineutrino emissivity enhancement 'more than two orders of magnitude' and luminosity 'by a factor of 50'; the introduction reports 'more than 4 orders of magnitude' for emissivity and '3 orders of magnitude' for luminosity; the conclusions state '~10^3 times' luminosity. Since these are the paper's headline quantitative claims, the authors must reconcile them and specify the exact thermodynamic epoch and baseline (e+ + n) used for each ratio. This inconsistency suggests the headline numbers are not yet stable.","section":"Abstract vs. Introduction vs. Conclusions"},{"comment":"The central quantitative enhancement is driven by FT-QRPA β-decay rates that are not benchmarked to experiment. Fig. S2 shows FT-QRPA rates exceed shell-model rates by up to three orders of magnitude for pf-shell nuclei, attributed to de-excitations from highly excited parent states (Ref. [42]). No measured β-decay half-lives or β-strength functions are presented to validate this de-excitation component. Since the emissivity in Eq. (4) is linear in the rates, an overcount of de-excitation strength would propagate directly into the claimed 4-order enhancement. The authors should provide experimental validation (e.g., half-lives at low T) and/or isolate the de-excitation contribution to demonstrate that it is not an artifact. Until then, the specific enhancement factors should be treated as conditional.","section":"Supplemental Fig. S2; Methods; Eq. (4)"},{"comment":"The antineutrino spectrum n_i(E) is a one-parameter ansatz with effective Q-value q_i fitted to reproduce the average energy. While the total luminosity is independent of the shape because Eq. (2) defines λ̄ν = ∫ E n_i dE, the differential emissivity and the claimed spectral extension to ~11 MeV (Fig. 1(a)) are not derived from the microscopic calculation. The authors should validate the shape against a microscopic spectral function for representative nuclei, or soften the high-energy-tail claim.","section":"Eq. (3) and Fig. 1"}],"minor_comments":[{"comment":"The β-decay rate table is only available 'upon request.' For reproducibility and to enable community use, the table should be deposited in a public repository (e.g., Zenodo) or included as supplementary material.","section":"Data availability"},{"comment":"The abstract describes the baseline as 'thermal emission alone,' but the baseline model is positron capture on free neutrons (e+ + n), not thermal emission. Please correct the wording.","section":"Abstract"},{"comment":"The effective Q-value q_i is chosen to reproduce the average antineutrino energy; this is a fit parameter. The authors should state how many parameters are in the spectral model and whether results are sensitive to the chosen ansatz.","section":"Methods, Eq. (3)"},{"comment":"The notation D3C* should be defined (it is the relativistic EDF parametrization). Also, the paper would benefit from a sentence explaining the de-excitation treatment in FT-QRPA, since Fig. S2 shows it is the main source of rate differences.","section":"Methods"},{"comment":"The dotted curve in Fig. 3 is labeled 'FT-QRPA νe' in the legend and 'FT-QRPA antineutrinos' in the caption; check that all labels are consistent.","section":"Fig. 3"},{"comment":"The label 'log10(¯ν emissivity)' in Fig. S3 is ambiguous; the text refers to antineutrino energy loss rate. Please unify the terminology.","section":"Supplemental Material, Fig. S3"}],"recommendation":"major_revision","confidential_remarks":"The paper's central quantitative claim is sensitive to unvalidated de-excitation strength in the FT-QRPA rates. I recommend requiring experimental benchmarking or a sensitivity analysis before publication. The internal inconsistency in the headline numbers must be fixed. If these points are addressed, the paper would be a valuable contribution to the field."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear Colleague,\n\nThe genuinely new thing here is the first CCSN simulation that plugs in a global microscopic beta-decay table: 902 nuclei from Ca to Zr, EDF+FT-QRPA, including de-excitations and first-forbidden transitions, alongside self-consistent EC rates from the same framework. That is solid work and a real extension of the group's electron-capture program. The qualitative conclusion is also credible: beta decay makes pre-bounce antineutrinos far brighter than positron capture alone, and the shell-model comparison supports the direction, not just the FT-QRPA model.\n\nThe quantitative size of the effect, though, is not stable. The abstract gives two orders of magnitude in emissivity and a factor of 50 in luminosity; the introduction gives four and three orders. That's a red flag that the headline numbers haven't converged. The rate table is only 'available upon request,' and there's no benchmark against measured beta-decay half-lives. In the overlap region with shell-model rates, FT-QRPA is systematically higher by up to three orders, and the paper attributes this to de-excitations from highly excited parent states without isolating and validating that component. If that strength is overcounted, the enhancement collapses to the shell-model level. The spectral ansatz with an effective Q-value fitted to the average energy is fine, but it's another reason to treat the specific spectrum as indicative.\n\nNone of this kills the paper. The direction is almost certainly right, and the N=50 neutron-rich nuclei contribution is a sensible, testable prediction. But the specific enhancement factors should be treated as conditional until the rates are benchmarked and the table is public.\n\nWho gets value? Supernova neutrino phenomenologists and nuclear weak-rate groups. It deserves a serious referee, not a desk rejection, because it's the first result of this kind. I'd send it out, expecting the referee to demand a public rate table, an experimental half-life benchmark, and one consistent set of numbers between abstract and body.\n\nBest, [name]","headline":"First global microscopic beta-decay table in a CCSN simulation; qualitative effect is credible, but the headline enhancement factors are not yet trustworthy because the rates are unbenchmarked and the paper's own numbers disagree internally.","tokens_in":13504,"tokens_out":2844,"would_cite":true,"duration_ms":28455,"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":"Including microscopic β-decay rates magnifies the pre-bounce antineutrino luminosity of a core-collapse supernova by a factor of roughly a thousand, with the spectrum reaching about 11 MeV.","keywords":["core-collapse supernovae","antineutrino emission","beta decay","weak-interaction rates","finite-temperature QRPA","energy density functional","N=50 nuclei","pre-supernova neutrinos"],"falsifier":"Measure the β-decay half-lives (or Gamow-Teller strength) of the neutron-rich N≈50 nuclei that dominate the emissivity—around zinc–germanium isotopes—and compare with FT-QRPA predictions; a systematic overprediction of the size seen in the shell-model comparison (~10^3) would falsify the enhancement. Alternatively, a Galactic pre-supernova antineutrino burst whose flux matches the shell-model expectation rather than the FT-QRPA one.","tokens_in":12579,"feed_emoji":"⚛️","tokens_out":5283,"duration_ms":53761,"temperature":0.7,"pith_summary":"The paper claims that β decay on neutron-rich nuclei—not just positron capture on free neutrons—produces the dominant electron-antineutrino emission in the seconds before a core-collapse supernova's bounce. Using a microscopic nuclear model (relativistic energy-density functional plus finite-temperature quasiparticle random-phase approximation), the authors compute β-decay rates for 902 nuclei and feed them into supernova simulations. They report that the pre-bounce antineutrino emissivity rises by more than four orders of magnitude and the luminosity by about three orders of magnitude compared with the standard baseline, with the spectrum extending to roughly 11 MeV. If correct, this makes the presupernova antineutrino signal a far stronger target for neutrino detectors and gives β decay a potential role in late-stage stellar evolution, even though the collapse dynamics themselves are unchanged.","feed_headline":"Beta decay lifts pre-supernova antineutrino luminosity 1000-fold","feed_subtitle":"Microscopic rate table shows neutron-rich nuclei near N=50 dominate the pre-bounce signal.","key_machinery":"The central object is the FT-QRPA β-decay rate table: β-decay and antineutrino-energy-loss rates for 902 nuclei from calcium to zirconium, computed at temperatures 1–30 GK and densities ρYe = 10^5–10^12 g/cm^3. It combines relativistic energy-density-functional theory (D3C* interaction) for the nuclear ground state with the finite-temperature quasiparticle random-phase approximation for transitions to excited states, including first-forbidden transitions and de-excitations from highly excited parent states. Coupled with the NSE composition and the neutrino-transport library, this table replaces the usual assumption of only positron capture on free neutrons as an antineutrino source.","core_discovery":"The discovery is that β decay, included self-consistently via FT-QRPA rates, enhances the pre-bounce antineutrino emissivity by more than four orders of magnitude and the luminosity by more than three orders relative to positron-capture-only estimates. The enhancement peaks at about 3 MeV and extends to about 11 MeV, with neutron-rich nuclei near N=50 contributing most of the high-energy tail. This is the first core-collapse supernova simulation to incorporate global β-decay rates from a microscopic theory, and it leaves the electron-neutrino channel and the collapse dynamics essentially unchanged.","pith_inferences":["The FT-QRPA rates exceed shell-model rates by up to three orders of magnitude in the overlap region (supplementary Fig. S2); a conservative reading would treat the shell-model result as a lower bound and the FT-QRPA result as an upper bound on the true enhancement.","If the de-excitation contributions are verified, the same mechanism could also enhance electron-neutrino emission at even higher densities and alter the presupernova antineutrino background relevant to early-warning systems.","A direct test is to compare FT-QRPA β-decay half-lives or Gamow-Teller strengths against terrestrial measurements on N≈50 nuclei (e.g., zinc and germanium isotopes); a systematic discrepancy on the scale seen in the shell-model comparison would overturn the enhancement claim.","The paper focuses on pre-bounce signals; extending these rates to post-bounce or to rotating/3D models could reveal larger dynamical effects than the spherical, non-rotating case considered here."],"forward_implications":["Pre-bounce antineutrino luminosity is roughly 10^3 times the positron-capture-only value, bringing it within about two orders of magnitude of the electron-neutrino luminosity.","The antineutrino spectrum extends to about 11 MeV, in the efficient detection window of water-Cherenkov and scintillator detectors, so a Galactic pre-supernova could yield a much larger burst than previously expected.","Reliable antineutrino predictions require rate tables covering medium-heavy neutron-rich nuclei (N≈50), beyond the pf-shell region accessible to shell-model calculations.","Although the new rates do not alter the collapse trajectory in these simulations, they set the stage for changes in presupernova stellar evolution, which could feed back into collapse initial conditions.","The same FT-QRPA framework gives, for the first time, average antineutrino and neutrino energies from a consistent rate set, enabling direct channel-by-channel comparison."],"fun_headline_variants":["β decay drives 1000× antineutrino luminosity jump in CCSNe","Microscopic β-decay rates boost CCSN antineutrinos 1000×","Self-consistent β decay makes pre-bounce antineutrinos shine 1000× brighter","First CCSN run with microscopic β decay reveals 1000× antineutrino boost","β decay raises CCSN antineutrino luminosity by 3 orders of magnitude"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The premise that the FT-QRPA β-decay rate table, especially the de-excitation contributions from highly excited parent states, is quantitatively accurate; if those rates are overstated, the orders-of-magnitude enhancement shrinks to the much smaller shell-model level.","fun_headline_variants_meta":{"raw":{"variants":["β decay drives 1000× antineutrino luminosity jump in CCSNe","Microscopic β-decay rates boost CCSN antineutrinos 1000×","Self-consistent β decay makes pre-bounce antineutrinos shine 1000× brighter","First CCSN run with microscopic β decay reveals 1000× antineutrino boost","β decay raises CCSN antineutrino luminosity by 3 orders of magnitude"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001396,"raw_usage":{"total_tokens":5475,"prompt_tokens":729,"completion_tokens":4746,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":473,"completion_tokens_details":{"reasoning_tokens":4643}},"tokens_in":473,"tokens_out":4746,"duration_ms":30982,"temperature":1.0,"reasoning_tokens":4643,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T19:57:18.509363+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the β-decay half-lives (or Gamow-Teller strength) of the neutron-rich N≈50 nuclei that dominate the emissivity—around zinc–germanium isotopes—and compare with FT-QRPA predictions; a systematic overprediction of the size seen in the shell-model comparison (~10^3) would falsify the enhancement. Alternatively, a Galactic pre-supernova antineutrino burst whose flux matches the shell-model expectation rather than the FT-QRPA one.","supporting_citations":[],"review_version":1}