Measurement of Electron Antineutrino Oscillation Amplitude and Frequency via Neutron Capture on Hydrogen at Daya Bay.

This Letter reports the first measurement of the oscillation amplitude and frequency of reactor antineutrinos at Daya Bay via neutron capture on hydrogen using 1958 days of data. With over 3.6 million signal candidates, an optimized candidate selection, improved treatment of backgrounds and efficiencies, refined energy calibration, and an energy response model for the capture-on-hydrogen sensitive region, the relative ν[over ¯]_{e} rates and energy spectra variation among the near and far detectors gives sin^{2}2θ_{13}=0.0759_{-0.0049}^{+0.0050} and Δm_{32}^{2}=(2.72_{-0.15}^{+0.14})×10^{-3} eV^{2} assuming the normal neutrino mass ordering, and Δm_{32}^{2}=(-2.83_{-0.14}^{+0.15})×10^{-3} eV^{2} for the inverted neutrino mass ordering. This estimate of sin^{2}2θ_{13} is consistent with and essentially independent from the one obtained using the capture-on-gadolinium sample at Daya Bay. The combination of these two results yields sin^{2}2θ_{13}=0.0833±0.0022, which represents an 8% relative improvement in precision regarding the Daya Bay full 3158-day capture-on-gadolinium result.

First Search for Dark-Trident Processes Using the MicroBooNE Detector.

