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.

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 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 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 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.

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 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.

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^{-}.

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.

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 Inclusive Muon Neutrino Charged Current Differential Cross Sections on Argon at E_{ν}∼0.8 GeV with the MicroBooNE Detector.

We report the first measurement of the double-differential and total muon neutrino charged current inclusive cross sections on argon at a mean neutrino energy of 0.8 GeV. Data were collected using the MicroBooNE liquid argon time projection chamber located in the Fermilab Booster neutrino beam and correspond to 1.6×10^{20} protons on target of exposure. The measured differential cross sections are presented as a function of muon momentum, using multiple Coulomb scattering as a momentum measurement technique, and the muon angle with respect to the beam direction. We compare the measured cross sections to multiple neutrino event generators and find better agreement with those containing more complete treatment of quasielastic scattering processes at low Q^{2}. The total flux integrated cross section is measured to be 0.693±0.010(stat)±0.165(syst)×10^{-38} cm^{2}.

Extraction of the

This Letter reports the first extraction of individual antineutrino spectra from ^{235}U and ^{239}Pu fission and an improved measurement of the prompt energy spectrum of reactor antineutrinos at Daya Bay. The analysis uses 3.5×10^{6} inverse beta-decay candidates in four near antineutrino detectors in 1958 days. The individual antineutrino spectra of the two dominant isotopes, ^{235}U and ^{239}Pu, are extracted using the evolution of the prompt spectrum as a function of the isotope fission fractions. In the energy window of 4-6 MeV, a 7% (9%) excess of events is observed for the ^{235}U (^{239}Pu) spectrum compared with the normalized Huber-Mueller model prediction. The significance of discrepancy is 4.0σ for ^{235}U spectral shape compared with the Huber-Mueller model prediction. The shape of the measured inverse beta-decay prompt energy spectrum disagrees with the prediction of the Huber-Mueller model at 5.3σ. In the energy range of 4-6 MeV, a maximal local discrepancy of 6.3σ is observed.

Search for Sterile Neutrinos in MINOS and MINOS+ Using a Two-Detector Fit.

A search for mixing between active neutrinos and light sterile neutrinos has been performed by looking for muon neutrino disappearance in two detectors at baselines of 1.04 and 735 km, using a combined MINOS and MINOS+ exposure of 16.36×10^{20} protons on target. A simultaneous fit to the charged-current muon neutrino and neutral-current neutrino energy spectra in the two detectors yields no evidence for sterile neutrino mixing using a 3+1 model. The most stringent limit to date is set on the mixing parameter sin^{2}θ_{24} for most values of the sterile neutrino mass splitting Δm_{41}^{2}>10^{-4} eV^{2}.

Measurement of the Electron Antineutrino Oscillation with 1958 Days of Operation at Daya Bay.

We report a measurement of electron antineutrino oscillation from the Daya Bay Reactor Neutrino Experiment with nearly 4 million reactor ν[over ¯]_{e} inverse ß decay candidates observed over 1958 days of data collection. The installation of a flash analog-to-digital converter readout system and a special calibration campaign using different source enclosures reduce uncertainties in the absolute energy calibration to less than 0.5% for visible energies larger than 2 MeV. The uncertainty in the cosmogenic ^{9}Li and ^{8}He background is reduced from 45% to 30% in the near detectors. A detailed investigation of the spent nuclear fuel history improves its uncertainty from 100% to 30%. Analysis of the relative ν[over ¯]_{e} rates and energy spectra among detectors yields sin^{2}2θ_{13}=0.0856±0.0029 and Δm_{32}^{2}=(2.471_{-0.070}^{+0.068})×10^{-3} eV^{2} assuming the normal hierarchy, and Δm_{32}^{2}=-(2.575_{-0.070}^{+0.068})×10^{-3} eV^{2} assuming the inverted hierarchy.