{"id":"9ee4482a-13e5-4c9e-85de-2b2f4a73626c","arxiv_id":"2501.09120","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":3.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"With current light-element abundances, hadronic decays of long-lived particles are excluded above a yield curve in lifetime, and for gravitinos this translates into reheating temperature upper limits as low as 5×10^5 GeV.","lead":"Long-lived particles that decay into quarks or gluons in the early universe can break apart the light nuclei formed during Big Bang nucleosynthesis. This paper redraws the exclusion limits using current abundance measurements and applies them to constrain how hot the universe could have been after inflation in supersymmetric gravitino models.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The headline reheating-temperature bound rests on an imported 2005 gravitino-yield fit (Eq. 7); an independent cross-check with modern production calculations is needed.","rationale":"The paper's contribution is not a new BBN code but an application of existing hadronic BBN constraints to the gravitino with updated observed abundances. The central numerical output is the reheating-temperature upper bound, and that output is obtained by inverting the BBN constraints on Y_chi using the analytic fitting formula for the thermal gravitino yield. That formula (Eq. 7) is the single point where the paper's calculation becomes dependent on an external, un-re-derived result; every TR limit in Figs. 4 and 5 scales with it. The reader's weakest-assumption analysis identified exactly this formula, and I agree. The paper itself flags the assumption of light gauginos, but does not test the sensitivity of the headline limit to the normalization or logarithmic structure of the fit. Modern calculations (Rychkov & Strumia 2007; Eberl et al. 2021) are cited but not used, so there is no evidence that the 2005 fit remains accurate at the 10-30% level that would matter for a one-significant-figure claim like TR < 5e5 GeV. A dedicated re-derivation or cross-check with the public modern result would settle the issue. The secondary inconsistency between the Section IV value (5e6 GeV for sub-TeV masses) and the conclusion (5e5 GeV) reinforces that the headline number lacks a precise provenance, but the root concern is the unvalidated yield formula. This does not warrant rejection; the qualitative BBN constraint on hadronic injection is well established and the dependence on TR is physically sound. It does warrant a conditional: the authors should validate Eq. (7) or soften the quoted bound.","tokens_in":11752,"tokens_out":10451,"duration_ms":100449,"concrete_test":"Evaluate the gravitino yield Y3/2 from Eq. (7) and from the modern thermal-production code of Eberl et al. (Ref. [49]) at TR = 1e5, 1e6, and 1e7 GeV, with all gauginos taken lighter than the gravitino. If the two yields differ by more than 30% at any of these temperatures, recompute the gluon-gluino TR limit in Fig. 4 using the modern yield; if the limit moves by more than the observational uncertainty, the headline bound TR < 5e5 GeV is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim is TR < 5e5 GeV for the gluon-gluino channel, obtained by combining the BBN limits of Figs. 3-4 with the gravitino yield formula of Eq. (7), taken verbatim from Ref. [20]. Eq. (7) is a fitting formula for thermal freeze-in production that assumes all gauginos are lighter than the gravitino; the paper itself notes that for heavy gauginos the yield must be rescaled by 3 m3/2^2 / m_gluino^2. Because the BBN limit fixes Y_chi at a given lifetime, the inferred TR is inversely proportional to the yield normalization: if the true thermal yield differs from Eq. (7), every curve in Figs. 4 and 5 shifts. The fit is from 2005 and predates the more recent calculations cited by the authors (Rychkov-Strumia 2007, Eberl et al. 2021, Refs. [48,49]); the paper does not quantify the difference. A 30% change in the yield changes TR by roughly 30% at fixed BBN limit, which could move the 5e5 GeV bound outside the stated precision. The claim also mixes scales: Section IV quotes TR <~ 5e6 GeV for sub-TeV gravitinos, while the conclusion quotes 5e5 GeV without specifying the mass at which this strongest limit occurs, making the headline number ambiguous. This is load-bearing because the paper's new contribution is precisely the numerical TR bound; if Eq. (7) is not validated, the headline result is not established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper revisits Big Bang Nucleosynthesis (BBN) constraints on hadronic energy injection from long-lived particles. Using central observed values of D/H, Yp, 3He/D, and 7Li/H, the authors draw exclusion contours in the (τχ, EvisYχ) plane for mχ = 1 TeV and Bh = 1, and then apply these contours to the gravitino by combining mass-lifetime relations, Eqs. (8)-(9), with the thermal gravitino yield formula, Eq. (7), imported from Ref. [20]. The central result is a reheating-temperature bound for the gluon-gluino channel, quoted in the conclusions as TR < 5×10^5 GeV, with a weaker bound for the photon-photino channel. The model-independent BBN analysis and the updated observational inputs are potentially useful, but the headline gravitino bound depends on a 2005 yield fit that is not re-derived or cross-checked against the modern calculations cited by the authors, and the BBN exclusion contours are based on single central values without uncertainty propagation.","tokens_in":12021,"tokens_out":6736,"duration_ms":71749,"significance":"If the results hold, the paper provides a useful update of hadronic-injection BBN constraints and sharpens the well-known gravitino bound on the reheating temperature. The explicit presentation of constraints in the (τχ, EvisYχ) plane is a strength, and the inclusion of 3He/D and 7Li/H alongside D/H and Yp adds completeness. The authors also state their assumptions clearly, e.g., mχ = 1 TeV and Bh = 1. The main limitation is that the final quantitative claim is only as robust as Eq. (7), which is adopted without re-derivation or comparison with the more recent production calculations cited in Refs. [48,49]; in addition, the absence of observational and nuclear uncertainty envelopes means that the contours and the derived TR limits are not statistical exclusion regions in the usual sense.","major_comments":[{"comment":"The central quantitative claim, TR < 5×10^5 GeV, is obtained by converting the BBN limits on Y3/2 into limits on TR using the fitting formula Eq. (7), taken verbatim from Ref. [20]. That 2005 formula assumes all supersymmetric particles are lighter than the gravitino and a particular thermal-history setup; the more recent thermal-production calculations cited by the authors, Refs. [48,49], are not used, and no estimate is given for the systematic error in Y3/2 from hadronic versus gauge contributions, heavier gauginos, or nonthermal production. At fixed lifetime the inferred TR scales inversely with the yield normalization, so even a 30% change in Y3/2 shifts the headline bound by roughly 30%. The authors should re-derive or at least cross-check Eq. (7) with modern inputs and quantify the sensitivity of the TR limits to the yield normalization.","section":"Sec. IV, Eq. (7), Figs. 4-5"},{"comment":"The exclusion contours are computed from single central values of the observed abundances and central values of nuclear cross sections; the text states that the abundances were 'estimated using the central values of cross sections and model parameters.' No 1σ or 95% intervals are propagated, and the quoted observational values are treated as hard cutoffs. Since the contours in Figs. 2 and 3 and the TR limits in Figs. 4 and 5 are derived from these lines, the sharp statement TR < 5×10^5 GeV is not robust against even small shifts in the adopted central values. Please propagate observational and nuclear uncertainties or, failing that, show how the contours move when the abundances are varied within their current uncertainties.","section":"Sec. III and Figs. 2-5"},{"comment":"The BBN constraints are computed only for mχ = 1 TeV and Bh = 1, but they are then applied to gravitino masses spanning roughly 10^2-10^5 GeV through the identification Evis = mχ. The paper does not justify that the hadronic cascade and photodissociation efficiencies are independent of the primary jet energy over this range. If the efficiencies change with Ejet, the (τχ, mχYχ) constraints are not directly transferable to all gravitino masses considered. Please quantify the mχ dependence, for example by repeating the BBN calculation at a low and a high mass, or state the mass interval over which the 1 TeV approximation is valid.","section":"Secs. III-IV, Figs. 3-5"},{"comment":"The conclusion states 'we find TR < 5×10^5 GeV' for the gluon-gluino channel, but Sec. IV gives a more nuanced result: for sub-TeV gravitinos the 3He bound gives TR <~ 5×10^6 GeV, while for masses between the weak scale and 7 TeV the D/H bound gives limits varying from 5×10^5 GeV to 2×10^7 GeV. As written, the headline number is ambiguous and appears to be a global bound when it is actually a mass-dependent limit. The conclusions should specify the mass interval on which the 5×10^5 GeV value applies.","section":"Sec. IV vs. Sec. V"},{"comment":"The 'updated' light-element abundances are attributed to Ref. [33], which is the authors' own earlier publication rather than an independent observational compilation. The numerical values used in Figs. 2-5, such as 3He/D = 0.4352, are not traced to primary measurements or to standard compilations such as Refs. [4,5,7,34]. Since the updated observational inputs are the main driver of the claimed improvement, the paper should cite the primary abundance determinations and state any corrections or re-analyses taken from Ref. [33].","section":"Sec. II, Refs. [33]"}],"minor_comments":[{"comment":"Ref. [35] duplicates Ref. [20] (the same paper by Kawasaki, Kohri, and Moroi is cited with different formatting); please merge them.","section":"References"},{"comment":"The text defines Yi = ni/nb for light elements but Yχ = nχ/s for the decaying particle; the two normalizations should be reconciled or their relationship clarified, since the conversion between nb and s enters the comparison of BBN abundances with Yχ.","section":"Sec. II, Eq. (1) and Eq. (4)"},{"comment":"The legend shows '3He/D' twice and it is not clear what the dashed versus solid purple lines represent; please clarify whether the two curves correspond to different abundance values or different epochs.","section":"Fig. 3"},{"comment":"The subsection heading 'Hadronic Injection' appears without a number, interrupting the numbering of Sections III.A and the following material; please fix the section numbering.","section":"Sec. III"},{"comment":"There are several small typos, including the duplicated comma after 'The gravitino is the supersymmetric partner of the graviton [23],,' and 'dark mater relics' in the conclusions; a careful proofreading pass is needed.","section":"Sec. I and Sec. V"}],"recommendation":"major_revision","confidential_remarks":"The paper is a useful update but the headline TR bound needs to be decoupled from the unvalidated 2005 yield formula. I would ask the authors to show the BBN-only contours as the primary model-independent result, and then present the gravitino conversion as a specific application with a clear sensitivity estimate. It would also help to state explicitly how the present constraints improve on Refs. [11,12] beyond the updated central abundances, given that Ref. [33] is a previous paper by the same group. These issues are fixable within the manuscript's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The useful part of this paper is straightforward: it recomputes hadronic BBN constraints on long-lived particles using current central values for D/H, Yp, 3He/D, and 7Li/H, and then translates those into upper bounds on the gravitino reheating temperature. The qualitative physics is sound and matches the established picture — early decays raise Yp, later decays overproduce D and 3He, and 7Li is destroyed. The figures are clear, and the authors correctly flag that the 7Li bounds should be treated cautiously. As an update of the Kawasaki-Kohri-Moroi framework, it does what it says.\n\nThe real soft spots are in the quantitative packaging. First, the bounds are computed from single central values of the observed abundances with no error propagation, and no code or data files are provided. That makes the exclusion contours hard to audit, and for 3He/D in particular the observational systematics are not negligible. Second, the headline result — TR < 5 × 10^5 GeV for the gluon-gluino channel — is ambiguous. Section IV gives a mass-dependent statement: the 5 × 10^5 GeV bound applies to gravitino masses between the weak scale and about 7 TeV, while sub-TeV masses give TR below roughly 5 × 10^6 GeV. The conclusion drops that nuance and states only the strongest limit. That should be fixed.\n\nThird, and most important, the paper imports Eq. (7), the 2005 fitting formula for thermal gravitino production, without cross-checking it against the more recent calculations it cites (Rychkov-Strumia, Eberl et al.). The stress-test note is right that a 30% change in the yield shifts TR by about 30% at fixed BBN limit. The authors do mention the rescaling for heavy gauginos, but they do not quantify how much their limits would move under modern production computations. Since the entire gravitino application hangs on that formula, this is a load-bearing caveat, not a corner case.\n\nThat said, the qualitative conclusions are unlikely to be wrong: hadronic gravitino decays still force reheating temperatures well below 10^6–10^7 GeV for sub-TeV gravitinos, consistent with a large existing literature. The paper is a legitimate, if incremental, contribution. It deserves a serious referee, but it needs revision before acceptance: add uncertainty bands or at least justify their neglect, clarify the mass dependence of the headline bound, and either validate Eq. (7) against modern thermal production calculations or state clearly how the limits would shift.","headline":"Useful refresh of hadronic BBN bounds with updated central abundances, but the gravitino headline is oversimplified and the analysis lacks uncertainty propagation.","tokens_in":12604,"tokens_out":1761,"would_cite":false,"duration_ms":20651,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["98.80.Cq","14.80.Ly"],"model":"deepseek-v4-flash","headline":"Updated abundances of deuterium, helium-3, helium-4, and lithium-7 tighten Big Bang Nucleosynthesis bounds on hadronic injection and, for a gluon-gluino gravitino, cap the reheating temperature at about 5×10^5 GeV.","keywords":["big bang nucleosynthesis","hadronic injection","gravitino problem","reheating temperature","long-lived particles","light-element abundances","photodissociation","hadrodissociation"],"falsifier":"Measure the primordial ${}^3\\mathrm{He}$ abundance independently of the central value adopted here, since it drives the low-mass bound; a value a factor of two lower would weaken the $T_R \\lesssim 5\\times10^5$ GeV constraint, while a value a factor of two higher would strengthen it. Alternatively, determine from colliders whether the gluino is heavier than the gravitino, which rescales the gravitino yield by a known factor and would shift all the quoted $T_R$ curves.","tokens_in":11520,"feed_emoji":"🌌","tokens_out":10229,"duration_ms":87270,"temperature":0.7,"pith_summary":"This paper argues that the newest measured abundances of the light elements—deuterium, helium-4, helium-3, and lithium-7—put stronger Big Bang Nucleosynthesis limits on long-lived particles that inject hadronic energy into the primordial plasma. Using those updated abundances, it maps model-independent exclusion contours in the plane of particle lifetime and energy injection, and converts them into an upper bound on the reheating temperature for the gravitino, a classic long-lived supersymmetric particle. For a gravitino decaying into a gluon and a gluino, the updated bounds require a reheating temperature below about $5\\times10^5$ GeV for gravitino masses below the weak scale, with the tightest constraint coming from helium-3 overproduction. Since the reheating temperature controls the thermal production of dark matter relics, these bounds constrain the early thermal history of the universe.","feed_headline":"Big Bang nucleosynthesis caps reheat temperature at 5×10^5 GeV","feed_subtitle":"Updated deuterium, helium-3, and lithium-7 abundance data tighten limits on late particle decays and the early thermal history.","key_machinery":"The analysis is carried out in the $(\\tau_\\chi,\\, E_{\\rm vis}Y_\\chi)$ plane, where $E_{\\rm vis}Y_\\chi = m_\\chi Y_\\chi$ is the visible energy injected by a particle of mass $m_\\chi$ and yield $Y_\\chi \\equiv n_\\chi/s$ at times $t \\ll \\tau_\\chi$; the paper fixes $m_\\chi = 1$ TeV and hadronic branching ratio $B_h = 1$. The physics engine is the hadronic shower: at $t \\lesssim 100$ s injected nucleons and mesons interconvert protons and neutrons, raising the $n/p$ ratio and hence ${}^4\\mathrm{He}$ and D production, while at $t \\gtrsim 100$ s high-energy nucleons destroy ${}^4\\mathrm{He}$ and produce D, ${}^3\\mathrm{He}$, ${}^6\\mathrm{Li}$, and ${}^7\\mathrm{Li}$. For the gravitino application, the key identity is the thermal yield-to-reheating relation $Y_{3/2} \\simeq 1.9\\times10^{-12}\\,(T_R/10^{10}\\,\\mathrm{GeV})\\,(1+0.045\\ln(T_R/10^{10}\\,\\mathrm{GeV}))\\,(1-0.28\\ln(T_R/10^{10}\\,\\mathrm{GeV}))$ imported from Ref. [20], combined with the lifetime scalings $\\tau_{3/2} \\propto (m_{3/2}/100\\,\\mathrm{GeV})^{-3}$ for the gluon-gluino and photon-photino channels, which converts the BBN exclusion curves into upper limits on $T_R$ as a function of $m_{3/2}$.","core_discovery":"The paper's central claim is that present-day measurements of D/H, $Y_p$, ${}^3\\mathrm{He}/\\mathrm{D}$, and ${}^7\\mathrm{Li}/\\mathrm{H}$, combined with a hadronic-decay treatment of long-lived particles, yield updated BBN constraints on the hadronic injection parameter $E_{\\rm vis}Y_\\chi$ across lifetimes $10^{-2}$ to $10^{12}$ s, and that these constraints translate directly into an upper bound on the reheating temperature in the gravitino scenario. The dominant constraints come from overproduction of ${}^3\\mathrm{He}$ and D, with ${}^4\\mathrm{He}$ ($Y_p$) and ${}^7\\mathrm{Li}$ playing secondary roles. Adopting the standard thermal yield $Y_{3/2}(T_R)$ and the gluon-gluino decay channel ($B_h=1$), the paper finds $T_R < 5\\times10^5$ GeV for gravitino masses below the weak scale, relaxing to about $10^{10}$ GeV for masses up to $\\sim100$ TeV; the photon-photino channel ($B_h \\approx 10^{-3}$) yields weaker bounds, with $T_R$ ranging from $5\\times10^5$ GeV to $7\\times10^8$ GeV below $\\sim2$ TeV. The authors state that BBN can thereby be used as a probe of the reheating temperature.","pith_inferences":["The constraint map is model-independent in ($\\tau_\\chi$, $E_{\\rm vis}Y_\\chi$), so the same updated contours would apply to other long-lived hadronically decaying species such as moduli, saxions, or dark-sector particles; the gravitino is the application the authors choose, and this extrapolation goes beyond their explicit scope.","If the gravitino yield were instead dominated by nonthermal production, for example from inflaton decay, the $T_R$ bounds derived here would not directly apply; the same BBN contours would instead bound the product of branching ratio and yield, so the quoted $T_R$ limit is contingent on the assumed thermal freeze-in efficiency.","A positive measurement of the gluino mass heavier than the gravitino would rescale the gravitino yield by a factor of order $3m_{3/2}^2/m_{\\tilde g}^2$, shifting all quoted $T_R$ limits and giving a collider-tested handle on these cosmological bounds.","The ${}^3\\mathrm{He}$ channel provides the strongest constraints at low gravitino mass, so a future measurement of primordial helium-3 (currently the least well-determined abundance used here) would most directly tighten or loosen the headline $T_R < 5\\times10^5$ GeV result."],"forward_implications":["A gravitino with hadronic decays and mass below the weak scale forces $T_R \\lesssim 5\\times10^5$ GeV in the gluon-gluino channel, which is below the typical temperature needed for thermal weakly interacting dark matter production, thereby constraining dark matter models.","The reheating temperature bounds are order-of-magnitude stronger for the gluon-gluino channel than for the photon-photino channel, so the gravitino's decay mode matters as much as its mass for cosmological viability.","Hadronic injection from long-lived particles can destroy ${}^7\\mathrm{Be}$ and reduce ${}^7\\mathrm{Li}$, and the paper's ${}^7\\mathrm{Li}/\\mathrm{H}$ constraints bracket the parameter region that could alleviate the cosmological lithium problem.","For gravitino masses between the weak scale and about 7 TeV, deuterium overproduction gives the strongest limits, with $T_R$ varying from $5\\times10^5$ GeV to $2\\times10^7$ GeV; above that, proton-neutron interconversion dominates and the bound relaxes toward $10^{10}$ GeV.","The bounds are presented as model-independent in the ($\\tau_\\chi$, $E_{\\rm vis}Y_\\chi$) plane, so the same updated contours constrain any long-lived particle that decays hadronically, not only the gravitino."],"supporting_citations":[{"why":"Supplies the hadronic-decay BBN formalism, the thermal gravitino yield formula of Eq. (7), and the gravitino lifetime scalings used to convert BBN bounds into $T_R$ limits.","marker":"[20]"},{"why":"Establishes the method for hadronic decay of late-decaying particles and the BBN constraints that this paper updates with new abundance data.","marker":"[11]"},{"why":"Source of the updated light-element abundances (D/H, $Y_p$, ${}^3\\mathrm{He}/\\mathrm{D}$, ${}^7\\mathrm{Li}/\\mathrm{H}$) that define the new exclusion contours.","marker":"[33]"},{"why":"Previous revisiting of BBN constraints on long-lived decaying particles that this work refines with updated abundances.","marker":"[12]"},{"why":"The prior BBN bound on the reheating temperature from an unstable gravitino that this paper updates.","marker":"[27]"},{"why":"Particle-data compilation from which the adopted observational light-element abundance values are drawn.","marker":"[34]"}],"fun_headline_variants":["BBN sets reheat cap at 5×10^5 GeV for gravitinos","Hadronic BBN bounds squeeze gravitino reheat temperature","Updated BBN bounds tighten the gravitino reheat limit","Gravitino problem: BBN puts reheat temperature below 5×10^5 GeV","Reheat temperature capped by BBN hadronic bounds"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire reheating-temperature bound rests on the assumed thermal gravitino yield formula, which presumes gravitinos are produced thermally during reheating with all supersymmetric gauginos lighter than the gravitino; if gravitino production is non-thermal or the sparticle spectrum differs, every quoted $T_R$ limit shifts.","fun_headline_variants_meta":{"raw":{"variants":["BBN sets reheat cap at 5×10^5 GeV for gravitinos","Hadronic BBN bounds squeeze gravitino reheat temperature","Updated BBN bounds tighten the gravitino reheat limit","Gravitino problem: BBN puts reheat temperature below 5×10^5 GeV","Reheat temperature capped by BBN hadronic bounds"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000812,"raw_usage":{"total_tokens":3582,"prompt_tokens":987,"completion_tokens":2595,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":603,"completion_tokens_details":{"reasoning_tokens":2497}},"tokens_in":603,"tokens_out":2595,"duration_ms":17203,"temperature":1.0,"reasoning_tokens":2497,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:10:37.778509+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the primordial ${}^3\\mathrm{He}$ abundance independently of the central value adopted here, since it drives the low-mass bound; a value a factor of two lower would weaken the $T_R \\lesssim 5\\times10^5$ GeV constraint, while a value a factor of two higher would strengthen it. Alternatively, determine from colliders whether the gluino is heavier than the gravitino, which rescales the gravitino yield by a known factor and would shift all the quoted $T_R$ curves.","supporting_citations":[{"cited_title":"Kawasaki, K","cited_arxiv_id":null,"evidence_quote":"Supplies the hadronic-decay BBN formalism, the thermal gravitino yield formula of Eq. (7), and the gravitino lifetime scalings used to convert BBN bounds into $T_R$ limits."},{"cited_title":"Updated Big Bang Nucleosynthesis Bounds on Long-lived Particles from Dark Sectors","cited_arxiv_id":"2311.07688","evidence_quote":"Source of the updated light-element abundances (D/H, $Y_p$, ${}^3\\mathrm{He}/\\mathrm{D}$, ${}^7\\mathrm{Li}/\\mathrm{H}$) that define the new exclusion contours."},{"cited_title":"Kohri, T","cited_arxiv_id":null,"evidence_quote":"The prior BBN bound on the reheating temperature from an unstable gravitino that this paper updates."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Particle-data compilation from which the adopted observational light-element abundance values are drawn."}],"review_version":1}