We present a first search for dark-trident scattering in a neutrino beam using a dataset corresponding to 7.2×10^{20} protons on target taken with the MicroBooNE detector at Fermilab. Proton interactions in the neutrino target at the main injector produce π^{0} and η mesons, which could decay into dark-matter (DM) particles mediated via a dark photon A^{'}. A convolutional neural network is trained to identify interactions of the DM particles in the liquid-argon time projection chamber (LArTPC) exploiting its imagelike reconstruction capability. In the absence of a DM signal, we provide limits at the 90% confidence level on the squared kinematic mixing parameter ϵ^{2} as a function of the dark-photon mass in the range 10≤M_{A^{'}}≤400 MeV. The limits cover previously unconstrained parameter space for the production of fermion or scalar DM particles χ for two benchmark models with mass ratios M_{χ}/M_{A^{'}}=0.6 and 2 and for dark fine-structure constants 0.1≤α_{D}≤1.

First Simultaneous Measurement of Differential Muon-Neutrino Charged-Current Cross Sections on Argon for Final States with and without Protons Using MicroBooNE Data.

We report the first double-differential neutrino-argon cross section measurement made simultaneously for final states with and without protons for the inclusive muon neutrino charged-current interaction channel. The proton kinematics of this channel are further explored with a differential cross section measurement as a function of the leading proton's kinetic energy that extends across the detection threshold. These measurements use data collected with the MicroBooNE detector from 6.4×10^{20} protons on target from the Fermilab booster neutrino beam with a mean neutrino energy of ∼0.8 GeV. Extensive data-driven model validation utilizing the conditional constraint formalism is employed. This motivates enlarging the uncertainties with an empirical reweighting approach to minimize the possibility of extracting biased cross section results. The extracted nominal flux-averaged cross sections are compared to widely used event generator predictions revealing severe mismodeling of final states without protons for muon neutrino charged-current interactions, possibly from insufficient treatment of final state interactions. These measurements provide a wealth of new information useful for improving event generators which will enhance the sensitivity of precision measurements in neutrino experiments.

Search for a Sub-eV Sterile Neutrino Using Daya Bay's Full Dataset.

This Letter presents results of a search for the mixing of a sub-eV sterile neutrino with three active neutrinos based on the full data sample of the Daya Bay Reactor Neutrino Experiment, collected during 3158 days of detector operation, which contains 5.55×10^{6} reactor ν[over ¯]_{e} candidates identified as inverse beta-decay interactions followed by neutron capture on gadolinium. The analysis benefits from a doubling of the statistics of our previous result and from improvements of several important systematic uncertainties. No significant oscillation due to mixing of a sub-eV sterile neutrino with active neutrinos was found. Exclusion limits are set by both Feldman-Cousins and CLs methods. Light sterile neutrino mixing with sin^{2}2θ_{14}≳0.01 can be excluded at 95% confidence level in the region of 0.01 eV^{2}≲|Δm_{41}^{2}|≲0.1 eV^{2}. This result represents the world-leading constraints in the region of 2×10^{-4} eV^{2}≲|Δm_{41}^{2}|≲0.2 eV^{2}.

Search for Heavy Neutral Leptons in Electron-Positron and Neutral-Pion Final States with the MicroBooNE Detector.

We present the first search for heavy neutral leptons (HNLs) decaying into νe^{+}e^{-} or νπ^{0} final states in a liquid-argon time projection chamber using data collected with the MicroBooNE detector. The data were recorded synchronously with the NuMI neutrino beam from Fermilab's main injector corresponding to a total exposure of 7.01×10^{20} protons on target. We set upper limits at the 90% confidence level on the mixing parameter |U_{µ4}|^{2} in the mass ranges 10≤m_{HNL}≤150 MeV for the νe^{+}e^{-} channel and 150≤m_{HNL}≤245 MeV for the νπ^{0} channel, assuming |U_{e4}|^{2}=|U_{τ4}|^{2}=0. These limits represent the most stringent constraints in the mass range 35

First Measurement of η Meson Production in Neutrino Interactions on Argon with MicroBooNE.

We present a measurement of η production from neutrino interactions on argon with the MicroBooNE detector. The modeling of resonant neutrino interactions on argon is a critical aspect of the neutrino oscillation physics program being carried out by the DUNE and Short Baseline Neutrino programs. η production in neutrino interactions provides a powerful new probe of resonant interactions, complementary to pion channels, and is particularly suited to the study of higher-order resonances beyond the Δ(1232). We measure a flux-integrated cross section for neutrino-induced η production on argon of 3.22±0.84(stat)±0.86(syst) 10^{-41} cm^{2}/nucleon. By demonstrating the successful reconstruction of the two photons resulting from η production, this analysis enables a novel calibration technique for electromagnetic showers in GeV accelerator neutrino experiments.

First Constraints on Light Sterile Neutrino Oscillations from Combined Appearance and Disappearance Searches with the MicroBooNE Detector.

We present a search for eV-scale sterile neutrino oscillations in the MicroBooNE liquid argon detector, simultaneously considering all possible appearance and disappearance effects within the 3+1 active-to-sterile neutrino oscillation framework. We analyze the neutrino candidate events for the recent measurements of charged-current ν_{e} and ν_{µ} interactions in the MicroBooNE detector, using data corresponding to an exposure of 6.37×10^{20} protons on target from the Fermilab booster neutrino beam. We observe no evidence of light sterile neutrino oscillations and derive exclusion contours at the 95% confidence level in the plane of the mass-squared splitting Δm_{41}^{2} and the sterile neutrino mixing angles θ_{µe} and θ_{ee}, excluding part of the parameter space allowed by experimental anomalies. Cancellation of ν_{e} appearance and ν_{e} disappearance effects due to the full 3+1 treatment of the analysis leads to a degeneracy when determining the oscillation parameters, which is discussed in this Letter and will be addressed by future analyses.

First Double-Differential Measurement of Kinematic Imbalance in Neutrino Interactions with the MicroBooNE Detector.

We report the first measurement of flux-integrated double-differential quasielasticlike neutrino-argon cross sections, which have been made using the Booster Neutrino Beam and the MicroBooNE detector at Fermi National Accelerator Laboratory. The data are presented as a function of kinematic imbalance variables which are sensitive to nuclear ground-state distributions and hadronic reinteraction processes. We find that the measured cross sections in different phase-space regions are sensitive to different nuclear effects. Therefore, they enable the impact of specific nuclear effects on the neutrino-nucleus interaction to be isolated more completely than was possible using previous single-differential cross section measurements. Our results provide precision data to help test and improve neutrino-nucleus interaction models. They further support ongoing neutrino-oscillation studies by establishing phase-space regions where precise reaction modeling has already been achieved.

Improved Measurement of the Evolution of the Reactor Antineutrino Flux and Spectrum at Daya Bay.

Reactor neutrino experiments play a crucial role in advancing our knowledge of neutrinos. In this Letter, the evolution of the flux and spectrum as a function of the reactor isotopic content is reported in terms of the inverse-beta-decay yield at Daya Bay with 1958 days of data and improved systematic uncertainties. These measurements are compared with two signature model predictions: the Huber-Mueller model based on the conversion method and the SM2018 model based on the summation method. The measured average flux and spectrum, as well as the flux evolution with the ^{239}Pu isotopic fraction, are inconsistent with the predictions of the Huber-Mueller model. In contrast, the SM2018 model is shown to agree with the average flux and its evolution but fails to describe the energy spectrum. Altering the predicted inverse-beta-decay spectrum from ^{239}Pu fission does not improve the agreement with the measurement for either model. The models can be brought into better agreement with the measurements if either the predicted spectrum due to ^{235}U fission is changed or the predicted ^{235}U, ^{238}U, ^{239}Pu, and ^{241}Pu spectra are changed in equal measure.

First Measurement of Quasielastic Λ Baryon Production in Muon Antineutrino Interactions in the MicroBooNE Detector.

We present the first measurement of the cross section of Cabibbo-suppressed Λ baryon production, using data collected with the MicroBooNE detector when exposed to the neutrinos from the main injector beam at the Fermi National Accelerator Laboratory. The data analyzed correspond to 2.2×10^{20} protons on target running in neutrino mode, and 4.9×10^{20} protons on target running in anti-neutrino mode. An automated selection is combined with hand scanning, with the former identifying five candidate Λ production events when the signal was unblinded, consistent with the GENIE prediction of 5.3±1.1 events. Several scanners were employed, selecting between three and five events, compared with a prediction from a blinded Monte Carlo simulation study of 3.7±1.0 events. Restricting the phase space to only include Λ baryons that decay above MicroBooNE's detection thresholds, we obtain a flux averaged cross section of 2.0_{-1.7}^{+2.2}×10^{-40} cm^{2}/Ar, where statistical and systematic uncertainties are combined.

Precision Measurement of Reactor Antineutrino Oscillation at Kilometer-Scale Baselines by Daya Bay.

We present a new determination of the smallest neutrino mixing angle θ_{13} and the mass-squared difference Δm_{32}^{2} using a final sample of 5.55×10^{6} inverse beta-decay (IBD) candidates with the final-state neutron captured on gadolinium. This sample is selected from the complete dataset obtained by the Daya Bay reactor neutrino experiment in 3158 days of operation. Compared to the previous Daya Bay results, selection of IBD candidates has been optimized, energy calibration refined, and treatment of backgrounds further improved. The resulting oscillation parameters are sin^{2}2θ_{13}=0.0851±0.0024, Δm_{32}^{2}=(2.466±0.060)×10^{-3} eV^{2} for the normal mass ordering or Δm_{32}^{2}=-(2.571±0.060)×10^{-3} eV^{2} for the inverted mass ordering.

First Measurement of High-Energy Reactor Antineutrinos at Daya Bay.

This Letter reports the first measurement of high-energy reactor antineutrinos at Daya Bay, with nearly 9000 inverse beta decay candidates in the prompt energy region of 8-12 MeV observed over 1958 days of data collection. A multivariate analysis is used to separate 2500 signal events from background statistically. The hypothesis of no reactor antineutrinos with neutrino energy above 10 MeV is rejected with a significance of 6.2 standard deviations. A 29% antineutrino flux deficit in the prompt energy region of 8-11 MeV is observed compared to a recent model prediction. We provide the unfolded antineutrino spectrum above 7 MeV as a data-based reference for other experiments. This result provides the first direct observation of the production of antineutrinos from several high-Q_{ß} isotopes in commercial reactors.

Search for an Excess of Electron Neutrino Interactions in MicroBooNE Using Multiple Final-State Topologies.

We present a measurement of ν_{e} interactions from the Fermilab Booster Neutrino Beam using the MicroBooNE liquid argon time projection chamber to address the nature of the excess of low energy interactions observed by the MiniBooNE Collaboration. Three independent ν_{e} searches are performed across multiple single electron final states, including an exclusive search for two-body scattering events with a single proton, a semi-inclusive search for pionless events, and a fully inclusive search for events containing all hadronic final states. With differing signal topologies, statistics, backgrounds, reconstruction algorithms, and analysis approaches, the results are found to be either consistent with or modestly lower than the nominal ν_{e} rate expectations from the Booster Neutrino Beam and no excess of ν_{e} events is observed.

Search for Neutrino-Induced Neutral-Current Δ Radiative Decay in MicroBooNE and a First Test of the MiniBooNE Low Energy Excess under a Single-Photon Hypothesis.

We report results from a search for neutrino-induced neutral current (NC) resonant Δ(1232) baryon production followed by Δ radiative decay, with a ⟨0.8⟩ GeV neutrino beam. Data corresponding to MicroBooNE's first three years of operations (6.80×10^{20} protons on target) are used to select single-photon events with one or zero protons and without charged leptons in the final state (1γ1p and 1γ0p, respectively). The background is constrained via an in situ high-purity measurement of NC π^{0} events, made possible via dedicated 2γ1p and 2γ0p selections. A total of 16 and 153 events are observed for the 1γ1p and 1γ0p selections, respectively, compared to a constrained background prediction of 20.5±3.65(syst) and 145.1±13.8(syst) events. The data lead to a bound on an anomalous enhancement of the normalization of NC Δ radiative decay of less than 2.3 times the predicted nominal rate for this process at the 90% confidence level (C.L.). The measurement disfavors a candidate photon interpretation of the MiniBooNE low-energy excess as a factor of 3.18 times the nominal NC Δ radiative decay rate at the 94.8% C.L., in favor of the nominal prediction, and represents a greater than 50-fold improvement over the world's best limit on single-photon production in NC interactions in the sub-GeV neutrino energy range.

Joint Determination of Reactor Antineutrino Spectra from

A joint determination of the reactor antineutrino spectra resulting from the fission of ^{235}U and ^{239}Pu has been carried out by the Daya Bay and PROSPECT Collaborations. This Letter reports the level of consistency of ^{235}U spectrum measurements from the two experiments and presents new results from a joint analysis of both data sets. The measurements are found to be consistent. The combined analysis reduces the degeneracy between the dominant ^{235}U and ^{239}Pu isotopes and improves the uncertainty of the ^{235}U spectral shape to about 3%. The ^{235}U and ^{239}Pu antineutrino energy spectra are unfolded from the jointly deconvolved reactor spectra using the Wiener-SVD unfolding method, providing a data-based reference for other reactor antineutrino experiments and other applications. This is the first measurement of the ^{235}U and ^{239}Pu spectra based on the combination of experiments at low- and highly enriched uranium reactors.

First Measurement of Energy-Dependent Inclusive Muon Neutrino Charged-Current Cross Sections on Argon with the MicroBooNE Detector.

We report a measurement of the energy-dependent total charged-current cross section σ(E_{ν}) for inclusive muon neutrinos scattering on argon, as well as measurements of flux-averaged differential cross sections as a function of muon energy and hadronic energy transfer (ν). Data corresponding to 5.3×10^{19} protons on target of exposure were collected using the MicroBooNE liquid argon time projection chamber located in the Fermilab booster neutrino beam with a mean neutrino energy of approximately 0.8 GeV. The mapping between the true neutrino energy E_{ν} and reconstructed neutrino energy E_{ν}^{rec} and between the energy transfer ν and reconstructed hadronic energy E_{had}^{rec} are validated by comparing the data and Monte Carlo (MC) predictions. In particular, the modeling of the missing hadronic energy and its associated uncertainties are verified by a new method that compares the E_{had}^{rec} distributions between data and a MC prediction after constraining the reconstructed muon kinematic distributions, energy, and polar angle to those of data. The success of this validation gives confidence that the missing energy in the MicroBooNE detector is well modeled and underpins first-time measurements of both the total cross section σ(E_{ν}) and the differential cross section dσ/dν on argon.

Search for a Higgs Portal Scalar Decaying to Electron-Positron Pairs in the MicroBooNE Detector.

We present a search for the decays of a neutral scalar boson produced by kaons decaying at rest, in the context of the Higgs portal model, using the MicroBooNE detector. We analyze data triggered in time with the Fermilab NuMI neutrino beam spill, with an exposure of 1.93×10^{20} protons on target. We look for monoenergetic scalars that come from the direction of the NuMI hadron absorber, at a distance of 100 m from the detector, and decay to electron-positron pairs. We observe one candidate event, with a standard model background prediction of 1.9±0.8. We set an upper limit on the scalar-Higgs mixing angle of θ<(3.3-4.6)×10^{-4} at the 95% confidence level for scalar boson masses in the range (100-200) MeV/c^{2}. We exclude, at the 95% confidence level, the remaining model parameters required to explain the central value of a possible excess of K_{L}^{0}→π^{0}νν[over ¯] decays reported by the KOTO collaboration. We also provide a model-independent limit on a new boson X produced in K→πX decays and decaying to e^{+}e^{-}.

New opportunities at the next-generation neutrino experiments I: BSM neutrino physics and dark matter.

The combination of the high intensity proton beam facilities and massive detectors for precision measurements of neutrino oscillation parameters including the charge-parity violating (CPV) phase will open the door to help make beyond the standard model (BSM) physics reachable even in low energy regimes in the accelerator-based experiments. Large-mass detectors with highly precise tracking and energy measurements, excellent timing resolution, and low energy thresholds will enable the searches for BSM phenomena from cosmogenic origin, as well. Therefore, it is also conceivable that BSM topics in the next-generation neutrino experiments could be the dominant physics topics in the foreseeable future, as the precision of the neutrino oscillation parameter and CPV measurements continue to improve.This paper provides a review of the current landscape of BSM theory in neutrino experiments in two selected areas of the BSM topics-dark matter and neutrino related BSM-and summarizes the current results from existing neutrino experiments to set benchmarks for both theory and experiment. This paper then provides a review of upcoming neutrino experiments throughout the next 10 to 15 year time scale and their capabilities to set the foundation for potential reach in BSM physics in the two aforementioned themes. An important outcome of this paper is to ensure theoretical and simulation tools exist to carry out studies of these new areas of physics, from the first day of the experiments, such as Deep Underground Neutrino Experiment in the U.S. and Hyper-Kamiokande Experiment in Japan.

Improved Constraints on Sterile Neutrino Mixing from Disappearance Searches in the MINOS, MINOS+, Daya Bay, and Bugey-3 Experiments.

Searches for electron antineutrino, muon neutrino, and muon antineutrino disappearance driven by sterile neutrino mixing have been carried out by the Daya Bay and MINOS+ collaborations. This Letter presents the combined results of these searches, along with exclusion results from the Bugey-3 reactor experiment, framed in a minimally extended four-neutrino scenario. Significantly improved constraints on the θ_{µe} mixing angle are derived that constitute the most constraining limits to date over five orders of magnitude in the mass-squared splitting Δm_{41}^{2}, excluding the 90% C.L. sterile-neutrino parameter space allowed by the LSND and MiniBooNE observations at 90% CL_{s} for Δm_{41}^{2}<13 eV^{2}. Furthermore, the LSND and MiniBooNE 99% C.L. allowed regions are excluded at 99% CL_{s} for Δm_{41}^{2}<1.6 eV^{2}.

First Measurement of Differential Charged Current Quasielasticlike ν_{µ}-Argon Scattering Cross Sections with the MicroBooNE Detector.

We report on the first measurement of flux-integrated single differential cross sections for charged-current (CC) muon neutrino (ν_{µ}) scattering on argon with a muon and a proton in the final state, ^{40}Ar (ν_{µ},µp)X. The measurement was carried out using the Booster Neutrino Beam at Fermi National Accelerator Laboratory and the MicroBooNE liquid argon time projection chamber detector with an exposure of 4.59×10^{19} protons on target. Events are selected to enhance the contribution of CC quasielastic (CCQE) interactions. The data are reported in terms of a total cross section as well as single differential cross sections in final state muon and proton kinematics. We measure the integrated per-nucleus CCQE-like cross section (i.e., for interactions leading to a muon, one proton, and no pions above detection threshold) of (4.93±0.76_{stat}±1.29_{sys})×10^{-38} cm^{2}, in good agreement with theoretical calculations. The single differential cross sections are also in overall good agreement with theoretical predictions, except at very forward muon scattering angles that correspond to low-momentum-transfer events.