{"id":"c289fc47-834a-4fb3-a7f7-63808bcef6b8","arxiv_id":"2411.13123","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Brink-Axel hypothesis fails badly for stellar beta-decay rates of fp/fpg-shell nuclei in a pn-QRPA model, with deviations up to three orders of magnitude.","lead":"Weak interaction rates that drive supernova evolution are often computed by assuming the Brink-Axel hypothesis, which says excited nuclei respond to transitions just like ground states. Using a microscopic quasiparticle random-phase approximation model, this paper finds that approximation badly underestimates or overestimates beta-decay and electron-capture rates in heavy nuclei, especially in beta decay at high temperature and density.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Unvalidated pn-QRPA excited-state GT strengths are the benchmark for the BA test; without an independent check, the three-order and neon-burning claims overstate a model-internal result.","rationale":"The paper's internal logic is straightforward and the numerical comparisons are consistent: within the pn-QRPA model, adopting BA changes the rates, particularly for BD at high temperature and density. This is a valid model-internal statement. The paper also compares ground-state GT strength functions with data (Figures 2 and 3) and state-by-state rates with LSSM (Tables 8 and 9), which is real support for the low-energy part of the model. However, the central claim—that BA is not appropriate for reliable stellar weak-rate calculations—depends on the excited-state strength functions being a trustworthy benchmark. The manuscript validates only ground-state strengths; it provides no experimental or independent theoretical benchmark for the excited-state distributions. The acknowledged zero-temperature, schematic nature of the model (Eq. 37 and the limitation statements in the Introduction and Section 2) makes this gap consequential: the same model that reproduces ground-state strengths could misplace the excitation-energy dependence of GT strength, which is exactly the quantity under test. The abstract and conclusion also overstate the onset and magnitude of deviations relative to Table 4, undercutting the rhetorical claim even though the direction is plausible. These issues are not fatal: the paper stands as a useful model-based warning. They do, however, justify a CONDITIONAL verdict rather than an unqualified one. A single external check using a finite-temperature QRPA or shell model for the test nuclei would settle whether the BA violation is physical or an artifact of the pn-QRPA excited-state spectroscopy.","tokens_in":31214,"tokens_out":8498,"duration_ms":84010,"concrete_test":"Recompute the GT strength functions from the first few excited states of 58Cr, 67Co, and 70Cu with an independent model, e.g., the finite-temperature QRPA of Dzhioev et al. (2020) or a large-scale shell-model calculation where feasible, and compare B(GT) distributions and the resulting λQRPA in Tables 4-6. If the independent excited-state strengths differ substantially from the pn-QRPA ones, the BA violation is model-dependent and the central claim must be downgraded to a model-specific result. A lighter check: rerun the pn-QRPA calculation with the GT force constants χGT and κGT varied over the range used in prior pn-QRPA work; if the deviation pattern changes qualitatively, the conclusion lacks robustness.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central comparison in Section 2 tests BA by replacing the model's own excited-state GT strength distributions with ground-state (or low-lying) distributions and comparing the resulting weak rates. The benchmark rates λQRPA therefore inherit all errors in the pn-QRPA excited-state strengths. Only ground-state strengths are validated against experiment (Figures 2 and 3); no excited-state strength is compared to data or to an independent calculation. The model itself is acknowledged to be a zero-temperature approximation with a schematic pairing-plus-quadrupole Hamiltonian and separable GT forces (Eq. 37), and the force constants χGT and κGT are never specified (Eqs. 5-8), so the excited-state distributions are both unanchored and irreproducible. If the true excited-state GT strength is closer to the ground-state strength than the model suggests, the apparent BA violation shrinks or disappears. Additionally, the abstract's 'as early as neon burning' and 'three orders of magnitude' claims outrun the tables: Table 4 shows ΔGBE ≈ 0 for 58Cr at T = 3 GK and low densities, and the three-order deviation occurs only at T = 30 GK and ρYe = 10^11 g cm^-3, an extreme corner where the zero-temperature approximation is least justified. The conclusion's own threshold 'exceeding 3 GK and densities beyond 10^6 g cm^-3' is also not comfortably identified with neon burning. Hence the central claim is plausible as a model-internal warning but is not yet a demonstrated general failure of BA for real nuclei.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript investigates whether the Brink-Axel (BA) hypothesis can be used to compute stellar weak rates of fp- and fpg-shell nuclei. Using a proton-neutron QRPA model, the authors compute state-by-state Gamow-Teller strength distributions for parent excited states and compare the resulting electron-capture (EC) and beta-decay (BD) rates with rates obtained by replacing excited-state strength functions with ground-state (or low-lying state) distributions according to global or local BA. They report deviations that grow with temperature and density, with BD rates differing by up to three orders of magnitude, and conclude that BA is not appropriate for high-temperature, high-density conditions. For one EC set (78Ge, 67Ni, 57Ni) and one BD set (58Cr, 67Co, 70Cu), the paper reports Brink-errors and standard deviations; it also compares selected rates with IPM, IPM-03, and large-scale shell-model rates.","tokens_in":31495,"tokens_out":5022,"duration_ms":50707,"significance":"If the central claim were robustly supported, the paper would provide useful quantitative guidance for astrophysical rate tabulations and would reinforce earlier work questioning BA. The study's strength is that it performs a fully state-by-state QRPA calculation over a broad temperature-density grid and carefully distinguishes global and local BA; the standard-deviation reduction from GBA to LBA is a clear, potentially useful result. Its principal weakness is that the benchmark itself -- the excited-state GT strength distributions -- is not validated against experiment or an independent model, and the model's force constants are not given. The paper is therefore best read as a model-internal warning; the general conclusion that BA is not appropriate outruns the presented evidence.","major_comments":[{"comment":"The BA test uses the pn-QRPA excited-state GT strength distributions as the benchmark, but only ground-state strengths are compared with data (Figs. 2-3). The manuscript itself acknowledges in Section 1 and Section 2 that the model is a zero-temperature approximation with a schematic pairing-plus-quadrupole Hamiltonian and separable GT forces (Eq. 37), and the force constants χGT and κGT appearing in Eqs. (5)-(8) are never specified in the text or a table. Without an independent check of excited-state strengths (e.g., a comparison with finite-temperature QRPA or shell-model results, or a sensitivity study in χGT and κGT), the reported deviations cannot be attributed to a failure of BA rather than to the specific excited-state response of this model. This is the load-bearing point of the paper.","section":"Section 2 (Figs. 2-3, Tables 4-6) and Section 1 (Eqs. 5-8, 37)"},{"comment":"The claim that deviations \"become significant as early as neon burning\" is not supported by Table 4. For 58Cr at T = 3 GK and ρYe = 10^6 g cm^-3, ΔGBE = 0.00, and at T = 5 GK it is only 0.01; large deviations appear only at T = 15 GK and above. The three-order-of-magnitude deviation occurs only at T = 30 GK and ρYe = 10^11 g cm^-3, an extreme corner of the grid where the zero-temperature approximation is least trustworthy and where the rates are of order 10^-100 s^-1. Please revise the abstract and conclusion to state the actual onset conditions and the magnitude of deviations at physically relevant presupernova conditions.","section":"Abstract and Section 3 (Table 4)"},{"comment":"The standard deviation σGBE is dominated by the largest |Δ| values rather than by typical behavior. For 58Cr, σGBE = 330.85 is driven by Δ ≈ -1814 at the T = 30 GK, ρYe = 10^11 g cm^-3 grid point; the median and quartile deviations are far smaller. Reporting percentiles of the Brink-error distribution, or evaluating deviations only in the astrophysically relevant window (T ≲ 10 GK, ρYe ≲ 10^10 g cm^-3), would give a fairer measure of BA's validity and would not leave the impression that three-order deviations are typical.","section":"Section 2, Eq. (47)"}],"minor_comments":[{"comment":"There is a typo in the first sentence: \"the the validity\" should be \"the validity.\"","section":"Abstract"},{"comment":"The affiliations contain spacing errors: \"Pakist an\" should be \"Pakistan.\"","section":"Author affiliations"},{"comment":"The text says k is the \"total number of temperature-density grid points,\" but the grid is not listed anywhere; please specify the grid points or provide a table, so that σGBE and σLBE are reproducible.","section":"Section 2, Eq. (47)"},{"comment":"The caption says \"total GT strength (arbitrary units),\" but B(GT) is normally dimensionless; please clarify the normalization or change the caption.","section":"Table 7 caption"},{"comment":"The sentence \"The pn-QRPA approach ... was able to calculate GT strength distributions from parent excited states in the computation of stellar rates\" is awkwardly phrased; consider rephrasing for clarity.","section":"Section 2, last paragraph"},{"comment":"The figure captions and axis labels appear garbled in the manuscript text; please ensure that the published figures have legible captions and clearly labeled axes.","section":"Figures 2-5"}],"recommendation":"major_revision","confidential_remarks":"The paper's internal arithmetic appears consistent, and I found no evidence of misconduct. The main issues are external validity and overclaiming. I would be comfortable with publication if the authors either (a) supply the missing force constants, add a sensitivity analysis, and explicitly reframe the conclusion as a model-internal statement, or (b) obtain independent confirmation from finite-temperature QRPA or shell-model calculations. The 'neon burning' statement in the abstract should be corrected regardless. The comparison with LSSM and IPM-03 is useful but does not resolve the excited-state validation problem, because those models also rely on BA for higher-lying states in some cases."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is a solid model-internal test of the Brink-Axel hypothesis for fp- and fpg-shell nuclei, and within the pn-QRPA framework the BA approximation indeed looks bad. But the abstract's 'neon burning' and 'three orders of magnitude' claims outrun what the tables actually show, and the benchmark is the model's own unvalidated excited-state strengths.\n\nWhat's new: the extension of the sd-shell test [27] to the fp/fpg region, plus the explicit GBA vs LBA classification. The tables are extensive and internally consistent. The observation that LBA reduces the deviations is a useful nuance, and the comparisons with LSSM and IPM rates are a plus. The authors also state plainly that their model is a zero-temperature approximation with a schematic pairing-plus-quadrupole Hamiltonian, which is honest.\n\nSoft spots, in proportion: First, the paper claims deviations become significant at neon burning, but Table 4 for 58Cr shows Delta=0.00 at T=3 GK and rhoYe=10^6 g/cm^3, and the three-order deviation appears only at T=30 GK and rhoYe=10^11 g/cm^3, an extreme corner where the zero-temperature approximation is least justified. The stated threshold 'above 3 GK and beyond 10^6' does not match the example. Second, the benchmark lambda_QRPA uses the same pn-QRPA excited-state GT strengths that are never compared to experiment for excited states; only ground-state strengths appear in Figures 2-3. So the result is a model-internal warning, not a demonstrated failure of BA for real nuclei. Third, the GT force constants chi_GT and kappa_GT are never given numerically, so the calculation is not reproducible from the paper alone. A sensitivity analysis over these constants would be the natural fix. These are fixable; the arithmetic is not the problem.\n\nThe intended reader is someone using BA-based weak rate tables in supernova or presupernova models; the paper is a useful caution flag from one of the standard QRPA groups. It deserves a serious referee, but only with a request for major revision: tone down the abstract, add the missing parameters or a sensitivity test, and frame the conclusion as a warning within the model. I would send it to review, not desk reject.","headline":"The BA test itself holds up internally, but the neon-burning and three-order claims outrun the tables and the benchmark is model-internal.","tokens_in":32067,"tokens_out":3155,"would_cite":true,"duration_ms":30191,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["23.40.-s","26.50.+x","21.60.Jz"],"model":"deepseek-v4-flash","headline":"Reusing ground-state Gamow-Teller strength for every excited state—the Brink-Axel shortcut—misestimates stellar beta-decay rates for fp- and fpg-shell nuclei by up to three orders of magnitude.","keywords":["Brink-Axel hypothesis","proton-neutron QRPA","Gamow-Teller strength","stellar weak rates","beta decay","electron capture","core-collapse supernovae","fp- and fpg-shell nuclei"],"falsifier":"A decisive check would be to compute the same six representative nuclides (78Ge, 67Ni, 57Ni, 58Cr, 67Co, 70Cu) with an independent state-by-state nuclear model whose excited-state strengths are not taken from the pn-QRPA, on the same temperature–density grid, and compare the BA-versus-microscopic rate ratios. If an independent model finds ratios within a factor of a few rather than up to three orders of magnitude, the claimed BA failure would be a property of the QRPA excited-state strengths rather than of the nuclei.","tokens_in":30951,"feed_emoji":"⚛️","tokens_out":16811,"duration_ms":149997,"temperature":0.7,"pith_summary":"The paper tests whether a standard shortcut used in nuclear astrophysics—the Brink-Axel (BA) hypothesis, which reuses the ground-state Gamow-Teller (GT) strength distribution for every excited parent state—is safe for computing allowed weak transitions in heavy nuclei. Because stellar $\\beta$-decay and electron-capture rates set the lepton fraction, entropy, and collapsing-core mass during presupernova evolution, unreliable rates directly change the dynamics and nucleosynthesis of core-collapse supernovae. Using the proton-neutron QRPA model to compute GT strength functions state by state for eighty nuclei with proton numbers 20–34 and mass numbers 48–83, the authors compare the resulting rates with rates obtained from the global and localized BA recipes across temperatures of 1–30 GK and densities of $10$–$10^{11}$ g cm$^{-3}$. They report that BA-based $\\beta$-decay rates deviate by up to three orders of magnitude, that the deviation starts as early as neon-burning conditions, and that a localized BA version reduces but does not remove the error. The paper concludes that BA is not an appropriate assumption for reliable stellar weak rates under these conditions.","feed_headline":"Beta-decay rate estimates off by up to 1000x under Brink-Axel","feed_subtitle":"The shortcut misses state-by-state rates in the conditions where supernova simulations depend on them.","key_machinery":"The load-bearing mechanism is the proton-neutron quasiparticle random-phase approximation (pn-QRPA), a microscopic nuclear-structure method that builds correlated multi-quasiparticle states and computes Gamow-Teller strength functions from each parent state individually, including highly excited states. This state-by-state output defines the reference rates $\\lambda_{\\mathrm{QRPA}}$ against which the Brink-Axel hypothesis—the assertion that a strength function depends only on the transition energy, not on which excited state decays—is tested. The paper quantifies the comparison with the Brink errors $\\Delta^{\\mathrm{GBA}}$ and $\\Delta^{\\mathrm{LBA}}=(\\lambda_{\\mathrm{QRPA}}-\\lambda_{\\mathrm{BA}})/\\lambda_{\\mathrm{QRPA}}$ for the global and localized BA recipes, plus standard deviations over the temperature–density grid. Those error measures, together with tables of total GT strengths and centroid energies, are what show that BA systematically misplaces strength and misestimates rates.","core_discovery":"On its own terms, the paper shows that the Brink-Axel hypothesis fails quantitatively for the weak rates that matter in pre-collapse and collapsing stellar cores. For each selected nucleus, the pn-QRPA model supplies GT transition strengths from individual parent states, including states well above the ground state; replacing those excited-state strength distributions with the ground-state profile (global BA), or with a profile built from the first one or few excited states (local BA), changes the integrated $\\beta$-decay and electron-capture rates in a way that grows with temperature and density. In the $\\beta$-decay direction the BA-based rates differ from the state-by-state rates by up to three orders of magnitude, and the standard deviation of the Brink error over the grid is more than a hundred times larger for $\\beta$ decay than for electron capture. Deviations emerge at temperatures above about 3 GK and densities above about $10^6$ g cm$^{-3}$, conditions that correspond roughly to neon burning in massive stars. The authors read these results as evidence that BA-based stellar weak rates are not reliable for core-collapse supernova conditions.","pith_inferences":["The paper deliberately leaves the source of the BA failure undissected: the larger total GT strengths and misplaced centroids are tabulated, but the three-order-of-magnitude beta-decay error is never broken into a phase-space contribution and a strength-shape contribution.","An adaptive recipe is a natural next step: compute excited states explicitly below some excitation threshold and use a renormalized ground-state strength above it; the local-BA results suggest this could recover much of the accuracy at a fraction of the cost.","Because only ground-state strengths are benchmarked against experiment, a measurement of excited-state GT strength for any one representative nucleus, such as 58Cr, would be the most direct external check of whether the reported violation is real or a QRPA artifact.","The same model can be run on the full set of newly proposed presupernova nuclei to produce BA-free weak-rate tables, which would let supernova simulators quantify how much of their output depends on the approximation."],"forward_implications":["Stellar weak-rate databases for fp- and fpg-shell nuclei that rely on the Brink-Axel shortcut will need to be remade with explicit excited-state strength functions if they are to be trusted in pre-supernova and collapse simulations.","Because the deviations set in around neon-burning conditions, the error enters existing models well before the final collapse, so it can alter the lepton fraction and core entropy over an extended evolutionary phase.","The localized BA recipe reduces the standard deviation by one to two orders of magnitude but still leaves beta-decay rates off by large factors; it is a partial remedy, not a solution.","Beta decay is the direction that breaks first and hardest; any future approximation that treats beta decay and electron capture symmetrically will misrepresent the dominant error."],"supporting_citations":[{"why":"Brink's original formulation that giant-resonance strength is insensitive to the details of the initial state; this is the hypothesis under test.","marker":"[17]"},{"why":"Axel's independent statement of the same hypothesis, which together with Brink gives the BA name and the generalized GT version used here.","marker":"[18]"},{"why":"The foundational stellar weak-rate tabulations that first adopted BA for excited-state GT strengths and serve as the IPM comparison rates.","marker":"[20]"},{"why":"The companion tabulations that extended BA to GT resonances and established the recipe the paper argues against.","marker":"[10]"},{"why":"Large-scale shell-model weak-rate tables that apply BA above low-lying states and provide the state-by-state comparison for 63Ni and 63Co.","marker":"[6]"},{"why":"The recent presupernova simulation that fixes the list of rate-relevant nuclei from which the 80-nucleus sample is taken.","marker":"[12]"},{"why":"The study that introduced the energy-localized BA variant, which motivates the paper's class 2 and class 3 recipes.","marker":"[45]"},{"why":"A finite-temperature QRPA study showing that GT strengths evolve with temperature as pairing weakens; the paper cites it as prior evidence that BA is questionable.","marker":"[43]"},{"why":"The authors' earlier sd-shell study whose state-by-state BA-testing method is directly extended here to heavier nuclei.","marker":"[27]"}],"fun_headline_variants":["Brink-Axel beta-decay rates off by 1000x","Brink-Axel fails for stellar weak rates","Supernova weak rates unreliable with Brink-Axel","State-by-state rates oppose Brink-Axel shortcut"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the pn-QRPA model's Gamow-Teller strength functions from highly excited parent states are accurate enough to serve as the benchmark, yet only ground-state strengths are checked against experiment; if those excited-state strengths are wrong, the reported BA violation could be a model artifact rather than a reflection of real nuclei.","fun_headline_variants_meta":{"raw":{"variants":["Brink-Axel beta-decay rates off by 1000x","Brink-Axel fails for stellar weak rates","Supernova weak rates unreliable with Brink-Axel","State-by-state rates oppose Brink-Axel shortcut"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000305,"raw_usage":{"total_tokens":1720,"prompt_tokens":884,"completion_tokens":836,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":500,"completion_tokens_details":{"reasoning_tokens":768}},"tokens_in":500,"tokens_out":836,"duration_ms":7582,"temperature":1.0,"reasoning_tokens":768,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T16:48:30.112820+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive check would be to compute the same six representative nuclides (78Ge, 67Ni, 57Ni, 58Cr, 67Co, 70Cu) with an independent state-by-state nuclear model whose excited-state strengths are not taken from the pn-QRPA, on the same temperature–density grid, and compare the BA-versus-microscopic rate ratios. If an independent model finds ratios within a factor of a few rather than up to three orders of magnitude, the claimed BA failure would be a property of the QRPA excited-state strengths rather than of the nuclei.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Brink's original formulation that giant-resonance strength is insensitive to the details of the initial state; this is the hypothesis under test."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Axel's independent statement of the same hypothesis, which together with Brink gives the BA name and the generalized GT version used here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The foundational stellar weak-rate tabulations that first adopted BA for excited-state GT strengths and serve as the IPM comparison rates."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The companion tabulations that extended BA to GT resonances and established the recipe the paper argues against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Large-scale shell-model weak-rate tables that apply BA above low-lying states and provide the state-by-state comparison for 63Ni and 63Co."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The recent presupernova simulation that fixes the list of rate-relevant nuclei from which the 80-nucleus sample is taken."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The study that introduced the energy-localized BA variant, which motivates the paper's class 2 and class 3 recipes."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"A finite-temperature QRPA study showing that GT strengths evolve with temperature as pairing weakens; the paper cites it as prior evidence that BA is questionable."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The authors' earlier sd-shell study whose state-by-state BA-testing method is directly extended here to heavier nuclei."}],"review_version":1}