{"id":"52a98e24-1f23-4766-b9a2-5f05f93047ac","arxiv_id":"hep-ph/0405172","paper_version":6,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Combining all 2004 neutrino oscillation data yields best-fit three-flavour parameters with sin^2(theta_13) below about 0.05 at 3 sigma and confirms that four-neutrino schemes cannot accommodate LSND alongside the rest.","lead":"A global fit of solar, atmospheric, reactor and accelerator neutrino data, updated with the 2004 KamLAND and K2K releases and new flux calculations, pins down the three-flavour oscillation parameters and tightens the bound on the last unknown mixing angle. It is the reference parameter set that the neutrino community used to plan the next generation of long-baseline and reactor experiments.","discovery_kind":"extension","skeptic_critique":{"model":"selfhosted/deepseek-v4.1-flash","headline":"The 3σ sin²θ₁₃ bound is read off a Δχ² profile whose minimum sits on the physical boundary, using a method the authors call naive; that uncalibrated step is exactly what tightens the global limit from 0.067 to 0.047.","rationale":"The reader's verdict of CONDITIONAL is the right call, and I would not move it. The paper is a careful, well-documented global fit whose numbers became a standard reference, and its statistical treatment is above the norm for the genre: exact pulls to all orders, correlated statistical errors for SNO salt, and an explicit cross-check against Super-K's own atmospheric result. The four-neutrino section is the most robust part, since the (2+2) exclusion follows from solar/atmospheric disagreement and does not depend on LSND being confirmed. I agree with the reader that the paper is unusually candid about its limits, with appendices that revise the numbers as ¹³C, MINOS and later KamLAND data arrive. Where I part company is on which soft input is load-bearing. The reader's weakest_assumption is the digitisation of the KamLAND spectrum from Fig. 3 of Ref. [15]. That is a genuine soft input, and the authors flag it in Appendix A.1, but the paper's own later appendices largely defuse it: the 2007 KamLAND sample, analysed the same way, gives sin²θ₁₃ ≤ 0.050 at 3σ versus 0.047, and Δm²₂₁ drifts only from 8.1 to 7.6 × 10⁻⁵ eV². A digitisation error large enough to matter would have shown up as a larger shift when four times the exposure arrived. The sharper unresolved issue is the one the reader mentions only in passing in the rationale: the θ₁₃ bound is obtained by the authors' own 'naive' Gaussian reading of Δχ² in a restricted physical region, with the KamLAND-only minimum at sin²θ₁₃ = −0.13. Footnote ‖ of Appendix A.2 acknowledges this and says a Monte Carlo is needed, but no such calibration is performed. Because the global best fit also sits at the boundary (Table 1, sin²θ₁₃ = 0.000), the Δχ² = 9 cut is expected to undercover, and the quoted 3σ limit is likely somewhat tight. The effect is bounded, not open-ended: Appendix D's 2007 update shows the global bound moving only from 0.040 to 0.050 at 3σ, and the authors there attribute part of the weakening to the solar+KamLAND preference for non-zero θ₁₃, which is consistent with the boundary/fluctuation effect being real but modest. So the concern lands on the precision of a quoted number, not on the qualitative conclusion that θ₁₃ is small and that solar+KamLAND now contributes materially. That is why the verdict stays CONDITIONAL rather than moving to REJECT, and why I would not ask the authors to re-prove anything cited from Ref. [15] — only to calibrate their own bound.","tokens_in":39771,"tokens_out":5001,"duration_ms":80207,"concrete_test":"Run the Monte Carlo the authors themselves propose in footnote ‖ of Appendix A.2. Generate synthetic KamLAND data sets at the global best-fit point (Δm²₂₁ = 8.1×10⁻⁵ eV², sin²θ₁₂ = 0.29, sin²θ₁₃ = 0) with the stated reactor-flux, thermal-power, fuel-composition and ¹³C/⁹Li/⁸He backgrounds, bin in 1/E_prompt exactly as in Appendix A.1, and for each set compute the 3σ upper limit on sin²θ₁₃ using the same Δχ² = 9 physical-region cut. Measure the fraction of sets whose true value lies below the quoted limit. If coverage falls below roughly 99.7%, the true 3σ limit exceeds 0.047; then repeat the global combination with the corrected solar+KamLAND profile and report both the shifted bound and the revised size of the solar+KamLAND contribution relative to CHOOZ.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim is Table 1's sin²θ₁₃ ≤ 0.047 at 3σ (1 d.o.f.), together with §4.2's assertion that solar+KamLAND data 'contribute notably' to it. Fig. 14 gives the arithmetic: CHOOZ+atm+K2K alone yields 0.067, solar+KamLAND alone 0.079, and the combination 0.047. The tightening therefore comes from combining two profiles, one of which — solar+KamLAND — is not a calibrated confidence interval. Appendix A.2 is explicit: without the sin²θ₁₃ ≥ 0 constraint the KamLAND-only best fit sits at sin²θ₁₃ = −0.13, with Δχ² = 3.5 relative to the physical point (1.75 after the ¹³C correction of Appendix B). Footnote ‖ then states 'we throughout use the naive method of calculating bounds by considering Δχ²-values as implied by Gaussian statistics and restricting the analysis to the physical region,' and notes that a reliable bound would require a Monte Carlo of synthetic data sets. That is the unresolved step. When the MLE lies on a parameter boundary, a Δχ² = 9 cut does not carry 3σ coverage; it undercovers, so 0.047 is plausibly too tight. The KamLAND resolution within the physical region quoted in Appendix A.2, σ ≈ 0.02, is not small compared with the 0.047 bound, so the boundary correction is not obviously negligible. This matters because the paper leans on the solar+KamLAND contribution precisely where CHOOZ weakens: §4.2 and Fig. 15 argue that for Δm²_atm ≲ 2×10⁻³ eV² the CHOOZ bound deteriorates and solar+KamLAND must hold the line. If that profile is over-tightened near zero, the global bound relaxes toward the CHOOZ+atm+K2K value of 0.067 and the 'contributes notably' claim weakens quantitatively. The reader's digitisation concern is real but largely defused internally: Appendix D's 2007 KamLAND sample (2881 ton-year) gives 0.050 at 3σ, essentially unchanged from 0.047, and Δm²₂₁ moves 8.1 → 7.9 → 7.6 × 10⁻⁵ eV², consistent with the real-data values.","agreement_with_reader":"partial"},"referee_report":{"model":"selfhosted/deepseek-v4.1-flash","summary":"This manuscript is a review-and-update of the authors' global oscillation analyses, combining solar (Homestake, SAGE, GALLEX/GNO, Super-K, SNO salt and day/night), KamLAND 2004, atmospheric Super-K (with the new three-dimensional Honda flux), K2K-I+II, and CHOOZ data. Section 2 gives two-flavour atmospheric+K2K regions and tests robustness against flavour-changing and non-universal non-standard interactions in Earth matter; Section 3 gives two-flavour solar+KamLAND regions and tests robustness against solar density fluctuations, convective-zone magnetic fields, and neutrino transition magnetic moments; Section 4 presents the five-parameter three-flavour fit (Table 1) and discusses the small parameters alpha and sin^2 theta_13; Section 5 updates the four-neutrino (3+1) and (2+2) interpretations of LSND using a parameter-goodness-of-fit statistic. Appendices A and B document the KamLAND reanalysis and the effect of the 13C(alpha,n)16O background, and Appendices C and D append 2006 and 2007 updates including MINOS, the final K2K sample, the BS05(OP) solar model, and the 2881 ton-year KamLAND sample.","tokens_in":40212,"tokens_out":7640,"duration_ms":140653,"significance":"The paper is a careful and useful synthesis. The statistical machinery is described in enough detail to be reproduced: pull parameters for systematics, an exact treatment of pulls to all orders rather than a first-order expansion (which lets the SNO neutral-current measurement and the SSM boron-flux prediction be fitted on the same footing), correlated statistical errors for the SNO salt data, and an explicit cross-check that the atmospheric fit agrees with the Super-K collaboration's own result. The robustness studies are genuine tests rather than assertions: the NSI fit of Fig. 4 shows the oscillation parameters are stable while bounding epsilon and epsilon-prime, and the solar-fluctuation and transition-moment analyses yield falsifiable limits (|Lambda| < 3.4 x 10^-10 mu_B, and 0.29 x 10^-10 mu_B after three years of Borexino). The parameter-goodness-of-fit treatment of Section 5 is a reproducible and appropriate way to reject (2+2) schemes without dilution by data-point counting. If the numbers hold, the paper fixes the size of the last undetermined leptonic mixing angle and hence the scale of CP-violating effects for the next generation of long-baseline experiments.","major_comments":[{"comment":"The headline bound sin^2 theta_13 <= 0.047 at 3 sigma (1 d.o.f.) is obtained by applying a Delta chi^2 = 9 cut to a profile whose minimum lies on the physical boundary. The authors state this plainly in Appendix A.2 (the unconstrained KamLAND best fit sits at sin^2 theta_13 = -0.13, with Delta chi^2 = 3.5 relative to the physical point, 1.75 after the 13C correction) and in the footnote, which calls the procedure naive and notes that a reliable bound would require a Monte Carlo of synthetic data sets. The limitation is acknowledged, but it is not carried through to the presentation: Eq. (17) and Table 1 quote 0.047 without qualification, and Section 4.2 and Fig. 15 then use the solar+KamLAND constraint as the argument that the bound is held where CHOOZ weakens (Delta m^2_atm below about 2 x 10^-3 eV^2). Because a boundary MLE undercovers, 0.047 is plausibly too tight, and the KamLAND resolution within the physical region quoted in Appendix A.2 (sigma about 0.02) is not small compared with the bound. The authors' own Appendix B demonstrates the sensitivity concretely: correcting the 13C background moves the 3 sigma value from 0.047 to 0.051 (Table B1). I ask that the bound be quoted with the calibration caveat alongside Eq. (17), or that the 13C-updated value be promoted to the primary result, or that the Monte Carlo calibration be supplied. The qualitative conclusion that theta_13 is small survives either way; the specific number does not currently carry the stated coverage.","section":"Section 4.2, Eq. (17), Table 1, and footnote of Appendix A.2","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The KamLAND spectral information is reconstructed by digitising the ratio plot of Fig. 3 of Ref. [15] into 13 bins of 180 km/E and re-binning in 1/E_prompt, because event-level energies are unavailable outside the collaboration. The authors justify the 1/E binning well, but they also show that the choice is not innocuous: the Delta chi^2 of the low-Delta m^2 minimum changes from 4.0 to 7.5 and that of the high-LMA minimum from 8.5 to 11.3 when going from E_prompt to 1/E_prompt binning. What is not shown is the effect of this choice on the two quantities that the rest of the paper leans on, namely the Delta m^2_sol error bar and the solar+KamLAND contribution to the theta_13 bound. Since the latter is one of the paper's headline results, a systematic check would be proportionate: quote the theta_13 bound and the Delta m^2_sol interval under both binnings, or state explicitly how much of the quoted precision is attributable to the digitisation and binning choice. This is a request for a sensitivity statement, not for a reanalysis.","section":"Appendix A.1 and Fig. A1","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The paper's stated purpose is to report the status of global fits, but the main text and abstract present the 2004 numbers as current while Appendices C and D (June 2006 and September 2007) supersede them. Table D1 gives sin^2 theta_13 <= 0.050 at 3 sigma, Delta m^2_31 = 2.4 x 10^-3 eV^2, and sin^2 theta_12 = 0.32, against Table 1's 0.047, 2.2 x 10^-3 eV^2, and 0.30. The abstract's phrase 'the latest KamLAND and K2K updates presented at Neutrino2004' is now factually stale given Appendix D. For a review whose title is 'Status of global fits', the reader is left to reconcile five tables spanning three years of data. I ask that the abstract and introduction be updated to point to the latest appendix, or that a short signposting paragraph be added at the start of Section 4 stating which table is the current result and which are historical. This is a currency and presentation issue rather than a physics error, but it bears directly on what the paper claims to deliver.","section":"Table 1 versus Tables B1, C1, and D1; Abstract and Section 1","authors_acknowledge":null,"authors_acknowledge_p":null}],"minor_comments":[{"comment":"The phrase 'the current worlds' global neutrino data sample' should read 'the current world's global neutrino data sample'.","section":"Abstract","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"Eq. (8) gives sin^2 theta_sol = 0.29 from solar+KamLAND, while Table 1 gives sin^2 theta_12 = 0.30 from the global three-flavour fit. A sentence explaining that the difference is the two- versus three-flavour treatment (and the inclusion of atmospheric, K2K, and CHOOZ data) would prevent the reader from reading this as an inconsistency.","section":"Section 3.3, Eq. (8) and Table 1","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The caption states that the regions delimited by the lines correspond to solar data, but does not say whether the same solar contours are overlaid in both the 2002 and 2004 panels. Please state this explicitly, since the visual comparison of the two panels is the point of the figure.","section":"Fig. 6 caption","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The legend lists 'Solar' and 'KamLAND' as separate entries alongside 'Total', but the surrounding text discusses 'solar+KamLAND' as a single data set. Please clarify whether the figure shows solar-only and KamLAND-only profiles in addition to their combination, and label them accordingly.","section":"Fig. 14 legend","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"Section 5.2 quotes chi^2_PG = 26.1 corresponding to 5.1 sigma, which is a one-degree-of-freedom test in eta_s, whereas Table 2 quotes chi^2_PG = 33.9 for (2+2) with four degrees of freedom and 4.9 sigma. The two numbers are consistent but the text does not distinguish them, and a reader may take them as contradictory. A clause noting that the first is the one-parameter eta_s test of Eq. (21) and the second the four-parameter test of Eq. (22) would remove the ambiguity.","section":"Section 5.2 and Table 2","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"Typographical: 'No significant improvement has been otained here' should read 'obtained'.","section":"Appendix D","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The NSI bounds -0.03 <= epsilon <= 0.02 and |epsilon'| <= 0.05 are taken from Ref. [50], which used the one-dimensional Bartol flux, while the rest of Section 2 now uses the three-dimensional Honda flux. The footnote acknowledges this, but a sentence in the main text stating that the NSI bounds are not on the same flux footing as the oscillation parameters would help the reader weigh them.","section":"Section 2.4, Eq. (5)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The 2.69 accidental background events added to the lowest energy bin are quoted without an uncertainty. Please state the uncertainty or explain why it is negligible.","section":"Appendix A.1","authors_acknowledge":null,"authors_acknowledge_p":null}],"recommendation":"minor_revision","confidential_remarks":"The paper is a well-executed review and the analysis machinery is sound. My reservation is narrow: the sin^2 theta_13 bound of 0.047 is quoted as a calibrated 3 sigma limit although the authors themselves describe the method as naive and show in Appendix B that a background correction moves it to 0.051. The fix is textual, not computational, since the corrected number already exists in Table B1. The appended 2006 and 2007 updates also leave the abstract describing 2004 data; the journal may wish to decide whether the published version should carry these appendices or whether they belong in a separate update. The self-citation to Refs. [30], [36], [121]-[123], and [125] is appropriate, as these are the published analyses being updated, but a short table mapping each section to the prior analysis it supersedes would help the reader."},"author_rebuttal":null,"desk_editor":{"model":"selfhosted/deepseek-v4.1-flash","letter":"Bottom line: this is the reference global fit for the 2004–2007 era, and the parameter tables are the numbers people actually cited. It is labelled a review, but the fits are new work, not a summary of prior results. Cite it, and cite the appendix version you mean — the paper revises itself in place four times (Tables 1, B1, C1, D1). That is honest but makes a bare 'Table 1' ambiguous, which is my only real usability complaint. Minor, but worth knowing before you put a number in a slide. What is genuinely new: the three-flavour fit with KamLAND 2004, K2K-I+II, three-dimensional atmospheric fluxes and BP04/BS05 solar fluxes; the decomposition of the sin²θ13 bound by data set; and the (2+2) four-neutrino exclusion at 4.9σ via parameter goodness of fit, which holds whether or not LSND is confirmed. That last section is the most robust thing in the paper and does not depend on the contested inputs. The statistical work is better than typical for this genre: pulls treated exactly to all orders rather than expanded, correlated statistical errors for SNO salt, and an explicit cross-check that the atmospheric fit reproduces Super-K's own result. The appendices are unusually candid about their own limits. On the soft spots, the stress-test lands on the mechanism but not on the conclusion. Appendix A.2 footnote ‖ admits the θ13 bound uses 'the naive method' of a Δχ² cut restricted to the physical region, while the KamLAND-only best fit sits at sin²θ13 = −0.13. At a boundary that cut does not carry 3σ coverage, so 0.047 is plausibly too tight and should not be quoted as a calibrated interval. But the 2007 KamLAND sample gives 0.050, so the physics conclusion — θ13 small, solar+KamLAND holding the line where CHOOZ weakens — survives intact. The digitisation of the KamLAND spectrum from a figure is disclosed and largely defused by Appendix D's stability check. Who should read it: anyone doing oscillation phenomenology or long-baseline design in this parameter range. It deserves a serious referee, though the referee's job is mostly to check the boundary treatment and the binning, not the framework.","headline":"The standard 2004–2007 global neutrino fit: tables are genuinely new work, the θ13 bound is a Δχ² cut rather than a calibrated interval, and the authors say so themselves.","tokens_in":41031,"tokens_out":1351,"would_cite":true,"duration_ms":23725,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"pith_extraction":{"msc":[],"pacs":["26.65.+t","13.15.+g","14.60.Pq","95.55.Vj"],"model":"selfhosted/deepseek-v4.1-flash","headline":"A single five-parameter fit to all neutrino oscillation data fixes the size of the last unknown mixing angle and with it the scale of CP violation in the lepton sector.","keywords":["neutrino oscillations","global fit","three-flavour mixing","theta_13","mass hierarchy ratio","KamLAND","solar neutrinos","sterile neutrinos"],"falsifier":"A direct reactor measurement of $sin^2theta_{13}$ above roughly 0.05, or a KamLAND reanalysis with event-level data that shifts the spectral distortion and hence $Delta m^2_{sol}$, would settle whether the quoted bound and its solar-plus-KamLAND attribution hold.","tokens_in":41066,"feed_emoji":"⚛️","tokens_out":2927,"duration_ms":50653,"temperature":0.7,"pith_summary":"Three-flavour neutrino oscillations are described by five numbers: two mass-squared differences, three mixing angles. This paper's job is to extract those five numbers, with honest error bars, from every oscillation data set available in mid-2004 — solar rates and spectra, the KamLAND reactor spectrum, Super-K atmospheric zenith distributions, K2K, and the CHOOZ reactor limit — and to say what the result implies for the two small quantities that control whether neutrinos can violate CP. Those quantities are the hierarchy ratio $alpha = Delta m^2_{sol}/Delta m^2_{atm}$ and the mixing angle $theta_{13}$. The paper reports $alpha = 0.035$ with $0.024 leq alpha leq 0.060$ at $3sigma$, and $sin^2theta_{13} leq 0.047$ at $3sigma$ for one degree of freedom. The second number is the one that matters for planning: it says the CP-violating term in the oscillation probability is small, so any future experiment that wants to see leptonic CP violation has to be built to reach it.","feed_headline":"Neutrino fits pin the last mixing angle below 0.05","feed_subtitle":"A five-parameter fit to all 2004 solar, atmospheric, reactor and accelerator data sets the scale of future CP violation.","key_machinery":"The load-bearing device is a five-dimensional $chi^2$ built from a pull-parameter treatment of every systematic uncertainty, minimised over $(sin^2theta_{12}, sin^2theta_{23}, sin^2theta_{13}, Delta m^2_{21}, Delta m^2_{31})$ and projected onto one- and two-dimensional subspaces. Two approximations make it tractable: $Delta m^2_{sol}$ is set to zero when fitting atmospheric and K2K data, and $Delta m^2_{atm}$ is set to infinity when fitting solar and KamLAND data, which also makes the Dirac CP phase drop out. The pull method is extended in two ways the paper flags as new: it handles correlated statistical errors, needed for the SNO salt data, and it treats each pull parameter exactly rather than to first order, which lets the fit choose its own compromise between the SNO neutral-current measurement and the Standard Solar Model prediction for the boron flux instead of expanding around the predicted value. A separate parameter-goodness-of-fit statistic, built from the sum of the individual data sets' $Delta chi^2$ functions rather than from the total $chi^2$, is what converts the solar-versus-atmospheric tension in four-neutrino schemes into a meaningful rejection level without being diluted by the large number of data points.","core_discovery":"The central claim is that a global three-flavour fit to the 2004 data set returns a consistent, tightly constrained parameter set, and that the bound on $sin^2theta_{13}$ is no longer a CHOOZ monopoly. The best-fit points are $Delta m^2_{21} = 8.1 times 10^{-5}$ eV$^2$, $Delta m^2_{31} = 2.2 times 10^{-3}$ eV$^2$, $sin^2theta_{12} = 0.30$, $sin^2theta_{23} = 0.50$, with $sin^2theta_{13}$ compatible with zero and bounded above at 0.047 ($3sigma$, 1 d.o.f.). The argument for the $theta_{13}$ bound is the paper's most consequential move: the authors show that solar plus KamLAND data alone now give $sin^2theta_{13} leq 0.079$ at $3sigma$, comparable to the CHOOZ-plus-atmospheric-plus-K2K bound of 0.067, and that the two constraints are complementary rather than redundant. Solar data disfavour the large $Delta m^2_{sol}$ that would let a nonzero $theta_{13}$ hide, while KamLAND's spectral distortion pins $Delta m^2_{sol}$ down and breaks the $sin^2theta_{13}$–$sin^2theta_{12}$ degeneracy. The paper also uses the same machinery to rule out four-neutrino schemes: (2+2) spectra fail at $4.9sigma$ from the disagreement between solar and atmospheric data on how much sterile neutrino participates, and (3+1) schemes are strongly disfavoured by short-baseline disappearance data, leaving the LSND anomaly without a satisfactory oscillation explanation.","pith_inferences":["The paper's own later appendices show how much the headline numbers move with new data: by the 2007 update the $3sigma$ bound on $sin^2theta_{13}$ relaxes from 0.047 to 0.050 and the best fit shifts to a slightly nonzero 0.007. A reader should treat the 2004 numbers as a snapshot of a fit that was still being re-run, not as a fixed constant of nature.","The claim that solar plus KamLAND contributes materially to the $theta_{13}$ bound depends on a spectral-shape argument in KamLAND rather than on a rate. That makes the bound sensitive to the shape of the reactor flux and the energy-scale calibration in a way that the CHOOZ bound is not, so the two constraints are less independent than their separate listing suggests.","If the fit prefers a slightly negative $sin^2theta_{13}$ when the physical region is not imposed — as Appendix A.2 reports, at $Delta chi^2 approx 3.5$ before the $^{13}$C background correction — then the quoted upper bound is being set partly by a downward fluctuation, and a future experiment with better statistics could move the central value up without contradicting anything here.","The pattern that emerges, with two large angles, one small angle, and a hierarchy ratio near 0.03, is exactly the regime in which the CP-violating asymmetry is largest relative to the matter-induced fake asymmetry, so the smallness of $theta_{13}$ is bad news for statistics but not necessarily for the cleanliness of a CP measurement."],"forward_implications":["If $sin^2theta_{13} lesssim 0.05$, the CP-violating term in the long-baseline $nu_e to nu_mu$ appearance probability is suppressed by roughly this factor, so a next-generation superbeam or beta-beam must be designed around a small $theta_{13}$, not a large one.","The lower bound on $Delta m^2_{atm}$ from K2K, combined with the $Delta m^2_{atm}$-independent solar-plus-KamLAND bound on $theta_{13}$, prevents the CHOOZ limit from degrading at the low $Delta m^2_{atm}$ values implied by the new three-dimensional atmospheric fluxes.","Because $alpha$ is measured to be about 0.03, the sub-leading three-flavour corrections to atmospheric and long-baseline oscillations are small but not negligible, and the combination $alpha sin 2theta_{12} = 0.032$ (range 0.022 to 0.054) is the quantity a long-baseline experiment actually constrains.","The (2+2) four-neutrino explanation of LSND is excluded independently of whether LSND itself is confirmed, because the rejection comes from solar and atmospheric data disagreeing about the sterile fraction.","The oscillation interpretation of solar data survives departures from the Standard Solar Model, from random radiative-zone density fluctuations, and from sub-leading spin-flavour precession, up to the magnetic moment bounds the paper derives."],"supporting_citations":[{"why":"The 2004 KamLAND reactor data, whose spectral distortion supplies the improved $Delta m^2_{sol}$ determination and the solar-plus-KamLAND contribution to the $theta_{13}$ bound.","marker":"[15]"},{"why":"The three-dimensional atmospheric neutrino flux calculation that shifts the atmospheric mass splitting downward relative to the older one-dimensional fluxes.","marker":"[32]"},{"why":"The BP04 Standard Solar Model update, used to demonstrate that the solar oscillation parameters are stable against changes in the flux predictions.","marker":"[31]"},{"why":"The SNO salt-phase neutral-current, charged-current and elastic-scattering fluxes that tighten the bound on a sterile component in solar oscillations.","marker":"[10]"},{"why":"The K2K-I plus K2K-II single-ring muon spectrum, which provides the independent accelerator confirmation of $Delta m^2_{atm}$ and the lower bound on it.","marker":"[17]"},{"why":"The CHOOZ reactor limit, the baseline constraint on $sin^2theta_{13}$ against which the new solar-plus-KamLAND bound is compared.","marker":"[105]"},{"why":"The pull-parameter $chi^2$ method that the paper extends to correlated statistical errors and to all-orders treatment of the pull parameters.","marker":"[61]"},{"why":"The parameter-goodness-of-fit statistic used to convert the disagreement between data sets in four-neutrino schemes into a rejection level.","marker":"[125]"},{"why":"The fit of atmospheric data with flavour-changing and non-universal neutral-current interactions, used to show the oscillation parameters are stable against non-standard physics.","marker":"[50]"}],"fun_headline_variants":["2004 global fit bounds sin^2 theta13 at 0.047 (3 sigma, 1 d.o.f.)","Solar plus KamLAND data rival CHOOZ on the theta13 bound","Four-neutrino LSND schemes disfavoured at 4.9 sigma by solar and atmospheric data","Global fit fixes five three-flavour parameters; theta13 stays consistent with zero","KamLAND spectrum breaks the theta13-theta12 degeneracy in the solar fit"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The KamLAND spectral information is read off a published ratio plot rather than from event-level data, so every KamLAND-derived number in the paper, including the $Delta m^2_{sol}$ best fit and the solar-plus-KamLAND share of the $theta_{13}$ bound, rests on that digitisation.","fun_headline_variants_meta":{"raw":{"variants":["2004 global fit bounds sin^2 theta13 at 0.047 (3 sigma, 1 d.o.f.)","Solar plus KamLAND data rival CHOOZ on the theta13 bound","Four-neutrino LSND schemes disfavoured at 4.9 sigma by solar and atmospheric data","Global fit fixes five three-flavour parameters; theta13 stays consistent with zero","KamLAND spectrum breaks the theta13-theta12 degeneracy in the solar fit"]},"model":"selfhosted/deepseek-v4.1-flash","effort":"low","attempts":1,"cost_usd":0.0,"raw_usage":{"total_tokens":1637,"prompt_tokens":1090,"completion_tokens":547,"prompt_tokens_details":null,"completion_tokens_details":{"reasoning_tokens":443}},"tokens_in":1090,"tokens_out":547,"duration_ms":8522,"image_pages":[],"temperature":1.0,"reasoning_budget":50,"reasoning_tokens":443,"transport_provider":"selfhosted","cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-10-29T16:46:00.885875+00:00","model_set":{"reader":"selfhosted/deepseek-v4.1-flash"},"falsifier":"A direct reactor measurement of $sin^2theta_{13}$ above roughly 0.05, or a KamLAND reanalysis with event-level data that shifts the spectral distortion and hence $Delta m^2_{sol}$, would settle whether the quoted bound and its solar-plus-KamLAND attribution hold.","supporting_citations":[{"cited_title":"Measurement of Neutrino Oscillation with KamLAND: Evidence of Spectral Distortion","cited_arxiv_id":"hep-ex/0406035","evidence_quote":"The 2004 KamLAND reactor data, whose spectral distortion supplies the improved $Delta m^2_{sol}$ determination and the solar-plus-KamLAND contribution to the $theta_{13}$ bound."},{"cited_title":"A New calculation of the atmospheric neutrino flux in a 3-dimensional scheme","cited_arxiv_id":"astro-ph/0404457","evidence_quote":"The three-dimensional atmospheric neutrino flux calculation that shifts the atmospheric mass splitting downward relative to the older one-dimensional fluxes."},{"cited_title":"What do we (not) know theoretically about solar neutrino fluxes?","cited_arxiv_id":"astro-ph/0402114","evidence_quote":"The BP04 Standard Solar Model update, used to demonstrate that the solar oscillation parameters are stable against changes in the flux predictions."},{"cited_title":"Measurement of the Total Active 8B Solar Neutrino Flux at the Sudbury Neutrino Observatory with Enhanced Neutral Current Sensitivity","cited_arxiv_id":"nucl-ex/0309004","evidence_quote":"The SNO salt-phase neutral-current, charged-current and elastic-scattering fluxes that tighten the bound on a sterile component in solar oscillations."},{"cited_title":"Limits on Neutrino Oscillations from the CHOOZ Experiment","cited_arxiv_id":"hep-ex/9907037","evidence_quote":"The CHOOZ reactor limit, the baseline constraint on $sin^2theta_{13}$ against which the new solar-plus-KamLAND bound is compared."},{"cited_title":"Getting the most from the statistical analysis of solar neutrino oscillations","cited_arxiv_id":"hep-ph/0206162","evidence_quote":"The pull-parameter $chi^2$ method that the paper extends to correlated statistical errors and to all-orders treatment of the pull parameters."},{"cited_title":"Testing the statistical compatibility of independent data sets","cited_arxiv_id":"hep-ph/0304176","evidence_quote":"The parameter-goodness-of-fit statistic used to convert the disagreement between data sets in four-neutrino schemes into a rejection level."},{"cited_title":"Probing neutrino non-standard interactions with atmospheric neutrino data","cited_arxiv_id":"hep-ph/0108043","evidence_quote":"The fit of atmospheric data with flavour-changing and non-universal neutral-current interactions, used to show the oscillation parameters are stable against non-standard physics."}],"review_version":1}