Table of Contents
Fetching ...

Low Energy Threshold Analysis of the Phase I and Phase II Data Sets of the Sudbury Neutrino Observatory

The SNO Collaboration

TL;DR

This work reports a joint low-energy threshold analysis of SNO Phase I and Phase II data down to $T_{ m eff}=3.5$ MeV, enabling a more precise measurement of the total $^8$B solar-neutrino flux via neutral-current interactions and a direct, model-independent reconstruction of the energy-dependent $ u_e$ survival probability. The analysis employs two independent signal-extraction approaches (binned histograms and kernel estimation) and an advanced energy estimator (FTK) to improve energy resolution, while leveraging calibration data to correct Monte Carlo PDFs and reduce systematics. The results yield a high-precision NC flux $ ext{Phi}_{ m NC}$ and a precise total $^8$B flux $ ext{Phi}_{^8 m B}$, with a two-flavor fit giving $ an^2 heta_{12}=0.457^{+0.040}_{-0.029}$ and $ riangle m^2_{21}=7.59^{+0.20}_{-0.21} imes 10^{-5}$ eV$^2$, and a three-flavor fit placing an upper bound on $ ext{sin}^2 heta_{13}<0.057$ (95% C.L.). These results provide the most precise solar-neutrino oscillation constraints to date, incorporating robust cross-checks against KamLAND and other solar experiments. The LE TA framework demonstrates the value of low-threshold, high-precision solar-neutrino measurements for testing MSW physics and constraining the solar neutrino parameter space.

Abstract

Results are reported from a joint analysis of Phase I and Phase II data from the Sudbury Neutrino Observatory. The effective electron kinetic energy threshold used is T_eff=3.5 MeV, the lowest analysis threshold yet achieved with water Cherenkov detector data. In units of 10^6 cm^{-2} s^{-1}, the total flux of active-flavor neutrinos from 8B decay in the Sun measured using the neutral current (NC) reaction of neutrinos on deuterons, with no constraint on the 8B neutrino energy spectrum, is found to be Phi_NC = 5.140 ^{+0.160}_{-0.158} (stat) ^{+0.132}_{-0.117} (syst). These uncertainties are more than a factor of two smaller than previously published results. Also presented are the spectra of recoil electrons from the charged current reaction of neutrinos on deuterons and the elastic scattering of electrons. A fit to the SNO data in which the free parameters directly describe the total 8B neutrino flux and the energy-dependent nu_e survival probability provides a measure of the total 8B neutrino flux Phi_8B = 5.046 ^{+0.159}_{-0.152} (stat) ^{+0.107}_{-0.123} (syst). Combining these new results with results of all other solar experiments and the KamLAND reactor experiment yields best-fit values of the mixing parameters of theta_{12}=34.06 ^{+1.16}_{-0.84} degrees and Delta m^2_{21}=7.59 ^{+0.20}_{-0.21} x 10^{-5} eV^2. The global value of Phi_8B is extracted to a precision of ^{+2.38}_{-2.95} %. In a three-flavor analysis the best fit value of sin^2θ_{13} is 2.00 ^{+2.09}_{-1.63} x 10^{-2}. Interpreting this as a limit implies an upper bound of sin^2θ_{13} < 0.057 (95% C. L.).

Low Energy Threshold Analysis of the Phase I and Phase II Data Sets of the Sudbury Neutrino Observatory

TL;DR

This work reports a joint low-energy threshold analysis of SNO Phase I and Phase II data down to MeV, enabling a more precise measurement of the total B solar-neutrino flux via neutral-current interactions and a direct, model-independent reconstruction of the energy-dependent survival probability. The analysis employs two independent signal-extraction approaches (binned histograms and kernel estimation) and an advanced energy estimator (FTK) to improve energy resolution, while leveraging calibration data to correct Monte Carlo PDFs and reduce systematics. The results yield a high-precision NC flux and a precise total B flux , with a two-flavor fit giving and eV, and a three-flavor fit placing an upper bound on (95% C.L.). These results provide the most precise solar-neutrino oscillation constraints to date, incorporating robust cross-checks against KamLAND and other solar experiments. The LE TA framework demonstrates the value of low-threshold, high-precision solar-neutrino measurements for testing MSW physics and constraining the solar neutrino parameter space.

Abstract

Results are reported from a joint analysis of Phase I and Phase II data from the Sudbury Neutrino Observatory. The effective electron kinetic energy threshold used is T_eff=3.5 MeV, the lowest analysis threshold yet achieved with water Cherenkov detector data. In units of 10^6 cm^{-2} s^{-1}, the total flux of active-flavor neutrinos from 8B decay in the Sun measured using the neutral current (NC) reaction of neutrinos on deuterons, with no constraint on the 8B neutrino energy spectrum, is found to be Phi_NC = 5.140 ^{+0.160}_{-0.158} (stat) ^{+0.132}_{-0.117} (syst). These uncertainties are more than a factor of two smaller than previously published results. Also presented are the spectra of recoil electrons from the charged current reaction of neutrinos on deuterons and the elastic scattering of electrons. A fit to the SNO data in which the free parameters directly describe the total 8B neutrino flux and the energy-dependent nu_e survival probability provides a measure of the total 8B neutrino flux Phi_8B = 5.046 ^{+0.159}_{-0.152} (stat) ^{+0.107}_{-0.123} (syst). Combining these new results with results of all other solar experiments and the KamLAND reactor experiment yields best-fit values of the mixing parameters of theta_{12}=34.06 ^{+1.16}_{-0.84} degrees and Delta m^2_{21}=7.59 ^{+0.20}_{-0.21} x 10^{-5} eV^2. The global value of Phi_8B is extracted to a precision of ^{+2.38}_{-2.95} %. In a three-flavor analysis the best fit value of sin^2θ_{13} is 2.00 ^{+2.09}_{-1.63} x 10^{-2}. Interpreting this as a limit implies an upper bound of sin^2θ_{13} < 0.057 (95% C. L.).

Paper Structure

This paper contains 69 sections, 28 equations, 40 figures, 28 tables.

Figures (40)

  • Figure 1: (Color online) The Monte Carlo-generated distributions of (a) energy ($T_{\rm eff}$), (b) radius cubed ($R^3$), (c) direction ($\cos\theta_{\odot}$), and (d) isotropy ($\beta_{14}$) for signal events. The same simulation was used to build multi-dimensional PDFs to fit the data. In calculating $R^3$, the radius $R$ is first normalized to the 600 cm radius of the AV. The CC and NC $\cos\theta_{\odot}$ distributions are scaled by a factor of 10 for clarity against the ES peak.
  • Figure 2: (Color online) The Monte Carlo-generated distributions of (a) energy ($T_{\rm eff}$) on a log scale, (b) radius cubed ($R^3$), (c) direction ($\cos\theta_{\odot}$), and (d) isotropy ($\beta_{14}$) for background events. The same simulation was used to build multi-dimensional PDFs to fit the background events. The backgrounds shown are internal $^{214}$Bi, internal $^{208}$Tl, and AV $^{208}$Tl.
  • Figure 3: (Color online) Comparison of new model of photomultiplier angular response to data and the old model for Phase I at 365$\,$nm.
  • Figure 4: Comparison of $^{16}$N simulation to data for (a) PMT hit time-of-flight residuals and (b) photoelectron charge spectra.
  • Figure 5: Fit of the $^{16}$N Compton-electron position distribution convolved with a Gaussian to the reconstructed $z$ position of $^{16}$N data events for a typical central run in Phase II.
  • ...and 35 more figures