{"id":"48e406fe-2144-440d-83ea-4d3ea0d4fcea","arxiv_id":"1908.10189","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Including neutrino oscillations and self-interactions in BBN calculations sets the reheating temperature above about 1.8 MeV for radiative decays and 4 to 5 MeV for hadronic decays of long-lived massive particles.","lead":"This paper computes how neutrinos were produced and mixed when long-lived heavy particles decayed and reheated the universe to a few MeV, then uses big-bang nucleosynthesis to set new lower limits on the reheating temperature. For the first time in this setting, both neutrino oscillations and neutrino self-interactions are included, giving limits of about 1.8 MeV for radiative decays and 4 to 5 MeV for hadronic decays.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Simplified QKE damping in Eq. (2.3) is load-bearing; the un-reconciled 4.1 MeV bound of Ref. [5] suggests it may underestimate flavor equilibration and shift the radiative T_RH limit.","rationale":"The reader's weakest assumption identifies the simplified QKE damping treatment after Eq. (2.3), and that is also the most load-bearing point in my reading. The paper's central radiative bound is numerically driven by how much neutrino oscillations and self-interactions equilibrate the νe and νx spectra, and the collision term's off-diagonal structure is precisely what sets the coherence damping. Because no code or data are released and because the paper gives no quantitative estimate of the neglected 'damping-like terms', the central bound is not yet independently verifiable. The unresolved tension with the previously published 4.1 MeV BBN bound from Ref. [5] makes this more than a formal concern: either the simplified treatment underestimates the oscillation effect, or the two analyses use different BBN inputs, and the paper does not tell the reader which. For the hadronic bounds, the effect is smaller because hadronic injection dominates, but the radiative bound is a headline result and the hadronic analysis inherits the same neutrino spectra. I therefore agree with the conditional verdict: the paper should provide the numerical implementation or a full-collision-term cross-check, and it should explicitly reconcile its 1.8 MeV radiative bound with the 4.1 MeV result of Ref. [5]. No additional independent reason to reject or elevate the verdict emerged from the manuscript text; the internal equations are consistent, the parameter choices are standard, and the BBN likelihood treatment is careful about the CMB prior issue (footnote 10).","tokens_in":24971,"tokens_out":22210,"duration_ms":203934,"concrete_test":"Reproduce the 100% radiative-decay case with a full three-flavor QKE solver that includes the complete collision terms, specifically the 'damping-like terms' of Ref. [11] (or the code behind Ref. [5]/[13]), using the same oscillation parameters, T_RH grid, BBN code, and χ² prescription as this paper. Compare the resulting 95% C.L. lower bound on T_RH with 1.8 MeV and with the 4.1 MeV result of Ref. [5]. If the bound moves by more than ~0.1 MeV, or if it rises toward 4.1 MeV, the simplified collision term is load-bearing and the missing reconciliation is resolved; if it stays near 1.8 MeV, the approximation is validated and the discrepancy must instead be traced to differences in BBN inputs or likelihood.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central radiative-decay bound, T_RH ≳ 1.8 MeV (Eq. 4.10), depends on the quantum kinetic equation whose collision term is written in Eq. (2.3) as diagonal production rates plus a purely off-diagonal damping term, −D(Px x̂ + Py ŷ). The paper explicitly states that it neglects additional 'damping-like terms' from Ref. [11], following Ref. [5]. This approximation directly controls how much νe ↔ νx flavor equilibration occurs near T ∼ 1–3 MeV, where neutrino production and oscillations overlap and where the collision rate D is comparable to or larger than the vacuum oscillation frequency. If the omitted terms change the coherence evolution materially, the final νe spectrum—and therefore the weak rates (4.1)–(4.3) that set the BBN neutron-to-proton ratio—will shift, moving the T_RH bound. The paper does not estimate the size of this shift. A concrete symptom of the risk is that the new radiative bound (1.8 MeV) is more than a factor of two weaker than the three-flavor BBN bound T_RH ≳ 4.1 MeV from Ref. [5], which the introduction cites but never reconciles. Since adding neutrino self-interactions should not relax the bound, the discrepancy suggests either the simplified two-flavor/damping treatment underestimates flavor equilibration or the BBN likelihood differs; in either case the quantitative central claim is not yet secure.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper computes lower bounds on the MeV-scale reheating temperature T_RH from big-bang nucleosynthesis, including for the first time both neutrino oscillations and neutrino self-interactions in the neutrino thermalization calculation. The authors solve momentum-dependent quantum kinetic equations for an effective two-flavor system, with production and damping terms, coupled to the background evolution with energy conservation. They feed the resulting neutrino spectra into a Monte Carlo BBN calculation and combine D/H and Yp in a chi-square analysis, scanning over the baryon-to-photon ratio eta_B and T_RH. For 100% radiative decays of the massive particle X they obtain T_RH ≳ 1.8 MeV at 95% C.L.; for 100% hadronic decays they obtain T_RH ≳ 4.1–4.9 MeV for m_X = 10 GeV–100 TeV at Br = 1 and T_RH ≳ 2.1–3.7 MeV at Br = 0.001. The paper concludes that neutrino oscillations and self-interactions enhance neutrino thermalization and tighten the radiative-decay bound compared to the case without these effects.","tokens_in":25284,"tokens_out":5927,"duration_ms":58936,"significance":"If the bounds are correct, this is an important update of BBN constraints on low reheating scenarios, and the inclusion of neutrino self-interactions and oscillations in the thermalization calculation is a genuine technical step forward. The paper is careful in several respects: it evolves full momentum-dependent QKEs with energy conservation, it propagates nuclear rate uncertainties through a Monte Carlo BBN calculation, and it presents a clear chi-square combination of D/H and Yp. The hadronic-decay bounds are broadly consistent with earlier work, while the radiative-decay bound is the main new quantitative result. The central claim, however, rests on a simplified treatment of collisional damping in the QKE, and the factor-of-two discrepancy with the three-flavor BBN bound of Ref. [5] is not reconciled in the manuscript.","major_comments":[{"comment":"The collision term in Eq. (2.3) keeps only the diagonal production rates R_nu_e, R_nu_x and an off-diagonal damping term -D rho_ex, and the text immediately after Eq. (2.3) states that 'damping-like terms' from Ref. [11] are neglected following Ref. [5]. This approximation directly controls the degree of flavor coherence near T ~ 1-3 MeV, where the production and oscillation rates overlap and where the final nu_e spectrum is set. Since the weak rates in Eqs. (4.1)-(4.3) depend on f_nu_e, the approximation is load-bearing for the central radiative-decay bound of Eq. (4.10). The paper should estimate the sensitivity of T_RH,min to this approximation, for example by comparing with the full collision-term expression of Ref. [11] or by varying the functional form of the damping term.","section":"Sec. II, Eq. (2.3)"},{"comment":"The Introduction cites Ref. [5] as obtaining T_RH > 4.1 MeV (95% C.L.) from BBN with three-flavor oscillations, whereas the present analysis, which adds neutrino self-interactions and should if anything strengthen the bound, obtains T_RH ≳ 1.8 MeV. The manuscript never reconciles this factor-of-two discrepancy. Because the difference could in principle arise from the likelihood setup (free eta_B scan versus a CMB prior on eta_B, different Y_p and D/H data, or different treatment of theoretical errors) rather than from the QKE approximation, the authors should provide a direct comparison with Ref. [5] under controlled assumptions, or otherwise quantify which input drives the weaker bound.","section":"Sec. IV.B, Eq. (4.10)"}],"minor_comments":[{"comment":"There are typos in the title ('osc illating') and in the caption of Fig. 12 ('in of hadronic decaythe case'); these should be corrected.","section":"Title and Fig. 12 caption"},{"comment":"The definition of R_dist uses T_nu,eff, which depends on the neutrino number density n_nu through the same ratio; the text should state explicitly that T_nu,eff is evaluated with the same n_nu appearing in the ratio, so that the definition is unambiguous.","section":"Sec. III, Eq. (3.6)"},{"comment":"The footnote explains why the CMB prior on eta_B is not used; it would be helpful to state whether the quoted bound from Ref. [5] used such a prior, since this is directly relevant to interpreting the comparison in Eq. (4.10).","section":"Sec. IV.B, footnote 10"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a solid, workmanlike paper. The genuinely new ingredient is that neutrino self-interactions are included in the quantum kinetic treatment of neutrino thermalization for low-reheating BBN, and that changes the final abundances at the O(10)% level in the radiative case. The machinery is standard: two-flavour QKE with production and damping, energy conservation, Monte Carlo BBN with nuclear-rate uncertainties, and a chi-square combination of D/H and Yp. They are explicit about what they drop — the two-flavour approximation, the damping-like terms from McKellar-Thomson, the 7Li constraint — and those choices are defensible.\n\nThe main soft spot is quantitative: their radiative-decay bound, T_RH ≳ 1.8 MeV, is more than a factor of two weaker than the three-flavour bound of Ref [5] (4.1 MeV), and the paper never explains why. My guess is the difference is mostly statistical — they explicitly avoid using a CMB eta_B prior (footnote 10), which Ref [5] may have used — but leaving it implicit is asking for trouble. A referee should push for a direct comparison: same data, same priors, with and without self-interactions. Without that, the 'newness' of the central number is hard to calibrate.\n\nThe stress-test worry about the simplified damping term in Eq. (2.3) is legitimate but, I think, not fatal. The authors flag the approximation and follow Ref [5], and the bound moves in the direction you would expect when you add equilibration. The more serious practical issue is that no code or data are released, so the numerical results cannot be checked. That is increasingly a fair requirement for a paper whose whole point is a number.\n\nWho is this for: cosmologists and particle phenomenologists who need the current BBN-only lower bound on T_RH for moduli/gravitino scenarios. They will read this, and they will cite it alongside Ref [5]. It deserves a serious referee and, with a request for a reconciliation plot and a code release, could be accepted.","headline":"A careful, internally consistent update of low-reheating BBN constraints; the new self-interaction physics is real, but the radiative bound's discrepancy with Ref [5] is left unexplained.","tokens_in":25849,"tokens_out":3587,"would_cite":true,"duration_ms":123178,"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":"BBN alone sets a 1.8 MeV floor on the reheating temperature for radiative decays.","keywords":["reheating temperature","big-bang nucleosynthesis","neutrino thermalization","neutrino oscillations","neutrino self-interactions","hadronic decay","radiative decay","quantum kinetic equations"],"falsifier":"A direct check would be to solve the full QKE including the damping-like terms of Ref. [11] (rather than neglecting them) and recompute the neutrino spectra for $T_{\\rm RH} \\sim 1$–$5$ MeV; if the resulting $Y_p$ and D/H abundances move outside the observational 2-$\\sigma$ regions at a given $T_{\\rm RH}$, the bounds shift. Alternatively, a future precision measurement of $Y_p$ with uncertainty below 0.1% would test the predicted sharp rise in $Y_p$ at $T_{\\rm RH} \\lesssim 2$ MeV.","tokens_in":24769,"feed_emoji":"🌌","tokens_out":2750,"duration_ms":231241,"temperature":0.7,"pith_summary":"This paper asks how low the cosmic reheating temperature can be without spoiling big-bang nucleosynthesis (BBN), under the assumption that the Universe was once matter-dominated by long-lived massive particles that decay either radiatively or hadronically. The authors compute neutrino thermalization during such a low-temperature reheating, including for the first time the combined effects of neutrino oscillations and neutrino self-interactions. They find a conservative 95% C.L. lower bound of $T_{\\rm RH} \\gtrsim 1.8$ MeV for purely radiative decays, and $T_{\\rm RH} \\gtrsim 4$–5 MeV for hadronic decays, for parent masses between 10 GeV and 100 TeV. This matters because many beyond-standard-model scenarios naturally predict reheating temperatures around the MeV scale, and these bounds are the current BBN-only constraints that such models must satisfy.","feed_headline":"Neutrino mixing raises reheating floor to 1.8 MeV","feed_subtitle":"First BBN bound on MeV-scale reheating that includes neutrino self-interactions and oscillations together.","key_machinery":"The argument runs on a momentum-dependent quantum kinetic equation (QKE) for the $2\\times2$ neutrino density matrix, written in terms of polarization vectors $\\mathbf{P}$ with collision terms for production, scattering, and neutrino self-interactions. The effective two-flavor mixing scheme treats $\\nu_e$ mixing with a degenerate $\\nu_x$ state, using the solar mixing parameters $(\\delta m^2_{12}, \\theta_{12})$, and includes neutrino self-interactions through the off-diagonal matter potential and collision terms. The reheating temperature is tied to the parent decay rate via $\\Gamma_X = 3H(T_{\\rm RH})$, and the resulting neutrino spectra feed into a BBN code that tracks the neutron-to-proton ratio and light-element abundances.","core_discovery":"The central claim is that including neutrino oscillations and self-interactions in the neutrino thermalization calculation raises the minimum reheating temperature allowed by BBN, compared with earlier calculations that omitted these effects. For 100% radiative decay of the massive parent, the lower bound becomes $T_{\\rm RH} \\gtrsim 1.8$ MeV at 95% C.L.; for 100% hadronic decay, the bound rises to $T_{\\rm RH} \\gtrsim 4$–5 MeV, depending on the particle mass between 10 GeV and 100 TeV. The mechanism is that both oscillations and self-interactions equilibrate neutrino flavors, enhancing the total neutrino production rate and thereby increasing both the effective number of relativistic species and the helium-4 and deuterium abundances. A stronger neutrino bath and altered neutron-proton exchange rates push BBN further from observations, tightening the bound on $T_{\\rm RH}$.","pith_inferences":["The BBN-only bound could be sharpened by future measurements of the primordial helium abundance, since $Y_p$ responds strongly to the neutrino spectra at $T_{\\rm RH} \\sim 2$ MeV.","A full three-flavor treatment, beyond the effective two-flavor scheme, might alter the quantitative bounds by a few percent at most, given the smallness of $\\theta_{13}$ effects shown here; this could still matter for precision cosmology.","The assumption that hadrons thermalize and interact only through the listed channels leaves room for additional hadronic reactions or meson injection scenarios that could either strengthen or weaken the hadronic bound, depending on the hadron spectrum.","The same QKE machinery could be applied to scenarios with direct neutrino decay channels (e.g., $X \\to \\nu\\bar{\\nu}$), which the paper explicitly leaves out; those would produce a very different neutrino spectrum and hence different bounds."],"forward_implications":["Low-reheating models with $T_{\\rm RH}$ below about 1.8 MeV are excluded by BBN alone when the decay is purely radiative, and below about 4–5 MeV when hadronic decay channels are appreciable.","The inclusion of neutrino self-interactions is not a minor correction: it changes the bound by up to a factor of a few in the radiative case, so future low-reheating analyses should treat the full collision terms rather than only oscillations.","The bound is nearly independent of the parent particle mass for radiative decays, while hadronic-decay bounds are stronger for lighter parents because the comoving abundance of the parent scales as $T_{\\rm RH}/m_X$.","The effective two-flavor treatment with $\\theta_{12}$ suffices for BBN constraints; the reactor angle $\\theta_{13}$ has a negligible effect regardless of mass ordering.","If the same neutrino thermalization is used as an input to CMB and large-scale-structure calculations, the combined constraints on $T_{\\rm RH}$ will be at least as strong as these BBN bounds, and likely stronger."],"supporting_citations":[{"why":"Provides the first BBN-based lower bound on the reheating temperature from incomplete neutrino thermalization, establishing the method this paper extends.","marker":"[2]"},{"why":"Computes BBN constraints for hadronic decays and introduces the distortion parameter $R_{\\rm dist}$ and hadronic neutron-proton interconversion effects used here.","marker":"[3]"},{"why":"First study to include neutrino oscillations in the thermalization calculation, yielding a $T_{\\rm RH} > 2$ MeV bound that is superseded by the present work.","marker":"[4]"},{"why":"Provides the simplified treatment of collisional damping in the QKE that the paper adopts, and also reports CMB-based low-reheating constraints.","marker":"[5]"},{"why":"Derives the full QKE with damping-like terms; the present paper neglects those terms following Ref. [5], so this reference defines the approximation made.","marker":"[11]"},{"why":"Establishes the general kinetic description of mixed neutrinos, the formal basis for Eq. (2.2).","marker":"[12]"},{"why":"Provides the full collision-term expressions and computational approach for neutrino self-interactions that the paper evaluates without simplifying assumptions.","marker":"[13]"},{"why":"Gives the hadronic cross sections for neutron-proton interconversion by pions and nucleons, including the 30% error bars used for conservative bounds.","marker":"[23]"},{"why":"Supplies the hadron yields and the scaling of emitted hadron number with $m_X$, used to compute hadronic decay effects on BBN.","marker":"[25]"},{"why":"Provides the observational value of the primordial helium abundance $Y_p$ used in the $\\chi^2$ analysis.","marker":"[30]"}],"fun_headline_variants":["Neutrino oscillations and self-interactions push BBN reheating bound higher","First BBN bound includes neutrino self-interactions and oscillations: 1.8 MeV","Oscillating neutrinos set stricter BBN reheating limit","Neutrino flavor equilibrium raises BBN reheating floor to 1.8 MeV","Neutrino mixing raises BBN reheating floor to 1.8-5 MeV"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes that the collisional damping of flavor coherence is fully captured by the single $D$-term in the quantum kinetic equations, and that the additional damping-like terms from the most general description are negligibly small at the MeV temperatures where neutrino production and oscillations both operate.","fun_headline_variants_meta":{"raw":{"variants":["Neutrino oscillations and self-interactions push BBN reheating bound higher","First BBN bound includes neutrino self-interactions and oscillations: 1.8 MeV","Oscillating neutrinos set stricter BBN reheating limit","Neutrino flavor equilibrium raises BBN reheating floor to 1.8 MeV","Neutrino mixing raises BBN reheating floor to 1.8-5 MeV"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001282,"raw_usage":{"total_tokens":5217,"prompt_tokens":902,"completion_tokens":4315,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":518,"completion_tokens_details":{"reasoning_tokens":4209}},"tokens_in":518,"tokens_out":4315,"duration_ms":26502,"temperature":1.0,"reasoning_tokens":4209,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:57:17.052098+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct check would be to solve the full QKE including the damping-like terms of Ref. [11] (rather than neglecting them) and recompute the neutrino spectra for $T_{\\rm RH} \\sim 1$–$5$ MeV; if the resulting $Y_p$ and D/H abundances move outside the observational 2-$\\sigma$ regions at a given $T_{\\rm RH}$, the bounds shift. Alternatively, a future precision measurement of $Y_p$ with uncertainty below 0.1% would test the predicted sharp rise in $Y_p$ at $T_{\\rm RH} \\lesssim 2$ MeV.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the first BBN-based lower bound on the reheating temperature from incomplete neutrino thermalization, establishing the method this paper extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Computes BBN constraints for hadronic decays and introduces the distortion parameter $R_{\\rm dist}$ and hadronic neutron-proton interconversion effects used here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the simplified treatment of collisional damping in the QKE that the paper adopts, and also reports CMB-based low-reheating constraints."},{"cited_title":"Therefore, we do not use the value i n this paper","cited_arxiv_id":null,"evidence_quote":"Gives the hadronic cross sections for neutron-proton interconversion by pions and nucleons, including the 30% error bars used for conservative bounds."}],"review_version":1}