{"id":"1bb514d9-af9b-4f09-a802-dd6564d63cb6","arxiv_id":"2412.18711","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Using the complete 3800-day RENO sample, the collaboration measures sin^2(2theta13)=0.0920+/-0.0060 and |Delta m^2_ee|=(2.57+/-0.12)x10^-3 eV^2, the final precision of the experiment.","lead":"The RENO experiment's final 3,800-day dataset yields its most precise reactor antineutrino oscillation parameters: sin^2(2theta13)=0.0920 and |Delta m^2_ee|=2.57x10^-3 eV^2. A generalist might read this because the precision of theta13 directly anchors upcoming searches for CP violation in neutrinos.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 5 MeV spectral excess is acknowledged but never stress-tested; a fit excluding 3.6–6.6 MeV should be reported to show it cannot bias the central oscillation parameters.","rationale":"The reader's verdict is CONDITIONAL, and my concrete concern reinforces the need for caution without changing that verdict. The paper's central physics claim is internally consistent: the candidate counts, background-subtracted rates, and chi^2/NDF = 64/66 are coherent, and the far-to-near ratio is a standard and generally reliable technique. I do not see a derivation gap in the quoted central values. My load-bearing concern is narrower: the paper explicitly acknowledges the 5 MeV excess as an unexplained deviation from the reactor model, yet the oscillation fit spans the affected energy range and no robustness test is reported. The far-to-near geometry does not guarantee that a common reactor-model error cancels, because the near and far detectors sample different reactor mixtures (Table VII), so an energy-dependent discrepancy can leave a residual that mimics or distorts the oscillation signal. This is a testable concern, and the test I propose would settle it. I also note the internal inconsistency in Section XI (precision improved to 6.5% versus 6.4% elsewhere) and the absence of released data or code, but these do not change the verdict beyond what the reader already recorded. My agreement with the reader is partial: we both flag a systematic-assumption risk, but the reader emphasizes detector-identity and background cancellation, while I focus on the explicitly acknowledged 5 MeV model discrepancy, which is the most direct threat to the fitted spectrum.","tokens_in":15718,"tokens_out":11042,"duration_ms":114050,"concrete_test":"Re-run the Section XI chi^2 fit with all prompt-energy bins in 3.6 < Ep < 6.6 MeV excluded, or equivalently add a free 5 MeV Gaussian-shaped nuisance term, and compare the best-fit sin^2(2theta13), |Delta m^2_ee|, and their quoted errors. If excluding that region shifts either parameter by more than about 0.5 sigma of the quoted total uncertainty, the 5 MeV model discrepancy is biasing the central result and the stated precision is not model-robust. If the shift is negligible, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central fit (Section XI, Eq. 2) determines sin^2(2theta13) and |Delta m^2_ee| from the far-to-near ratio over 1.2–8.0 MeV using a reactor antineutrino model. Figure 7 shows a large, unmodeled spectral deviation near 5 MeV in both detectors, which the paper attributes to a flaw in the reactor model but does not otherwise constrain or remove. A purely common multiplicative flux error would cancel in a far-to-near ratio, but cancellation is not automatic: the near and far detectors weight the six reactors very differently (Table VII: Reactor 3 contributes 32.4% near versus 15.3% far; Reactors 5+6 contribute 21.9% near versus 39.3% far), so any energy-dependent, reactor-mix-dependent component of the 5 MeV excess will not cancel. The chi^2 as written has no nuisance parameter for this bump, and the paper reports no stability test with the 3.6–6.6 MeV region removed. If a residual 5 MeV discrepancy enters the fit, the oscillation amplitude and frequency could be pulled beyond the quoted 0.0044/0.0041 (sin^2(2theta13)) and 0.10/0.05 (|Delta m^2_ee|) errors, because those uncertainties do not include a model-validity term for the acknowledged spectral failure. This is the most load-bearing point: both central parameters are extracted from the energy spectrum in the region where the reactor model is explicitly shown to disagree with data.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript reports the final RENO measurement of reactor antineutrino oscillation parameters using the complete data sample of 3800 days, with 1,211,995 (144,667) IBD candidates in the near (far) detector. The far-to-near spectral ratio is fit with a pull-parameter chi-square that accounts for background, reactor flux, efficiency, and energy-scale systematics. The best-fit values are sin^2(2theta_13) = 0.0920 (+0.0044/-0.0042 stat.) (+0.0041/-0.0041 syst.) and |Delta m^2_ee| = (2.57 (+0.10/-0.11 stat.) (+0.05/-0.05 syst.)) x 10^-3 eV^2, with chi^2/NDF = 64.0/66. The paper also reports improved estimates of backgrounds, a reduced reactor-flux uncertainty, and a merged 9Li/8He spectrum from both detectors. The main caveat is the acknowledged 5 MeV spectral excess, which is shown in Fig. 7 but not stress-tested in the central fit.","tokens_in":15895,"tokens_out":4476,"duration_ms":45914,"significance":"If the result is robust, this is a valuable final measurement from RENO, improving the precision of sin^2(2theta_13) and |Delta m^2_ee| beyond the previous 2200-day result and reaching the experiment's statistical goal. The paper benefits from the full 3800-day dataset, explicit background estimates, a far-to-near ratio analysis that cancels many correlated systematic uncertainties, and a toy-MC treatment of reactor correlations. The central values are consistent with previous RENO results and with the global reactor picture. However, because the central fit includes the spectral region where the reactor antineutrino model is explicitly shown to disagree with data, the claimed systematic uncertainties do not yet include a model-validity term for that discrepancy. The manuscript would be strengthened by a robustness test excluding the 3.6-6.6 MeV region or by an additional nuisance parameter.","major_comments":[{"comment":"The chi-square fit of Eq. (2) uses the full 1.2-8.0 MeV range, while Fig. 7 demonstrates a large unmodeled spectral deviation near 5 MeV in both detectors. Because the near and far detectors weight the six reactors very differently (Reactor 3 is 32.4% near vs 15.3% far; Reactors 5+6 are 21.9% near vs 39.3% far, Table VII), any reactor-mix-dependent component of the 5 MeV excess will not automatically cancel in the far-to-near ratio. The paper reports no fit excluding 3.6-6.6 MeV and no nuisance parameter for the bump. I request a stability test with the 3.6-6.6 MeV region removed, or an explicit additional systematic term; without this, the quoted uncertainties omit the model-validity question in the very region where the model is known to fail.","section":"Section XI, Eq. (2); Fig. 7; Table VII"},{"comment":"The 9Li/8He background is the largest background rate uncertainty, and the analysis merges the near and far detector spectra into a single shape, assuming identical detector responses. The correlated shape uncertainties then cancel in the far-to-near ratio by construction. The paper should show the far-only background fit or quantify how a detector-dependent energy-scale difference in the merged spectrum changes the extrapolation from the 7-12 MeV normalization region to 1.2-8 MeV. Otherwise the cancellation of the 9Li/8He shape uncertainty is an assumption rather than a demonstrated behavior.","section":"Section VIII, Fig. 4"},{"comment":"The definition of T_i^{F/N} in Eq. (2) references the reactor antineutrino model, but Fig. 8 and Fig. 10 construct the no-oscillation expectation from the measured near-detector spectrum. The text should state explicitly whether the model-based prediction or the near-spectrum-based prediction is used in the chi-square fit; if both are used in different places, the difference in their treatment of the 5 MeV excess should be quantified.","section":"Section XI, Eq. (2)"}],"minor_comments":[{"comment":"The abstract and Section XII refer to '3800 live days', but Table VI reports 3307.25 live days for the near detector and 3737.85 for the far detector. Please clarify which quantity the '3800 live days' label refers to and whether it is the far-detector live time or the total elapsed time.","section":"Abstract and Table VI"},{"comment":"The text states that the precision of sin^2(2theta_13) improves from 7.5% to 6.5%, while the abstract and Table XI state 6.4%. These numbers should be made consistent.","section":"Section XI, closing paragraph"},{"comment":"The shaded band is described as the 'uncertainty of the expected spectrum' but the text does not specify which sources of uncertainty are included. Please define the band (e.g., reactor flux, energy scale, backgrounds, or a combination) so that the reader can judge the significance of the 5 MeV deviation.","section":"Fig. 7"},{"comment":"The number of energy bins used in the fit, the bin width, and the exact treatment of the two data-taking periods before and after the 252Cf contamination should be given explicitly alongside Eq. (2).","section":"Section XI"},{"comment":"The text says that, because of the 5 MeV excess, the expected L_eff/E_nu distribution is taken from the measured near-detector spectrum instead of the IBD Monte Carlo spectrum. It should be clarified that this choice is used for the display of the survival probability, not for the central chi-square fit of Eq. (2).","section":"Fig. 10"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a standard final-result paper from RENO and the measurement is internally consistent. The main issue is the lack of a robustness test for the 5 MeV spectral excess in the central oscillation fit; adding such a test is feasible within the manuscript's scope and should be required before publication. The 9Li/8He merging assumption also deserves a brief sensitivity check. I do not see a reason to question the data integrity or the general analysis framework."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The headline is simple: this is RENO's final word on theta13 and Delta m^2_ee, and it is a competent, internally consistent measurement. The new content is the complete 3800-day dataset, a toy-MC treatment of reactor-flux correlations that shrinks an old systematic, and a combined near/far 9Li/8He background spectrum. The precision improvement from 7.5% to 6.4% on sin^2(2theta13) and 5.2% to 4.5% on |Delta m^2_ee| is real and roughly what you would expect from the added statistics and the reduced systematics. The central values agree with the previous RENO results and with Daya Bay, so nothing here is going to shake the three-flavor picture.\n\nThe one soft spot worth a referee's attention is the 5 MeV excess. The paper shows it clearly in both detectors and attributes it to the reactor model, but then proceeds to fit the full 1.2–8.0 MeV region without reporting a stability test that excludes 3.6–6.6 MeV. The stress-test note worries that the near and far detectors weight the six reactors differently, so a reactor-mix-dependent component of the excess might not cancel in the far-to-near ratio. That is a fair worry in principle. In practice, the fit quality is good (chi2/NDF = 64/66) and the excess appears to scale with total IBD rate, which suggests it mostly cancels in the ratio. But the paper does not demonstrate that explicitly, and it would cost little to add a fit with the bump region removed. That is a minor revision, not a fatal flaw.\n\nTwo smaller issues. The text says 6.5% in one sentence in Section XI and 6.4% in the abstract and Table XI; obviously a typo, but it should be fixed. And the paper, like most collaboration papers, does not release data or code, so independent reanalysis is impossible; that is standard but limits how much we can verify beyond the reported numbers.\n\nOverall, the measurement itself is sound. The far-to-near ratio analysis is well executed, the systematic tables are detailed, and the background estimates are careful. The reader's conditional verdict is about right. The paper deserves serious peer review; the main request should be a robustness check for the 5 MeV region and a cleanup of the precision inconsistency. I would take it to a reading group to discuss the toy-MC correlation treatment, and I would cite it as the final RENO result.","headline":"The final RENO 3800-day result is a solid, mature measurement: no surprises, modest precision gain, and one legitimate robustness request about the 5 MeV excess that the paper should address before publication.","tokens_in":16753,"tokens_out":2375,"would_cite":true,"duration_ms":25591,"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":"With 3,800 live days of data, RENO's final analysis fixes sin^2(2theta13)=0.0920 and |Delta m^2_ee|=(2.57+/-0.12)x10^-3 eV^2.","keywords":["reactor antineutrino","neutrino oscillation","theta13","inverse beta decay","far-to-near ratio","mass-squared difference","RENO experiment"],"falsifier":"Re-fit the far-detector spectrum alone, without using the near detector as a normalization, modeling the reactor flux and detector response absolutely; if the resulting $\\sin^{2}2\\theta_{13}$ or $|\\Delta m^{2}_{ee}|$ differs from the ratio-based result by more than the combined quoted uncertainty, the cancellation assumption is falsified. Alternatively, an analysis that removes the 9Li/8He background constraint or replaces the merged background spectrum with independently measured near and far spectra would reveal whether the background treatment biases the fitted parameters.","tokens_in":15399,"feed_emoji":"⚛️","tokens_out":6445,"duration_ms":57354,"temperature":0.7,"pith_summary":"The paper reports the final oscillation measurement from the RENO reactor neutrino experiment, using the complete 3,800-day data sample collected before the detectors were switched off. The central claim is that the full dataset determines the neutrino mixing angle $\\sin^{2}2\\theta_{13}=0.0920$ with a total precision of 6.4% and the effective mass-squared difference $|\\Delta m^{2}_{ee}|=(2.57\\pm0.12)\\times10^{-3}~\\mathrm{eV}^{2}$ with 4.5% precision. This matters because $\\theta_{13}$ is the smallest neutrino mixing angle and a key input for future searches for CP violation in long-baseline experiments. The analysis reaches the statistical precision the experiment was designed for, and the improvement comes from a 40% increase in statistics plus a reduced reactor-flux uncertainty.","feed_headline":"RENO's final run pins theta13 to 6.4 percent","feed_subtitle":"Complete 3,800-day dataset also fixes the neutrino mass splitting to 4.5 percent precision.","key_machinery":"The load-bearing observable is the spectral ratio of far-detector to near-detector inverse-$\\beta$-decay candidates, binned in prompt energy. Because the two detectors are nearly identical, all correlated uncertainties, including detection efficiency, reactor flux normalization, and energy-scale correlations, cancel in this ratio, leaving only uncorrelated uncertainties plus background-rate and pull parameters. A $\\chi^{2}$ function with pull terms for reactor flux, detection efficiency, energy scale, and backgrounds is minimized over the oscillation parameters, and the 9Li/8He cosmogenic background spectrum is updated with a merged near-plus-far sample to reduce shape errors.","core_discovery":"Using 1,211,995 inverse-$\\beta$-decay candidates in the near detector and 144,667 in the far detector, the collaboration observes a clear energy- and baseline-dependent disappearance of reactor antineutrinos. A fit to the far-to-near ratio of prompt spectra yields $\\sin^{2}2\\theta_{13}=0.0920^{+0.0044}_{-0.0042}(\\mathrm{stat.})^{+0.0041}_{-0.0041}(\\mathrm{syst.})$ and $|\\Delta m^{2}_{ee}|=\\left[2.57^{+0.10}_{-0.11}(\\mathrm{stat.})^{+0.05}_{-0.05}(\\mathrm{syst.})\\right]\\times10^{-3}~\\mathrm{eV}^{2}$. This is the experiment's final word on neutrino oscillation amplitude and frequency, improving precision from 7.5% to 6.4% on the mixing angle and from 5.2% to 4.5% on the mass splitting compared with the 2,200-day result. The paper also reports that the statistical error has reached the experiment's design goal, so additional running would barely improve the measurement.","pith_inferences":["Going beyond the paper: a comparison of this result with the other leading reactor experiment's $\\theta_{13}$ measurement would test the consistency of the three-flavor oscillation framework at the few-percent level, since the two experiments use different baselines, detector technologies, and reactor cores.","Going beyond the paper: if the far-to-near cancellation is as complete as claimed, the same 3,800-day dataset could be re-analyzed with an absolute-rate fit to extract the reactor antineutrino flux normalization, providing a cross-check of the reactor flux anomaly.","Going beyond the paper: the energy-scale uncertainty dominates the $|\\Delta m^{2}_{ee}|$ error, so a future analysis using a precisely calibrated monoenergetic source or a different neutron-capture channel could reduce the dominant systematic and sharpen the mass-splitting measurement.","Going beyond the paper: the 5 MeV excess, if due to a reactor spectrum model defect, may affect the oscillation fit at low energies; testing the fit with the 5 MeV region masked would quantify its influence."],"forward_implications":["If correct, $\\theta_{13}$ is known to 6.4% total precision and $|\\Delta m^{2}_{ee}|$ to 4.5%, making the reactor measurement the most precise determination of the smallest mixing angle.","The reduced reactor-flux uncertainty, from 0.9% to 0.04–0.05% uncorrelated, means future reactor experiments can rely on the far-to-near technique to suppress flux systematics.","The persistent spectral excess near 5 MeV, correlated with reactor thermal power, indicates that current reactor antineutrino spectrum models need revision.","Statistical precision is saturated; any further improvement in $\\sin^{2}2\\theta_{13}$ and $|\\Delta m^{2}_{ee}|$ from this experiment would require better control of detector energy scale and background systematics."],"supporting_citations":[{"why":"Previous RENO 2,200-day result that this measurement extends and compares its precision against.","marker":"[1]"},{"why":"Earlier RENO 500-day measurement establishing the measurement chain and baseline for oscillation parameters.","marker":"[10]"},{"why":"Defines the effective mass-squared difference $\\Delta m^{2}_{ee}$ used in the survival probability formula.","marker":"[12]"},{"why":"Documents the detector simulation, energy calibration, reconstruction, and 9Li/8He spectrum methods used here.","marker":"[15]"},{"why":"Reference reactor antineutrino spectrum model used for the expected rates and spectral shape.","marker":"[26]"},{"why":"Companion reactor antineutrino spectrum model used in constructing the expected far-to-near ratio.","marker":"[27]"},{"why":"Source of the $\\chi^{2}$ pull-term fitting method for propagating systematic uncertainties into oscillation parameters.","marker":"[28]"},{"why":"Measures cosmogenic 9Li and 8He production rates that underpin the updated background spectrum.","marker":"[34]"}],"fun_headline_variants":["RENO's 3800-day run pins theta13 to 6.4% precision","RENO complete dataset: theta13 at 6.4%, mass splitting at 4.5%","After 3800 days, RENO achieves precise theta13 and mass splitting","RENO's final measurement sharpens theta13 and mass splitting","RENO's complete sample yields precise oscillation amplitude and frequency"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The analysis assumes the near and far detectors behave identically within the quoted uncorrelated uncertainties, so every correlated detector, reactor, and efficiency effect cancels in the far-to-near ratio; if a difference between the detectors is missed, both fitted oscillation parameters would shift.","fun_headline_variants_meta":{"raw":{"variants":["RENO's 3800-day run pins theta13 to 6.4% precision","RENO complete dataset: theta13 at 6.4%, mass splitting at 4.5%","After 3800 days, RENO achieves precise theta13 and mass splitting","RENO's final measurement sharpens theta13 and mass splitting","RENO's complete sample yields precise oscillation amplitude and frequency"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001661,"raw_usage":{"total_tokens":6691,"prompt_tokens":1141,"completion_tokens":5550,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":757,"completion_tokens_details":{"reasoning_tokens":5442}},"tokens_in":757,"tokens_out":5550,"duration_ms":35273,"temperature":1.0,"reasoning_tokens":5442,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T04:33:25.194799+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-fit the far-detector spectrum alone, without using the near detector as a normalization, modeling the reactor flux and detector response absolutely; if the resulting $\\sin^{2}2\\theta_{13}$ or $|\\Delta m^{2}_{ee}|$ differs from the ratio-based result by more than the combined quoted uncertainty, the cancellation assumption is falsified. Alternatively, an analysis that removes the 9Li/8He background constraint or replaces the merged background spectrum with independently measured near and far spectra would reveal whether the background treatment biases the fitted parameters.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier RENO 500-day measurement establishing the measurement chain and baseline for oscillation parameters."},{"cited_title":"Nunokawa et al., Phys","cited_arxiv_id":null,"evidence_quote":"Defines the effective mass-squared difference $\\Delta m^{2}_{ee}$ used in the survival probability formula."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Documents the detector simulation, energy calibration, reconstruction, and 9Li/8He spectrum methods used here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reference reactor antineutrino spectrum model used for the expected rates and spectral shape."},{"cited_title":"Huber, Phys","cited_arxiv_id":null,"evidence_quote":"Companion reactor antineutrino spectrum model used in constructing the expected far-to-near ratio."},{"cited_title":"Abe et al","cited_arxiv_id":null,"evidence_quote":"Source of the $\\chi^{2}$ pull-term fitting method for propagating systematic uncertainties into oscillation parameters."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Measures cosmogenic 9Li and 8He production rates that underpin the updated background spectrum."}],"review_version":1}